DSB信的仿真分析

xlabel('已调信号频谱')

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% yn=awgn(y,4); %加入高斯白噪声 znn=demod(y,fc,fs,'am'); %无噪声已调信号解调 zn=demod(yn,fc,fs,'am'); %加噪声已调信号解调 fft3=fft(znn,N); mag3=abs(fft3); f3=(0:length(fft3)-1)'*fs/length(fft3); figure(2); subplot(3,1,1); plot(t,zn); xlabel('加噪声解调信号波形') subplot(3,1,2); plot(t,znn); xlabel('无噪声解调信号波形') subplot(3,1,3); plot(f3,mag3); axis([0 500 0 500]); xlabel('解调信号频谱')

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

yn1=awgn(y,8); yn2=awgn(y,12); yn3=awgn(y,16); yn4=awgn(y,20);

zn1=demod(yn1,fc,fs,'am'); zn2=demod(yn2,fc,fs,'am'); zn3=demod(yn3,fc,fs,'am'); zn4=demod(yn4,fc,fs,'am'); dyi=yn-y; %高斯白噪声 s_ni=var(y)/var(dyi); %输入信噪比 dyo=zn-znn; %解调后噪声 s_no=var(znn)/var(dyo); %输出信噪比 dyi1=yn1-y; s_ni1=var(y)/var(dyi1); dyo1=zn1-znn; s_no1=var(znn)/var(dyo1);

dyi2=yn2-y; s_ni2=var(y)/var(dyi2); dyo2=zn2-znn; s_no2=var(znn)/var(dyo2); dyi3=yn3-y; s_ni3=var(y)/var(dyi3); dyo3=zn3-znn; s_no3=var(znn)/var(dyo3); dyi4=yn4-y; s_ni4=var(y)/var(dyi4); dyo4=zn4-znn; s_no4=var(znn)/var(dyo4); in=[s_ni,s_ni1,s_ni2,s_ni3,s_ni4]; out=[s_no,s_no1,s_no2,s_no3,s_no4]; figure(3); plot(in,out,'*')

hold on plot(in,out)

xlabel('输入信噪比'); ylabel('输出信噪比') ×载波幅度

调用程序，程序中A=0.8。 调制信号、已调信号的波形、频谱如图2-1所示： 图2-1 调制信号、已调信号的波形、频谱图 解调信号的波形、频谱如图2-2所示： 图2-2解调信号的波形、频谱图 输入输出信噪比关系曲线如图2-3所示： 图2-3 输入输出信噪比关系曲线 调用函数，函数中A=1。调制信号、已调信号的波形、频谱如图2-4所示： 图2-4调制信号、已调信号的波形、频谱图 解调信号的波形、频谱如图2-5所示：

图2-5解调信号的波形、频谱图

输入输出信噪比关系曲线如图2-6所示：