數(shù)字通信系統(tǒng)課程設(shè)計(jì)

1、成都學(xué)院（成都大學(xué)）課程設(shè)計(jì)報告數(shù)字通信系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)摘要：本設(shè)計(jì)為掌握利用matlab加深對2dps瞰字頻帶通信系統(tǒng)的理解與掌握，理解運(yùn)用所學(xué)的知識。我主要是用模塊和程序設(shè)計(jì)出數(shù)字通信系統(tǒng)中的信源、信源編碼、調(diào)制器、噪聲、信道、解調(diào)器、信源譯碼、信宿。通過 matlab真平臺，運(yùn)用所學(xué)的理論和方法進(jìn)行仿真、調(diào)試、波形眼圖分析，最終成功實(shí)現(xiàn)了 2dps敝字通信系統(tǒng)。關(guān)鍵詞：matlab 2dpsk調(diào)制；解調(diào)3目錄第1章弓i言 11.1 背景 11.2 選題的目的和意義 11.3 本課程設(shè)計(jì)的主要內(nèi)容 2第2章2dpsk基本原理 32.1 2dpsk言號的原理 32.2 2dpsk信號的調(diào)制原

2、理 32.2.1 2dpsk 調(diào)制 32.2.2 模擬調(diào)制法 32.2.3 鍵控法 42.3 2dpsk信號的解調(diào)原理 42.3.1 2dpsk 解調(diào) 42.3.2 極性比較法 62.3.3 相位比較法 6第3章2dpsk系統(tǒng)模塊設(shè)計(jì)仿真 73.1 模擬調(diào)制法和極性比較法構(gòu)成的 2dpskk統(tǒng) 73.2 模擬調(diào)制法和相位比較法構(gòu)成的 2dpskk統(tǒng) 73.3 鍵控法和相位比較法構(gòu)成的 2dpsk!（統(tǒng) 83.4 模擬調(diào)制法和極性比較法模塊分析 83.4.1 模擬調(diào)制法模塊 83.4.2 鍵控法調(diào)制模塊 93.4.3 模擬信道模塊 103.4.4 極性比較法模塊 103.4.5 相位比較法模塊

3、113.4.6 誤碼率模塊 123.4.7 延時器和觀測模塊 123.4.8 眼圖模塊 133.5 模塊調(diào)試 143.5.1 模擬調(diào)制法極性調(diào)解法仿真 143.5.2 模擬調(diào)制法極性調(diào)解法仿真分析 16第4章2dpsk系統(tǒng)程序設(shè)計(jì)仿真 184.1 matla程序 184.2 仿真波形 204.3 波形分析 21第5章結(jié)束語 225.1 綜合總結(jié) 225.2 心得 22參考文獻(xiàn) 23成都學(xué)院（成都大學(xué)）課程設(shè)計(jì)報告第1章引言a a1.1 目樂數(shù)字通信是通信技術(shù)和計(jì)算機(jī)技術(shù)相結(jié)合而產(chǎn)生的一種新的通信方式。要在兩地間傳輸信息必須有傳輸信道，根據(jù)傳輸媒體的不同有有線數(shù)字通信與無線數(shù)據(jù)通信之分。但它們都

4、是通過傳輸信道將數(shù)據(jù)終端與計(jì)算機(jī)聯(lián)結(jié)起來，而使不同地點(diǎn)的數(shù)字終端實(shí)現(xiàn)軟、硬件和信息資源的共享。為使數(shù)字信號能在帶通信道中傳輸，必須用數(shù)字信號對載波進(jìn)行調(diào)制，其調(diào)制方式與模擬信號調(diào)制類似。根據(jù)數(shù)字信號控制波的參量不同也分為調(diào)幅、調(diào)頻、調(diào)相三種方式。因數(shù)字信號對載波參數(shù)的調(diào)制通常采用數(shù)字信號的離散值對載波進(jìn)行鍵控，故這三種數(shù)字調(diào)制方式被稱為幅移鍵控頻移鍵控和相移鍵控。matlab是由mathworks公司于1984年推出的一種面向科學(xué)與工程的設(shè)計(jì)的計(jì)算機(jī)軟件，它將 不同的領(lǐng)域的計(jì)算用函數(shù)的形式提供給給用戶；用戶在使用時，只需要用這些函數(shù)并賦予實(shí)際參數(shù)就能解決實(shí)際問題，它涉及數(shù)值分析、自動控制、數(shù)字

