第一篇：通信原理仿真

通信原理仿真實(shí)驗(yàn)提綱

1.任意產(chǎn)生一個(gè)調(diào)制信號(hào)，畫出其波形及其頻譜；

2.產(chǎn)生一個(gè)余弦載波信號(hào)，畫出其波形及其頻譜；

3.分別采用AM（幅度），DSB（雙邊），SSB（單邊）的方式對(duì)調(diào)

制信號(hào)進(jìn)行調(diào)制，畫出已調(diào)信號(hào)的波形及頻譜；

4.采用適當(dāng)?shù)姆绞?，分別對(duì)3中得到的已調(diào)信號(hào)進(jìn)行解調(diào)，畫

出解調(diào)信號(hào)的波形；

5.產(chǎn)生一個(gè)高斯白噪聲，疊加在已調(diào)信號(hào)上，然后進(jìn)行解調(diào)，畫出解調(diào)信號(hào)的波形；

6.比較4和5中的結(jié)果；

7.編寫A律13折線PCM編碼的程序，能夠?qū)θ我廨斎胄盘?hào)輸

出其PCM編碼；

8.產(chǎn)生一個(gè)隨機(jī)數(shù)字信號(hào)，分別進(jìn)行ASK，F(xiàn)SK，PSK調(diào)制解調(diào)，畫出解調(diào)前后的波形

第二篇：通信原理仿真作業(yè)2

通信原理仿真作業(yè)

第五章 模擬調(diào)制

1.AM、DSB調(diào)制及解調(diào)

用matlab產(chǎn)生一個(gè)頻率為1Hz，功率為1的余弦信源m(t)，設(shè)載波頻率?c?10Hz，m0?2,試畫出：

?

?

?

? AM及DSB調(diào)制信號(hào)的時(shí)域波形； 采用相干解調(diào)后的AM及DSB信號(hào)波形； AM及DSB已調(diào)信號(hào)的功率譜； 調(diào)整載波頻率及m0，觀察分的AM的過調(diào)與DSB反相點(diǎn)現(xiàn)象。? 在接收端帶通后加上窄帶高斯噪聲，單邊功率譜密度n0?0.1，重新解調(diào)。

2.SSB調(diào)制及解調(diào)

用matlab產(chǎn)生一個(gè)頻率為1Hz，功率為1的余弦信源，設(shè)載波頻率?c?10Hz，,試畫出：

? SSB調(diào)制信號(hào)的時(shí)域波形；

? 采用相干解調(diào)后的SSB信號(hào)波形；

? SSB已調(diào)信號(hào)的功率譜；

? 在接收端帶通后加上窄帶高斯噪聲，單邊功率譜密度n0?0.1，重新解調(diào)。

3.VSB調(diào)制及解調(diào)（參照作業(yè)題5－4）

用matlab產(chǎn)生一個(gè)頻率分別為5Hz、5/2 Hz的余弦和正弦疊加信號(hào)作為信源

m(t)，兩個(gè)頻率分量功率相同，總信號(hào)功率為2，設(shè)載波頻率為20Hz，試畫出：

? 殘留邊帶為0.2fm的VSB調(diào)制信號(hào)時(shí)域波形；

? 采用相干解調(diào)后的VSB信號(hào)波形；

? 調(diào)制信號(hào)的功率譜密度；

? 在接收端帶通后加上窄帶高斯噪聲，單邊功率譜密度n0?0.1，重新解調(diào)。

4.FM調(diào)制及解調(diào)

設(shè)輸入信號(hào)為m(t)?cos2?t，載波中心頻率為fc?10Hz，VCO的壓控振蕩系數(shù)為5Hz/V，載波平均功率為1W。試畫出：

? 已調(diào)信號(hào)的時(shí)域波形；

? 已調(diào)信號(hào)的振幅譜；

? 用鑒頻器解調(diào)該信號(hào)，并與輸入信號(hào)比較。

第三篇：通信原理課程設(shè)計(jì)----基帶信號(hào)眼圖仿真

需求分析

1．設(shè)基帶傳輸系統(tǒng)響應(yīng)是α＝1的升余弦滾降系統(tǒng)，畫出在接收端的基帶數(shù)字信號(hào)波形及其眼圖。

2．設(shè)定二進(jìn)制數(shù)字基帶信號(hào) an∈｛＋1,－1｝，g(t)= 1 0≤t≤Ts；t為其他值時(shí)g(t)= 0。系統(tǒng)加性高斯白噪聲的雙邊功率譜密度為0。畫出：(1)經(jīng)過理想低通H(f)= 1 │f│≤5/(2 Ts)后的眼圖。(2)經(jīng)過理想低通H(f)= 1 │f│≤1/ Ts后的眼圖。(3)比較分析上面圖形。

在該部分中敘述：對(duì)題目中要求的功能進(jìn)行的簡單的敘述分析，把題目內(nèi)容給介紹一下，還需要介紹分工情況。概要設(shè)計(jì)

1、基帶傳輸特性

基帶系統(tǒng)的分析模型如圖3-1所示，要獲得良好的基帶傳輸系統(tǒng)，就應(yīng)該

?a??t?nT?nsn基帶傳輸H(?)?ah?t?nT?nsn抽樣判決

圖2-1 基帶系統(tǒng)的分析模型

抑制碼間干擾。設(shè)輸入的基帶信號(hào)為?an??t?nTs?，Ts為基帶信號(hào)的碼元周期，n則經(jīng)過基帶傳輸系統(tǒng)后的輸出碼元為?anh?t?nTs?。其中

n1h(t)?2??????H(?)ej?td?

（3-1）

理論上要達(dá)到無碼間干擾，依照奈奎斯特第一準(zhǔn)則，基帶傳輸系統(tǒng)在時(shí)域應(yīng)滿足：

k?0?1，h(kTs)??0,k為其他整數(shù)?頻域應(yīng)滿足：

(3-2)

??T，???sTs

H(?)???0,其他??1

(3-3)

H(?)Ts??Ts0?Ts?

圖2-2 理想基帶傳輸特性

此時(shí)頻帶利用率為2Baud/Hz,這是在抽樣值無失真條件下，所能達(dá)到的最高頻率利用率。由于理想的低通濾波器不容易實(shí)現(xiàn)，而且時(shí)域波形的拖尾衰減太慢，因此在得不到嚴(yán)格定時(shí)時(shí)，碼間干擾就可能較大。在一般情況下，只要滿足：

???2?i?2??2??H?????H?????H(?)?H?????Ts,?TTTis?s?s???????Ts(3-4)基帶信號(hào)就可實(shí)現(xiàn)無碼間干擾傳輸。這種濾波器克服了拖尾太慢的問題。從實(shí)際的濾波器的實(shí)現(xiàn)來考慮，采用具有升余弦頻譜特性H(?)時(shí)是適宜的。

