【正文】 PSK的調(diào)制框圖如圖513所示。圖513 產(chǎn)生OQPSK信號的框圖k 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 +1 +1 +1 +1 1 1 +1 +1 1 1 +1 1 +1 +1 +1 1圖514 OQPSK的同相及正交支路基帶波形圖514中的OQPSK信號表示為： 式(511)其中， 式(512) 式(513)是成形濾波器沖擊響應(yīng)，若是不歸零矩陣脈沖，可表示為： 式(514)關(guān)于OQPSK及QPSK信號波形示于圖515及圖516。OQPSK最佳接收的平均誤比特率與QPSK相同： 式(515) 程序仿真圖519 QPSK及OQPSK經(jīng)帶通濾波器后波形及功率譜密度 經(jīng)過帶通濾波器后，QPSK信號的包絡(luò)起伏明顯變大，而OQPSK信號的包絡(luò)起伏相對小得多，再經(jīng)過非線性功率放大器后，此時QPSK信號的功率譜旁瓣比OQPSK信號大。在這次設(shè)計(jì)中，我仿真出了QPSK在加性高斯白噪聲下的最佳接收，比較并分析了QPSK及OQPSK的性能。在數(shù)字頻帶通信系統(tǒng)中，QPSK調(diào)制解調(diào)技術(shù)已經(jīng)比較成熟，在實(shí)際應(yīng)用中也發(fā)揮著極為重要的作用。使我能夠在短時間內(nèi)有很大的進(jìn)步，順利完成畢業(yè)設(shè)計(jì)的各項(xiàng)工作。G ml=2。 SNR=0。 datanrz=data.*21。 for q=1:nb data0((q1)/delta_T+1:q/delta_T)=data(q)。% QPSK信號的調(diào)制ich=zeros(1,nb/delta_T/2)。 qch=zeros(1,nb/2/delta_T)。 s=idata1+qdata1。Ray_idata=idata1.*ray_ich。 % 通過高斯信道之后的信號s11=abs(fft(s1))。 % 這里面其實(shí)隱藏了一個串并轉(zhuǎn)換的過程qdata2=s1.*b。 if sum(idata2((n1)/delta_T+1:n/delta_T))=0 idata3(n)=1。 else qdata3(n)=0。 if Awgn_ichsum(n)=0 idata4(n)=1。 else qdata4(n)=0。 % 存放偶數(shù)位%為了顯示，將它變成波形信號（即傳輸一個1代表單位寬度的高電平）demodata1=zeros(1,nb/delta_T)。 % Ray_idata3=zeros(1,nb/2)。 % end% if Ray_sum(Ray_qdata2((n1)/delta_T+1:n/delta_T))=0% Ray_qdata3(n)=1。for n=1:nb/2 Ray_ichsum(n)=sum(idata2((n1)/delta_T+1:n/delta_T))*delta_T。 if Ray_qchsum(n)=0 Ray_qdata4(n)=1。 % 存放奇數(shù)位Ray_demodata(2:ml:nb)=Ray_qdata4。和cm_sm33(）函數(shù)%%聲明： 函數(shù)聲明在另外倆個M文件中%%作用： cm_sm32()用于瑞利信道誤碼率的計(jì)算%% cm_sm33()用于高斯信道誤碼率的計(jì)算%% ecoh on/off 作用在于決定是否顯示指令內(nèi)容%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%SNRindB1=0:1:8。 smld_symbol_ray_err_prb(i)=ps。% 高斯信道 echo on。 disp([ps1,pb1])。for i=1:length(SNRindB2), SNR=exp(SNRindB2(i)*log(10)/10)。end。fs39。subplot(211)。axis([0 20000 2 2])。,fs),title(39。plot(s),title(39。figure(3)subplot(212)。解調(diào)輸出功率譜密度39。解調(diào)輸出39。psd(h,s1,39。)。)。高斯噪聲曲線39。r+39。hold on。實(shí)際值（接收端）39。*39。r39。hold on。hold on。*39。)。 legend(39。,39。)。M=4。N_num=100。py1f=zeros(1,length(t))。 gt=ones(1,fc*N_sample)。 b2=conv(s2,gt)。 s2_delay=[ones(1,N_sample*fc/2) s2(1:endN_sample*fc/2)]。 [t y1]=bpf(f,y1f,fc1/Ts,fc+1/Ts)。xlabel(39。QPSK波形39。經(jīng)過帶通后的波形39。xlabel(39。OQPSK波形39。經(jīng)過帶通后的波形39。 [f y1f]=T2F(t,y1)。 %OQPSK不同段信號功率譜密度相加endpy1f=py1f/100 %QPSK100段功率譜密度平均py2f=py2f/100 %OQPSK100段功率譜密度平均subplot(223)。)。title(39。subplot(224)。)。title(39。figure。)。OQPSK功率譜密度(dB/Hz)39。xlabel(39。)。QPSK功率譜密度(dB/Hz)39。xlabel(39。 py1f=py1f+abs(y1f).^2/T。 %經(jīng)過非線性電路 y1=*tanh(2*y1)。 axis([5 15 ])。)。 subplot(222)。 axis([5 15 ])。)。 subplot(221)。 %經(jīng)過帶通后，在經(jīng)過非線性電路 [f y1f]=T2F(t,st_qpsk)。 s2=b2(1:length(s2))。 s2=sigexpand(d2,fc*N_sample)。 %功率譜密度2for PL=1:100 %輸入100段N_num個碼字的波形，為了使功率譜密度看起來更加平滑， %可以取這100段信號功率譜密度的平均 d1=sign(randn(1,N_num))。t=0:dt:N_num*Tsdt。fc=10。 OQPSKclear all。,39。,39。BER39。 xlabel(39。hold on。r*39。誤碼率曲線39。hold on。axis([2 2 2 2])。理論值（發(fā)射端）39。QPSK信號星座圖（Awgn）39。axis([0 2000 5 5])。subplot(311)。plot(s1),title(39。,fs),title(39。axis([0 20000 2 2])。subplot(211)。fs39。)。)。psd(h,s,39。基帶信號39。基帶信號功率譜密度39。 % 類似于C語言的宏定義，方便以下的調(diào)用 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 第一部分（理想）figure(1)subplot(212)。 % 高斯噪聲理論誤碼率 theo_err_ray_prb(i)=*(11/sqrt(1+1/SNR))。end。 smld_bit_awgn_err_prb(i)=pb1。 echo off。% 瑞利衰落信道 for i=1:length(SNRindB1), [pb,ps]=cm_sm32(SNRindB1(i))。 % 創(chuàng)建一個1*nb/delta_T的零矩陣for q=1:nb Ray_demodata1((q1)/delta_T+1:q/delta_T)=Ray_demodata(q)。 endend % 將判決之后的數(shù)據(jù)存放進(jìn)數(shù)組Ray_demodata=zeros(1,nb)。 else Ray_idata4(n)=0。% end% end % 為了顯示星座圖,將信號進(jìn)行處理Ray_idata4=zeros(1,nb/2)。% 抽樣判決的過程，與0作比較，data=0,則置1，否則置0% for n=1:nb/2% if Ray_sum(Ray_idata2((n1)/delta_T+1:n/delta_T))=0% Ray_idata3(n)=1。 % 將極性碼變成對應(yīng)的波形信號end % 累計(jì)誤碼數(shù)% abs(demodatadata)求接收端和發(fā)射端% 數(shù)據(jù)差的絕對值，累計(jì)之后就是誤碼個數(shù)Awgn_num_BER=sum(abs(demodatadata)) %%%%%%%%%%%%%%%%%%%% 解調(diào)部分（瑞利+高斯）Ray_idata2=Ray_s.*a。demodata(1:ml:(nb1))=idata3。 end Awgn_qchsum(n)=sum(qdata2((n1)/delta_T+1:n/delta_T))*delta_T。qdata4=zeros(1,nb/2)。 end% A2(n)=sum(qdata2((n1)/delta_T+1:n/delta_T))。 % 建立1*nb數(shù)組，以存放解調(diào)之后的信號qdata3=zeros(1,nb/2)。 % 高斯噪聲曲線%%%%%%%%%%%%%%%%%%%%%%%%%%Awgn_s=awgn(Ray_s,SNR)。Ray_s=Ray_idata+Ray_qdata。 % 瑞利衰落信道ray_ich=raylrnd(,1,nb/2/delta_T)。endfor jj=1:N/2 b(jj)=sqrt(2/T)*sin(2*pi*fc*t(jj))。endfor ii=1:N/2 a(ii)=sqrt(2/T)*cos(2*pi*fc*t(ii))。% 串并轉(zhuǎn)換，將奇偶位數(shù)據(jù)分開idata=datanrz(1:ml:(nb1))。 for q=1:nb data1((q1)/delta_T+1:q/delta_T)=datanrz(q)。 N=length(t)。 delta_T=T/200。T=1。 感謝在大學(xué)四年中所有關(guān)心我的任課老師，同學(xué)和朋友們，感謝我的父母，因?yàn)橛兴麄兊闹С郑也拍懿粩嘞蚯斑~進(jìn)。致謝首先，我要感謝北京郵電大學(xué)世紀(jì)學(xué)院給我提供這么一個實(shí)踐的機(jī)會，為我們在校本科生提供如圖書館之類的各種學(xué)習(xí)場所，及可以完成實(shí)踐的學(xué)習(xí)工作條件。OQPSK的性能優(yōu)于QPSK，結(jié)果與預(yù)想中的一樣。本論文的題目是基于MATLAB數(shù)字頻帶通信系統(tǒng)仿真。 OQPSK信號的最佳接收及誤比特率圖517 在加性白高斯噪聲干擾下OQPSK的最佳接收框圖（1）OQPSK信號的平均功率譜密度由于OQPSK也可以看為是由同相及正交支路的2PSK信號的疊加，所以O(shè)QPSK的功率譜與QPSK的相同。它與QPSK不同的是，他的正交支路碼元與同相支路碼元在時間上要偏移一個比特間隔。恒包絡(luò)技術(shù)所產(chǎn)生的已調(diào)波經(jīng)過發(fā)送帶限后，當(dāng)通過非線性部件時，只產(chǎn)生很小的頻譜擴(kuò)展。圖512 QPSK信號限帶濾波前后的波形若將此包絡(luò)起伏很大的限帶QPSK信號再進(jìn)行硬限幅或非線性功率放大，雖然此已調(diào)信號包絡(luò)的起伏可減弱，但卻使非線性功率放大后的信號功率譜旁瓣增生，頻率又被擴(kuò)展，其旁瓣將會干擾鄰近頻道的信號，這是不希望出現(xiàn)的現(xiàn)象。瑞利衰落信道下的誤碼率比高斯信道下的誤碼率高。并且可以看出調(diào)制解調(diào)的結(jié)果沒有誤碼。于是，得到I支路及Q支路的平均錯判概率為： 式(57)其中，是平均比特能量。下面，采用第二種方法來進(jìn)行計(jì)算。圖55表示在2PSK及QPSK的二進(jìn)制信息速率相同時，QPSK信號的平均功率譜密度的主瓣寬度是2PSK平均功率譜主瓣寬度的一半。 QPSK信號的平均功率譜密度由于QPSK信號時由兩正交載波的2PSK線性疊加而成，所以QPSK信號的平均功率譜密度是同相支路及正交支路2PSK信號平均功率譜密度的線性疊加。k 1 2 3 4 5 6 7 8 9 10 11 12 13 14n +1 +1 1 +1 +1 1 +1 1 +1 +1 1 1 1 +12n +1 +1 1 1 +1 1 +12n1 +1 1 +1 +1 +1 1 1