一种采样率差异估计与校正系统及其处理方法

专利名称：一种采样率差异估计与校正系统及其处理方法
技术领域：
本发明涉及一种多通道信号处理方法，更具体的说，涉及到自适应噪声消除 (ANC)、盲源信号分离(BSS)等多通道信号处理应用中的一种采样率差异估计与校正 系统及其处理方法。
背景技术：
当音频流从一个设备传输到另一个设备，在播放或者录制的过程中，存在一个 问题在两个设备的数模转换模块中，缺乏一个共同的时钟信号。这可能会导致信 号的采样率出现不一致的情况。导致采样率出现不一致现象的原因如下首先、产 生时钟信号的晶振存在一定频率范围内的容忍度，在一些商用的数字信号处理设备 中，这个容忍度可以从几十个PPM (parts per million)到上万个PPM;其次，晶振 的频率可能收到温度的影响；另外，在一些便携设备中，其所需要的频率可能是从 一个更高的频率通过分频得到的。由于诸如以上的这些原因，例如一个标称为8000Hz 采样的信号，其实际采样率可能为8002Hz，或者其它频率。
在一些基本应用中，比如说VOIP，可以采用一些技巧来克服采样率差异的影响， 例如在那些不为人耳感知的语音段插入一些零值样点，或者在这些语音段删除一些 样点，来保证音频信号流的同步。
另外，在一些更为复杂的应用中，比如涉及到多路信号处理的应用中，如盲源 信号分离，自适应噪声消除等应用，仅仅保证信号的同步并不能满足算法的需要。 当音频信号来自多个不同的采样设备时，在经过算法处理之前，必须保证多路信号 的采样率是一致的。这时候必须进行采样率差异的估计和校正。现有的信号分离或
者自适应噪声消除系统，并没有考虑采样率差异的影响。

