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