【正文】 %已調(diào)信號(hào)與載波信號(hào)相乘N=length(ask)。 %長度tau=(0:N1)/(M1)/f。 %時(shí)間Tmin=min(tau)。 %最小時(shí)刻Tmax=max(tau)。 %最大時(shí)刻T=ones(size(carry1))。 %一個(gè)數(shù)字信號(hào)‘1’dsig=kron(m,T)。 %數(shù)字信號(hào)波形 Fask=fft(ask)。 %對(duì)信號(hào)進(jìn)行傅里葉變換Fask=abs(Fask(1:length(Fask)/2+1))。 %已調(diào)信號(hào)頻譜f1=(0:length(Fask)1)*f/length(Fask)/2。 %已調(diào)信號(hào)頻譜帶寬、%Butterworth 模擬低通濾波器設(shè)計(jì)%獲得Butterworth 模擬低通原型濾波器的階數(shù)及3dB 截止頻率[Ord,omega_c]=buttord(4*pi*f*,4*pi*f*,2,30,39。s39。)。%由原型濾波器向?qū)嶋H濾波器轉(zhuǎn)換，獲得濾波器的分子，分母系數(shù)[num,den]=butter(Ord,omega_c,39。s39。)。 h=tf(num,den)。 %獲得濾波器傳遞函數(shù)%濾波x=lsim(h,sig_mul,tau)。 %運(yùn)用模擬濾波器對(duì)信號(hào)進(jìn)行濾波Fx=fft(x)。 %對(duì)信號(hào)進(jìn)行傅里葉變換Fx=abs(Fx(1:length(Fx)/2+1))。 %解調(diào)信號(hào)頻譜f1=(0:length(Fx)1)*f/length(Fx)/2。 %解調(diào)信號(hào)頻譜帶寬subplot(3,2,2) %子圖分割plot(tau,x),title(39。解調(diào)輸出波形39。) %畫出濾波后的波形grid on %添加網(wǎng)格axis([Tmin Tmax ]) %設(shè)置坐標(biāo)范圍subplot(323)。plot(f1,Fx),title(39。解調(diào)信號(hào)頻譜39。)。 %解調(diào)信號(hào)頻譜grid on %添加網(wǎng)格axis([Tmin Tmax ]) %設(shè)置坐標(biāo)范圍th=。 %抽樣判決的閾值設(shè)置t_judge=(0:Lm1)*Lcarry1+Lcarry1/2。 %抽樣判決點(diǎn)的選取y=(x(t_judge))39。 %抽樣判決時(shí)刻時(shí)的信號(hào)值y_judge=1*(y=th)+0*(y=th)。 %抽樣判決信號(hào)值的0 階保持y_value=kron(y_judge,ones(size(carry1)))。 %抽樣判決后的數(shù)字信號(hào)波形n_tau=tau+。 %抽樣判決后的信號(hào)對(duì)應(yīng)的時(shí)間subplot(3,2,4) %子圖分割plot(n_tau,y_value), title(39。抽樣判決后的數(shù)字信號(hào)波形39。)%畫出抽樣判決%后%%的數(shù)字%%信號(hào)波形axis([min(n_tau) max(n_tau) ]) %設(shè)置坐標(biāo)范圍grid on subplot(3,2,5) %子圖分割plot(tau,dsig),title(39。原始信號(hào)波形與解調(diào)后的信號(hào)作對(duì)比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。已調(diào)信號(hào)與載波信號(hào)相乘波形39。)%已調(diào)信號(hào)與%載波信號(hào)相乘波形grid on axis([Tmin Tmax *A *A]) %設(shè)置坐標(biāo)范圍運(yùn)行結(jié)果：分析結(jié)果：解調(diào)信號(hào)的頻譜沒出來，檢查代碼終究還是沒找到問題所在。解調(diào)抽樣后恢復(fù)原二進(jìn)制序列。2FSK調(diào)制解調(diào)原理與系統(tǒng)仿真、2FSK調(diào)制原理 2FSK是利用數(shù)字基帶信號(hào)控制在波的頻率來傳送信息。例如，1碼用頻率f1來傳輸，0碼用頻率f2來傳輸，而其振幅和初始相位不變。故其表示式為： 式中，假設(shè)碼元的初始相位分別為和；和為兩個(gè)不同的碼元的角頻率；幅度為A為一常數(shù)，表示碼元的包絡(luò)為矩形脈沖。2FSK信號(hào)的產(chǎn)生方法有兩種：（1） 模擬法，即用數(shù)字基帶信號(hào)作為調(diào)制信號(hào)進(jìn)行調(diào)頻。如圖11（a）所示。（2） 鍵控法，用數(shù)字基帶信號(hào)及其反相分別控制兩個(gè)開關(guān)門電路，以此對(duì)兩個(gè)載波發(fā)生器進(jìn)行選通。如圖11（b）所示。這兩種方法產(chǎn)生的2FSK信號(hào)的波形基本相同，只有一點(diǎn)差異，即由調(diào)頻器產(chǎn)生的2FSK信號(hào)在相鄰碼元之間的相位是連續(xù)的，而鍵控法產(chǎn)生的2FSK信號(hào)，則分別有兩個(gè)獨(dú)立的頻率源產(chǎn)生兩個(gè)不同頻率的信號(hào)，故相鄰碼元的相位不一定是連續(xù)的。 (a) (b) 2FSK信號(hào)產(chǎn)生原理圖 由鍵控法產(chǎn)生原理可知，一位相位離散的2FSK信號(hào)可看成不同頻率交替發(fā)送的兩個(gè)2ASK信號(hào)之和，即 其中是脈寬為的矩形脈沖表示的NRZ數(shù)字基帶信號(hào)。 其中，為的反碼，即若，則；若，則。2FSK信號(hào)的頻譜特性： 由于相位離散的2FSK信號(hào)可看成是兩個(gè)2ASK信號(hào)之和，所以，這里可以直接應(yīng)用2ASK信號(hào)的頻譜分析結(jié)果，比較方便，即 2FSK信號(hào)帶寬為 式中，是基帶信號(hào)的帶寬。、2FSK解調(diào)原理2FSK的解調(diào)方法有包絡(luò)檢波法、相干檢波法、過零檢測(cè)法、差分檢測(cè)法，本實(shí)驗(yàn)用到相干檢測(cè)法。相干檢測(cè)法原理：相干檢測(cè)的具體解調(diào)電路是同步檢波器，原理框圖如下圖所示。圖中兩個(gè)帶通濾波器的作用同于包絡(luò)檢波法，起到分路作用。它們的輸出分別與相應(yīng)的同步相干載波相乘，再分別經(jīng)低通濾波器濾掉二倍頻信號(hào)，取出含基帶數(shù)字信息的低頻信號(hào)，抽樣判決器在抽樣脈沖到來時(shí)對(duì)兩個(gè)低頻信號(hào)的抽樣值vv2進(jìn)行比較判決（判決規(guī)則同于包絡(luò)檢波法），即可還原出基帶信號(hào)。假設(shè)信道噪聲為高斯白噪聲。2FSK調(diào)制解調(diào)代碼：clear allclose alli=10。 %基帶信號(hào)碼元數(shù)j=5000。a=round(rand(1,i))。 %產(chǎn)生隨機(jī)序列t=linspace(0,5,j)。f1=10。 %載波1頻率f2=5。 %載波2頻率fm=i/5。 %基帶信號(hào)頻率B1=2*f1。 %載波1帶寬B2=2*f2。 %載波2帶寬%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%產(chǎn)生基帶信號(hào)st1=t。for n=1:10 if a(n)1。 for m=j/i*(n1)+1:j/i*n st1(m)=0。 end else for m=j/i*(n1)+1:j/i*n st1(m)=1。 end endendst2=t。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%基帶信號(hào)求反for n=1:j。 if st1(n)=1。 st2(n)=0。 else st2(n)=1。 endend。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%載波信號(hào)s1=cos(2*pi*f1*t)s2=cos(2*pi*f2*t)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%調(diào)制F1=st1.*s1。 %加入載波1F2=st2.*s2。 %加入載波2fsk=F1+F2。 %fsk信號(hào)nosie=rand(1,j)。 %高斯噪聲fsk1=fsk+nosie。 %fsk信號(hào)加高斯白噪聲Ffsk1=fft(fsk1)。 %對(duì)信號(hào)進(jìn)行傅里葉變換%%%%%%%%%%%%%%%%%%%%已調(diào)信號(hào)加噪聲后的頻譜Ffsk1=abs(Ffsk1(1:length(Ffsk1)/2+1))。%%%%%%%%%%%%%已調(diào)信號(hào)加噪聲信號(hào)后的頻譜帶寬f3=(0:length(Ffsk1)1)*fm/length(Ffsk1)/2 subplot(421)。plot(t,st1)。title(39?；鶐盘?hào)39。)。%顯示基帶波形axis([0,5,1,2])。subplot(422)。plot(t,st2)。title(39。基帶信號(hào)反碼39。)。axis([0,5,1,2])。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%載波信號(hào)s1=cos(2*pi*f1*t)