1、<p><b>　　課程設(shè)計說明書</b></p><p>　　題目：利用MATLAB實現(xiàn)QPSK調(diào)制及解調(diào)</p><p>　　姓 名： </p><p>　　院 系： 電氣信息工程學(xué)院 </p><p>　　專業(yè)班級： 電子信息工程09-1 </p

2、><p>　　學(xué) 號： </p><p>　　指導(dǎo)教師： </p><p>　　成 績： </p><p>　　時間： 2012 年 6 月 18 日至 2012 年 6 月 22 日</p><p&

3、gt;　　課 程 設(shè) 計 任 務(wù) 書</p><p>　　題目 利用MATLAB實現(xiàn)QPSK調(diào)制及解調(diào) </p><p>　　專業(yè)班級 電子信息工程09級 1班 學(xué)號 54 姓名 </p><p>　　主要內(nèi)容、基本要求、主要參考資料等：</p>&

4、lt;p><b>　　主要內(nèi)容：</b></p><p>　　已知數(shù)字信號1011000101101011，碼元速率為2400波特，載波頻率為1200Hz，利用MATLAB畫出QPSK調(diào)制波形，并畫出調(diào)制信號經(jīng)過高斯信道傳輸后解調(diào)波形及接收誤碼率，將其與理論值進行比較。</p><p><b>　　基本要求：</b></p>

5、<p>　　通過本課程設(shè)計，鞏固通信原理QPSK調(diào)制的有關(guān)知識；</p><p>　　熟悉QPSK產(chǎn)生原理；</p><p>　　熟悉高斯信道的建模及QPSK解調(diào)原理；</p><p>　　熟悉誤碼率的蒙特卡羅仿真；</p><p>　　學(xué)會用MATLAB來進行通信系統(tǒng)仿真。</p><p><b>

6、;　　主要參考資料：</b></p><p><b>　　主要參考資料：</b></p><p>　　1、王秉鈞等. 通信原理[M].北京：清華大學(xué)出版社，2006.11</p><p>　　2、陳懷琛.數(shù)字信號處理教程----MATLAB釋義與實現(xiàn)[M].北京：電子工業(yè)出版社,2004.</p><p>　

7、　完 成 期 限： 2012.6.18—2012.6.23 </p><p>　　指導(dǎo)教師簽名： </p><p>　　課程負責(zé)人簽名： </p><p>　　2012年 6月 16日</p><p><b>　　目 錄</b><

8、;/p><p><b>　　一 前言4</b></p><p>　　1.1QPSK系統(tǒng)的應(yīng)用背景簡介4</p><p>　　1.2 QPSK實驗仿真的意義4</p><p>　　1.3 實驗平臺和實驗內(nèi)容4</p><p>　　1.3.1實驗平臺4</p><p>　

9、　1.3.2實驗內(nèi)容4</p><p>　　二、系統(tǒng)實現(xiàn)框圖和分析5</p><p>　　2.1、QPSK調(diào)制部分，5</p><p>　　2.2、QPSK解調(diào)部分6</p><p>　　三、實驗結(jié)果及分析7</p><p>　　3.1、理想信道下的仿真7</p><p>　　3.2

10、、高斯信道下的仿真8</p><p>　　3.3、先通過瑞利衰落信道再通過高斯信道的仿真9</p><p><b>　　參考文獻：11</b></p><p><b>　　附錄12</b></p><p>　　基于MATLAB的QPSK仿真設(shè)計與實現(xiàn)</p><p>

11、<b>　　一 前言</b></p><p>　　1.1QPSK系統(tǒng)的應(yīng)用背景簡介</p><p>　　QPSK是英文Quadrature Phase Shift Keying的縮略語簡稱，意為正交相移鍵控，是一種數(shù)字調(diào)制方式。在19世紀(jì)80年代初期,人們選用恒定包絡(luò)數(shù)字調(diào)制。這類數(shù)字調(diào)制技術(shù)的優(yōu)點是已調(diào)信號具有相對窄的功率譜和對放大設(shè)備沒有線性要求,不足之處是其頻譜