??Ts???(1??)?(1??)????1?sin?(??)?,Ts?TsTs?2????(1??)?

H(?)??1,0???Ts???(1??)0,???Ts??這里?稱為滾降系數(shù)，0???1。所對(duì)應(yīng)的其沖激響應(yīng)為：

sinh(t)?

(3-5)

?t

(3-6)

Tscos(??tTs)

?t1??4?2t2Ts2?Ts此時(shí)頻帶利用率降為2/(1??)Baud/Hz，這同樣是在抽樣值無失真條件下，所能 達(dá)到的最高頻率利用率。換言之，若輸入碼元速率Rs'?1/Ts，則該基帶傳輸系統(tǒng)輸出碼元會(huì)產(chǎn)生碼間干擾。

2、眼圖

所謂眼圖就是將接收濾波器輸出的，未經(jīng)再生的信號(hào)，用位定時(shí)以及倍數(shù)作為同步信號(hào)在示波器上重復(fù)掃描所顯示的波形（因傳輸二進(jìn)制信號(hào)時(shí)，類似人的眼睛）。干擾和失真所產(chǎn)生的畸變可以很清楚的從眼圖中看出。眼圖反映了系統(tǒng)的最佳抽樣時(shí)間，定時(shí)的靈敏度，噪音容限，信號(hào)幅度的畸變范圍以及判決門限電平，因此通常用眼圖來觀察基帶傳輸系統(tǒng)的好壞。

圖2-3

眼圖示意圖

部分的算法設(shè)計(jì)說明（可以是描述算法的流程圖）；每個(gè)程序中使用的存儲(chǔ)結(jié)構(gòu)設(shè)計(jì)說明（如果題目已經(jīng)指定了數(shù)據(jù)存儲(chǔ)的，按照指定的設(shè)計(jì)，并且寫出該存儲(chǔ)結(jié)構(gòu)的定義）。也就是簡單介紹一下你程序中自定義的函數(shù)（函數(shù)名，參數(shù)，以及功能）等。運(yùn)行環(huán)境

Windows 7 開發(fā)工具和編程語言

MATLAB V6.0 詳細(xì)設(shè)計(jì)

clear all;close all;N_eye = 7;%眼圖中碼元個(gè)數(shù) N0 = 0;%高斯白噪聲功率譜設(shè)為0 Ts = 2;%碼元寬度設(shè)為1 N_sample = 10;%每個(gè)碼元抽樣10個(gè)點(diǎn) dt = Ts/N_sample;gt = ones(1,N_sample);%矩形不歸零脈沖 N_code = 1000;%碼元個(gè)數(shù)

t = 0:dt:N_code*N_sample*dt-dt;data = sign(randn(1,N_code));data1 = sigexpand(data,N_sample);st = conv(data1,gt);%產(chǎn)生NRZ基帶信號(hào)

h1t =5/Ts*sinc(5*(t-5)/Ts);%濾波器1的沖激響應(yīng) h2t = 2/Ts*sinc(2*(t-5)/Ts);%濾波器2的沖激響應(yīng) h3t = 1/Ts*sinc(1*(t-5)/Ts);%濾波器3的沖激響應(yīng) h4t = 0.5/Ts*sinc(0.5*(t-5)/Ts);

r1t = conv(h1t,st);r2t = conv(h2t,st);r3t = conv(h3t,st);r4t = conv(h4t,st);

r1t = [0 r1t*dt];r2t = [0 r2t*dt];r3t = [0 r3t*dt];r4t = [0 r4t*dt];

sigma1 = sqrt(2.5*N0/Ts);sigma2 = sqrt(N0/Ts);sigma3 = sqrt(0.5*N0/Ts);sigma4 = sqrt(0.25*N0/Ts);

n1 = sigma1*randn(1,length(r1t));n2 = sigma2*randn(1,length(r2t));n3 = sigma2*randn(1,length(r3t));%加入白噪聲，該白噪聲是通過濾波器后的帶限白噪聲

n4 = sigma2*randn(1,length(r4t));

r1t = r1t + n1;r2t = r2t + n2;r3t = r3t + n3;r3t = r3t + n3;%至此得到基帶信號(hào)分別通過濾波器1和濾波器2 濾波器3濾波器4的輸出波形

ttt = 0:dt:N_eye*N_sample*dt-dt;figure(1)subplot(411)%畫出濾波器1輸出的眼圖 for n = 9:1:100 ss = r1t(n*N_sample*N_eye+1:(n+1)*N_sample*N_eye);drawnow;plot(ttt,ss);hold on;end ylabel('B=2.5Hz');subplot(412)%畫出濾波器2輸出的眼圖 for n = 9:1:100 ss = r2t(n*N_sample*N_eye+1:(n+1)*N_sample*N_eye);drawnow;plot(ttt,ss)hold on;end ylabel('B=1Hz')subplot(413)%畫出濾波器3輸出的眼圖 for n = 9:1:100 ss = r3t(n*N_sample*N_eye+1:(n+1)*N_sample*N_eye);drawnow;plot(ttt,ss)hold on;end ylabel('B=0.5Hz')subplot(414)%畫出濾波器4輸出的眼圖 for n = 9:1:100 ss = r4t(n*N_sample*N_eye+1:(n+1)*N_sample*N_eye);drawnow;plot(ttt,ss)hold on;end ylabel('B=0.25Hz')6 調(diào)試分析

剛開始的時(shí)候，程序出現(xiàn)了問題

??? Undefined command/function 'sigexpand'.之后添加了一個(gè)函數(shù)

function[out]=sigexpand(d,M)%將輸入的序列擴(kuò)展成間隔為N-1個(gè)0的序列 N=length(d);out=zeros(M,N);%例如，M=8，N=1000，產(chǎn)生一個(gè)8行1000列的矩陣 out(1,:)=d;%d的每一個(gè)元素賦值給out矩陣每列的第一個(gè)元素，因?yàn)閐中共1000個(gè)元素（0，-1或者1），且out矩陣共有1000列，所以當(dāng)然成立

out=reshape(out,1,M*N);%把out這個(gè)矩陣重新變成一個(gè)1行M*N列的矩陣，即1行8000列的有8000個(gè)元素的矩陣，其中每隔7個(gè)0就是一個(gè)原來的d中的元素

之后再調(diào)用sigexpand()函數(shù)的到時(shí)候就就沒有問題了。測試結(jié)果

程序運(yùn)行結(jié)果如圖7-1所示

圖7-1

程序運(yùn)行結(jié)果圖

