FM调制解调系统设计与仿真linWord文档下载推荐.docx

1、2.3 FM调制仿真结果 （62.4 FM解调模型的建立 （62.5 解调过程分析 （72.6高斯白噪声信道特性 （82.7信噪比分析 （92.8调频系统的抗噪声性能分析 （10三.仿真实现 （123.1 MATLAB源代码 （123.2MATLAB仿真结果及分析 （12四.心得体会 （14五.参考文献 （14摘要调制在通信系统中有十分重要的作用。通过调制,不仅可以进行频谱搬移,把调制信号的频谱搬移到所希望的位置上,从而将调制信号转换成适合于传播的已调信号,而且它对系统的传输有效性和传输的可靠性有着很大的影响,调制方式往往决定了一个通信系统的性能。FM信号的调制属于频谱的非线性搬移,它的解调也

2、有相干和非相干解调两种方式。本课程设计使用的仿真软件为MATLAB,利用MATLAB集成环境下的M文件,编写程序来实现FM信号的仿真分析,并分别绘制出基带信号、载波信号、已调信号的时域波形和频域波形;再进一步分别绘制出对已调信号叠加噪声后信号、同步解调前信号和解调后基带信号的时域波形;最后绘出FM基带信号调制和解调系统后的输入输出信噪比的关系,并通过与理论结果波形对比来分析该仿真调制与解调过程的正确性及噪声对FM信号解调的影响。在课程设计中,系统开发平台为Windows7,在该平台运行MATLAB程序完成了对FM调制和解调以及对叠加噪声后解调结果的观察和分析。关键词:FM 调制解调MATLAB

3、 噪声AbstractModulation has a very important role in the communication system . Modulation, can be used to not only move the spectrum , the spectrum of the modulated signal is moved to a desired position so as to convert the modulated signal suitable for transmission of modulated signals,but also it s

4、ignificantly affect the validity of the transmission system and the transmission reliability . Often determines the performance of the modulation scheme of a communication system . FM modulated signal is Non-linear movement of spectrum . it also has the coherent demodulation and non-coherent demodul

5、ation of two ways.Simulation Software This course is designed for use as MATLAB, the use of M-files MATLAB integrated environment , write a program to achieve the simulation of FM signals , and were drawn out of the baseband signal , the carrier signal, modulated signal domain waveform and frequency

6、 domain waveforms ; were drawn further modulated signal superimposed on the noise signal, a synchronization signal and demodulates the first time domain waveform of the demodulated baseband signal ; final output SNR plot FM input baseband signal after the modulation and demodulation system of relati

7、ons and, through comparison with the theoretical results to analyze waveforms and noise affect the correctness of the simulation process of modulation and demodulation of the FM signal demodulation . In the course design , system development platform for Windows7, the platform running MATLAB program

8、 completed and FM modulation and demodulation noise superimposed on the demodulation results after observation and analysis .Keywords : FM Modulation Demodulation MATLAB Noise一.通信系统简介通信的目的是传输信息。通信系统的作用就是将信息从信息源发送到一个或多个目的地。对于任何一个通信系统,均可视为由发送端、信道和接收端三大部分组成（如图1所示。图1 通信系统一般模型信息源（简称信源的作用是把各种信息转换成原始信号。发送设

9、备的作用产生适合传输的信号,即使发送信号的特性和信道特性相匹配,具有抗噪声的能力,并且具有足够的功率满足原距离传输的需求。信息源和发送设备统称为发送端。发送端将信息直接转换得到的较低频率的原始电信号称为基带信号。通常基带信号不宜直接在信道中传输。因此,在通信系统的发送端需将基带信号的频谱搬移（调制到适合信道传输的频率范围内进行传输。这就是调制的过程。信号通过信道传输后,具有将信号放大和反变换功能的接收端将已调制的信号搬移（解调到原来的频率范围,这就是解调的过程。信号在信道中传输的过程总会受到噪声的干扰,通信系统中没有传输信号时也有噪声,噪声永远存在于通信系统中。噪声对于信号的传输是有害的,它能

