psk調(diào)制與解調(diào)武漢理工

1、課程設(shè)計任務(wù)書學(xué)生姓名： 王xx 專業(yè)班級： 通信1004班 指導(dǎo)教師： 許建霞 工作單位： 信息工程學(xué)院 題 目: 設(shè)計一個4PSK調(diào)制解調(diào)系統(tǒng)初始條件：本設(shè)計基于數(shù)字信號處理技術(shù)基礎(chǔ)實驗，通過自行設(shè)計程序并在電腦上利用MATLAB軟件進(jìn)行仿真。要求完成的主要任務(wù): （包括課程設(shè)計工作量及其技術(shù)要求，以及說明書撰寫等具體要求）1) 4PSK信號波形的載頻和相位參數(shù)應(yīng)隨機置或者可有幾組參數(shù)組合供選擇2) 系統(tǒng)中要求加入高斯白噪聲3) 4PSK解調(diào)方框圖采用相干接收形式4) 分析誤碼率參考書目:1 謝自美.電子線路設(shè)計·實驗·測試(第三版).武漢：華中科技大學(xué)出版社2 康華光

2、. 電子技術(shù)基礎(chǔ)模擬部分.高等教育出版社，20053 康華光. 電子技術(shù)基礎(chǔ) 數(shù)字部分.高等教育出版社，20054 樊昌信. 通信原理（第五版）.北京：國防工業(yè)出版社，2005時間安排：第1周，安排任務(wù)（鑒主15樓實驗室）第1-17周，仿真設(shè)計（鑒主13樓計算機實驗室）第18周，完成（答辯，提交報告，演示）指導(dǎo)教師簽名: 年 月 日系主任簽名: 年 月 日目 錄摘 要3Abstract41引言51.1 背景介紹51.2 設(shè)計要求52 4PSK調(diào)制解調(diào)的基本原理62.1 2PSK數(shù)字調(diào)制原理62.2 4PSK的調(diào)制和解調(diào)73 4PSK調(diào)制解調(diào)系統(tǒng)仿真103.1 MATLAB軟件介紹103.2 2

3、PSK調(diào)制解調(diào)系統(tǒng)仿真113.3 4PSK調(diào)制解調(diào)系統(tǒng)仿真124 4PSK誤碼率分析154.1 4PSK誤碼率的計算154.2 4PSK誤碼率的仿真165總結(jié)17參考文獻(xiàn)1826摘 要在數(shù)字信號的調(diào)制方式中4PSK是目前最常用的一種數(shù)字信號調(diào)制方式,它具有較高的頻譜利用率、較強的抗干擾性、在電路上實現(xiàn)也較為簡單。調(diào)制技術(shù)是通信領(lǐng)域里非常重要的環(huán)節(jié)，一種好的調(diào)制技術(shù)不僅可以節(jié)約頻譜資源而且可以提供良好的通信性能。4PSK調(diào)制是一種具有較高頻帶利用率和良好的抗噪聲性能的調(diào)制方式，在數(shù)字移動通信中已經(jīng)得到了廣泛的應(yīng)用。本次設(shè)計在理解4PSK調(diào)制解調(diào)原理的基礎(chǔ)上應(yīng)用MATLAB語言來完成仿真，仿真出了

4、4PSK的調(diào)制以及解調(diào)的仿真圖，包括已調(diào)信號的波形，解調(diào)后的信號波形，眼圖和誤碼率。在仿真的基礎(chǔ)上分析比較了各種調(diào)制方法的性能，并通過比較仿真模型與理論計算的性能，證明了仿真模型的可行性。關(guān)鍵字：4PSK 調(diào)制解調(diào) MATLAB 分析與仿真AbstractIn digital signal modulation method of 4PSK is the most commonly used a digital signal modulation mode, it has high frequency spectrum use efficiency, strong anti-jamming,

5、in circuit implementation are relatively simple. Modulation technique is communication field very important link, a good modulation technology can not only save spectrum resources and can provide good communication performance. 4PSK modulation is a kind of high efficiency and good resistance to band

6、 noise performance of the modulation mode, in digital mobile communication has been widely used. This design in understanding 4PSK demodulation based on the principle of MATLAB language to complete the simulation, the simulation out 4PSK modulation and demodulation of the simulation diagram, includi

7、ng the signal waveform has, after the demodulation signal waveform, the eye chart and the bit error rate. On the basis of simulation analyses and compares the performance of various modulation method, and by comparing the simulation model and theoretical calculation of performance, proved the feasib