比較上圖中的1和3可以發(fā)現(xiàn)，3的性能比1的性能。1中的失真比較嚴(yán)重。

參考文獻(xiàn)

[1] 孫屹，MATLAB通信仿真開發(fā)手冊(cè)，北京：國防工業(yè)出版社

[2] 李建新，現(xiàn)代通信系統(tǒng)分析與仿真－MATLAB通信工具箱，西安：西安電子科技大學(xué)出版社

[3] 曹志剛，現(xiàn)代通信原理，北京：清華大學(xué)出版社 [4] 教學(xué)用“通信原理”教材

心得體會(huì)

通過這次的課程設(shè)計(jì)，我更加理解了眼圖的用途，通過眼圖可以迅速的定性評(píng)價(jià)系統(tǒng)的性能?？梢杂^察到麻將干擾和信道噪聲等影響因素的情況蟲兒快速的判斷出系統(tǒng)性能的優(yōu)劣程度。最佳抽煙時(shí)刻是‘眼睛’張開最大的時(shí)刻，定時(shí)誤差靈敏度是眼圖斜邊的斜率，斜率越大，對(duì)應(yīng)的定時(shí)誤差越敏感。圖的陰影區(qū)的垂直高度便是抽樣時(shí)刻上信號(hào)受噪音干擾的畸變程度。圖中央的橫軸位置對(duì)應(yīng)于判決門限電平，抽樣時(shí)刻時(shí)上下兩陰影區(qū)的間隔距離之半為噪聲容限。此外，又學(xué)習(xí)了用matlab仿真通信原理實(shí)驗(yàn)的有關(guān)方法，及常用的相關(guān)函數(shù)。通過對(duì)程序錯(cuò)誤的調(diào)試，更加熟練的應(yīng)用軟件。通過仿真，更容易理解了通信原理的相關(guān)內(nèi)容。

通信原理是一門很有用的課程，但是內(nèi)容繁雜，需要記憶的很多，所以我們應(yīng)該加深的理解而不是單純的記憶，另外如果增加實(shí)踐的內(nèi)容，我們將學(xué)得更多，更快！路漫漫其修遠(yuǎn)兮，吾將上下而求索！

第四篇：通信原理課程設(shè)計(jì)_(基于MATLAB的_2PSK_2DPSK仿真)

江西農(nóng)業(yè)大學(xué)

通信原理課程設(shè)計(jì)報(bào)告

題 目 基于Matlab的相移鍵控仿真設(shè)計(jì)

專 業(yè) 電子信息工程

學(xué)生姓名 曾凡文

學(xué) 號(hào) 20121206

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告 2015年6月

基于Matlab的2PSK,2DPSK仿真

摘要：現(xiàn)代通信系統(tǒng)要求通信距離遠(yuǎn)、通信容量大、傳輸質(zhì)量好，作為其關(guān)鍵技術(shù)之一的調(diào)制技術(shù)

一直是研究的一個(gè)重要方向。本設(shè)計(jì)主要敘述了數(shù)字信號(hào)的調(diào)制方式，介紹了2PSK數(shù)字調(diào)制方式的

基本原理，功率譜密度，并運(yùn)用MATLAB軟件對(duì)數(shù)字調(diào)制方式2PSK進(jìn)行了編程仿真實(shí)現(xiàn)，在MATLAB平

臺(tái)上建立2PSK和2DPSK調(diào)制技術(shù)的仿真模型。進(jìn)一步學(xué)習(xí)了MATLAB編程軟件，將MATLAB與通信系統(tǒng)

中數(shù)字調(diào)制知識(shí)聯(lián)系起來，為以后在通信領(lǐng)域?qū)W習(xí)和研究打下了基礎(chǔ)在計(jì)算機(jī)上，運(yùn)用MATLAB軟件

來實(shí)現(xiàn)對(duì)數(shù)字信號(hào)調(diào)制技術(shù)的仿真。

關(guān)鍵詞：數(shù)字調(diào)制與解調(diào)；MATLAB；2PSK；2DPSK；

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

第1章 緒論

1.1 調(diào)制方式

數(shù)字通信系統(tǒng), 按調(diào)制方式可以分為基帶傳輸和帶通傳輸。數(shù)字基帶信號(hào)的功率一般處于從零開始到某一頻率（如0～6M）低頻段，因而在很多實(shí)際的通信（如無線信道）中就不能直接進(jìn)行傳輸，需要借助載波調(diào)制進(jìn)行頻譜搬移，將數(shù)字基帶信號(hào)變換成適合信道傳輸?shù)臄?shù)字頻帶信號(hào)進(jìn)行傳輸，這種傳輸方式，稱為數(shù)字信號(hào)的頻帶傳輸或調(diào)制傳輸、載波傳輸。所謂調(diào)制，是用基帶信號(hào)對(duì)載波波形的某參量進(jìn)行控制，使該參量隨基帶信號(hào)的規(guī)律變化從而攜帶消息。對(duì)數(shù)字信號(hào)進(jìn)行調(diào)制可以便于信號(hào)的傳輸；實(shí)現(xiàn)信道復(fù)用；改變信號(hào)占據(jù)的帶寬；改善系統(tǒng)的性能。

數(shù)字基帶通信系統(tǒng)中四種基本的調(diào)制方式分別稱為振幅鍵控（ASK，Amplitude-Shift keying）、移頻鍵控（FSK，F(xiàn)requency-Shift keying）、移相鍵控(PSK，Phase-Shift keying)和差分移相鍵（DPSK，Different Phase-Shift keying）。本次課程設(shè)計(jì)對(duì)PSK,DPSK這兩種調(diào)制方式進(jìn)行了仿真。

1.2 設(shè)計(jì)要求 1.2.1 設(shè)計(jì)內(nèi)容

用MATLAB完成對(duì)2PSK、2DPSK的調(diào)制與解調(diào)仿真電路設(shè)計(jì)，并對(duì)仿真結(jié)果進(jìn)行分析，可編寫程序，也可硬件設(shè)計(jì)框圖

1.2.2 設(shè)計(jì)參數(shù)（參數(shù)可以自行設(shè)置）

1、傳輸基帶數(shù)字信號(hào)（15位）碼元周期T=0.01S

2、載波頻率：15KHz 1.2.3 設(shè)計(jì)儀器

計(jì)算機(jī)和MATLAB軟件

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

第2章 2PSK,2DPSK原理

2.1 2PSK原理 2.1.1 2PSK基本原理

二進(jìn)制移相鍵控，簡記為2PSK或BPSK。2PSK信號(hào)碼元的“0”和“1”分別用兩個(gè)不同的初始相位“0”和“?”來表示，而其振幅和頻率保持不變.因此，2PSK信號(hào)的時(shí)域表達(dá)式為:

