課程設(shè)計--基于matlab的數(shù)字頻帶通信系統(tǒng)仿真

1、<p><b>　　通信專業(yè)課程設(shè)計二</b></p><p>　　課 程 設(shè) 計（論 文）</p><p>　　設(shè)計(論文)題目：基于MATLAB的數(shù)字頻帶通信系統(tǒng)仿真</p><p>　　姓 名 </p><p>　　學(xué) 號 </p><p

2、>　　班 級 通信 </p><p>　　學(xué) 院 電子信息工程學(xué)院 </p><p>　　指導(dǎo)教師 </p><p>　　2013年 1 月 8 日</p><p>　　課程設(shè)計（論文）任務(wù)書</p><p>　　學(xué)院（直屬系）：電子信息工程學(xué)院

3、 時間： 2013年1月8日</p><p>　　基于MATLAB的數(shù)字頻帶通信系統(tǒng)仿真</p><p><b>　　摘 要</b></p><p>　　本論文研究了數(shù)字頻帶傳輸?shù)男阅堋@基本結(jié)構(gòu)發(fā)送濾波器、調(diào)制器、信道、解調(diào)器、接收濾波器等部分進行仿真。在分析數(shù)字頻帶傳輸基本原理的基礎(chǔ)上，從兩方面理論分析其傳輸系統(tǒng)的基本性能。<

4、;/p><p>　　第一個方面使用MATLAB中M語言完成QPSK的仿真，得出在加性高斯白噪聲的信道下，傳輸比特錯誤率，得出關(guān)系曲線，并將比特錯誤率與理論值相比較，分析其傳輸性能。</p><p>　　第二個方面是通過仿真得到QPSK和OQPSK信號經(jīng)過帶通后的波形。對比兩個信號的功率譜密度圖，理論分析其傳輸性能。</p><p>　　關(guān)鍵詞：仿真 MATLAB Q

5、PSK OQPSK</p><p><b>　　目 錄</b></p><p><b>　　摘 要</b></p><p>　　第1章 緒論……………………………………………………………………………………1</p><p>　　1.1數(shù)字頻帶傳輸系統(tǒng)的簡介……………………………………………………

6、…………1</p><p>　　1.2本論文的選題意義………………………………………………………………………2</p><p>　　第2章 數(shù)字信號的頻帶傳輸…………………………………………………………………4</p><p>　　2.1 數(shù)字頻帶傳輸系統(tǒng)………………………………………………………………………4</p><p>　　2.2 四

7、相移相鍵控（QPSK）…………………………………………………………………5</p><p>　　2.2.1 QPSK信號的產(chǎn)生……………………………………………………………………5</p><p>　　2.2.2 QPSK信號的平均功率譜密度………………………………………………………8</p><p>　　2.2.3 QPSK信號的接收及其平均誤比特率…………………

8、……………………………9</p><p>　　2.3偏移四相移相鍵控（OQPSK）…………………………………………………………11</p><p>　　2.3.1 OQPSK調(diào)制…………………………………………………………………………11</p><p>　　2.3.2 OQPSK信號的最佳接收及誤比特率………………………………………………13</p>

9、<p>　　第3章仿真及其結(jié)果分析……………………………………………………………………15</p><p>　　3.1 QPSK仿真………………………………………………………………………………15</p><p>　　3.2 OQPSK仿真………………………………………………………………………………16</p><p>　　第4章 結(jié)論………………………

10、……………………………………………………………18</p><p><b>　　參考文獻(xiàn)19</b></p><p><b>　　附 錄20</b></p><p><b>　　第1章 緒論</b></p><p>　　在現(xiàn)代信息化的社會里，通信是很重要的，尤其是移動通信，

11、它在現(xiàn)代化進程中起到的巨大作用是不容質(zhì)疑的。移動通信可以說從M·G·馬可尼（1987）所完成的無線通信實驗之日就產(chǎn)生了，當(dāng)時的試驗是在固定站與一艘拖船之間進行的，距離為18海里。移動通信的發(fā)展走過了一段很漫長的路程。通信信號處理方式也從模擬向數(shù)字通信轉(zhuǎn)變。同時通信系統(tǒng)已經(jīng)由單一的硬件系統(tǒng)轉(zhuǎn)向與軟件的結(jié)合。隨著通信技術(shù)的發(fā)展，信號處理方面硬件設(shè)計與專業(yè)軟件設(shè)計結(jié)合日趨緊密。比較優(yōu)秀而且得到廣大技術(shù)人員認(rèn)可的有MTALA

12、B。MATLAB是Matrix Labotratory(矩陣實驗室)的縮寫，最初由美國Cleve Moler博士在20世紀(jì)70年代末講授矩陣?yán)碚摵蛿?shù)據(jù)分析等課程編寫的軟件包Linkpack與Eispack組成，旨在使應(yīng)用人員免去大量重復(fù)的矩陣運算和基本數(shù)字運算等繁瑣的編程工作。MATLAB作為一種仿真平臺，提供了強大的工程，設(shè)計者可以自己設(shè)計函數(shù)，同時也可以調(diào)用Toolbox里的函數(shù)調(diào)用。</p><p>　　M

13、ATLAB等優(yōu)秀軟件使仿真技術(shù)得到很好的應(yīng)用。通過對通信過程的仿真，我們就可以在低成本的條件下檢測某一個方案是否能夠?qū)崿F(xiàn)，是否有更好的方案可以代替原來的方案，這樣對通信的研究就站在了一個更高的起點，使通信技術(shù)的發(fā)展日新月異，近十年來手機的普及率的迅速提高就從側(cè)面反映移動通信技術(shù)的發(fā)展。</p><p>　　移動通信中，數(shù)字調(diào)制解調(diào)技術(shù)是關(guān)鍵技術(shù)，其中數(shù)字調(diào)制信號具有數(shù)字通信的諸多優(yōu)點，在數(shù)字移動通信中，廣泛使用它

14、來傳送各種控制信息。</p><p>　　1.1 數(shù)字頻帶傳輸系統(tǒng)的簡介</p><p>　　人類的社會活動和發(fā)展是與信息及其傳遞——通信密切相關(guān)的。當(dāng)今世界已進入信息時代。通信已深入到社會各個領(lǐng)域，對人類的生活和社會的發(fā)展更加具有特殊重要的意義，對于遠(yuǎn)程通信來說，數(shù)字頻帶傳輸技術(shù)的發(fā)展尤為重要。</p><p>　　隨著數(shù)字技術(shù)的飛速發(fā)展與數(shù)字器件的廣泛使用數(shù)字信

