基于MATLAB的FSK调制解调

1、目录 1 引言 .3 1.1 课程设计目的 .3 1.2 课程设计要求 .3 1.3 课程设计步骤 .3 1.4 设计平台简介 .4 2 设计原理 .5 2.1 FSK 基本原理.5 2.2 FSK 调制原理.5 2.3 FSK 解调原理.6 3 仿真实现过程 .7 3.1 FSK 信号的产生.7 3.2 FSK 的解调.8 3.3 叠加噪声的 FSK 解调 .9 4 遇到问题及解决办法.12 4.1 编程中问题解答 .12 4.2 程序解调过程中的问题 .12 结束语.13 参考文献.14 基于 MATLAB 的 FSK 调制解调 摘摘 要：要： 本课程设计利用 MATLAB 集成环境下的

2、M 文件，编写程序来实现 FSK 的调制解调，并绘制出解调前后的时域和频域波形及叠加噪声时解调前后 的时频波形，并观察解调前后频谱有何变化以加深对 FSK 信号解调原理的理解。 对信号叠加噪声，并进行解调，绘制出解调前后信号的时频波形，改变噪声功 率进行解调，根据运行结果和波形来分析该解调过程的正确性及信道对信号传 输的影响。完成整个 FSK 的调制解调过程。程序开发平台为 MATLAB7.1，使 用其自带的 M 文件实现。运行平台为 Windows 2000。 关键词关键词：程序设计；FSK ；调制解调 ；MATLAB7.1；M 文件 1 引言 本课程设计是利用 MATLAB 集成环境下的

3、M 文件，编写程序来实现 FSK 的调制解调，并绘制出解调前后的时域和频域波形及叠加噪声时解调前后的时 频波形，根据运行结果和波形来分析该解调过程的正确性及信道对信号传输的 影响。 1.1 课程设计目的 此次课程设计的目的是熟悉 MATLAB 中 M 文件的使用方法，编写 M 文件 实现 FSK 的调制和解调，绘制出 FSK 信号解调前后在时域和频域中的波形， 观察调解前后频谱的变化，再对信号进行噪声叠加后解调同样绘制解调前后的 信号时频波形，最后改变噪声功率进行调解，分析噪声对信号传输造成的影响， 加深对 FSK 信号解调原理的理解。 1.2 课程设计要求 熟悉 MATLAB 中 M 文件的

4、使用方法，并在掌握 FSK 调制解调原理的基础 上，编写出 FSK 调制解调程序。在 M 文件环境下运行程序绘制出 FSK 信号解 调前后在时域和频域中的波形，观察波形在解调前后的变化，对其作出解释， 同时对信号加入噪声后解调，得到解调后的时频波形，分析噪声对信号传输造 成的影响。解释所得到的结果。 1.3 课程设计步骤 本课程设计采用 M 文件编写的方法实现二进制的 FSK 的调制与解调，然后 在信号中叠加高斯白噪声。一，调用 dmode 函数实现 FSK 的解调，并绘制出 FSK 信号调制前后在时域和频域中的波形，两者比较。二，调用 ddemod 函数解 调，绘制出 FSK 信号解调前后在

5、时域和频域中的波形，两者比较。三，调用 awgn 函数在新海中叠加不同信噪比的噪声，绘制在各种噪声下 的时域频域图。 最后分析结果。 1.4 设计平台简介 Matlab 是美国 MathWorks 公司开发的用于概念设计,算法开发,建模仿真, 实时实现的理想的集成环境。是目前最好的科学计算类软件。 作为和 Mathematica、Maple 并列的三大数学软件。其强项就是其强大的矩阵计 算以及仿真能力。Matlab 的由来就是 Matrix + Laboratory = Matlab，这个软件在 国内也被称作矩阵实验室 。Matlab 提供了自己的编译器：全面兼容 C+以 及 Fortran

6、两大语言。Matlab 7.1 于 2005.9 最新发布-完整版,提供了 MATLAB，SIMULINK 的升级以及其他最新的 75 个模块的升级。7.1 版本提 高了产品质量，同时也提供了新的用于数据分析、大规模建模、固定点开发、 编码等。所以 Matlab 是工程师，科研工作者手上最好的语言，最好的工具和环 境。 2 设计原理 2.1 FSK 基本原理 频移键控是利用载波的频率变化来传递数字信息。在 2FSK 中载波的频率 随二进制基带信号在 f1 和 f2 两个频率点间变化。故其表达式为 e2fsk(t)=Acos(1t+n) 发送“1”时；e2fsk(t)=Acos(2t+n) 发送

