利用matlab实现2FSK调制解调的仿真论文.docx

1、利用matlab实现2FSK调制解调的仿真论文潍坊学院专业课综合课程设计说明书利用 matlab 实现 2FSK调制解调的仿真系 部： 信息控制与工程学院专 业： 电子信息工程班 级：学生: 学 号:指导教师 ：2012年 12 月 01 日. . .1 MATLAB 软件简介 12 理论分析 . . 22.1 2FSK信号的产生 22.2 2FSK信号的解调方式 33 MATLAB 实现程序 53.1 生成 2FSK的程序 53.2 2FSK调制解调系统进行设计和仿真源程序 75心得体会 . . 166参考文献 . . 17. . .1 MATLAB 软件简介MATLAB 是目前国际上流行的

2、进行科学研究、工程计算的软件。它起源于矩阵运算，并已经发展成为一种高度集成的计算机语言。 MATLAB具有强大的数学运算能力、方便实用的绘图功能及语言的高度集成性。 除具备卓越的数值计算能力之外，它还提供了专业水平的符号计算、文字处理、可视化建模仿真、实时控制等功能。在通信领域 MATLAB更是优势明显，因为通信领域中很多问题是研究系统性能的，传统的方法只有构建一个实验系统， 采用各种方法进行测量， 才能得到所需的数据，这样不仅需要花费大量的资金用于实验系统的构建， 而且系统构建周期长，系统参数的调整也十分困难。而 MATLAB的出现使得通信系统的仿真能够用计算机模拟实现， 免去构建实验系统的

3、不便， 而且操作十分简便， 只需要输入不同的参数就能得到不同情况下系统的性能， 而且在结构的观测和数据的存储方面也比传统的方式有很多优势。因而 MATLAB在通信仿真领域得到越来越多的应用。. . .2理论分析2.1 2FSK 信号的产生2FSK 是利用数字基带信号控制在波的频率来传送信息。例如， 1 码用频率f1 来传输， 0 码用频率 f2 来传输，而其振幅和初始相位不变。故其表示式为_s2FSK (t) an g(t nTs ) cos( 1tn ) an g(tnTs ) cos( 2t n )nn式中，假设码元的初始相位分别为n 和 n ；12 f1 和 22 f2 为两个不同的码元

4、的角频率；幅度为 A 为一常数，表示码元的包络为矩形脉冲。2FSK 信号的产生方法有两种：模拟法，即用数字基带信号作为调制信号进行调频。键控法，用数字基带信号 g(t) 及其反 g(t ) 相分别控制两个开关门电路，以此对两个载波发生器进行选通。这两种方法产生的 2FSK信号的波形基本相同，只有一点差异，即由调频器产生的 2FSK信号在相邻码元之间的相位是连续的，而键控法产生的 2FSK信号，则分别有两个独立的频率源产生两个不同频率的信号， 故相邻码元的相位不一定是连续的。. . .ak1011001s1 (t)s2 (t)cos (w1 t+ n)cos (w2 t+ n)s1(t) cos

5、(w1t+ n)s2 (t) cos (w2 t+ n)2 FSK 信号2.2 2FSK 信号的解调方式ttttttt非相干解调经过调制后的 2FSK数字信号通过两个频率不同的带通滤波器 f1 、f2 滤出不需要的信号，然后再将这两种经过滤波的信号分别通过包络检波器检波， 最后将两种信号同时输入到抽样判决器同时外加抽样脉冲， 最后解调出来的信号就是调制前的输入信号。其原理图如下图所示：. . .带通滤波器F1 包络检波器输入 输出抽样脉冲 抽样判决器带通滤波器F2 包络检波器非相干方式原理图相干解调根据已调信号由两个载波 f1 、f2 调制而成，则先用两个分别对 f1 、f2 带通的滤波器对已

