通信工程实验教学中心 通信系统原理实验报告
z1=z+noise; subplot(5,2,5) plot(z1);
title('加入噪声后信号波形');
%对加噪信号进行频谱分析 ba=fft(z1,512); ba=abs(ba); subplot(5,2,6) plot(ba);
title('加噪信号频谱') axis([-200 600 0 150])
%设计一个低通滤波器 Wp =50/80; Ws = 70/80;
[n,Wn] = buttord(Wp,Ws,1,5) [b,a] = butter(n,Wn);
%对加入噪声的信号进行滤波 x_fir=filter(b,1,z1);
%观察滤波之后的信号波形 subplot(5,2,7) plot(x_fir);
title('滤波之后的信号') axis([0 3000 -2 2])
%相干解调 x2=x_fir.*y; subplot(5,2,8); plot(x2);
title('与相干载波相乘波形') axis([0 3000 -0.5 2])
%对加入噪声的信号进行滤波 x3=filter(b,1,x2); subplot(5,2,9) plot(x3);
title('与相干载波相乘后滤波波形 ') axis([0 3000 -0.5 2.5])
%抽样判决 for i=1:N
if abs(x3(i*100-20))>=0.5;
xn2(i*100-99:i*100)=ones(1,100); else
xn2(i*100-99:i*100)=zeros(1,100); end end
subplot(5,2,10);
plot(xn2);
title('恢复波形') axis([0 3000 -1 2])
11