5、信號處理、圖像處理、小波分析及神經(jīng) 元網(wǎng)絡(luò)等十幾個領(lǐng)域的計(jì)算和圖形顯示，而且隨著新出版的推出，涉及的領(lǐng)域更多，功能強(qiáng)大。 matla棍供實(shí)現(xiàn)動態(tài)系統(tǒng)建模和仿真的軟件包，它讓用戶把精力從編程轉(zhuǎn)向模型的構(gòu)造，為用戶省 去了許多重復(fù)的代碼編寫工作；用戶只須知道模塊的輸入、輸出以及模塊的功能，而不必管模塊內(nèi) 部是怎么實(shí)現(xiàn)的，于是留給用戶的事情就是如何利用這些模塊來建立模型以完成自己的仿真任務(wù)， 它被廣泛的應(yīng)用在信號仿真中。本課程設(shè)計(jì)報告主要介紹了用用窗函數(shù)法設(shè)計(jì)一個線性相位fir低通濾波器和用雙線性變化法設(shè)計(jì)一個butterworth低通濾波器，同時還用了 matlab軟件進(jìn)行仿真設(shè)計(jì)。1.2 選題的

6、目的和意義在傳輸信號中，2dpskw 2psk信號和2ask及2fsk信號相比，具有較好的誤碼率性能， 但2fsk對相位不敏感，為了保證 2psk的優(yōu)點(diǎn)，又不會產(chǎn)生誤碼，將2psk體制改進(jìn)為二進(jìn)制差分相移鍵控（2dpsk,及相對相移鍵控。2dpsk方式即是利用前后相鄰碼元的相對相位值去表示數(shù)字信息的一種方式。所以我們進(jìn)行數(shù)字通信系統(tǒng)設(shè)計(jì)可以讓我們在設(shè)計(jì)中獲得科學(xué)信息，培養(yǎng)自己的邏輯能力及個人設(shè)計(jì)能力，同時我還可以熟悉使用matlaba件，練習(xí) wordc檔的用法，及函數(shù)的編譯，流程圖的設(shè)計(jì)等等，所以這次課程設(shè)計(jì)很有意義。matla的真平臺，運(yùn)用所本設(shè)計(jì)為掌握利用計(jì)算機(jī)來加深對所學(xué)知識的理解和掌

7、握，通過 學(xué)的理論和方法進(jìn)行仿真、解決問題。1.3 本課程設(shè)計(jì)的主要內(nèi)容設(shè)計(jì)2dps/字通信系統(tǒng)分別設(shè)計(jì)各模塊和編寫代碼程序，進(jìn)行仿真、調(diào)試、分析波形、評價系統(tǒng)、心得體會。通信系統(tǒng)基本模型如圖1-1圖1-1通信系統(tǒng)基本模型3成都學(xué)院(成都大學(xué))課程設(shè)計(jì)報告第2章2dpsks本原理2.1 2dpsk信號的原理2dpsk是利用前后相鄰碼元的載波相對相位變化傳遞數(shù)字信息，載波相鄰兩碼元的相位差定義 為q-1(2-1)%、q-1分別表示第n及n-1個碼元的載波初相。 鴛=0通常表示數(shù)字信息“ 0"，中n=n 通常表示數(shù)字信息“ 1 ”。2.2 2dpsk信號的調(diào)制原理2dpsk信號有兩種調(diào)制

8、方法，分別是模擬調(diào)制法和鍵控法。2.2.1 2dpsk 調(diào)制2dpsk的基本原理和2ask是一樣的，只是把輸入的數(shù)字信息進(jìn)行碼變換，s2dpsk(t)等于調(diào)制信號s(t)碼變換后的信號s乘以載波信號cos2n fct。所以sdpsk的數(shù)學(xué)表達(dá)式為 's2dpsk (t) = s(t)，cos2h fct(2-2 )根據(jù)功率譜公式可以算出s(t)'的功率譜，用p(f)來表示。根據(jù)頻移定理得到s(t)信號的功率譜a2 -1f2dpsjf)=ps(f fc)+ps(f + fc)(2-3)4 -其中，ps (f )為雙極性全占空矩形脈沖序列s(t)的功率譜。2.2.2模擬調(diào)制法模擬調(diào)