8、ility of the simulation model.Key word: 4PSK demodulation MATLAB analysis and simulation1 引言1.1 背景介紹隨著現(xiàn)代電子技術(shù)的發(fā)展，通信技術(shù)正向著數(shù)字化、網(wǎng)絡(luò)化、智能化和寬帶化的方向發(fā)展。隨著科學(xué)技術(shù)的進(jìn)步，人們對通信的要求越來越高，各種技術(shù)會不斷地應(yīng)用于通信領(lǐng)域，各種新的通信業(yè)務(wù)將不斷地被開發(fā)出來。到那時人們的生活將越來越離不開通信。本文中提到的調(diào)制方式大都是可以實用的，已經(jīng)采用多年，并且至今仍然被采用著。但是，這些調(diào)制方法還不是很完善，有許多值得改進(jìn)之處。因此，在這些基本的數(shù)字調(diào)制方法基礎(chǔ)上，多年來

9、不斷研究出新的或改進(jìn)的調(diào)制方法。實際上，在基本的和先進(jìn)的調(diào)制方法之間并沒有明確的界限。這些方法都是不間斷地發(fā)展出來的，后來者自然比原有者更先進(jìn)。此外，隨著技術(shù)的進(jìn)步，特別是超大規(guī)模集成電路和數(shù)字信號處理技術(shù)的發(fā)展，使得復(fù)雜的電路設(shè)計得以用少量的幾塊即成電路模塊實現(xiàn)，有些硬件電路的功能還可以用軟件代替實現(xiàn)。因此使得一些較復(fù)雜的調(diào)制技術(shù)能夠容易地實現(xiàn)并投入使用。這方面的條件使得新的更復(fù)雜的調(diào)制體制迅速地不斷涌現(xiàn)。 目前，改進(jìn)的數(shù)字調(diào)制方式主要有偏置正交相移鍵控， p/4正交差分相移鍵控，最小頻移鍵控，高斯最小頻移鍵控，正交頻分復(fù)用，網(wǎng)格編碼調(diào)制等，這里對最小頻移鍵控作一介紹。1.2 設(shè)計要求1)

10、4PSK信號波形的載頻和相位參數(shù)應(yīng)隨機置或者可有幾組參數(shù)組合供選擇2) 系統(tǒng)中要求加入高斯白噪聲3) 4PSK解調(diào)方框圖采用相干接收形式4) 分析誤碼率2 4PSK調(diào)制解調(diào)的基本原理2.1 2PSK數(shù)字調(diào)制原理2PSK信號用載波相位的變化來表征被傳輸信息的狀態(tài)，通常規(guī)定0相位載波和相位載波分別表示傳“1”和傳“0”。2PSK碼元序列的波形與載頻和碼元持續(xù)時間之間的關(guān)系有關(guān)。當(dāng)一個碼元中包含有整數(shù)個載波周期時，在相鄰碼元的邊界處波形是不連續(xù)的，或者說相位是不連續(xù)的。當(dāng)一個碼元中包含的載波周期數(shù)比整數(shù)個周期多半個周期時，則相位連續(xù)。當(dāng)載波的初始相位差90度時，即余弦波改為正弦波時，結(jié)果類似。以上說

11、明，相鄰碼元的相位是否連續(xù)與相鄰碼元的初始相位是否相同不可混為一談。只有當(dāng)一個碼元中包含有整數(shù)個載波周期時，相鄰碼元邊界處的相位跳變才是由調(diào)制引起的相位變化。2PSK信號的產(chǎn)生方法主要有兩種。第一種叫相乘法，是用二進(jìn)制基帶不歸零矩形脈沖信號與載波相乘，得到相位反相的兩種碼元。第二種方法叫選擇法，是用此基帶信號控制一個開關(guān)電路，以選擇輸入信號，開關(guān)電路的輸入信號是相位相差的同頻載波。這兩種方法的復(fù)雜程度差不多，并且都可以用數(shù)字信號處理器實現(xiàn)如圖1所示。 碼變換相乘S（t）載波eo（t）雙極性不歸零圖1 2PSK及2DPSK的調(diào)制方框2.2 4PSK的調(diào)制和解調(diào)四進(jìn)制絕對相移鍵控(4PSK)直接利