6、调信号进行滤波，然后再分别将滤波后的信号与相应的载波 f1 、f2 相乘进行相干解调，再分别低通滤波、用抽样信号进行抽样判决器即可。原理图如下：cos2 f1t带通滤波器 F1相乘器低通滤波器输入抽样脉冲抽样判决器输出带通滤波器 F2相乘器低通滤波器cos2 f2t相干方式原理图. . .3 Matlab 实现程序3.1 生成 2FSK的程序源程序代码：clear allclose alli=10;% 基带信元数j=5000;a=round(rand(1,i);% 产生随机序列t=linspace(0,5,j);f1=10;%载波 1频率f2=5;%载波 2频率fm=i/5;% 基带信号频率B

7、1=2*f1;%载波 1带宽B2=2*f2;%载波 2带宽%产生基带信号 st1=t;for n=1:10if a(n)=1;st2(n)=0;elsest2(n)=1;endend;figure(1);subplot(411);plot(t,st1);title( 基带信号 );axis(0,5,-1,2);subplot(412);plot(t,st2);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

8、),plot(s2);title( 载波信号 2);%调制F1=st1.*s1;% 加入载波 1F2=st2.*s2;% 加入载波 2. . .figure(2);subplot(311);plot(t,F1);title(s1*st1);subplot(312);plot(t,F2);title(s2*st2);e_fsk=F1+F2;subplot(313);plot(t,e_fsk);title(2FSK 信号 )3.2 2FSK 调制解调系统进行设计和仿真源程序Fc=10; % 载频Fs=35; % 系统采样频率Fd=1; % 码速率N=Fs/Fd;df=10;numSymb=25;%

9、 进行仿真的信息代码个数M=2; % 进制数SNRpBit=60;% 信噪比SNR=SNRpBit/log2(M);%60seed=12345 54321;numPlot=15;x=randsrc(numSymb,1,0:M-1);% 产生 25 个二进制随机码figure(1)stem(0:numPlot-1,x(1:numPlot),bx);% 显 15 个码元，杆图，从 x 的前十五个随机数中选取title( 二进制随机序列 )xlabel(Time);ylabel(Amplitude);%调制y=dmod(x,Fc,Fd,Fs,fsk,M,df);% 数字带通调制 numModPlot

