麦克风算法仿真及实验台搭建

1、麦克风算法仿真及平台搭建基础知识1.1阵列信号系统模型123M-1 M图1.1：麦克风等间距线阵示意图我们考虑等间距摆放的麦克风线性阵列，其示意图见图1.1。考 虑到语音信号为宽带信号，我们在窄带信号模型的基础上，引出宽带 信号模型。窄带信号模型在信号处理中，窄带信号与宽带信号是相对而言的。设信号s(t)的 带宽为B，中心频率为f。，若信号s(t)满足相对带宽条件，即： TOC o 1-5 h z 匕口 1，一般取匕 10(1)BB则信号s(t)为窄带信号，反之为宽带信号。其中寿定义为J+ w | S(w) |2 dw(2)0I S (w) |2 dwB定义为i+ (w - f )2 I S

2、(w) |2 dw/c、B =一0 (3)VJ+81 S (w) |2 dw-s式中S(w)为s(t)的频谱。中心频率为f的窄带信号s(t)，其实数形式的表达式为：0s(t) = u(t )cos2 丸 ft+(t)(4)0为了推导方便写成复数形式，即s(t) = u(t )ej2 兀 f0+中(t)(5)式中，u (t)为基带信号，其最高频率远小于正弦信号的频率f0， cos2丸ft+顿t)为载波信号，其中中(t)是正弦信号的初相位，假定为一 0个常数。为了简化模型，通常我们做如下假设：麦克风阵列位于信号源的远场，则接收到的信号为平面波；传输介质是无损的、线性的、非扩散性的、均匀的且各向同

3、性；麦克风的孔径远小于入射波的波长，而且麦克风无指向性， 可近似认为麦克风是点元；麦克风阵列的间距远大于麦克风孔径，各麦克风之间的相互 影响可以忽略不计。现在考虑当信号源为单个窄带信号时，对于间距为d的均匀线阵 (如图1.1)，窄带信号s (t)的入射角度为0，麦克风数为M个，声速 为c = 340m/s。由于信号到达各个麦克风的时间不同，第m个麦克风 接收到的信号为(6)x (t) = s(t -t ) + n (t) m = 0,1,., M -1侦 生羿m为信号到达第m个麦克风相对其到达参考麦克风的时间 延迟，n (t)为加性噪声。由式(6 )可知，信号到达第m个麦克风相对于到达参考麦克

4、风只发生了一个相移ej矶0 则式(6 )可变为:七(t) = e -j 2 兀 f0T m(0)s (t) + 七(t)将麦克风接收到的信号写成M x 1维的向量形式，则麦克风输出为：x (t)e -j2 兀 fT 0(0)0: -=s (t):.x(t)L M-1_ e j fT m -1(0)_x(t)=+n(t)(8)设a(f ,0 ) = e-j2兀)e-j2兀醇M 10小为方向矢量，则 0 x(t) = a(f ,0)S(t) + n(t) 0对于D个窄带入射信号，设入射角为0d,d =1,2,D，则方向矢量为如f0,0 )=回 f 吓 a( f 02)，.,a(九乩) 则麦克风输出

5、可表示为：(10)x(t) = A( f，0 )s(t) + n(t) 0其 中，s(t) = *(t),t),., SD(t)Tx(t X.t)x,t ( ) x.T. 为麦)克(11)D个窄带入射信号。接收到的信号。n(t) = n0(t),n1(t),., n (t)t 为加性噪声向量。宽带信号模型根据阵列信号处理理论，如果用窄带波束形成方法来处理宽带语音信号会导致信号畸变，且信号的带宽越大，畸变越严重。因此，需要用宽带波束形成方法来进行语音增强。假定宽带信号W)是由/个窄带信号叠加而成，即s(t)二乙(t)(12)fm m=1式中七(t)是以fm为中心频率的窄带信号。假设收到的为单个宽