12、用載波的四種不同相位來表示數(shù)字信息。如圖2所示參考相位000o11 180o01 270o10 90o45o 11135o 0100225o10315o參考相位圖2 4PSK信號相位n矢量圖由于每一種相位代表兩個比特信息，因此每個四進(jìn)制碼元可以用兩個二進(jìn)制碼元的組合來表示。兩個二進(jìn)制碼元中的前一比特用a來表示，后一比特用b表示，則雙比特ab與載波相位的關(guān)系入下表1：表1 雙比特ab與載波相位的關(guān)系雙比特碼元載波相位（n） abA方式B方式000o225o1090o315 o11180o45 o01270o135 o四進(jìn)制信號可等效為兩個正交載波進(jìn)行雙邊帶調(diào)制所得信號之和。這樣，就把數(shù)字調(diào)相和線

13、性調(diào)制聯(lián)系起來，為四相波形的產(chǎn)生提供依據(jù)。4PSK信號調(diào)制和解調(diào)（1）4PSK調(diào)制原理： 4PSK的調(diào)制方法有正交調(diào)制方式（雙路二相調(diào)制合成法或直接調(diào)相法）、相位選擇法、插入脈沖法等。這里我們采用正交調(diào)制方式。4PSK的正交調(diào)制原理如圖3：串/并變換單/雙極性換單/雙極性換×移相/2載波震蕩+acosctsinct-+輸入4PSK輸出b×圖3 4PSK正交調(diào)制原理方框圖它可以看成是由兩個載波正交的2PSK調(diào)制器構(gòu)成的。圖中串/并變換器將輸入的二進(jìn)制序列分為速度減半的兩個并行雙極性序列a和b（a,b碼元在事件上是對齊的），再分別進(jìn)行極性變換，把極性碼變?yōu)殡p極性碼（0-1，1+

14、1）然后分別調(diào)制到cosct和sinct兩個載波上，兩路相乘器輸出的信號是相互正交的抑制載波的雙邊帶調(diào)制（DSB）信號，其相位與各路碼元的極性有關(guān)，分別由a和b碼元決定。經(jīng)相加電路后輸出兩路的合成波形，即是4PSK信號。圖中兩個乘法器，其中一個用于產(chǎn)生0o與180o兩種相位狀態(tài)，另一個用于產(chǎn)生90o與270o兩種相位狀態(tài)，相加后就可以得到45o，135o，225o，和315o四種相位（2）4PSK解調(diào)原理4PSK信號是兩個載波正交的2PSK信號的合成。所以，可以仿照2PSK相干檢測法，用兩個正交的相干載波分別檢測兩個分量 a和b，然后還原成二進(jìn)制雙比特串行數(shù)字信號。此法稱作極性比較法（相干解調(diào)

15、加碼反變換器方式或相干正交解調(diào)發(fā)）如圖4帶通濾波器××低通濾波器抽樣判決抽樣判決并/串變換正交載波源4PSK輸入yi(t)yB(t)cosctsinctyA(t)zB(t)xA(t)zA(t)xB(t)ab位定時低通濾波器圖4 4PSK信號解調(diào)器原理方圖在不考慮噪聲及傳輸畸變時，接收機輸入的4PSK信號碼元可表示為 （2.2.1）表2 抽樣判決器的判決準(zhǔn)則輸入相位的極性的極性判決器輸出+11-+01-00+-10判決器是按極性來判決的。即正抽樣值判為1，負(fù)抽樣值判為0.兩路抽樣判決器輸出a、b，經(jīng)并/串變換器就可將并行數(shù)據(jù)恢復(fù)成串行數(shù)據(jù)如表2所示。3 4PSK調(diào)制解調(diào)系統(tǒng)仿

16、真3.1 MATLAB軟件介紹MATLAB 軟件是美國 Math works 公司的產(chǎn)品，MATLAB 是英文 MATrix LABoratory(矩陣實驗室)的縮寫。MATLAB軟件系列產(chǎn)品是一套高效強大的工程技術(shù)數(shù)值運算和系統(tǒng)仿真軟件，廣泛應(yīng)用于當(dāng)今的航空航天、汽車制造、半導(dǎo)體制造、電子通信、醫(yī)學(xué)研究、財經(jīng)研究和高等教育等領(lǐng)域，被譽為“巨人肩膀上的工具”。研發(fā)人員借助MATLAB軟件能迅速測試設(shè)想構(gòu)想，綜合評測系統(tǒng)性能，快速設(shè)計更好方案來確保更高技術(shù)要求。同時MATLAB也是國家教委重點提倡的一種計算工具。MATLAB主要由C語言編寫而成，采用LAPACK 為底層支持軟件包。MATLAB的

