2FSK调制解调及其仿真.pdf

2FSK

一、题目

1.2FSK调制解调及其仿真。

2.相关调制解调的原理图如

ω1

带通滤相乘器低通滤

波器波器输出
输入
抽样判
Cosω1t
抽样脉冲决器

带通滤
相乘器低通滤
波器
波器

ω2
Cosω2t

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)中有反码的存在，则需

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

的表达式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.经过相乘器输出的波形再通过低通滤波器，设置低通滤波器的参
数，用一维数字滤波韩式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;

%随机信号
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)%FSK

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'

sw1=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]);

figure(4)

st1=filter(bn,1,sw1);

st2=filter(bn,1,sw2);

subplot(211);

plot(t,st1);

'经过低通滤波器后的波形')
title(sw