【正文】 。加入噪聲三角波解調(diào)信號時域圖 39。)。 plot(t,yyn_tri) xlabel(39。)。 title(39。time39。 plot(t,yn_tri) axis([0 ])。)。 title(39。time39。 %***********三角波時域圖 *********************************** figure(9) subplot(3,1,1)。加入 噪聲正弦波解調(diào)信號時域圖 39。)。 plot(t,yyn_sin) %plot(t,yyn_sin1) xlabel(39。)。 title(39。time39。 plot(t,yn_sin) plot(t,yn_sin1) axis([0 ])。)。 title(39。time39。) %***********正弦波時域圖 *********************************** figure(8) subplot(3,1,1)。) ylabel(39。) hold on plot(in,out) xlabel(39。 %**解調(diào)器輸入輸出信噪比關(guān)系 ********************************* figure(7) plot(in,out,39。三角波已調(diào)信號頻譜圖 39。)。 xlabel(39。 figure(6) subplot(2,1,2)。三角波調(diào)制信號頻譜圖 39。)。 xlabel(39。 %*****************顯示三角波調(diào)制解調(diào)過程頻譜圖 ************************ figure(6) subplot(2,1,1)。正弦波已調(diào)信號頻譜圖 39。)。 xlabel(39。 figure(5) subplot(2,1,2)。正弦波調(diào)制信號頻譜圖 39。)。 xlabel(39。 %*****************顯示正弦波調(diào)制解調(diào)過程頻譜圖 ****************************** figure(5) subplot(2,1,1)。加入高斯白噪聲三角波解調(diào)信號時域圖 39。)。 plot(t,yyn_tri) xlabel(39。)。 title(39。time39。 %*********顯示加入高斯白噪聲解調(diào)信號 ********************* figure(4) subplot(2,1,1)。)。 title(39。time39。 figure(3) subplot(2,1,2)。)。 《武漢理工大學 MATLAB 課程設計報告》 25 title(39。time39。 %***********加入高斯白噪聲已調(diào)波的顯示 *************************** figure(3) subplot(2,1,1)。)。 title(39。time39。 figure(2) subplot(2,1,2)。)。 title(39。time39。 %***********已調(diào)波的顯示 *************************** figure(2) subplot(2,1,1)。載波時域圖 39。)。 plot(t,ct) xlabel(39。)。 title(39。time39。 subplot(3,1,2)。正弦波調(diào)制信號時域圖 39。)。 plot(t,mt_sin) 《武漢理工大學 MATLAB 課程設計報告》 24 xlabel(39。 out=[snr2,snr21,snr22,snr23,snr24]。 % snr24=var(y_sin)/var(dy24)。 dy24=yyn_sin4y_sin。 snr14=var(sfm_sin)/var(dy14)/Fs/10。 snr23=var(y_sin)/var(dy23)/Fs*400。 dy23=yyn_sin3y_sin。 %snr13=var(sfm_sin)/var(dy13)。 snr22=var(y_sin)/var(dy22)/Fs*400。 dy22=yyn_sin2y_sin。 %snr12=var(sfm_sin)/var(dy12)。 《武漢理工大學 MATLAB 課程設計報告》 23 %snr21=var(y_sin)/var(dy21)。 dy21=yyn_sin1y_sin。 snr11=var(sfm_sin)/var(dy11)/Fs/10。 %snr2=var(y_sin)/var(dy2)。 dy2=yyn_siny_sin。 snr1=var(sfm_sin)/var(dy1)/Fs/10。)。 yyn_tri4=demod(yn_tri4,30000,Fs,39。fm39。)。 yyn_tri2=demod(yn_tri2,30000,Fs,39。fm39。)。 yyn_sin4=demod(yn_sin4,30000,Fs,39。fm39。)。 yyn_sin2=demod(yn_sin2,30000,Fs,39。fm39。*Fs/length(ff4)。 mag4=abs(ff4)。*Fs/length(ff3)。 《武漢理工大學 MATLAB 課程設計報告》 22 mag3=abs(ff3)。*Fs/length(ff2)。 mag2=abs(ff2)。*Fs/length(ff1)。 mag1=abs(ff1)。)。 yyn_tri=demod(yn_tri,fc,Fs,39。fm39。)。 y_tri=demod(sfm_tri,fc,Fs,39。fm39。 yn_tri4=awgn(sfm_sin,40)。 yn_tri3=awgn(sfm_tri,80)。 yn_tri2=awgn(sfm_tri,60)。 yn_tri1=awgn(sfm_tri,10)。 yn_tri=awgn(sfm_tri,50)。 sfm_tri=am*cos(2*pi*fc*t+mf*2*pi*fm*int_mt_tri)。 for i=1:length(t)1 %三角波積分 int_mt_tri(i+1)=int_mt_tri(i)+mt_tri(i)*dt。 %正弦波積分 for i=1:length(t)1 int_mt_sin(i+1)=int_mt_sin(i)+mt_sin(i)*dt。 %三角波 mf=3。 %調(diào)制信號頻率 mt_sin=am*cos(2*pi*fm*t)。 %載波信號 am=1。 %***************信號設定 *********************** fc=30000。 %抽樣頻率 Fs=100000hz,滿足抽樣定理 n=0:N1。 %積分步長 N=1000。 北京：中國水利水電出版社。北京：機械工業(yè)出版社。 北京：國防工業(yè)出版社。這次課程設計終于順利完成了，在設計中遇到了很多問題，最后 在不懈的努力下 終于 迎刃而解 。通過這次課程設計使我懂得了理論與實際相結(jié)合是很重要的，只有理論 知識是遠遠不夠的，只有把所學的理論知識與實踐相結(jié)合起來，從理論中得出結(jié)論，才能真正為社會服務，從而提高自己的實際動手能力和獨立思考的能力。結(jié)果表明：信道噪聲很大程度上影響著調(diào)制解調(diào)系統(tǒng)性能，實際運用中應盡量降低信道噪聲；為了少失真地得到原始信號，在提高系統(tǒng) 性能的基礎上，還應增大輸入信噪比。 （ 5）輸入輸出信噪比呈正向線性關(guān)系，輸入信噪比越大，輸出信噪比也越大。 （ 4）在加入較強的高斯白噪聲時，調(diào)制系統(tǒng)輸出不太理想的 FM 已調(diào)信號，已調(diào)信號的頻率與原信號的幅度大小關(guān)系不明顯。 （ 3）在加入較弱的高斯白噪聲時，調(diào)制系統(tǒng)輸出較為理 想的 FM 已調(diào)信號，已調(diào)信號的頻率隨著原信號的幅度大小而變但并不是理想的。 （ 2）在未加入高斯白噪聲時，調(diào)制系統(tǒng)輸出理想的 FM 已調(diào)信號，已調(diào)信號的頻率隨著原信號的幅度大小而變。這里針對部分仿真結(jié)果進行分析。 圖 8 程序流程圖 MATLAB 源代碼 根據(jù)設計要求所編寫的源代碼已在在附頁中給出，詳情請看附頁 I?？梢?，加大調(diào)頻指數(shù) fm ，可使調(diào)頻系統(tǒng)的抗噪聲性能迅速改善。 設調(diào)頻信號為 ( ) c o s( ( ) )tF M c fS t A t K m d? ? ????? ? 故其輸入功率為 22i AS ? 輸入噪聲功率為 i o FMN n B? 因此輸入信噪比為 2