QPSK通信系統(tǒng)的設計與仿真

1、 摘要 介紹了數(shù)字通信中的qpsk調(diào)制解調(diào)的原理,通過用matlab編寫腳本程序?qū)psk通信系統(tǒng)的發(fā)射和接收過程的具體實現(xiàn)進行模擬仿真,繪出信號在理想信道和加噪信道中模擬傳輸時的時域圖,并對各模塊進行了頻譜分析,所得到的結果與理論基本相符,對于理解qpsk系統(tǒng)的性能并在系統(tǒng)的實際應用上作進一步的設計,提供了有效的參考依據(jù)。 通過利用matlab軟件sumlink實現(xiàn)了qpsk通信系統(tǒng)的仿真，完成了qpsk通信系統(tǒng)的調(diào)制解調(diào)過程的仿真實現(xiàn)，使接收端能夠準確地接收到來自發(fā)放的信息。qpsk調(diào)制方式在通信工程中的應用十分廣泛，其誤碼率隨信噪比的增加而減少并最終可能為零。在qpsk調(diào)制方式以后，還會

2、出現(xiàn)進制更多的調(diào)制方式。而我們著重要解決的問題也從如何提高相位譜利用率轉(zhuǎn)變?yōu)槿绾螠p少誤差以及提高傳送速率。闡述qpsk調(diào)制解調(diào)的實現(xiàn)過程，并運用軟件實現(xiàn)手段對信號變換過程加以分析，希望有所收獲。關鍵詞：數(shù)字通信 四相相位鍵控(quadrature phase shift keying,qpsk) matlab目錄前言1一、 qpsk簡介22.1 qpsk調(diào)制原理32.2 qpsk解調(diào)原理5三、 matlab仿真結果63.1 qpsk的產(chǎn)生及加噪仿真63.2 qpsk解調(diào)前后信號仿真63.3 qpsk碼元恢復73.4誤碼率分析9四、 程序清單10五、 總結17六、 參考文獻18七、致謝19前言

3、隨著數(shù)字技術的飛速發(fā)展與數(shù)字器件的廣泛使用,數(shù)字信號處理在通信系統(tǒng)中的應用已經(jīng)越來越重要。數(shù)字信號傳輸系統(tǒng)分為基帶傳輸系統(tǒng)和頻帶傳輸系統(tǒng),頻帶傳輸系統(tǒng)也叫數(shù)字調(diào)制系統(tǒng),該系統(tǒng)對基帶信號進行調(diào)制,使其頻譜搬移到適合在信道(一般為帶通信道)上傳輸?shù)念l帶上。數(shù)字調(diào)制和模擬調(diào)制一樣都是正弦波調(diào)制,即被調(diào)制信號都為高頻正弦波。數(shù)字調(diào)制信號又稱為鍵控信號,數(shù)字調(diào)制過程中處理的是數(shù)字信號,而載波有振幅、頻率和相位3個變量,且二進制的信號只有高低電平兩個邏輯量即1和0,所以調(diào)制的過程可用鍵控的方法由基帶信號對載頻信號的振幅、頻率及相位進行調(diào)制,最基本的方法有3種:正交幅度調(diào)制(qam) 、頻移鍵控( fsk)

4、 、相移鍵控( psk) 。根據(jù)所處理的基帶信號的進制不同分為二進制和多進制調(diào)制(m進制) 。多進制數(shù)字調(diào)制與二進制相比,其頻譜利用率更高。在數(shù)字通信、數(shù)字視頻廣播、數(shù)字衛(wèi)星廣播等領域中,廣泛采用m進制的調(diào)制方式 1 。理論分析指出,在恒參信道中利用psk方式可以獲得最佳接收性能,在接收機噪聲的作用下, psk誤碼率最低。本文用matlab軟件對m = 4的四相相位調(diào)制解調(diào)(qpsk)進行模擬,并對各模塊進行了仿真分析。一、 qpsk簡介 在多相調(diào)制中，可以用多種相位差來表示數(shù)字數(shù)據(jù)信息。把要傳送的二進制信息序列每k比特編為一組，那么便構成了k比特碼元，每一k比特碼元都有2k種不同的狀態(tài)，要用

