信号分析处理设计书(编辑修改稿)

信号分析处理设计书(编辑修改稿)内容摘要：

%抽样频率 fs=100000。 t=[0:1/Fs:]。 %抽样间隔 Fc=30000。 %载波频率 a=1。 %设 a=1 m=a*cos(300*2*pi*t)。 %调制信号 X=fft(m)。 %对 m 进行傅里叶变换 X=abs(X(1:length(X)/2+1))。 %调制信号频谱 frqX=(0:length(X)1)*Fs/length(X)/2 S = modulate(m,Fc,Fs,39。 amssb39。 )。 %对信号进行调制 武汉理工大学《信号分析处理 》课程设计说明书 5 Y=fft(S)。 %对 S 进 行傅里叶变换 Y=abs(Y(1:length(Y)/2+1))。 frqY=(0:length(Y)1)*Fs/length(Y)/2。 %已调信号频谱 set(gcf,39。 color39。 ,39。 w39。 ) %颜色设置 figure(1) subplot(421) %绘制调制信号时域图 plot(t,m) title(39。 调制信号 时域图 39。 ) subplot(422) %绘制调制信号频谱图 plot(frqX,X) axis([0 3000 0 max(X)]) title(39。 调制信号频谱 图 39。 ) subplot(423) %绘制 已调信号 时域图 plot(t,S) title(39。 已调信号 时域图 39。 ) subplot(424) plot(frqY,Y) %绘制 已调信号频谱 图 axis([0 60000 0 max(Y)]) title(39。 已调信号频谱 图 39。 ) sn=awgn(S,4)。 %加入高斯白噪声 sn1=awgn(S,10)。 sn2=awgn(S,15)。 sn3=awgn(S,20)。 sn4=awgn(S,25)。 Y1=demod(S,Fc,Fs,39。 amssb39。 )。 %无噪声已调信号解调 YYN=demod(sn,Fc,Fs,39。 amssb39。 )。 %加噪声已调信号解调 YYN1=demod(sn1,Fc,Fs,39。 amssb39。 )。 YYN2=demod(sn2,Fc,Fs,39。 amssb39。 )。 YYN3=demod(sn3,Fc,Fs,39。 amssb39。 )。 武汉理工大学《信号分析处理 》课程设计说明书 6 YYN4=demod(sn4,Fc,Fs,39。 amssb39。 )。 J1=fft(sn)。 J1=abs(J1(1:length(J1)/2+1)) frqJ1=(0:length(J1)1)*Fs/length(J1)/2。 %加噪声后已调 信号频谱 J2=fft(YYN) J2=abs(J2(1:length(J2)/2+1)) frqJ2=(0:length(J2)1)*Fs/length(J2)/2 %加噪声后解调信号频谱 set(gcf,39。 color39。 ,39。 w39。 ) figure(1) subplot(425) %绘制 加噪声解调信号 时域图 plot(t,YYN)。 title(39。 加噪声解调信号 时域图 39。 ) subplot(426) %绘制 加噪声解调信号频谱 图 plot(frqJ2,J2)。 axis([0 3000 0 max(J2)]) title(39。 加噪声解调信号频谱 图 39。 ) subplot(427) %绘制 无噪声解调信号 时域图 plot(t,Y1) title(39。 无噪声解调信号 时域图 39。 ) subplot(428) %绘制 无噪声解调信号频谱 图 plot。