【正文】 tle(39。) figure,subplot(211),plot(R_Y),title(39。) figure,subplot(211),plot(Sig),title(39。 end if i==0 figure,subplot(211),plot(Num),title(39。 %加上高斯白噪聲 %解調(diào) RYc=R_Y.*cos(2*pi*Fc/Fs*n)。 %濾波 y=x.*cos(2*pi*Fc/Fs*n)。 Num=rcosfir(Alfa,N_Filter,Rate,Tc,39。 %升余弦濾波器參數(shù) Tc=1/sigv。sigv=16。) % %子函數(shù) runup，完成調(diào)制與解調(diào)，返回解調(diào)信號 function R_X=runup(Sig,SNR,i) NSig=length(Sig)。) ylabel(39。 %畫出誤碼率曲線 gtext(39。2psk39。 plot(t,pe(1,:),39。) axis([0,100,]) subplot(414),plot(abs(fft(x1))),title(39。) axis([0,100,]) subplot(412),plot(abs(fft(O_Y))),title(39。 if O_a(n)==1 O_Y(n)=1。O_a(1)=xor(R_X(1),1)。 end end R_X=runup(Sig,SNR,i)。b(1)=1。) axis([0,100,inf,inf]) else i=1。) axis([0,100,inf,inf]) subplot(413),plot(x1),title(39。 end end if SNR==0 figure,subplot(411),plot(R_X),title(39。%2PSK 信號 else Sig(n)=1。原信號波形頻譜 39。%產(chǎn)生隨機信源，信源符號 1和 0的概率分別為 ， if SNR==0 i=0。 %用于存儲不同信噪比不同調(diào)制方式下的錯誤符號個數(shù) pe=zeros(2,21)。 仿真參數(shù)設(shè)計 抽樣速率 fs=256Hz 載波速率 fc=64Hz 碼元速率 sigv=16B 碼元周期 Tc=1/sigv。)。 %向功率為 0dB 的信號 x 添加高斯白噪聲，輸出信號 y的信噪比 SNR為參數(shù) snr (單位 dB)。sqrt39。%計算誤碼率 Xor : b(n+1)=xor(x1(n),b(n))。 圖 310 2DPSK解調(diào)低通濾波信號 SNR=0dB時及 2DPSK 解調(diào)濾波信號頻譜 SNR=0dB 時 圖 311 a圖 為 2DPSK抽樣判決后信號波形 SNR=0dB b圖為 2DPSK抽樣判決后信號頻譜 SNR=0dB c圖為原信號波形 d圖為 原信號頻譜 解調(diào)時加有噪聲的 2DPSK信號與同步載波在頻域相乘后經(jīng)過低通濾波后，再經(jīng)過抽樣判決得到相對碼，然后進行碼反變換得到絕對碼，從而恢復(fù)出原始數(shù)字信號，與源信號相比 是一致的，說明沒有產(chǎn)生誤碼，各時間波形和頻譜如圖 311所示。 圖 37 2DPSK信號和 2DPSK信號頻譜 源信號是隨機產(chǎn)生二進制碼元，載波頻為 fc=256Hz，幅度為１的余弦波，源信號經(jīng)過碼變換（差分編碼）產(chǎn)生碼變換后信號，各時間波形和頻譜如圖 38所示。 圖 35 2PSK 解調(diào)低通濾波后信號及 2PSK解調(diào)濾波后信號頻譜 圖 36 a圖為 2PSK抽樣判決后信號波形 b圖為 抽樣判決后信號頻譜， c圖為原信號波形 d圖為 原信號頻譜 解調(diào)時加噪聲的 2PSK 信號與同步載波在頻域相乘再經(jīng)過低通濾波后，再經(jīng) 過抽樣判決后恢復(fù)出原始數(shù)字信號，各時間波形和頻譜如圖 36所示，抽樣判決 后 信號 時域波形 和頻譜與 源信號相比 是一致的，說明沒有產(chǎn)生誤碼。 圖２－ 6 2DPSK 差分相干解調(diào)器原理框圖 e2DPSK (t) 相乘器 低通濾波器 抽樣判決器 定時脈沖 輸出 帶通濾波器 延遲 Ts 圖 31原信號波形及原信號波形頻譜 其中源信號是隨機產(chǎn)生的二進制碼元，載波頻為 fc=256Hz，幅度為１的余弦波。 2DPSK系統(tǒng)是一種實用的數(shù)字調(diào)相系統(tǒng)，但其抗加性白噪聲性能比 2PSK的要差。 2DPSK信 號的另一種解調(diào)方法是差分相干解調(diào)（相位比較法），其原理框圖如圖２－ 6所示。 式（２ .3－ 2）稱為差分編碼（碼變換），即把絕對碼變換為相對碼；其逆過程稱為差分譯碼（碼反變換），即 1??? nnn bba （ － 3） 2DPSK信號的解調(diào)方法之一是相干解調(diào)（極性比較法）加碼反變換法。 2DPSK信號的產(chǎn)生方法：先對二進制數(shù)字基帶信號進行差分編碼，即把表示數(shù)字信息的序列的絕對碼變換成相對碼（差分碼），然后再根據(jù)相對碼絕對調(diào)相，從而產(chǎn)生二進制差分相移鍵控信號。因此，在相鄰碼元之間必定有相位突跳。這樣， 2DPSK方式雖然解決了載波相位不確定性問題，但是碼元的定時問題仍沒有解決。 為了更直觀地說明信號碼元的相位關(guān)系，我們可以用矢量圖 來表述。為了克服此缺點，提出了二進制差分相移鍵控（ 2DPSK）方式。在相干解調(diào)中，如何得到與接收的 2PSK信號同頻同相的相干載波是個關(guān)鍵問題。這種以載波的不同相位直接去表示相應(yīng)二進制數(shù)字信號的調(diào)制方式，稱為二進制絕對相移方式。這樣，可以用載波的某些離散狀態(tài)來表示數(shù)字基帶信號的離散狀態(tài)。然而，在實際信道中，大多數(shù)信道具有帶通傳輸特性，數(shù)字基帶信號不能直接在這種帶通傳輸特性的信道中傳輸。 計算機仿真軟件在通信系統(tǒng)工程設(shè)計中發(fā)揮著越來越重要的作用。因此，大部分現(xiàn)代通信系統(tǒng)都使用數(shù)字調(diào)制技術(shù)。通過仿真更深刻地理解了數(shù)字調(diào)制解調(diào)系統(tǒng)基本原理 。通過改變載波幅度、頻率、 相位，來傳輸數(shù)字基帶信號，所以帶通傳 輸也 叫做載波傳輸。 目前， MATLAB 集科學(xué)計算(putation) 、可視化 (visualization)、編程 (programming)于一身，并提供了豐富的 Windows 圖形界面設(shè)計方法。最后， 對兩 種調(diào)制解調(diào)系統(tǒng)的性能進行了 比較。 數(shù)字調(diào)制是通信系統(tǒng)中最為重要的環(huán)節(jié)之一，數(shù)字調(diào)制技術(shù)的改進也是通信系統(tǒng)性能提高的重要途徑。利用 MATLAB作為編程工具，設(shè)計了相移鍵控系統(tǒng)的模型，并且對模型的方針流程以及仿真結(jié)果都給出具體詳實的分析，為實際系統(tǒng)的構(gòu)建提供了很好的依據(jù)。通過仿真更深刻地理解了數(shù)字調(diào)制解調(diào)系統(tǒng)基本原理 。 MATLAB 的基本數(shù)據(jù)單位是矩陣，它的指令表達(dá)式與數(shù)學(xué)、工程中常用的形式十分相似，雇傭 MATLAB 可以進行矩陣、控制設(shè)計、信號處理與通信、圖像處理、信號檢測等領(lǐng)域。 在接收端通過解調(diào)器把帶通信號還原成數(shù)字基帶信號的過程稱為數(shù)字解調(diào) ，而包括調(diào)制和解調(diào)的過程數(shù)字傳輸 系統(tǒng) 叫做數(shù)字帶通傳輸系統(tǒng)。通過仿真， 分析 了 這 兩 種 調(diào)制解調(diào)過程中各環(huán)節(jié)時域和頻域的波形 ，并考慮了信道噪聲的影響。為了進行長途傳輸，必須對數(shù)字信號進行載波調(diào)制，將信號頻譜搬移到高頻處才能在信道中傳輸。通過對調(diào)制系統(tǒng)的仿真，我 們可以更加直觀的了解數(shù)字調(diào)制系統(tǒng)的性能及影響性能的因素，從而便于改進系統(tǒng)，獲得更佳的傳輸性能。 相關(guān)理論 數(shù)字調(diào)制與解調(diào)的相關(guān)原理 在數(shù)字基帶傳輸系統(tǒng)中，為了使 數(shù)字基帶信號能夠在信道中傳輸，要求信道應(yīng)具有低通形式的傳輸特性。但是，數(shù)字基帶信號具有與模擬基帶信號不同的特點，其取值是有限的離散狀態(tài)。因此， 2P