(t)=Acos其中，表示第n個(gè)符號(hào)的絕對(duì)相位：

t+)

=因此，上式可以改寫為:

這種以載波的不同相位直接表示相應(yīng)二進(jìn)制數(shù)字信號(hào)的調(diào)制方式,稱為二進(jìn)制移相鍵控方式。二進(jìn)制移相鍵控信號(hào)的典型時(shí)間波形如圖2-1。

10011tTs圖2-1 二進(jìn)制相移鍵控信號(hào)的時(shí)間波形

2.1.2 2PSK調(diào)制原理

在二進(jìn)制數(shù)字調(diào)制中，當(dāng)正弦載波的相位隨二進(jìn)制數(shù)字基帶信號(hào)離散變化時(shí)，則產(chǎn)生二進(jìn)制移相鍵控(2PSK)信號(hào)。2PSK信號(hào)調(diào)制有兩種方法，即模擬調(diào)制法和鍵控法。通常用已調(diào)信號(hào)載波的 0°和 180°分別表示二進(jìn)制數(shù)字基帶信號(hào)的 1 和 0，模擬調(diào)

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

制法用兩個(gè)反相的載波信號(hào)進(jìn)行調(diào)制。2PSK以載波的相位變化作為參考基準(zhǔn)的，當(dāng)基帶信號(hào)為0時(shí)相位相對(duì)于初始相位為0°，當(dāng)基帶信號(hào)為1時(shí)相對(duì)于初始相位為180°。

鍵控法，是用載波的相位來攜帶二進(jìn)制信息的調(diào)制方式。通常用0°和180°來分別代表0和1。其時(shí)域表達(dá)式為：

??e2PSK???ang(t?nTs)?cos?ct

?n?其中，2PSK的調(diào)制中an必須為雙極性碼。兩種方法原理圖分別如圖2-2和圖2-3所示。

圖2-2 模擬調(diào)制

原理圖

圖 2-3 鍵控法原理

圖

2.1.3 2PSK解調(diào)原理

由于2PSK的幅度是恒定的，必須進(jìn)行相干解調(diào)。經(jīng)過帶通濾波的信號(hào)在相乘器中與本地載波相乘，然后用低通濾波器濾除高頻分量，在進(jìn)行抽樣判決。判決器是按極性來判決的。即正抽樣值判為1，負(fù)抽樣值判為0。2PSK信號(hào)的相干解調(diào)原理圖如圖2-4所示，各點(diǎn)的波形如圖2-5所示。

由于2PSK信號(hào)的載波回復(fù)過程中存在著180°的相位模糊，即恢復(fù)的本地載波與所需 3

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

相干載波可能相同，也可能相反，這種相位關(guān)系的不確定性將會(huì)造成解調(diào)出的數(shù)字基帶信號(hào)與發(fā)送的基帶信號(hào)正好相反，即“1”變成“0”嗎“0”變成“1”，判決器輸出數(shù)字信號(hào)全部出錯(cuò)。這種現(xiàn)象稱為2PSK方式的“倒π”現(xiàn)象或“反相工作”。

e2PSK(t)帶通濾波器a相乘器c低通濾波器d抽樣判決器定時(shí)脈沖e輸出

cos?ct

b

圖 2-4 2PSK的相干解調(diào)原理圖

edb10011atTstctt10011t圖 2-5 相干解調(diào)中各點(diǎn)波形圖

2.2 2DPSK原理 2.2.1 2DPSK基本原理

二進(jìn)制差分相移鍵控常簡稱為二相相對(duì)調(diào)相，記為2DPSK。它不是利用載波相位的絕對(duì)數(shù)值傳送數(shù)字信息，而是用前后碼元的相對(duì)載波相位值傳送數(shù)字信息。所謂相對(duì) 4

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

載波相位是只本碼元初相與前一碼元初相之差。

傳輸系統(tǒng)中要保證信息的有效傳輸就必須要有較高的傳輸速率和很低的誤碼率。在傳輸信號(hào)中，2PSK信號(hào)和2ASK及2FSK信號(hào)相比，具有較好的誤碼率性能，但是，在2PSK信號(hào)傳輸系統(tǒng)中存在相位不確定性，并將造成接收碼元“0”和“1”的顛倒，產(chǎn)生誤碼。為了保證2PSK的優(yōu)點(diǎn)，又不會(huì)產(chǎn)生誤碼，將2PSK體制改進(jìn)為二進(jìn)制差分相移鍵控（2DPSK），及相對(duì)相移鍵控。

2DPSK方式即是利用前后相鄰碼元的相對(duì)相位值去表示數(shù)字信息的一種方式。現(xiàn)假設(shè)用Φ表示本碼元初相與前一碼元初相之差，并規(guī)定：Φ＝0表示0碼，Φ＝π表示1碼。則數(shù)字信息序列與2DPSK信號(hào)的碼元相位關(guān)系可舉例表示如2PSK信號(hào)是用載波的不同相位直接去表示相應(yīng)的數(shù)字信號(hào)而得出的，在接收端只能采用相干解調(diào)，它的時(shí)域波形圖見圖2-6。

圖2-6 2DPSK信號(hào)波形圖(a)絕對(duì)碼(b)相對(duì)碼10參考100011011(c)2DPSKt

2.2.2 2DPSK調(diào)制原理

二進(jìn)制差分相移鍵控。2DPSK方式是用前后相鄰碼元的載波相對(duì)相位變化來表示數(shù)字信息。假設(shè)前后相鄰碼元的載波相位差為??，可定義一種數(shù)字信息與??之間的關(guān)系為：

?0（數(shù)字信息“0”）????

?（數(shù)字信息“1 ??為前一碼元的相位。

實(shí)現(xiàn)二進(jìn)制差分相移鍵控的最常用的方法是：先對(duì)二進(jìn)制數(shù)字基帶信號(hào)進(jìn)行差分編碼，然后對(duì)變換出的差分碼進(jìn)行絕對(duì)調(diào)相即可。2DPSK調(diào)制原理圖如圖2-7所示。

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

絕對(duì)碼Dn相對(duì)碼BnCnS2dpsk(t)+延時(shí)Ts波形變換×Coswc(t)

圖2-7 2DPSK調(diào)制原理框圖

2.2.3 2DPSK解調(diào)原理

2DPSK信號(hào)解調(diào)有相干解調(diào)方式和差分相干解調(diào)。用差分相干解調(diào)這種方法解調(diào)時(shí)不需要恢復(fù)本地載波，只要將DPSK信號(hào)精確地延遲一個(gè)碼元時(shí)間間隔，然后與DPSK信號(hào)相乘，相乘的結(jié)果就反映了前后碼元的相對(duì)相位關(guān)系，經(jīng)低通濾波后直接抽樣判決即可恢復(fù)出原始的數(shù)字信息，而不需要在進(jìn)行差分解碼。

