【正文】 16 。 。,39。 legend 39。*39。 %Eb/N0 in ratio thBER *erfc sqrt EbN0_R/2 。 。 ylabel 39。Eb/N0 in dB39。 。%和 while 連用 end end end semilogy myEbN0,myBER,39。 myBER index tempBER myEbN0 index EbN0 index index+1。%根據(jù)計算出來的誤碼率來計算實際需要的測試次數(shù)下限（ thresholdTestLength） else thresholdTestLength TestLength。 %控制 while 循環(huán) end %接下來判斷要不要繼續(xù)測試下去 tempBER errorCount/testCount。 end testCount testCount+1。ave2 1/2 errorCount errorCount+1。ave2 1/2 errorCount errorCount+1。%抽樣值的絕對值的平均值，大于1/2 為‘ 1’，小于 1/2 為‘ 0’ ave2 sum abs efsk22 /length mm 。%低通濾波 efsk22 filtfilt ones 1,10 ,efsk2 。%頻率為 w1 的相干解調(diào) efsk2 Noise_fsk.*cos w2*mm 。%高斯白噪聲矩陣 Noise_fsk cos w1*mm *src1+cos w2*mm *src2+noise 。 src2 1src。 %抽樣判決在 t 軸上的抽樣矩陣 while 1 % 以下為源信號經(jīng)過 FSK 調(diào)制解調(diào) src rand 1 。 %信號‘ 1’載波角頻率 w2 2021 。 %誤碼的符號個數(shù) testCount 0。 %noise density 噪聲密度 noiseRoot sqrt noiseVar 。 N0 Eb*10^ EbN0/10 。 %作為每個信噪比的下標(biāo) for EbN0 0:1:10 %信噪比由 1 到十 disp 39。 。%signal amplitude Eb A*A。 % test length should at least reach this number TestLength 50000。,sum 100*clock 。 %Reset the generator of uniformly distributed random numbers randn 39。state39。 clear all。PropD 0。 %濾除高斯白噪聲 Delay 3。 title 39。 [errorSym ratioSym] symerr x,z2 。 title 39。 %誤碼率統(tǒng)計 [errorSym ratioSym] symerr x,z1 。Amplitude39。 。 xlabel 39。,39。 legend 39。 title 39。 hold off。ro39。 hold on。bx39。非相干解調(diào)后的信號的眼圖 39。,M,df 。 %非相干解調(diào) figure 6 z2 ddemod y,Fc,Fd,Fs,39。Amplitude39。 。 xlabel 39。,39。 legend 39。 title 39。 hold off。ro39。 hold on。bx39。相干解調(diào)后的信號的眼圖 39。,M,df 。 %相干解調(diào) figure 4 z1 ddemod y,Fc,Fd,Fs,39。Amplitude39。 。 xlabel 39。 title 39。 。%在已調(diào)信號中加入高斯白噪聲 figure 3 plot t,y 1:length t ,39。dB39。measured39。,seed 2 。 %在已調(diào)信號中加入高斯白噪聲 randn 39。Amplitude39。 。 xlabel 39。 title 39。 。 figure 2 plot t,y 1:length t ,39。 numModPlot numPlot*Fs。fsk39。 。 ylabel 39。Time39。二進(jìn)制隨機序列 39。 。 %產(chǎn)生 25 個二進(jìn)制隨機碼 figure 1 stem [0:numPlot1],x 1:numPlot ,39。 numPlot 15。 %信噪比 SNR SNRpBit/log2 M 。 %進(jìn)行仿真的信 息代碼個數(shù) M 2。 df 10。 %系統(tǒng)采樣頻率 Fd 1。 參考文獻(xiàn) [2] 吳翊，吳盂達(dá)，成禮智．?dāng)?shù)學(xué)建模的理論與實踐．北京：國防科技大學(xué)出版社， 1999 [3] 樊昌信 , 曹麗娜 ,通信原理 北京 :國防工業(yè)出版社 2021 [4] 孫祥，徐流美，吳清 .MATLAB ：清華大學(xué)出版社 2021. [5] 陳懷琛 吳大正 高西全 .MATLAB 及在電子信息課程中的應(yīng)用 [M].北京：電子工業(yè)出 ,2021 [6] 薛年喜 .MATLAB 在數(shù)字信號處理中的應(yīng)用 [M].清華大學(xué)出版社， 2021 附 錄 附件 1 function FSK Fc 10。其中，對MATLAB 軟件的學(xué)習(xí)也讓我感慨良多。但是由于時間比較倉促，加上本人所學(xué)知識和經(jīng)驗的不足，因此，在分析問題、解決問題時顯得不夠嚴(yán)密、完善，還需要在以后的學(xué)習(xí)工作中不斷地改進(jìn)和完善。 致謝 經(jīng)過兩個多月的構(gòu)思與設(shè)計，在指導(dǎo)老師肖永隆老師的精 心引導(dǎo)下，本論文已經(jīng)基本完畢。因此對其產(chǎn)生了更大的興趣，很有感覺。在將理論用各種方法實現(xiàn)的同時，我們也在不斷的搜集資料，不斷的學(xué)習(xí)，獲得更多的相關(guān)知識。實現(xiàn)了所學(xué) FSK 調(diào)制解調(diào)的仿真，對FSK 的原理更加熟悉了，體會到理論和實際是有好大不 同的，實踐離不開理論，理論只有應(yīng)用于實踐才能發(fā)揮其作用。 其三是對圖形函數(shù) stairs 的用法不時很清楚，總是有一個碼元是畫不出來的，可以看到在程序開始時， n 8，其實產(chǎn)生了 8 個碼元，但只畫出了 7 個碼元，用限定坐標(biāo)的方法解決了此問題，但仍需改進(jìn)。后來經(jīng)過參考各方面的資料 ，發(fā)現(xiàn)用 plot 就可以解決，而且效果很好，鑒于程序結(jié)構(gòu)不適合更改為 plot，并未應(yīng)用。當(dāng)然理論與實際還會有很大的出入，在設(shè)計時還要考慮各種干擾和噪聲等因素。在仿真過程中，構(gòu)建 FSK 的調(diào)制模型、相干和非相干解調(diào)模型、過零檢測和差分檢波模型，充分發(fā)揮了 SIMULINK 功能強大，建模簡單，參數(shù)易于調(diào)整的特點。并用 stairs 命令畫出了原始信號的波形圖 src1 與 src2。 src2 ones 1,n src。 end 而濾波在 MATLAB 中我們用 filtfilt 命令實現(xiàn) ,對已調(diào)信號進(jìn)行濾波后的碼元波形見圖 39 . 圖３ 9 已調(diào)信號經(jīng)低通濾波后時域圖 最后 ,經(jīng)過抽樣判決 ,可以得到解調(diào)后的數(shù)字信號序列 ,其波形見圖 310 . 圖 310 原始信號與解調(diào)信號對比圖 注：這里，我們用代碼 src rand 1,n P。%分支 2 抽樣值的絕對值的平均值， if ave1 ave2 final i 1。因此前者絕對值的的