《MATLAB环境下ISI信道仿真及自适应均衡器设计程序说明》由会员分享，可在线阅读，更多相关《MATLAB环境下ISI信道仿真及自适应均衡器设计程序说明（6页珍藏版）》请在装配图网上搜索。
1、MATLAB环境下ISI信道仿真及自适应均衡器设计程序说明一、系统模型二、ISI信道仿真及LSM算法自适应均衡器原理1、发送端和接收端滤波器的级联和在采样瞬间时的信道可用等效的离散时间 FIR 信道 滤波器来表示，Xn二0.05 -0.063 0.088 -0.126 -0.25 0.9047 0.25 0 0.126 0.038 0.088， n二-5,-4,，5。2、基于MSE （均方准则）的均衡器抽头系数的自适应算法为：其中Ck代表抽头系数向量的估值，A为迭代过程中的步长参数，e为误差信号，y代表在瞬kkk时k包含均衡器中2k+1接收信号值的行向量。误差信号e表示为：ek=ak-zk；
2、zk为均衡器输出，a为已知信号序列。最初用一已 kk k k kk知伪随机序列 ak 在信道上将这个自适应均衡器进行训练。在解调器端，均衡器用这个已k知序列去调整它的系数，一旦初始调节完成，自适应均衡器就从一个训练模式切换到直接判A 决模式，这时：e - ak-z，式中a是检测器的输出。k k kk为了确保收敛 和 在慢变化信道中好的跟踪能力，选择步长参数的一种经验公式是15(2 k + 叫式中P代表接收到的信号加噪声的功率，它可以从接收信号中估计出。R三、仿真结果图四、结论分析从结果图中我们可以看出，在信噪比逐渐增大的过程中，未经均衡器均衡的差错率没 有明显改善，可知系统中始终存在码间干扰造
3、成的误码；经均衡器均衡后的差错率则有明显 改善。但我们同时也可以看到在信噪比较低情况下，均衡器均衡之后的误码率并没有明显改 善，甚至没有未均衡的差错率低，这主要是因为噪声为随机信号，功率大时对源信号影响较 大，而且均衡器不易跟踪；当我们把均衡器的步长调低后，跟踪能力增强，差错率降低。附源程序代码：main_plot.mclear;clc;echo off;%指定信号序列长度%产生二进制信号序列%AWGN信道信噪比dB(j); %通过既有码间干扰又有白噪声信道close all;N=10000; info=random_binary(N);SNR_ in _dB=8:1:18;for j=1:l
4、ength(SNR _in_ dB)y, len =channel(info, SNR _innumoferr=0;for i=len+1:N+len,if (y(i)0),decis=-1;elsedecis=1;end;if (decis=info(i-5), numoferr=numoferr+1;end;end;Pe(j)=numoferr/N;end;semilogy(SNR_in_dB,Pe,red*-);hold on;delta_1=0.11;delta_2=0.09;for j=1:length(SNR_in_dB) y=channel(info,SNR_in_dB(j);
5、z=lms_equalizer(y,info,delta_1); numoferr=0;for i=1:N,if (z(i)0),decis=-1;elsedecis=1;end;if (decis=info(i), numoferr=numoferr+1;end;end;Pe(j)=numoferr/N;end;semilogy(SNR _in _ dB, Pe ,blacko-) hold on;for j=1:length(SNR_in_dB) y=channel(info,SNR_in_dB(j); z=lms_equalizer(y,info,delta_2); numoferr=0
6、;for i=1:N,if (z(i)0),%初始误码统计数%从第len个码元开始为真实信元%判决译码%判断是否误码，统计误码码元个数%未经均衡器均衡，得到的误码率%未经均衡器，误码率结果图%指定自适应均衡器的步长%指定自适应均衡器的步长%通过信道%通过自适应均衡器，并设置步长为0.11%经自适应均衡器均衡后，得到的误码率%自适应均衡器均衡之后，误码率结果图%通过信道%通过自适应均衡器，并设置步长为0.09decis=-1;elsedecis=1;end;if (decis=info(i),numoferr=numoferr+1;end;end;Pe(j)=numoferr/N; %经自适应均
7、衡器均衡后，得到的误码率end;semilogy(SNR_in_dB,Pe,blue.-); %自适应均衡器均衡之后，误码率结果图 hold on;xlabel(SNR in dB);ylabel(Pe);title(ISI信道自适应均衡系统仿真)；legend(未经均衡器均衡,经自适应均衡器均衡，步长detla=0.11,.经自适应均衡器均衡，步长detla=0.09)；random_binary.m%产生二进制信源随机序列function info=random_binary(N)if nargin = 0,%如果没有输入参数，则指定信息序列为1 0 0 0 0个码元N=10000；end
8、；for i=1:N,temp=rand；if (temp0.5),info(i)=-1；% 1/2的概率输出为-1elseinfo(i)=1；% 1/2的概率输出为1endend；channel.m%模拟既有码间干扰又有高斯白噪声的信道function y,len=channel(x,snr_in_dB)SNR=exp(snr_in_dB*log(10)/10)；%信噪比真值转换sigma=1/sqrt(2*SNR)；%高斯白噪声的标准差%指定信道的ISI参数，可以看出此信道质量还是比较差的actual_isi=0.05 -0.063 0.088 -0.126 -0.25 0.9047 0.
9、25 0 0.126 0.0380.088;len_ actual _isi=(length(actual_isi)-1)/2; len=len_actual_isi;y=conv(actual_isi,x); %信号通过信道，相当于信号序列与信道模型序列作卷积 %需要指出，此时码元序列长度变为N+len-1,译码时我们从第len个码元开始到N+len个结束 for i=1:2:size(y,2),noise(i) noise(i+1)=gngauss(sigma); %产生噪声end;y=y+noise;%叠加噪声gngauss . m%产生高斯白噪声function gsrv1,gsrv2
10、=gngauss(m,sgma)if nargin = 0,%如果没有输入实参，则均方为0，标准差为1m=0; sgma=1;elseif nargin = 1,%如果输入实参为1个参数，则标准差为输入实参，均值为0sgma=m; m=0;end;u=rand; z=sgma*(sqrt(2*log(1/(1-u);u=rand;gsrv1=m+z*cos(2*pi*u); gsrv2=m+z*sin(2*pi*u);lm _equalizer .m%LSM算法自适应滤波器实现function z=lms_equalizer(y,info,delta) estimated_c=0 0 0 0 0 1 0 0 0 0 0;%初始抽头系数K=5;for k=1:size(y,2)-2*K,y_k=y(k:k+2*K);%获取码元，一次11个z_k=estimated_c*y_k;%各抽头系数与码元相乘后求和e_k=info(k)-z_k;%误差估计estimated_c=estimated_c+delta*e_k*y_k; %计算校正抽头系数 z(k)=z_k;%均衡后输出的码元序列end;
MATLAB环境下ISI信道仿真及自适应均衡器设计程序说明