10、使模拟信号失真。在本仿真的过程中我们假设信道为高斯白噪声信道。调制在通信系统中具有十分重要的作用。一方面,通过调制可以把基带信号的频谱搬移到所希望的位置上去,从而将调制信号转换成适合于信道传输或便于信道多路复用的已调信号。另一方面,通过调制可以提高信号通过信道传输时的抗干扰能力,同时,它还和传输效率有关。具体地讲,不同的调制方式产生的已调信号的带宽不同,因此调制影响传输带宽的利用率。可见,调制方式往往决定一个通信系统的性能。在本仿真的过程中我们选择用调频调制方法进行调制。信息源 发送设备 信 道接受设备 信息源 噪声源 发送端接收端 信道调制过程是一个频谱搬移的过程,它是将低频信号的频谱搬移到

11、载频位置。而解调是将位于载频的信号频谱再搬回来,并且不失真地恢复出原始基带信号。在本仿真的过程中我们选择用相干解调方法进行解调。角度调制及解调电路不同于频谱线性搬移电路。它是用低频信号去调制高频振荡的相角,或是从已调波中解出调制信号所进行的频谱变换,这种变换不是线性变换,而是非线性变换。因此,我们把角度调制及调角波的解调电路称为频谱非线性变换电路。FM在通信系统中的使用非常广泛。FM广泛应用于高保真音乐广播、电视伴音信号的传输、卫星通信和蜂窝电话系统等。本次课设中使用功能强大的MATLAB仿真软件对FM信号进行仿真分析。6max（K dt t m f 二.FM 调制解调系统设计2.1 MATB

12、LAB 简介MATLAB 的名称源自Matrix Laboratory ,它的首创者是在数值线性代数领域颇有影响的Cleve Moler 博士,他也是生产经营MATLAB 产品的美国Mathworks 公司的创始人之一。MATLAB 是一种科学计算软件,专门以矩阵的形式处理数据。MATLAB 将高性能的数值计算和可视化集成在一起,并提供了大量的内置函数,从而使其被广泛地应用于科学计算、控制系统、信息处理等领域的分析、仿真和设计工作中。MATLAB 与其它计算机高级语言如C ,C+等相比,MATLAB 语言编程要简洁得多,编程语句更加接近数学描述,可读性好,其强大的图形功能和可视化数据处理能力也

13、是其它高级语言望尘莫及的。对于具有任何一门高级语言基础的学生来说,学习MATLAB 十分容易。MATLAB 使得人们摆脱了常规计算机编程的繁琐,让人们能够将大部分精力投入到研究问题的数学建模上。可以说,应用MATLAB 这一数学计算和系统仿真的强大工具,可以使科学研究的效率得以成百倍的提高。2.2 FM 调制模型的建立从频率调制的相位与频率关系可以看出,调频信号可通过直接调频和间接调频两种方法得到,所谓间接调频就是先对调制信号积分再调相而得到。同样,调相信号也可以通过直接调相和间接调相两种方法得到,间接调相就是先对调制信号进行微分再进行频率调制。根据调制后已调信号的瞬时相位偏移的大小,可将角度

14、调制分为宽带调制（宽带调频和宽带调相和窄带调制（窄带调频和窄带调相。如果调频信号或调相信号的最大瞬时相位偏移保持在很小的范围内,一般小于30即满足条件:时,则称为窄带调频或窄带调相。当上述条件不满足时,就称为宽带调频或宽带调相。 图2 FM 间接调制模型其中,m（t为基带调制信号,设调制信号为m（t=cos（2*pi*fm*t设正弦载波为c（t=cos（2*pi*fc*t在调制时,调制信号的频率去控制载波的频率的变化,载波的瞬时频偏随调制信号m（t成正比例变化,即:dt t m k t t tf +=00（式中,f K 为调频灵敏度。 这时相位偏移为f （00+=则可得到调频信号为（cos（u

15、 00FM dt t m k t U t tf c +=FM 的频谱的计算:可以看出FM 的频谱与（f n m J 的值有关。mf m =1204202cos （2（4cos （22cos （2（sin cos（+=+=+=n f n f f f f f tn m J m J t m J t m J m J t m12（53112sin（25sin （23sin （2sin （2sin sin（+=+=+=n f n f f f f tn m J t m J t m J t m J t mtn m J A t V k t A t c n f n f c FM （cos （sin cos （U

16、_+=+=+=其信号带宽为 +=1（2f m B 调制信号产生的部分M 文件如下:fc=30000; %载波频率 ct=cos（2*pi*fc*t; %载波信号 am=1; %调制信号幅度 fm=300; %调制信号频率 mt_sin=am*cos（2*pi*fm*t; %正弦波 mt_tri=sawtooth（2*pi*fm*t,0.5; %三角波 mf=1; %调制指数,可调int_mt_sin（1=0; %正弦波积分 for i=1:length（t-1int_mt_sin（i+1=int_mt_sin（i+mt_sin（i*dt; endint_mt_tri（1=0;for i=1:l

