1、探2009級通信工程專業*通信原理課程設計*探通信原理課程設計報告書16QAM調制與解調課題名稱的MATLA實現及調制性能分析姓 名學 號學 院通信與電子工程學院專 業通信工程指導教師夢醒2012年01月01日一、設計任務及要求:設計任務：利用MATLA設計一個16QAMS制與解調系統，并對其進行性能分析。要 求：1. 設計一個16QA碉制與解調系統。2. 設計程序時必須使得程序盡可能的簡單。3. 利用MATLA進行程序編寫并對系統進行仿真分析。指導教師簽名：20 年 月曰、指導教師評語:指導教師簽名：2010 年 月 日二、成績驗收蓋章2010年 月 日16QA碉制與解調的MATLA實現及調

2、制性能分析1設計目的(1) 掌握16QAM調制與解調的原理。(2) 掌握星座圖的原理并能熟悉星座圖的應用。(3) 熟悉并掌握MATLAB勺使用方法。(4) 通過對16QAM調制性能的分析了解16QA調制相對于其它調制方式的優缺 點。2設計原理正交振幅調制(Quadrature Amplitude Modulation,QAM )是一種振幅和相位聯合鍵控。雖然MPSI和MDPS等相移鍵控的帶寬和功率方面都具有優勢，即 帶寬占用小和比特噪聲比要求低。但是由圖1可見，在MPS處制中，隨著圖1 8PSK信號相位M的增大，相鄰相位的距離逐漸減小，使噪聲容限隨之減小，誤碼率難于保證。為了改善在M大時的噪聲

3、容限，發展出了 QAM體制。在QAM體制中，信號的振幅和相位作為兩個獨立的參量同時受到調制。這種信號的一個碼元可以表示為Sk(t)Ak cos( otk)kT t(k1)T(2 1)式中：k=整數；Ak和k分別可以取多個離散值。式(2 1)可以展開為沁)Ak cos k cos 0tAkSin ksin°t(22)令Xk=Akcos k，Yk=-A kSink則式(2 1)變為Sk(t)XkCos 0tYk sin0t(2 3)Xk和Yk也是可以取多個離散的變量。從式(2 3)看出，Sk(t)可以看作是兩個正交的振幅鍵控信號之和。在式（2 1）中，若 k值僅可以取 /4和-/4 ,

4、Ak值僅可以取+A和-A,則此QAM言號就成為QPSK言號，如圖2所示：所以，QPSI信號就是一種最簡單的QAM信號。有代表性的QAM言號是16進制的, 記為16QA M它的矢量圖示于下圖中：圖3 16QAM信號矢量圖圖中用黑點表示每個碼元的位置，并且示出它是由兩個正交矢量合成的。類似地,有64QAM和 256QAM等 QAM言號，如圖4、圖5所示。它們總稱為 MQAM圖4 64 QAM信號矢量圖圖5 256QAM信號矢量圖調制。由于從其矢量圖看像是星座，故又稱星座調制。16QAMW號的產生方法主要有兩種。第一種是正交調幅法，即用兩路獨立的 正交4ASK信號疊加，形成16QAMW號，如圖6所示

5、。第二種方法是復合相圖6正交調幅法移法，它用兩路獨立的圖7復合相移法虛線大圓上的4個大黑點表示一個QPSK言號矢量的位置。在這4個位置上可以 疊加上第二個QPSK矢量，后者的位置用虛線小圓上的4個小黑點表示。3設計過程3.1設計思路由設計原理中可知 MQA調制又稱為星座調制，故我們在設計 16QAMM制系 統時就可以星座圖來進行編程。下面我們就借用如圖8所示的星座圖設計一個16QAM調制系統。圖8 16QAM星座在圖中共有16個點，每個點用4個比特表示，代表調制以后的一個矢量位 置(這個點擁有唯一的振幅與相位)。因此我們可以把橫軸看作是實軸，縱軸看 作虛軸。由于每個點與跟它相鄰的四個點是等距，

