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。2、基于MSE准则的均衡器抽头系数的自适应算法为：其中代表抽头系数向量的估值，为迭代过程中的步长参数，为误差信号，代表在瞬时k包含均衡器中2k+1接收信号值的行向量。误差信号表示为：；为均衡器输出，为已知信号序列。

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

3、有未均衡的差错率低，这主要是因为噪声为随机信号，功率大时对源信号影响较大，而且均衡器不易跟踪；当我们把均衡器的步长调低后，跟踪能力增强，差错率降低。附源程序代码：main_plot.mclear;clc;echo off;close all;N=10000; %指定信号序列长度info=random_binary(N); %产生二进制信号序列SNR_in_dB=8:1:18; %AWGN信道信噪比for j=1:length(SNR_in_dB) y,len=channel(info,SNR_in_dB(j); %通过既有码间干扰又有白噪声信道 numoferr=0; %初始误码统计数 for

4、 i=len+1:N+len, %从第len个码元开始为真实信号码元 if (y(i)<0), %判决译码 decis=-1; else decis=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; %指定自

5、适应均衡器的步长for j=1:length(SNR_in_dB) y=channel(info,SNR_in_dB(j); %通过信道 z=lms_equalizer(y,info,delta_1); %通过自适应均衡器，并设置步长为0.11 numoferr=0; for i=1:N, if (z(i)<0), decis=-1; else decis=1; end; if (decis=info(i), numoferr=numoferr+1; end; end; Pe(j)=numoferr/N; % 经自适应均衡器均衡后，得到的误码率end;semilogy(SNR_in_dB

6、,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); %通过自适应均衡器，并设置步长为0.09 numoferr=0; for i=1:N, if (z(i)<0), decis=-1; else decis=1; end; if (decis=info(i), numoferr=numoferr+1; end; end; Pe(j)=numoferr/

7、N; % 经自适应均衡器均衡后，得到的误码率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%产生二进制信源随机序列func

8、tion info=random_binary(N) if nargin = 0, %如果没有输入参数，则指定信息序列为10000个码元 N=10000;end;for i=1:N, temp=rand; if (temp<0.5), info(i)=-1; % 1/2的概率输出为-1 else info(i)=1; % 1/2的概率输出为1 endend;channel.m%模拟既有码间干扰又有高斯白噪声的信道function y,len=channel(x,snr_in_dB)SNR=exp(snr_in_dB*log(10)/10); %信噪比真值转换sigma=1/sqrt(2*

9、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=(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)

10、 noise(i+1)=gngauss(sigma); %产生噪声end;y=y+noise; %叠加噪声gngauss.m%产生高斯白噪声function gsrv1,gsrv2=gngauss(m,sgma) if nargin = 0, %如果没有输入实参，则均方为0，标准差为1 m=0; sgma=1;elseif nargin = 1, %如果输入实参为1个参数，则标准差为输入实参，均值为0 sgma=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);lms_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: