通信信号处理仿真大作业.docVIP

第五章计算机仿真大作业 采用计算机编程实现图1中的自适应均衡器： 图1 信号传输的系统模型 图1中为频率为10Hz采样频率为1000Hz，经过上述信道的输出信号与高斯白噪声叠加后作为自适应均衡器的输入信号，为自适应均衡器的输出信号。其中图1中所示的自适应均衡器为一N=31阶FIR滤波器，抽头系数为其结构如图2所示： 图2自适应均衡器结构图 按要求分析回答下列问题，并给出分析结果和波形图： 1.在一个图中用子图的形式（subplot）画出图1中： (1)信号； (2)经信道传输后的信号； (3)当加的信噪比SNR(dB）为20dB时均衡器的输入序列的波形图； 对上述波形进行对比分析和说明。 分析说明：s(n)通过具有码间干扰特性的信道h(n)，由于信道存在一定的误差和码间干扰使系统的性能下降，x(n)的波形密度减小了，但整体波形没有发生变化。加入噪声后，y（n）的幅值没有变化，但整个波形由于受到噪声干扰出现”现象，波形不在平滑 %第一题：在一个图中用子图的形式画出s(n),x(n),y(n) fs=1000; %采样频率 f0=10; %模拟频率 t=0:1:fs; h=[-0.005,0.009,-0.024,0.854,-0.218,0.049,-0.0323]; s=sin(2*pi*f0/fs*t); x=filter(h,1,s); y=awgn(x,20,measured); %添加高斯白噪声 figure(1); subplot(311),plot(t,s,LineWidth,2),grid; title(正弦信号s(n),FontSize,17);set(gca,FontSize,14); subplot(312),plot(t,x,LineWidth,2),grid; title(x(n)序列,FontSize,17);set(gca,FontSize,14); subplot(313),plot(t,y,LineWidth,2),grid; title(y(n)序列,FontSize,17);set(gca,FontSize,14); 2.为了确保收敛和在慢变化信道中有好的跟踪能力，通常选择步长的一种经验公式是：，公式中的N为滤波器阶数，代表接收信号加噪声的功率。如果图2中的自适应均衡器采用的自适应算法为LMS算法，进行下列仿真实验： 当加的信噪比SNR(dB）=10dB时，在归一化功率分别为：① 0.158；②0.202；③ 0.404时的不同迭代因子下，画出： 均衡器的输出信号，并将均衡器的输出结果与均衡器输入信号y(n)和期望信号s(n)进行分析对比和说明，给出相应的结论； 分析说明：加入10dB噪声后，输入信号y(n)的波形出现严重的干扰，但在通过均衡器后输出信号z（n）在一定程度上恢复了原来信号的状态，接近期望信号是s（n），这是因为自适应均衡器能对信道的响应进行自动调节，以适应信道的时变性。 以迭代次数为横轴，在一个图中用不同线型画出不同迭代因子下的均方误差mse（Mean Squared Error, MSE，即估计值与真值之差平方 分析说明：在不同的归一化频率下，随迭代次数的增加均方误差发生一些变化，收敛速度也不一样。整体上而言，Pr越小，均方误差越小，收敛速度越快。 %第二题：SNR=10dB下，不同u下的LMS算法 clc;clear all; fs=1000;%采样频率 f0=10;%模拟频率 t=0:1:fs; h=[-0.005,0.009,-0.024,0.854,-0.218,0.049,-0.0323]; s=sin(2*pi*f0/fs*t);% x=filter(h,1,s);%x=conv(info,h); M=31; %抽头个数 N=1000; %迭代次数 Pr(1)=0.158; Pr(2)=0.202; Pr(3)=0.404; %归一化功率pr SNR(1)=10; SNR(2)=20; SNR(3)=30; %信噪比 for j=1:3 u=1/(5*Pr(j)*M); %迭代因子 y=awgn(x,SNR(1),measured); % u=0.03;