【正文】 仿真2ASK有兩種基本解調(diào)方法：相干解調(diào)法（同步檢測(cè)法）和非相fc干解調(diào)法（包絡(luò)檢波法）。figure(2)subplot(111),plot(tau,ask1) %畫出調(diào)制加噪波形title(39。) %子圖分割grid on axis([Tmin Tmax *A *A]) %設(shè)置坐標(biāo)范圍subplot(224)。) %畫出數(shù)字信號(hào)的波形圖grid on axis([Tmin Tmax ]) %設(shè)置坐標(biāo)范圍subplot(2,2,2) %子圖分割plot(tau,carry2),title (39。 %最大時(shí)刻T=ones(size(carry1))。 %已調(diào)信號(hào)頻譜帶寬N=length(ask)。,0,V)。 %載波信號(hào)ask=kron(m,carry1)。 %一個(gè)正弦周期內(nèi)的采樣點(diǎn)數(shù)t=(0:M1)/M/f。 %數(shù)字信號(hào)的帶寬f=600。相應(yīng)的調(diào)制如圖23：乘法器圖23模擬相乘法調(diào)制代碼：clear。2ASK是利用代表數(shù)字信息“0”或“1”的基帶矩形脈沖去鍵控一個(gè)連續(xù)的載波，使載波時(shí)斷時(shí)續(xù)地輸出?！?MATLAB的基本數(shù)據(jù)單位是矩陣，它的指令表達(dá)式與數(shù)學(xué)、工程中常用的形式十分相似，故用MATLAB來解算問題要比用C，F(xiàn)ORTRAN等語言完成相同的事情簡(jiǎn)捷得多，并且MATLAB也吸收了像Maple等軟件的優(yōu)點(diǎn)，使MATLAB成為一個(gè)強(qiáng)大的數(shù)學(xué)軟件。任務(wù)和要求、課程設(shè)計(jì)目的?通過課程設(shè)計(jì)加深對(duì)《通信原理》課程基本內(nèi)容的理解?學(xué)習(xí)與掌握通信系統(tǒng)仿真的基本方法、課程設(shè)計(jì)要求1）完成2ASK、2FSK、2PSK、4PSK、16QAM五種調(diào)制系統(tǒng)的仿真，要求有調(diào)制、信道、解調(diào)三個(gè)部分；2）信道為高斯信道；3）系統(tǒng)參數(shù)可自行設(shè)置；4）仿真結(jié)果要求給出已調(diào)信號(hào)時(shí)域波形和頻譜，解調(diào)后信號(hào)的時(shí)域波形與頻譜。為了使數(shù)字信號(hào)在帶通信道中傳輸，必須使用數(shù)字基帶信號(hào)對(duì)載波進(jìn)行調(diào)制，以使信號(hào)與信道的特性相匹配。這種用數(shù)字基帶信號(hào)控制載波，把數(shù)字基帶信號(hào)變換為數(shù)字帶通信號(hào)的過程稱為數(shù)字調(diào)制。、MATLAB仿真平臺(tái)簡(jiǎn)介MATLAB是矩陣實(shí)驗(yàn)室（Matrix Laboratory）的簡(jiǎn)稱，是美國(guó)MathWorks公司出品的商業(yè)數(shù)學(xué)軟件，用于算法開發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計(jì)算的高級(jí)技術(shù)計(jì)算語言和交互式環(huán)境，主要包括MATLAB和Simulink兩大部分。在新的版本中也加入了對(duì)C，F(xiàn)ORTRAN，C++，JAVA的支持。有載波輸出時(shí)表示發(fā)送“1”，無載波輸出時(shí)表示發(fā)送“0”。 %清除空間變量m=[1 0 1 1 0 1 1 0]。 %正弦載波信號(hào)的頻率A=1。 %一個(gè)正弦信號(hào)周期內(nèi)的時(shí)間carry1=repmat(A*sin(2*pi*f*t),1,Q)。 %調(diào)制后的信號(hào)V=。 %加高斯白噪聲后的已調(diào)信號(hào)Fask=fft(ask)。 %長(zhǎng)度tau=(0:N1)/(M1)/f。 %一個(gè)數(shù)字信號(hào)‘1’dsig=kron(m,T)。載波波形39。plot(f1,Fask),title(39。調(diào)制后加噪聲波形39。相干解調(diào)需要將載頻位置的已調(diào)信號(hào)頻譜重新搬回原始基帶位置，因此用相乘器與載波相乘來實(shí)現(xiàn)。經(jīng)過各個(gè)模塊后波形變化如圖7所示。 %數(shù)字信號(hào)序列Lm=length(m)。 %載波的幅度Q=f/F。 %一個(gè)碼元寬度內(nèi)的正弦載波信號(hào)Lcarry1=length(carry1)。ask1=imnoise(ask,39。 %已調(diào)信號(hào)與載波信號(hào)相乘N=length(ask)。 %最大時(shí)刻T=ones(size(carry1))。 %已調(diào)信號(hào)頻譜f1=(0:length(Fask)1)*f/length(Fask)/2。%由原型濾波器向?qū)嶋H濾波器轉(zhuǎn)換，獲得濾波器的分子，分母系數(shù)[num,den]=butter(Ord,omega_c,39。 %獲得濾波器傳遞函數(shù)%濾波x=lsim(h,sig_mul,tau)。 %解調(diào)信號(hào)頻譜帶寬subplot(3,2,2) %子圖分割plot(tau,x),title(39。解調(diào)信號(hào)頻譜39。 %抽樣判決點(diǎn)的選取y=(x(t_judge))39。 %抽樣判決后的信號(hào)對(duì)應(yīng)的時(shí)間subplot(3,2,4) %子圖分割plot(n_tau,y_value), title(39。) %畫出原始%%信號(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。2FSK調(diào)制解調(diào)原理與系統(tǒng)仿真、2FSK調(diào)制原理 2FSK是利用數(shù)字基帶信號(hào)控制在波的頻率來傳送信息。如圖11（a）所示。 (a) (b) 2FSK信號(hào)產(chǎn)生原理圖 由鍵控法產(chǎn)生原理可知，一位相位離散的2FSK信號(hào)可看成不同頻率交替發(fā)送的兩個(gè)2ASK信號(hào)之和，即 其中是脈寬為的矩形脈沖表示的NRZ數(shù)字基帶信號(hào)。相干檢測(cè)法原理：相干檢測(cè)的具體解調(diào)電路是同步檢波器，原理框圖如下圖所示。2FSK調(diào)制解調(diào)代碼：clear allclose alli=10。f1=10。 %載波1帶寬B2=2*f2。 end else for m=j/i*(n1)+1:j/i*n st1(m)=1。 st2(n)=0。 %加入載波1F2=st2.*s2。 %fsk信號(hào)加高斯白噪聲Ffsk1=fft(fsk1)。title(39。subplot(422)。)。載波信號(hào)139。載波信號(hào)239。 %加入載波2subplot(425)。)。s2*st239。 %對(duì)信號(hào)進(jìn)行傅里葉變換Ffsk=abs(Ffsk(1:length(Ffsk)/2+1))。title(39。)。title(39。)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解調(diào)st11=fsk1.*s1。plot(t,st11)。st22=fsk1.*s2。plot(t,st22)。for m=0:i1。endelsefor j=m*500+1:(m+1)*500。plot(t,at)。)subplot(326)。)。 %關(guān)閉窗口fs=8e5。 fc=2e5。 A=1。 x(x0)=am。 axis([0 2e4 2 2])。 grid on %畫網(wǎng)格car=sin(2*pi*fc*t)。 %已調(diào)信號(hào)頻譜f=(0:length(Fpsk)1)*fs/length(Fpsk)/2。 title(39。 subplot(313)。)。 2PSK信號(hào)屬于DSB信號(hào)，它的解調(diào)，不再能采用包絡(luò)檢測(cè)的方法，只能進(jìn)行相干解調(diào)。2PSK信號(hào)經(jīng)過帶通濾波器得到有用信號(hào)，經(jīng)相乘器與本地載波相乘再經(jīng)過低通濾波器得到低頻信號(hào)v(t)，再經(jīng)抽樣判決得到基帶信號(hào)。 fs=8e5。 fc=2e5。 A=1。 x(x0)=am。 %載波調(diào)制vn=。 grid on。 axis([0 .2e3 1 1])。%已調(diào)信號(hào)加噪聲后的頻譜%已調(diào)信號(hào)加噪聲信號(hào)后的頻譜帶寬f1=(0:length(Fpskn)1)*fs/length(Fpskn)/2 subplot(212)。加噪后信號(hào)39。subplot(321)。)。 w=2*pi*f/fs。%BW=2(1a)/sqrt(a) p=(j^2*a^2)。%avoid log(0) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%濾波器系數(shù) a=[1 0 ]。 plot(t,fpskn)。)。 plot(t,cm)。通過相乘器后輸出39。%gain=(1p)/2 Hz1=gain1*(z+1)./(z(p))。 so=so*10。 axis([0 6e4 ])。 grid on。 %設(shè)立比較標(biāo)準(zhǔn) error=0。 end end Vo=Vs。), axis([0 2e4 5 5]) grid on。 %坐標(biāo)軸標(biāo)記 ylabel(39。 plot (t,x)。 axis([0 2e4 2 2])。)。運(yùn)行結(jié)果： 分析可知：解調(diào)后恢復(fù)原始信號(hào)。%%%%%%%%%