12、利用率低于線性調(diào)制技術(shù)。19世紀(jì)80年代中期以后,四相絕對移相鍵控(QPSK)技術(shù)以其抗干擾性能強、誤碼性能好、頻譜利用率高等優(yōu)點,廣泛應(yīng)用于數(shù)字微波通信系統(tǒng)、數(shù)字衛(wèi)星通信系統(tǒng)、寬帶接入、移動通信及有線電視系統(tǒng)之中。</p><p>　　1.2 QPSK實驗仿真的意義</p><p>　　通過完成設(shè)計內(nèi)容， 復(fù)習(xí)QPSK調(diào)制解調(diào)的基本原理，同時也要復(fù)習(xí)通信系統(tǒng)的主要組成部分，了解調(diào)制解調(diào)方

13、式中最基礎(chǔ)的方法。了解QPSK的實現(xiàn)方法及數(shù)學(xué)原理。并對“通信”這個概念有個整體的理解，學(xué)習(xí)數(shù)字調(diào)制中誤碼率測試的標(biāo)準(zhǔn)及計算方法。同時還要復(fù)習(xí)隨機信號中時域用自相關(guān)函數(shù)，頻域用功率譜密度來描述平穩(wěn)隨機過程的特性等基礎(chǔ)知識，來理解高斯信道中噪聲的表示方法，以便在編程中使用。 </p><p>　　理解QPSK調(diào)制解調(diào)的基本原理，并使用MATLAB編程實現(xiàn)QPSK信號在高斯信道和瑞利衰落信道下傳輸，以及該方式的誤碼率

14、測試。復(fù)習(xí)MATLAB編程的基礎(chǔ)知識和編程的常用算法以及使用MATLAB仿真系統(tǒng)的注意事項，并鍛煉自己的編程能力，通過編程完成QPSK調(diào)制解調(diào)系統(tǒng)的仿真，以及誤碼率測試，并得出響應(yīng)波形。在完成要求任務(wù)的條件下，嘗試優(yōu)化程序。</p><p>　　通過本次實驗，除了和隊友培養(yǎng)了默契學(xué)到了知識之外，還可以將次實驗作為一種推廣，讓更多的學(xué)生來深入一層的了解QPSK以至其他調(diào)制方式的原理和實現(xiàn)方法?？梢苑奖銓W(xué)生進行測試和

15、對比。足不出戶便可以做實驗。</p><p>　　1.3 實驗平臺和實驗內(nèi)容</p><p><b>　　1.3.1實驗平臺</b></p><p>　　本實驗是基于Matlab的軟件仿真，只需PC機上安裝MATLAB 6.0或者以上版本即可。</p><p>　?。ū緦嶒灨綆Щ贛atlab Simulink （模塊化

16、）仿真，如需使用必須安裝simulink 模塊）</p><p><b>　　1.3.2實驗內(nèi)容</b></p><p>　　1.構(gòu)建一個理想信道基本QPSK仿真系統(tǒng),要求仿真結(jié)果有</p><p>　　a.基帶輸入波形及其功率譜 </p><p>　　b.QPSK信號及其功率譜 </p>&l

17、t;p>　　c.QPSK信號星座圖 </p><p>　　2.構(gòu)建一個在AWGN（高斯白噪聲）信道條件下的QPSK仿真系統(tǒng)，要求仿真結(jié)果有</p><p>　　a.QPSK信號及其功率譜 </p><p>　　b.QPSK信號星座圖</p><p>　　c.高斯白噪聲信道條件下的誤碼性能以及高斯白噪聲的理論曲線，

18、要求所有誤碼性能曲線在同一坐標(biāo)比例下繪制</p><p>　　3驗可選做擴展內(nèi)容要求：</p><p>　　構(gòu)建一個先經(jīng)過Rayleigh（瑞利衰落信道），再通過AWGN（高斯白噪聲）信道條件下的條件下的QPSK仿真系統(tǒng)，要求仿真結(jié)果有</p><p>　　a.QPSK信號及其功率譜 </p><p>　　b.通過瑞利衰落信道之前和之后