6、且設為2,則每個點都可用一個虛數進行表示。例如點0000可用-1-j表示，這個虛數的模就是相當于16QAM 信號的振幅，相角就相當于16QAM言號的相位。所以16QAM勺調制過程就可以用 如下語句進行描述：if A(1,b:c)=0 0 0 0B(k)=-1-1i;elseif A(1,b:c)=0 0 0 1B(k)=-3-1i;elseif A(1,b:c)=0 0 1 0B(k)=-1-3i;elseif A(1,b:c)=0 0 1 1B(k)=-3-3i;elseif A(1,b:c)=0 1 0 0B(k)=1-1i;elseif A(1,b:c)=0 1 0 1B(k)=1-3i

7、;elseif A(1,b:c)=0 1 1 0B(k)=3-1i;elseif A(1,b:c)=0 1 1 1B(k)=3-3i;elseif A(1,b:c)=1 0 0 0B(k)=-1+1i;elseif A(1,b:c)=1 0 0 1B(k)=-1+3i;elseif A(1,b:c)=1 0 1 0B(k)=-3+1i;elseif A(1,b:c)=1 0 1 1 B(k)=-3+3i;elseif A(1,b:c)=1 1 0 0B(k)=1+1i;elseif A(1,b:c)=1 1 0 1B(k)=3+1i;elseif A(1,b:c)=1 1 1 0B(k)=1+