5、2k 種不同的相位或相位差來表示。由于k比特碼元包括的信息量是二進制碼元所含信息量的k倍，所以，多相調(diào)制系統(tǒng)與二相調(diào)制系統(tǒng)相比，在碼元傳輸速率相同時，信息傳輸速率更高。例如，四相調(diào)制系統(tǒng)的信息傳輸速率是二相調(diào)制系統(tǒng)的兩倍。另外，在系統(tǒng)信息傳輸速率的情況下，由于多相調(diào)制系統(tǒng)的碼元速率比二相調(diào)制要低，因而多相調(diào)制時，信號碼元持續(xù)時間比二相調(diào)制時長。碼元時間長，則對應的帶寬就窄，因此壓縮了信號的頻帶，可以叫囂由于新到特性引起的碼間串擾的影響。由此可見，多相調(diào)制是提高數(shù)據(jù)通信有效性、可靠性的一種途徑。 由于四相絕對移相信號可以看成兩個正交2psk信號只能采用相干檢測法解調(diào)， 4psk信號也只能采用兩個

6、正交的相干載波分別來檢測這兩個分量a和b。在四相絕對移相調(diào)制系統(tǒng)中，實現(xiàn)起來最困難的問題在于接收端如何產(chǎn)生一標準的相干載波。通常采用的方法是四倍頻-四分頻法，然后經(jīng)過四次分頻，便可以得到所有的相干載波輸出。這種方法也會由于分頻起始點不同而使提取的相干載波的相位出現(xiàn)四種不確定性。因此，在實際中，通常采用四相相對移相調(diào)制。 二、qpsk調(diào)制解調(diào)原理2.1 qpsk調(diào)制原理 四相相移鍵控信號簡稱“qpsk”。它分為絕對相移和相對相移兩種。由于絕對相移方式存在相位模糊問題，所以在實際中主要采用相對移相方式qdpsk。它具有一系列獨特的優(yōu)點，目前已經(jīng)廣泛應用于無線通信中，成為現(xiàn)代通信中一種十分重要的調(diào)制

7、解調(diào)方式。在數(shù)字信號的調(diào)制方式中qpsk四相移鍵控是目前最常用的一種衛(wèi)星數(shù)字信號調(diào)制方式,它具有較高的頻譜利用率、較強的抗干擾性、在電路上實現(xiàn)也較為簡單。 qpsk四相絕對相移鍵控，用載波信號的4個初始相位對應4進制碼元。因此，對于輸入的二進制數(shù)字序列每2比特分為一組，稱為雙比特碼元ab，然后用4種不同的載波相位分別表征這4種數(shù)字碼元。按ab分組與載波相位的對應關系有/2體系(也稱a方式)與/4體系(也稱b方式)，具體定義如表所示：雙比特碼元ab載波相位ab/2體系（a方式）/4體系(b方式)000（）5/4（225°）10/2（）7/4（315°）11（）/4（45

8、76;）013/2（）3/4（135°） 圖2-1 qpsk信號雙比特碼元與相位矢量關系圖 圖2-2 四相相移調(diào)制是利用載波的四種不同相位差來表征輸入的數(shù)字信息，是四進制移相鍵控。qpsk是在m=4時的調(diào)相技術，它規(guī)定了四種載波相位，分別為45°，135°，225°，315°，調(diào)制器輸入的數(shù)據(jù)是二進制數(shù)字序列，為了能和四進制的載波相位配合起來，則需要把二進制數(shù)據(jù)變換為四進制數(shù)據(jù)，這就是說需要把二進制數(shù)字序列中每兩個比特分成一組，共有四種組合，即00，01，10，11，其中每一組稱為雙比特碼元。每一個雙比特碼元是由兩位二進制信息比特組成，它們分別

9、代表四進制四個符號中的一個符號。qpsk中每次調(diào)制可傳輸2個信息比特，這些信息比特是通過載波的四種相位來傳遞的。解調(diào)器根據(jù)星座圖及接收到的載波信號的相位來判斷發(fā)送端發(fā)送的信息比特。數(shù)字調(diào)制用“星座圖”來描述，星座圖中定義了一種調(diào)制技術的兩個基本參數(shù)：（1）信號分布；（2）與調(diào)制數(shù)字比特之間的映射關系。星座圖中規(guī)定了星座點與傳輸比特間的對應關系，這種關系稱為“映射”，一種調(diào)制技術的特性可由信號分布和映射完全定義，即可由星座圖來完全定義。同時qpsk信號可以看作兩個載波正交2psk信號的合成，下圖表示qpsk正交調(diào)制器。圖2-3 qpsk調(diào)制系統(tǒng)原理圖2.2 qpsk解調(diào)原理 數(shù)字調(diào)制用“星座圖”

