基于MATLAB的FSK調(diào)制

1、 2FSK 一個信號可以看成是兩個不同載波的2ASK信號的疊加。其解調(diào)和解調(diào)方法和ASK差不多。2FSK信號的頻譜可以看成是f1和f2的兩個2ASK頻譜的組合1FSK程序：2FSK圖形(包含誤碼率分析)：close allclear allti=0;fpefsk=;startn=-6;endn=18;for ti=startn:endnn=1000;f1=18000000;f2=6000000;bitRate=1000000;N=50;noise=ti;signal=source(n,N);transmittedSignal=fskModu(signal,bitRate,f1,f2,N);si

2、gnal1=gussian(transmittedSignal,noise);configueSignal=demoFSK(signal1,bitRate,f1,f2,N);configueSignal;P=CheckRatePe(signal,configueSignal,n)fpefsk=fpefsk,P;endfigure(8);semilogy(startn:length(fpefsk)+startn-1,fpefsk);grid on;title('Bit Error Rate Of FSK');xlabel('r/dB');ylabel('P

3、eFSK');load PeRatesave PeRate.mat fpefsk fpeask1.隨機(jī)信號的產(chǎn)生2.FSK信號調(diào)制3.信道加噪聲后信號4.接收信號的解調(diào)5.解調(diào)出的二進(jìn)制信號誤碼率分析附錄：程序FSKclear all close all i=10;%基帶信號碼元數(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;%基帶信號頻率 B1=2*f1;%載波1帶寬 B2=2*f2;%載波2帶寬 %產(chǎn)生基帶信號 st1=t; for n=1:10 if a(

4、n)<1; for m=j/i*(n-1)+1:j/i*n st1(m)=0; end else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end end end st2=t; %基帶信號求反 for n=1:j; if st1(n)>=1; st2(n)=0; else st2(n)=1; end end; figure(1); subplot(411); plot(t,st1); title('基帶信號'); axis(0,5,-1,2); subplot(412); plot(t,st2); title('基帶信號反碼'

5、;); axis(0,5,-1,2); %載波信號 s1=cos(2*pi*f1*t) s2=cos(2*pi*f2*t) subplot(413),plot(s1); title('載波信號1'); subplot(414),plot(s2); title('載波信號2'); %調(diào)制 F1=st1.*s1;%加入載波1 F2=st2.*s2;%加入載波2 figure(2); subplot(411); plot(t,F1); title('s1*st1'); subplot(412); plot(t,F2); title('s2*st

6、2'); e_fsk=F1+F2; subplot(413); plot(t,e_fsk); title('2FSK信號') nosie=rand(1,j); fsk=e_fsk+nosie; subplot(414); plot(t,fsk); title('加噪聲信號') %相干解調(diào) st1=fsk.*s1;%與載波1相乘 f,sf1 = T2F(t,st1);%通過低通濾波器 t,st1 = lpf(f,sf1,2*fm); figure(3); subplot(311); plot(t,st1); title('與載波1相乘后波形'); st2=fsk.*s2;%與載波2相 f,sf2 = T2F(t,st2);%通過低通濾波器 t,st2 = lpf(f,sf2,2*fm); subplot(312); plot(t,st2); title('與載波2相乘后波形'); for m=0:i-1; %抽樣判決 if st1(1,m*500+250)<0.25; for j=m*500+1:(m+1)*500; at(1,j)=0; end else for j=m*5