8、3i;elseif A(1,b:c)=1 1 1 1 B(k)=3+3i;end當調制以后的信號經過信道加噪以后，我們必須對其進行解調。由于加噪了的緣故，調制以后的信號不再是原來的信號，而應該有不同。因此在解調時不能簡單的將上述過程逆轉，即不能由-1-j就判斷為0000。而應該對虛數的實部和 虛部設定一個圍后再進行判斷。這個圍邊界的選取原理我們可以借用量化的概 念，取相鄰虛數的實部的平均數和虛部的平均數。以下為16QAM軍調過程的程序語句：if (real(D( n) )<-2)&&(imag(D( n) )<-2)C(1,d:e)=0 0 1 1;elseif (

9、real(D( n)v-2)&&(imag(D( n) )<0)C(1,d:e)=0 0 0 1;elseif (real(D( n)v-2)&&(imag(D( n) )<2)C(1,d:e)=1 0 1 0;elseif (real(D( n)v-2)&&(imag(D( n)>=2)C(1,d:e)=1 0 1 1;elseif (real(D( n)<0)&&(imag(D( n) )v-2)C(1,d:e)=0 0 1 0;elseif (real(D (n )<0)&&(i

10、mag(D( n)<0)C(1,d:e)=0 0 0 0;elseif (real(D (n )<0)&&(imag(D( n) )<2)C(1,d:e)=1 0 0 0;elseif (real(D( n)<0)&&(imag(D( n) )>=2)C(1,d:e)=1 0 0 1;elseif (real(D( n)<2)&&(imag(D( n) )v-2)C(1,d:e)=0 1 0 1;elseif (real(D (n )v2)&&(imag(D( n) )v0)C(1,d:e)=0

11、 1 0 0;elseif (real(D (n )v2)&&(imag(D( n) )v2)C(1,d:e)=1 1 0 0;elseif (real(D( n)v2)&&(imag(D( n) )>=2)C(1,d:e)=1 1 1 0;elseif (real(D( n) )>=2)&&(imag(D( n) )v-2) C(1,d:e)=0 1 1 1;elseif (real(D( n)>=2)&&(imag(D( n) )<0) C(1,d:e)=0 1 1 0;elseif (real(D (

12、n )>=2)&&(imag(D( n)<2) C(1,d:e)=1 1 0 1;elseif(real(D( n) )>=2)&&(imag(D( n)>=2) C(1,d:e)=1 1 1 1;end3.2設計總程序clear all;close all;N=40000;K=4*N;%言息長度L=7*N;W=7*N/4;E=ra ndsrc(1,K,0,1);%信源B=zeros(1,W);%16QAM調制后的信號C=zeros(1,L);%16QAM 解調后的信號 num=zeros(20,1);% 誤比特數 ber=zeros(2

13、0,1);%誤比特率for SNR=1:1:20%漢明編碼A=en code(E,7,4,'ham min g/b in ary');%16QA調制過程for k=1:Wb=4*k-3;c=4*k;if A(1,b:c)=0 0 0 0B(k)=-1-1i;elseif A(1,b:c)=0 0 0 1B(k)=-3-1i;elseif A(1,b:c)=0 0 1 0 B(k)=-1-3i;elseif A(1,b:c)=0 0 1 1 B(k)=-3-3i;elseif A(1,b:c)=0 1 0 0 B(k)=1-1i;elseif A(1,b:c)=0 1 0 1 B

14、(k)=1-3i;elseif A(1,b:c)=0 1 1 0B(k)=3-1i;elseif A(1,b:c)=0 1 1 1B(k)=3-3i;elseif A(1,b:c)=1 0 0 0B(k)=-1+1i;elseif A(1,b:c)=1 0 0 1 B(k)=-1+3i;elseif A(1,b:c)=1 0 1 0 B(k)=-3+1i;elseif A(1,b:c)=1 0 1 1 B(k)=-3+3i;elseif A(1,b:c)=1 1 0 0 B(k)=1+1i;elseif A(1,b:c)=1 1 0 1 B(k)=3+1i;elseif A(1,b:c)=1

15、1 1 0 B(k)=1+3i;elseif A(1,b:c)=1 1 1 1 B(k)=3+3i;endend%言道加噪D=awg n( B,SNR);%16QA解調過程for n=1:Wd=4* n-3;e=4* n;if (real(D( n) )<-2)&&(imag(D( n) )<-2)C(1,d:e)=0 0 1 1;elseif (real(D( n)<-2)&&(imag(D( n) )<0) C(1,d:e)=0 0 0 1;elseif (real(D( n)<-2)&&(imag(D( n)

16、)<2) C(1,d:e)=1 0 1 0;elseif (real(D( n)v-2)&&(imag(D( n)>=2) C(1,d:e)=1 0 1 1;elseif (real(D( n)<0)&&(imag(D( n) )v-2) C(1,d:e)=0 0 1 0;elseif (real(D (n )<0)&&(imag(D( n) )<0) C(1,d:e)=0 0 0 0;elseif (real(D (n )<0)&&(imag(D( n) )<2) C(1,d:e)=1

17、0 0 0;elseif (real(D (n )<0)&&(imag(D( n) )>=2) C(1,d:e)=1 0 0 1;elseif (real(D (n )<2)&&(imag(D( n) )v-2) C(1,d:e)=0 1 0 1;elseif (real(D (n )<2)&&(imag(D( n) )<0) C(1,d:e)=0 1 0 0;elseif (real(D (n )<2)&&(imag(D( n) )<2) C(1,d:e)=1 1 0 0;elseif

18、(real(D (n )<2)&&(imag(D( n)>=2) C(1,d:e)=1 1 1 0;elseif (real(D( n) )>=2)&&(imag(D( n) )v-2) C(1,d:e)=0 1 1 1;elseif (real(D (n )>=2)&&(imag(D( n)<0) C(1,d:e)=0 1 1 0;elseif (real(D( n)>=2)&&(imag(D( n) )<2) C(1,d:e)=1 1 0 1;elseif (real(D( n) )&

19、gt;=2)&&(imag(D( n)>=2) C(1,d:e)=1 1 1 1;endend%漢明譯碼F=decode(C,7,4,'ham min g/b in ary');%求誤比特率nu m(SNR,1),ber(SNR,1)=biterr(F,E);end%誤比特數圖plot (nu m);figure%誤比特率圖plot(ber);3.3系統仿真結果圖0.04在不同信噪比下的誤比特率0.0350.030.0250.0150.010.00524681214161810SNR0P 0.0220圖9不同SNF下的誤比特率4系統性能分析三種在AWGNf

20、S連中的謨碼性龍比較S W'1010,5十EPSKl真課碼率 BPSkiii 碼率 FSk|f論課碼率 FEK仿其課碼率 iSQAM理論渓碼率 0 祐CAM真逞碼率曲線15E/NfdB)圖10三種調制方式在 AWG信道中的誤碼性能比較分析：將QPSK BPSK 16QAM FSK四種調制方式，包括理論值與實際值，放在同一個圖下，對他們進行對比，可以很清晰地發現，QPSI在信噪比較小時，仿真值和理論值就有了偏離，且兩者數值都比較大，當信噪比越來越大時，仿真 值成直線幾乎沒變化，而實際值的 Pe值逐漸變小，這種調制方式不是很可取； 16QAM勺性能跟QPSKS比，在低信噪比時，Pe值較大（還要