9、制法如圖2-1所示，其中碼變換過程為將絕對碼變換為相對嗎；碼型變換的過程為將輸 入的單極性不歸零碼轉(zhuǎn)換為雙極性不歸零碼型；乘法器過程是將雙極性不歸零信號與載波相乘得到圖2-1模擬調(diào)制法2.2.3鍵控法鍵控法如圖2-2所示，其中差分變換功能同圖1的碼變換；選項(xiàng)開關(guān)的作用是輸入“0”時接相位0,輸入“ 1”時接相位兀。圖2-2鍵控法2.3 2dpsk信號的解調(diào)原理2dpsk信號最常用的解調(diào)方法有兩種，一種是極性比較法(相干解調(diào)法)，另一種是相位比較法(非相干解調(diào)法)。2.3.1 2dpsk 解調(diào)當(dāng)采用2psk解調(diào)時，設(shè)調(diào)制采用“ 1”變“0”不變規(guī)則。當(dāng)發(fā)送端“ 1”時，收到的2psk信 號為s2

10、psk = c0s2n fct(2-4)帶通濾波器的輸出是信號加窄帶噪聲cos2n fct 十 ni (t) = 一1 + ni (t)】cos2n fct - nq (t)sin2 n fct(2-5)上式與相干載波cos2n fct相乘，得i - cos2 .fctni(t)lcos2 二fct-1-1ni (t)bos2 2二fct-nq(t)sin2 二 fctcos2二fct5成都學(xué)院（成都大學(xué)）課程設(shè)計(jì)報告1 .1 1 , 一=一-1 ni(t)l- 1-1 ni(t) icos4r： fct 22-rq (t)sin 4n fct cos2n fct(2-6)式(2-6)所示信號

11、經(jīng)低通濾波器后得x(t) = 1 +ni (t)(2-7)顯然，x(t)的瞬時值是均值為-1、方差為6： = r0b2psk =2n0f的高斯隨機(jī)變量。所以，x(t)的 取樣值的概率密度函數(shù)為(x -1)21一 tti(x) = e 注(2-8),2二y同理，發(fā)端發(fā)“ 0”時，收到的2psk信號為s2psk = cos2n fct(2-9)帶通濾波器的輸出是信號加窄帶噪聲cos2冗 fct + 丹(t) = 1 + ni (t) cos2n fct 一 nq (t)sin2 冗 fct(2-10)上式與相干載波cos2n fct相乘，得cos2二 fctni (t) lcos2二fct- 1n

12、i (t) lcos2 2二fct-nq(t)sin2二fctcos2二fct=11 ni(t) ,1 -1 1 ni (t) lcos4二 fct22nq(t)sin 4n fct cos2n fct(2-11)式(2-11)所示信號經(jīng)低通濾波器后得x(t) =1 +ni (t)(2-12)顯然，x(t)的瞬時值是均值為1、方差為6： =n0b2psk =2n0f的高斯隨機(jī)變量。所以，x(t)的 取樣值的概率密度函數(shù)為(x4)2127f0(x) = j r e(2-13)二2:、當(dāng)“ 1”、“o”等概率時，最彳i判決門限為0。發(fā)“ 1”錯判成“ 0”的概率為p(0/1) =f1 (x)dx

13、=1 erfc( jr)(2-14)02根據(jù)r = p(0)p(1/0) + p(1)p(0/1)得解調(diào)器平均誤碼率為1 ，1 廠r =-erfc(v7)【p(0) + p(1)】 = erfc(v7)(2-15)2 2式中，=注。 、：2psk的反向工作問題：二分頻電路恢復(fù)的載波有時與發(fā)光載波相同，有時反相。當(dāng)本地載波反相，變?yōu)閏0s(2nfct+ n )時，則相乘器以后的輸出波形都和載波同頻同相時的情況相反，判決器輸出的數(shù)字信號全錯，與發(fā)送數(shù)碼完全相反，這種情況稱為反向工作。2dpsk只是解決了 2psk的反向工作問題，但是在 2dpsk差分碼中有一個出錯時會引起兩個相鄰碼元錯誤。所以 2

