武汉理工大学《MATLAB课程设计》报告

2.1.3调制信号幅度=1.5*载波幅度

调用程序，程序中A=1.5。调制信号、已调信号的波形、频谱如图2-7所示：

图2-7调制信号、已调信号的波形、频谱图

解调信号的波形、频谱如图2-8所示：

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

14

武汉理工大学《MATLAB课程设计》报告

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

图2-9输入输出信噪比关系曲线

2.2矩形波调制

用频率300HZ矩形波调制频率30KHZ的正弦波，采用同步解调，观察调制信号、已调信号、解调信号的波形、频谱以及解调器输入输出信噪比的关系。

MATLAB源程序如下：

fc=30000; %载波频率 fs=100000; %抽样频率 N=10000; n=0:N-1; t=n/fs;

x=A*square(2*pi*300*t,50);

y=modulate(x,fc,fs,'am'); %抑制双边带振幅调制 fft1=fft(x,N); %傅里叶变换 mag1=abs(fft1); %取模 f1=(0:length(fft1)-1)'*fs/length(fft1);

15

武汉理工大学《MATLAB课程设计》报告

fft2=fft(y,N); mag2=abs(fft2);

f2=(0:length(fft2)-1)'*fs/length(fft2); figure(1);

subplot(2,2,1); plot(t,x); axis([0 0.1 -2 2]); xlabel('调制信号波形')

subplot(2,2,2); plot(f1,mag1); axis([0 5000 0 8000]); xlabel('调制信号频谱')

subplot(2,2,3); plot(t,y); axis([0 0.01 -2 2]); xlabel('已调信号波形')

subplot(2,2,4); plot(f2,mag2); axis([0 50000 0 8000]); 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 5000 0 4000]); xlabel('解调信号频谱')

16

武汉理工大学《MATLAB课程设计》报告

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

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; %高斯白噪声 %输入信噪比 解调后噪声

%输出信噪比 17

%