1、系数可变的自适应均衡器可以分为两类：基于导频的估计方法和盲估计方法。第一种方法利用数据序列中的已知数据 （可以是离散的或连续的）得到导频 位置处的信道响应，然后利用有关内差算法得到整个频域信道的响应，这种方法 简单，运算量小，但需要发送已知的导频信息，降低了系统效率。而盲估计和跟踪方法利用了接收数据的统计特性来实现信道的估计和跟踪， 如利用子空间分解算法等，相对于基于导频的估计和跟踪算法， 盲算法提高了系统效率，但极大地增加了运算量。盲均衡是一种在信道畸变相当严重的条件下， 不借助训练序列，仅根据接受 到的信号序列本身对信道进行自适应均衡的方法。 与普通的均衡器相比，盲均衡 具有收敛域大，应用
2、范围广的特点。1.3盲均衡算法与分类1.3.1盲均衡概述含有盲均衡功能的接收系统如图所示。其中信道包括收发部分的滤波器以及 空间传播媒体，其时变冲激响应序列九未知。信道输出信号形式为：r（n） hks（n k） n（n） hn s（n） n（n）, n 0, 1, 2,k为了保证无噪信道输出u（n） hn s（n）方差不变，通常采用自动增益控制技 术，使得 h2 1。令wi为一个理想逆滤波器的冲激响应序列，他与信道冲激响应序列hn之间满足逆关系，即Wg i ni这样，在发射信号通过信道传输后，首先接入这个逆滤波器，其输出为（先 不考虑噪声因素）：wir（n i） s（n l） hl i s（n
3、 l） l s（n）。i l i l在实际应用中，理想逆滤波器w通常采用长度为2L+1的有限抽头，这样 滤波器输出为Ly（n） W?s（n i）i L这就是众所周知的用横向滤波器实现逆滤波器的形式。由于逆滤波器截断， 必然会带来残余码间干扰，进一步分析可知： y（n） s（n） v（n）,其中v（n） W（n） w（n）s（n i）,? 0 i L称为卷积噪声，也就是残余码间干扰。以此作为误差信号去调节逆滤波器就得 到盲均衡器。1.3.2盲均衡算法分类考虑一个有2N+1抽头的线性均衡器如下图所示。其中NZn Cimyn i，式中m和n取整数,yn为第NT时刻均衡器的输出参数，i N为第m次高速
4、后第i个抽头的增益系数，T为发送端信号的符号周期。算法的一般形式为Cim1 Cim oynif（Zn），这里是迭代步长，f（）是起误差控制的函数，其选取关系到算法的收敛性。Sato提出的盲均衡算法表达式为2f（Zn） Zn r.sign（乙），其中 r E（an）/E|an|;Godard给出的盲均衡算法表达式为K 2 K 2 K kf（Zn） Zn|Zn| （| Zn | Rk ）,其中 Rk E | 务 | / E | an | ；Serra给出的盲均衡算法表达式为f（Zn） | Zn |k .sig n（Zn） RK.Zn，其中 R E | an |K 1/E（an2）；Ben ve n
5、i sete-Goursat 提出的均衡算法表达式为f （Zn） K1en K2 | en | en，其中 en Zn an ；以上各种算法的盲均衡器总的要求是快速跟踪信道的变化，快速收敛，且收敛以后的剩余误差要小2. CMA算法2.1CMA算法的原理利用自适应滤波算法，合理的人工制造一个“期望响应”来代替缺失的“期 望响应”。其实，人工制造一个期望响应的思想，在非盲均衡器的应用中已经被 采用，即训练序列，但训练序列只在初始系统训练阶段存在，一旦训练结束，训 练序列不再存在，通信系统将传输用户的有用数据，期望响应也不再存在，自适 应滤波器切换成一个固定系数滤波器，对于平稳信道来讲这样做是可以接
6、受的， 但对于性能不稳定的信道，接收机性能将会显著下降。对原理加以改进，在训练序列传输结束后，通过人造一个期望响应，使得自 适应滤波过程能够继续，以保证自适应均衡器跟踪信道的变换。 人造“期望响应” 的方法是，在训练结束后，将均衡器输出送入判决器，判决器的输出作为期望响 应，与滤波器输出相减构成误差量用于调整自适应均衡器系数。 由于判决器运算是一种非线性运算，因此训练结束后，利用人造期望响应的自适应均衡算法不再 是线性自适应滤波器，而是非线性自适应滤波器。下图表示了 CMA盲均衡算法的框图。图2.1 CMA盲均衡算法框图在通信系统中，角度调制是常用的调制形式，它包括频率调制（ FM ）和相位调
7、制（PM），这些调制信号满足包络是常数的性质， 利用这个性质，构造一类盲自适应均 衡算法，即CMA算法。传输信号满足恒模性，即|s（n）|2 R2，因为接收到的 信号经过信道引起了畸变并且混入了干扰噪声， 已不满足恒模性，当接收到的信号通过均衡器后，如果性能得到改善，误差函数g （n） = _ R、会下降，理想的均衡器是误差函数下降到零。定义 （y（d）=（ （n） ）2 = （ lv（n） | - Rj 2使（y（n）最小，利用LMS算法的基本思路，可以导出 CAM算法如下dv （n I 1） =w （n） / 屮（v（n）补尉）=w （ti） （ v（m） 2 R ? ） y （n） x
8、（n）对于复信号和复系统，权更新算法为w n-Fl） =w （n） / （ v（h） 2 Ri ） y Hi） x* （n）2.2CMA 算法的MATLAB 程序实现先以4QAM调制为例。第一步：初始化。取1000个数据，调制方式为4QAM，从星座可知，其模为常数，步长为0.02，信道冲激响应随机生成，为复信道。第二步：生成信道噪声。第三步：通过CMA均衡器处理。第四步：计算SER。程序如下:% QAM的CMA算法实现%初始化T=1000;dB_max=30;dBn ter=3;N=5;Lh=5;Ap=4;h=ra ndn（ Ap,Lh+1）+sqrt（-1）*ra ndn（ Ap,Lh+1）
9、;for i=1:Ap, h（i,:）=h（i,:）/no rm（h（i,:）;ends=rou nd（ra nd（1,T）*2-1;s=s+sqrt（-1）*（rou nd（ra nd（1,T）*2-1）;SER=zeros（1,dB_max）;for dB=0:dBn ter:dB_max%产生信道噪声x=zeros（Ap,T）;SNR=zeros（1,Ap）;Apx（i,:）=filter（h（i,:）,1,s）;vn=randn （1,T）+sqrt（-1）*ra ndn （1,T）;vn=v n/norm（v n）*10A（-dB/20）* norm（x（i,:SNR（i）=20*l
10、og10（ norm（x（i,:）/norm（v n）; x（i,:）=x（i,:）+v n;% CMA盲均衡器Lp=T-N;X=zeros（N+1）*Ap,Lp）;Lpfor j=1:X（j-1）*（N+1）+1:j*（N+1）,i）=x（j, i+N:-1:i）.;e=zeros（1,Lp）;f=zeros（N+1）*Ap,1）; f（N*Ap/2+3）=1;R2=2;mu=0.001;e（i）=abs（f*X（:,i）A2-R2;f=f-mu*2*e（i）*X（:,i）*X（:,i）*f;sb=f*X;%计算SERH=zeros（N+1）*Ap,N+Lh+1）; temp=O;N+1te
11、mp=temp+1;H（temp,i:i+Lh）=h（j,:）;fh=f*H;temp=0;temp=fi nd（abs（fh）=max（abs（fh）;sb1=zeros（1,size（sb）;sb1=sb./（fh（temp）;sb仁sig n（ real（sb1）+sqrt（-1）*sig n（imag（sb1）;start=N+1-temp;sb2=sb1（10:le ngth（sb1）-s（start+10:start+le ngth（sb1）;SER（dB+1）=le ngth（fi nd（sb2=0）/le ngth（sb2）;%画图if 1figure（1）;subplot（2
12、21）,plot（s,ogrid,title（Transmitted symbols xlabel（Real）,ylabel（lmage）axis（-2 2 -2 2）subplot（222）,plot（x,grid, title（Received samples）, ylabel（Imagesubplot（223）,plot（sb,Equalized symbols）, xlabel（figure（2）;plot（abs（e）;C on verge nee nError e（n）figure （3）;i=O:dB_max;semilogy（i,SER（i+1）,gp-grid;lege nd
13、（SGDCMAylabel（误码率xlabel（信噪比 dB figure（4）;h=resh ape（h,1,（Ap*（Lh+1）;ii=1:（N+1）*Ap;stem（ii,h（ii）;cha nnel impluse resp on sestem（ii,f（ii）;equalization coefficienee 生成的星座对比图如下: . . .jZ1 .-._ XjL?口i-j a :1 1 B Transmitted symbols1iu4-1Received sampleso6E-2 0 2RealEqualized symbole2 -10 12从这张对比图可以看出，当采用C
14、MA盲均衡以后，盲均衡输出汇聚到四个星座点上，这样在判决的时候将极大提高判决准确率。Convergence54 3c=ot山O0 oon _yl no8 on71 no6 no6 on .II on20c no-该图表示了 QAM经过盲均衡处理器以后的收敛曲线2.3CMA算法和LMS算法的性能比较LMS算法是一种线性自适应滤波算法。LMS算法包括两个过程：一个是滤 波过程，一个是自适应过程。在滤波过程中，自适应滤波器计算其对输入的响应， 并且通过与期望响应比较，得到估计的误差信号。在自适应过程中，系统估计误 差自动调整滤波器的参数。M对于FIR横向滤波器y（n） WjX（n i 1） wTx（
15、n），使用最小均方误差i 1T 2（LMS）作为代价函数J（w） E y（n） wTx（n），在最小均方误差意义下的最佳权向量wopt argmin J（w） R 1r，其中R Ex（n）xT（n）,r Ex（n）y（n）。动态系统中，加权向量应该根据观测信息自适应调整， 应用最广的是下降算法，即卩w（n） w（n 1） （n）v（n）， （n）为更新步长，v（n）为更新方向向量。基本LMS算法：又称为最陡下降法，更新方向向量为 n-1次迭代代价函数的负梯度，v（n） Jw（n 1），为了简化梯度计算量，通常用估计值 ?（n） 2e（n）x（n），其中误差信号定义为期望输出与滤波器实际输出之间
16、的误差 e（n） d（n） wT（n 1）x（n）。如果期望信号未知，也可以用y（n）代替d（n）。根 据更新步长的不同又分为三种情况：（1） （n） constant为基本 LMS 算法；（2） （n）xT（ n）x（n）（0,2）, 0为归一化LMS算法;率归一化LMS算法， 为遗忘因子，M为加权向量维数。基本LMS算法是由Windrow 在60年代初提出的。时域解相关LMS算法：上述LMS算法收敛速度慢，而解相关可以显著加快收敛速度。定义x（n）和x（n1）在n时刻的相关系数a（n）总噱?1），更新方向向量为v（n） x（n）a（n）x（n 1），更新步长（n） e（n）/xT（n）v（
17、n），0 1 为修正因子。此法是 Doherty 在 1997 年提出的。变换域解相关LMS算法：通过对输入数据进行酉变换，在不增加计算复杂度的前提下，提高收敛速度。首先给定一个酉变换矩阵 SSH I，u（n） Sx（n） ； e（n） y（n） WH（n 1）u（n） ;W?（n） W（n 1） （n）u（n）e（n）。下面给出程序：%psk的盲均衡分别用CMA和LMSclear allM=4; k=log2（M）; n=5000; %u=0.05; u1=0.001; u2=0.0001;m=500;%h=0.05 -0.063 0.088;%-0.126; -0.25;h=1 0.3 -
18、0.3 0.1 -0.1;L=7;mse1_av=zeros（1, n-L+1）;mse2_av=mse1_av;ma=ra ndi nt（1, n,M）;a仁pskmod（a,M）;m1= abs（a1）.A4;m2=abs（a1）.A2;r1=mea n（m1）;r2=mea n（ m2）;R2=r1/r2;%R2=sqrt（2%）;s=filter（h,1,a1）;sn r=15;x=awg n（ s,s nr,measuredc1=0 0 0 1 0 0 0 ;c2=c1;y=zeros（ n-L+1,2）; n-L+1y=x（i+L-1:i）; z1（i）=c1*yz2（i）=c2*
19、ye仁 R2-（abs（z1（i）A2）;e2=a1（i）-z2（i）;c1= c1+u1*e1*y*z1（i）;c2=c2+u2*e2*y;mse1（i）=e1A2;%u（i）=0.2*（1-exp（-（0.3*abs（e（i）; mse2（i）=abs（e2F2;end;mse1_av=mse1_av+mse1;mse2_av=mse2_av+mse2;mse1_av=mse1_av/m;mse2_av=mse2_av/m;figureplot（1: n-L+1,mse1_av,r,1: n-L+1,mse2_av,baxis（0,5100,0 2.8）;scatterplot（a1,1,
20、0,r*hold onscatterplot（x,1,0,g*scatterplot（z1（2300:4800）,1,0,hold offscatterplot（z2（2300:1.5U.5-0.5Scatter plot-1.6 -1 -0.5 0 0.6 1 1.5-1.5.5O.巴 nlEpenQ-1 s5 L H I i d 可-1.5 1 -0.6 0 0.5 1 1.5IrrP hase上面两张图表示了分别米用 CMA算法和LMS算法的星座图。从两张图可以看出，CMA比LMS性能要好，在判决的时候更不容易出现判决错误参考资料：1. 通信中的自适应信号处理.邱天爽 魏东兴等编著.电子工业出版社。2 . 信号处理的自适应理论谢胜利 何昭水等编著.科学出版社。3.盲信号处理及应用.张发启等编著.西安电子科技大学出版社。4.DSE-CMA: 种新的常数模盲均衡算法.朱小刚等.上海交通大学电子 工程系。5.基于混沌的通信系统的信道盲均衡.王世元等.西南师范大学电子信息 工程学学院。