14、dpsk誤碼率為2psk誤碼率的兩倍1pe =2父erfc(tt)p(0) +p(1)=erfc"r)(2-16)22.3.2 極性比較法解調(diào)原理：對2dpsk言號進(jìn)行相干解調(diào)，恢復(fù)出想相對碼，再經(jīng)碼反變換器變換為絕對碼，從 而恢復(fù)出發(fā)送的二進(jìn)制數(shù)字信息。圖2-3極性比較法解調(diào)2.3.3 相位比較法解調(diào)原理：對接收到的 2dpsk言號延時一個碼元間隔 ts,然后與2dpsk言號本身相乘，相乘結(jié) 果反映了前后碼元的相位差，經(jīng)過低通濾波器后再抽樣判決，可直接恢復(fù)出原始數(shù)字信息。圖2-4相位比較法解調(diào)9第3章2dpsk系統(tǒng)模塊設(shè)計(jì)仿真3.1模擬調(diào)制法和極性比較法構(gòu)成的2dpsk（統(tǒng)如圖3-

15、1所示。他出幅通buimrbamaulli日置mull曜 rmatmprmudfiimhjaiar3ag fttaeu0rlium maysreajnhg fter defgrtri®eied subreemtiucwdria typa coniskfl（冷卻 犀明 bdgutrirpubegfliwicrun pd 5r i.：cofiwtflrztfn.onhfhoidem瞄看擊酒店痂電由金語圖3-1模擬調(diào)制和極性比較法構(gòu)成的2dpsk系統(tǒng)3.2 模擬調(diào)制法和相位比較法構(gòu)成的 2dpsk（統(tǒng)如圖3-2所示。j-unpotar tubipdsr dmh!erterpipdiktll

16、um glx ta bh polar conweterbuttershe wsvbrelay1小gau-saafianatag rber de/rhgaussian noise genmturtransport ddayfitar de 現(xiàn) ntrinfld subsyuqebuu4rn.iztroorderhddll 但 zefo-onfer hio w圖3-2模擬調(diào)制和相位比較法構(gòu)成的2dpsk系統(tǒng)成都學(xué)院(成都大學(xué))課程設(shè)計(jì)報告3.3 鍵控法和相位比較法構(gòu)成的2dpsk（統(tǒng)11圖3-3鍵控法和相位比較法構(gòu)成的2dpsk!（統(tǒng)3.4 模擬調(diào)制法和極性比較法模塊分析3.4.1 模擬調(diào)制法模塊

17、如圖3-4所示。sine wave圖3-4模擬調(diào)制法調(diào)制模塊模塊功能說明：通過伯努利二進(jìn)制發(fā)生器模塊(bernoulli binary generator)產(chǎn)生二進(jìn)制序列模擬數(shù)字基帶信號，然后使用differential encoder模塊對該基帶信號進(jìn)行差分編碼，unipolar tobipolar模塊將前面所得的單極性差分碼轉(zhuǎn)換成雙極性差分碼，之后使用相乘模塊( product )把雙 極性差分碼與sine wave模塊產(chǎn)生的載波信號相乘，輸出即是2dpsk調(diào)整信號。參數(shù)設(shè)置(沒提及的都默認(rèn))：伯努利二進(jìn)制發(fā)生器模塊 (bernoulli binary generator)-抽樣時間 sa

18、mple time:0.001differential encoder模塊-m-arry number:2sine wave 模塊-ferquency:3000*pi3.4.2 鍵控法調(diào)制模塊如圖3-5所示。sine wave圖3-5鍵控法調(diào)制模型模塊模塊功能：通過伯努利二進(jìn)制發(fā)生器模塊(bernoulli binary generator)產(chǎn)生二進(jìn)制序列模擬數(shù)字基帶信號，transport delay模塊對sine wave模塊產(chǎn)生的正弦波延時半個周期，相當(dāng)于將正弦波反相，switch模塊為一個邏輯開關(guān)，通過接受differential encoder輸出的單極性差分碼的控制，以決定接受正弦

