【正文】 k:39。 end。 xn=xn39。 I1=find(yn= amp。 yn=。 end。 n=(.6:1:N)*m。 [Q,Qn]=two2four(Q,4*m)。解調(diào)信號(hào)頻譜39。 Pe(j)=numoferr/N。)。 y=20*log10(y)。 dt=1/fs。然后經(jīng)過(guò)低通濾波器，低通濾波器LPF濾除乘法器產(chǎn)生的高頻分量，獲得有用信號(hào)，低通濾波器LPF輸出經(jīng)抽樣判決可恢復(fù)出電平信號(hào)。subplot(2,1,2)。axis([0 100 2 2])。title(39。for i=1:20:20000 if sum(data_recover(i:i+19))0。QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc)。data_I1=repmat(data_I39。 %坐標(biāo)范圍grid on。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%一下為顯示波形subplot(325)。%add gain so=somean(so)。 grid on。 z=exp(w*j)。 plot(t,pskn)。 x(x0)=1。信號(hào)經(jīng)過(guò)信道，受到信道噪聲干擾，假設(shè)噪聲為高斯白噪聲。PSK信號(hào)39。 x(x0)=1。axis([0,5,1,2])。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽樣判決if st11(1,m*500+250)。 %與載波1相乘[f,sf1] = t2f(t,st11)。 %已調(diào)信號(hào)頻譜f=(0:length(Ffsk)1)*fm/length(Ffsk)/2。)?；鶐盘?hào)39。 end endendst2=t。圖中兩個(gè)帶通濾波器的作用同于包絡(luò)檢波法，起到分路作用。已調(diào)信號(hào)與載波信號(hào)相乘波形39。解調(diào)輸出波形39。 %一個(gè)數(shù)字信號(hào)‘1’dsig=kron(m,T)。 %頻率比即一個(gè)碼元寬度中的正弦周期個(gè)數(shù),為適配下面濾波器參數(shù)選%取M=500。)。 %時(shí)間Tmin=min(tau)。 %載波的幅度Q=f/F。　MATLAB是由美國(guó)mathworks公司發(fā)布的主要面對(duì)科學(xué)計(jì)算、可視化以及交互式程序設(shè)計(jì)的高科技計(jì)算環(huán)境。它將數(shù)值分析、矩陣計(jì)算、科學(xué)數(shù)據(jù)可視化以及非線性動(dòng)態(tài)系統(tǒng)的建模和仿真等諸多強(qiáng)大功能集成在一個(gè)易于使用的視窗環(huán)境中，為科學(xué)研究、工程設(shè)計(jì)以及必須進(jìn)行有效數(shù)值計(jì)算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計(jì)語(yǔ)言（如C、Fortran）的編輯模式，代表了當(dāng)今國(guó)際科學(xué)計(jì)算軟件的先進(jìn)水平。 %頻率比，即一個(gè)碼元寬度中的正弦周期個(gè)數(shù),為適配下面濾波器參數(shù)%%選%%取，推薦將Q 設(shè)為Q=1/3M=500。 %最小時(shí)刻Tmax=max(tau)。grid on axis([Tmin Tmax *A *A]) %設(shè)置坐標(biāo)范圍運(yùn)行結(jié)果：信號(hào)經(jīng)過(guò)信道后，會(huì)受到噪聲的干擾，假設(shè)信道噪聲為高斯白噪聲。 %一個(gè)正弦周期內(nèi)的采樣點(diǎn)數(shù)t=(0:M1)/M/f。 %數(shù)字信號(hào)波形 Fask=fft(ask)。) %畫(huà)出濾波后的波形grid on %添加網(wǎng)格axis([Tmin Tmax ]) %設(shè)置坐標(biāo)范圍subplot(323)。)%已調(diào)信號(hào)與%載波信號(hào)相乘波形grid on axis([Tmin Tmax *A *A]) %設(shè)置坐標(biāo)范圍運(yùn)行結(jié)果：分析結(jié)果：解調(diào)信號(hào)的頻譜沒(méi)出來(lái)，檢查代碼終究還是沒(méi)找到問(wèn)題所在。它們的輸出分別與相應(yīng)的同步相干載波相乘，再分別經(jīng)低通濾波器濾掉二倍頻信號(hào)，取出含基帶數(shù)字信息的低頻信號(hào)，抽樣判決器在抽樣脈沖到來(lái)時(shí)對(duì)兩個(gè)低頻信號(hào)的抽樣值vv2進(jìn)行比較判決（判決規(guī)則同于包絡(luò)檢波法），即可還原出基帶信號(hào)。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%基帶信號(hào)求反for n=1:j。)。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%調(diào)制F1=st1.*s1。 %已調(diào)信號(hào)頻譜帶寬subplot(427)。 %通過(guò)低通濾波器[t,st11] = lpf(f,sf1,2*fm)。for j=m*500+1:(m+1)*500。運(yùn)行結(jié)果：調(diào)制輸出：加噪后解調(diào)輸出：2PSK調(diào)制解調(diào)原理與系統(tǒng)仿真、2PSK調(diào)制原理與仿真 當(dāng)相移常數(shù)Kp=π時(shí)， 當(dāng)基帶數(shù)字信號(hào)采用幅度為1寬度為TS的矩形脈沖的雙極性非歸零碼表示時(shí)，時(shí)域表示式為 2PSK信號(hào)波形 2PSK信號(hào)的調(diào)制原理框圖2PSK調(diào)制代碼：clear all。 subplot(311)。)。解調(diào)代碼：clear all。car=sin(2*pi*fc*t)。 axis([0 200e6 2 2])。 BW=2*pi*fBW/fs。 cm=fpskn.*car。%removes DC ponent subplot(324)。 plot (t,Vo), title(39。xlabel(39。,20,1)。QPSK_rc_n0=QPSK_rc+n0。 data_recover_a(i:i+19)=1。4psk信號(hào)波形39。title(39。plot(t2,data_recover_a)。16QAM正交相干解調(diào)如圖2所示。 t=0:dt:Tdt。 f1=m/n。 SNR_in_dB=8:2:24。 %統(tǒng)計(jì)誤碼率 end。)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%一下為需要用到的M文件%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%產(chǎn)生二進(jìn)制信源隨機(jī)序列 function [info]=random_binary(N) if nargin == 0, %如果沒(méi)有輸入?yún)?shù)，則指定信息序列為10000個(gè)碼元 N=10000。 if Kbase==2。 n=fix(n)。 if nargin4。 y=yn。 yn)。 xn=xn(:)。%如果無(wú)輸出，則作圖 if nargout1。,[0 0],[ ],39。 plot([ ],[0 0],39。 i1=i1+m。 end。 yn(I0)=zeros(size(I0))。 end。 alfa=。 N=length(y)/m。 Q=x(2:2:nn)。title(39。 end。16QAM調(diào)制方式對(duì)應(yīng)的星座圖39。 y(q)=1e04。nn=length(info)。 、16QAM解調(diào)原理16QAM信號(hào)采取正交相干解調(diào)的方法解調(diào)，解調(diào)器首先對(duì)收到的16QAM信號(hào)進(jìn)行正交相干解調(diào)，一路與tcωcos相乘，一路與tcωsin相乘。)。plot(n0)。axis([0 100 2 2])。bit_recover=[]。Q_rc=data_Q2_rc.*sin(2*pi*f1*t1)。data_Q=2*bit_Q+1。 axis([0 2e4 2 2])。 end end Vo=Vs。 so=so*10。)。 w=2*pi*f/fs。%已調(diào)信號(hào)加噪聲后的頻譜%已調(diào)信號(hào)加噪聲信號(hào)后的頻譜帶寬f1=(0:length(Fpskn)1)*fs/length(Fpskn)/2 subplot(212)。 x(x0)=am。2PSK信號(hào)經(jīng)過(guò)帶通濾波器得到有用信號(hào)，經(jīng)相乘器與本地載波相乘再經(jīng)過(guò)低通濾波器得到低頻信號(hào)v(t)，再經(jīng)抽樣判決得到基帶信號(hào)。 title(39。 x(x0)=am。)。for m=0:i1。)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解調(diào)st11=fsk1.*s1。 %對(duì)信號(hào)進(jìn)行傅里葉變換Ffsk=abs(Ffsk(1:length(Ffsk)/2+1))。載波信號(hào)239。title(39。 end else for m=j/i*(n1)+1:j/i*n st1(m)=1。相干檢測(cè)法原理：相干檢測(cè)的具體解調(diào)電路是同步檢波器，原理框圖如下圖所示。) %畫(huà)出原始%%信號(hào)波形與%解調(diào)后的信號(hào)作對(duì)比axis([Tmin Tmax ]) %設(shè)置坐標(biāo)范圍grid on %添加網(wǎng)格subplot(321),plot(tau,sig_mul) title(39。 %解調(diào)信號(hào)頻譜帶寬subplot(3,2,2) %子圖分割plot(tau,x),title(39。 %最大時(shí)刻T=ones(size(carry1))。 %載波的幅度Q=f/F。調(diào)制后加噪聲波形39。 %長(zhǎng)度tau=(0:N1)/(M1)/f。 %正弦載波信號(hào)的頻率A=1。、MATLAB仿真平臺(tái)簡(jiǎn)介MATLAB是矩陣實(shí)驗(yàn)室（Matrix Laboratory）的簡(jiǎn)稱，是美國(guó)MathWorks公司出品的商業(yè)數(shù)學(xué)軟件，用于算法開(kāi)發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計(jì)算的高級(jí)技術(shù)計(jì)算語(yǔ)言和交互式環(huán)境，主要包括MATLAB和Simulink兩大部分?！?MATLAB