19、的信號星座圖，前后進行比較</p><p>　　c.在瑞利衰落信道和在高斯白噪聲條件下的誤碼性能曲線，并和二.2.c中所要求的誤碼性能曲線在同一坐標(biāo)比例下繪制</p><p>　　二、系統(tǒng)實現(xiàn)框圖和分析</p><p>　　2.1、QPSK調(diào)制部分，</p><p><b>　　原理框圖如圖1所示</b></p&g

20、t;<p><b>　　1（t）＝</b></p><p><b>　　2（t）＝ </b></p><p><b>　　圖1</b></p><p><b>　　原理分析：</b></p><p><b>　　基本原理及系統(tǒng)結(jié)構(gòu)&

21、lt;/b></p><p>　　QPSK與二進制PSK一樣，傳輸信號包含的信息都存在于相位中。的別的載波相位取四個等間隔值之一，如л/4, 3л/4,5л/4,和7л/4。相應(yīng)的，可將發(fā)射信號定義為</p><p><b>　　0≤t≤T</b></p><p><b>　　Si（t） ＝</b></p&g

22、t;<p>　　0。， 其他</p><p>　　其中，i＝1，2，2，4；E為發(fā)射信號的每個符號的能量，T為符號持續(xù)時間，載波頻率f等于nc/T，nc為固定整數(shù)。每一個可能的相位值對應(yīng)于一個特定的二位組。例如，可用前述的一組相位值來表示格雷碼的一組二位組：10，00，01，11。</p><p>　　下面介紹QPSK信號的產(chǎn)生和檢測。如果a為典型的

23、QPSK發(fā)射機框圖。輸入的二進制數(shù)據(jù)序列首先被不歸零（NRZ）電平編碼轉(zhuǎn)換器轉(zhuǎn)換為極性形式，即負號1和0分別用和－表示。接著，該二進制波形被分接器分成兩個分別由輸入序列的奇數(shù)位偶數(shù)位組成的彼此獨立的二進制波形，這兩個二進制波形分別用a1（t），和a2（t）表示。容易注意到，在任何一信號時間間隔內(nèi)a1（t），和a2（t）的幅度恰好分別等于Si1和 Si2，即由發(fā)送的二位組決定。這兩個二進制波形a1（t），和a2（t）被用來調(diào)制一對正交載波

24、或者說正交基本函數(shù)：1（t）＝，2（t）＝。這樣就得到一對二進制PSK信號。1（t）和2（t）的正交性使這兩個信號可以被獨立地檢測。最后，將這兩個二進制PSK信號相加，從而得期望的QPSK。</p><p>　　2.2、QPSK解調(diào)部分</p><p>　　，原理框圖如圖2所示：</p><p>　　1（t） 同相信道 門限＝0</p>

25、;<p>　　2（t） 正交信道 門限＝0 </p><p><b>　　圖2</b></p><p><b>　　原理分析：</b></p><p>　　QPSK接收機由一對共輸入地相關(guān)器組成。這兩個相關(guān)器分別提供本地產(chǎn)生地相干參考信號1（t）和2（t）。相關(guān)器接收信號x（t），相關(guān)器輸出地x

26、1和x2被用來與門限值0進行比較。如果x1>0,則判決同相信道地輸出為符號1；如果x1<0 ,則判決同相信道的輸出為符號0。；類似地。如果正交通道也是如此判決輸出。最后同相信道和正交信道輸出這兩個二進制數(shù)據(jù)序列被復(fù)加器合并，重新得到原始的二進制序列。在AWGN信道中，判決結(jié)果具有最小的負號差錯概率。</p><p><b>　　三、實驗結(jié)果及分析</b></p>&

27、lt;p>　　根據(jù)圖1和圖2的流程框圖設(shè)計仿真程序，得出結(jié)果并且分析如下：</p><p>　　3.1、理想信道下的仿真，實驗結(jié)果如圖3所示</p><p><b>　　圖3</b></p><p><b>　　實驗結(jié)果分析：</b></p><p>　　如圖上結(jié)果顯示，完成了QPSK信號在理

