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