相干解調(diào)碼變換法及相干解調(diào)法的解調(diào)原理是，先對(duì)2DPSK信號(hào)進(jìn)行相干解調(diào)，恢復(fù)出相對(duì)碼，再通過碼反變換器變換為絕對(duì)碼，從而恢復(fù)出發(fā)送的二進(jìn)制數(shù)字信息。

在解調(diào)過程中，若相干載波產(chǎn)生180?相位模糊，解調(diào)出的相對(duì)碼將產(chǎn)生倒置現(xiàn)象，但是經(jīng)過碼反變換器后，輸出的絕對(duì)碼不會(huì)發(fā)生任何倒置現(xiàn)象，從而解決了載波相位模糊的問題。本次設(shè)計(jì)采用相干解調(diào)。兩種解調(diào)方式的原理圖如圖2-8和圖2-9所示。

2DPSK相干解調(diào)各點(diǎn)波形圖如圖 2-10所示。

圖 2-8 2DPSK差分相干解調(diào)原理圖

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

圖 2-9 2DPSK相干解調(diào)原理圖

第3章 實(shí)驗(yàn)過程

3.1 2PSK仿真部分 3.1.1 2PSK仿真圖

用MATLAB搭建好的2PSK仿真圖如下：

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

圖3-1PSK仿真圖

3.1.2 2PSK模塊的參數(shù)設(shè)置： 1）相乘模塊

圖3-2 相乘器參數(shù)設(shè)置

2）低通濾波器模塊

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

圖3-3 濾波器其參數(shù)設(shè)置

3）抽樣判決模塊

圖3-4 pulse generator 參數(shù)設(shè)置

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

3.2 2DPSK仿真部分 3.2.1 2DPSK仿真圖

用MATLAB搭建好的2DPSK仿真圖如下：

圖3-5 2DPSK仿真圖

2.2.2 2DPSK模塊的參數(shù)設(shè)置： 1）載波模塊

圖3-6 載波參數(shù)設(shè)置

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

2）乘法器模塊

圖3-7 乘法器參數(shù)設(shè)置

3）基帶模塊

圖3-8 基帶信號(hào)參數(shù)設(shè)置

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

4）Unipolar to Bipolar Converte模塊

圖3-9 Unipolar to Bipolar Converter參數(shù)設(shè)置

5）碼變換模塊

圖3-10 Logical Operator參數(shù)設(shè)置

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

圖3-11 Unit Delay參數(shù)設(shè)置

圖3-12 Data Type Conversion參數(shù)設(shè)置

6）濾波器模塊

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

圖3-13 帶通濾波器參數(shù)設(shè)置

圖3-14 低通濾波器參數(shù)設(shè)置

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

第4章 仿真結(jié)果

4.1 2PSK仿真結(jié)果

圖4-1 2PSK電路仿真波形

4.2 2DPSK仿真結(jié)果

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

圖4-2 2DPSK電路仿真波形

附錄：

通過編寫M文件程序： 2PSK調(diào)制解調(diào)程序及注釋 clear all close all i=10;j=5000;fc=4;%載波頻率 fm=i/5;%碼元速率 B=2*fm;t=linspace(0,5,j);a=round(rand(1,i));%隨機(jī)序列,基帶信號(hào) figure(3);stem(a);st1=t;16

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

for n=1:10 if a(n)<1;for m=j/i*(n-1)+1:j/i*n st1(m)=0;end else for m=j/i*(n-1)+1:j/i*n st1(m)=1;end end end figure(1);subplot(411);plot(t,st1);title('基帶信號(hào)st1');axis([0,5,-1,2]);%由于PSK中的是雙極性信號(hào)，因此對(duì)上面所求單極性信號(hào)取反來與之一起構(gòu)成雙極性碼 st2=t;

for k=1:j;

if st1(k)>=1;

st2(k)=0;

else

st2(k)=1;

end end;subplot(412);plot(t,st2);title('基帶信號(hào)反碼st2');axis([0,5,-1,2]);st3=st1-st2;subplot(413);plot(t,st3);title('雙極性基帶信號(hào)st3');axis([0,5,-2,2]);s1=sin(2*pi*fc*t);subplot(414);plot(s1);title('載波信號(hào)s1');e_psk=st3.*s1;figure(2);subplot(511);plot(t,e_psk);title('e_2psk');noise=rand(1,j);psk=e_psk+noise;

%加入噪聲 subplot(512);plot(t,psk);title('加噪后波形');psk=psk.*s1;

%與載波相乘 subplot(513);plot(t,psk);title('與載波s1相乘后波形');[f,af] = T2F(t,psk);

%通過低通濾波器

[t,psk] = lpf(f,af,B);subplot(514);plot(t,psk);title('低通濾波后波形');for m=0:i-1;

if psk(1,m*500+250)<0;

for j=m*500+1:(m+1)*500;

psk(1,j)=0;

end

else

for j=m*500+1:(m+1)*500;

psk(1,j)=1;

end

end end subplot(515);plot(t,psk);axis([0,5,-1,2]);title('抽樣判決后波形')

2DPSK調(diào)制解調(diào)程序及注釋 clear all close all i=10;j=5000;fc=4;%載波頻率 fm=i/5;%碼元速率B=2*fm;t=linspace(0,5,j);a=round(rand(1,i));figure(4);stem(a);st1=t;for n=1:10 if a(n)<1;for m=j/i*(n-1)+1:j/i*n st1(m)=0;end else for m=j/i*(n-1)+1:j/i*n st1(m)=1;end end end 江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

figure(1);subplot(321);plot(t,st1);title('絕對(duì)碼');axis([0,5,-1,2]);b=zeros(1,i);%全零矩陣 b(1)=a(1);for n=2:10 if a(n)>=1;if b(n-1)>=1 b(n)=0;else b(n)=1;end else b(n)=b(n-1);end end st1=t;for n=1:10 if b(n)<1;for m=j/i*(n-1)+1:j/i*n st1(m)=0;end else for m=j/i*(n-1)+1:j/i*n st1(m)=1;end end end subplot(323);plot(t,st1);title('相對(duì)碼st1');axis([0,5,-1,2]);st2=t;for k=1:j;if st1(k)>=1;st2(k)=0;else st2(k)=1;江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告 end end;subplot(324);plot(t,st2);title('相對(duì)碼反碼st2');axis([0,5,-1,2]);s1=sin(2*pi*fc*t);subplot(325);plot(s1);title('載波信號(hào)s1');s2=sin(2*pi*fc*t+pi);subplot(326);plot(s2);title('低通濾波后波形');st=zeros(1,i);

