【正文】 。N=60。仿真結(jié)果4：(16QAM）正交調(diào)制及相干解調(diào)原理框圖正交調(diào)制原理框圖相干解調(diào)原理框圖16QAM調(diào)制與解調(diào)基于MATLAB的仿真：仿真程序：clc。ps=numofsymbolerror/N。 end。 end。 if(decis2~=dsource2(i)) numofbiterror=numofbiterror+1。 symbolerror=1。 symbolerror=0。decis2=1。decis2=0。decis2=1。decis2=0。 c_max=max([c00 c01 c10 c11])。 c10=dot(r,s10)。 c00=dot(r,s00)。 else r=s11+n。 elseif((dsource1(i)==1)amp。 elseif((dsource1(i)==0)amp。 if((dsource1(i)==0)amp。numofbiterror=0。end。 dsource2(i)=1。 dsource2(i)=0。 dsource2(i)=1。 dsource2(i)=0。for i=1:N temp=rand。s11=[1 0]。s00=[1 0]。SNR=10^(SNRindB/10)。function[pb,ps]=cm_sm32(SNRindB)N=10000。hold off。理論比特誤碼率39。仿真比特誤碼率39。semilogy(SNRindB2,theo_err_prb)。o39。)。% semilogy(SNRindB1,smld_symbol_err_prb,39。axis([0 10 10e8 1])。*39。)。title(39。 theo_err_prb(i)=Qfunct(sqrt(2*SNR))。end。 smld_bit_err_prb(i)=pb。SNRindB2=0::10。axis([0 20 6 6])。subplot(4,1,4)。plot(t2,data_recover)。axis([0 20 6 6])。subplot(4,1,2)。plot(t3,I)。for i=1:1e4 ddd2(i)=ddd1(i)。ddd1=repmat(ddd39。end。data_recover=[]。t2=0::。I=I_recover(11:10010)。I_recover=conv(I_demo,xrc)。% 解調(diào)I_demo=QPSK_rc.*cos(2*pi*f1*t1)。Q_rc=data_Q2_rc.*sin(2*pi*f1*t1)。t1=0::1e3+。data_Q2_rc=conv(data_Q2,xrc)/。xrc=+*cos(pi*f)。t=0::。 data_Q2(i)=data_Q1(i)。,20,1)。,20,1)。data_Q = 2*bit_Q+1。bit_Q = bit_in(2:2:1e3)。 仿真結(jié)果2：% 設(shè)定 T=1, 不加噪聲clear allclose all% 調(diào)制bit_in = randint(1e3, 1, [0 1])。解調(diào)后序列39。axis([0 100 2 2])。subplot(2,1,2)。原序列39。axis([0 100 2 2])。figure(1)subplot(2,1,1)。endfor i=1:1e3 if bit_recover(i)~=ddd(i) error=error+1。ddd1=repmat(ddd,20,1)。ddd=[dd39。 endenderror=0。 else data_recover_a(i:i+19)=1。for i=1:20:20000 if sum(data_recover(i:i+19))0 data_recover_a(i:i+19)=1。end。% 抽樣判決data_recover=[]。t2=0::。I=I_recover(11:10010)。% 低通濾波I_recover=conv(I_demo,xrc)。% 解調(diào)I_demo=QPSK_rc_n0.*cos(2*pi*f1*t1)。QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc)。I_rc=data_I2_rc.*cos(2*pi*f1*t1)。t1=0::1e3+。data_Q2_rc=conv(data_Q2,xrc)/。xrc=+*cos(pi*f)。end。for i=1:1e4 data_I2(i)=data_I1(i)。data_Q1=repmat(data_Q39。data_I1=repmat(data_I39。仿真結(jié)果1：data_I = 2*bit_I+1。加入噪聲39。axis([1 8 1 1])。subplot(4,1,4)。合成序列39。axis([1 8 1 1])。subplot(4,1,3)。b序列39。axis([1 8 1 1])。subplot(4,1,2)。a序列39。axis([1 8 1 1])。figure(1)subplot(4,1,1)。QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc)。I_rc=i_rc(101:800)。i_rc=y1.*cos(2*pi*f1*t3)。QPSK_n=(sqrt(1/2).*I+sqrt(1/2).*Q)+n0。Q=q(1:700)。q=x2.*sin(2*pi*f1*t1)。f1=1。y2=conv(x2,xrc)/。xrc=+*cos(pi*f)。 else x2(i)=1。 t1(i)=8)。for i=1:tt1 if (t1(i)=0 amp。tt1=length(t1)。t3=1::。 endendt1=[0::]。 x1(i)=1。 t(i)=1) | (t(i)=5amp。x1=ones(1,800)。4PSK基于MATLAB的仿真% 調(diào)相法clear allclose allt=[1::]。和270176。、90176。相位模糊”現(xiàn)象。4PSK信號(hào)的相干解調(diào)若解調(diào)4PSK信號(hào)（A方式），只需適當(dāng)改變相移網(wǎng)絡(luò)。圖829圖中兩個(gè)相互正交的相干載波分別檢測(cè)出兩個(gè)分量和，然后，經(jīng)并/串變換器還原成二進(jìn)制雙比特串行數(shù)字信號(hào)，從而實(shí)現(xiàn)二進(jìn)制信息恢復(fù)。（2）4PSK信號(hào)的解調(diào)由于4PSK信號(hào)可以看作是兩個(gè)載波正交的2PSK信號(hào)的合成，因此，對(duì)4PSK信號(hào)的解調(diào)可以采用與2PSK信號(hào)類似的解調(diào)方法進(jìn)行。顯然其為B方式4PSK相位配置情況。B方式4PSK時(shí)的原理方框圖如圖828（a）所示。要想形成A方式的4PSK信號(hào)，只需調(diào)整四相載波發(fā)生器輸出的載波相位即可。例如，B方式情況下，雙比特碼元為11時(shí)，輸出相位為45?的載波；雙比特碼元為01時(shí)，輸出相位為135?的載波等。輸入的二進(jìn)制數(shù)碼經(jīng)串/并變換器輸出雙比特碼元。因此，可以用相位選擇法產(chǎn)生4PSK信號(hào)，其原理如圖827所示。（1）4PSK信號(hào)的產(chǎn)生多相制信號(hào)常用的產(chǎn)生方法有：直接調(diào)相法及相位選擇法。4PSK利用載波的四種不同相位來(lái)表征數(shù)字信息。)仿真結(jié)果： (4PSK）4PSK信號(hào)的產(chǎn)生與解調(diào)在進(jìn)制數(shù)字相位調(diào)制中，四進(jìn)制絕對(duì)移相鍵控（4PSK，又稱QPSK）和四進(jìn)制差分相位鍵控（4DPSK，又稱QDPSK）用的最為廣泛。), ylabel(39。 xlabel(39。解調(diào)后輸出信號(hào)39。 subplot(325)。 else Vs(ii)=Low。 len1=length(so)。 vt=0。 %Comparator High=。)。 title(39。 plot(t,so)。%add gain so=somean(so)。%gain*[1 1] so=filter(b1,a1,cm)。 %濾波器系數(shù) a1=[1 ]。)。 title(39。%avoid log(0) plot(f,20*log10(abs(Hz1)))。 subplot(323)。 gain1=。)。 title(39。 axis([0 100e6 2 2])。%解調(diào) subplot(322)。 grid on。通過(guò)帶通濾波后輸出39。 axis([0 100e6 2 2])。 figure(2) subplot(321)。%[1 0 p] b=[ 0 ]。 axis([1e5 3e5 3 1])。Receiver 3dB Filter Response39。 grid on。%avoid log(0) subplot(326)。 grid on。帶通濾波器39。 plot(f,abs(Hz))。 Hz=gain*(z+1).*(z1)./(z.^2(p))。%BW=2(1a)/sqrt(a) p=(j^2*a^2)。 BW=2*pi*fBW/fs。 w=2*pi*f/fs。 %帶通濾波 fBW=40e3。)。 title(39。 plot(t,askn)。 askn=(ask+noise)。)。 title(39。 plot(t,noise)。 noise=vn*(randn(size(t)))。 grid on。PSK信號(hào)39。 axis([0 200e6 2 2])。%載波調(diào)制 subplot(322)。 grid on car=sin(2*pi*fc*t)?；鶐盘?hào)39。 axis([0 2e4 2 2])。 figure(1) subplot(321)。 x(x0)=am。 x = A * cos(twopi_fc_t + phi)。 A=1。 Fn=fs/2。 fc=2e5。%基帶頻率 n=2*(6*fs/fm)。 fs=8e5。2PSK基于MATLAB的仿真：2PSK基于MATLAB的程序代碼：clear all。這種因?yàn)楸镜貐⒖驾d波倒相，而在接收端發(fā)生錯(cuò)誤恢復(fù)的現(xiàn)象稱為“倒”現(xiàn)象或“反向工作”現(xiàn)象。由于2PSK信號(hào)實(shí)際上是以一個(gè)固定初相的末調(diào)載波為參考的，因此，解調(diào)時(shí)必須有與此同頻同相的同步載波。2PSK接收系統(tǒng)各點(diǎn)波形如圖814所示。（825）經(jīng)低通濾波器濾除高頻分量，得解調(diào)器輸出為時(shí)，代表數(shù)字“0”；時(shí)，代表數(shù)字“1”。 工作原理簡(jiǎn)要分析如下。而就鍵控法來(lái)說(shuō)，用數(shù)字基帶信號(hào)控制開關(guān)電路，選擇不同相位的載波輸出，這時(shí)為單極性NRZ或雙極性NRZ脈沖序列信號(hào)均可。圖（a）是產(chǎn)生2PSK信號(hào)的模擬調(diào)制法框圖；圖（b）是產(chǎn)生2PSK信號(hào)的鍵控法框圖。 （821）式中，是高度為1，寬度為的門函數(shù)； （820）這里，與2ASK及2FSK時(shí)不同，為雙極性數(shù)字基帶信號(hào)，即2PSK已調(diào)信號(hào)的時(shí)域表達(dá)式為)仿真結(jié)果： (2PSK）2PSK原理2PSK的一般原理及實(shí)現(xiàn)方法絕對(duì)相移是利用載波的相位（指初相）直接表示數(shù)字信號(hào)的相移方式。plot(f2,abs(fftshift(s))),%axis([400 400 0 300])xlabel(39。)s=fft(y7,n2)。plot(f2,abs(fftshift(s))),%axis([400 400 0 300])xlabel(39。)s=fft(y1,n2)。plot(f2,abs(fftshift(s))),%axis([400 400 0 300])xlabel(39。)s=fft(y2,n2)。plot(f2,abs(fftshift(s))),%axis([400 400 0 300])xlabel(39。s=fft(y5,n2)。)n2=1024。plot(f1,abs(fftshift(s))),%axis([30 30 0 7])xlabel(39。s=fft(y3,n1)。)figure(2)n1=1024。xlabel(39。)subplot(3,2,6)。xlabel(39。)subplot(3,2,3)。xlabel(39。)subplot(3,2,2)。)xlabel(39。title(39。)subplot(3,2,1)。xlabel(39。)subplot(3,2,4)。xlabel(39。 %引入示波器的數(shù)據(jù)subplot(3,2,5)。t7=ScopeData3(:,1)。t6=ScopeData7(:,1)。t5=ScopeData1(:,1)。t3=ScopeData(:,1)。t2=ScopeData2(:,1)。t1=ScopeData5(:,1)。dm39。因此，實(shí)際中很少采