查閱數(shù)字通信的相關教材,完成OQPSK信號通信系統(tǒng)的仿真

1、OQPSK調制解調匯報人：組員：7.查閱數(shù)字通信的相關教材，完成OQPSK信號通信系統(tǒng)的仿真。OQPSK研究背景及意義 偏移正交相移鍵控(OQPSK: Offset Quadrature Phase Shift Keying)是一種調制技術，它是一種改良的QPSK，具有較高的頻譜效率和良好的光譜特性。相移鍵控調制解調技術在軟件無線電的數(shù)據傳輸系統(tǒng)中，在衛(wèi)星通信、無線環(huán)境中，有著越來越高的研究價值。 OQPSK相對于其他數(shù)字調制性能更好，更適合非線性信道下的軟件定義無線電(SDR)傳輸系統(tǒng)的實際應用。OQPSK通過模塊化的硬件平臺和軟件編程來實現(xiàn)調制和解調功能，其主要問題集中在數(shù)字解調同步算法的

2、研究上。OQPSK的基本原理OQPSK的調制方法與QPSK類似 QPSK調制原理： QPSK信號可以用正交調制的方法得到： QPSK調制可以看成兩個BPSK調制構成，輸入信號經過串并變換后分成兩個速率變?yōu)樵瓉硪话氲男蛄校缓蠼涍^極性轉換變成I(t)和Q(t)雙極性信號。I(t)和Q(t)兩路信號相互正交。QPSK解調原理：在QPSK解調中，常采用的解調方式是相干解調 已調信號為同相路相乘可得到下式正交支路為:經過低通濾波得到： 因為同相路和正交路是經過極性轉換的，即:1對應二進制數(shù)據中的1，-1對應二進制數(shù)據中的0，由上式，經過判決電路后就可以得到如下表所示的結果。QPSK的星座圖 QPSK信

3、號可以用星座圖表示，也成為矢量圖。它表明了各個符號(11,00,01，10)間的相位和幅度關系，四個雙比特符號分別表示QPSK信號的四個相位，相鄰兩個相位相互正交。在這些相位偏移關系中，01對應11, 00對應01等，相位偏移90度.或者是00對11相位偏移180度。OQPSK調制解調原理： 與QPSK相同的是相位關系，即：把輸入信號分為兩路，然后正交調制。所不同的是，同相和正交兩支路的碼流在時間上錯開了半個周期。因為兩支路碼元上偏移了半個周期，每次只能有一路可能發(fā)生極性翻轉。所以，OQPSK信號能跳變的相位只能是0度、正負90度，不會出現(xiàn)180度的跳變，OQPSK與QPSK相比，信號的包絡波

