资源描述：
《2fsk调制解调及仿真通信原理matlab》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、2FSK调制解调及其仿真1.2FSK调制解调及其仿真。2.相关调制解调的原理图如带通滤波器带通滤波器相乘器相乘器低通滤波器低通滤波器抽样判决器Cosω1tCosω2tω1ω2抽样脉冲输出输入3.输入的信号为：S（t）=[∑аn*g(t-nTs)]cosω1t+[ān*g(t-nTs)]cosω1t；ān是аn的反码。二、仿真思路1.首先要确定采样频率fs和两个载波频率的值f1，f2。2.写出输入已经信号的表达式S(t)。由于S(t)中有反码的存在，则需要将信号先反转后在从原信号和反转信号中进行抽样。写出已调信
2、号的表达式S(t)。3.在2FSK的解调过程中，如上图原理图，信号首先通过带通滤波器，设置带通滤波器的参数，后用一维数字滤波函数filter对信号S(t)的数据进行滤波处理。输出经过带通滤波器后的信号波形。由于已调信号中有两个不同的载波（ω1,ω2）,则经过两个不同频率的带通滤波器后输出两个不同的信号波形H1,H2。4.经过带通滤波器后的2FSK信号再经过相乘器（cosω1，cosω2），两序列相乘的MATLAB表达式y=x1.*x2→SW=Hn.*Hn，输出得到相乘后的两个不同的2FSK波形h1,h2。5.
3、经过相乘器输出的波形再通过低通滤波器，设置低通滤波器的参数，用一维数字滤波韩式filter对信号的数据进行新的一轮的滤波处理。输出经过低通滤波器后的两个波形（sw1,sw2）。6.将信号sw1和sw2同时经过抽样判决器，分别输出st1,st2。其抽样判决器输出的波形为最后的输出波形st。对抽样判决器经定义一个时间变量长度i，当st1(i)>=st2(i)时，则st=0，否则st=st2(i).其中st=st1+st2。三、仿真程序程序如下：fs=2000;%采样频率dt=1/fs;f1=20;f2=120;%
4、两个信号的频率a=round(rand(1,10));%随机信号g1=ag2=~a;%信号反转，和g1反向g11=(ones(1,2000))'*g1;%抽样g1a=g11(:)';g21=(ones(1,2000))'*g2;g2a=g21(:)';t=0:dt:10-dt;t1=length(t);fsk1=g1a.*cos(2*pi*f1.*t);fsk2=g2a.*cos(2*pi*f2.*t);fsk=fsk1+fsk2;%产生的信号no=0.01*randn(1,t1);%噪声sn=fsk+no;
5、subplot(311);plot(t,no);%噪声波形title('噪声波形')ylabel('幅度')subplot(312);plot(t,fsk);title('产生的波形')ylabel('幅度')subplot(313);plot(t,sn);title('将要通过滤波器的波形')ylabel('幅度的大小')xlabel('t')figure(2)%FSK解调b1=fir1(101,[10/80020/800]);b2=fir1(101,[90/800110/800]);%设置带通参数H1=f
6、ilter(b1,1,sn);H2=filter(b2,1,sn);%经过带通滤波器后的信号subplot(211);plot(t,H1);title('经过带通滤波器f1后的波形')ylabel('幅度')subplot(212);plot(t,H2);title('经过带通滤波器f2后的波形')ylabel('幅度')xlabel('t')sw1=H1.*H1;sw2=H2.*H2;%经过相乘器figure(3)subplot(211);plot(t,sw1);title('经过相乘器h1后的波形')yl
7、abel('幅度')subplot(212);plot(t,sw2);title('经过相乘器h2后的波形')ylabel('·幅度')xlabel('t')bn-fir1(101,[2/80010/800]);%经过低通滤波器figure(4)st1=filter(bn,1,sw1);st2=filter(bn,1,sw2);subplot(211);plot(t,st1);title('经过低通滤波器sw1后的波形')ylabel('幅度')subplot(212);plot(t,st2);title('
8、经过低通滤波器sw2后的波形')ylabel('幅度')xlabel('t')%判决fori=1:length(t)if(st1(i)>=st2(i))st(i)=0;elsest(i)=st2(i);endendfigure(5)st=st1+st2;subplot(211);plot(t,st);title('经过抽样判决器后的波形')ylabel('幅度')subplot(212);plot(
