综合课程设计基于Matlab的自适应均衡器设计

1、电子信息系综合课程设计基于matlab的自适应均衡器设计专业名称通信工程班级学号学生姓名指导教师设计时间2010.12.202011.1.7课程设计任务书专业： 学号： 学生姓名（签名）： 设计题目：基于matlab的自适应均衡器设计一、设计实验条件实验室，matlab软件二、设计任务及要求1. 课题要求系统学习时域均衡原理，掌握理论知识；2. 首先进行时域均衡原理和算法设计，再在所用的仿真软件matlab上对设计进行仿真分析，最后写实验报告；3. 对整个系统设计进行回顾，总结心得。三、设计报告的内容1. 设计题目与设计任务（设计任务书）2. 前言（绪论）(设计的目的、意义等)3. 设计主体（

2、各部分设计内容、分析、结论等）4. 结束语（设计的收获、体会等）5. 参考资料四、设计时间与安排1、设计时间： 3周2、设计时间安排： 熟悉实验设备、收集资料: 4天设计图纸、实验、计算、程序编写调试: 7天编写课程设计报告: 3天答辩: 1天基于matlab的自适应均衡器设计一、 设计目的及意义：通过本学期通信原理课程的学习，主要对数字信号系统的通信原理、传输机制等有了系统深入的了解。而实践性的课程设计能够起到提高综合运用能力，加强理论知识的学习，提高实验技术，启发创造新思想的效果。此次课程设计是自适应均衡器设计。我们按照查找资料、软件选择、系统设计、仿真实现、结果优化这一流程进行。不仅使我

3、们进一步巩固了课程知识，也提高了我们分析问题、解决问题的能力。二、设计主体：1 、设计原理数字信号经过这样的信道传输以后，由于受到了信道的非理想特性的影响，在接收端就会产生码间干扰(isi)，使系统误码率上升，严重情况下使系统无法继续正常工作。理论和实践证明，在接收系统中插入一种滤波器，可以校正和补偿系统特性，减少码间干扰的影响。这种起补偿作用的滤波器称为均衡器。时域均衡是利用均衡器产生的时间波形去直接校正已畸变的波形，使包括均衡器在内的整个系统的冲击响应满无码间串扰条件。频域均衡是从校正系统频率特性出发，使包括均衡器的基带系统的总特性满足无失真传输条件；频域均衡在信道特性不变，且在传输低速数

4、据时是适用的。而时域均衡可以根据信道特性的变化进行调整，能够有效地减小码间串扰，故在高速数据传输中得以广泛应用。图1数字基带传输系统在实际中，当数字基带传输系统（如图1）的特性不满足奈奎斯特第一准则时，就会产生有码间串扰的响应波形。现在我们来证明：如果在接收滤波器和抽样判决之间插入一个称之为横向滤波器的可调滤波器，其冲激响应为其中， 完全依赖于h()，那么，理论上就可消除抽样时刻上的码间串扰。由 上式可以看出，这里的是图2所示网络的单位冲激响应。该网络是有无限多的按横向排列的迟延单元和抽头加权系数组成的，因此成为横向滤波器。它的功能是利用它产生的无限多个响应波形之和，将接收滤波器输出端抽样时刻

5、上有码间串扰的响应波形变换成抽样时刻上无码间串扰的响应波形。由于很想滤波器的均衡原理是建立在响应波形上的，故把这种均衡成为时域均衡。图2横向滤波器不难看出，横向滤波器的特性将取决于各抽头系数。如果是可调整的，则图3所示的滤波器是通用的；特别是当可自动调整时，则它能够适应信道特性的变化，可以动态校正系统的事件响应。图3有限长横向滤波器及其输入和输出波形理论上，无限长的横向滤波器可以完全消除抽样时刻上的码间串扰，但实际中是不可能实现的。因为，不仅均衡器的长度受限制，并且系数的调整准确度也受到限制。如果的调整准确度可不到保证，即使增加长度也不不会获得显著的效果。因此，有必要进一步讨论有限长横向滤波器