15、號處理在通信系統(tǒng)中的應(yīng)用已經(jīng)越來越重要。數(shù)字信號傳輸系統(tǒng)分為基帶傳輸系統(tǒng)和頻帶傳輸系統(tǒng)。頻帶傳輸系統(tǒng)也叫數(shù)字調(diào)制系統(tǒng)，數(shù)字頻帶傳輸系統(tǒng)是一種利用調(diào)制器對傳輸信號進行頻率交換的傳輸方式，該系統(tǒng)對基帶信號進行調(diào)制，使其頻譜搬移到適合在信道(一般為帶通信道)上傳輸?shù)念l帶上。信號調(diào)制的目的是為了更好的適應(yīng)信號傳輸通道的頻率特性，傳輸信號經(jīng)過調(diào)制處理也能克服基帶傳輸同頻帶過寬的缺點，提高線路的利用率，一舉兩得。研究數(shù)字頻帶通信系統(tǒng)是研究數(shù)字通信的

16、基礎(chǔ)，現(xiàn)代大多數(shù)數(shù)字通信系統(tǒng)都是頻帶通信系統(tǒng)。數(shù)字頻帶傳輸在應(yīng)用上要比基帶傳輸更廣泛，應(yīng)用范圍更大，如：衛(wèi)星通信、移動通信、光纖通信、均是在所規(guī)定信道頻帶內(nèi)傳輸頻帶信號。</p><p>　　通信信號的調(diào)制識別成為研究熱點之一，國內(nèi)外都有相關(guān)方面的研究，并且取得很好的結(jié)果。80年代以來，通信和信號處理系統(tǒng)越來越復(fù)雜，各種新技術(shù)的發(fā)展對通信系統(tǒng)的實現(xiàn)起著重大的影響。通信系統(tǒng)復(fù)雜性的增加使得分析與設(shè)計所需的時間和費用

17、也迅速上升，為了節(jié)約人力、物力、財力和時間，就需要進行系統(tǒng)仿真。近10年來，數(shù)字通信系統(tǒng)得到了迅速的發(fā)展，尤其是人們對數(shù)字通信系統(tǒng)的需求，已超過現(xiàn)今技術(shù)所能支持?jǐn)?shù)據(jù)速率幾個數(shù)量級，而可利用的無線頻譜有限，只有使通信頻譜利用率得到顯著提高，才能滿足通信容量的要求。盡管通過增加發(fā)射功率和帶寬來提高系統(tǒng)輸出信噪比可以改善系統(tǒng)性能，但這與移動通信的發(fā)展目標(biāo)相違背。一個使數(shù)據(jù)速率增長得以實現(xiàn)的重要技術(shù)突破就是研究數(shù)字頻帶通信系統(tǒng)。</p&g

18、t;<p>　　過去，數(shù)字信號處理一直是用模擬設(shè)備完成，近代，數(shù)字計算機的出現(xiàn)和大規(guī)模集成技術(shù)的高度發(fā)展，為數(shù)字頻帶通信系統(tǒng)的處理提供了強有力的手段。在電子各個領(lǐng)域，正日益廣泛地用數(shù)字方法實現(xiàn)數(shù)字頻帶通信系統(tǒng)。</p><p>　　數(shù)字通信在其發(fā)展過程中表現(xiàn)出了強大的生命力，它沖破了傳統(tǒng)模擬通信方式的統(tǒng)治，逐步地發(fā)展、完善。隨著通信事業(yè)的發(fā)展，特別是各種寬帶傳輸技術(shù)(例如光纖傳輸、數(shù)字微波等)、綜合

19、業(yè)務(wù)數(shù)字網(wǎng)(ISDH)的實用化，全數(shù)字化的通信方式必將逐步取代模擬通信方式而得到蓬勃發(fā)展。數(shù)字調(diào)制傳輸在現(xiàn)在通信中發(fā)揮著越來越重要的作用。為使數(shù)字信號能在帶通信道中傳輸，必須用數(shù)字信號對載波進行調(diào)制，其調(diào)制方式與模擬信號調(diào)制相類似。根據(jù)數(shù)字信號控制波的參量不同也分為調(diào)幅、調(diào)頻和調(diào)相三種方式。因數(shù)字信號對載波參數(shù)的調(diào)制通常采用數(shù)字信號的離散值對載波進行鍵控，故這三種數(shù)字調(diào)制方式被稱為幅移鍵控（ASK）、頻移鍵控（FSK）和相移鍵控（PSK

20、）。基于頻帶系統(tǒng)仿真是通信電子類專業(yè)教學(xué)和科研經(jīng)常使用的一種方法。我們在計算機嘗試仿真，在計算機上顯示不同系統(tǒng)的工作原理進行波形觀察、頻譜分析和性能分析等，為數(shù)字頻帶傳輸系統(tǒng)以及通信系統(tǒng)設(shè)計和研究提供強有力的指導(dǎo)。</p><p>　　1.2 本論文的選題意義</p><p>　　當(dāng)要測試一個傳輸系統(tǒng)的性能時，往往都是把它做出來，再對其進行測量計算。這樣做的話，需要浪費大量時間。而我們現(xiàn)在

21、利用計算機，把需要設(shè)計的方案進行計算機仿真，模擬測試出來傳輸系統(tǒng)，最后計算出來和他相關(guān)的誤碼率，誤比特率等數(shù)據(jù)，得出這個傳輸系統(tǒng)性價比。MATLAB是一種功能強大，效率高，交互性好的數(shù)值計算和可視化計算機高級語言，它將數(shù)值的分析，矩陣計算，信號處理和圖形顯示有機地融合為一體，形成了一個極其方便，用戶界面有好的操作環(huán)境。隨著其自身版本的不斷提高，MATLAB得到了各個領(lǐng)域?qū)＜覍W(xué)者的廣泛關(guān)注，其強大的擴展功能可以為各個領(lǐng)域的應(yīng)用提供基礎(chǔ)。目

22、前，MATLAB已廣泛應(yīng)用于科學(xué)研究工程計算，仿真，自動控制等領(lǐng)域，是目前工程上流行最廣泛的科學(xué)語言。數(shù)字頻帶系統(tǒng)的應(yīng)用廣泛，所以研究數(shù)字頻帶通信系統(tǒng)并用軟件仿真，具有現(xiàn)實意義。</p><p>　　我首先要復(fù)習(xí)頻帶傳輸調(diào)制解調(diào)的基本原理，同時也要復(fù)習(xí)通信系統(tǒng)的主要組成部分，了解調(diào)制解調(diào)方式中最基礎(chǔ)的方法及數(shù)學(xué)原理。并對“通信”這個概念有個整體的理解，學(xué)習(xí)數(shù)字調(diào)制中誤碼率測試的標(biāo)準(zhǔn)及計算方法。同時還要復(fù)習(xí)隨機信號