7、“ 0”时。一个 2FSK 可以看成两个不同载频的 2ASK 信号的叠加。因此 2FSK 信 号时域表达式又可以写成 e2fsk(t)= ang(t-nTs)cos(1t+n)+ng(t-nTs)cos(2+n) n n a (2-1) 式中：g(t)为单个矩形脉冲，脉宽为 Ts；an=1 时，概率为 P; an=0 时，概率为 1-P;n 是 an 的反码，若 an=1，则n=0；若若 an=0，则n=1，于是aaa n =1 时，概率为 1-P; n =0 时，概率为 P;aa n 和 n 分别是第 n 个信号码元（1 或 0）的初始相位。在移频键控中 n 和 n 不携带信息，通常和令 n

8、 和 n 为零。因此 2FSK 信号表达式可简化为: e2fsk(t)=s1(t)cos1t+ s2(t)cos2t ( 2-2) s1(t)= ang(t-nTs) (2-3) n s2(t)= ng(t-nTs)(2-4) n a 2.2 FSK 调制原理 在二进制频移键控中，幅度恒定不变的载波信号的频率随着输入码流的变化而 切换（称为高音和低音，代表二进制的 1 和 0） 。产生 FSK 信号最简单的方法 是根据输入的数据比特是还是，在两个独立的振荡器中切换。采用这种方 法产生的波形在切换的时刻相位是不连续的，因此这种 FSK 信号称为不连续 F SK 信号。 ）其实现如图一 所示： 振

9、荡器 FH 振荡器 FL 放大输出 输入 图 2-1 非连续相位 FSK 的调制方式 2.3 FSK 解调原理 对于 FSK 信号的解调方式很多：相干解调、滤波非相干解调、正交相乘非相干 解调。而 FSK 的非相干解调一般采用滤波非相干解调，解调原理是将 2FSK 信 号分为上下两路 2ASK 信号分别进行解调，然后判决，这里的抽样判决是直接 比较两路信号抽样值的大小，可以不专门设置门限。判决则应与调制规则相呼 应，调制时若规定“1”符号为对应载波频率 f1,则接受时上支路的样值大小， 应判为“1”反之则判为“0”。 带通滤波 FH 带通滤波 FL 包络检 波 包络检 波 判决 FS K 图

10、2-2 2FSK 信号非相干解调原理图 3 仿真实现过程 3.1 FSK 信号的产生 调用 domde 函数实现 FSK 的调制： y=dmod(x,Fc,Fd,Fs,fsk,M,df); (3-1) 其中 y 为已调信号；x 为调制信号；Fc=10 为载波频率；Fs=40 系统采样频率； Fd=1 码元速率。 此程序输入的是二进制的 FSK 信号且为数字信号，所谓为连续的二进制图 形如图三所示： 图 3-1 输入的二进制随即序列 对该信号进行频谱分析在 M 文件下得到了调制后的频域和时域图设定好各个的 参数 图 3-2 调制后的信号 图 3-3 调制后的频谱图 通过调制，从所得的图形可知输入

11、的二进制随机信号没受到噪声干扰的时候是 很规则的信号调制后的信号 3.2 FSK 的解调 对调制好的 FSK 信号进行非相干解调。输入的 FSK 中频信号分别经过中心频 率为 fH、fL 的带通滤波器，然后分别经过包络检波，包络检波的输出在 t=kTb。时抽样（其中 k 为整数） ，并且将这些值进行比较。根据包络检波器输 出的大小，比较器判决数据比特是还是。 图 3-4 解调前后的频谱图比较 从所得到的结果来看，经过解调后，在没有噪声的干扰下，解调出来的信号在 频谱上和原信号十分吻合，基本没有误差。 3.3 叠加噪声的 FSK 解调 在信号中如果存在了不需要的电信号，我们称之为噪声。而在通信系

12、统中， 我们一般需要处理的是高斯白噪声。 一般在物理上把它翻译成白噪声 （white noise）。 白噪声是指功率谱密度在整个频域内均匀分布的噪声。 所有频率具有 相同能量的随机噪声称为白噪声。从我们耳朵的频率响应听起来它是非常明 亮的“咝”声。白噪声或白杂讯，是一种功率频谱密度为常数的随机信号或 随机过程。换句话说，此信号在各个频段上的功率是一样的，由于白光是由 各种频率（颜色）的单色光混合而成，因而此信号的这种具有平坦功率谱的 性质被称作是 “白色的”，此信号也因此被称作白噪声。相对的，其他不具 有这一性质的噪声信号被称为有色噪声。 理想的白噪声具有无限带宽，因而其能量是无限大，这在现实