28、想信道上的調(diào)制，傳輸，解調(diào)的過程，由于調(diào)制過程中加進了載波，因此調(diào)制信號的功率譜密度會發(fā)生變化。并且可以看出調(diào)制解調(diào)的結(jié)果沒有誤碼。</p><p>　　3.2、高斯信道下的仿真，結(jié)果如圖4所示：</p><p><b>　　圖4</b></p><p><b>　　實驗結(jié)果分析：</b></p><p

29、>　　由圖4可以得到高斯信道下的調(diào)制信號，高斯噪聲，調(diào)制輸出功率譜密度曲線和QPSK信號的星座圖。</p><p>　　在高斯噪聲的影響下，調(diào)制信號的波形發(fā)生了明顯的變化，其功率譜密度函數(shù)相對于圖1中的調(diào)制信號的功率譜密度只發(fā)生了微小的變化，原因在于高斯噪聲是一個均值為0的白噪聲，在各個頻率上其功率是均勻的，因此此結(jié)果是真確的。星座圖反映可接收信號早高斯噪聲的影響下發(fā)生了誤碼，但是大部分還是保持了原來的特性

30、。</p><p>　　3.3、先通過瑞利衰落信道再通過高斯信道的仿真。實驗結(jié)果如圖5所示：</p><p><b>　　圖5</b></p><p><b>　　實驗結(jié)果分析：</b></p><p>　　由圖5可以得到瑞利衰落信道前后的星座圖，調(diào)制信號的曲線圖及其功率譜密度。最后顯示的是高斯信道

31、和瑞利衰落信道的誤碼率對比。由圖可知瑞利衰落信道下的誤碼率比高斯信道下的誤碼率高。</p><p>　　至此，仿真實驗就全部完成。</p><p><b>　　參考文獻：</b></p><p>　　1、《MATLAB 寶典》 陳杰等編著 電子工業(yè)出版社</p><p>　　2、《MATLAB信號處理》 劉波, 文忠,

32、 曾涯編著 北京電子工業(yè)出版社</p><p>　　3、《數(shù)字信號處理的MATLAB實現(xiàn)》 萬永革編著 北京科學(xué)出版社</p><p><b>　　4、網(wǎng)上資料</b></p><p><b>　　附錄</b></p><p><b>　　MATLAB程序</b></p&

33、gt;<p><b>　　% 調(diào)相法</b></p><p><b>　　clear all</b></p><p><b>　　close all</b></p><p>　　t=[-1:0.01:7-0.01];</p><p>　　tt=length(t);

34、</p><p>　　x1=ones(1,800);</p><p>　　for i=1:tt</p><p>　　if (t(i)>=-1 & t(i)<=1) | (t(i)>=5& t(i)<=7);</p><p><b>　　x1(i)=1;</b></p>

35、<p>　　else x1(i)=-1;</p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　t1=[0:0.01:8-0.01];</p><p>　　t2=0:0.01:7-0.01;</p><p>　

36、　t3=-1:0.01:7.1-0.01;</p><p>　　t4=0:0.01:8.1-0.01;</p><p>　　tt1=length(t1);</p><p>　　x2=ones(1,800);</p><p>　　for i=1:tt1</p><p>　　if (t1(i)>=0 & t1

37、(i)<=2) | (t1(i)>=4& t1(i)<=8);</p><p><b>　　x2(i)=1;</b></p><p>　　else x2(i)=-1;</p><p><b>　　end</b></p><p><b>　　end</b>

38、;</p><p>　　f=0:0.1:1;</p><p>　　xrc=0.5+0.5*cos(pi*f);</p><p>　　y1=conv(x1,xrc)/5.5;</p><p>　　y2=conv(x2,xrc)/5.5;</p><p>　　n0=randn(size(t2));</p>&

39、lt;p><b>　　f1=1;</b></p><p>　　i=x1.*cos(2*pi*f1*t);</p><p>　　q=x2.*sin(2*pi*f1*t1);</p><p>　　I=i(101:800);</p><p>　　Q=q(1:700);</p><p>　　QPSK