6、的抽头增益调整问题。设一个具有2n+1个抽头的横向滤波器，如图3（a）所示，其单位冲激响应为，则有有设它的输入为，是被均衡的对象，并设它没有附加噪声，如图3（b）所示，则均衡后输出波形为在抽样时刻 (设系统无延时时)上，有上式说明，均衡器在第k个抽样时刻上的样值将有2n+1个与乘积之和来决定。显然，其中除以外的所有都属于波形失真引起的码间串扰。当输入波形给定，及各种可能的确定是，通过调整使指定的等于零时容易办到的，但同时要求所有的（除k=0外）都等于零却是一件很难的事。这也说明，利用有限长的横向滤波器减小码间串扰是可能的，但完全消除是不可能的。为了反映这些失真的大小，我们可采用最小均方和迫零算

7、法准则来衡量。本实验采用最小均方误差算法（lms）。lms（least mean square）算法最早由widrow于1960年建立。采用最小均方算法的均衡器比迫零算法均衡器要稳定一些，它所采用的准则是使均衡器的期望输出值和实际输出值之间的均方误差(mse)最小化的准则。由wiener滤波理论可知，滤波器的最优抽头系数矢量可表示为式中是输入信号矢量的互相关函数，与期望信号矢量的互相关函数。直接求解式此式的运算量非常大，尤其当滤波器的抽头系数个数比较多时。考虑如图4所示的自适应fir滤波器：图4 自适应fir滤波器 令表示滤波器在n时刻的误差，并定义均方误差 为代价函数。由梯度的定义可以得出，

8、代价函数相对于滤波器抽头系数矢量w梯度为定义梯度矢量为 以及输入矢量和抽头系数矢量： 则此式可以写作以下矢量形式： 最广泛使用的自适应算法形式为“下降算法”： 式中是第n 步迭代的系数矢量，是第n步迭代的更新步长，而v(n)是第n 步迭代的更新方向矢量。最常用的下降法为梯度下降法，又称最陡下降法。在这种算法里，更新方向矢量v(n)取作第n-1 次迭代的代价函数j（n-1）的负梯度，即最陡下降法的统一形式为 系数为了使得到的更新公式更加简单可得抽头系数矢量的更新公式为 从上式中可以看出，滤波器抽头系数矢量的更新迭代与和有关，而实际的信号处理过程中，和都不是能够先验知道的，因此通常用梯度的估计值来

9、代替. 可得 其中为滤波器输出与期望信号的误差函数。上式就是lms算法抽头系数的更新公式，容易验证。的无偏估计，通过上面的讨论可以总结出lms算法的具体步骤如下：步骤l ：初始化抽头系数矢量；步骤2 ：当n=n+1时计算滤波器输出，计算误差函数，更新滤波器抽头系数矢量关于上面的lms算法有以下几点注释：注释l ：若取常数，则称为基本lms 算法。注释2 ：若取则的到归一化lms算法。注释3 ：在功率归一化算法中，取递推计算，这里为遗忘因子，由0<a<确定，而m是滤波器阶数.注释4 ：当期望信号未知时，步骤2中的d(n)可直接用判决器的实际输出代替。本实验中选取滤波器阶数order=

10、63，调整步长u=0.02，均为经验值。2 、matlab仿真过程1.定义一个随机信道c=0.4;x1=c*rand();x2=c*rand();x3=c*rand();%信道参数isi=1，x1，x2，x3;snr=10; %awgn信道信噪比2.定义滤波器各参数order=63; %滤波器阶数len=1000; %训练序列长度u=0.02; %调整步长m=(order-1)/2;n=len+length(isi)-1;e=zeros(1，n);error=e;y=zeros(1，n); %fir滤波器输出3.信道均衡过程for i=1:100 x=sign(rand(1，len)-0.5)

11、;%采用pn码作为训练序列 noise=randn(1，n)/10.(snr/10);%awgn y=conv(isi，x)+noise;%输出 wk=zeros(1，order);%横向滤波器系数 for n=order:n-m+1 %lms算法 y1=y(n+m-1:-1:n-m-1); d1=wk*y1' e(n)=x(n)-d1; wk=wk+u*e(n)*y1; e(n)=10*log10(abs(e(n); endq error=error+e;end4.得到横向滤波器系数wk输出结果:wk = columns 1 through 9 -0.0058 0.0297 0.02