17、ength（t-1 %三角波积分 int_mt_tri（i+1=int_mt_tri（i+mt_tri（i*dt;sfm_sin=am*cos（2*pi*fc*t+mf*2*pi*fm*int_mt_sin; sfm_tri=am*cos（2*pi*fc*t+mf*2*pi*fm*int_mt_tri;2.3 FM 调制仿真结果图3 FM 调制2.4 FM 解调模型的建立由前述公式知道输入调频信号为相干载波为乘法器的作用是把调频信号变成有多种频率的波的混合,乘法器输出为2c o s 1（22s i n 21（t dt t m K t t s c f c p -*+-= 经低通滤波器后取出器低频

18、分量为dt t m K t s f d =（2（ 在经过微分器,即得出解调出的基带信号:（0t m K t m f =相干解调可以恢复出原来的基带信号,而且要求本地载波与调制载波同步,否则会使解调信号失真。2.5 解调过程分析在这次的解调设计中,采用MATLAB 自带的解调函数demod 对已调信号进行解调。该函数使用形式如下所示: Z = demod（Y ,Fc,Fs,Method其中,Y 为输入的调频波函数,Fc 、Fs 分别为载波频率和抽样频率,而且要满足Fc 2Fs ,Method 为调制方式。FM 解调的M 文件部分代码如下: y_sin=demod（sfm_sin,fc,Fs,fm

19、; y_tri=demod（sfm_tri,fc,Fs, yyn_sin=demod（yn_sin,fc,Fs, yyn_tri=demod（yn_tri,fc,Fs,通过M 文件绘制出大、小信噪比解调的输出波形如下:图4小信噪比时的FM 解调图5 大信噪比时的FM 解调2.6 高斯白噪声信道特性设正弦波通过加性高斯白噪声信道后的信号为（cos（c r t A t n t =+其中,白噪声（n t 的取值的概率分布服从高斯分布。MATLAB 本身自带了标准高斯分布的内部函数awgn 。awgn 函数产生的随机序列服从均值为0m =,方差21=的高斯分布。正弦波通过加性高斯白噪声信道后的信号为故

20、其有用信号功率为2A S =噪声功率为N =信噪比S N 满足公式1010log （S B N则可得到公式2210210B A =我们可以通过这个公式方便的设置高斯白噪声的方差。2.7信噪比分析在本仿真过程中,选择了1,3,5三种不同调制指数以示区别,并绘制输入输出信噪比关系曲线,如下图6所示。mf=1mf=3mf=5图6 正弦波输入输出信噪比2.8 调频系统的抗噪声性能分析从前面的分析可知,调频信号的解调有相干解调和非相干解调两种。相干解调仅适用于窄带调频信号,且需同步信号;而非相干解调适用于窄带和宽带调频信号,而且不需同步信号,因而是FM 系统的主要解调方式,所以这里仅仅讨论非相干解调系统

21、的抗噪声性能,其分析模型如图所示。图7 调频系统抗噪声性能分析模型图中带通滤波器的作用是抑制信号带宽以外的噪声。（n t 是均值为零,单边功率谱密度为0n 的高斯白噪声,经过带通滤波器后变为窄带高斯噪声（i n t 。限幅器是为了消除接收信号在幅度上可能出现的畸变。设调频信号为（cos（tFM c f S t A t K m d -=+故其输入功率为i A S =输入噪声功率为i o FM N n B =因此输入信噪比为2i i FMS A N B = 在大信噪比条件下,信号和噪声的相互作用可以忽略,这时可以把信号和噪声分开来算,这里,我们可以得到解调器的输出信噪比武汉理工大学MATLAB 课

22、程设计报告222233（8ofoo mA K m t S N n f =上式中,A 为载波的振幅,f K 为调频器灵敏度,m f 为调制信号（m t 的最高频率,o n 为噪声单边功率谱密度。我们如若考虑（m t 为单一频率余弦波时的情况,可得到解调器的制度增益为232A oFMi o mi S N G m S n f N =考虑在宽带调频时,信号带宽为2（12（FM f m m B m f f f =+则可以得到23（1FM f f G m m =+可以看出,大信噪比时宽带调频系统的信噪比增益是很高的,它与调频指数的立方成正比。可见,加大调频指数f m ,可使调频系统的抗噪声性能迅速改善。三