19、波或反相正弦波，完成 2psk調(diào)制，綜合整個過程，最后輸出所需要的2dpsk調(diào)制信號。參數(shù)設(shè)置(沒提及的都默認(rèn))：伯努利二進(jìn)制發(fā)生器模塊(bernoulli binary generator)-抽樣時間 sample time:0.001differential encoder模塊-m-arry number:2sine wave 模塊-ferquency:3000*pitransport delay模塊-time delay:1/33.4.3模擬信道模塊如圖3-6所示。oaussign noise generator圖3-6模擬信道模型模塊功能：通過對 2dpsk調(diào)制信號輸出加入一個gaus

20、sian噪聲信號，可以來模擬信號在信道中的傳輸過程，通過加法器模塊加入gaussian噪聲信號3.4.4 極性比較法模塊圖3-7極性比較法解調(diào)模塊模塊功能：analog filter design模塊設(shè)計(jì)為一個帶通濾波器，接受信道輸出的信號，基本濾去與所需信號不相關(guān)的雜波信號，然后通過相乘器 product將帶通輸出的信號與調(diào)制時所用的正弦 波信號進(jìn)行相乘，是信號波形產(chǎn)生變化，并進(jìn)行頻譜搬移，所得的信號含高頻和低頻部分，再 經(jīng)過analog filter design模塊設(shè)計(jì)的低通濾波器，將高頻信號過濾，同時達(dá)到對波形整形的目的。pulsegenerator,triggered subsyst

21、em 和relay三個模塊構(gòu)成抽樣判決器，對低通輸出的信號進(jìn)行抽樣成都學(xué)院（成都大學(xué)）課程設(shè)計(jì)報告判決，還原出差分變換得到的相對碼，最后通過由 logical operator 和unit delay 模塊構(gòu)成的逆 差分變換模塊組，將判決出的相對碼變換出和基帶信號相同的絕對碼。data type conversion為一功能模塊，作用是對模塊輸出的數(shù)字類型進(jìn)行轉(zhuǎn)換，使變換成simulink仿真時能識別的double類型。參數(shù)設(shè)置（沒提及的都默認(rèn)）：anglog filter design帶通模塊-filter type:bandpass lower:8*pi upper:8000*piangl

22、og filter design低通模塊-filter type:lowpass passban:3000*pipulse generator模塊-period : 0.001 pulse:5logical operator模塊-simple time:0.001unit delay 模塊-simple time:0.001data type conversion模塊-output data type:double3.4.5 相位比較法模塊模塊功能：analog filter design模塊設(shè)計(jì)為一個帶通濾波器，接受信道輸出的信號，基本濾去與所需信號不相關(guān)的雜波信號，transport de

23、lay模塊的作用是對帶通輸出的信號進(jìn)行延時一個周期，得到差分波形，然后與延時前的波形相乘，形成差分相干過程，之后通過低通和進(jìn)行抽樣判決器,將高頻信號過濾，同時達(dá)到對波形整形的目的。pulse generator,triggered subsystem 和relay三個模塊構(gòu)成抽樣判決器。通過相位比較法抽樣判決出的碼型即為絕對碼，無需進(jìn)行碼反變換過程。如圖3-8所示。21transportrelay圖3-8相位比較法解調(diào)模塊參數(shù)設(shè)置（沒提及的都默認(rèn)）：anglog filter design帶通模塊-filter type:bandpass lower:8*pi upper:8000*piang

24、log filter design低通模塊-filter type:lowpass passban:3000*pipulse generator模塊-period : 0.001 pulse:5transport delay模塊-time delay:0.0013.4.6 誤碼率模塊如圖3-9所示。zero-order hold圖3-9誤碼率模塊模塊功能：zero-order-hold模塊對采樣時間進(jìn)行零階保持，對輸入進(jìn)行固定頻率取樣,輸出離散信號。error-rate calculation是誤碼率計(jì)算模塊，計(jì)算的誤碼率顯示在display模塊上。3.4.7延時器和觀測模塊如圖3-10所示。