发明内容
为了解决上述问题，本发明的目的在于提供一种采样率差异估计与校正系统及 其处理方法。其通过一种改进的相位变换算法，以及sinc函数内插的精细方法，用 于估计相关信号间的采样率差异，并利用polyphase滤波器组进行重采样校正。该 方法，解决了 ANC/BSS等多通道信号处理算法在实际应用中存在的由于设备或者
7环境的原因，导致标称为同一采样率的多通道信号出现采样率不一致的问题。
为了实现上述目的，本发明的采样率差异估计与校正系统，包括采样率差异估 计模块和采样率校正模块。所述采样率差异估计模块，通过改进的相位变换算法 (IPHAT)，以及sinc函数内插的精细方法，用于两路相关信号之间的采样率差异估 计；所述采样率校正模块，用于对两路存在采样率差异的相关信号进行调整，利用 抽取和内插滤波器组进行重采样校正，使得它们具有一致的采样率。本发明的采样 率差异估计与校正系统，首先利用采样率差异估计模块，对两相关信号的采样率差 异进行估计，然后利用采样率校正模块，将两路信号的采样率调整为一致。
另外，本发明的采样率差异估计与校正系统的处理方法，包括如下步骤
1) 选取来自同一源的两路信号中的其中一路，将它从标称采样率向目标范围内 的采样率重新采样，这里认为两路信号的采样率与标称采样率的差异在/M以内，即
目标范围为标称采样率的正负/M范围。
在本发明的方法中，两路信号的采样率与标称采样率的差异/M的范围/^^20Hz 的情况都能处理，例如当/^为20Hz时，目标范围为标称采样率的正负20Hz范围。
2) 每次重采样之后，与另一路信号求互相关函数，记录每次互相关函数的最大 值，最大值对应的目标采样率即为信号的真实采样率。
这里，利用改进的相位变换算法IPHAT,计算两段相关信号的互相关函数的最大 值及其对应坐标，IPHAT算法的完整描述如下
A(咖F7T0(")) %2(w) = ,iTCy("》 G,2(w)""w)X;(w)
12、 2;r丄IG一)卩
其中，x(w)，:Kn)是等长的两段相关信号，X,(w),&(w)是对应的频谱，GJw)是
对应的互相关功率谱，C7代表信号的信噪比，/l是与信噪比fT对应的一个比例因子，
它们之间满足如下的关系
C7<C70
-义
—"0
II
义3) 利用IPHAT算法计算出x(n)与y(n)的互相关函数/^(w)以后，由下面的公式
计算x(n)与y(n)之间的整数点采样率差异A/—int: /or △/ = -/M
及砂=xcwr(W); em/ = max((i);
采样率差异为
△/_int = zWe;c 1;
其中，/M是采样率差异估计的范围，"^/7^/e是重采样函数，表示将信号X
的采样率调整/4//倍，；ccwr是互相关运算函数，max是取最大值函数，/"&x表示
最大值对应序列的坐标，A/一int表示计算出的整数点采样率差异值。
4) 在求出整数点的采样率差异以后，如果需要更精确地确定两段信号之间的采 样率差异，进一步采用基于sine函数内插的方法，对IPHAT算法求出的一系列i 12(")
的最大值进行连续时间域的重构，通过找出重构信号的最大值，从而确定更精确的 最大值坐标，来找到分数阶的采样率差异4/v^，计算公式如下
假设上述利用IPHAT算法求出的各个采样率下，两信号相关函数的最大值组成 的序列为^W，则
= Z d(w) sin c(f - w)
— - "o
式中，首先通过sine核函数内插的方法，将最大值序列重构成连续时间信号，
"。为^(n)中最大值对应的坐标，/mfec为连续时间信号中最大值对应的坐标，A/_/rac
表示分数阶的采样率差异。
5) 将步骤3)中计算出的整数点的采样率差异A/Jnt与步骤4)中计算出的分 数阶采样率差异A/—》。e相加，即得到两路信号之间的更为精确的采样率差异值A/s:int + △/—加c o
6)在计算出两路信号之间的采样率差异A,以后，以乂为基准，对其中一路信 号进行重采样校正，该重采样步骤是利用抽取和内插滤波器实现，最终得到经校正 后的信号。
在采样率校正模块中，首先，根据估计出来的采样率差异值4/;，计算出内插滤 波器和抽取滤波器中内插因子Z和抽取因子M的大小。
其中，内插和抽取因子的计算与采样率差异值之间满足如下的关系式
￡/M = (/s+4/;)//s 其中，Z、 M对应最小的整数，例如，(乂+ 为1.2时，L=6， M=5。
接着，在计算出抽取因子L和内插因子M以后，首先将所选的其中一路信号经 过一个内插滤波器，实现一个向上L倍的内插，然后将内插后的信号通过低通滤波 器，该低通滤波器的作用是避免出现频率混叠，最后将低通滤波器输出的信号通过 抽取滤波器，实现向下M倍的抽取，从而完成将所选的某一路信号的采样率变化 L/M倍，得到采样率转换之后的信号。
本发明的采样率差异估计与校正系统及其处理方法的有益效果在于通过一种 改进的相位变换算法以及sine函数内插的精细方法，估计相关信号间的采样率差异， 并利用polyphase滤波器组进行重采样校正，能够解决ANC/BSS等多通道信号处理 算法在实际应用中存在的由于设备或者环境的原因导致标称为同一采样率的多通道 信号出现采样率不一致的问题。


