【正文】 。程序的界面可以做得更加人性化，這樣使用戶使用更加方便，畫圖函數(shù)也還可以改進(jìn)，例如：在使用 figure 分圖顯示的時候，我盡量選擇部分圖，都在一個 figure 中顯示，這樣既符合了，老師的要求 —— 盡量在一個圖中完成設(shè)計，同時也使整個設(shè)計的顯得干凈利落。 （四）要學(xué)會總結(jié)，提高。大學(xué)中的老師是每個大學(xué)生的好資源，他的博學(xué)、 重教，向他們請教會有很大收獲。 在編程過程中，我遇到了許多問題，通過查找資料也未能找到原因，每次我想老師請教，老師總能不厭其煩地一點點講解。由于 課堂時間有限， MATLAB 課程只能靠我們自學(xué)，然而 guide 畫圖功能我們必須用到，因此只能通過自習(xí)，實踐也證明我們是有自學(xué)的潛能的，通過自學(xué)課本，不理解的知識，通過查找資料，請教老師，最終掌握知識，因此，這次課程設(shè)計時間培養(yǎng)了我們的自學(xué)能力。 （二）我們在大學(xué)需要學(xué)習(xí)的不僅是基礎(chǔ)知識、專業(yè)知識，更重要的是一種學(xué)習(xí)能力。經(jīng)過一次又一次調(diào)試、修改又修改，一點一點發(fā)現(xiàn)問題并改正，我才真正 發(fā)現(xiàn)編程遠(yuǎn)沒有想象中的簡單。 當(dāng)然成 功編寫一個程序絕非易事，之前，我總以為程序能夠正常運行，就代表著編程成功，后來我才發(fā)現(xiàn)我大錯特錯了。這是由于自己對 MATLAB的模塊設(shè)計不夠理解，同時對 MATLAB 的基本語句一無所知，不過通過請教老師懂得了首先要設(shè)計 guide 目錄，再根據(jù)用戶輸入 ,執(zhí)行 swtich 語句，在每個 case后調(diào)用一個函數(shù)，來實現(xiàn)要求的功能，這樣一下子豁然開朗，掌握了基本設(shè)計思路之后，后面的編程就順利多了。 以下是我的幾點切身感受： （一）編寫程序需要一個清醒的頭腦，明確的思路，同時也要有耐心毅力。 通過這兩周的通信原理課程設(shè)計實踐，我復(fù)習(xí)了 MATLAB 編程語言的基本概念、 語法、語義和數(shù)據(jù)類型的使用特點，加深了對課堂所學(xué)理論知識的理解，掌握了運用結(jié)構(gòu)化程序設(shè)計的基本思想和方法，更重要的是培養(yǎng)了自己的自學(xué)能力。)。title(39。 subplot(3, 1, 3)。MFSK后的波形 39。 plot(t, s)。 axis([0, 1, , ])。發(fā)送端二進(jìn)制波形 39。 plot(t, m)。 Pss = S.* conj(S) / 512。 %Power Spectrum f = 1000*(0:256)/512。 m = 2*a(ceil(10*t+))3。 %Source amp。 M = 8 df = 10。 MFSK 的調(diào)制與解調(diào) %Time t = 0 : 1/1e3 : 。功率譜密度圖 39。 plot(f, Pss(1:257))。)。title(39。 subplot(3, 1, 2)。)。title(39。 subplot(3, 1, 1)。 S = fft(s,512)。 s = cos(2*pi*100*t+m/M)。 m = randint(10, 1, M)39。 MPSK 的調(diào)制與解調(diào) M = 8。MASK后的波形 39。 plot(t, s)。)。title(39。 subplot(2, 1, 1)。 g = [g, g, g, g, g, g, g, g, g, g]。 a = (2*a 1 M)*d。 t = 0 : 1/1e3 : 。 （ 2） M=8； MASK 的調(diào)制與解調(diào) M = 8。功率譜密度圖 39。 plot(f, Pss(1:257))。)。title(39。 subplot(3, 1, 2)。)。title(39。 %Drawing subplot(3, 1, 1)。 S = fft(s,512)。 s = cos(2*pi*(100+m*df).*t)。 4FSK a = randint(1, 10, M)。 M = 4。 %Frequency offset amp。)。title(39。 subplot(3, 1, 3)。MPSK后的波形 39。 plot(t, s)。 axis([0 1 min(m) max(m)+])。發(fā)送端二進(jìn)制波形 39。 plot(t, m)。 Pss = S.* conj(S) / 512。 f = 1000*(0:256)/512。 m = m(ceil(10*t+))。 t = 0::。)。title(39。 subplot(2, 1, 2)。發(fā)送的二進(jìn)制波形 39。 plot(t, a(ceil(10*t+)))。 s = a(ceil(10*t+)).*g.*cos(2*pi*100*t)。 g = ones(1, 100)。 a = randint(1, 10, M)。 d = 1。 （ 2）輸出提示語句未顯示 錯誤：提示語與程序中設(shè)置顯示語言發(fā)生沖突。 七 、調(diào)試分析 調(diào)試方法：運行程序，按要求測試各功能能否實現(xiàn)，如果不能實現(xiàn)，則尋找問題原因，修改程序后再次運行，再次測試功能能否實現(xiàn)，直到成功實現(xiàn)所有功能。 grid 六 、 調(diào)試過程及結(jié)論 基于 MATLAB 的 2ASK 調(diào)制解調(diào)仿真過程及結(jié)論 仿真過程中的各點波形 圖 7 發(fā)送端二進(jìn)制波形 圖 8 已調(diào) 2ASK 波形 圖 9 2ASK 頻譜分析 圖 10 加入 高斯白噪聲后的 2ASK 波形 圖 11 解調(diào)過程中各點的輸出波形 圖 13 解調(diào)后的波形 基于 MATLAB 的 2FSK 調(diào)制解調(diào)仿真過程及結(jié)論 圖 14 發(fā)送端二進(jìn)制波形 圖 15 已調(diào) 2FSK 波形 圖 16 FSK 仿真頻譜分析圖 圖 17 圖 18 圖 19 解調(diào)后的波形 基于 MATLAB 的 2PSK 調(diào)制解調(diào)仿真過程及結(jié)論 圖 20 接收端二進(jìn)制波形 圖 21 2PSK 調(diào)制后的波形 圖 22 2PSK 仿真頻譜分析圖 圖 23 加入 高斯白噪聲后的 2PSK 波形 圖 24 圖 25 接收端二進(jìn)制波形 基于 MATLAB 的 2DPSK 調(diào)制解調(diào)仿真過程及結(jié)論 圖 26 圖 27 圖 28 基于 MATLAB 的 HDB3 碼的編碼與譯碼 仿真過程及結(jié)論 調(diào)試過程及結(jié)論 按照設(shè)計的調(diào)制解調(diào)框圖，根據(jù)每部分在整個系統(tǒng)的作用寫好程序，運行尋找錯誤并修改，最后得到可以滿足本設(shè)計的程序。逆碼變換后的波形 39。 axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), 3*(A / 2), 3*(A / 2) ] )。 else sign_result( (I1)*Num_Unit + 1 : I*Num_Unit) = Low_Level。 )。 title ( 39。 subplot(2,1,2) plot ( t, dt )。 [bz,az]=impinvar(b,a,fs)。s39。)。 grid [N,Wn] = buttord( 2*pi*50, 2*pi*150,3,25,39。相乘后的波形 39。 axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), 3*(A / 2), 3*(A / 2) ] )。 grid dt = st .* cos ( 2 * pi * w * t )。調(diào)制后的信號 39。 axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), 3*(A / 2), 3*(A / 2) ] )。 grid subplot ( 2, 1, 2 )。原始信號 39。 axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), (A / 2), A + (A / 2) ] )。 else st( (I1)*Num_Unit + 1 : I*Num_Unit) = A * cos ( 2 * pi * w * t( (I1)*Num_Unit + 1 : I*Num_Unit ) )。 else sign_orign( (I1)*Num_Unit + 1 : I*Num_Unit) = Low_Level。 t = 0 : 1/fs : Time_Hold_On * Lenth_Of_Sign 1/fs。 sign_orign = zeros ( 1, Num_Unit * Lenth_Of_Sign )。 Sign_Set = [0,1,1,0,1,0,0,1] Lenth_Of_Sign = length ( Sign_Set )。 w = 300。 High_Level = ones ( 1, Num_Unit )。 Time_Hold_On = 。)。 ylabel(39。r/dB39。)。 title(39。 semilogy(startn:length(fpepsk)+startn1,fpepsk)。 P=CheckRatePe(signal,configueSignal,n) fpepsk=[fpepsk,P]。 configueSignal=demoBPSK(signal1,bitRate,fc,n,N)。 transmitted