23、.仿真实现本次设计是在MATLAB 平台上,在前述FM 通信系统的分析基础上,进行程序代码的编写,从而得到相关仿真波形,由此进行分析讨论。图8 程序流程图3.1 MATLAB 源代码根据设计要求所编写的源代码已在在附页中给出,详情请看附页I 。3.2MATLAB 仿真结果及分析在MATLAB 平台上运行程序得到各仿真图,由于篇幅有限,已在附页II 给出。这里针对部分仿真结果进行分析。（1系统生成满足条件的信号波形,包括正弦波和三角波。（2在未加入高斯白噪声时,调制系统输出理想的FM 已调信号,已调信号的频率随着原信号的幅度大小而变。解调系统也能理想地解调出原始输入信号,解调后的波形光滑,与FM

24、 调制程序添加高斯白噪声程序FM 解调程序绘制频谱图程序显示程序初始化程序原信号一致。（3在加入较弱的高斯白噪声时,调制系统输出较为理想的FM已调信号,已调信号的频率随着原信号的幅度大小而变但并不是理想的。解调系统也能正确地解调出原始输入信号,解调后的波形稍带毛刺,与原信号基本一致。（4在加入较强的高斯白噪声时,调制系统输出不太理想的FM已调信号,已调信号的频率与原信号的幅度大小关系不明显。解调后的波形有失真现象。（5输入输出信噪比呈正向线性关系,输入信噪比越大,输出信噪比也越大。总述:FM调制解调系统设计成功。结果表明:信道噪声很大程度上影响着调制解调系统性能,实际运用中应尽量降低信道噪声;

25、为了少失真地得到原始信号,在提高系统性能的基础上,还应增大输入信噪比。四.心得体会课程设计是培养学生综合运用所学知识,发现,提出,分析和解决实际问题,锻炼实践能力的重要环节,是对学生实际工作能力的具体训练和考察过程. 回顾起此次通信原理课程设计,至今我仍感慨颇多,的确,从选题到定稿,从理论到实践,在整整一星期的日子里,可以说得是苦多于甜,但是可以学到很多很多的的东西,同时不仅可以巩固了以前所学过的知识,而且学到了很多在书本上所没有学到过的知识。通过这次课程设计使我懂得了理论与实际相结合是很重要的,只有理论知识是远远不够的,只有把所学的理论知识与实践相结合起来,从理论中得出结论,才能真正为社会服

26、务,从而提高自己的实际动手能力和独立思考的能力。在设计的过程中遇到问题,可以说得是困难重重,这毕竟第一次做的,难免会遇到过各种各样的问题,同时在设计的过程中发现了自己的不足之处,对以前所学过的知识理解得不够深刻,掌握得不够牢固通过这次课程设计之后,一定把以前所学过的知识重新温故。这次课程设计终于顺利完成了,在设计中遇到了很多问题,最后在不懈的努力下终于迎刃而解。感谢这次旅程,我又成长不少!五.参考文献1通信原理（第六版樊昌信等。北京:国防工业出版社。2MATLAB7.0在数字信号处理中的应用罗军辉等。机械工业出版社。3MATLAB程序设计教程刘卫国等。中国水利水电出版社。附页I仿真波形图附页I

27、I源代码%*本程序主要实现FM调制与解调系统设计* %*以正弦波和三角波进行设计*%*并考虑信道中高斯白噪声的影响*%*初始化*dt=0.00001; %积分步长N=1000; %采样点数Fs=1/dt; %抽样频率Fs=100000hz,满足抽样定理n=0:N-1;t=n/Fs;%*信号设定* %载波频率ct=cos（2*pi*fc*t; %载波信号am=1; %调制信号幅度fm=300; %调制信号频率mt_sin=am*cos（2*pi*fm*t; %正弦波mt_tri=sawtooth（2*pi*fm*t,0.5; %三角波mf=3; %调制指数%*调制模块* %正弦波积分endlength（t-1 %三角波积分int_mt_tri（i+1=int_mt_tri（i+mt_tri（i*dt;sfm_tri=am*cos（2*pi*fc*t+mf*2*pi*fm*int_mt_tri;%*