6、带信 源，对于由M个麦克风组成的麦克风阵列，其输出为(13)x(t)室 a( f , 0 )s (t) + n(t)jJj式中a(f ,0) = 32f (0) 32f M1(0勺，其中T (0)为宽带信号到达第m个 jm麦克风相对其到达参考麦克风的时间延迟。当有D(M D 1)个远场信号，分别来自0 ,0 ,.,0方向，则第m麦克风阵列输出为x (t) = 2s (t -t (0 ) + n (t)(14)mi m imi=1式中，t (0 )是第i个信号到达第m个麦克风相对其到达参考麦克风的 m i时间延迟，N (t)为加性噪声。对式(14)进行傅里叶变换，得到(15)x (f ) = %

7、 e-J2兀 fj (d/c)sin(0j )(m-1)s (f ) + n (f ) m ji j m ji=1s( f ) = s (f ), s (f ),., sj 1 j 2 j D(f )r表示D个远场宽带信号的傅里叶变换，n(f ) = N (f ),N (f ),.,N(f )t表示加性噪声信号的傅里叶变换，j0 j 1 jM -1 jA(f ) = a( f ,0 ),a( f ,0 ),., a( f ,0 )表示D个远场信号的方向矢量。则麦 jj 1 j 2j D(16)克风阵列输出宽带频域模型的矩阵形式表示为:x( f ) = A(f )s( f ) + n( f )

8、j j j j在麦克风阵列语音增强中，噪声来源取决于实际的应用环境，噪 声场是一个很重要的概念，不同的方法可能适用于不同类型的噪声场。语音信号在封闭环境传输时，由于反射、衍射的影响，会使麦克 风阵列接收的语音来自多个路径，这种现象称为混响。混响对语音的 影响表现为两个方面：首先，混响使语音信号产生失真；其次，混响 改变了噪声场的空间特性。混响对噪声场空间特性的改变体现为，单一的噪声源被扩展为多 个不同方位的镜像源，噪声空间接近于散射噪声场。这一特性借助模 平方相干函数来界定。设位于,和尸的两麦克风接收到的噪声分别为 和n，相应的自功率谱及互功率谱分别(r, w)、Q(r , w)和中(p ,

9、w)。1 i1ia于是，定义噪声的空时域模平方相干函数为C (w) =1 r (w)|2= Q(piw)|2(17)司 司(r, w)(r , w)根据C (w)的值把噪声场分为相干噪声场、非相干噪声场和散射噪 声场。相干噪声场相干噪声场是指噪声信号直接从噪声源传播到麦克风阵列，而没 有经过任何因传播环境而产生的反射、散射和耗散。在相干噪声场中， 阵列中的不同麦克风接收到的噪声信号具有强相干性，即ir (w)i2归1。ij在实际环境中，相干噪声场存在于开放的空气环境中，在这种环境的 周围没有主要的障碍物影响声音传播，而热湍流的影响也很小。非相干噪声场非相干噪声场是指在给定的空间位置的噪声信号与

10、其它所有位置的噪声信号都不相关，也就是说ir(讪也.0。这样理想的非相干噪 万声场在实际环境中不存在。而在麦克风阵列应用中，麦克风阵列中的 电学噪声通常被认为是一种非相干噪声。非相干噪声也被叫做空间白 噪声。散射噪声场在散射噪声场，相同能量的噪声信号在各个方向同时进行传播。 因此，散射噪声场中的麦克风接收到的噪声信号是弱相关的，并且具 有近似相同的能量。许多实际的噪声环境都可以被看作是散射噪声场。 在散射噪声场中，任何两点的噪声信号之间的相干性可以用麦克风之 间的距离函数来表示(f)|= sin c(色K)(18)式(18)中d是麦克风i和j之间的距离。，j语音增强算法分析及仿真语音增强技术是

11、信号处理领域中的主要研究内容，在现代多媒 体通信体技术以及人工智能系统中应用广泛。语音增强方法自提出至 今，经历了单麦克风语音增强技术到阵列语音增强技术的转化，因为 麦克风阵列语音处理技术采用“电子瞄准”的方法，可以从声源位置 获取较高品质的语音因而得到广泛应用。其中波束形成是麦克风阵列 语音增强的典型算法，它主要分为三类4：固定波束形成、后置滤波波束形成和自适应波束形成。固定波束形成语音增强由于采用的滤波器的权值固定不变，因而 方法简单，运算量低，容易实现，能够抑制背景噪声，但是不具备自 适应抑制强干扰的能力。因此对于可变噪声场效果一般。后置滤波语音增强方法可以用来去除非相干噪声，但在相干噪

12、声 场情况下性能退化，甚至不能应用。自适应波束形成方法，例如Frost提出的线性约束自适应波束形 成法5,是在某种约束条件下，使输出信号能量达到最小。这是用优 化手段实现语音增强的一种方法，这种算法在相干噪声场情况下，可 以得到较高的信噪比改善，但在弱相干噪声场或在散射噪声场中，性 能甚至不如传统的固定波束形成方法。2.1固定波束形成固定波束形成方法是麦克风阵列语音增强算法中最简单易于实 现的一种方法，也是最经典的麦克风阵列语音增强的算法之一。固定 波束形成的含义是，是其波束形成器的权值固定不变，与麦克风阵列 的接收信号无关。固定波束形成算法通过延时控制来补偿从声源到每 个麦克风的延时，对每个

13、麦克风接收到的信号进行延迟补偿，然后使 麦克风阵列波束指向有最大输出功率的方向，即波束对准相应空间位 置的声源信号，其原理如图2.1所示。时延补偿x ；(n) ox ；( n)Q.X M 图2.1固定波束形成原理图该固定波束形成器是典型的延迟-求和波束形成器，它包含三个 部分：时延估计部分，时延补偿部分和加权求和部分。通过麦克风接 收到的数据来估计各个麦克风的时延，然后通过时延补偿使麦克风信 号同步，再加权相加得到波束形成器的输出。设麦克风接收到的信号为x:.，权系数为巧，则延迟-求和波束 形成输出为：y(t) = M 1 w(i)x: (t -t ,)(19)i=0其中M为麦克风的数日，t，

14、表示由时延估计得到的时延。若转换 成数字形式时延为T = d -Sin (i-1)f，其中f为采样频率。i css该方法适合用来消除非相干噪声或者弱相干噪声，对于相干噪声 几乎没有消噪能力。如果用来消除非相干噪声或弱相干噪声，它的消 噪能力主要取决于麦克风的数日，麦克风的数日越多，去噪效果越明 显。因此，该方法获得去噪效果是以增加麦克风数日为代价的，在实际中较少单独使用。2.2后置滤波波束形成后置滤波可以用来去除声学环境中的非相干噪声，人们提出将波 束形成器的输出信号通过后置滤波器以进一步提高输出信号信噪比。 1977年Allen利用波束形成的方向性增益和维纳滤波器对非相干噪声 抑制能力，提出

15、将后置滤波用于语音增强。1988年Zelinski将其扩展 到麦克风阵列中，对延时一求和波束形成的输出端使用自适应维纳滤 波器，从而提出基于后置滤波的麦克风阵列语音增强方法，进一步改 善了噪声抑制能力。鉴于实际声学环境下噪声类型和噪声场的复杂性， 许多学者将性能互补的方法相结合，以改善语音增强方法的噪声抑制 能力。1996年，Fischer和Simmer提出了 GSC与 Wiener滤波相结合 的频域麦克风阵列语音增强方法，以抑制声学环境中的相干和非相 干噪声。2003年Gannot和Cohen提出了基于传递函数GSC和后置 滤波的频域麦克风阵列语音增强方法口爬,以去除噪声场中可能存在 的相干

16、噪声、非相干噪声及非平稳噪声。下面具体介绍一下基本后置滤波算法的原理。后置维纳滤波是后置滤波和波束形成技术相结合的产物。滤波求 和波束形成器的输出由维纳滤波器进行后置滤波处理，其目的是进一 步减少己经由固定波束形成器去除不相干的噪声后的残余噪声成分。 经典的后置滤波方法是R.Zelinski在1988年提出，其系统结构如图 2.2所示。延时补偿图2.2后置自适应滤波系统框图其中，固定波束形成器累加之前每个通道的输出为：j (n) = w (n)x (n)(20)因此，固定波束形成器的输出为j(n) = M 1 w (n)x (n)(2Di=0其中M为麦克风的数目。信号经过固定波束形成器后，进行

17、维纳滤波的后置处理。维纳滤波方法就是按照最小均方误差的准则对期望信号s(,)进行估计，即选取爪)的估计值血)，使均方误差如式(22)最小。 = Es(t) - s(t )2 = Ee2(t)(22)也即要设计一个数字滤波器h(t)，当输入为x(t)时(x(t) = s(t) + n(t)，(23)滤波器输出：s(t) = x(t) * h(t)满足式(22)。根据正交性原理，最佳仰)必须满足对所有k有下式成立：Es(t) - s(t )(t- k) = 0(24)将式(23)代入式(24)，并对式子两边取傅里叶变换，得到我们常说的维纳-霍夫方程H (o)=R (o)R*(o)(25)式中R(o

18、 )为X(t)的功率谱密度，R (o )为s(t)与x(t)的互功率谱密度。由于s(t)和n(t)不相关，即R (o) = 0，则可得snR (o) = R (o)(26)R (o) = R (o) + R (o)(27)于是，式(25)变为H (o)=R (o)R (o j+ R (o)(28)式(28)即为维纳滤波器。采用维纳滤波最大的好处是增强后的残留噪声类似于白噪声，而 不是有节奏起伏的音乐噪声。但是维纳滤波只是在平稳条件下才能保 证在最小均方误差意义下的最优估计。而语音是非平稳的，只在短时 间内近似平稳，实际环境中的噪声也是非平稳的，这是维纳滤波的缺 陷所在。由于维纳滤波器可以实现理

19、论上的最优滤波，因此，如果把波束 形成方法与维纳滤波相结合，可以进一步减少固定波束形成器去除不 相干噪声后残留的噪声成分，使得处理的效果得到进一步的加强。图2.2为具有后置滤波的语音增强系统的结构图。该系统首先对各个麦克风接收到的信号进行延时估计和补偿，并进行了 FFT变换； 随后计算接收信号的自功率谱和互功率谱，由这些短时功率谱可以得 到自适应维纳滤波系数；最后将相加平均后的信号输入到维纳滤波器 中，从而得到所求语音信号的估计值。这种算法要求假定麦克风阵列 的间距比较大，同时要求所求语音信号源于麦克风的距离不小于各噪 声源与麦克风的距离，这时可认为语音和噪声，各通道噪声之间是彼 此不相干的。

20、每个通道的信号经过延时估计和补偿后的信号为X () = s() + n ()(29)即该信号包括所期望的有用信号s(s)和噪声信号n (w)。h(.)为后置滤波的传递函数，其最优解可以表示为件)=5(30)ssnn其中，七(w)是所期望的有用信号的自功率谱密度，虬(w)是噪声信号的自功率谱密度。对于随机过程来说，由于功率谱和相关函数互为傅里叶变换对， 于是通过计算，我们有第，和第;声道上信号的自功率谱和互功率谱密度分别如下：(31)(32)(w) =e(w) +e (w) +e(w) +e(w)e(w) =e(w) +e(w) +e (w) +e (w)XiXjSSjSSljniS假定：信号与

21、噪声是不相关的(即e (w) = 0,Vi)；每个麦克风上的噪声信号是不相关的(e (w) = 0,Vi* j)；n ni j所有麦克风上的噪声信号有着相同的功率谱(。(o)=。(),Vi)。i i那么式子(13 )和式(14)可简化为：(33)(34)。(O)二。()+。(o) TOC o 1-5 h z 矿(o) =8(o)罕jss将式(33)、(34)代入式(30)中，则可以得到后置滤波的传递函数。通过对所有麦克风组合的谱密度估计进行平均，就可以得到更加准确 的估计值:(35)点 U Re(8 (o) )=NN - D i=i J=i+1立、1 为(o)Nxxi=1其中Reg表示取实部的

22、运算符，由于8 (o)实际上代表期望信号的自i j功率谱，所以它应该是一个实数值。因此取8 (o)的实部，以减少估 i j计误差。另外对于随机过程来说，由于功率谱和相关函数互为傅里叶变换对，所以有8 (o) = X (o)X*(o)， 8 (o) = X (o)X*(o) =| X (o)|2。 xxijx.x.iiii j后置维纳滤波算法需要假定阵列中麦克风间的距离比较大，同时 要求所求语音信号源与麦克风的距离小于各噪声源与麦克风的距离，这时可以认为语音和噪声、各通道噪声2间彼此是不相干的。后置维 纳滤波算法的性能与噪声的相干场有关，对于相干噪声场，后置维纳 滤波器儿乎没有消噪能力，而对于非

23、相干场则效果明显。2.3自适应波束形成：线性约束最小方差波束形成算法线性约束最小方差(linearly constrained minimum variance, LCMV)波束形成算法可以实时地对所求信号进行处理，同时又可有效抑制噪 声干扰。这种算法对宽带麦克风阵列处理系统的滤波系数不断进行迭 代更新，在保持期望信号方向上的频响不变的条件下，使输出的噪声 功率达到最小。该算法实质上是一种带约束的LMS算法。这种算法 需要知道所求语音信号的方向，并对其事先设定一个固定的频响。在 阵列作完时间延迟补偿后，对于所求信号而言，各麦克风在对应延迟 线上的波形是完全相同的，其总的处理效果相当于一个单独的

24、阶数为 K的滤波器，其系数等于对应M个延迟抽头系数之和。所以，这MK 个系数的选择，应当满足事先给定的滤波系数值，即需要满足K个 线性约束等式。这K个约束同时也将所求信号的功率输出保持在一 定量值上。这样，对于这些系数的选择，自由度还有MK-K，它们用 来使输出总功率最小。由于事先保证了所求信号的功率，总功率最小 也就意味着系统输出的噪声功率达到最小，所以它对其它方向上的干 扰噪声能够进行有效消除。这种算法的一个主要优点在于它具有自纠正功能，不会由于程序 在运行过程中的舍入误差或截取误差的累积而偏离约束值。该处理系 统见图2.3。延时补偿x 2(n) :x ) 时延估计11TWW(n -1)

25、i(n - 2)* TWMTT1T1( n - (K -1)T (K-1)M2 MM +12M +1+ TITSI T zTy (n)K-1)M+1+.f图2.3： LCMV算法-基本宽带波束形成系统框图麦克风阵列接收信号为： TOC o 1-5 h z x (n) = s( n) + n( n)(36)其中：x(n) = x (n), x (n), , x (n - (K - 1)r(37)s(n) = s(n), , s(n),.,s(n - (K -1), , s(n - (K - 1)r(38)M个M个 n(n) = n (n顼Jn), , n/(39)假定所求信号和干扰噪声的均值均为

26、0，并且期望方向所求信号与其它方向上的干扰噪声互不相关，即Es(n)nr (n) = 0(40)阵列的输出y(n) = wt x(n) = xt (n)w(41)其中w为权系数向量w = w , w , , w 进而输出功率可以通过下式来计算(42) TOC o 1-5 h z Ey2(n) = Ewtx(n)xt(n)w = wtR w(43)约束条件要求麦克风阵列中延时相同的抽头上的权系数之和满足事先给定的滤波系数，即ctw = f , k = 0,1,., K -1(44)这里MK维向量ck表示如下c = 0,.,0,.,0,.,0,1,.,1 ,0,.,0,.,0,.,0 T(45)k

27、MM第k组(M个1)MM定义约束矩阵CMKx广C0，C，CK-1(46)及滤波系数向量F = f0, f1, , fK-1T(47)式(43)中的线性约束条件可以表示为CTw = F(48)常见的一种约束形式(即无失真约束)如下f = 1, f = 0(k = 1,., K -1)(49)0 k由于阵列注视方向上的频响由K个约束完全确定，因此所求信号的 输出功率也完全确定。这样，减少非注视方向的干扰噪声功率等同于 要求使全部输出功率最小，所以这里优化问题的目标就是使式(43)达 到最小。最后，关于MK个权系数的求解问题就变为带约束的LMS 优化问题min wtR w s.t Ctw = F(5

28、0)XXw根据拉格朗日乘子法可以得到上述优化问题的解为 TOC o 1-5 h z w = R CCtR-1C-1 F(51)opt xxxx这样就构造出一个在空间上和频域上的滤波器。将上式代入式(41), 得到信号估计值为y (n) = wt x(n)(52)在上述过程中，输入相关矩阵R假定是已知的，而实际情况下 它只能通过自适应学习方法获得。如果直接代入自相关矩阵的估计值， 由于需要求解它的逆矩阵，这就使得运算量很大，所以一般利用梯度 下降方法来求解带约束的LMS问题。(53)对于滤波器的权系数向量初始值，首先应该满足约束条件w(0) = CCt C-1 F利用梯度下降法，迭代公式应为w(

29、n +1) = w(n) RR w(n) + 人(n)C(54)xx其中人(n)为拉格朗日乘子，h为步长。新的权系数向量仍需满足约束条件，从而可得w(n +1) = w(n) - h I - C(CtC)-1 Ct R w(n) + C(CtC)-1F 一 Ctw(n)(55)xx上式右边最后一项中的因子F-Ct w(n)在理想情况下等于零，但在应 用中由于存在舍入误差或截取误差，这一项实际上并不等于零。加入 该项可以避免误差累积，使权系数向量在迭代后不致偏离优化解。令:(56)n = C(Ct C)-1 fP = I - C(Ct C)-1 Ct实际迭代中，利用每次输入信号向量的外积来作为自

30、相关矩阵的估计值。最后该算法归结为：(57)w(0)= nw(n +1) = Pw(n) - h y(n)x(n) + CCC-iF式（57）即为LCMV算法求解权矢量的表达式。实验证明，线性约束自适应波束形成方法对于相干噪声场有比较 好的消噪性能，而对于非相干噪声场，去除噪声的效果并不比延迟- 求和波束形成的效果好。2.4仿真验证仿真条件:8单元均匀线阵，阵元间距d=7cm，期望信号来向38， 干扰信号（噪声语音信号）来向5，SNR=-3dB，INR=5dB，干扰 信号和期望信号同时到达第一个阵兀（参考阵兀）。纯净语音信号（期望信号）、噪声语音信号（干扰信号）及噪声时域波形见图2.4。纯净语

31、音信号10度 幅-1 00.5110度 幅1.522.533.5Time/sed噪声语音信号-1 00.511.522.533.5Time/sed图2.4：各信号时域波形阵列接收到的信号，延时-求和法、后置滤波法及LCMV法处理后置滤波处理后的语音信号01234Time/sed延时-求和处理后的语音信号度幅0度幅012340Time/sedLCMV处理后的信号后的语音信号如图2.5所示。其中LCMV算法中滤波器阶数为J=24, 步长为r = 0.009。1度0幅-101234Time/sed图2.5：各算法处理后的语音信号由上图可见：延时-求和算法结构简单，能有效消除相干噪声， 但不具有自适应

32、干扰抑制能力，且对非相干噪声消除作用有限；后置 滤波法在延时-求和的基础上进行滤波处理，效果好于延时-求和法， 能有效消除非相干噪声，但经过播放处理后的语音，可以发现信号有 一定的失真；LCMV自适应算法，适合时变的声学环境，对相干噪声 有很好的消除能力，对非相干噪声的消除作用有限（这里滤波器的阶 数比较大J=24，如果变小性能会变差，小到一定程度时，消噪效果还不如后置滤波法)。3.麦克风实验平台的搭建1自适应天线是当前信号处理和通信领域最热门的技术之一，它的 应用可以充分利用无线资源的空间可分割性，有效地抑制干扰，提高 无线通信系统对无线资源的利用率，从根本上提高系统容量和服务质 量。为了使

33、这项技术能够尽快实用化，进行自适应天线实物的研究已 经迫在眉睫。因此，在对自适应天线算法进行理论分析和仿真的基础 上，有必要进行一些实际的测试。为此，可以通过搭建一种自适应天 线系统的实验平台，以验证各种智能天线阵列和控制算法的有效性和 实用性。但是，鉴于自适应天线系统的天线阵元加工难度较大、成本 较高，还要制作相应的射频单元，同时还需要根据实验情况适时地调 整阵元和算法，因此，在实际的电磁环境下实验和测试各种天线阵列、 自适应天线控制算法的性能存在较大的难度。为了验证自适应算法的 有效性，我们可以通过麦克风阵列来验证。3.1麦克风阵验证DBF算法的可行性分析下面提出了一种用工作原理与“自适应

34、天线系统”相同的“麦克 风阵”来构建实验平台的想法。因为自适应天线是采用一定阵元数的 阵列天线，各阵元的空间位置不同，同一无线电信号到达各个阵元波 程差也不同，因而可用于波达估计(DOA)。而不同方向的信号通过阵 列天线后产生不同的阵列响应，自适应天线接收系统可采取一定的算 法（可在射频、中频或基带实现）把不同方向的信号区分开来，从而降 低干扰，提高系统性能。因此，在自适应天线系统中，天线阵列结构 和控制算法对整个系统的性能起着重要的影响，同时也是自适应天线 系统的研究核心。麦克风阵系统也是基于麦克风阵列结构和自适应算法两大核心 内容的系统。虽然自适应天线系统与麦克风阵系统接收的信号分别是 射

35、频信号和语音信号，两者存在较大的差异，但是，对于两个系统所 研究的核心内容却是相同的，不管是射频信号还是语音信号，在系统 中都只是当作一个信号来进行分析和处理，其最终日的是为了选择最 好的阵列和最优的算法。因此，可以用麦克风阵列来代替自适应天线 阵列，在计算机平台中进行算法处理，来测试不同的自适应天线系统 中的天线阵列和控制算法。自适应天线是由多个天线单元组成的，自适应的含义就是其通过 反馈控制方式恰当的改变调整加权系数，从而连续调整自适应天线本 身方向图。有时为了使其同时具有空域和时域的处理能力，在结构上， 可以在自适应天线阵的每个天线后接一个延时抽头加权网。而自适应 麦克风阵在结构上与自适

36、应天线类似，也是通过连续调整自适应麦克 风阵的加权系数，实现辨识信号的方向和降低干扰。所以，完全可以 用麦克风阵元来代替天线阵元，利用自适应麦克风阵作为自适应天线 的实验平台，基于麦克风阵列来寻求较佳的自适应天线的天线阵列结 构和自适应算法。此外，用麦克风阵列替代天线阵列，可以大大减少了研制费用、 周期，降低制作难度，同时还不需要射频部分。一旦要进行自适应天 线系统的实验只需更改天线阵列、接入射频部分即可，也可以根据需 要在计算机中修改算法，具有通用性和灵活性。3.2实验平台的搭建下面主要介绍基于麦克风阵的自适应天线实验平台的原理、结构、 功能及主要技术指标。实验系统的原理与实现基于麦克风的自

37、适应天线实验平台的信号处理结构图见3.1。声源图3.1基于麦克风的自适应天线实验平台的信号处理结构图系统各模块功能：麦克风阵列。在构建系统时，本系统拟选用8个阵元的阵列， 并且采用间距为0.07m的均匀线列阵，较为适合于语音频段的空间采 样。放大滤波电路。该电路用来提高系统的灵敏度，滤除部分噪 声以及提供采集卡正常工作所需的阻抗和信号，包括阻抗匹配和动态 范围的匹配。采集卡。选用8通道的数据采集卡，将放大滤波后的信号进 行采集转换，并通过一定的方式传至处理系统，然后在处理系统进行 语音增强。对于采集卡有一个采样率设置问题。为了使采样后输出的 离散时间序列信号能无失真地复现原输入信号，由采样定理

38、可知采 样频率f至少应为输入信号最高有效频率/_的2倍，否则会出现 频率混淆误差。实际系统为了保证数据采样精度，一般有下列关系：f = (710) f x Nsmax式中N为多通道数采集系统的通道数。考虑到实际语音信号90%95% 的频率集中在300Hz3kHz频段，则对于16路的语音采集系统可设 定采样频率为216kHz。处理系统。处理系统完成语音信号的处理。接收到数据采 集卡传来的数据后，利用一定的自适应算法，对数据进行处理(如语 音增强等)，验证算法及天线结构的有效性。对于处理系统可以采用 计算机作为数据处理平台。系统以计算机作为数据处理平台即处理系 统，可以根据不同自适应天线实验的需要

39、，采用不同的阵列结构，选 用不同的自适应算法，能够达到通用平台的效果。另外，考虑到数据传输的稳定性及传输速度，考虑采用计算机的 PCI总线作为数据输入接口，由进行具体编程的数据采集卡来实现。由以上对实验系统各个模块的分析，用麦克风阵列实验平台验证 DBF算法的实验过程可表述为：设置信号源（声源）和干扰源（语 音信号），图3.1中的噪声为背景噪声；然后声音信号经过传播，到 达麦克风阵列。麦克风将声音信号转换成模拟的电信号，该模拟电信 号经过放大滤波后传至数据采集部分，数据采集部分又将模拟电信号 转换成计算机平台可以处理的数字信号，然后该数字信号经过输出送 至计算机平台，计算机平台采用一定的自适应

40、算法对信号进行处理， 验证自适应算法的有效性。具体实验中，可设置声源来自38方向，干扰源来自5方向， 其中干扰源也为语音信号。选用DBF中常用的自适应算法LCMV算 法对语音数据采集卡采集到的语音信号进行处理，以验证LCMV算 法。实验平台各模块的设计由上节知该实验平台主要包括4个模块：麦克风阵列、放大器滤 波、数据采集卡和计算机处理平台。简化的实验系统结构如图3.2所 示。本节对各个模块的设计和功能进行详细介绍。声源s(h, r)图3.2麦克风阵实验系统框图干扰源I (0., l)IP计算机平台(1)麦克风阵麦克风阵完成语音数据的采集，其功能等同于自适应天线系统的 天线阵元，可以灵活地采用不

41、同的结构。一般采用416天线阵元结 构，相邻阵元间距取接收信号中心频率波长的1/2。阵元间距过小， 降低接收信号相关度；阵元间距过大，将在方向图上引起不必要的波 瓣，因此阵元半波长间距通常是优选的。天线阵列形式主要有直线形 阵列、圆环形阵列或者平面阵。不同的阵列结构，它们的空间位置(方 位或仰角)也是不相同的。当然，在实际的自适应天线实验中，该阵 列结构的选择将根据所使用的自适应天线的天线阵元来决定，对于不 同的实验可以灵活地调整麦克风阵的阵列结构来进行测试，从而使该 实验平台能对不同的自适应天线系统进行实验，具有很高的通用性和 灵活性。在这里，我们使用全方向性驻极体传声器 松下公司的WM-6

42、1。一个底面直径为6mmS 3.4mm的圆柱体。其灵敏度达到-35分贝，信噪比高达62分贝。并且，在20Hz-200Hz这一频段内,WM-61的灵敏度非常稳定，确保了语音信号转化为电信号的准确性。(2)前置放大器前置放大的作用在于提高系统的灵敏度以及提供采集卡正常工 作所需的阻抗和信号，包括阻抗匹配和动态范围的匹配。前置放大以两块LM324集成运放为主体。每块LM324包含4个 放大单元，将每个放大单元设计为反相放大器，由2个放大单元对一 路麦克风信号进行放大，增益设计为20dB。对于16路语音信号，需 要 8 块 LM324。数据采集卡数据采集单元主要是将阵结构从各个方位收到的信号采集到计

43、算机平台中去，以便进一步的处理。由于输入的数据流量非常大，故 采用计算机的PCI总线作为计算机的数据输入接口。可以采用NI公司的多功能数据采集卡NI PCI-6220，如图3.3所 示。它具有16位的AD采样精度，16路模拟输入，总带宽速度为250k / s。如同时采集8路声音信号。每个通道分配的AD采样速率达30k /s。完全能够满足语音信号采样的要求。另外，数据采集卡需要通过软件对其进行设置。注册采集卡，让 计算机识别采集卡并获得采集卡号；设置采样频率采样的通道数、采 集卡号以及输入动态范围等。设置完成后开始采集数据，并将将采集 到的数据以二进制形式写入预定好的文件。PCI-6220 （低

44、价位）250 kS/s, 16 位、16 路16路模拟输入、16位分辨率250 kS/s采样速率24路数字I/O(8路高速可达1 MHz)2个32位80 MHz计数器/定时器图3.3 N1PCI-6220数据采集卡计算机平台部分计算机平台主要完成两部分工作，一是对采集的数据选用好的自 适应算法进行自适应处理，来控制和调整麦克风阵列加权网的加权系 数，改变阵列的方向性，实现对信源的定位和对信号的增强，；二是 可以将数据生成语音文件，完成对信号的复原。自适应天线系统的核心是自适应算法，因此重要的是选择较好算 法实现波束的自适应控制。日前已提出很多著名的自适应天线算法， 可以分为两类。一是在时域获得天线的最优加权，有最小均方算法 (LMS)，取样协方差矩阵的直接求逆(DMI)，递归最小均方误差算法 (RLS)等；二是在空间域对频谱进行分析，以获得DOA的估计。其 中DFT和最大熵估计(MEM)可以直接从空间频谱中估计到达角度， 而多信号分类法(MUSIC)算法和旋转不变技术参数估计法(ESPRIT) 是基于特征