40、=sqrt(1/2).*I+sqrt(1/2).*Q;</p><p>　　QPSK_n=(sqrt(1/2).*I+sqrt(1/2).*Q)+n0;</p><p>　　n1=randn(size(t2));</p><p>　　i_rc=y1.*cos(2*pi*f1*t3);</p><p>　　q_rc=y2.*sin(2*pi*f

41、1*t4);</p><p>　　I_rc=i_rc(101:800);</p><p>　　Q_rc=q_rc(1:700);</p><p>　　QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc);</p><p>　　QPSK_rc_n1=QPSK_rc+n1;</p><p>&

42、lt;b>　　figure(1)</b></p><p>　　subplot(4,1,1);plot(t3,i_rc);axis([-1 8 -1 1]);ylabel('a序列');</p><p>　　subplot(4,1,2);plot(t4,q_rc);axis([-1 8 -1 1]);ylabel('b序列');</p&

43、gt;<p>　　subplot(4,1,3);plot(t2,QPSK_rc);axis([-1 8 -1 1]);ylabel('合成序列');</p><p>　　subplot(4,1,4);plot(t2,QPSK_rc_n1);axis([-1 8 -1 1]);ylabel('加入噪聲');</p><p><b>　　

44、效果圖：</b></p><p>　　% 設(shè)定 T=1,加入高斯噪聲</p><p><b>　　clear all</b></p><p><b>　　close all</b></p><p><b>　　% 調(diào)制</b></p><p>

45、;　　bit_in = randint(1e3, 1, [0 1]);</p><p>　　bit_I = bit_in(1:2:1e3);</p><p>　　bit_Q = bit_in(2:2:1e3);</p><p>　　data_I = -2*bit_I+1;</p><p>　　data_Q = -2*bit_Q+1;<

46、;/p><p>　　data_I1=repmat(data_I',20,1);</p><p>　　data_Q1=repmat(data_Q',20,1);</p><p>　　for i=1:1e4</p><p>　　data_I2(i)=data_I1(i);</p><p>　　data_Q2(i

47、)=data_Q1(i);</p><p><b>　　end;</b></p><p>　　f=0:0.1:1;</p><p>　　xrc=0.5+0.5*cos(pi*f);</p><p>　　data_I2_rc=conv(data_I2,xrc)/5.5;</p><p>　　data

48、_Q2_rc=conv(data_Q2,xrc)/5.5;</p><p><b>　　f1=1;</b></p><p>　　t1=0:0.1:1e3+0.9;</p><p>　　n0=rand(size(t1));</p><p>　　I_rc=data_I2_rc.*cos(2*pi*f1*t1);</p&

49、gt;<p>　　Q_rc=data_Q2_rc.*sin(2*pi*f1*t1);</p><p>　　QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc);</p><p>　　QPSK_rc_n0=QPSK_rc+n0;</p><p><b>　　% 解調(diào)</b></p><

50、;p>　　I_demo=QPSK_rc_n0.*cos(2*pi*f1*t1);</p><p>　　Q_demo=QPSK_rc_n0.*sin(2*pi*f1*t1);</p><p><b>　　% 低通濾波</b></p><p>　　I_recover=conv(I_demo,xrc); </p><p

51、>　　Q_recover=conv(Q_demo,xrc);</p><p>　　I=I_recover(11:10010);</p><p>　　Q=Q_recover(11:10010);</p><p>　　t2=0:0.05:1e3-0.05;</p><p>　　t3=0:0.1:1e3-0.1;</p><

52、;p><b>　　% 抽樣判決</b></p><p>　　data_recover=[];</p><p>　　for i=1:20:10000</p><p>　　data_recover=[data_recover I(i:1:i+19) Q(i:1:i+19)];</p><p><b>　　en

53、d;</b></p><p>　　bit_recover=[];</p><p>　　for i=1:20:20000</p><p>　　if sum(data_recover(i:i+19))>0</p><p>　　data_recover_a(i:i+19)=1;</p><p>　　bit_

