信號(hào)分析處理設(shè)計(jì)書(編輯修改稿)
2025-06-29 18:50 本頁面
【文章內(nèi)容簡介】 %抽樣頻率 fs=100000。 t=[0:1/Fs:]。 %抽樣間隔 Fc=30000。 %載波頻率 a=1。 %設(shè) a=1 m=a*cos(300*2*pi*t)。 %調(diào)制信號(hào) X=fft(m)。 %對(duì) m 進(jìn)行傅里葉變換 X=abs(X(1:length(X)/2+1))。 %調(diào)制信號(hào)頻譜 frqX=(0:length(X)1)*Fs/length(X)/2 S = modulate(m,Fc,Fs,39。amssb39。)。 %對(duì)信號(hào)進(jìn)行調(diào)制 武漢理工大學(xué)《信號(hào)分析處理 》課程設(shè)計(jì)說明書 5 Y=fft(S)。 %對(duì) S 進(jìn) 行傅里葉變換 Y=abs(Y(1:length(Y)/2+1))。 frqY=(0:length(Y)1)*Fs/length(Y)/2 。 %已調(diào)信號(hào)頻譜 set(gcf,39。color39。,39。w39。) %顏色設(shè)置 figure(1) subplot(421) %繪制調(diào)制信號(hào)時(shí)域圖 plot(t,m) title(39。調(diào)制信號(hào) 時(shí)域圖 39。) subplot(422) %繪制調(diào)制信號(hào)頻譜圖 plot(frqX,X) axis([0 3000 0 max(X)]) title(39。調(diào)制信號(hào)頻譜 圖 39。) subplot(423) %繪制 已調(diào)信號(hào) 時(shí)域圖 plot(t,S) title(39。已調(diào)信號(hào) 時(shí)域圖 39。) subplot(424) plot(frqY,Y) %繪制 已調(diào)信號(hào)頻譜 圖 axis([0 60000 0 max(Y)]) title(39。已調(diào)信號(hào)頻譜 圖 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。)。 %無噪聲已調(diào)信號(hào)解調(diào) YYN=demod(sn,Fc,Fs,39。amssb39。)。 %加噪聲已調(diào)信號(hào)解調(diào) YYN1=demod(sn1,Fc,Fs,39。amssb39。)。 YYN2=demod(sn2,Fc,Fs,39。amssb39。)。 YYN3=demod(sn3,Fc,Fs,39。amssb39。)。 武漢理工大學(xué)《信號(hào)分析處理 》課程設(shè)計(jì)說明書 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。 %加噪聲后已調(diào) 信號(hào)頻譜 J2=fft(YYN) J2=abs(J2(1:length(J2)/2+1)) frqJ2=(0:length(J2)1)*Fs/length(J2)/2 %加噪聲后解調(diào)信號(hào)頻譜 set(gcf,39。color39。,39。w39。) figure(1) subplot(425) %繪制 加噪聲解調(diào)信號(hào) 時(shí)域圖 plot(t,YYN)。 title(39。加噪聲解調(diào)信號(hào) 時(shí)域圖 39。) subplot(426) %繪制 加噪聲解調(diào)信號(hào)頻譜 圖 plot(frqJ2,J2)。 axis([0 3000 0 max(J2)]) title(39。加噪聲解調(diào)信號(hào)頻譜 圖 39。) subplot(427) %繪制 無噪聲解調(diào)信號(hào) 時(shí)域圖 plot(t,Y1) title(39。無噪聲解調(diào)信號(hào) 時(shí)域圖 39。) subplot(428) %繪制 無噪聲解調(diào)信號(hào)頻譜 圖 plot
研究報(bào)告相關(guān)推薦
鄂ICP備17016276號(hào)-1