4、動幅度有限，經過限幅放大后的頻帶范圍也要小，所以性能也更加優(yōu)良。其星座映射圖如下： OQPSK信號的產生方法相對于QPSK的產生方法僅僅在Q支路上增加了一個Ts/2(半個周期)的延時電路。 OQPSK信號的解調同QPSK解調原理基本相同，因為在調制時Q支路信號在時間上延遲了半個周期，所以在解調端的Q支路上也應該延時Ts/2再進行抽樣判決，從而保證對兩個支路能進行正確的交錯抽樣，然后經過并/串變換來恢復原始信號。OQPSK仿真產生二進制信源m=randn(1,N); %產生1到n的隨機數(shù)d=sign(m); % 將大于0的變?yōu)?小于0的變?yōu)?1dd=sigexpand(d,fc*N_sample

5、); %序列d的周期變?yōu)門sgt=ones(1,fc*N_sample); d_NRZ=conv(dd,gt);%卷積產生基帶信號figure(5);subplot(2,2,1); plot(t,dd(1:Lt);axis(tx1,tx2,ty1,ty2);figure(1);subplot(2,4,1);plot(t,d_NRZ(1:Lt);axis(tx1,tx2,ty1,ty2);xlabel(時間(S);ylabel(幅度);title(基帶信號時域波形圖); grid;f,d_NRZf=T2F(t,d_NRZ(1:Lt);% 進行傅里葉變換 figure(1);subplot(2,4

6、,5);plot(f,10*log10(abs(d_NRZf).2/T); axis(fx1,fx2,fy1,fy2); xlabel(頻率(Hz);ylabel(功率譜密度(dB/Hz); title(基帶信號功率譜圖); grid;figure(4);subplot(2,2,1);plot(t,d_NRZ(1:Lt);axis(tx1,tx2,ty1,ty2); xlabel(時間(S); ylabel(幅度);title(基帶信號時域波形圖); grid;figure(4);subplot(2,2,2);plot(f,10*log10(abs(d_NRZf).2/T);axis(fx1,

7、fx2,fy1,fy2);xlabel(頻率(Hz);ylabel(功率譜密度(dB/Hz); title(基帶信號功率譜圖); grid; 串并轉換d1=; d2=; for i=1:Nif rem(i,2)=1d1(i+1)/2)=d(i); elsed2(i/2)=d(i); end enddd1=sigexpand(d1,2*fc*N_sample); %功能同上 gt1=ones(1,2*fc*N_sample); d_NRZ1=conv(dd1,gt1);figure(1);subplot(2,4,2);plot(t,d_NRZ1(1:Lt); axis(tx1,tx2,ty1,t

8、y2); xlabel(時間(S); ylabel(幅度);title(上支路基帶信號時域波形圖); grid;f1,d_NRZ1f=T2F(t,d_NRZ1(1:Lt); figure(1);subplot(2,4,6);plot(f1,10*log10(abs(d_NRZ1f).2/T); axis(fx1,fx2,fy1,fy2); xlabel(頻率(Hz);ylabel(功率譜密度(dB/Hz);title(上支路基帶信號功率譜圖); grid;figure(3);subplot(2,4,1);plot(t,d_NRZ1(1:Lt); axis(tx1,tx2,ty1,ty2); x

9、label(時間(S); ylabel(幅度);title(上支路基帶信號時域波形圖); grid;dd2=sigexpand(d2,2*fc*N_sample); gt1=ones(1,2*fc*N_sample); d_NRZ2=conv(dd2,gt1); d_NRZ2DLY=-ones(1,N_sample*fc),d_NRZ2(1:end-N_sample*fc); %進行延時，在前面添-1figure(1);subplot(2,4,3);plot(t,d_NRZ2DLY(1:Lt); axis(tx1,tx2,ty1,ty2); xlabel(時間(S);ylabel(幅度);ti

10、tle(下支路基帶信號時域波形圖); grid;f2,d_NRZ2f=T2F(t,d_NRZ2DLY(1:Lt); figure(1);subplot(2,4,7);plot(f2,10*log10(abs(d_NRZ2f).2/T); axis(fx1,fx2,fy1,fy2); xlabel(頻率(Hz);ylabel(功率譜密度(dB/Hz);title(下支路基帶信號功率譜圖); grid;figure(3);subplot(2,4,5);plot(t,d_NRZ2DLY(1:Lt); axis(tx1,tx2,ty1,ty2); xlabel(時間(S); ylabel(幅度);ti

11、tle(下支路基帶信號時域波形圖); grid; 載波h1t=A*cos(2*pi*fc*t); h2t=A*sin(2*pi*fc*t);figure(1);subplot(2,4,4); plot(t,h1t);axis(tx1,tx2,ty1,ty2); xlabel(時間(S); ylabel(幅度);title(載波信號時域波形圖); grid;f3,h1tf=T2F(t,h1t); figure(1);subplot(2,4,8);plot(f3,10*log10(abs(h1tf).2/T);% p=2/T *10*log10(abs(h1tf)為求功率譜的公式 axis(fx1

12、,fx2,fy1,fy2); xlabel(頻率(Hz);ylabel(功率譜密度(dB/Hz); title(載波信號功率譜圖);grid; 生成OQPSK信號s_qpsk1=d_NRZ1(1:Lt).* h1t; %上下支路分別調制 s_qpsk2=d_NRZ2DLY(1:Lt).* h2t; figure(2);subplot(2,2,1); plot(t,s_qpsk1);axis(tx1,tx2,ty1,ty2); xlabel(時間(S); ylabel(幅度);title(上支路頻帶信號時域波形圖); grid;f4,s_qpsk1f=T2F(t,s_qpsk1); figure

13、(2);subplot(2,2,3);plot(f4,10*log10(abs(s_qpsk1f).2/T); axis(fx1,fx2,fy1,fy2); xlabel(頻率(Hz);ylabel(功率譜密度(dB/Hz);title(上支路頻帶信號功率譜圖); grid;figure(2);subplot(2,2,2); plot(t,s_qpsk2);axis(tx1,tx2,ty1,ty2); xlabel(時間(S); ylabel(幅度);title(下支路頻帶信號時域波形圖); grid;f5,s_qpsk2f=T2F(t,s_qpsk2); figure(2);subplot(

14、2,2,4);plot(f5,10*log10(abs(s_qpsk2f).2/T); axis(fx1,fx2,fy1,fy2); xlabel(頻率(Hz);ylabel(功率譜密度(dB/Hz);title(下支路頻帶信號功率譜圖); grid;s_oqpsk=s_qpsk1+s_qpsk2; %兩路信號相加得到調制后的信號 figure(8);subplot(2,2,1);plot(t,s_oqpsk);axis(tx1,tx2,ty1,ty2);xlabel(時間(S);ylabel(幅度);title(已調信號時域波形圖);grid;f6,s_oqpskf=T2F(t,s_oqps

15、k);figure(8);subplot(2,2,3);plot(f6,10*log10(abs(s_oqpskf).2/T);axis(fx1,fx2,fy1,fy2);xlabel(頻率(Hz);ylabel(功率譜密度(dB/Hz);title(已調信號功率譜圖);grid; 信道加入高斯白噪聲進行接收解調 產生高斯白噪聲 m=1; p1=-10; noise = wgn(m,Lt,p1);接收信號y_qpsk = s_oqpsk + noise;figure(8);subplot(2,2,2);plot(t,y_qpsk);axis(tx1,tx2,ty1,ty2);xlabel(時間

16、(S);ylabel(幅度);title(接收信號時域波形圖 );grid;f7,y_qpskf=T2F(t,y_qpsk);figure(8);subplot(2,2,4);plot(f7,10*log10(abs(y_qpskf).2/T);axis(fx1,fx2,fy1,fy2);xlabel(頻率(Hz);ylabel(功率譜密度(dB/Hz);title(接收信號功率譜圖 Pn=-10dB);grid;%t,ny_qpsk=bpf(f7,y_qpskf,1,8); 相干解調通過乘法器1r_qpsk1 = y_qpsk .* h1t;figure(3);subplot(2,4,2);

17、plot(t,r_qpsk1);axis(tx1,tx2,ty1,ty2);xlabel(時間(S);ylabel(幅度);title(上支路通過乘法器信號時域波形圖);grid;%通過低通濾波器f8,r_qpsk1f=T2F(t,r_qpsk1);B1=1;t1,r_qpsk11=lpf(f8,r_qpsk1f,B1);figure(3);subplot(2,4,3);plot(t1,r_qpsk11)axis(tx1,tx2,ty1,ty2);xlabel(時間(S);ylabel(幅度);title(上支路通過低通濾波器信號時域波形圖);grid;%抽樣判決dd11=r_qpsk11(f

18、c*N_sample:2*fc*N_sample:end);dd22=sign(dd11); %判決dd222=sigexpand(dd22,2*fc*N_sample);d_NRZ11=conv(dd222,gt1); %得到上支路信號d_NRZ11DLY=-ones(1,N_sample*fc),d_NRZ11(1:end-N_sample*fc); %路信號延時figure(3);subplot(2,4,4);plot(t,d_NRZ11DLY(1:Lt);axis(tx1,tx2,ty1,ty2);xlabel(時間(S);ylabel(幅度);title(上支路抽樣判決后信號時域波形

19、圖);grid;通過乘法器2r_qpsk2 = y_qpsk .* h2t;figure(3);subplot(2,4,6);plot(t,r_qpsk2);axis(tx1,tx2,ty1,ty2);xlabel(時間(S);ylabel(幅度);title(下支路通過乘法器信號時域波形圖);grid;%通過低通濾波器f9,r_qpsk2f=T2F(t,r_qpsk2);t2,r_qpsk21=lpf(f9,r_qpsk2f,B1);figure(3);subplot(2,4,7);plot(t2,r_qpsk21)axis(tx1,tx2,ty1,ty2);xlabel(時間(S);yla

20、bel(幅度);title(下支路通過低通濾波器信號時域波形圖);grid;%抽樣判決dd33=r_qpsk21(fc*N_sample:2*fc*N_sample:end);dd44=sign(dd33); %判決dd444=sigexpand(dd44,2*fc*N_sample);d_NRZ21=conv(dd444,gt1)figure(3);subplot(2,4,8);plot(t,d_NRZ21(1:Lt)axis(tx1,tx2,ty1,ty2);xlabel(時間(S);ylabel(幅度);title(下支路抽樣判決后信號時域波形圖);grid;dd22DLY=d_NRZ11DLY(fc*N_sample:2*fc*N_sample:end);并串轉換ddd=;for s=1:N/2ddd(2*s-1)=dd22DLY(s);ddd(2*s)=dd44(s);endrddd=sigexpand(ddd,fc*N_sample);%解調信號輸出r_qpsk=conv(rddd,gt);qr_qpsk