17、編程非常簡單，它有著比其他任何計算機高級語言更高的編程效率、更好的代碼可讀性和移植性，以致被譽為“第四代”計算機語言，MATLAB是所有MathWorks公司產(chǎn)品的數(shù)值分析和圖形基礎(chǔ)環(huán)境。此外MATLAB 還擁有強大的2D和3D甚至動態(tài)圖形的繪制功能，這樣用戶可以更直觀、更迅速的進(jìn)行多種算法的比較，從中找出最好的方案。從通信系統(tǒng)分析與設(shè)計、濾波器設(shè)計、信號處理、小波分析、神經(jīng)網(wǎng)絡(luò)到控制系統(tǒng)、模糊控制等方面來看，MATLAB提供了大量的面向?qū)I(yè)領(lǐng)域的工具箱。通過工具箱，以往需要復(fù)雜編程的算法開發(fā)任務(wù)往往只需一個函數(shù)就能實現(xiàn)，而且工具箱是開放的可擴展集，用戶可以查看或修改其中的算法，甚至開發(fā)自己的

18、算法。目前， MATLAB已經(jīng)廣泛地應(yīng)用于工程設(shè)計的各個領(lǐng)域，如電子、通信等領(lǐng)域；它已成為國際上最流行的計算機仿真軟件設(shè)計工具?，F(xiàn)在的MATLAB不再僅僅是一個矩陣實驗室，而是一種實用的、功能強大的、不斷更新的高級計算機編程語言?，F(xiàn)在從電子通信、自動控制圖形分析處理到航天工業(yè)、汽車工業(yè)，甚至是財務(wù)工程。MATLAB都憑借其強大的功能獲得了極大的用武之地。廣大學(xué)生可以使用MATLAB來幫助進(jìn)行信號處理、通信原理、線性系統(tǒng)、自動控制等課程的學(xué)習(xí)；科研工作者可以使用MATLAB進(jìn)行理論研究和算法開發(fā)；工程師可以使用MATLAB進(jìn)行系統(tǒng)級的設(shè)計與仿真。3.2 2PSK調(diào)制解調(diào)系統(tǒng)仿真系統(tǒng)方框圖如圖5所

19、示圖5 系統(tǒng)方框圖運行結(jié)果如下圖6所示：圖6 運行結(jié)果3.3 4PSK調(diào)制解調(diào)系統(tǒng)仿真系統(tǒng)框圖如圖7所示圖7系統(tǒng)框圖運行結(jié)果如下圖8所示：圖8 運行結(jié)果效果圖如圖9所示圖9 解調(diào)效果圖解調(diào)序列圖如圖10所示圖10 解調(diào)序列圖加入高斯白噪聲后效果圖如圖11所示圖11 假如噪聲效果圖4 4PSK誤碼率分析4.1 4PSK誤碼率的計算 對于4PSK其有以下圖12表示兩種形式參考相位000o11 180o01 270o10 90o45o 11135o 0100225o10315o參考相位圖12 4PSK信號相位矢量圖對于第二種4PSK的表達(dá)式為： （4.1.1）即： （4.1.2）已知PSK信號數(shù)字表

20、達(dá)式為： （4.1.3）計算可得2PSK誤碼率為： （4.1.4）故4PSK系統(tǒng)誤碼率為： （4.1.5）4.2 4PSK誤碼率的仿真4PSK誤碼率仿真圖如圖13所示圖13 誤碼率仿真圖5 總結(jié)在通信和信息傳輸系統(tǒng)、工業(yè)自動化或電子工程技術(shù)中，調(diào)制和解調(diào)應(yīng)用最為廣泛。本設(shè)計研究了2PSK和4PSK的調(diào)制和解調(diào)原理，以及利用MATLAB對其調(diào)制和解調(diào)進(jìn)行了編程和編譯仿真，得到的結(jié)論和理論上是一致的。簡單而且快捷。同時利用MATLAB中的SIMYULINK對2PSK和4PSK的通信系統(tǒng)進(jìn)行了仿真研究了其傳輸?shù)奶匦?，及傳輸中噪聲對系統(tǒng)的影響。而調(diào)制和解調(diào)的基本原理是利用信號與系統(tǒng)的頻域分析和傅里葉變