25、unit oeilayw mlpower spectral density圖3-10延遲器和觀測模塊power spectral density 模塊模塊功能：unit delay 為延遲器模塊，scope2為示波器模塊, 用來觀測信號功率在頻譜上的分布。3.4.8眼圖模塊如圖3-11所示。discrete-time eye diagram scope圖3-11眼圖模塊模塊功能：輸出代觀測點(diǎn)信號的眼圖參數(shù)設(shè)置：由最佳觀測條件決定，在調(diào)試中設(shè)置。3.5模塊調(diào)試3.5.1 模擬調(diào)制法極性調(diào)解法仿真圖3-12調(diào)制過程波形圖3-14解調(diào)過程波形圖3-15輸入輸出波形比較-1mkf isqutncy l

26、iadml)富1q+time holorypumir sf»d/jl daftsiypmwi spectral 0e«lyle*i35?l醯6 o.n d 切2 o fh?？谄?p一二衛(wèi)- 15 0pswtr sfnciiai dtnsiiy &115 s 鼻與braqumcy e口由聞 10 po»wsp«trji densityijpham-)3i 1.5225fe 串 mty卜圖3-16調(diào)制過程功率譜變化。英6 0 硒 o.m m32 mmtime 4fke|k 也看po*ei spectral dqnmiy1142ff«4jei

27、k¥ iradsfi吟fuvtf印mu" 口白曲面神jo.s 1!14221 mf-imutncy iradu'smlj元/powur sp«£irjl mi.iyi聲3鐘-100bf>pqu?ncy jisdvw)fguecy lied或)q.s 11 £2 2s 3rnatyjancy 2uwt卜徜 對r?t5eh%hw圖3-17解調(diào)過程功率譜變化10d珈7d0砸圖3-18沒失真眼圖圖3-19有失真眼圖3.5.2模擬調(diào)制法極性調(diào)解法仿真分析圖3-12波形表示基帶信號經(jīng)過差分編碼再經(jīng)過極性變換然后與載波相乘得到2dpsk言號。圖

28、3-13波形表示2dpsk言號經(jīng)過信道再經(jīng)過帶通濾波器然后與載波相乘得到相乘后的信號。圖3-14波形表示相乘后的信號經(jīng)低通濾波器經(jīng)過判決器后得到解調(diào)差分信號再得到解調(diào)輸出信號。圖3-15波形是基帶信號和系統(tǒng)輸出信號的對比，考慮到有延時情況，從波形可以看出系統(tǒng)解調(diào)出了原基帶信號，所以系統(tǒng)設(shè)計(jì)成功。圖3-16從左到右分別表示基帶信號功率譜、2dpsk言號功率譜和信道輸出信號功率譜。圖3-17從左到右分別表示信道輸出信號經(jīng)過帶通濾波器后的功率譜、低通濾波器輸出信號功率譜、系統(tǒng)解調(diào)輸出信號功率譜。從功率譜變化可以看出幸好從低頻搬移到高頻傳輸然后又從高頻搬移到低頻，這些完全符合數(shù)字信號的調(diào)制解調(diào)原理。圖

29、3-18表示在噪聲simpletime為1時低通濾波器輸出的眼圖。圖3-19表示在噪聲simpletime為0.001時低通濾波器輸出的眼圖。成都學(xué)院(成都大學(xué))課程設(shè)計(jì)報告第4章2dpsk系統(tǒng)程序設(shè)計(jì)仿真4.1 matlab 程序clear,close allbit=1000;n=16;p=0.6；signal=rand(1,n)<=p;% receive=0; % j=1； while j<=n %產(chǎn)生n位隨機(jī)二進(jìn)制信號 y = rand(m,n)或y = rand(m n) 返回一個m x n的隨機(jī)矩陣差分編碼if signal(j)=receive(j)%=不等于bi=1;