图1是本发明的采样率差异估计与校正系统的构成示意图。
图2是表示采样率存在差异时的自适应噪声消除算法的性能表现，图中，从上
至下依次为原始噪声、采样率差异为OHz、 lHz、 2Hz、 5Hz时的自适应噪声消除输
出结果。
图3是本发明的采样率差异估计与校正系统的处理方法的流程框图。 图4是本发明的采样率差异估计与校正系统的处理方法中的采样率校正部分的 详细流程图。
图5是利用本发明的方法对两段采样率差异为1Hz的信号进行采样率差异估计 的结果示意图，其中，图5(a)表示两段相关信号， 一段的信噪比为一10dB;图5(b) 是采样率差异估计结果图。
10
具体实施例方式
下面结合附图和具体实施例对本发明的采样率差异估计与校正系统及其处理方 法进行详细的说明。
图1是本发明的采样率差异估计与校正系统的构成示意图。如图1所示，在两 路具有采样率差异的BSS/ANC系统中，两路信号具有标称的同一采样率乂，而实际 上，两者的真实采样率却不是,，在这种情况下，进行多通道信号处理的时候，算 法的性能就会受到影响。
由此，如图1所示，本发明的采样率差异估计与校正系统，包括采样率差异 估计模块和采样率校正模块。采样率差异估计模块完成两路相关信号之间的采样率 差异估计；采样率校正模块完成对其中一路信号的采样率调整，使得它们具有一致 的采样率。
图2是表示采样率存在差异时的自适应噪声消除算法的性能表现，其中，从上 至下依次为原始噪声、以及采样率差异为OHz、 lHz、 2Hz、 5Hz时的自适应噪声消 除输出结果，从实验结果来看，当存在采样率差异时，自适应算法的性能明显受到 影响，表现为出现较大的残余噪声，而且，采样率差异越大，算法表现越糟糕。
在采样率差异估计中，本发明基于的原理是两路相干信号，当它们的采样率 一致时，具有最大的相关性。为此，本发明采取下述策略选取两路信号中的其中
一路，将它从标称的采样率向某个范围的采样率重新釆样，每次重采样之后，与另 一路信号求互相关函数，记录每次互相关函数的最大值。在这些最大值构成的序列 中，最大值对应的目标采样率即为信号的真实采样率。
图3是本发明的采样率差异估计与校正系统的处理方法的流程框图。图3中， ;c(w)和:K")是两相关信号，由两个不同的设备录制同一音源获得，其标称的的采样
率均为/"而实际的采样率可能并不一致，而且，jc(")和;K")中还可能混入了干扰
噪声，表达式如下式所示。
;c(w)-》(w) + ^(M) _y(w) = s2(") + "2(")
其中，A(")，^(")是同一音源到达不同设备被录制的信号，巧(")，"2(")是到达录 音设备的干扰信号。
由于噪声的影响，直接计算x(n)和y(n)时域的互相关函数，其最大值往往不能准确反映《(W)， ^(")的相关性，尤其在噪声巧(n), "2(")也具有较强的相关性时，时域相 关函数的最大值可能反应的是噪声之间的相关性。
由此，如图3所示，本发明的采样率差异估计与校正系统的处理方法，包含采
样率差异估计过程和采样率校正过程，具体包括如下步骤
1) 选取来自同一源的两路信号中的其中一路，将它从标称采样率向目标范围内
的采样率重新采样，两路信号的采样率与标称采样率的差异/w的范围为/M《20Hz。
例如，认为两路信号的采样率与标称采样率的差异/M在20Hz以内，即目标范 围为标称采样率的正负20Hz范围。
2) 每次重采样之后，与另一路信号求互相关函数，记录每次互相关函数的最大 值，最大值对应的目标采样率即为信号的真实采样率。
这里，利用改进的相位变换算法IPHAT，计算两段相关信号的互相关函数的最大 值及其对应坐标，IPHAT算法的完整描述如下
x,(咖FFr(x("))
X2(w) = /^TO;(")) G,2(w)-X,(w)X;0)
物二了^LlM^ w
其中，x(w)，:KW是等长的两段相关信号，A(w),^(MO是对应的频谱，C^(w)是
对应的互相关功率谱，C7代表信号的信噪比，义是与信噪比O"对应的一个比例因子，
它们之间满足如下的关系
-(CT-CT》+ ^ CT。〈0"〈0",
3)利用IPHAT算法计算出x(n)与y(n)的互相关函数i^(")以后，由下面的公式 计算x(n)与y(n)之间的整数点采样率差异A/一itit 。
12/or A/ = -/M
A 二画msww_p/e(jc，/,州；
及砂=xcwr(x,，力；
"(A/ + /M+l) = max("; 二 max④；
整数点采样率差异为
A/一int -画zWex - /M —1;
其中，/M是采样率差异估计的范围，"sfl附p/e是重采样函数，表示将信号x 的采样率由乂调整到/， ；ccwr是互相关运算函数，max是取最大值函数，/"A;c表 示最大值对应序列的坐标，A/一int表示计算出的整数点采样率差异值。
这里，假设两路信号的采样率与标称采样率的差异/A/在20Hz以内，即目标范 围为标称采样率的正负20Hz范围。则x(n)与y(n)之间的整数点采样率差异A/jnt为
，A/ = -20:1:20 /
A = resam—(x，/, ;
牟/ + 20 + 1) = max(&); em/ = max(d); 」nt 二 -20—1;
4)在求出整数点的采样率差异以后，如果需要更精确地确定两段信号之间的采 样率差异，进一步采用基于sine函数内插的方法，对IPHAT算法求出的一系列及12(") 的最大值进行连续时间域的重构，通过找出重构信号的最大值，从而确定更精确的 最大值坐标，来找到分数阶的采样率差异，计算公式如下
假设上述利用IPHAT算法求出的各个采样率下，两信号相关函数的最大值组成 的序列为""入则d(0 = Z (w) sin c(卜w) = max(d(O) = 一 "o
式中，首先通过sinc核函数内插的方法，将最大值序列d(w)重构成连续时间信号， 为d(")中最大值对应的坐标，/mfec为连续时间信号中最大值对应的坐标， 表示分数阶的采样率差异。
5) 将上述步骤3)中计算出的整数点的采样率差异与步骤4)中计算出的分数
阶采样率差异相加，即得到两路信号之间的更为精确的采样率差异值4/;:
A/s =A/—int + A/—>ac 。
6) 在计算出两路信号之间的采样率差异以后，以fs为基准，对其中一路信号进 行重采样校正，该重采样步骤是利用抽取和内插滤波器实现，最终得到经校正后的 信号。
在采样率校正模块中，首先根据估计出来的采样率差异值4/;，计算出内插滤波
器和抽取滤波器中内插因子Z和抽取因子M的大小。其中，内插和抽取因子的计算
与采样率差异值之间满足如下的关系式
其中，Z、 M对应最小的整数，例如，(X+A/;)/X为1.2时，L=6， M=5。
接着，在计算出抽取因子丄和内插因子M以后，首先将所选的其中一路信号经 过一个内插滤波器完成信号向上L倍的内插，然后将内插后的信号通过低通滤波器， 低通滤波器的作用是避免出现频率混叠，最后将低通滤波器输出的信号通过抽取滤 波器，完成向下M倍的抽取，从而完成将所选的某一路信号的采样率变化L/M倍， 得到采样率转换之后的信号。
图4是本发明的釆样率差异估计与校正系统的处理方法中的采样率校正部分的 详细流程图。
如图4所示，抽取和内插滤波器合二为一的采样率校正系统图，完成将信号x(n) 的采样率变化L/M倍，首先，将信号经过一个内插滤波器，实现一个向上L倍的内 插，然后，将内插后的信号通过一个低通滤波器以防止频率混叠，最后，将低通滤 波器的输出信号通过一个抽取滤波器，实现向下M倍的抽取。其中，抽取和内插滤 波器的实现采用了多相分解的方法。性能分析
在实际应用中，两路信号的采样率可能都不是标称的采样率々，这里假设两个
采样率分别为/,/2，两者之间的采样率差异为4/;=/2-/，以々为基准，根据估 计出的采样率差异(这里默认估计准确)，把信号x(n)进行采样率转换，转换后的采
样率为
,=乂*鹿=/;*(,+4/;)//,
此时，重采样之后的x(n)与y(n)之间的采样率差异A/为
△/=/2—,'=/2—/*(/s+A/;)a =(/2—y;)/x 《/广y;。
由此可见，经过本发明的估计与校正系统的处理，最终所获得的采样率差异远
小于两路信号原始的采样率差异，从而能够解决ANC/BSS等多通道信号处理算法在
实际应用中存在的由于设备或者环境的原因导致标称为同一采样率的多通道信号出 现采样率不一致的问题。
性能试验
为考察本发明方法的性能，通过下述考量方法进行了试验
图5是利用本发明的方法对两段采样率差异为1Hz的信号进行采样率差异估计 的结果示意图，其中，图5 (a)表示用于采样率差异估计的两段相关信号， 一段的信 噪比为-10dB;图5(b)是采样率差异估计结果图。
在没有cross-talk的情况下，定义一个变量ENR (error to noise ratio)
￡A^ = 10*logl0(^^) 五A^值反应了残余噪声与初始噪声的能量之比，上式中，e代表残余噪声，n代
表原始噪声。ENR值越小，自适应噪声消除算法性能越好。
下表中给出针对语音信号，采样率从0Hz变为10Hz时，ENR的变化。
0Hz1Hz2Hz3Hz4Hz5Hz6Hz7Hz8Hz9Hz
E廳(dB)-37.28-16.08-11.18-8.11-6.18-3.39-2.57-2.29-1.91-1.47
ENR2(dB)-37.28-37.34-37.09-37.26-37.26-37.16-37.03-37.14-37.03-37.15
15其中，ENR1代表未进行采样率差异调整时的值，ENR2代表进行采样率差异调 整以后的值。
在实际环境中，带有cross-talk的情况下，考査ANC算法的性能，两路采用两 种设备录制， 一个为电脑声卡， 一个为爱国者录音笔，如图5 (b)所示，结果如下: 校正之前ENR = -12.516dB
校正之后，ENR = -16.338dB
由此可见，采样率差异问题确实存在于自适应噪声消除或者盲源分离的应用中， 并且严重影响了自适应算法的性能，而本发明所提出的相关信号间采样率差异估计 和校正方法，则较好地解决了这一问题。
权利要求
1、一种采样率差异估计与校正系统，其特征在于，包括采样率差异估计模块和采样率校正模块，所述采样率差异估计模块，用于通过改进的相位变换算法IPHAT以及sinc函数内插的精细方法估计两路相关信号之间的采样率差异；所述采样率校正模块，用于对两路存在采样率差异的相关信号进行调整，利用抽取和内插滤波器组进行重采样校正，使得它们具有一致的采样率，该采样率差异估计与校正系统，首先利用所述采样率差异估计模块，对两相关信号的采样率差异进行估计，然后利用所述采样率校正模块，将两路信号的采样率调整为一致。
2、 如权利要求1所述的采样率差异估计与校正系统，其特征在于，所述采样率 差异估计模块，是基于相关性原理，计算信号间互相关函数时，采取一种改进的相 位变换算法IPHAT， IPHAT算法的具体过程如下Z2(w) = FiTCK")) G,2 0)-《0)X;0)'八 2;r JjG12 其中，x("),:K")是等长的两段相关信号，《(w)，^^(nO是对应的频谱，G"w)是对应的互相关功率谱，C7代表信号的信噪比，；t是与信噪比^T对应的一个比例因子， 它们之间满足如下的关系<formula>formula see original document page 2</formula>
3、如权利要求2所述的采样率差异估计与校正系统，其特征在于，所述采样率 差异估计模块，在利用IPHAT算法计算出x(n)与y(n)的互相关函数/^(")以后，由下面的公式计算x(n)与y(n)之间的整数点采样率差异A/_int:<formula>formula see original document page 3</formula>及砂=xcwrO',力；<formula>formula see original document page 3</formula>其中，/M是采样率差异估计的范围，m^wp/e是重采样函数，表示将信号x的采样率调整/"/倍，；ccwr是互相关运算函数，max是取最大值函数，/"&x表示最大值对应序列的坐标，A/—int表示计算出的整数点采样率差异值；在求出整数点的采样率差异之后，采用基于sine函数内插的方法进一步确定两 段信号之间的采样率差异，对IPHAT算法求出的一系列及12(")的最大值进行连续时 间域的重构，通过找出重构信号的最大值，从而确定更精确的最大值坐标，来找到 分数阶的采样率差异，计算公式如下假设上述利用IPHAT算法求出的两信号相关函数的最大值组成的序列为W"人 求分数阶采样率差异A/—/rae的公式表述如下，<formula>formula see original document page 3</formula>式中，首先通过sinc核函数内插的方法，将最大值序列d(")重构成连续时间信号，"。为d(w)中最大值对应的坐标，/mtec为连续时间信号中最大值对应的坐标，表示分数阶的采样率差异；然后，将上述计算出的整数点的采样率差异与分数阶采样率差异相加，获取两 路信号之间的更为精确的采样率差异值A/s: 4^=A/_int + A/一加c 0
4、 一种采样率差异估计与校正系统的处理方法，其特征在于，包括如下步骤1) 选取来自同一源的两路信号中的其中一路，将它从标称采样率向目标范围内 的采样率重新采样，这里认为两路信号的采样率与标称采样率的差异在/m以内；2) 每次重采样之后，与另一路信号求互相关函数，记录每次互相关函数的最大 值，在这些最大值构成的峰值序列中，最大值对应的目标采样率即为信号的真实采 样率，这里，利用改进的相位变换算法IPHAT，计算两段相关信号的互相关函数的最大 值及其对应坐标，IPHAT算法的完整描述如下<formula>formula see original document page 0</formula>12、 2;r丄IG,)卩 其中，x(")，K")是等长的两段相关信号，《(w),X2(w)是对应的频谱，( 12 )是对应的互相关功率谱，ff代表信号的信噪比，A是与信噪比cr对应的一个比例因子，它们之间满足如下的关系<formula>formula see original document page 0</formula>3)利用IPHAT算法计算出x(n)与y(n)的互相关函数&(")以后，由下面的公式 计算x(n)与y(n)之间的整数点采样率差异A/Jnt:= ms"w/ /e(;c,/， ; "(A/ + /M+l) = max(及砂)；end[max—va/we, zWex] = max(t/); V—int = /w^x - j^ -l;其中，/M是采样率差异估计的范围，"^附;^是重采样函数，表示将信号x 的采样率调整/"/倍，jccwr是互相关运算函数，max是取最大值函数，/mfec表示最大值对应序列的坐标，A/—int表示计算出的整数点采样率差异值；4) 在求出整数点的采样率差异之后，采用基于sinc函数内插的方法进一步确定 两段信号之间的采样率差异，对IPHAT算法求出的一系列及12(")的最大值进行连续 时间域的重构，通过找出重构信号的最大值，从而确定更精确的最大值坐标，来找 到分数阶的采样率差异，计算公式如下假设上述利用IPHAT算法求出的两信号相关函数的最大值组成的序列为W^， 求分数阶采样率差异的公式表述如下，<formula>formula see original document page 5</formula>式中，首先通过sinc核函数内插的方法，将最大值序列d(")重构成连续时间信 号， 为c/(")中最大值对应的坐标，/mfec为连续时间信号中最大值对应的坐标，A/l/^表示分数阶的采样率差异；5) 将步骤3)中计算出的整数点的采样率差异与步骤4)中计算出的分数阶采样率差异相加，即得到两路信号之间的更为精确的采样率差异值4/;:6) 在计算出两路信号之间的采样率差异A/;以后，以々为基准，对其中一路信号进行重采样校正，该重采样利用抽取和内插滤波器实现，最终得到经校正后的信 号。
5、 如权利要求4所示的采样率差异估计与校正系统的处理方法，其特征在于， 所述两路信号的采样率与标称采样率的差异/M^20Hz 。
6、 如权利要求4所示的采样率差异估计与校正系统的处理方法，其特征在于，所述步骤6)中，首先，根据估计出来的采样率差异值4/;,计算出内插滤波器和抽取滤波器中内插因子L和抽取因子M的大小，其中，内插因子i和抽取因子M的计 算与采样率差异值之间满足如下的关系式其中，￡、 M对应最小的整数；接着，在计算出抽取因子Z和抽取因子M以后，将所选的其中一路信号经过一 个内插滤波器，实现向上丄倍的内插，然后，将内插后的信号通过一个低通滤波器 以防止频率混叠，最后，将低通滤波器输出的信号通过一个抽取滤波器，实现向下M倍的抽取，从而完成将所选的某一路信号的采样率变化L/M倍，得到采样率转换 之后的信号。
7、 如权利要求6所示的采样率差异估计与校正系统的处理方法，其特征在于， 所述抽取和内插滤波器的实现采用多相分解的方法。
8、 如权利要求6所示的采样率差异估计与校正系统的处理方法，其特征在于， 所述用于防止频率混叠的低通滤波器采用有限冲激响应低通滤波器。
全文摘要
本发明提供一种采样率差异估计与校正系统及其处理方法。该系统包括采样率差异估计模块和采样率校正模块。首先通过采样率差异估计模块，利用改进的相位变换算法IPHAT，求两路信号的互相关函数；然后，计算它们之间的整数点采样率差异；接着进一步采用基于sinc函数内插的方法，获得分数阶的采样率差异；将整数点的采样率差异与分数阶采样率差异相加，即得到两路信号之间的更为精确的采样率差异值；在计算出两路信号之间的采样率差异以后，通过采样率校正模块，利用抽取和内插滤波器对其中一路信号进行重采样校正。该方法可解决ANC/BSS等多通道信号处理算法在实际应用中存在的由于设备或者环境的原因导致标称为同一采样率的多通道信号出现采样率不一致的问题。
文档编号G10L21/00GK101645273SQ20091008873
公开日2010年2月10日 申请日期2009年7月10日 优先权日2009年7月10日
发明者强 付, 波 覃, 颜永红 申请人:中国科学院声学研究所;北京中科信利技术有限公司