【正文】 y=fft(st2)。 plot(t,y)。 title(39?；鶐盘柗创a的頻譜39。)。 xlabel(39。f(Hz)39。)。 ylabel(39。幅度39。)。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%載波信號 s1=cos(2*pi*f1*t)。 s2=cos(2*pi*f2*t)。 figure(2) subplot(221) plot(s1)。 title(39。載波信號139。)。 xlabel(39。t39。)。 ylabel(39。幅度39。)。 subplot(222) y=fft(s1)。 plot(t,y) title(39。s1的頻譜圖39。)。 xlabel(39。f(Hz)39。)。 ylabel(39。幅度39。)。 subplot(223) plot(s2)。 title(39。載波信號239。)。 xlabel(39。t39。)。 ylabel(39。幅度39。)。 subplot(224) y=fft(s2)。 plot(t,y)。 title(39。s2的頻譜圖39。)。 xlabel(39。f(Hz)39。)。 ylabel(39。幅度39。)。 figure(3) fsk1=st1.*s1。 fsk2=st2.*s2。 fsk=fsk1+fsk2。 no=.*randn(1,t1)。 sn=fsk+no。 subplot(411) plot(t,fsk)。 title(39。fsk波形39。)。 ylabel(39。幅度39。) subplot(412) plot(t,sn)。 title(39。加入噪聲之后的波形39。)。 xlabel(39。t39。)。 ylabel(39。幅度39。)。 subplot(413) y=fft(fsk)。 plot(t,y)。 title(39。fsk的頻譜圖39。)。 xlabel(39。f(Hz)39。)。 ylabel(39。幅度39。)。 subplot(414) y=fft(sn)。 plot(t,y)。 title(39。sn的頻譜圖39。)。 xlabel(39。f(Hz)39。)。 ylabel(39。幅度39。)。 figure(4) b1=fir1(80,[15300 17300])。 b2=fir1(80,[18300 90300]) %設置帶寬參數 H1=filter(b1,1,sn)。 %b1為分子，1為分母，sn為濾波器輸入序列 H2=filter(b2,1,sn)。 %噪聲信號同時通過兩個濾波器 subplot(221) plot(t,H1)。 title(39。經過帶通濾波器H1后的波形39。)。 xlabel(39。t39。)。 ylabel(39。幅度39。)。 subplot(222)。 plot(t,H2)。 title(39。經過帶通濾波器H2后的波形39。)。 xlabel(39。t39。)。 ylabel(39。幅度39。)。 sw1=H1.*H1。 sw2=H2.*H2。 subplot(223) plot(t,sw1)。 title(39。經過乘法器一后的波形39。)。 xlabel(39。t39。)。 ylabel(39。幅度39。)。 subplot(224) plot(t,sw2)。 title(39。經過乘法器二的波形39。)。 xlabel(39。t39。)。 ylabel(39。幅度39。)。 bn=fir1(8,[1300 9300])。 %經過低通濾波器 figure(5) sr1=filter(bn,1,sw1)。 sr2=filter(bn,1,sw2)。 subplot(221) plot(t,sr1) title(39。經過低通濾波器后的波形39。)。 xlabel(39。t39。)。 ylabel(39。幅度39。)。 subplot(222) plot(t,sr2)。 title(39。經過低通濾波器后的波形39。)。 xlabel(39。t39。)。 ylabel(39。幅度39。)。 %判決for i=1:length(t) if(sr1(i)=sr2(i)) sr(i)=1。 else sr(i)=0。 end endsr=sr1+sr2。subplot(223)plot(t,sr)。title(39。經過抽樣判決器的波形39。)。xlabel(39。t39。)。ylabel(39。幅度39。)。subplot(224)a=SNR_dB=6::20。SNR=10.^()。SNR2=a.^2.(2*SNR)。fsk_pe0=*erfc(sqrt(a.^2.(4*SNR2)))。semilogy(SNR_dB,fsk_pe0,39。39。)。legend(39。相干2FSK39。)title(39。相干解調的誤碼率39。)。xlabel(39。SNR_dB39。)。ylabel(39。Pe39。)。參考文獻[1] 曹志剛. 現代通信原理. 清華大學出版社. [2] 劉樹棠. 信號與系統(tǒng). 西安交通大學出版社. [3] 劉衛(wèi)國. MATLAB程序設計教程. 中國水利水電出版社. [4] 樓順天. .[5] 梁虹、梁潔、.[6] 樊昌信,2006