30、elsebi=0;endreceive=receive bi;j=j+1；enddifference=receive(2:n+1);% m=0:1/bit:(n-1)/bit;figure(1)%subplot(3,1,1) %除去差分碼參考位創(chuàng)建圖形畫出基帶原碼stairs(m,signal),axis(0,n/bit,-0.5,1.5);%title( '基帶原碼'),xlabel( 'time/sec' ),ylabel( grid on%添加網(wǎng)格線subplot(3,1,2) %畫出差分碼stairs(m,difference),axis(0,n/bit

31、,-0.5,1.5);stairs函數(shù)繪制階梯狀圖axis函數(shù)設(shè)置坐標(biāo) 幅值)；title( '差分碼),xlabel( 'time/sec' ),ylabel( '幅值');grid onfor i=1:n;%將差分碼變成雙極性差分碼biploarif difference(i)=0biploar(i)=-1;elsebiploar(i)=1;endendt=linspace(0,16/1000,3200); % carrier=sin(2*pi*1500*t); % biploar=repmat(biploar,200,1);%載波時間切片在0到16

32、/1000 平均分為3200點(diǎn)載波頻率為1500復(fù)制平鋪biploar=reshape(biploar,1,numel(biploar);%reshape重新調(diào)整行數(shù)列數(shù)modulate=biploar*carrier; % numel計(jì)算矩陣元素個數(shù)，模擬調(diào)制，產(chǎn)生2dpsk調(diào)制波25subplot(3,1,3),plot(t,modulate),axis(0,n/bit,-1.5,1.5);title( '2dpsk 調(diào)制波),xlabel( 'time/sec' ),ylabel( '幅值');grid onfigure(2)%&&

33、&&&&&&&&&&&&&&&&&&&&&&& 模擬信道傳輸 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&channel=awgn(modulate,10);%awgn隨

34、機(jī)無線噪聲信噪比為 10dbsubplot(3,1,1),plot(t,channel)axis(0,n/bit,-1.5,1.5);title( '信道輸出),xlabel( 'time/sec' ),ylabel( '幅值)；grid on%&&&&&&&&&&&&&&&&&&&&&&&&& 相干解調(diào) &&&&&&

35、&&&&&&&&&&&&&&&&&&&&&&&&&&& coherent=channel*carrier;subplot(3,1,2)plot(t,coherent);axis(0,n/bit,-1.5,1.5);title( '相乘輸出'),xlabel( 'time/sec' ),ylabel( '幅值');grid on%&am

36、p;&&&&&&&&&&&&&&&&&&&&&&&&& 彳氐通橢圓濾波器 &&&&&&&&&&&&&&&&&&&&&&&&&&& fp=60;fs=150;fs=4000;rp=3;r

37、s=60;wp=fp/fs;ws=fs/fs;n,wn=buttord(wp,ws,rp,rs);% bz,az=butter(n,wn); % low-pass=filter(bz,az,coherent);% subplot(3,1,3) plot(t,low-pass),axis(0,n/bit,-1.5,1.5);計(jì)算濾波器階數(shù)和截止頻率計(jì)算數(shù)字濾波器系統(tǒng)函數(shù)的分子分母多項(xiàng)式系數(shù)y=filter (b,a,x ) b/a為濾波器系數(shù)，b為分子a為分母title( '低通輸出'),xlabel( 'time/sec' ),ylabel( '幅值&#

38、39;);grid on judge=;bridge=;%&&&&&&&&&&&&&&&&&&&&&&&&&& 抽樣判決&&&&&&&&&&&&&&&&&&&&&&&&&&&&

39、amp;&&&for b=0:200:3000;if b=0c=1;judge=(1&&low-pass(c)=0);elsebridge=1&&(low-pass(b)>0);endjudge=judge bridge;end m=0:1/bit:(n-1)/bit;%judge=1*low-pass>0);%&&&&&&&&&&&&&&&&&&&&&&

40、;&&& 繪制輸入輸出碼型圖以比較&&&&&&&&&&&&&&&&&&&&&figuer(3)subplot(4,1,1)stairs(m,signal,'r' ).axis(0,n/bit,-1.5,1.5);title( '基帶原碼'),xlabel( 'time/sec' ),ylabel( '幅值)；grid onsubplot(4,1,2)

41、stairs(m,difference,'r' ).axis(0,n/bit,-1.5,1.5);title( '差分碼),xlabel( 'time/sec' ),ylabel( '幅值');grid onsubplot(4,1,3)stairs(m,judge, 'r' ).axis(0,n/bit,-1.5,1.5);title( '抽樣判決'),xlabel( 'time/sec' ),ylabel( '幅值)； grid on%&&&&&&&&&&&&&&&&&&&&&&&&&& 碼型差分逆變換 &&&&&&&&&&&&&&&&&&&&&am