21、換的基本性質(zhì)，將信號的頻譜進(jìn)行搬移，使之滿足一定需要，從而完成信號的傳輸或處理。調(diào)制與解調(diào)又分模擬和數(shù)字兩種，在現(xiàn)代通信中，調(diào)制器的載波信號幾乎都是正弦信號，數(shù)字基帶信號通過調(diào)制器改變正弦載波信號的幅度、頻率或相位，產(chǎn)生幅度鍵控（ASK）、相位鍵控（PSK）、頻率鍵控（FSK）信號，或同時改變正弦載波信號的幾個參數(shù)，產(chǎn)生復(fù)合調(diào)制信號。本課程設(shè)計主要介紹基于Matlab對2PSK 和4PSK進(jìn)制的調(diào)制仿真實現(xiàn).本研究具有可對比性，對比2PSK和4PSK的通信原理和星座圖可發(fā)現(xiàn)其中的不同點，但是頻譜圖近似相同。通信中信道的信噪比設(shè)置越大信噪傳輸越理想，與理論上是相符合的。2PSK和4PSK的傳輸系

22、統(tǒng)也具有對比性，本研究在文中列出了仿真過程中每個元件的仿真參數(shù)的設(shè)置。比較其中不同點我們發(fā)現(xiàn)其中參數(shù)基本相似。也說明了他們的傳輸原理基本相同，都利用了相位的不同表示了不同的碼元傳輸參考文獻(xiàn)1 Vinay K.Ingle,數(shù)字信號處理及其MATLAB實現(xiàn).電子工業(yè)出版社，1998 2 Sanjit K.Miltra, Digital Signal Processing Laboratory Using Matlab. McGraw-Hill,20003 王紫婷.電子技術(shù)試驗教程.成都：西南交大出版社，20014 樊昌信. 通信原理（第五版）.北京：國防工業(yè)出版社，20055 孫屹，吳磊. Sim

23、ulink通信仿真開發(fā)手冊通信仿真開發(fā)手冊.北京：國防工業(yè)出版社，20046 鐘麟，王峰. MATLAB仿真技術(shù)與應(yīng)用教程.北京：國防工業(yè)出版社，2003附錄：% 調(diào)相法clear allclose allt=-1:0.01:7-0.01;tt=length(t);x1=ones(1,800);for i=1:tt if (t(i)>=-1 & t(i)<=1) | (t(i)>=5& t(i)<=7); x1(i)=1; else x1(i)=-1; endendt1=0:0.01:8-0.01;t2=0:0.01:7-0.01;t3=-1:0.01:

