基于MATLAB的2FSK调制由刀豆文库小编整理,希望给你工作、学习、生活带来方便,猜你可能喜欢“基于matlab的2fsk调制”。
基于MATLAB的2FSK信号调制
2FSK 一个FSK信号可以看成是两个不同载波的2ASK信号的叠加。其解调和解调方法和ASK差不多。2FSK信号的频谱可以看成是f1和f2的两个2ASK频谱的组合1FSK程序:
2FSK图形(包含误码率分析): close all clear all ti=0;fpefsk=[];startn=-6;endn=18;for ti=startn:endn n=1000;f1=18000000;f2=6000000;bitRate=1000000;N=50;基于MATLAB的2FSK信号调制
noise=ti;signal=source(n,N);transmittedSignal=fskModu(signal,bitRate,f1,f2,N);signal1=guian(transmittedSignal,noise);configueSignal=demoFSK(signal1,bitRate,f1,f2,N);configueSignal;P=CheckRatePe(signal,configueSignal,n)fpefsk=[fpefsk,P];end figure(8);semilogy(startn:length(fpefsk)+startn-1,fpefsk);grid on;title('Bit Error Rate Of FSK');xlabel('r/dB');ylabel('PeFSK');load PeRate save PeRate.mat fpefsk fpeask 基于MATLAB的2FSK信号调制
1.随机信号的产生
2.FSK信号调制 基于MATLAB的2FSK信号调制
3.信道加噪声后信号 基于MATLAB的2FSK信号调制
4.接收信号的解调
基于MATLAB的2FSK信号调制
5.解调出的二进制信号
误码率分析 基于MATLAB的2FSK信号调制
附录:程序
FSK clear all close all i=10;%基带信号码元数 j=5000;a=round(rand(1,i));%产生随机序列 t=linspace(0,5,j);f1=10;%载波1频率 f2=5;%载波2频率 fm=i/5;%基带信号频率 B1=2*f1;%载波1带宽 B2=2*f2;%载波2带宽
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%产生基带信号 基于MATLAB的2FSK信号调制
st1=t;for n=1:10 if a(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);基于MATLAB的2FSK信号调制
title('基带信号反码');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');%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%调制
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*st2');e_fsk=F1+F2;subplot(413);plot(t,e_fsk);title('2FSK信号')nosie=rand(1,j);fsk=e_fsk+nosie;subplot(414);基于MATLAB的2FSK信号调制
plot(t,fsk);title('加噪声信号')%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解调
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)
end end end;subplot(313);plot(t,at);axis([0,5,-1,2]);title('抽样判决后波形')