12、83 0.0247 -0.0147 -0.0118 0.0385 -0.0348 0.0095 columns 10 through 18 0.0422 0.0321 0.0236 0.0126 -0.0068 -0.0310 0.0094 0.0207 -0.0136 columns 19 through 27 0.0354 0.0016 -0.0281 0.0181 0.0014 0.0073 -0.0104 0.0038 0.0422 columns 28 through 36 0.0149 0.0396 0.0150 0.9751 -0.2937 -0.2446 0.1236 0.08

13、25 -0.0315 columns 37 through 45 -0.0098 -0.0055 0.0035 -0.0426 -0.0316 0.0045 0.0007 -0.0574 0.0115 columns 46 through 54 0.0346 -0.0578 0.0389 -0.0503 -0.0110 0.0133 0.0500 -0.0114 -0.0143 columns 55 through 63 -0.0055 -0.0214 0.0119 0.0058 -0.0181 -0.0184 -0.0392 0.0115 -0.02725.计算自适应均衡器误差收敛性erro

14、r=error(order:n-m+1)/100;%误差值t=1:length(error);plot(t，error);title('误差收敛曲线');ylabel('误差db');xlabel('n');图5误差收敛曲线6. 随机信号通过理想均衡信道m = 4; msg = randint(1500，1，m); %产生一个1500 x 1 在0，m-1区域内的随机整数序列modmsg = qammod(msg，m); %进行qam调制的映射trainlen = 500; %定义训练序列的帧长chan = .986; .845; .237; .

15、123+.31i; %定义信道参量filtmsg = filter(chan，1，modmsg); % 模拟信道变化eq1 = lineareq(8， lms(0.01);%生成均衡器对象eq1.sigconst = qammod(0:m-1，m); % 设置星座图symbolest，yd = equalize(eq1，filtmsg，modmsg(1:trainlen); h = scatterplot(filtmsg，1，trainlen，'bx'); hold on; %以蓝星画出未经过均衡的信号scatterplot(symbolest，1，trainlen，'

16、g.'，h); %在原图的基础上以绿色画出经过均衡后的信号scatterplot(eq1.sigconst，1，0，'k*'，h); %在原图的基础上以黄色标出理想星座图legend('filtered signal'，'equalized signal'， 'ideal signal constellation');hold off;% 计算误码率demodmsg_noeq = qamdemod(filtmsg，m); %解调未均衡的码字demodmsg = qamdemod(yd，m); %解调已均衡的码字nnoeq，

17、rnoeq = symerr(demodmsg_noeq(trainlen+1:end)，.对比未均衡信号和样本信号的误码和误码率 msg(trainlen+1:end);neq，req = symerr(demodmsg(trainlen+1:end)，. msg(trainlen+1:end); %对比均衡信号和样本信号的误码和误码率disp('symbol error rates with and without equalizer:')disp(req rnoeq)图6信号星座图输出结果: symbol error rates with and without equa

18、lizer:0 0.32303、实验结果分析输入信号通过参数可调的数字滤波器后产生输出信号，将其与期望相应进行比较，形成误差信号，并以此通过lms算法对滤波器参数进行逐步调整，最终使误差均方值最小，并得到自适应均衡器各阶参数。随机信号通过理想信道均衡器可完全消除码间串扰。三、实验心得通过分析课题和动手实践，培养了对抽象的实际问题进行逻辑抽象，以确定输入输出及其关系进而进行分析的能力。了解并掌握正确运用matlab各种函数在数字信号处理中的作用。做课程设计同时也是对课本知识的巩固和加强。由于我主要负责系统的设计和仿真，在对课题理解和整体设计中对课本知识有了更深一步的了解。通过在图书馆细心地查找，也寻找到了很多有关书籍文献，对我们的设计有很大帮助，从中再一次证明了书籍永远是我们生命的源泉，增强了我们的自学能力。生活就是这样，汗水预示着结果也见证了收获。劳动是人类生存生活永恒不变的话题。通过实习，我才真正领略到“艰苦奋斗”这一词的真正含义。我想说，设计确实有些辛苦，但苦中也有乐，在如今单一的理论学习中，很少有机会能有实践的机会，但我们可以，而且设计也是一个团队的任务，一起的工作可以让我们有说有笑，相互帮助，配合默契，多少人间欢乐在这里洒下，大学里一年的相处还赶不上这十来天的合作，我感觉我和同学们之间的距离更加近了；我想说，确实很累，但当我们看到自己所做的成果时，心中也不免产生