10、=numPlot*Fs; %15*40t=0:numModPlot-1./Fs;% 数组除法（仿真时间）figure(2)plot(t,y(1:length(t),b-);axis(min(t) max(t) -1.5 1.5);. . .title( 调制后的信号 )xlabel(Time);ylabel(Amplitude);%在已调信号中加入高斯白噪声randn(state,seed(2); % 生成 -2 到+2 之间的随机数矩阵y=awgn(y,SNR-10*log10(0.5)-10*log10(N),measured,dB);% 在 已调信号中加入高斯白噪声figure(3)pl

11、ot(t,y(1:length(t),b-);% 画出经过信道的实际信号axis(min(t) max(t) -1.5 1.5);title( 加入高斯白噪声后的已调信号 )xlabel(Time);ylabel(Amplitude);%相干解调figure(4)z1=ddemod(y,Fc,Fd,Fs,fsk/eye,M,df);title( 相干解调后的信号的眼图 )%带输出波形的相干 M元频移键控解调 figure(5) stem(0:numPlot-1,x(1:numPlot),bx); hold on;stem(0:numPlot-1,z1(1:numPlot),ro);hold o

12、ff;axis(0 numPlot -0.5 1.5);title( 相干解调后的信号原序列比较 )legend( 原输入二进制随机序列 , 相干解调后的信号 )xlabel(Time);ylabel(Amplitude);%非相干解调figure(6)z2=ddemod(y,Fc,Fd,Fs,fsk/eye/noncoh,M,df);title( 非相干解调后的信号的眼图 )%带输出波形的非相干 M元频移键控解调 figure(7)stem(0:numPlot-1,x(1:numPlot),bx);hold on;stem(0:numPlot-1,z2(1:numPlot),ro);hold

13、 off;axis(0 numPlot -0.5 1.5);title( 非相干解调后的信号 )legend( 原输入二进制随机序列 , 非相干解调后的信号 )xlabel(Time);ylabel(Amplitude);%误码率统计errorSym ratioSym=symerr(x,z1);. . .figure(8)simbasebandex(0:1:5);title( 相干解调后误码率统计 )errorSym ratioSym=symerr(x,z2);figure(9)simbasebandex(0:1:5);title( 非相干解调后误码率统计 )%滤除高斯白噪声Delay=3;R

14、=0.5;PropD=0; % 滞后 3syf,tf=rcosine(Fd,Fs,fir,R,Delay); % 升余弦函数yo2,to2=rcosflt(y,Fd,Fs,filter,yf);%加入高斯白噪声后的已调信号和经过升余弦滤波器后的已调信号 t=0:numModPlot-1./Fs;figure(10)plot(t,y(1:length(t),r-);hold on;plot(to2,yo2,b-);% 滤出带外噪声hold off;axis(0 30 -1.5 1.5);xlabel(Time);ylabel(Amplitude);legend( 加入高斯白噪声后的已调信号 ,

15、经过升余弦滤波器后的已调信号)title( 升余弦滤波前后波形比较 )eyediagram(yo2,N);% 眼图title( 加入高斯白噪声后的已调信号的眼图 )4仿真结果. . .s1*st110-100.511.522.533.544.55s2*st210-10 0.5 1 1.5 2 2.5 3 3.5 4 4.5 52FSK 信 号10-10 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5基带信号210-100.511.522.533.544.55基带信号反码210-100.511.522.533.544.55载波信号110-10500100015002000250030

16、003500400045005000载波信号210-10 500 1000 1500 2000 2500 3000 3500 4000 4500 5000. . .二进制随机序列10.90.80.7e 0.6dutilp 0.5mA0.40.30.20.100 2 4 6 8 10 12 14Time调制后的信号1.510.5edutilp 0mA-0.5-1-1.50 2 4 6 8 10 12 14Time. . .加入高斯白噪声后的已调信号1.510.5edutilp 0mA-0.5-1-1.50 2 4 6 8 10 12 14Time相干解调后的信号的眼图10.50-0.50 0.1

17、 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9. . .相干解调后的信号原序列比较1.5原输入二进制随机序列相干解调后的信号1edutilp 0.5mA0-0.5510150Time非相干解调后的信号的眼图10.50-0.50 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9. . .edutilpmA非相干解调后的信号1.5原输入二进制随机序列非相干解调后的信号10.50-0.50 5 10 15010Time相干解调后误码率统计Theoretical SERTheoretical BERSimulated SERSimulated BERREBdna

18、RES-110-210-3100 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5EbNo (dB). . .非相干解调后误码率统计010Theoretical SERTheoretical BERSimulated SERSimulated BER-110REBdnaRES-210-3100 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5EbNo (dB)升余弦滤波前后波形比较1.5加入高斯白噪声后的已调信号经过升余弦滤波器后的已调信号10.5edutilp 0mA-0.5-1-1.50 5 10 15 20 25 30Time. . .加入高斯白噪声后的已调信号的眼图

19、1.510.5edutilp 0mA-0.5-1-1.5-0.4 -0.2 0 0.2 0.4Time5心得体会两周的课程设计在忙忙碌碌中一晃而过。 经历了一次次的困惑， 却积累了一. . .定的知识。在整个课程设计中的过程中遇到的问题主要有以下三点： 第一，基础知识的不牢固， 主要表现在一些常用的电路的形式和功能不清楚， 对书本上的容理解不够透彻； 第二，对一些常用的应用软件缺少应用， 体现在华电路图和系统仿真的时候，对这些软件操作的不熟练，浪费了很多时间。第三，相关知识掌握不够全面，缺少系统设计额仿真的经验。这次课程设计进一步端正了我的学习态度，学会了实事求是，严谨的作风，提高了动手能力。

20、 对自己要严格要求， 不能够一知半解要力求明明白白。 急于求成是不好的，不仅会浪费时间还会适得其反。在我看来，懂得少并不可怕，可怕的是不向别人虚心学习。 没有人生下就知道什么， 也没有人生下来就很聪明。 即使是天才，也要通过后天的努力，才获得成功的。我觉得动手之前，头脑里必须清楚怎么做，通过我的不懈努力，在这方面我总会得到提高，这一点我坚信。因为别人能做到的，我也一定能做到。在此次的课程设计中我最大的体会就是进一步认识到了理论联系实践的重要性。一份耕耘一分收获。 通过一个星期的实习， 让我明白科学的思维方法和学习方法是多么重要，只有这样才能够让自己工作更完美。总而言之，在此次课程设计让我学到了

21、好多平时在课堂上学不到得东西， 增加了我的知识运用能力，为我走向社会奠定了一个号的基础。通过对 matlab 这个软件的学习，使我对通信原理又有了进一步的认识。在以往的学习中我多是注重理论知识没注重实践， 这次实践使我对课本知识有了新的理解。6参考文献1保锁 付晓梅 侯永宏编著、 现代通信原理试验 、天津大学电信学院通信试验室2曹志刚 , 钱亚生 . 现代通信原理 . 北京 : 清华大学 ,1992:254-256.3晁冰 , 东生 , 雍爱霞 . 最小频移键控系统实现技术的仿真研究 . 现代电子. . .技术 ,2002,(12):88-89.4东生 , 雍爱霞 , 左洪浩 .System View 系统设计及仿真入门与应用 . 北京 :电子工业 ,2002:120-156.5青松 , 程岱松 , 武建华 . 数字通信系统的 System View 仿真与分析 . 北京 :电子工业 ,2001:53-135.6 东升 . 信号与电子系统原理及 EDA仿真 M. 北京 : 中国科技大学 , 2000.6. . .