基于matlab的2FSK调制解调系统仿真设计

在日常的生活中，通信是人们用来传递信息的方式。随着数字系统的飞速发展，对数字系统的性能和调制解调技术要求也越来越高。同时，由于计算技术的发展，通信系统的仿真已日益普遍，已逐渐成为今天设计和分析通信系统的主要工具。 本次设计将使用MATLAB软件设计函数对2FSK调制解调技术进行仿真和研究。 本文在第一章中介绍了通信系统的组成、MATLAB的使用。第二章深入分析了2FSK的调制解调原理理论知识，熟悉了原理后，在第三章中用MATLAB编程进行仿真和研究。本设计主要实现2FSK调制解调过程的仿真。最后一章对数字调制与解调作了一个总结

fs=2000;   %抽样频率 
dt=1/fs; 
f1=20;     %定义两列载波的频率 
f2=100; 
a=round（rand（110））;     %产生二进制随机序列 
g1=a; 
g2=~a; 
g11=（ones（12000））‘*g1;  %产生方波信号 
g1a=g11（:）‘;  
g21=（ones（12000））‘*g2; 
g2a=g21（:）‘; 
t=0:dt:10-dt; 
t1=length（t）; 
tuf1=cos（2*pi*f1.*t）;
tuf2=cos（2*pi*f2.*t）; 
 
subplot（311） 
n=0:9; 
x=square（150）; 
stem（[0:9]a*x）; 
grid; 
xlabel（‘二进制随机序列‘） 
ylabel（‘幅度‘） 
 
subplot（312）; 
plot（ttuf1）; 
title（‘频率为f1的余弦波‘） 
ylabel（‘幅度‘）  
subplot（313）; 
plot（ttuf2）; 
title（‘频率为f2的余弦波‘） 
ylabel（‘幅度‘） 
 
 
figure（2） 
fsk1=g1a.*tuf1; 
fsk2=g2a.*tuf2; 
fsk=fsk1+fsk2; 
no=0.01*randn（1t1）;%噪声 
sn=fsk+no; 
subplot（211）; 
plot（tfsk）; 
title（‘2fsk波形‘） 
ylabel（‘幅度‘） 
 
subplot（212）; 
plot（tsn）; 
title（‘加入高斯噪声后的2fsk波形‘） 
ylabel（‘幅度的大小‘） 
xlabel（‘t‘） 
 
 
figure（3）   %FSK解调 
b1=fir1（101[10/800 20/800]）; 
b2=fir1（101[90/800 110/800]）;      %设置带宽参数 
H1=filter（b11sn）;       %b1为分子1为分母sn为滤波器输入序列 
H2=filter（b21sn）;       %噪声信号同时通过两个滤波器 
subplot（211）; plot（tH1）; 
 
title（‘经过带通滤波器H1后的波形‘）      %画出经过H1滤波器后的波形 
ylabel（‘幅度‘）; 
subplot（212）; 
plot（tH2）;                            %画出经过滤波器二后的波形 
title（‘经过带通滤波器H2后的波形‘） 
 
ylabel（‘幅度‘） 
xlabel（‘t‘） 
sw1=H1.*H1;                        %相干解调乘以同频同相的载波 
sw2=H2.*H2;%经过相乘器 
  
 
figure（4） 
subplot（211）; 
plot（tsw1）; 
title（‘经过相乘器h1后的波形‘）      %画出乘以同频同相载波后的波形 
ylabel（‘幅度‘） 
subplot（212）; 
plot（tsw2）; 
title（‘经过相乘器h2后的波形‘） 
ylabel（‘.幅度‘） 
xlabel（‘t‘） 
bn=fir1（101[2/800 10/800]）;             %经过低通滤波器 
 
 
figure（5） 
st1=filter（bn1sw1）; 
st2=filter（bn1sw2）; 
subplot（211）; plot（tst1）; 
title（‘经过低通滤波器sw1后的波形‘）      %经过低通滤波器滤出频率 
ylabel（‘幅度‘）                        %为f1f2的基带调制信号波形 
subplot（212）; 
plot（tst2）; 
title（‘经过低通滤波器sw2后的波形‘） 
ylabel（‘幅度‘） 
xlabel（‘t‘） 
%判决 
for i=1:length（t） 
 if（st1（i）>=st2（i）） 
   st（i）=0; 
  else st（i）=st2（i）; 
 end  
end 
 
 
figure（6） 
st=st1+st2; 
subplot（211）; 
plot（tst）; 
title（‘经过抽样判决器后的波形‘）         %画出经过抽样判决的波形     
ylabel（‘幅度‘） 
subplot（212）; 
plot（tsn）; 
title（‘原始的波形‘） 
ylabel（‘幅度‘） 
xlabel（‘t‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-10-10 15:39  源码\
     文件     1047552  2016-04-21 17:05  源码\基于matlab的2FSK调制解调系统仿真设计（修改）.doc
     目录           0  2016-04-03 09:49  源码\程序\
     文件        2420  2016-03-31 14:24  源码\程序\FSK.m
     文件       10180  2016-03-31 15:38  源码\程序\FSK_GUI.m
     目录           0  2018-10-10 15:39  演示\
     文件     3317612  2016-03-31 19:20  演示\9742.exe
     文件          64  2018-10-10 16:34  说明.txt
     目录           0  2018-10-10 15:38  需求\
     文件       11903  2016-03-29 14:51  需求\功能说明.docx