54、recover=[bit_recover 1];</p><p><b>　　else</b></p><p>　　data_recover_a(i:i+19)=-1;</p><p>　　bit_recover=[bit_recover -1];</p><p><b>　　end</b><

55、/p><p><b>　　end</b></p><p><b>　　error=0;</b></p><p>　　dd = -2*bit_in+1;</p><p>　　ddd=[dd'];</p><p>　　ddd1=repmat(ddd,20,1);</p&

56、gt;<p>　　for i=1:2e4</p><p>　　ddd2(i)=ddd1(i);</p><p><b>　　end</b></p><p>　　for i=1:1e3</p><p>　　if bit_recover(i)~=ddd(i)</p><p>　　erro

57、r=error+1;</p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　p=error/1000;</p><p><b>　　figure(1)</b></p><p>　　subplot(2,

58、1,1);plot(t2,ddd2);axis([0 100 -2 2]);title('原序列');</p><p>　　subplot(2,1,2);plot(t2,data_recover_a);axis([0 100 -2 2]);title('解調(diào)后序列');</p><p><b>　　效果圖：</b></p>

59、<p>　　% 設(shè)定 T=1, 不加噪聲</p><p><b>　　clear all</b></p><p><b>　　close all</b></p><p><b>　　% 調(diào)制</b></p><p>　　bit_in = randint(1e3, 1

60、, [0 1]);</p><p>　　bit_I = bit_in(1:2:1e3);</p><p>　　bit_Q = bit_in(2:2:1e3);</p><p>　　data_I = -2*bit_I+1;</p><p>　　data_Q = -2*bit_Q+1;</p><p>　　data_I

61、1=repmat(data_I',20,1);</p><p>　　data_Q1=repmat(data_Q',20,1);</p><p>　　for i=1:1e4</p><p>　　data_I2(i)=data_I1(i);</p><p>　　data_Q2(i)=data_Q1(i);</p>&

62、lt;p><b>　　end;</b></p><p>　　t=0:0.1:1e3-0.1;</p><p>　　f=0:0.1:1;</p><p>　　xrc=0.5+0.5*cos(pi*f);</p><p>　　data_I2_rc=conv(data_I2,xrc)/5.5;</p>&l

63、t;p>　　data_Q2_rc=conv(data_Q2,xrc)/5.5;</p><p><b>　　f1=1;</b></p><p>　　t1=0:0.1:1e3+0.9;</p><p>　　I_rc=data_I2_rc.*cos(2*pi*f1*t1);</p><p>　　Q_rc=data_Q2

64、_rc.*sin(2*pi*f1*t1);</p><p>　　QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc);</p><p><b>　　% 解調(diào)</b></p><p>　　I_demo=QPSK_rc.*cos(2*pi*f1*t1);</p><p>　　Q_demo=QPS

65、K_rc.*sin(2*pi*f1*t1);</p><p>　　I_recover=conv(I_demo,xrc);</p><p>　　Q_recover=conv(Q_demo,xrc);</p><p>　　I=I_recover(11:10010);</p><p>　　Q=Q_recover(11:10010);</p&g