13、世界是不 可能存在的。实际上，我们常常将有限带宽的平整讯号视为白噪音，因为这 让我们在数学分析上更加方便。然而，白噪声在数学处理上比较方便，因此 它是系统分析的有力工具。一般，只要一个噪声过程所具有的频谱宽度远远 大于它所作用系统的带宽，并且在该带宽中其频谱密度基本上可以作为常数 来考虑，就可以把它作为白噪声来处理。 我们要分析的是在加了高斯白噪 声后的 FSK 信号。 调用 awgn 函数来实现给 FSK 信号的噪声叠加： yl=awgn(y,20); (3-2) yl 为加了噪声后的信号，y 为需要叠加噪声的信号。 调用前面已经产生的调制信号，我们设定噪声比为 20，则可以得到加入噪声后

14、的已调信号和原二进制随机序列的比较。如图 6 所示 图 3-5 原始信号调制后信号和加入高斯白噪声后的已调信号比较 从图中可知，加入噪声之后信号的波动明显增大，幅度超出了“1” ， “-1”的范 围，信号的失真情况很明显，但对信号的波动没有很大的改变。由此说明，叠 加噪声之后信号传输效果会受到影响，但正确性不会改变。 图 3-6 调制后的频谱图 3-7 加入高斯白噪声解调后的频谱图（snr=10） 3-8 加入高斯白噪声解调后的频谱图（snr=3） 加入噪声后，信号的变化很明显。改变噪声的大小：当信噪比为 10 时， 失真不是非常明显，与原波形基本上没大的变化。当信噪比为 3 时，噪声大了 之

15、后对信号的影响非常大。说明 FSK 信号具有一定的抗噪声能力。 图 3-9 加入噪声后非相干解调后信号 加入噪声后，对信号进行非相干解调，与原信号没有变化，任然能正确的 解调出原始信号，这是因为虽然加入了噪声，但数字信号的抗噪性能较模拟信 号而言更强。并且我们取的数据组数有限，不能很好的反应其真实的抗噪能力。 通过上述的操作之后，信道中噪声对信号的传输有一定的影响，噪声越大 影响越明显，掩盖了原始信号，所以我们在平常信号传输中，尽量减少信道中 噪声的干扰，保证信号的传输效率。 4 遇到问题及解决办法 4.1 编程中问题解答 在此次课程实际过程中，碰到了很多的问题，而解决这些问题使我学得了 很多

16、的新知识，同时也重拾以前遗忘的知识，加深了对现有知识的巩固。在同 学们和老师的帮助下，将这些问题一一解决。现将问题做如下总结。 首先是在编写程序时，虽然在网上找到了相关程序，但这些程序是网上程 序员编写的，其中有许多的函数是我们的 MATLAB 虽然有，但是我们所谓涉 及过的，因而无法读懂这些函数 对其作出解释。而我要做的是将这些函数变为 自己现有的程序在 MATLAB 的函数库中现有的我们接触过的能自己组织语言 描述的。例如： yl=awgn(y,SNR-10*log10(0.5)-10*log10(N),measured,dB);此函数在 以前的课程中并未看到过，在问过老师之后，才真正明白

17、其含义：次句程序的 意思就和我们的 yl=awgn(y,20); 一样，在信号中添加高斯白噪声，对于这一系列的函数和程序用语问题，不懂 的经过老师和同学的解答，都得到了很好的解释。也是我懂得在平时学习中不 能只局限于课本知识，对于课后的各类文章也要去多多留意。 4.2 程序解调过程中的问题 对于程序中做出 FSK 解调后的频谱图的程序我们可以单独作为一个 M 文 件在 MATLAB 中运行，这样就减少了程序在各个 M 文件中存在的问题，同样 也就较少了程序过于繁琐的问题。 对于矩阵 n=0:1:length(y)-1;w=0:1:1000-1*pi/500;）在没有用 1000-1 之前，程序

18、一直运行有错误。在请教了我们班的同学之后，按照矩阵乘 法的规则第一个矩阵的列数要和第二个矩阵的行数相等，否则不能相乘。所以 在定义 1000 个点，与前一个矩阵相对应。 结束语 紧张气氛并存的机房中，度过了很多个日子。终于完成了通信原理的课程 设计在此我要谢谢在整个课程设计过程中帮助过的我老师和同学们，在同学们 的详细讲解和不断指导下，使得我完成此次课程设计的过程相对较轻松，特别 是周博同学，在自己还有课程设计需要完成时，抽出时间来，对我的整个程序 做出了详细的讲解，在程序运行的过程中，帮我查找错误，努力思考，对每个 错误存在的地方做了详细的讲解和更正。同样，也要感谢吴老师，在此次课程 设计过程中耐心的为我们的每个同学作辅导，关心每个同学的课程设计进程， 对我们的提问都给出了详细和让每个人满意的答复。 最后通过自己的努力和大家的帮助，在这次的课程设计中我懂得了 FSK 的 调制和解调原理，以及在 MATLAB 环境下运行的过程。对各种情况下的时域 和频域图形能做出合理的解释。理解了加入高斯白噪声后信号的变化。通过这 些，我才能完成以上的课程设计