nst=awgn(s_2ask,70); %解调部分 nst=nst.*ht; figure(2)

subplot(311); plot(t,nst);

axis([0 50 -0.2 1.2]);

title('乘以相干载波后的信号') %低通滤波器设计 wp=2*pi*2*fc*0.5; ws=2*pi*2*fc*0.9; Rp=2; As=45;

[N1,wc]=buttord(wp,ws,Rp,As,'s'); [B,A]=butter(N1,wc,'s'); %低通滤波

h=tf(B,A); %转换为传输函数 dst=lsim(h,nst,t); subplot(312); plot(t,dst);

axis([0 50 -0.2 1.2]);

title('经过低通滤波器后的信号'); %判决器 k=0.125;

pdst=1*(dst>0.125); subplot(313); plot(t,pdst);

axis([0 50 -0.2 1.2]);

title('经过抽样判决后的信号') %频谱观察

%调制信号频谱 T=t(end); df=1/T;

N3=length(s_2ask); f=(-N3/2:N3/2-1)*df;

sf=fftshift(abs(fft(s_2ask))); figure(3) subplot(411) plot(f,sf)

title('调制信号频谱') %信源频谱

mf=fftshift(abs(fft(d_NRZ(1:length(t))))); subplot(412) plot(f,mf)