%全零矩陣for m=0:i-1;

if dpsk(1,m*500+250)<0;

st(m+1)=0;

for j=m*500+1:(m+1)*500;

dpsk(1,j)=0;

end

else

for j=m*500+1:(m+1)*500;

st(m+1)=1;

dpsk(1,j)=1;

end

end end subplot(413);plot(t,dpsk);axis([0,5,-1,2]);title('抽樣判決后波形')dt=zeros(1,i);

%全零矩陣 dt(1)=st(1);for n=2:10;

if(st(n)-st(n-1))<=0&&(st(n)-st(n-1))>-1;

dt(n)=0;

else

dt(n)=1;

end end st=t;for n=1:10

if dt(n)<1;

for m=j/i*(n-1)+1:j/i*n

st(m)=0;

end

else

for m=j/i*(n-1)+1:j/i*n

st(m)=1;

end

end end 江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

subplot(414);plot(t,st);axis([0,5,-1,2]);title('碼反變換后波形')21

第五篇：通信原理課程設(shè)計(jì)_(基于MATLAB的_2PSK_2DPSK仿真)

江西農(nóng)業(yè)大學(xué)

通信原理課程設(shè)計(jì)報(bào)告

題 目 基于Matlab的相移鍵控仿真設(shè)計(jì)

專 業(yè) 電子信息工程

學(xué)生姓名 曾凡文

學(xué) 號(hào) 20121206

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告 二 0 一五 年 六 月

基于Matlab的2PSK,2DPSK仿真

摘要：現(xiàn)代通信系統(tǒng)要求通信距離遠(yuǎn)、通信容量大、傳輸質(zhì)量好，作為其關(guān)鍵技術(shù)之一的調(diào)制技術(shù)一直是研究的一個(gè)重要方向。本設(shè)計(jì)主要敘述了數(shù)字信號(hào)的調(diào)制方式，介紹了2PSK數(shù)字調(diào)制方式的基本原理，功率譜密度，并運(yùn)用MATLAB軟件對(duì)數(shù)字調(diào)制方式2PSK進(jìn)行了編程仿真實(shí)現(xiàn),在MATLAB平臺(tái)上建立2PSK和2DPSK調(diào)制技術(shù)的仿真模型。進(jìn)一步學(xué)習(xí)了MATLAB編程軟件，將MATLAB與通信系統(tǒng)中數(shù)字調(diào)制知識(shí)聯(lián)系起來，為以后在通信領(lǐng)域?qū)W習(xí)和研究打下了基礎(chǔ)在計(jì)算機(jī)上，運(yùn)用MATLAB軟件來實(shí)現(xiàn)對(duì)數(shù)字信號(hào)調(diào)制技術(shù)的仿真。

課程設(shè)計(jì)目的：通過課程設(shè)計(jì)，鞏固已學(xué)過的*****知識(shí)，加深對(duì)其理解和應(yīng)用，學(xué)會(huì)應(yīng)用Matlab Simulink工具對(duì)通信系統(tǒng)仿真。

關(guān)鍵詞：數(shù)字調(diào)制與解調(diào)；MATLAB；2PSK；2DPSK；

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

第1章 基本工作原理

1.1 2PSK原理

1.1.1 2PSK基本原理

二進(jìn)制移相鍵控，簡記為2PSK或BPSK。2PSK信號(hào)碼元的“0”和“1”分別用兩個(gè)不同的初始相位“0”和“?”來表示，而其振幅和頻率保持不變.因此，2PSK信號(hào)的時(shí)域表達(dá)式為:

(t)=Acos其中，表示第n個(gè)符號(hào)的絕對(duì)相位：

t+)

=因此，上式可以改寫為:

這種以載波的不同相位直接表示相應(yīng)二進(jìn)制數(shù)字信號(hào)的調(diào)制方式,稱為二進(jìn)制移相鍵控方式。二進(jìn)制移相鍵控信號(hào)的典型時(shí)間波形如圖1-1。

10011tTs江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

圖1-1 二進(jìn)制相移鍵控信號(hào)的時(shí)間波形

1.1.2 2PSK調(diào)制原理

在二進(jìn)制數(shù)字調(diào)制中，當(dāng)正弦載波的相位隨二進(jìn)制數(shù)字基帶信號(hào)離散變化時(shí)，則產(chǎn)生二進(jìn)制移相鍵控(2PSK)信號(hào)。2PSK信號(hào)調(diào)制有兩種方法，即模擬調(diào)制法和鍵控法。通常用已調(diào)信180°分別表示信號(hào)的 1 和 兩個(gè)反相的載制。2PSK以載

號(hào)載波的 0°和 二進(jìn)制數(shù)字基帶0，模擬調(diào)制法用波信號(hào)進(jìn)行調(diào)波的相位變化作為參考基準(zhǔn)的，當(dāng)基帶信號(hào)為0時(shí)相位相對(duì)于初始相位為0°，當(dāng)基帶信號(hào)為1時(shí)相對(duì)于初始相位為180°。鍵控法，是用載二進(jìn)制信息的調(diào)制和180°來分別代表達(dá)式為：

波的相位來攜帶方式。通常用0°表0和1。其時(shí)域

??e2PSK???ang(t?nTs)?cos?ct

?n?其中，2PSK的調(diào)制中an必須為雙極性碼。兩種方法原理圖分別如圖1-2和圖1-3所示。

圖1-2 模擬調(diào)制原理圖

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

圖 1-3 鍵控法原理圖

帶通濾波器ae2PSK(t)相乘器c低通濾波器d抽樣判決器定時(shí)脈沖e輸出

cos?ct

1.1.3 2PSK解調(diào)原理

b由于2PSK的幅度是恒定的，必須進(jìn)行相干解調(diào)。經(jīng)過帶通濾波的信號(hào)在相乘器中與本地載波相乘，然后用低通濾波器濾除高頻分量，在進(jìn)行抽樣判決。判決器是按極性來判決的。即正抽樣值判為1，負(fù)抽樣值判為0。2PSK信號(hào)的相干解調(diào)原理圖如圖1-4所示，各點(diǎn)的波形如圖1-5所示。

由于2PSK信號(hào)的載波回復(fù)過程中存在著180°的相位模糊，即恢復(fù)的本地載波與所需相干載波可能相同，也可能相反，這種相位關(guān)系的不確定性將會(huì)造成解調(diào)出的數(shù)字基帶信號(hào)與發(fā)送的基帶信號(hào)正好相反，即“1”變成“0”嗎“0”變成“1”，判決器輸出數(shù)字信號(hào)全部出錯(cuò)。這種現(xiàn)象稱為2PSK方式的“倒π”現(xiàn)象或“反相工作”。

