【正文】 與仿真2ASK有兩種基本解調(diào)方法：相干解調(diào)法（同步檢測(cè)法）和非相fc干解調(diào)法（包絡(luò)檢波法）。為確保無(wú)失真還原信號(hào)，必須在接收端提供一個(gè)與調(diào)制載波嚴(yán)格同步的本地載波，這是整個(gè)解調(diào)過程能否順利完好進(jìn)行的關(guān)鍵。S( t )定時(shí)脈沖 帶通 相乘 低通 抽樣判決eo( t )cosωc t2fs圖5 相干解調(diào)帶通 整流 低通 抽樣判決判決定時(shí)脈沖 S( t ) eo( t ) abcd包絡(luò)檢波器通常由整流器和低通濾波器組成，可以直接從已調(diào)波中提取原始基帶信號(hào)，結(jié)構(gòu)簡(jiǎn)單，如圖6所示。圖6 非相干解調(diào)圖7 非相干解調(diào)過程的時(shí)間波形本實(shí)驗(yàn)采用非相干解調(diào)方式。 %清除空間變量m=[1 0 1 1 0 1 1 0]。 %序列的長(zhǎng)度F=150。 %正弦載波信號(hào)的頻率A=1。 %頻率比即一個(gè)碼元寬度中的正弦周期個(gè)數(shù),為適配下面濾波器參數(shù)選%取M=500。 %一個(gè)正弦信號(hào)周期內(nèi)的時(shí)間carry1=repmat(A*sin(2*pi*f*t),1,Q)。 %一個(gè)碼元寬度內(nèi)的信號(hào)長(zhǎng)度carry2=kron(ones(size(m)),carry1)。 %調(diào)制后的信號(hào) V=。gaussian39。 %加高斯白噪聲后的已調(diào)信號(hào)sig_mul=ask.*carry2。 %長(zhǎng)度tau=(0:N1)/(M1)/f。 %最小時(shí)刻Tmax=max(tau)。 %一個(gè)數(shù)字信號(hào)‘1’dsig=kron(m,T)。 %對(duì)信號(hào)進(jìn)行傅里葉變換Fask=abs(Fask(1:length(Fask)/2+1))。 %已調(diào)信號(hào)頻譜帶寬、%Butterworth 模擬低通濾波器設(shè)計(jì)%獲得Butterworth 模擬低通原型濾波器的階數(shù)及3dB 截止頻率[Ord,omega_c]=buttord(4*pi*f*,4*pi*f*,2,30,39。)。s39。 h=tf(num,den)。 %運(yùn)用模擬濾波器對(duì)信號(hào)進(jìn)行濾波Fx=fft(x)。 %解調(diào)信號(hào)頻譜f1=(0:length(Fx)1)*f/length(Fx)/2。解調(diào)輸出波形39。plot(f1,Fx),title(39。)。 %抽樣判決的閾值設(shè)置t_judge=(0:Lm1)*Lcarry1+Lcarry1/2。 %抽樣判決時(shí)刻時(shí)的信號(hào)值y_judge=1*(y=th)+0*(y=th)。 %抽樣判決后的數(shù)字信號(hào)波形n_tau=tau+。抽樣判決后的數(shù)字信號(hào)波形39。原始信號(hào)波形與解調(diào)后的信號(hào)作對(duì)比39。已調(diào)信號(hào)與載波信號(hào)相乘波形39。解調(diào)抽樣后恢復(fù)原二進(jìn)制序列。例如，1碼用頻率f1來(lái)傳輸，0碼用頻率f2來(lái)傳輸，而其振幅和初始相位不變。2FSK信號(hào)的產(chǎn)生方法有兩種：（1） 模擬法，即用數(shù)字基帶信號(hào)作為調(diào)制信號(hào)進(jìn)行調(diào)頻。（2） 鍵控法，用數(shù)字基帶信號(hào)及其反相分別控制兩個(gè)開關(guān)門電路，以此對(duì)兩個(gè)載波發(fā)生器進(jìn)行選通。這兩種方法產(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ù)的。 其中，為的反碼，即若，則；若，則。、2FSK解調(diào)原理2FSK的解調(diào)方法有包絡(luò)檢波法、相干檢波法、過零檢測(cè)法、差分檢測(cè)法，本實(shí)驗(yàn)用到相干檢測(cè)法。圖中兩個(gè)帶通濾波器的作用同于包絡(luò)檢波法，起到分路作用。假設(shè)信道噪聲為高斯白噪聲。 %基帶信號(hào)碼元數(shù)j=5000。 %產(chǎn)生隨機(jī)序列t=linspace(0,5,j)。 %載波1頻率f2=5。 %基帶信號(hào)頻率B1=2*f1。 %載波2帶寬%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%產(chǎn)生基帶信號(hào)st1=t。 for m=j/i*(n1)+1:j/i*n st1(m)=0。 end endendst2=t。 if st1(n)=1。 else st2(n)=1。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%載波信號(hào)s1=cos(2*pi*f1*t)s2=cos(2*pi*f2*t)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%調(diào)制F1=st1.*s1。 %加入載波2fsk=F1+F2。 %高斯噪聲fsk1=fsk+nosie。 %對(duì)信號(hào)進(jìn)行傅里葉變換%%%%%%%%%%%%%%%%%%%%已調(diào)信號(hào)加噪聲后的頻譜Ffsk1=abs(Ffsk1(1:length(Ffsk1)/2+1))。plot(t,st1)?；鶐盘?hào)39。%顯示基帶波形axis([0,5,1,2])。plot(t,st2)?；鶐盘?hào)反碼39。axis([0,5,1,2])。title(39。)。title(39。)。 %加入載波1F2=st2.*s2。 %顯示載波1與基帶信號(hào)相乘波形plot(t,F1)。s1*st139。subplot(426)。title(39。)。Ffsk=fft(fsk)。 %已調(diào)信號(hào)頻譜f=(0:length(Ffsk)1)*fm/length(Ffsk)/2。plot(t,fsk)。2FSK信號(hào)39。已調(diào)信號(hào)頻譜39。figure(2)。plot(t,fsk1)。fsk信號(hào)加高斯白噪聲輸出信號(hào)39。調(diào)制信號(hào)加高斯噪聲后頻%譜%39。 %與載波1相乘[f,sf1] = t2f(t,st11)。subplot(323)。title(39。)。 %與載波2相乘[f,sf22] = t2f(t,st22)。subplot(324)。title(39。)。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽樣判決if st11(1,m*500+250)。at(1,j)=0。at(1,j)=1。subplot(325)。axis([0,5,1,2])。抽樣判決后波形39。plot(t,st1)。原基帶信號(hào)39。axis([0,5,1,2])。 %清除close all。 %抽樣頻率 fm=20e3。 final=(1/fs)*(n1)。 %載波頻率 t=0:1/fs:(final)。 %耐奎斯特頻率 %用正弦波產(chǎn)生方波 twopi_fc_t=2*pi*fm*t。 %方波幅度phi=0。 % 方波 am=1。 x(x0)=1。 %顯示基帶信號(hào)圖像plot(t,x)。 %坐標(biāo)范圍title(39。)。 %載波 psk=x.*car。 %對(duì)信號(hào)進(jìn)行傅里葉變換Fpsk=abs(Fpsk(1:length(Fpsk)/2+1))。 %已調(diào)信號(hào)頻譜帶寬subplot(312)。 axis([0 200e6 2 2])。PSK信號(hào)39。 grid on。 %顯示psk調(diào)制信號(hào)頻譜plot(f,Fpsk)。PSK信號(hào)頻譜39。 grid on。而就鍵控法來(lái)說，用數(shù)字基帶信號(hào)s(t)控制開關(guān)電路，選擇不同相位的載波輸出，這時(shí)s(t)為單極性NRZ或雙極性NRZ脈沖序列信號(hào)均可。2PSK解調(diào)原理框圖由于PSK信號(hào)的功率譜中五載波分量，所以必須采用相干解調(diào)的方式。只有對(duì)PSK信號(hào)進(jìn)行非線性變換，才能產(chǎn)生載波分量。信號(hào)經(jīng)過信道，受到信道噪聲干擾，假設(shè)噪聲為高斯白噪聲。 close all。 %抽樣頻率 fm=20e3。final=(1/fs)*(n1)。 % 載波頻率 t=0:1/fs:(final)。 %耐奎斯特頻率 %用正弦波產(chǎn)生方波 twopi_fc_t=2*pi*fm*t。 %方波幅度phi=0。 % 方波 am=1。 x(x0)=1。 %載波 psk=x.*car。 noise=vn*(randn(size(t)))。 plot(t,noise)。 title(39。)。 pskn=(psk+noise)。 %對(duì)信號(hào)進(jìn)行傅里葉變換Fpskn=abs(Fpskn(1:length(Fpskn)/2+1))。 plot(t,pskn)。 t