23、中時域用自相關(guān)函數(shù)，頻域用功率譜密度來描述平穩(wěn)隨機過程的特性等基礎(chǔ)知識，來理解高斯信道中噪聲的表示方法，以便在編程中使用。 </p><p>　　在理解頻帶傳輸調(diào)制解調(diào)的基本原理后，使用MATLAB編程實現(xiàn)頻帶信號在高斯信道下傳輸，以及該方式的誤碼率測試。復(fù)習(xí)MATLAB編程的基礎(chǔ)知識，編程的常用算法以及使用MATLAB仿真系統(tǒng)的注意事項，并鍛煉自己的編程能力，通過編程完成QPSK調(diào)制解調(diào)系統(tǒng)的仿真，以及誤碼率測

24、試，并得出響應(yīng)波形。</p><p>　　本論文要針對數(shù)字頻帶傳輸?shù)男阅軄磉M行研究，在分析數(shù)字頻帶傳輸基本原理的基礎(chǔ)上，理論分析其傳輸系統(tǒng)的基本性能。通過使用MATLAB軟件利用M程序得到比較符合實際情況的高斯噪聲下QPSK傳輸?shù)恼`碼率曲線，QPSK和OQPSK信號經(jīng)過帶通后的波形和功率譜圖，理論分析其傳輸性能。</p><p>　　第2章 數(shù)字信號的頻帶傳輸</p><

25、;p>　　2.1數(shù)字頻帶傳輸系統(tǒng)</p><p>　　數(shù)字通信的基本特征是，它的消息或信號具有“離散”或“數(shù)字”的特征，從而使數(shù)字通信具有許多特殊的問題。</p><p>　　數(shù)字通信中存在以下突出問題：第一，數(shù)字信號傳輸時，信道噪聲或干擾所造成的差錯，原則上是可以控制的。這是通過所謂的差錯控制編碼來實現(xiàn)的。于是就需要在發(fā)送端增加一個編碼器，而在接收端相應(yīng)需要一個解碼器。第二，當(dāng)需要

26、實現(xiàn)保密通信時，可對數(shù)字基帶信號進行人為“擾亂”（加密），此時在接收端就必須進行解密。第三，由于數(shù)字通信傳輸?shù)氖且粋€接一個按一定節(jié)拍傳送的數(shù)字信號，因而接收端必須有一個與發(fā)送相同的節(jié)拍，否則，就會因收發(fā)步調(diào)不一致而造成混亂。另外，為了表述消息內(nèi)容，基帶信號都是按消息特征進行編組的，于是，在收發(fā)之間一組組的編碼的規(guī)律也必須一致。否則接收時消息的真正內(nèi)容將無法修復(fù)。在數(shù)字通信系統(tǒng)中，稱節(jié)拍一致位“位同步”或“碼元同步”，而稱編組一致為“群同

27、步”或“幀同步”，故數(shù)字通信中還必須有“同步”這個重要問題。</p><p>　　綜上所述，點對點的數(shù)字通信模型一般可用圖2.1表示。</p><p>　　圖2.1數(shù)字頻帶通信系統(tǒng)模型</p><p>　　圖中調(diào)制器/解調(diào)器、加密器/解密器、編碼器/譯碼器等環(huán)節(jié)，在具體通信系統(tǒng)中是否全部采用，這要取決于具體設(shè)計條件和要求。但在一個系統(tǒng)中，如果發(fā)端有調(diào)制/加密/編碼，

28、則收端必須有解調(diào)/解密/譯碼。通常把有調(diào)制器/解調(diào)器的數(shù)字通信系統(tǒng)稱為數(shù)字頻帶傳輸通信系統(tǒng)。</p><p>　　在實際通信中，不少信道都不能直接傳送基帶信號，必須用基帶信號對載波波形的某些參量進行控制，即對基帶信號進行調(diào)制。調(diào)制即是按照調(diào)制信號（基帶信號）的變化規(guī)律去改變載波某些參數(shù)的過程。調(diào)制不僅可以對被調(diào)制信號進行頻譜搬移和擴頻，而且對系統(tǒng)的傳輸有效性和傳輸可靠性有很大的影響。因此，調(diào)制方式往往決定了一個通

29、信系統(tǒng)的性能。在傳輸脈沖時，傳輸系統(tǒng)會產(chǎn)生噪聲，如果其幅度為反極性并超過脈沖的幅度，則產(chǎn)生誤碼。在以脈沖調(diào)制載波時，也有同樣現(xiàn)象產(chǎn)生。但是，在一定的噪聲情況下，根據(jù)對載波的調(diào)制方法和解調(diào)方法的不同，誤碼率也不同。一般，數(shù)字頻帶傳輸中，調(diào)制方式的選擇往往是頻帶利用率、誤比特率、(或S/N)和設(shè)備實現(xiàn)復(fù)雜性等因素綜合考慮的結(jié)果，必須根據(jù)具體使用條件進行比較才能做出判斷。</p><p>　　數(shù)字調(diào)制方式有三種：幅度鍵

30、控、頻移鍵控和相移鍵控。它們分別對應(yīng)于用正弦波的幅度、頻率和相位來傳遞數(shù)字基帶信號。當(dāng)調(diào)制信號為二進制數(shù)字信號時，該調(diào)制稱為二進制數(shù)字調(diào)制。而當(dāng)調(diào)制信號為多進制數(shù)字信號時，則稱為多進制數(shù)字調(diào)制。相移鍵控有很多方式，如PSK、DPSK、QPSK、OQPSK、pi/4-DPSK等。本設(shè)計重點對QPSK和OQPSK進行研究和仿真。</p><p>　　2.2 四相移相鍵控（QPSK）</p><p&

31、gt;　　2.2.1 QPSK信號的產(chǎn)生</p><p>　　四相移相鍵控(QPSK)又名四進制移相鍵控，該信號的正弦載波有4個可能的離散相位狀態(tài)，每個載波相位攜帶2個二進制符號，其信號表示式為：</p><p><b>　　式(2-1)</b></p><p>　　為四進制符號間隔，為正弦載波的相位，有4種可能狀態(tài)。</p>

32、<p>　　若，則為0、、、，此初始相位為0的QPSK的矢量圖于圖2.2（a）。</p><p>　　若，則為、、、，此初始相位為的矢量圖示于圖2.2（b）。</p><p>　　圖2.2 QPSK信號矢量圖</p><p>　　下面，著重分析圖2.2(a)的QPSK信號的產(chǎn)生及解調(diào)。將式2-1寫成：</p><p><b&g

33、t;　　式(2-2)</b></p><p>　　若為、、、，則；。于是，式2-2可寫成：</p><p>　　； 式(2-3)</p><p>　　根據(jù)式2-3可得到圖所示的正交調(diào)制框圖：</p><p>　　圖2.3產(chǎn)生QPSK信號的正交調(diào)制原理圖</p><p>　　從圖2.3看出，信息速

34、率為的二進制序列(取值為+1或-1)，串并變換后分成兩路速率減半的二進制序列，得到基帶信號波形I(t)及Q(t)，這兩路碼元在時間上是對齊的，稱這兩支路為同相支路及正交支路，將它們分別對正交載波及進行2PSK調(diào)制，再將這兩支路的2PSK信號相加即可得到QPSK信號。</p><p>　　k 1 2 3 4 5 6 7 8 9 10 11 12 13 14&l

35、t;/p><p>　　n +1 +1 -1 +1 +1 -1 +1 -1 +1 +1 -1 -1 -1 +1</p><p>　　2n +1 +1 -1 -1 +1 -1 +1</p><p>　　2n-1 +1 -1 +

36、1 +1 +1 -1 -1</p><p>　　圖2.4 QPSK串并變換及I(t)、Q(t)波形圖</p><p>　　從圖5-3及圖5-4看出，將二進制序列進行串并變換，其實質(zhì)是完成二進制與四進制符號的變換，稱并行支路所對應(yīng)的每一對比特及為雙比特碼元，其符號間隔為四進制符號間隔（），每一組雙比特碼元與四進制符號之一相對應(yīng)，而每個四進制符號()又與

37、QPSK信號的載波相位相對應(yīng)，稱QPSK的載波相位與它所攜帶的雙比特碼元之間的關(guān)系。用上述正交調(diào)制法所產(chǎn)生的QPSK的相位邏輯關(guān)系克通過表2-1作進一步解釋。</p><p>　　表2-1 QPSK信號載波相位與雙比特碼元的關(guān)系</p><p>　　從表2-1看出，QPSK的載波相位與雙比特碼之間的關(guān)系正好符合格雷碼的相位邏輯，即相鄰四進制符號所對應(yīng)的雙比特碼之間僅相差一個二進制信號。&l

38、t;/p><p>　　采用格雷碼的相位邏輯的優(yōu)點是：若QPSK信號在信道傳輸中受到加性噪聲干擾，在噪聲不太大時，接收到的載波相位有可能接近相鄰的載波相位，在解調(diào)時，會發(fā)生錯判為相鄰四進制符號的現(xiàn)象，在從四進制符號譯為雙比特二進制碼時，若采用格雷碼邏輯關(guān)系，則在兩個比特符號中僅錯一個比特符號，這樣可減小誤比特率，所以希望QPSK的相位邏輯符合格雷編碼的關(guān)系。</p><p>　　2.2.2 QP

39、SK信號的平均功率譜密度</p><p>　　由于QPSK信號時由兩正交載波的2PSK線性疊加而成，所以QPSK信號的平均功率譜密度是同相支路及正交支路2PSK信號平均功率譜密度的線性疊加。</p><p>　　2PSK信號的平均功率譜密度公式如式2-4所示。</p><p><b>　　式(2-4)</b></p><p&

40、gt;<b>　　式(2-5)</b></p><p>　　假設(shè)QPSK與2PSK的信息速率及發(fā)送信號的平均功率均一致，則QPSK信號每個支路的2PSK信號的振幅為，符號間隔，QPSK的總功率譜是兩支路功率譜之和，于是，得到QPSK功率譜如式2-5所示。</p><p>　　根據(jù)式2-4及式2-5，畫出了2PSK信號及QPSK信號的平均功率譜密度，如圖2-5所示。&l

41、t;/p><p>　　圖2-5表示在2PSK及QPSK的二進制信息速率相同時，QPSK信號的平均功率譜密度的主瓣寬度是2PSK平均功率譜主瓣寬度的一半。</p><p>　　圖2.5 在給定信息速率為條件下，2PSK及QPSK雙邊功率譜密度</p><p>　　2.2.3 QPSK信號的接收及其平均誤比特率</p><p>　?。?）在加性高斯白

42、噪聲信道條件下QPSK最佳接收</p><p>　　在加性高斯白噪聲干擾的信道條件下，利用匹配濾波器或相關(guān)解調(diào)器進行QPSK的解調(diào)，圖2.6為QPSK最佳接收框圖。</p><p>　　圖2.6 QPSK信號的最佳接收框圖</p><p>　　由于QPSK信號可看為同相及正交支路2PSK的疊加，所以在解調(diào)時可對兩路信號分別進行2PSK的解調(diào)，然后進行并串變換，得到所

43、傳輸?shù)臄?shù)據(jù)。</p><p>　　計算QPSK解調(diào)的誤比特率有兩種方法：一是先計算誤符率(平均錯判四進制信號的概率)，然后再根據(jù)誤符率計算從四進制符號譯為二進制符號的誤比特率；另一種計算方法是沿用2PSK匹配濾波器解調(diào)的誤比特率計算公式。下面，采用第二種方法來進行計算。</p><p>　　在加性高斯白噪聲信道下，2PSK最佳接收的平均誤比特率為：</p><p>

44、<b>　　式(2-6)</b></p><p>　　對于QPSK而言，在QPSK與2PSK的輸入二進制信息速率相同，二者的發(fā)送功率相同，加性噪聲的單邊功率譜密度相同的條件下，QPSK與2PSK的平均誤比特率是相同的。</p><p>　　在給定二進制信息速率的條件下，QPSK的同相支路及正交支路的四進制符號速率是二進制信息速率的一半，即。在給定信號總發(fā)送功率的條件下

45、，QPSK同相支路或正交支路的信號功率是總發(fā)送功率的一半。于是，得到I支路及Q支路的平均錯判概率為：</p><p><b>　　式(2-7)</b></p><p>　　其中，是平均比特能量。</p><p>　　由于QPSK發(fā)短信源輸出的二進制符號“1”和“0”等概出現(xiàn)，二進制碼元經(jīng)串并變換后在同相支路及正交支路也是等概分布的，所以在收端的

46、同相及正交支路解調(diào)的輸出經(jīng)并串變換后的數(shù)據(jù)，其總的平均誤比特率與I支路及Q支路的平均誤判概率是相同的，即：</p><p><b>　　式(2-8)</b></p><p>　　其中，及分別是總的二進制碼元出現(xiàn)在I支路及Q支路的概率，，因而QPSK的平均誤比特率為：</p><p><b>　　式(2-9)</b><

47、/p><p>　?。?）在理想限帶及加性高斯白噪聲信道條件下QPSK最佳接收</p><p>　　在理想限帶及加性高斯白噪聲加性高斯白噪聲信道條件下的QPSK最佳頻帶傳輸系統(tǒng)的框圖如5-7所示，他在實際限帶通信系統(tǒng)中獲廣泛應(yīng)用。</p><p>　　圖2.7在理想限帶及加性高斯白噪聲信道條件下的QPSK最佳頻帶傳輸系統(tǒng)</p><p>　　若發(fā)端

48、信源的“1”和“0”等概率出現(xiàn)，其最佳接收的平均誤比特率為：</p><p><b>　　式(5-10)</b></p><p>　　綜上所述，將QPSK與2PSK相比較，在兩者的信息速率、信號發(fā)送功率、噪聲功率譜密度相同的條件下，QPSK與2PSK的平均誤比特率是相同的，而QPSK的功率譜主瓣寬度比2PSK的窄一半。</p><p>　　2.

49、3偏移四相移相鍵控（OQPSK）</p><p>　　在QPSK數(shù)字調(diào)制中，若將二進制雙極性不歸零矩陣脈沖序列串并變換后在進行正交載波調(diào)制，則所得到的QPSK信號包絡(luò)是恒定的。在實際數(shù)字通信中，往往要求信道帶寬是有限的，為了對QPSK信號帶寬進行限制，經(jīng)常在QPSK數(shù)字調(diào)制器中，先將基帶雙極性矩陣不歸零脈沖序列經(jīng)過基帶成形濾波器限帶，再進行正交。載波調(diào)制，將限帶的基帶信號功率譜搬移到載頻上，成為限帶的QPSK信號

50、，但此時的限帶QPSK信號包絡(luò)不再恒定，并且在相鄰四進制符號的載波相位發(fā)生相移突變處，會出現(xiàn)包絡(luò)為零的現(xiàn)象，如圖2.12所示。</p><p>　　若將此包絡(luò)起伏很大的限帶QPSK信號再進行硬限幅或非線性功率放大，雖然此已調(diào)信號包絡(luò)的起伏可減弱，但卻使非線性功率放大后的信號功率譜旁瓣增生，頻率又被擴展，其旁瓣將會干擾鄰近頻道的信號，這是不希望出現(xiàn)的現(xiàn)象。因而對限帶QPSK信號的功率放大，只能采用線性功放，但線性功

51、放的功率轉(zhuǎn)換效率低，為此人們要探討能應(yīng)用非線性功率放大的限帶數(shù)字調(diào)制方式，其中之一是QPSK的修正形式：偏移四相移相鍵控OQPSK(或稱為參差QPSK)的調(diào)制方式。</p><p>　　2.3.1 OQPSK調(diào)制</p><p>　　OQPSK是在QPSK基礎(chǔ)上發(fā)展起來的一種恒包絡(luò)數(shù)字調(diào)制技術(shù)。這里，所謂恒包絡(luò)技術(shù)是指已調(diào)波的包絡(luò)保持為恒定，它與多進制調(diào)制是從不同的兩個角度來考慮調(diào)制技術(shù)的

52、。恒包絡(luò)技術(shù)所產(chǎn)生的已調(diào)波經(jīng)過發(fā)送帶限后，當(dāng)通過非線性部件時，只產(chǎn)生很小的頻譜擴展。這種形式的已調(diào)波具有兩個主要特點，其一是包絡(luò)恒定或起伏很??；其二是已調(diào)波頻譜具有高頻快速滾降特性，或者說已調(diào)波旁瓣很小，甚至幾乎沒有旁瓣。OQPSK的調(diào)制框圖如圖5-13所示。</p><p>　　將雙極性二進制序列串并變換后成為速率減半的同相支路及正交支路序列，其符號間隔。它與QPSK不同的是，他的正交支路碼元與同相支路碼元在時

53、間上要偏移一個比特間隔。同相及正交序列經(jīng)過成形濾波后再分別進行正交載波調(diào)制，其中I(t)及Q(t)的基帶波形如圖2.14所示。</p><p>　　圖2.13 產(chǎn)生OQPSK信號的框圖</p><p>　　k 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15</p><p>　　+1

54、+1 +1 +1 -1 -1 +1 +1 -1 -1 +1 -1 +1 +1 +1 -1</p><p>　　圖2.14 OQPSK的同相及正交支路基帶波形</p><p>　　圖2.14中的OQPSK信號表示為：</p><p><b>　　式(2-11)</b></p><p><b

55、>　　其中，</b></p><p><b>　　式(2-12)</b></p><p><b>　　式(2-13)</b></p><p>　　是成形濾波器沖擊響應(yīng)，若是不歸零矩陣脈沖，可表示為：</p><p><b>　　式(2-14)</b></

56、p><p>　　關(guān)于OQPSK及QPSK信號波形示于圖2.15及圖2.16。</p><p>　　圖2.15 OQPSK信號波形</p><p>　　圖2.16 QPSK信號波形</p><p>　　從圖2.15及圖2.16看出，QPSK信號在每隔時間，其信號載波相位可能發(fā)生的相位突變現(xiàn) OQPSK信號在每隔時間，其信號載波相位有可能發(fā)生相位變化

57、，不會發(fā)生的相位突變現(xiàn)象，因而此限帶的OQPSK信號的包絡(luò)起伏小，限帶OQPSK信號包絡(luò)的最大值與最小值之比約為，該信號經(jīng)線性功放后，不會引起功率譜旁瓣有大的增生，所以它適于在限帶非線性信道中使用。</p><p>　　2.3.2 OQPSK信號的最佳接收及誤比特率</p><p>　　圖2.17 在加性白高斯噪聲干擾下OQPSK的最佳接收框圖</p><p>　　

58、（1）OQPSK信號的平均功率譜密度</p><p>　　由于OQPSK也可以看為是由同相及正交支路的2PSK信號的疊加，所以O(shè)QPSK的功率譜與QPSK的相同。</p><p>　?。?）OQPSK的最佳接收及其誤比特率</p><p>　　在加性白高斯噪聲信道條件下，OQPSK最佳接收框圖如圖2.17所示。</p><p>　　OQPSK

59、最佳接收的平均誤比特率與QPSK相同：</p><p><b>　　式(2-15)</b></p><p>　　第3章 仿真及其結(jié)果分析</p><p>　　3.1 QPSK仿真</p><p>　　圖3.1 理想信道下的信號及對應(yīng)的功率譜密度</p><p>　　圖3.2 高斯噪聲曲線，高

60、斯信道下的調(diào)制信號，調(diào)制輸出功率譜密度曲線</p><p>　　由圖2.8上可知，完成了QPSK信號在理想信道上的調(diào)制，傳輸，解調(diào)的過程，由于調(diào)制過程中加進了載波，因此調(diào)制信號的功率譜密度會發(fā)生變化。并且可以看出調(diào)制解調(diào)的結(jié)果沒有誤碼。</p><p>　　由圖2.9，圖2.10可知，在高斯噪聲的影響下，調(diào)制信號的波形發(fā)生了明顯的變化，其功率譜密度函數(shù)相對于圖5-8中的調(diào)制信號的功率譜密度

61、只發(fā)生了微小的變化，原因在于高斯噪聲是一個均值為0的白噪聲，在各個頻率上其功率是均勻的，因此此結(jié)果是真確的。星座圖反映可接收信號在高斯噪聲的影響下發(fā)生了誤碼，但是大部分還是保持了原來的特性。</p><p>　　圖3.3 QPSK的星座圖</p><p>　　圖3.4 高斯信道和瑞利衰落信道的誤碼率對比</p><p>　　由圖2.11可知實際通過高斯信道下的信號誤

62、比特率和理論值基本相似。瑞利衰落信道下的誤碼率比高斯信道下的誤碼率高。</p><p>　　3.2 OQPSK仿真</p><p>　　圖3.5 QPSK及OQPSK經(jīng)帶通濾波器后波形及功率譜密度</p><p>　　經(jīng)過帶通濾波器后，QPSK信號的包絡(luò)起伏明顯變大，而OQPSK信號的包絡(luò)起伏相對小得多，再經(jīng)過非線性功率放大器后，此時QPSK信號的功率譜旁瓣比OQP

63、SK信號大。因此OQPSK信號性能優(yōu)于QPSK信號。</p><p><b>　　第4章 結(jié)論</b></p><p>　　本論文的題目是基于MATLAB數(shù)字頻帶通信系統(tǒng)仿真。在這次畢業(yè)設(shè)計中，我在老師的指導(dǎo)下，通過由不能熟練應(yīng)用MATLAB這款軟件，到自學(xué)其中的知識，并且在網(wǎng)絡(luò)上查找窗口中的界面的介紹，結(jié)合老師所給的參考的資料書目內(nèi)容，自學(xué)這款軟件，并獨立完成全部的

64、仿真的內(nèi)容。</p><p>　　在這次設(shè)計中，我仿真出了QPSK在加性高斯白噪聲下的最佳接收，比較并分析了QPSK及OQPSK的性能。仿真后的得到的QPSK在加性高斯白噪聲下的誤比特率與理論值接近。OQPSK的性能優(yōu)于QPSK，結(jié)果與預(yù)想中的一樣。</p><p>　　我通過自己獨立思考，統(tǒng)籌課本的知識點，然后一步一步著手實踐，最后完成了全部報告。</p><p>

65、;　　在數(shù)字頻帶通信系統(tǒng)中，QPSK調(diào)制解調(diào)技術(shù)已經(jīng)比較成熟，在實際應(yīng)用中也發(fā)揮著極為重要的作用。而在通信技術(shù)高速發(fā)展的今天，單一的數(shù)字調(diào)制技術(shù)已經(jīng)無法滿足通信的要求，根據(jù)不同通信方式而采取不同的調(diào)制方法，有時也會出現(xiàn)混合調(diào)制的方式，因此對基本的QPSK的研究與學(xué)習(xí)將為其他復(fù)雜的調(diào)制方法提供扎實的基礎(chǔ)。</p><p><b>　　參考文獻(xiàn)</b></p><p>　

66、　[1] 樊昌信，詹道庸，徐炳祥，吳成柯．通信原理[M]．北京：國防工業(yè)出版社，1980：171~260．</p><p>　　[2] 薛鵬騫．電子與通信電路計算機仿真[M]．北京：煤炭工業(yè)出版社，2003：113~144．</p><p>　　[3] 鄧華．MATLAB通信仿真及應(yīng)用實例詳解[M]．北京：人民郵電出版社，2003：99~115．</p><p>　

67、　[4] 陳懷琛．MATLAB及其在理工課程中的應(yīng)用指南[M]．西安：西安電子科技大學(xué)出版社，2002：146~180．</p><p>　　[5] 張謙．通信系統(tǒng)中MATLAB基礎(chǔ)與仿真應(yīng)用[M]．西安：西安電子科技大學(xué)出版社，2010：85~109．</p><p>　　[6] 孫屹．MATLAB通信仿真開發(fā)手冊[M]．北京：國防工業(yè)出版社，2005：138~168．</p>

68、;<p>　　[7] 王立寧，樂光新，詹菲．MATLAB與通信仿真[M]．北京：人民郵電出版社，2000：74~86．</p><p>　　[8] 約翰·G·普羅克斯．現(xiàn)代通信系統(tǒng)的使用MATLAB[M]．西安：西安交通學(xué)出版社，2001：95~121．</p><p>　　[9] 陳萍．現(xiàn)代通信實驗系統(tǒng)的計算機仿真[M]．北京：國防工業(yè)出版社，2003：

69、171~260．</p><p>　　[10] 孫亮．MATLAB語言與控制系統(tǒng)仿真[M]．北京：北京工業(yè)大學(xué)出版社，2001：64~83．</p><p>　　[11] William A Shay．Understanding Data Communications and Networks[M]：Wadsworth Publishing Company，1998：162~164．<

70、;/p><p>　　[12] William Stallings．Data and Computer Communications，Seventh Edition[M]: Prentice Hall，2003：201~207．</p><p><b>　　附錄</b></p><p>　　MATLAB仿真程序：</p><p&g

71、t;<b>　　一．QPSK仿真</b></p><p><b>　　clear;</b></p><p><b>　　T=1; </b></p><p><b>　　fc=10/T; </b></p><p><b>　　ml=2; </

72、b></p><p><b>　　nb=100; </b></p><p>　　delta_T=T/200; </p><p>　　fs=1/delta_T; </p><p><b>　　SNR=0; </b></p><p>　　t=0:delta_T:nb*T-d

73、elta_T; </p><p>　　N=length(t); </p><p>　　data=randn(1,nb)>0.5; </p><p>　　datanrz=data.*2-1; </p><p>　　data1=zeros(1,nb/delta_T); </p><p>　　for q=1:nb<

74、;/p><p>　　data1((q-1)/delta_T+1:q/delta_T)=datanrz(q); </p><p><b>　　end </b></p><p>　　% 將基帶信號變換成對應(yīng)波形信號</p><p>　　data0=zeros(1,nb/delta_T); </p><p>

75、;　　for q=1:nb</p><p>　　data0((q-1)/delta_T+1:q/delta_T)=data(q); % </p><p><b>　　end </b></p><p><b>　　% 發(fā)射的信號</b></p><p>　　data2=abs(fft(data1));

76、</p><p>　　% 串并轉(zhuǎn)換，將奇偶位數(shù)據(jù)分開</p><p>　　idata=datanrz(1:ml:(nb-1)); </p><p>　　qdata=datanrz(2:ml:nb);</p><p>　　% QPSK信號的調(diào)制</p><p>　　ich=zeros(1,nb/delta_T/2);

77、 </p><p>　　for i=1:nb/2</p><p>　　ich((i-1)/delta_T+1:i/delta_T)=idata(i);</p><p><b>　　end</b></p><p>　　for ii=1:N/2</p><p>　　a(ii)=sqrt(2/T)*c

78、os(2*pi*fc*t(ii)); </p><p><b>　　end</b></p><p>　　idata1=ich.*a; </p><p>　　qch=zeros(1,nb/2/delta_T);</p><p>　　for j1=1:nb/2</p><p&g

79、t;　　qch((j1-1)/delta_T+1:j1/delta_T)=qdata(j1);</p><p><b>　　end</b></p><p>　　for jj=1:N/2</p><p>　　b(jj)=sqrt(2/T)*sin(2*pi*fc*t(jj));</p><p><b>　　end

80、</b></p><p>　　qdata1=qch.*b; </p><p>　　s=idata1+qdata1; </p><p>　　ss=abs(fft(s)); </p><p><b>　　% 瑞利衰落信道</b></p><p>　　ray_ich=

81、raylrnd(0.8,1,nb/2/delta_T);</p><p>　　ray_qch=raylrnd(0.8,1,nb/2/delta_T);</p><p>　　Ray_idata=idata1.*ray_ich;</p><p>　　Ray_qdata=qdata1.*ray_qch;</p><p>　　Ray_s=Ray_id

82、ata+Ray_qdata;</p><p>　　%%%%%%%%%%%%%%%%%%%%%%%%%%</p><p><b>　　% 高斯信道 </b></p><p>　　s1=awgn(s,SNR); % 通過高斯信道之后的信號</p><p>　　s11=abs(fft(s1));

83、 % 快速傅里葉變換得到頻譜 </p><p>　　s111=s1-s; % 高斯噪聲曲線</p><p>　　%%%%%%%%%%%%%%%%%%%%%%%%%%</p><p>　　Awgn_s=awgn(Ray_s,SNR); % 通過高斯信道再通過瑞利衰落信道</p&g

84、t;<p>　　% QPSK 解調(diào)部分</p><p>　　% 解調(diào)部分（高斯信道）</p><p>　　idata2=s1.*a; % 這里面其實隱藏了一個串并轉(zhuǎn)換的過程</p><p>　　qdata2=s1.*b; % 對應(yīng)的信號與正余弦信號相乘</p><p>

85、　　idata3=zeros(1,nb/2); % 建立1*nb數(shù)組，以存放解調(diào)之后的信號</p><p>　　qdata3=zeros(1,nb/2);</p><p>　　% 抽樣判決的過程，與0作比較，data>=0,則置1，否則置0</p><p>　　for n=1:nb/2</p><p>　　% A1

86、(n)=sum(idata2((n-1)/delta_T+1:n/delta_T));</p><p>　　if sum(idata2((n-1)/delta_T+1:n/delta_T))>=0</p><p>　　idata3(n)=1;</p><p>　　else idata3(n)=0; </p><p><

87、b>　　end</b></p><p>　　% A2(n)=sum(qdata2((n-1)/delta_T+1:n/delta_T));</p><p>　　if sum(qdata2((n-1)/delta_T+1:n/delta_T))>=0</p><p>　　qdata3(n)=1;</p><p>

88、　　else qdata3(n)=0;</p><p><b>　　end</b></p><p><b>　　end </b></p><p>　　% 為了顯示星座圖,將信號進行處理</p><p>　　idata4=zeros(1,nb/2);</p><p>　　qd

89、ata4=zeros(1,nb/2);</p><p>　　for n=1:nb/2</p><p>　　Awgn_ichsum(n)=sum(idata2((n-1)/delta_T+1:n/delta_T))*delta_T;</p><p>　　if Awgn_ichsum(n)>=0</p><p>　　idata4(n)=1;

90、</p><p>　　else idata4(n)=0;</p><p><b>　　end</b></p><p>　　Awgn_qchsum(n)=sum(qdata2((n-1)/delta_T+1:n/delta_T))*delta_T;</p><p>　　if Awgn_qchsum(n)>=0<

91、/p><p>　　qdata4(n)=1;</p><p>　　else qdata4(n)=0;</p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　% 將判決之后的數(shù)據(jù)存放進數(shù)組</p><p>

92、　　demodata=zeros(1,nb);</p><p>　　demodata(1:ml:(nb-1))=idata3; % 存放奇數(shù)位</p><p>　　demodata(2:ml:nb)=qdata3; % 存放偶數(shù)位</p><p>　　%為了顯示，將它變成波形信號（即傳輸一個1代表單位寬度的高電平）</p><p>　

93、　demodata1=zeros(1,nb/delta_T); % 創(chuàng)建一個1*nb/delta_T的零矩陣</p><p>　　for q=1:nb</p><p>　　demodata1((q-1)/delta_T+1:q/delta_T)=demodata(q); % 將極性碼變成對應(yīng)的波形信號</p><p><b>　　end </b

94、></p><p><b>　　% 累計誤碼數(shù)</b></p><p>　　% abs(demodata-data)求接收端和發(fā)射端</p><p>　　% 數(shù)據(jù)差的絕對值，累計之后就是誤碼個數(shù)</p><p>　　Awgn_num_BER=sum(abs(demodata-data)) </p>&

95、lt;p>　　%%%%%%%%%%%%%%%%%%%</p><p>　　% 解調(diào)部分（瑞利+高斯）</p><p>　　Ray_idata2=Ray_s.*a; </p><p>　　Ray_qdata2=Ray_s.*b; </p><p>　　% Ray_idata3=ze

96、ros(1,nb/2); </p><p>　　% Ray_qdata3=zeros(1,nb/2);</p><p>　　% 抽樣判決的過程，與0作比較，data>=0,則置1，否則置0</p><p>　　% for n=1:nb/2</p><p>　　% if Ray_sum(Ray_idata2((n-1)/delta

97、_T+1:n/delta_T))>=0</p><p>　　% Ray_idata3(n)=1;</p><p>　　% else Ray_idata3(n)=0; </p><p><b>　　% end</b></p><p>　　% if Ray_sum(Ra

98、y_qdata2((n-1)/delta_T+1:n/delta_T))>=0</p><p>　　% Ray_qdata3(n)=1;</p><p>　　% else Ray_qdata3(n)=0;</p><p><b>　　% end</b></p><p><b>

99、;　　% end </b></p><p>　　% 為了顯示星座圖,將信號進行處理</p><p>　　Ray_idata4=zeros(1,nb/2);</p><p>　　Ray_qdata4=zeros(1,nb/2);</p><p>　　for n=1:nb/2</p><p>　　Ray_ic

100、hsum(n)=sum(idata2((n-1)/delta_T+1:n/delta_T))*delta_T;</p><p>　　if Ray_ichsum(n)>=0</p><p>　　Ray_idata4(n)=1;</p><p>　　else Ray_idata4(n)=0;</p><p><b>　　end&l

101、t;/b></p><p>　　Ray_qchsum(n)=sum(qdata2((n-1)/delta_T+1:n/delta_T))*delta_T;</p><p>　　if Ray_qchsum(n)>=0</p><p>　　Ray_qdata4(n)=1;</p><p>　　else Ray_qdata4(n)=0;

102、</p><p><b>　　end</b></p><p><b>　　end </b></p><p>　　% 將判決之后的數(shù)據(jù)存放進數(shù)組</p><p>　　Ray_demodata=zeros(1,nb);</p><p>　　Ray_demodata(1:ml:(n

103、b-1))=Ray_idata4; % 存放奇數(shù)位</p><p>　　Ray_demodata(2:ml:nb)=Ray_qdata4; % 存放偶數(shù)位</p><p>　　%為了顯示，將它變成波形信號（即傳輸一個1代表單位寬度的高電平）</p><p>　　Ray_demodata1=zeros(1,nb/delta_T); % 創(chuàng)建一個1*nb/

104、delta_T的零矩陣</p><p>　　for q=1:nb</p><p>　　Ray_demodata1((q-1)/delta_T+1:q/delta_T)=Ray_demodata(q); % 將極性碼變成對應(yīng)的波形信號</p><p><b>　　end </b></p><p><b>　　%

105、累計誤碼數(shù)</b></p><p>　　% abs(demodata-data)求接收端和發(fā)射端</p><p>　　% 數(shù)據(jù)差的絕對值，累計之后就是誤碼個數(shù)</p><p>　　Ray_num_BER=sum(abs(Ray_demodata-data)) </p><p>　　%%%%%%%%%%%%%%%%%%%%%%%%%

106、%%%%%%%%%%%%%%%%%%%%%%%%% </p><p>　　%% 誤碼率計算</p><p>　　%% 調(diào)用了cm_sm32();和cm_sm33(）函數(shù)</p><p>　　%%聲明： 函數(shù)聲明在另外倆個M文件中</p><p>　　%%作用： cm_sm32()用于瑞利信道誤碼

107、率的計算</p><p>　　%% cm_sm33()用于高斯信道誤碼率的計算</p><p>　　%% ecoh on/off 作用在于決定是否顯示指令內(nèi)容</p><p>　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</p><p>　　SNRindB1=0:1:8;&

108、lt;/p><p>　　SNRindB2=0:0.1:8;</p><p><b>　　% 瑞利衰落信道 </b></p><p>　　for i=1:length(SNRindB1),</p><p>　　[pb,ps]=cm_sm32(SNRindB1(i)); % 比特誤碼率</p>&l

109、t;p>　　smld_bit_ray_err_prb(i)=pb;</p><p>　　smld_symbol_ray_err_prb(i)=ps;</p><p>　　disp([ps,pb]);</p><p>　　echo off; </p><p><b>　　end;</b></p>&l

110、t;p><b>　　% 高斯信道 </b></p><p><b>　　echo on;</b></p><p>　　for i=1:length(SNRindB1),</p><p>　　[pb1,ps1]=cm_sm33(SNRindB1(i)); </p><p>

111、;　　smld_bit_awgn_err_prb(i)=pb1;</p><p>　　smld_symbol_awgn_err_prb(i)=ps1;</p><p>　　disp([ps1,pb1]);</p><p><b>　　echo off;</b></p><p><b>　　end;</b&

112、gt;</p><p><b>　　% 理論曲線</b></p><p><b>　　echo on;</b></p><p>　　for i=1:length(SNRindB2),</p><p>　　SNR=exp(SNRindB2(i)*log(10)/10);

113、 % 信噪比</p><p>　　theo_err_awgn_prb(i)=0.5*erfc(sqrt(SNR)); % 高斯噪聲理論誤碼率</p><p>　　theo_err_ray_prb(i)=0.5*(1-1/sqrt(1+1/SNR)); % 瑞利衰落信道理論誤碼率</p><p><b>　　echo off;</

114、b></p><p><b>　　end;</b></p><p>　　%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</p><p>　　h = spectrum.welch; % 類似于C語言的宏定義，方便以下的調(diào)用 </p><p>　　%%%%%%%%%%%

115、%%%%%%%%%%%%%%%%%%%%%%%%%</p><p>　　% 第一部分（理想）</p><p><b>　　figure(1)</b></p><p>　　subplot(212);</p><p>　　psd(h,data0,'fs',fs),title('基帶信號功率譜密度

116、9;);</p><p>　　subplot(211);</p><p>　　plot(data0),title('基帶信號');</p><p>　　axis([0 20000 -2 2]);</p><p><b>　　figure(2)</b></p><p>　　subpl

117、ot(212);</p><p>　　psd(h,s,'fs',fs),title('調(diào)制信號功率譜密度');</p><p>　　subplot(211);</p><p>　　plot(s),title('調(diào)制信號');</p><p>　　axis([0 500 -3 3]);</p

118、><p><b>　　figure(3)</b></p><p>　　subplot(212);</p><p>　　psd(h,demodata1,'fs',fs),title('解調(diào)輸出功率譜密度');</p><p>　　subplot(211);</p><p>