2FSK调制解调和仿真设计.pdf

2FSK调制解调及其仿真
一、题目

1.2FSK调制解调及其仿真。
2.相关调制解调的原理图如

31

32
3
Cos2t

3.输入的信号为：

S(t)=[Ean*g(t-nTs)]cos31t+[m*g(t-nTs)]cos31t;

m是an的反码。

二、仿真思路

1.首先要确定采样频率fs和两个载波频率的值f1,f2。

2.写出输入已经信号的表达式s(t)o由于s(t)中有反码的存在，贝y

需要将信号先反转后在从原信号和反转信号中进行抽样。写出已调信

号的表达式S(t)o

3.在2FSI的解调过程中，如上图原理图，信号首先通过带通滤波器，

.专业资料.

设置带通滤波器的参数，后用一维数字滤波函数filter对信号S(t)

.专业资料.

的数据进行滤波处理。输出经过带通滤波器后的信号波形。由于已调

信号中有两个不同的载波(31,32),则经过两个不同频率的带通滤波

器后输出两个不同的信号波形H1,H2

4.经过带通滤波器后的2FS信号再经过相乘器(C0S31,C0S32),

两序列相乘的MATLAB达式y=x1.*x2—SW=Hn.*Hn输出得到相乘

后的两个不同的2FSI波形h1,h2。

5.经过相乘器输出的波形再通过低通滤波器，设置低通滤波器的参数，

用一维数字滤波韩式filter对信号的数据进行新的一轮的滤波处理。输出

经过低通滤波器后的两个波形(sw1,sw2。

6.将信号sw和sw2同时经过抽样判决器，分另I」输出st1,st2。其抽样

判决器输出的波形为最后的输出波形st。对抽样判决器经定义一个时间变

量长度i,当st1(i)>=st2(i)时，贝Ust=O,否则st=st2(i).其中

st=st1+st2。
三、仿真程序

程序如下：

fs=2000;%采样频率
dt=1/fs;
f仁20;
f2=120;%两个信号的频率
a=round(rand(1,10));%随机信号
g1=a
g2=~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;沪生的信号
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)%FS解调
b1=fir1(101,[10/80020/800]);
b2=fir1(101,[90/800110/800]);%设置带通参数
H1=filter(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')
sw仁H1.*H1;
sw2=H2.*H2;%经过相乘器
figure(3)subplot(211);
plot(t,sw1);
title('经过相乘器h1后的波形’)ylabel('幅度')subplot(212);
plot(t,sw2);
title('经过相乘器h2后的波形’)
ylabel('•幅度')
xlabel('t')
bn-fir1(101,[2/80010/800]);%S过低通滤波器
figure(4)
st1=filter(bn,1,sw1);
st2=filter(bn,1,sw2);
subplot(211);
plot(t,st1);
title('经过低通滤波器swl后的波形’)ylabel('幅度')subplot(212);
plot(t,st2);
title('经过低通滤波器sw2后的波形’)ylabel('幅度')xlabel('t')

.专业资料.

%判决
fori=1:length(