24、7.1-0.01;t4=0:0.01:8.1-0.01;tt1=length(t1);x2=ones(1,800);for i=1:tt1 if (t1(i)>=0 & t1(i)<=2) | (t1(i)>=4& t1(i)<=8); x2(i)=1; else x2(i)=-1; endendf=0:0.1:1;xrc=0.5+0.5*cos(pi*f);y1=conv(x1,xrc)/5.5;y2=conv(x2,xrc)/5.5;n0=randn(size(t2);f1=1;i=x1.*cos(2*pi*f1*t);q=x2.*sin(2*pi*

25、f1*t1);I=i(101:800);Q=q(1:700);QPSK=sqrt(1/2).*I+sqrt(1/2).*Q;QPSK_n=(sqrt(1/2).*I+sqrt(1/2).*Q)+n0;n1=randn(size(t2);i_rc=y1.*cos(2*pi*f1*t3);q_rc=y2.*sin(2*pi*f1*t4);I_rc=i_rc(101:800);Q_rc=q_rc(1:700);QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc);QPSK_rc_n1=QPSK_rc+n1;figure(1)subplot(4,1,1);plot(t3,

26、i_rc);axis(-1 8 -1 1);ylabel('a序列');subplot(4,1,2);plot(t4,q_rc);axis(-1 8 -1 1);ylabel('b序列');subplot(4,1,3);plot(t2,QPSK_rc);axis(-1 8 -1 1);ylabel('合成序列');subplot(4,1,4);plot(t2,QPSK_rc_n1);axis(-1 8 -1 1);ylabel('加入噪聲');% 設(shè)定 T=1,加入高斯噪聲clear allclose all% 調(diào)制bit_in

27、= randint(1e3, 1, 0 1);bit_I = bit_in(1:2:1e3);bit_Q = bit_in(2:2:1e3);data_I = -2*bit_I+1;data_Q = -2*bit_Q+1;data_I1=repmat(data_I',20,1);data_Q1=repmat(data_Q',20,1);for i=1:1e4 data_I2(i)=data_I1(i); data_Q2(i)=data_Q1(i);end;f=0:0.1:1;xrc=0.5+0.5*cos(pi*f);data_I2_rc=conv(data_I2,xrc)/5

28、.5;data_Q2_rc=conv(data_Q2,xrc)/5.5;f1=1;t1=0:0.1:1e3+0.9;n0=rand(size(t1);I_rc=data_I2_rc.*cos(2*pi*f1*t1);Q_rc=data_Q2_rc.*sin(2*pi*f1*t1);QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc);QPSK_rc_n0=QPSK_rc+n0;% 解調(diào)I_demo=QPSK_rc_n0.*cos(2*pi*f1*t1);Q_demo=QPSK_rc_n0.*sin(2*pi*f1*t1);% 低通濾波I_recover=conv(

29、I_demo,xrc); Q_recover=conv(Q_demo,xrc);I=I_recover(11:10010);Q=Q_recover(11:10010);t2=0:0.05:1e3-0.05;t3=0:0.1:1e3-0.1;% 抽樣判決data_recover=;for i=1:20:10000 data_recover=data_recover I(i:1:i+19) Q(i:1:i+19);end;bit_recover=;for i=1:20:20000 if sum(data_recover(i:i+19)>0 data_recover_a(i:i+19)=1;

30、bit_recover=bit_recover 1; else data_recover_a(i:i+19)=-1; bit_recover=bit_recover -1; endenderror=0;dd = -2*bit_in+1;ddd=dd'ddd1=repmat(ddd,20,1);for i=1:2e4 ddd2(i)=ddd1(i);endfor i=1:1e3 if bit_recover(i)=ddd(i) error=error+1; endendp=error/1000;figure(1)subplot(2,1,1);plot(t2,ddd2);axis(0 10

31、0 -2 2);title('原序列');subplot(2,1,2);plot(t2,data_recover_a);axis(0 100 -2 2);title('解調(diào)后序列');% 設(shè)定 T=1, 不加噪聲clear allclose all% 調(diào)制bit_in = randint(1e3, 1, 0 1);bit_I = bit_in(1:2:1e3);bit_Q = bit_in(2:2:1e3);data_I = -2*bit_I+1;data_Q = -2*bit_Q+1;data_I1=repmat(data_I',20,1);data_

32、Q1=repmat(data_Q',20,1);for i=1:1e4 data_I2(i)=data_I1(i); data_Q2(i)=data_Q1(i);end;t=0:0.1:1e3-0.1;f=0:0.1:1;xrc=0.5+0.5*cos(pi*f);data_I2_rc=conv(data_I2,xrc)/5.5;data_Q2_rc=conv(data_Q2,xrc)/5.5;f1=1;t1=0:0.1:1e3+0.9;I_rc=data_I2_rc.*cos(2*pi*f1*t1);Q_rc=data_Q2_rc.*sin(2*pi*f1*t1);QPSK_rc=(

33、sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc);% 解調(diào)I_demo=QPSK_rc.*cos(2*pi*f1*t1);Q_demo=QPSK_rc.*sin(2*pi*f1*t1);I_recover=conv(I_demo,xrc);Q_recover=conv(Q_demo,xrc);I=I_recover(11:10010);Q=Q_recover(11:10010);t2=0:0.05:1e3-0.05;t3=0:0.1:1e3-0.1;data_recover=;for i=1:20:10000 data_recover=data_recover I(i:1:i+

34、19) Q(i:1:i+19);end;ddd = -2*bit_in+1;ddd1=repmat(ddd',10,1);for i=1:1e4 ddd2(i)=ddd1(i);endfigure(1)subplot(4,1,1);plot(t3,I);axis(0 20 -6 6);subplot(4,1,2);plot(t3,Q);axis(0 20 -6 6);subplot(4,1,3);plot(t2,data_recover);axis(0 20 -6 6);subplot(4,1,4);plot(t,ddd2);axis(0 20 -6 6);% QPSK誤碼率分析SNR

35、indB1=0:2:10;SNRindB2=0:0.1:10;for i=1:length(SNRindB1) pb,ps=cm_sm32(SNRindB1(i); smld_bit_err_prb(i)=pb; smld_symbol_err_prb(i)=ps;end;for i=1:length(SNRindB2) SNR=exp(SNRindB2(i)*log(10)/10); theo_err_prb(i)=Qfunct(sqrt(2*SNR);end;title('QPSK誤碼率分析');semilogy(SNRindB1,smld_bit_err_prb,'*');axis(0 10 10e-8 1);hold on;% semilogy(SNRindB1,smld_symbol_err_prb,'o');semilogy(SNRi