圖 1-4 2PSK的相干解調(diào)原理圖

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

edb10011atTstctt10011t

圖 1-5 相干解調(diào)中各點(diǎn)波形圖

1.2 2DPSK原理

1.2.1 2DPSK基本原理

二進(jìn)制差分相移鍵控常簡稱為二相相對(duì)調(diào)相，記為2DPSK。它不是利用載波相位的絕對(duì)數(shù)值傳送數(shù)字信息，而是用前后碼元的相對(duì)載波相位值傳送數(shù)字信息。所謂相對(duì)載波相位是只本碼元初相與前一碼元初相之差。

傳輸系統(tǒng)中要保證信息的有效傳輸就必須要有較高的傳輸速率和很低的誤碼率。在傳輸信號(hào)中，2PSK信號(hào)和2ASK及2FSK信號(hào)相比，具有較好的誤碼率性能，但是，在2PSK信號(hào)傳輸系統(tǒng)中存在相位不確定性，并將造成接收碼元“0”和“1”的顛倒，產(chǎn)生誤碼。為了保證2PSK的優(yōu)點(diǎn)，又不會(huì)產(chǎn)生誤碼，將2PSK體制改進(jìn)為二進(jìn)制差分相移鍵控（2DPSK），及相對(duì)相移鍵控。

2DPSK方式即是利用前后相鄰碼元的相對(duì)相位值去表示數(shù)字信息的一種方式?，F(xiàn)假設(shè)用Φ表示本碼元初相與前一碼元初相之差，并規(guī)定：Φ＝0表示0碼，Φ＝π表示1碼。則數(shù)字信息序列與2DPSK信號(hào)的碼元相位關(guān)系可舉例表示如2PSK信號(hào)是用載波的不同相位直接去表示相應(yīng)的數(shù)字信號(hào)而得出的，在接收端只能采用相干解調(diào)，它的時(shí)域波形圖見圖1-6。

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

圖1-6 2DPSK信號(hào)波形圖

(a)絕對(duì)碼(b)相對(duì)碼10參考100011011(c)2DPSKt1.2.2 2DPSK調(diào)制原理

二進(jìn)制差分相移鍵控。2DPSK方式是用前后相鄰碼元的載波相對(duì)相位變化來表示數(shù)字信息。假設(shè)前后相鄰碼元的載波相位差為??，可定義一種數(shù)字信息與??之間的關(guān)系為：

?0（數(shù)字信息“0”）????

??（數(shù)字信息“1 ?為前一碼元的相位。

實(shí)現(xiàn)二進(jìn)制差分相移鍵控的最常用的方法是：先對(duì)二進(jìn)制數(shù)字基帶信號(hào)進(jìn)行差分編碼，然后對(duì)變換出的差分碼進(jìn)行絕對(duì)調(diào)相即可。2DPSK調(diào)制原理圖如圖1-7所示。

絕對(duì)碼Dn相對(duì)碼BnCnS2dpsk(t)+延時(shí)Ts波形變換×Coswc(t)

圖1-7 2DPSK調(diào)制原理框圖

1.2.3 2DPSK解調(diào)原理

2DPSK信號(hào)解調(diào)有相干解調(diào)方式和差分相干解調(diào)。用差分相干解調(diào)這種方法解調(diào)時(shí)不需要恢復(fù)本地載波，只要將DPSK信號(hào)精確地延遲一個(gè)碼元時(shí)間間隔，然后與DPSK信號(hào)相乘，相乘的結(jié)果就反映了前后碼元的相對(duì)相位關(guān)系，經(jīng)低通濾波后直接抽樣判決即可恢復(fù)出原始的數(shù)字信息，而不需要在進(jìn)行差分解碼。

第二章 設(shè)計(jì)系統(tǒng)

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告 2.1框圖

兩種解調(diào)方式的原理框圖如圖1-8和圖1-9所示。

圖 1-8 2DPSK差分相干解調(diào)原理框圖

圖 1-9 2DPSK相干解調(diào)原理框圖

2.2工作原理

相干解調(diào)碼變換法及相干解調(diào)法的解調(diào)原理是，先對(duì)2DPSK信號(hào)進(jìn)行相干解調(diào)，恢復(fù)出相對(duì)碼，再通過碼反變換器變換為絕對(duì)碼，從而恢復(fù)出發(fā)送的二進(jìn)制數(shù)字信息。

在解調(diào)過程中，若相干載波產(chǎn)生180?相位模糊，解調(diào)出的相對(duì)碼將產(chǎn)生倒置現(xiàn)象，但是經(jīng)過碼反變換器后，輸出的絕對(duì)碼不會(huì)發(fā)生任何倒置現(xiàn)象，從而解決了載波相位模糊的問題。本次設(shè)計(jì)采用相干解調(diào)。

2.3設(shè)定參數(shù) 如附錄1

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

第三章 Matlab仿真

3.1顯示系統(tǒng)不同部分的信號(hào)波形

3.1.1PSK如圖3-1和圖3-2所示

圖3-1

圖3-2

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

3.1.2 2DPSK如圖3-

3、圖3-4和圖3-5所示

圖3-3

圖3-4

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告

圖3-5 3.2各種相移系統(tǒng)的比較和分析

例如“倒π”現(xiàn)象

對(duì)于相同的數(shù)字信號(hào)基帶序列，由于初始相位不同，2DPSK信號(hào)的相位并不直接代表基帶信號(hào)，而前后碼元相對(duì)相位的差才唯一決定信號(hào)的符號(hào)。2PSK信號(hào)載波恢復(fù)過程中，存在著180°的相位模糊即恢復(fù)的本地載波與與所需的想干載波可能同相也可能反相，這種相位關(guān)系的不確定性將會(huì)造成解調(diào)出來的數(shù)字基帶信號(hào)與發(fā)送的數(shù)字基帶信號(hào)正好相反，即“1”變成“0”，“0”變成“1”，判決器輸出的數(shù)字信號(hào)全部出錯(cuò), 這種現(xiàn)象稱為2PSK的“倒π”現(xiàn)象或“反相工作”。本地載波與發(fā)送端載波反向時(shí)，2PSK的解調(diào)波形與2DPSK完全相反

3.3不同方式解調(diào)下PSK、DPSK的誤碼率

誤碼率是指接收的碼元數(shù)在傳輸總碼元數(shù)中所占的比例，即：

誤碼率錯(cuò)誤碼元數(shù)傳輸總碼元數(shù) PSK相干解調(diào)：

誤碼率（r/2）DPSK相干解調(diào)：

