【正文】 于4PSK信號可以看作是兩個載波正交的2PSK信號的合成，因此，對4PSK信號的解調可以采用與2PSK信號類似的解調方法進行。圖829是B方式4PSK信號相干解調器的組成方框圖。圖中兩個相互正交的相干載波分別檢測出兩個分量和，然后，經(jīng)并/串變換器還原成二進制雙比特串行數(shù)字信號，從而實現(xiàn)二進制信息恢復。此法也稱為極性比較法。圖8294PSK信號的相干解調若解調4PSK信號（A方式），只需適當改變相移網(wǎng)絡。在2PSK信號相干解調過程中會產(chǎn)生“倒”即“180176。相位模糊”現(xiàn)象。同樣，對于4PSK信號相干解調也會產(chǎn)生相位模糊問題，并且是0176。、90176。、180176。和270176。四個相位模糊。4PSK基于MATLAB的仿真% 調相法clear allclose allt=[1::]。tt=length(t)。x1=ones(1,800)。for i=1:tt if (t(i)=1 amp。 t(i)=1) | (t(i)=5amp。 t(i)=7)。 x1(i)=1。 else x1(i)=1。 endendt1=[0::]。t2=0::。t3=1::。t4=0::。tt1=length(t1)。x2=ones(1,800)。for i=1:tt1 if (t1(i)=0 amp。 t1(i)=2) | (t1(i)=4amp。 t1(i)=8)。 x2(i)=1。 else x2(i)=1。 endendf=0::1。xrc=+*cos(pi*f)。y1=conv(x1,xrc)/。y2=conv(x2,xrc)/。n0=randn(size(t2))。f1=1。i=x1.*cos(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。figure(1)subplot(4,1,1)。plot(t3,i_rc)。axis([1 8 1 1])。ylabel(39。a序列39。)。subplot(4,1,2)。plot(t4,q_rc)。axis([1 8 1 1])。ylabel(39。b序列39。)。subplot(4,1,3)。plot(t2,QPSK_rc)。axis([1 8 1 1])。ylabel(39。合成序列39。)。subplot(4,1,4)。plot(t2,QPSK_rc_n1)。axis([1 8 1 1])。ylabel(39。加入噪聲39。) 。仿真結果1：data_I = 2*bit_I+1。data_Q = 2*bit_Q+1。data_I1=repmat(data_I39。,20,1)。data_Q1=repmat(data_Q39。,20,1)。for i=1:1e4 data_I2(i)=data_I1(i)。 data_Q2(i)=data_Q1(i)。end。f=0::1。xrc=+*cos(pi*f)。data_I2_rc=conv(data_I2,xrc)/。data_Q2_rc=conv(data_Q2,xrc)/。f1=1。t1=0::1e3+。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。% 解調I_demo=QPSK_rc_n0.*cos(2*pi*f1*t1)。Q_demo=QPSK_rc_n0.*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::。t3=0::。% 抽樣判決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))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]。 endenderror=0。dd = 2*bit_in+1。ddd=[dd39。]。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。 endendp=error/1000。figure(1)subplot(2,1,1)。plot(t2,ddd2)。axis([0 100 2 2])。title(39。原序列39。)。subplot(2,1,2)。plot(t2,data_recover_a)。axis([0 100 2 2])。title(39。解調后序列39。)。 仿真結果2：% 設定 T=1, 不加噪聲clear allclose all% 調制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_I39。,20,1)。data_Q1=repmat(data_Q39。,20,1)。for i=1:1e4 data_I2(i)=data_I1(i)。 data_Q2(i)=data_Q1(i)。end。t=0::。f=0::1。xrc=+*cos(pi*f)。data_I2_rc=conv(data_I2,xrc)/。data_Q2_rc=conv(data_Q2,xrc)/。f1=1。t1=0::1e3+。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)。% 解調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::。t3=0::。data_recover=[]。for i=1:20:10000 data_recover=[data_recover I(i:1:i+19) Q(i:1:i+19)]。end。ddd = 2*bit_in+1。ddd1=repmat(ddd39。,10,1)。for i=1:1e4 ddd2(i)=ddd1(i)。endfigure(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 6 6])。 % QPSK誤碼率分析SNRindB1=0:2:10。SNRindB2=0::10。for i=1:length(SNRindB1) [pb,ps]=cm_sm32(SNRindB1(i))。 smld_bit_err_prb(i)=pb。 smld_symbol_err_prb(i)=ps。end。for i=1:length(SNRindB2) SNR=exp(SNRindB2(i)*log(10)/10)。 theo_err_prb(i)=Qfunct(sqrt(2*SNR))。end。title(39。QPSK誤碼率分析39。)。semilogy(SNRindB1,smld_bit_err_prb,39。*39。)。axis([0 10 10e8 1])。hold on。% semilogy(SNRindB1,smld_symbol_err_prb,39。o39。)。仿真結果3：% semilogy(SNRindB1,smld_symbol_err_prb,39。o39。)。semilogy(SNRindB2,theo_err_prb)。legend(39。仿真比特誤碼率39。,39。理論比特誤碼率39。)。hold off。function[y]=Qfunct(x)y=(1/2)*erfc(x/sqrt(2))。function[pb,ps]=cm_sm32(SNRindB)N=10000。E=1。SNR=10^(SNRindB/10)。sgma=sqrt(E/SNR)/2。s00=[1 0]。s01=[0 1]。s11=[1 0]。s10=[0 1]。for i=1:N temp=rand。 if (temp) dsource1(i)=0。 dsource2(i)=0。 elseif (temp) dsource1(i)=0。 dsource2(i)=1。 elseif (temp) dsource1(i)=1。 dsource2(i)=0。 else dsource1(i)=1。 dsource2(i)=1。 end。end。numofsymbolerror=0。numofbiterror=0。for i=1:N n=sgma*randn(size(s00))。 if((dsource1(i)==0)amp。(dsource2(i)==0)) r=s00+n。 elseif((dsource1(i)==0)amp。(dsource2(i)==1)) r=s01+n。 elseif((dsource1(i)==1)amp。(dsource2(i)==0)) r=s10+n。 else r=s11+n。 end。 c00=dot(r,s00)。 c01=dot(r,s01)。 c10=dot(r,s10)。 c11=dot(r,s11)。 c_max=max([c00 c01 c10 c11])。 if (c00==c_max) decis1=0。decis2=0。 elseif(c01==c_max) decis1=0。decis2=1。 elseif(c10==c_max) decis1=1。decis2=0。 else decis1=1。decis2=1。 end。 symbolerror=0。 if(decis1~=dsource1(i)) numofbiterror=numofbiterror+1。 symbolerror=1。 end。 if(decis2~=dsource2(i)) numofbiterror=numofbiterror+1。 symbolerror=1。 end。 if(symbolerror==1) numofsymbolerror=numofsymbolerror+1。 end。end。ps=numofsymbolerror/N。pb=numofbiterror/(2*N)。仿真結果4：(16QAM）正交調制及相干解調原理框圖正交調制原理框圖相干解調原理框圖16QAM調制與解調基于MATLAB的仿真：仿真程序：clc。close all。N=60。 %設定碼元數(shù)量 fb=1。