应用于双麦克风小型手持设备的语音降噪方法

应用于双麦克风小型手持设备的语音降噪方法
【专利摘要】本发明公开了一种应用于双麦克风小型手持设备的语音降噪方法，该方法可以根据用户所使用模式的不同来选取对应的降噪方法；在正常手持通话模式下，基于功率水平差异的方法适用于该模式下的近场信号模型，且能有效利用主次麦克风接收语音信号在能量上的较大差异信息；在使用扬声器模式下，利用基于双麦克风之间相关性的方法不需要进行噪声功率谱的估计，且能够避免两个麦克风接收噪声不一致对于基于功率水平差异方法的影响；此外，相比于传统的基于功率水平差异的方法，本发明通过构造的维纳滤波能有效避免噪声功率谱估计错误导致的音乐噪声。
【专利说明】
应用于双麦克风小型手持设备的语音降噪方法
技术领域
[0001] 本发明涉及语音降噪技术领域，尤其涉及一种应用于双麦克风小型手持设备的语 音降噪方法。
【背景技术】
[0002] 在现有的小型手持设备中，常采用双麦克风以及三麦克风。现有技术利用主麦克 风接收期望语音及背景噪声，次麦克风接收背景噪声，并且假设这两个麦克风接收到的背 景噪声是一致的，从而利用两个麦克风接收信号的差值作为期望信号输出。
[0003] 这种方法在背景噪声不一致的情况下会产生音乐噪声，并且当次麦克风包含期望 语音成分的时候会产生语音失真。此外，利用相位差(时延估计）的双通道滤波技术会残留 音乐噪声，且由相位差计算声源方向的可行性是建立在远场条件下的，并不适用于手持近 场的环境。

【发明内容】

[0004] 本发明的目的是提供一种应用于双麦克风小型手持设备的语音降噪方法，可以尽 可能无失真地提取出期望语音，消除背景噪声。
[0005] 本发明的目的是通过以下技术方案实现的：
[0006] -种应用于双麦克风小型手持设备的语音降噪方法，其特征在于，包括：
[0007] 接收双麦克风的语音信号，并进行通话模式的判别；
[0008] 若当前为手持通话模式，则采用基于双麦克风功率水平差异的语音降噪方法来获 取降噪后的语音信号；
[0009] 若当前为扬声器通话模式，则采用基于双麦克风之间相干性的语音降噪方法来获 取降噪后的语音信号。
[0010] 2、根据权利要求1所述的方法，其特征在于，所述采用基于双麦克风功率水平差异 的语音降噪方法来获取降噪后的语音信号包括：
[0011] 主麦克风与次麦克风接收到的语音信号如下：
[0012] xi(m) =hi(m)*s(m)+m(m);
[0013] X2(m) =h2(m)*s(m)+n2(m);
[0014] 其中，m为采样点，Xl(m)为主麦克风接收到的语音信号，X2(m)为次麦克风接收到的 语音信号;hi(m)为声波传播模型冲击响应，m(m)为噪声，i = l，2;s(m)为目标语音，*为卷积 操作；
[0015] 对主麦克风与次麦克风接受到的语音信号分别做短时傅里叶变换，得到：
[0016] Xi(n,k)=Hi(n,k)S(n,k)+Ni(n,k)；
[0017] X2(n，k)=H2(n，k)S(n，k)+N2(n，k);
[0018] 其中，η和k分别表示时间点和频率点;将上述两个式子改写为：
[0019] Xi(n,k) = Si(n,k)+Ni(n,k)；
[0020] X2(n，k)=Hi2(n，k)Si(n，k)+N2(n，k);
[0021] 其中，SKhk)表示出(11，1〇5(11，1〇，从而有
：根据短时傅里叶变 换结果分别计算主麦克风与次麦克风带噪音的功率谱密度PSD，得到：
[0022]
[0023]
[0024]
[0025]
[0026]
[0027]
[0028] 利用估计的语音信号PSD和噪声信号PSD可以构造维纳滤波器G^(n，k);其中利用 估计的主麦克风语音信号PSD& (? j)和噪声信号PSD^, (?,幻进行构造时：
[0029]
[0030] 其中，G^(n，k)的下标ΔΡ表示该滤波器是基于功率水平差得到的;将| APx(n，k) 表达式带入上式，得到：
[0031]
[0032]在上式中加入一个自由参数α，则有：
[0033] 十'
一 ' · · -I
[0034] 其中，主麦克风噪声的PSD采用语音信号开始的Τ纯噪声帧计算，公式如下：
[0035]
[0036] 式中，λΝ为噪声遗忘因子;ΧΚη,?Ο表示主麦克风接收信号的时频域值；
[0037]主麦克风与次麦克风带噪音的互功率谱密度CPSD为：
[0038]
[0039] 其中幻为两个麦克风接收噪声信号的CPSD，通过下式进行估计：
[0040]
[0041 ]及(》，幻通过主麦克风接收信号PSD与估计噪声PSD的差值进行估计：
[0042]
[0043]从而得到估计的冲击响应函数H12(n，k):
[0044]
[0045] 两个麦克风接收信号的PSD和CPSD，即\ (/α)和巧1:1.2 ，采用如下的递归平均 方法来估计：
[0046]
[0047]
[0048] 其中λχ为带噪语音遗忘因子。将主麦克风以及次麦克风做短时傅里叶变换后的语 音信号与维纳滤波器相乘，并做离散傅里叶逆变换处理及叠接相加，从而得到降噪之后的 时域语音信号。
[0049] 3、根据权利要求1所述的方法，其特征在于，所述采用基于双麦克风之间相干性的 语音降噪方法来获取降噪后的语音信号包括：
[0050] 主麦克风与次麦克风接收到的语音信号如下：
[0051] xi(m) = si(m)+m(m)，i = 1，2
[0052] 其中，m为采样点，X1(m)为主麦克风接收到的语音信号，x2(m)为次麦克风接收到的 语音信号，m(m)为噪声， Sl(m)为目标语音；
[0053]进行短时傅里叶变换后得到：
[0054] Xi(n,k) = Si(n,k)+Ni(n,k), i = 1,2；
[0055] 其中，n和k分别表示时间点和频率点；
[0056] 定义主麦克风与次麦克风接收到的语音信号的相干性函数为：
[0057]
[0058] 式中，幻分别为主麦克风与次麦克风带噪音的功率谱密度PSD， 八为主麦克风与次麦克风带噪语音的互功率谱密度；
[0059] 相干件函教与Φ砉克风与汝砉克风届部信嚙bh SNR!和3爾9有加下关系：
[0060]
[0061] 其中，和分别表示两个麦克风接收的目标语音相干性函数和噪声相干性 函数;4
則上式改写为：
[0062]
[0063] 当使用扬声器通话模式时，假设目标语音信号源在双麦克的正前方，选取两麦克 风位置的中心作为阵列参考点，则目标语音方向即为〇°，背景噪声信号源等效为从Θ方向入 射，贝主麦克风与次麦克风接收到的信号Ul和U2之间的相干性为：
[0064] Γ响 ;
[0065]式中，fs是采样率，d是主麦克风与次麦克风的间距，c为声速;山和1]2表示主麦克风 与次麦克风接收到的语音信号或者噪声信号；
[0066]将上面两个表达式相结合，得到：
[0067]
[0068] 其中，τ = fs (d/c)，分别取出实部和虚部：
[0069]
[0070]
[0071 ] 其中，α= ω TSin0，分别对实部和虚部表达式变型，得到：
[0074]由于上述两式相等，则有：[0075] Icosa = (R-l )sina+I ；[0076] 根据 cos2a+sin2a = 1，得到：[0077] (I2+(l-R)2)sin2a+2I(R-l)sina = 〇；[0078] 其根为：
[0072]
[0073]
[0079]
[0080] 其中，sina = 〇是平凡解，忽略;则：
[0081]
[0082] 求得sina后，可求得：
[0083]
[0084] 从而构建非限制维纳滤波器G?h(n，k):
[0085]
[0086] 将主麦克风以及次麦克风做短时傅里叶变换后的语音信号与维纳滤波器相乘，并 做离散傅里叶逆变换处理及叠接相加，从而得到降噪之后的时域语音信号。
[0087] 由上述本发明提供的技术方案可以看出，可以根据用户所使用模式的不同来选取 对应的降噪方法;在正常手持通话模式下，基于功率水平差异的方法适用于该模式下的近 场信号模型，且能有效利用主次麦克风接收语音信号在能量上的较大差异信息;在使用扬 声器模式下，利用基于双麦克风之间相关性的方法不需要进行噪声功率谱的估计，且能够 避免两个麦克风接收噪声不一致对于基于功率水平差异方法的影响;此外，相比于传统的 基于功率水平差异的方法，本发明通过构造的维纳滤波能有效避免噪声功率谱估计错误导 致的音乐噪声。
【附图说明】
[0088] 为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用 的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本 领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他 附图。
[0089] 图1为本发明实施例提供的一种应用于双麦克风小型手持设备的语音降噪方法的 流程图；
[0090] 图2为本发明实施例提供的基于双麦克风功率水平差异的语音降噪方法的流程 图；
[0091 ]图3为本发明实施例提供的基于双麦克风之间相干性的语音降噪方法的流程图；
[0092] 图4为本发明实施例提供的基于双麦克风之间相干性的语音降噪技术信号模型。
【具体实施方式】
[0093] 下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整 地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本 发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施 例，都属于本发明的保护范围。
[0094] 如图1所示，为本发明实施例提供的一种应用于双麦克风小型手持设备的语音降 噪方法的流程图，其主要包括如下步骤：
[0095] 步骤11、接收双麦克风的语音信号，并进行通话模式的判别；
[0096] 步骤12、若当前为手持通话模式，则采用基于双麦克风功率水平差异的语音降噪 方法来获取降噪后的语音信号；
[0097] 步骤13、若当前为扬声器通话模式，则采用基于双麦克风之间相干性的语音降噪 方法来获取降噪后的语音信号。
[0098] 本发明实施例中，根据用户所使用模式的不同来选取对应的降噪方法。在正常手 持通话模式下，基于功率水平差异的方法适用于该模式下的近场信号模型，且能有效利用 主次麦克风接收语音信号在能量上的较大差异信息；在使用扬声器模式下，利用基于双麦 克风之间相关性的方法不需要进行噪声功率谱的估计，且能够避免两个麦克风接收噪声不 一致对于基于功率水平差异方法的影响。
[0099] 下面主要针对基于双麦克风功率水平差异的语音降噪方法以及基于双麦克风之 间相干性的语音降噪方法进行详细的介绍。
[0100] 1、基于双麦克风功率水平差异的语音降噪方法。
[0101] 在正常的手持通话模式下，我们默认的是双麦克风与说话人嘴巴构成一个近场环 境。在该环境下的声波模型为球形波，传统的基于相位的降噪技术并不适用这一模型。而利 用功率水平差异(Power Level Difference，PLD)，则避免了这一问题。
[0102 ]其具体的处理流程如图2所示，主要为：
[0103] 主麦克风与次麦克风接收到的语音信号如下：
[0104] XI(m) =hi(m)*s(m)+ni(m)；
[0105] X2(m) =h2(m)*s(m)+n2(m);
[0106] 其中，m为采样点，X1(m)为主麦克风接收到的语音信号，X2(m)为次麦克风接收到的 语音信号;hi(m)为声波传播模型冲击响应，m(m)为噪声，i = l，2;s(m)为目标语音，*为卷积 操作；
[0107] 对主麦克风与次麦克风接受到的语音信号分别做短时傅里叶变换将其转换到时 频域，得到：
[0108] Xi(n,k)=Hi(n,k)S(n,k)+Ni(n,k)；
[0109] X2(n，k)=H2(n，k)S(n，k)+N2(n，k);
[0110] 其中，η和k分别表示时间点和频率点;不失一般性，将上述两个式子改写为：
[0111] Xi(n,k) = Si(n,k)+Ni(n,k)；
[0112] X2(n，k)=Hi2(n，k)Si(n，k)+N2(n，k);
[0113] 其中，SKhk)表示出(11，1〇3(11，1〇，从而有
4艮据短时傅里叶变 换结果分别计算主麦克风与次麦克风带噪音的功率谱密度(PSD)，得到：
[0114]
[0115]
[0116] 其中，& 表示主麦克风接收到的语音信号的PSD，仏幻表示主麦克风接收 到的噪声信号的PSD，f v:(?A)表示次麦克风接收到的噪声信号的PSD。
[0117] 将 ' ("，々: (/7，幻相减得到：
[0118]
^
[0119]令碼("，幻-(仏々)，Δ/\ = Ρ、.（"，Α·)-Λ,」(/α);另外，一般假设两个麦克 风接收到的背景噪声可以认为差异不大，故A ΡΝ可以忽略不计，即Δ ΡΝ~〇，则有：
[0120]
[0121] 利用估计的语音信号PSD和噪声信号PSD可以构造维纳滤波器G^(n，k)，这里我们 利用估计的主麦克风语音信号PSD & 幻和噪声信号PSD^ (?，幻进行构造：
[0122]
[0123] 其中，GAP(n，k)的下标ΔΡ表示该滤波器是基于功率水平差得到的。&(?，幻不能直 接估计得到，故可将I ΔΡχ(η，1〇 I表达式带入上式，得到：
[0124]
[0125] 为了不使语音产生太大失真，在上式中加入一个自由参数α，则有G^(n，k)的估计 式：
[0126]
[0127] 其中，主麦克风噪声的PSD采用语音信号开始的T纯噪声帧计算，公式如下：
[0128]
[0129]式中，λΝ为噪声遗忘因子，由经验值设定;ΧΚη,?Ο表示主麦克风接收信号的时频域 值；
[0130] 则现在G^P(n，k)估计式中仅有冲击响应函数H12(n，k)为未知的；考虑互功率谱密 度(CPSD)，主麦克风与次麦克风接收信号的CPSD为：
[0131]
[0132] 其中，幻为两个麦克风接收噪声信号的CPSD，类似于&(?,幻的计算，可以 通过下式进行估计：
[0133]
[0134] 可以通过主麦克风接收信号PSD与估计噪声PSD的差值进行估计：
[0135]
[0136] 从而得到估计的冲击响应函数//12(〇，幻：
[0137] -1 -.
r. ."1 · /
[0138] 本发明实施例中，两个麦克风接收信号的PSD和CPSD(即Λ ("，幻和&；^0,幻)我们 采用如下的递归平均方法来估计：
[0139]
[0140]
[0141] 其中λχ为带噪语音遗忘因子。此时G^(n，k)表达式中的参数均求出，则获得了用于 增强语音的维纳滤波器，将主麦克风以及次麦克风做短时傅里叶变换后的语音信号与维纳 滤波器相乘，并做离散傅里叶逆变换处理（IFFT)及叠接相加，从而得到降噪之后的时域语 首?目号（5'〗、幻:）。：
[0142] 2、基于双麦克风之间相干性的语音降噪方法。
[0143] 在使用扬声器通话的模式下，可以认为两个麦克风与说话人嘴巴构成一个远场环 境。此时，本发明采用基于双麦克风接收信号相干性的语音降噪技术。这个方法不需要进行 噪声功率谱的估计，且能够避免两个麦克风接收噪声不一致对于基于功率水平差异方法的 影响。
[0144] 在双麦克风条件下，两个麦克风接收到的语音信号之间是有相干性的。这个相干 性可以用来估计带噪语音信号的每个时频点的SNR，从而利用维纳滤波的思想就能达到语 音降噪的目的。
[0145] 其具体的处理流程如图3所示，主要为：
[0146] 主麦克风与次麦克风接收到的语音信号如下：
[0147] xi(m) = si(m)+m(m)，i = 1，2
[0148] 其中，m为采样点，X1(m)为主麦克风接收到的语音信号，X2(m)为次麦克风接收到的 语音信号，m(m)为噪声， Sl(m)为目标语音；
[0149] 进行短时傅里叶变换后得到：
[0150] Xi(n,k) = Si(n,k)+Ni(n,k), i = 1,2；
[0151 ]其中，η和k分别表示时间点和频率点；
[0152] 定义主麦克风与次麦克风接收到的语音信号的相干性函数为：
[0153]
[0154] 式中幻与仏("，幻分别为主麦克风与次麦克风带噪音的功率谱密度PSD， 巧^2(?，幻为主麦克风与次麦克风带噪语音的互功率谱密度;同样可以采用前文所述的递归 平均方法求得。
[0155] 相干性函数与主麦克风与次麦克风局部信噪比SNRjPSNR2有如下关系：
[0156]
[0157] 其中，1^2和1^2分别表示两个麦克风接收的目标语音相干性函数和噪声相干性 函数。如果两个麦克风非常接近(例如2cm)，则SNRi^SNfc成立，但是如果两个麦克风相隔比 较远（例如15cm)，则上述假设不一定成立。但是可以认为在两个麦克风处(或 1-rSNR 1. )是相等的。
[0158] 贝lj
则上式改写为：
[0159] r%-rSiSG + r；.B2(l-(；);
[0160] 如图4所示，当使用扬声器通话模式时，假设目标语音信号源在双麦克的正前方， 选取两麦克风位置的中心作为阵列参考点，则目标语音方向即为〇°，背景噪声信号源等效 为从Θ方向入射。根据阵列信号处理理论，对于两个同源信号山和1] 2(同一个目标语音或噪声 分别被两个麦克风接收的信号），主麦克风与次麦克风接收到的信号之间的相干性可以表 示为：
[0161] Γ?(ι"；(?) = ^Λ(，//〇8?π0ι
[0162] 式中，fs是采样率，d是主麦克风与次麦克风的间距，c为声速；
[0163] 将上面两个表达式相结合，得到：
[0164]
[0165] 其中，i = fs(d/c)，分别取出实部和虚部：
[0166]
[0167]
[0168] 其中，α= ω TSin0，分别对实部和虚部表达式变型，得到：
[0169]
[0170]
[0171] 由于上述两式相等，则有：
[0172] Icosa = (R-l )sina+I ；
[0173] 根据 cos2a+sin2a = 1，得到：
[0174] (I2+(l-R)2)sin2a+2I(R-l)sina = 〇；
[0175] 其根为：
[0176]
[0177] 其中，sina = 〇是平凡解，忽略;则：
[0178]
(I _ /<厂+厂
[0179] 求得sina后，可求得：
[0180]
[0181] 从而构建非限制维纳滤波器G?h(n，k)(下标coh表示该滤波器是基于相干性得到 的）：
[0182]
[0183] 将主麦克风以及次麦克风做短时傅里叶变换后的语音信号与维纳滤波器相乘，并 做离散傅里叶逆变换处理及叠接相加，从而得到降噪之后的时域目标语音信号（U 2 >。：
[0184] 本发明实施例上述方案，相比于传统的手持设备双麦克风语音降噪技术，可以根 据用户所使用模式的不同来选取对应的降噪方法。在正常手持通话模式下，利用基于功率 水平差异的方法能避免近场效应的影响；在使用扬声器模式下，利用基于双麦克风之间相 关性的方法不需要进行噪声功率谱的估计，且能够避免两个麦克风接收噪声不一致对于基 于功率水平差异方法的影响。此外，相比于传统的基于功率水平差异的方法，通过构造的维 纳滤波能有效避免噪声功率谱估计错误导致的音乐噪声。
[0185] 通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例可 以通过软件实现，也可以借助软件加必要的通用硬件平台的方式来实现。基于这样的理解， 上述实施例的技术方案可以以软件产品的形式体现出来，该软件产品可以存储在一个非易 失性存储介质（可以是CD-R0M，U盘，移动硬盘等）中，包括若干指令用以使得一台计算机设 备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法。
[0186] 以上所述，仅为本发明较佳的【具体实施方式】，但本发明的保护范围并不局限于此， 任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可轻易想到的变化或替换， 都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书的保护范 围为准。
【主权项】
1. 一种应用于双麦克风小型手持设备的语音降噪方法，其特征在于，包括： 接收双麦克风的语音信号，并进行通话模式的判别； 若当前为手持通话模式，则采用基于双麦克风功率水平差异的语音降噪方法来获取降 噪后的语音信号； 若当前为扬声器通话模式，则采用基于双麦克风之间相干性的语音降噪方法来获取降 噪后的语音信号。2. 根据权利要求1所述的方法，其特征在于，所述采用基于双麦克风功率水平差异的语 音降噪方法来获取降噪后的语音信号包括： 主麦克风与次麦克风接收到的语音信号如下： xi(m) =hi(m)*s(m)+ni(m); X2(m) =h2(m)*s(m)+n2(m); 其中，m为采样点，xi(m)为主麦克风接收到的语音信号，X2(m)为次麦克风接收到的语音 信号;hi(m)为声波传播模型冲击响应，m(m)为噪声，i = l，2;s(m)为目标语音，*为卷积操 作； 对主麦克风与次麦克风接受到的语音信号分别做短时傅里叶变换，得到： Xi(n,k) =Hi(n,k)S(n,k)+化（n,k); X2(n,k)=出（n,k)S(n,k)+化（n,k); 其中，η和k分别表示时间点和频率点;将上述两个式子改写为： Xi(n,k) =Si(n,k)+化(n,k); X2(n,k) =Hi2(n,k)Si(n,k)+N2(n,k); 其中，Si(n，k)表示出(11，1〇5(11，1〇，从而有：根据短时傅里叶变换结 果分别计算主麦克风与次麦克风带噪音的功率谱密度PSD，得到：利用估计的语音信号PSD和噪声信号PSD可W构造维纳滤波器G&P(n，k);其中利用估计 的主麦克风语音信号PSD与,的，W和噪声信号PSD 进行构造时：其中，GAP(n，k)的下标ΔΡ表示该滤波器是基于功率水平差得到的;将I APx(n，k)|表达 式带入上式，得到：两个麦克风接收信号的？50和〔?50,即与-,的，<^1和&一-, (η，W，采用如下的递归平均方法来 估计：其中λχ为带噪语音遗忘因子。将主麦克风W及次麦克风做短时傅里叶变换后的语音信 号与维纳滤波器相乘，并做离散傅里叶逆变换处理及叠接相加，从而得到降噪之后的时域 语音信号。3.根据权利要求1所述的方法，其特征在于，所述采用基于双麦克风之间相干性的语音 降噪方法来获取降噪后的语音信号包括： 主麦克风与次麦克风接收到的语音信号如下： xi(m) = si(m)+ni(m), 1 = 1,2 其中，m为采样点，xi(m)为主麦克风接收到的语音信号，X2(m)为次麦克风接收到的语音 信号，m(m)为噪声，si(m)为目标语音； 进行短时傅里叶变换后得到： Xi(n,k) =Si(n,k)+Ni(n,k), 1 = 1,2; 其中，η和k分别表示时间点和频率点； 定义主麦克风与次麦克风接收到的语音信号的相干性函数为：式中，A-,(n，A)与分别为主麦克风与次麦克风带噪音的功率谱密度PSD, Λ., (?，A)为主麦克风与次麦克风带噪语音的互功率谱密度； 相干性函数与主麦克风与次麦克风局部信噪比SNRi和SNR2有如下关系：其中，Γ,,,,和:T",，b分别表示两个麦克风接收的目标语音相干性函数和噪声相干性函数； 令，则上式改写为：当使用扬声器通话模式时，假设目标语音信号源在双麦克的正前方，选取两麦克风位 置的中屯、作为阵列参考点，则目标语音方向即为〇°，背景噪声信号源等效为从Θ方向入射， 则主麦克风与次麦克风接收到的信号化和化之间的相干性为：式中，fs是采样率，d是主麦克风与次麦克风的间距，C为声速;Ui和化表示主麦克风与次 麦克风接收到的语音信号或者噪声信号； 将上面两个表达式相结合，得到：其中，s ina = ο是平凡解，忽略;则：将主麦克风W及次麦克风做短时傅里叶变换后的语音信号与维纳滤波器相乘，并做离 散傅里叶逆变换处理及叠接相加，从而得到降噪之后的时域语音信号。
【文档编号】G10L21/0216GK105976826SQ201610286545
【公开日】2016年9月28日
【申请日】2016年4月28日
【发明人】叶中付, 鲍光照, 罗友
【申请人】中国科学技术大学