誤碼率 DPSK差分非相干解調(diào)：

誤碼率1/2e^(-r)r為信噪比。

誤碼率是衡量一個(gè)數(shù)字通信系統(tǒng)性能的重要指標(biāo)。在信道高斯白噪聲的干擾下，各種二進(jìn)制數(shù)字調(diào)制系統(tǒng)的誤碼率取決于解調(diào)器輸入信噪比，而誤碼率表達(dá)式的形式則取決于解調(diào)方式。對(duì)于所有的數(shù)字調(diào)制系統(tǒng)誤碼率與信噪比的關(guān)系的圖表來看，所有的曲線呈減函數(shù)的下降曲線，即隨著信噪比的增大，誤碼率降低。橫向比較來看，對(duì)于同一種調(diào)制方式，當(dāng)信噪比相同時(shí)，采用相干解調(diào)方式的誤碼率低于非相干解調(diào)方式的誤碼率；縱向比較來看，對(duì)2PSK,2DPSK兩種調(diào)制方式若采用同一種解調(diào)方式

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告（相干解調(diào)或非相干解調(diào)），則2PSK的誤碼率最低，2DSPK的誤碼率次之。當(dāng)信噪比一定時(shí)，誤碼率由低到高依次是：2PSK的相干解調(diào)，2DPSK的相干解調(diào)，2DPSK的差分解的非相干調(diào)。

附錄1 2PSK調(diào)制解調(diào)程序及注釋： clear all close all i=10;j=5000;fc=4.6;%載波頻率 fm=i/5;%碼元速率 B=2*fm;t=linspace(0,5,j);a=round(rand(1,i));%隨機(jī)序列,基帶信號(hào) figure(3);stem(a);st1=t;for n=1:10 if a(n)<1;for m=j/i*(n-1)+1:j/i*n st1(m)=0;end else for m=j/i*(n-1)+1:j/i*n st1(m)=1;end end end figure(1);subplot(411);plot(t,st1);title('基帶信號(hào)st1');axis([0,5,-1,2]);%由于PSK中的是雙極性信號(hào)，因此對(duì)上面所求單極性信號(hào)取反來與之一起構(gòu)成雙極性碼

st2=t;for k=1:j;if st1(k)>=1;st2(k)=0;else

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告 st2(k)=1;end end;subplot(412);plot(t,st2);title('基帶信號(hào)反碼st2');axis([0,5,-1,2]);st3=st1-st2;subplot(413);plot(t,st3);title('雙極性基帶信號(hào)st3');axis([0,5,-2,2]);s1=sin(2*pi*fc*t);subplot(414);plot(s1);title('載波信號(hào)s1');e_psk=st3.*s1;figure(2);subplot(511);plot(t,e_psk);title('e_2psk');noise=rand(1,j);psk=e_psk+noise;%加入噪聲 subplot(512);plot(t,psk);title('加噪后波形');psk=psk.*s1;%與載波相乘 subplot(513);plot(t,psk);title('與載波s1相乘后波形');[f,af] = T2F(t,psk);%通過低通濾波器 [t,psk] = lpf(f,af,B);subplot(514);plot(t,psk);title('低通濾波后波形');for m=0:i-1;if psk(1,m*500+250)<0;for j=m*500+1:(m+1)*500;psk(1,j)=0;end else for j=m*500+1:(m+1)*500;psk(1,j)=1;end

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告 end end subplot(515);plot(t,psk);axis([0,5,-1,2]);title('抽樣判決后波形')

2DPSK調(diào)制解調(diào)程序及注釋: clear all close all i=10;j=5000;fc=4.6;%載波頻率 fm=i/5;%碼元速率 B=2*fm;t=linspace(0,5,j);a=round(rand(1,i));figure(4);stem(a);st1=t;for n=1:10 if a(n)<1;for m=j/i*(n-1)+1:j/i*n st1(m)=0;end else for m=j/i*(n-1)+1:j/i*n st1(m)=1;end end end figure(1);subplot(321);plot(t,st1);title('絕對(duì)碼');axis([0,5,-1,2]);b=zeros(1,i);%全零矩陣 b(1)=a(1);for n=2:10 if a(n)>=1;

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告 if b(n-1)>=1 b(n)=0;else b(n)=1;end else b(n)=b(n-1);end end st1=t;for n=1:10 if b(n)<1;for m=j/i*(n-1)+1:j/i*n st1(m)=0;end else for m=j/i*(n-1)+1:j/i*n st1(m)=1;end end end subplot(323);plot(t,st1);title('相對(duì)碼st1');axis([0,5,-1,2]);st2=t;for k=1:j;if st1(k)>=1;st2(k)=0;else st2(k)=1;end end;subplot(324);plot(t,st2);title('相對(duì)碼反碼st2');axis([0,5,-1,2]);s1=sin(2*pi*fc*t);subplot(325);plot(s1);title('載波信號(hào)s1');s2=sin(2*pi*fc*t+pi);subplot(326);plot(s2);

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告 title('載波信號(hào)s2');d1=st1.*s1;d2=st2.*s2;figure(2);subplot(411);plot(t,d1);title('st1*s1');subplot(412);plot(t,d2);title('st2*s2');e_dpsk=d1+d2;subplot(413);plot(t,e_dpsk);title('調(diào)制后波形');noise=rand(1,j);dpsk=e_dpsk+noise;%加入噪聲 subplot(414);plot(t,dpsk);title('加噪聲后信號(hào)');dpsk=dpsk.*s1;%與載波s1相乘 figure(3);subplot(411);plot(t,dpsk);title('與載波相乘后波形');[f,af]=T2F(t,dpsk);%通過低通濾波器 [t,dpsk]=lpf(f,af,B);subplot(412);plot(t,dpsk);title('低通濾波后波形');st=zeros(1,i);%全零矩陣 for m=0:i-1;if dpsk(1,m*500+250)<0;st(m+1)=0;for j=m*500+1:(m+1)*500;dpsk(1,j)=0;end else for j=m*500+1:(m+1)*500;st(m+1)=1;dpsk(1,j)=1;end end end subplot(413);

江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告 plot(t,dpsk);axis([0,5,-1,2]);title('抽樣判決后波形')dt=zeros(1,i);%全零矩陣 dt(1)=st(1);for n=2:10;if(st(n)-st(n-1))<=0&&(st(n)-st(n-1))>-1;dt(n)=0;else dt(n)=1;end end st=t;for n=1:10 if dt(n)<1;for m=j/i*(n-1)+1:j/i*n st(m)=0;end else for m=j/i*(n-1)+1:j/i*n st(m)=1;end end end subplot(414);plot(t,st);axis([0,5,-1,2]);title('碼反變換后波形');