【正文】 ，通常位于每個(gè)碼元的中央位置，其重復(fù)周期等于碼元的寬度。經(jīng)抽樣、判決后將碼元再生，即可恢復(fù)出數(shù)字序列。相干檢測就是同步解調(diào)，要求接收機(jī)產(chǎn)生一個(gè)與發(fā)送載波同頻同相的本地載波信號，稱其為同步載波或相干載波。低通濾波器的截止頻率與基帶數(shù)字信號的最高頻率相等。雖然2ASK信號中確實(shí)存在著載波分量，原則上可以通過窄帶濾波器或鎖相環(huán)來提取同步載波，但這會給接收設(shè)備增加復(fù)雜性。2ASK調(diào)制與解調(diào)基于MATLAB仿真：2ASK基于MATLAB的程序代碼：function ask1()sim(39。)。y1=ScopeData5(:,2)。y2=ScopeData2(:,2)。y3=ScopeData(:,2)。y5=ScopeData1(:,2)。y6=ScopeData7(:,2)。y7=ScopeData3(:,2)。plot(t1,y1),grid。濾波后的波形39。plot(t2,y2),grid。相干解調(diào)后的波形39。stairs(t3,y3),axis([0 20 ]),grid。 ASK調(diào)制與解調(diào)各環(huán)節(jié)波形39。源信號39。plot(t5,y5),grid。調(diào)制后的波形39。plot(t6,y6),grid。經(jīng)過高斯信道后的波形39。stairs(t7,y7),axis([0 20 ]),grid。抽樣判決后的波形39。f1=(n1/2:1:n1/21)。subplot(3,2,1)。源信號頻譜39。f2=(n2/2:1:n2/21)。subplot(3,2,2)。調(diào)制后的頻譜39。subplot(3,2,3)。相干解調(diào)后的頻譜39。subplot(3,2,4)。濾波后的頻譜39。subplot(3,2,5)。抽樣判決后的頻譜39。二進(jìn)制相移鍵控中，通常用相位0和來分別表示“0”或“1”。 （822）因此，在某一個(gè)碼元持續(xù)時(shí)間內(nèi)觀察時(shí)，有，或（823）當(dāng)碼元寬度為載波周期的整數(shù)倍時(shí)，2PSK信號的典型波形如圖811所示圖8112PSK信號的典型波形2PSK信號的調(diào)制方框圖如圖812示。圖8122PSK調(diào)制器框圖就模擬調(diào)制法而言，與產(chǎn)生2ASK信號的方法比較，只是對要求不同，因此2PSK信號可以看作是雙極性基帶信號作用下的DSB調(diào)幅信號。2PSK信號屬于DSB信號，它的解調(diào)，不再能采用包絡(luò)檢測的方法，只能進(jìn)行相干解調(diào)，其方框圖如圖813。圖813 2PSK信號接收系統(tǒng)方框圖不考慮噪聲時(shí)，帶通濾波器輸出可表示為 （824）式中為2PSK信號某一碼元的初相。與同步載波相乘后，輸出為 （826）根據(jù)發(fā)端產(chǎn)生2PSK信號時(shí)（0或）代表數(shù)字信息（“1”或“0”）的規(guī)定，以及收端與的關(guān)系的特性，抽樣判決器的判決準(zhǔn)則為 （828）其中為在抽樣時(shí)刻的值。圖814 2PSK解調(diào)各點(diǎn)波形可見，2PSK信號相干解調(diào)的過程實(shí)際上是輸入已調(diào)信號與本地載波信號進(jìn)行極性比較的過程，故常稱為極性比較法解調(diào)。如果同步載波的相位發(fā)生變化，如0相位變?yōu)橄辔换蛳辔蛔優(yōu)?相位，則恢復(fù)的數(shù)字信息就會發(fā)生“0”變“1”或“1”變“0”，從而造成錯(cuò)誤的恢復(fù)。絕對移相的主要缺點(diǎn)是容易產(chǎn)生相位模糊，造成反向工作。 close all。%抽樣頻率 fm=20e3。 final=(1/fs)*(n1)。 % 載波頻率 t=0:1/fs:(final)。%耐奎斯特頻率 %用正弦波產(chǎn)生方波 twopi_fc_t=2*pi*fm*t。 phi=0。 % 方波 am=1。 x(x0)=1。 plot(t,x)。 title(39。)。%載波 ask=x.*car。 plot(t,ask)。 title(39。)。 vn=。%產(chǎn)生噪音 subplot(323)。 grid on。噪音信號39。 axis([0 .2e3 1 1])。%調(diào)制后加噪 subplot(324)。 axis([0 200e6 2 2])。加噪后信號39。 grid on。 f=[0:3e3:4e5]。 z=exp(w*j)。 a=.8547。 gain=.135。 subplot(325)。 title(39。)。 Hz(Hz==0)=10^(8)。 plot(f,20*log10(abs(Hz)))。 title(39。)。 %濾波器系數(shù) a=[1 0 ]。%gain*[1 0 1] faskn=filter(b,a,askn)。 plot(t,faskn)。 title(39。)。 cm=faskn.*car。 plot(t,cm)。 grid on。通過相乘器后輸出39。 %低通濾波器 p=。%gain=(1p)/2 Hz1=gain1*(z+1)./(z(p))。 Hz1(Hz1==0)=10^(8)。 grid on。LPF 3dB response39。 axis([0 5e4 3 1])。%(z(p)) b1=[ ]。 so=so*10。%removes DC ponent subplot(324)。 axis([0 8e4 ])。通過低通濾波器后輸出39。 grid on。 Low=。%設(shè)立比較標(biāo)準(zhǔn) error=0。 for ii=1:len1 if so(ii) = vt Vs(ii)=High。 end end Vo=Vs。 plot (t,Vo), title(39。), axis([0 2e4 5 5]) grid on。時(shí)間 (s)39。幅度(V)39。下面著重介紹多進(jìn)制數(shù)字相位調(diào)制的這兩種形式。由于每一種載波相位代表兩個(gè)比特信息，故每個(gè)四進(jìn)制碼元又被稱為雙比特碼元，習(xí)慣上把雙比特的前一位用代表，后一位用代表。1）相位選擇法由式（854）可以看出，在一個(gè)碼元持續(xù)時(shí)間內(nèi)，4PSK信號為載波四個(gè)相位中的某一個(gè)。圖中，四相載波發(fā)生器產(chǎn)生4PSK信號所需的四種不同相位的載波。按照輸入的雙比特碼元的不同，邏輯選相電路輸出相應(yīng)相位的載波。圖827相位選擇法產(chǎn)生4PSK信號（B方式）方框圖圖中產(chǎn)生的是B方式的4PSK信號。2）直接調(diào)相法由式（856）可以看出，4PSK信號也可以采用正交調(diào)制的方式產(chǎn)生。它可以看成是由兩個(gè)載波正交的2PSK調(diào)制器構(gòu)成，分別形成圖828（b）中的虛線矢量，再經(jīng)加法器合成后，得圖（b）中實(shí)線矢量圖。圖828直接調(diào)相法產(chǎn)生4PSK信號方框圖若要產(chǎn)生4PSK的A方式波形，只需適當(dāng)改變振蕩載波相位就可實(shí)現(xiàn)。圖829是B方式4PSK信號相干解調(diào)器的組成方框圖。此法也稱為極性比較法。在2PSK信號相干解調(diào)過程中會產(chǎn)生“倒”即“180176。同樣，對于4PSK信號相干解調(diào)也會產(chǎn)生相位模糊問題，并且是0176。、180176。四個(gè)相位模糊。tt=length(t)。for i=1:tt if (t(i)=1 amp。 t(i)=7)。 else x1(i)=1。t2=0::。t4=0::。x2=ones(1,800)。 t1(i)=2) | (t1(i)=4amp。 x2(i)=1。 endendf=0::1。y1=conv(x1,xrc)/。n0=randn(size(t2))。i=x1.*cos(2*pi*f1*t)。I=i(101:800)。QPSK=sqrt(1/2).*I+sqrt(1/2).*Q。n1=randn(size(t2))。q_rc=y2.*sin(2*pi*f1*t4)。Q_rc=q_rc(1:700)。QPSK_rc_n1=QPSK_rc+n1。plot(t3,i_rc)。ylabel(39。)。plot(t4,q_rc)。ylabel(39。)。plot(t2,QPSK_rc)。ylabel(39。)。plot(t2,QPSK_rc_n1)。ylabel(39。) 。data_Q = 2*bit_Q+1。,20,1)。,20,1)。 data_Q2(i)=data_Q1(i)。f=0::1。data_I2_rc=conv(data_I2,xrc)/。f1=1。n0=rand(size(t1))。Q_rc=data_Q2_rc.*sin(2*pi*f1*t1)。QPSK_rc_n0=QPSK_rc+n0。Q_demo=QPSK_rc_n0.*sin(2*pi*f1*t1)。 Q_recover=conv(Q_demo,xrc)。Q=Q_recover(11:10010)。t3=0::。for i=1:20:10000 data_recover=[data_recover I(i:1:i+19) Q(i:1:i+19)]。bit_recover=[]。 bit_recover=[bit_recover 1]。 bit_recover=[bit_recover 1]。dd = 2*bit_in+1。]。for i=1:2e4 ddd2(i)=ddd1(i)。 endendp=error/1000。plot(t2,ddd2)。title(39。)。plot(t2,data_recover_a)。title(39。)。bit_I = bit_in(1:2:1e3)。data_I = 2*bit_I+1。data_I1=repmat(data_I39。data_Q1=repmat(data_Q39。for i=1:1e4 data_I2(i)=data_I1(i)。end。f=0::1。data_I2_rc=conv(data_I2,xrc)/。f1=1。I_rc=data_I2_rc.*cos(2*pi*f1*t1)。QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc)。Q_demo=QPSK_rc.*sin(2*pi*f1*t1)。Q_recover=conv(Q_demo,xrc)。Q=Q_recover(11:10010)。t3=0::。for i=1:20:10000 data_recover=[data_recover I(i:1:i+19) Q(i:1:i+19)]。ddd = 2*bit_in+1。,10,1)。endfigure(1)subplot(4,1,1)。axis([0 20 6 6])。plot(t3,Q)。subplot(4,1,3)。axis([0 20 6 6])。plot(t,ddd2)。 % QPSK誤碼率分析SNRindB1=0:2:10。for i=1:length(SNRindB1) [pb,ps]=cm_sm32(SNRindB1(i))。 smld_symbol_err_prb(i)=ps。for i=1:length(SNRindB2) SNR=exp(SNRindB2(i)*log(10)/10)。end。QPSK誤碼率分析39。semilogy(SNRindB1,smld_bit_err_prb,39。)。hold on。o39。仿真結(jié)果3：% semilogy(SNRindB1,smld_symbol_err_prb,39。)。legend(39。,39。)。function[y]=Qfunct(x)y=(1/2)*erfc(x/sqrt(2))。E=1。sgma=sqrt(E/SNR)/2。s01=[0 1]。s10=[0 1]。 if (temp) dsource1(i)=0。 elseif (temp) dsource1(i)=0。 elseif (temp) dsource1(i)=1。 else dsource1(i)=1。 end。numofsymbolerror=0。for i=1:N n=sgma*randn(size(s00))。(dsource2(i)==0)) r=s00+n。(dsource2(i)==1)) r=s01+n。(dsource2(i)==0)) r=s10+n。 end。 c01=dot(r,s01)。 c11=dot(r,s11)。 if (c00==c_max) decis1=0。 elseif(c01==c_max) decis1=0。 elseif(c10==c_max) decis1=1。 else decis1=1。 end。 if(decis1~=dsource1(i)) numofbiterror=numofbiterror+1。 end。 symbolerror=1。 if(symbolerror==1) numofsymbolerror=numofsymbolerror+1。end。pb=numofbiterror/(2*N)。close all。 %設(shè)定碼元數(shù)量 fb=