基于MATLAB的2ASK、2FSK和2PSK的調(diào)制仿真

2#波，最后將兩種信號同時輸入到抽樣判決器同時外加抽樣脈沖，最后解調(diào)出來的信號就是調(diào)制前的輸入信號。其原理圖如下圖所示：圖5(b圖5(b)非相干方式3、二進制相移鍵控(2PSK)相移鍵控是利用載波的相位變化來傳遞數(shù)字信息，而振幅和頻率保持不變。在2PSK中通常用初始相位0和兀分別表示二進制“1”和“0”。因此，2PSK信號的時域表達式為e(t)=Acos(wt+9)2FSKcn其中，9表示第n個符號的絕對相位：n{0{0兀發(fā)送“0”時發(fā)送“1”時因此，式子可以改寫為{Acoswt概率為P—Acoswt概率為1-Pc典型波形如圖所示。

一個雙極性全占空矩形脈沖序列與一個正弦載波的相乘，即e(t)=S(t)COSwt2FSKcs(s(t)=工ag(t-T)nnn這里，g(t)是脈寬為Ts的單個矩形脈沖，而an的統(tǒng)計特性為概率為P1-1概率為1-1概率為1-P2PSK信號的調(diào)制原理框圖如圖所示。2ASK信號的產(chǎn)生方法比較知識對s(t)的要求不同，在2ASK中s(t)是單極性的，而在2PSK中s(t)是雙極性的基帶信號。雙極性不歸開關(guān)電路(b)雙極性不歸開關(guān)電路(b)2PSK信號的解調(diào)通常采用相干解調(diào)法，解調(diào)器原理框圖如圖所示。四、仿真源程序和代碼clearall;closeall;A=1;fc=2;%載波頻率2Hz;N_sample=8;N=500;%碼元數(shù)Ts=1;%1baud/sdt=Ts/fc/N_sample;%波形采樣間隔t=0:dt:N*Ts-dt;T=dt*length(t);Lt=length(t);%產(chǎn)生二進制信源d=sign(randn(1,N));dd=sigexpand((d+1)/2,fc*N_sample);gt=ones(1,fc*N_sample);%NRZ波形figure(1)subplot(221);%輸入NRZ信號波形d_NRZ=conv(dd,gt);plot(t,d_NRZ(1:length(t)));axis([010001.2]);ylabel('輸入信號');subplot(222);%輸入NRZ頻譜[f,d_NRZf]=T2F(t,d_NRZ(1:length(t)));plot(f,10*log10(abs(d_NRZf).八2/T));axis([-22-5010]);ylabel('輸入信號功率譜密度(dB/Hz)');%2ASK信號ht=A*cos(2*pi*fc*t);s_2ask=d_NRZ(1:Lt).*ht;subplot(223);plot(t,s_2ask);axis([010-1.21.2]);ylabel('2ASK');[f,s_2askf]=T2F(t,s_2ask);subplot(224)plot(f,10*log10(abs(s_2askf).八2/T));axis([-fc-4fc+4-5010]);ylabel('2ASK功率譜密度(dB/Hz)');figure(2)%2PSK信號d_2psk=2*d_NRZ-1;s_2psk=d_2psk(1:Lt).*ht;subplot(221)plot(t,s_2psk);axis([010-1.21.2]);ylabel('2PSK');subplot(222);[f,s_2pskf]=T2F(t,s_2psk);plot(f,10*log10(abs(s_2pskf).八2/T));axis([-fc-4fc+4-5010]);ylabel('2PSK功率譜密度(dB/Hz)');%2FSK信號sd_2fsk=2*d_NRZ-1;s_2fsk=A*cos(2*pi*fc*t+2*pi*sd_2fsk(1:length(t)).*t);subplot(223)plot(t,s_2fsk);axis([010-1.21.2]);xlabel('t');ylabel('2FSK');subplot(224);[f,s_2fskf]=T2F(t,s_2fsk);plot(f,10*log10(abs(s_2fskf).八2/T));axis([-fc-4fc+4-5010]);xlabel('f');ylabel('2FSK功率譜密度(dB/Hz)');%兩個自定義函數(shù)如下：function[out]=sigexpand(d,M)%將輸入的序列擴展成間隔M-1個0的序列N=length(d);out=zeros(M,N);out(1,:)=d;out=reshape(out,1,M*N);function[f,sf]=T2F(t,st)%計算一個信號的傅里葉變換dt=t(2)-t(1);T=t(end);df=1/T;N=length(st);f=-N/2*df:df:N/2*df-df;sf=fft(st);sf=T/N*fftshift(sf);五、仿真結(jié)果0000號信入輸￡6備密譜率功號信入輸KSA2)zH/Bd(度密譜率功KSA2M—密譜率功KSP2KSF209876540000號信入輸￡6備密譜率功號信入輸KSA2)zH/Bd(度密譜率功KSA2M—密譜率功KSP2KSF20987654320密譜率功KSF2五、結(jié)果分析（1）2ASK信號的功率譜由連續(xù)