MATLAB环境下ISI信道仿真及自适应均衡器设计程序说明书

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 o2、基于MSE准则的均衡器抽头系数的自适应算法为：AAck i ckek ykA其中Ck代表抽头系数向量的估值，为迭代过程中的步长参数，q为误差信号，yk代表在瞬时k包含均衡器中2k+1接收信号值的行向量。误差信号ek表示为：Q ak

2、z< ； Zk为均衡器输出，ak为已知信号序列。最初用一已 知伪随机序列ak在信道上将这个自适应均衡器进行训练。 在解调器端，均衡器用这个已知 序列去调整它的系数，一旦初始调节完成，自适应均衡器就从一个训练模式切换到直接判决AA模式，这时：e ak Zk，式中ak是检测器的输出。为了确保收敛和在慢变化信道中好的跟踪能力，选择步长参数的一种经验公式是15(2k 1)Pr式中Pr代表接收到的信号加噪声的功率，它可以从接收信号中估计出三、仿真结果图四、结论分析从结果图中我们可以看出，在信噪比逐渐增大的过程中，未经均衡器均衡的差错率没有明显改善，可知系统中始终存在码间干扰造成的误码； 经均衡器均

3、衡后的差错率则有明显但我们同时也可以看到在信噪比较低情况下，均衡器均衡之后的误码率并没有明显改 善，甚至没有未均衡的差错率低，这主要是因为噪声为随机信号，功率大时对源信号影响较 大，而且均衡器不易跟踪；当我们把均衡器的步长调低后，跟踪能力增强，差错率降低。附源程序代码：main_plot.mclear;clc;echo off;close all;N=10000;in fo=ra ndom_bi nary(N);SNR_in_dB=8:1:18;for j=1:le ngth(SNR_in_dB)y,le n=cha nn el(i nfo,SNR_in_dB(j); numoferr=0;f

4、or i=le n+1:N+le n,if (y(i)<0),decis=-1;elsedecis=1;en d;if (decis=in fo(i-5),num oferr =num oferr+1;en d;en d;Pe(j)=num oferr/N;%指定信号序列长度%产生二进制信号序列%AWGN信道信噪比%通过既有码间干扰又有白噪声信道%初始误码统计数%从第len个码元开始为真实信号码元%判决译码%判断是否误码，统计误码码元个数%未经均衡器均衡，得到的误码率end;semilogy(SNR_in_dB,Pe,'red*-');hold on;delta_ 仁 0

5、.11;delta_2=0.09;for j=1:le ngth(SNR_in_dB)y=cha nn el(i nfo,SNR_in_dB(j); z=lms_equalizer(y,i nfo,delta_1); numoferr=0;for i=1:N,if (z(i)<0),decis=-1;elsedecis=1;en d;if (decis=i nfo(i),num oferr= num oferr+1; en d;en d;Pe(j)=num oferr/N;en d;semilogy(SNR_in_dB,Pe,'blacko-');%未经均衡器，误码率结果

6、图%指定自适应均衡器的步长%指定自适应均衡器的步长%通过信道%通过自适应均衡器，并设置步长为%经自适应均衡器均衡后，得到的误码率%自适应均衡器均衡之后，误码率结果图0.11%通过信道%通过自适应均衡器，并设置步长为0.09%经自适应均衡器均衡后，得到的误码率%自适应均衡器均衡之后，误码率结果图hold on;for j=1:le ngth(SNR_in_dB) y=cha nn el(i nfo,SNR_in_dB(j); z=lms_equalizer(y,i nfo,delta_2); numoferr=0;for i=1:N,if (z(i)<0), decis=-1;elsede

7、cis=1;en d;if (decis=i nfo(i),num oferr= num oferr+1; en d;en d;Pe(j)=num oferr/N;en d; semilogy(SNR_in_dB,Pe,'blue.-'); hold on;xlabel('SNR in dB');ylabel('Pe');title('ISI信道自适应均衡系统仿真');legend('未经均衡器均衡','经自适应均衡器均衡，步长detla=0.11','经自适应均衡器均衡，步长detla=0

8、.09');ran dom_bi nary.m%产生二进制信源随机序列fun ctio n in fo=ra ndom_b in ary(N)if nargin = 0,%如果没有输入参数，则指定信息序列为10000个码元N=10000;en d;for i=1:N,temp=ra nd;if (temp<0.5),in fo(i)=-1;% 1/2的概率输出为-1elsein fo(i)=1;% 1/2的概率输出为1end en d;cha nn el.m%模拟既有码间干扰又有高斯白噪声的信道function y,len=channel(x,snr_in_dB)SNR=exp(

9、s nr_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.25 0 0.126 0.038 0.088;len _actual_isi=(le ngth(actual_isi)-1)/2;len=le n_actual_isi;y=co nv(actual_isi,x);%信号通过信道，相当于信号序列与信道模型序列作卷积%需要指出，此时码元序列长度变为N+len-1 ，译码时

10、我们从第len个码元开始到N+len个结束for i=1:2:size(y,2),no ise(i) no ise(i+1)=g ngauss(sigma);% 产生噪声en d;y=y+noise;%叠加噪声gn gauss.m%产生高斯白噪声function gsrv1,gsrv2=g ngauss(m,sgma)if n argin = 0,%如果没有输入实参，则均方为0，标准差为1m=0; sgma=1;elseif n argin = 1,%如果输入实参为1个参数，则标准差为输入实参，均值为 0sgma=m; m=0;en d;u=ra nd;z=sgma*(sqrt(2*log(1/(1-u);u=ra nd;gsrv仁m+z*cos(2*pi*u);gsrv2=m+z*s in( 2*pi*u);K=5;for k=1:size(y,2)-2*K,y_k=y(k:k+2*K);z_k=estimated_c*y_k'e_k=i nfo(k)-z_k;estimated_c=estimated_c+delta*e_k*y_k;z(k)=