66、t;<p>　　t2=0:0.05:1e3-0.05;</p><p>　　t3=0:0.1:1e3-0.1;</p><p>　　data_recover=[];</p><p>　　for i=1:20:10000</p><p>　　data_recover=[data_recover I(i:1:i+19) Q(i:1:

67、i+19)];</p><p><b>　　end;</b></p><p>　　ddd = -2*bit_in+1;</p><p>　　ddd1=repmat(ddd',10,1);</p><p>　　for i=1:1e4</p><p>　　ddd2(i)=ddd1(i);<

68、;/p><p><b>　　end</b></p><p><b>　　figure(1)</b></p><p>　　subplot(4,1,1);plot(t3,I);axis([0 20 -6 6]);</p><p>　　subplot(4,1,2);plot(t3,Q);axis([0 20

69、-6 6]);</p><p>　　subplot(4,1,3);plot(t2,data_recover);axis([0 20 -6 6]);</p><p>　　subplot(4,1,4);plot(t,ddd2);axis([0 20 -6 6]);</p><p><b>　　效果圖：</b></p><p>

70、;　　% QPSK誤碼率分析</p><p>　　SNRindB1=0:2:10;</p><p>　　SNRindB2=0:0.1:10;</p><p>　　for i=1:length(SNRindB1)</p><p>　　[pb,ps]=cm_sm32(SNRindB1(i));</p><p>　　smld

71、_bit_err_prb(i)=pb;</p><p>　　smld_symbol_err_prb(i)=ps;</p><p><b>　　end;</b></p><p>　　for i=1:length(SNRindB2)</p><p>　　SNR=exp(SNRindB2(i)*log(10)/10);<

72、/p><p>　　theo_err_prb(i)=Qfunct(sqrt(2*SNR));</p><p><b>　　end;</b></p><p>　　title('QPSK誤碼率分析');</p><p>　　semilogy(SNRindB1,smld_bit_err_prb,'*'

73、);</p><p>　　axis([0 10 10e-8 1]);</p><p><b>　　hold on;</b></p><p>　　% semilogy(SNRindB1,smld_symbol_err_prb,'o');</p><p>　　semilogy(SNRindB2,theo_er

74、r_prb);</p><p>　　legend('仿真比特誤碼率','理論比特誤碼率');</p><p><b>　　hold off;</b></p><p>　　function[y]=Qfunct(x)</p><p>　　y=(1/2)*erfc(x/sqrt(2));<

75、/p><p>　　function[pb,ps]=cm_sm32(SNRindB)</p><p><b>　　N=10000;</b></p><p><b>　　E=1;</b></p><p>　　SNR=10^(SNRindB/10);</p><p>　　sgma=sq

76、rt(E/SNR)/2;</p><p>　　s00=[1 0];</p><p>　　s01=[0 1];</p><p>　　s11=[-1 0];</p><p>　　s10=[0 -1];</p><p><b>　　for i=1:N</b></p><p>　　

77、dsource1(i)=[1 0 1 1 0 0 0 1 0 1 1 0 1 0 1 1];</p><p>　　numofsymbolerror=0;</p><p>　　numofbiterror=0;</p><p><b>　　for i=1:N</b></p><p>　　n=sgma*randn(size(s

78、00));</p><p>　　if((dsource1(i)==0)&(dsource2(i)==0))</p><p><b>　　r=s00+n;</b></p><p>　　elseif((dsource1(i)==0)&(dsource2(i)==1))</p><p><b>　　r

79、=s01+n;</b></p><p>　　elseif((dsource1(i)==1)&(dsource2(i)==0))</p><p><b>　　r=s10+n;</b></p><p><b>　　else</b></p><p><b>　　r=s11+n

80、;</b></p><p><b>　　end;</b></p><p>　　c00=dot(r,s00);</p><p>　　c01=dot(r,s01);</p><p>　　c10=dot(r,s10);</p><p>　　c11=dot(r,s11);</p>

81、<p>　　c_max=max([c00 c01 c10 c11]);</p><p>　　if (c00==c_max)</p><p>　　decis1=0;decis2=0;</p><p>　　elseif(c01==c_max)</p><p>　　decis1=0;decis2=1;</p><p&

82、gt;　　elseif(c10==c_max)</p><p>　　decis1=1;decis2=0;</p><p><b>　　else</b></p><p>　　decis1=1;decis2=1;</p><p><b>　　end;</b></p><p>　　

83、symbolerror=0;</p><p>　　if(decis1~=dsource1(i))</p><p>　　numofbiterror=numofbiterror+1;</p><p>　　symbolerror=1;</p><p><b>　　end;</b></p><p>　　i

84、f(decis2~=dsource2(i))</p><p>　　numofbiterror=numofbiterror+1;</p><p>　　symbolerror=1;</p><p><b>　　end;</b></p><p>　　if(symbolerror==1)</p><p>

85、　　numofsymbolerror=numofsymbolerror+1;</p><p><b>　　end;</b></p><p><b>　　end;</b></p><p>　　ps=numofsymbolerror/N;</p><p>　　pb=numofbiterror/(2*N)