10、來描述，星座圖中定義了一種調(diào)制技術的兩個基本參數(shù)：（1）信號分布；（2）與調(diào)制數(shù)字比特之間的映射關系。星座圖中規(guī)定了星座點與傳輸比特間的對應關系，這種關系稱為“映射”，一種調(diào)制技術的特性可由信號分布和映射完全定義，即可由星座圖來完全定義。四相相移調(diào)制是利用載波的四種不同相位差來表征輸入的數(shù)字信息，是四進制移相鍵控。qpsk是在m=4時的調(diào)相技術，它規(guī)定了四種載波相位，分別為45°，135°，225°，315°，調(diào)制器輸入的數(shù)據(jù)是二進制數(shù)字序列，為了能和四進制的載波相位配合起來，則需要把二進制數(shù)據(jù)變換為四進制數(shù)據(jù)，這就是說需要把二進制數(shù)字序列中每兩個比特分

11、成一組，共有四種組合，即00，01，10，11，其中每一組稱為雙比特碼元。每一個雙比特碼元是由兩位二進制信息比特組成，它們分別代表四進制四個符號中的一個符號。qpsk中每次調(diào)制可傳輸2個信息比特，這些信息比特是通過載波的四種相位來傳遞的。解調(diào)器根據(jù)星座圖及接收到的載波信號的相位來判斷發(fā)送端發(fā)送的信息比特。由qpsk信號的調(diào)制可知，對它的解調(diào)可以采用與2psk信號類似的解調(diào)方法進行解調(diào)。解調(diào)原理圖如下所示，同相支路和正交支路分別采用相干解調(diào)方式解調(diào)，得到和，經(jīng)過抽樣判決和并/串交換器，將上下支路得到的并行數(shù)據(jù)恢復成串行數(shù)據(jù)。圖2-4 qpsk解調(diào)系統(tǒng)原理圖3、 matlab仿真結果3.1 qps

12、k的產(chǎn)生及加噪仿真 圖3-1 通過matlab產(chǎn)生兩路相互正交的2psk信號，信號疊加后的到合成序列，對合成序列進行加噪處理。3.2 qpsk解調(diào)前后信號仿真圖3-23.3 qpsk碼元恢復然后相干解調(diào)中，正交路和同相路分別設置兩個相關器(或匹配濾波器)，得到i(t)和q(t)，經(jīng)電平判決和并一串變換后即可恢復原始信息。當然，如果調(diào)制端是差分編碼的，那么解調(diào)中并串變換后還需一個差分解碼。假如已調(diào)信號分別為同相路和正交路，為載波頻率，那么相干解調(diào)后，同相路相乘可得 正交路為： 經(jīng)過低通濾波后，可得： 在完成通過濾波器之后，將碼元再通過并串轉(zhuǎn)換，就可以恢復出原始信號的碼元序列。圖3-33.4誤碼率

13、分析圖3-4有圖示可知，qpsk調(diào)制解調(diào)有比較低的誤碼率，所以因其抗干擾性能強、誤碼性能好、頻譜利用率高等優(yōu)點,廣泛應用于數(shù)字微波通信系統(tǒng)、數(shù)字衛(wèi)星通信系統(tǒng)、寬帶接入、移動通信及有線電視系統(tǒng)之中。4、 程序清單clear allclose all t=-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; endend t1=0:0.01:8-0.01;t2=0:0

14、.01:7-0.01;t3=-1:0.01: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; endend f=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

15、(2*pi*f1*t);q=x2.*sin(2*pi*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; f

16、igure(1)subplot(4,1,1);plot(t3,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('加入噪聲');clear allclos

17、e 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_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

18、_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; 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

19、*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+19) q(i:1:i+19);end; bit_recover=;for i=1:20:20000 if sum(data_recover(i:i+19)&g

20、t;0 data_recover_a(i:i+19)=1; bit_recover=bit_recover 1; else data_recover_a(i:i+19)=-1; bit_recover=bit_recover -1; endend error=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; endend p=error/1000; figure(1)subpl

21、ot(2,1,1);plot(t2,ddd2);axis(0 100 -2 2);title('原序列');subplot(2,1,2);plot(t2,data_recover_a);axis(0 100 -2 2);title('解調(diào)后序列');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(d

22、ata_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;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.*s

23、in(2*pi*f1*t1);qpsk_rc=(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_rec

24、over=data_recover i(i:1:i+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);end figure(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 -