基于matlab的FM调制

1、Matlab课程设计任务书学生姓名：专业班级：通信工程1403指导教师：韩一工作单位：信息工程学院题目：FM信号的仿真分析初始条件：理论方面电路分析基础、模拟电子、高频电子线路、通信原理、数字信号处理等专业知识设备方面微型计算机软件方面matlab仿真平台场地方面鉴主15楼通信实验室1要求完成的主要任务：调制信号：分别为300Hz正弦信号和三角波信号；载波频率：30kHz；解调方式：同步解调；要求：画出以下三种情况下调制信号、已调信号、解调信号的波形、频谱以及解调器输入和输出信噪比的关系曲线；1）调制指数=0.5；2）调制指数=1；3）调制指数=3。时间安排：2016.12.22至2016.1

2、.1仿真设计2016.1.2至2016.1.4写报告指导教师签名：系主任（或责任教师）签名2016年12月05日2016年12月05日目录摘要1Abstract21绪论32设计原理42.1MATLAB仿真软件42.2FM调制原理42.3FM解调原理62.3.1两种解调方式62.3.2解调过程分析72.4抗噪声性能分析73仿真实现93.1FM调制仿真93.1.1正弦波的情况93.1.2三角波的情况133.2 FM解调仿真163.3 抗噪声性能204结果分析214.1 仿真结果分析214.2 仿真设计的完整代码225总结23参考文献24武汉理工大学matlab原理与应用课程设计报告摘要频率调制广泛

3、应用于调频广播、电视伴音、微波通信、锁相电路和扫频仪等方面，在我们的日常生产生活中起着重要的作用。它是使载波的瞬时频率按照所需传递信号的变化规律而变化的调制方法。它是一种使受调波瞬时频率随调制信号而变的调制方法。由调频方法产生的无线电波叫调频波,其基本特征是载波的振荡幅度保持不变，振荡频率随调制信号而变。调频（FM）,就是高频载波的频率不是一个常数，是随调制信号而在一定范围内变化的调制方式，其幅值则是一个常数。本课程设计采用仿真软件matlab对不同信号进行FM调制，通过matlab软件对信号的FM调制进行分析，并分析其频谱特性以及不同调制系数下的调制增、人益。关键词：FM调制解调输出信噪比输

4、出信噪比28AbstractFrequencymodulationiswidelyusedinFMradio,TVsound,microwavecommunication,PLLandfrequencysweepmeteretc,playsanimportantroleinourdailylife.Itisamodulationmethodwhichmakestheinstantaneousfrequencyofthecarrierchangeaccordingtothechangingruleofthedesiredsignal.Itisamodulationmethodwhichmakes

5、themodulatedwaveinstantaneousfrequencychangewiththemodulationsignal.TheradiowavegeneratedbytheFMmethodiscalledFMwave,whosebasiccharacteristicisthattheamplitudeofthecarrieroscillationremainsunchanged,andtheoscillationfrequencychangeswiththemodulationsignal.Frequencymodulation(FM),thatis,thefrequencyo

6、fhigh-frequencycarrierisnotaconstant,withthemodulationsignalinacertainrangeofmodulation,theamplitudeisaconstant.ThiscourseisdesignedtousethesimulationsoftwareMATLABtocarryontheFMmodulationtothedifferentsignal,carriesontheanalysisthroughtheMATLABsoftwaretothesignalFMmodulation,andanalyzesitsspectrumc

7、haracteristicaswellasthemodulationgainunderdifferentmodulationcoefficient.Keywords:FM,modulation,demodulation,outputsignaltonoiseratio(SNR),outputsignaltonoiseratio(SNR)1绪论在通信系统中，从消息变换过来的原始信号所占的有效频带往往具有频率较低的频谱分量（例如语音信号），如果将这种信号直接在信道中进行传输，则会严重影响信息传送的有效性和可靠性，因此这种信号在许多信道中均是不适宜直接进行传输的。在通信系统的发射端通常需要有调制过程

8、，将调制信号的频谱搬移到所希望的位置上，使之转换成适于信道传输或便于信道多路复用的已调信号；而在接收端则需要有解调过程，以恢复原来有用的信号。调制解调方式常常决定了一个通信系统的性能。调制在通信系统中具有重要的作用。通过调制，不仅可以进行频谱搬移，把调制信号的频谱搬移到所希望的位置上，从而将调制信号转换成适合于信道传输或便于信道多路复用的已调信号，而且它对系统的传输有效性和传输可靠性有着很大的影响。调制方式往往决定了一个通信系统的性能。调制技术是指把基带信号变换成传输信号的技术。基带信号是原始的电信号，一般是指基本的信号波形，在数字通信中则指相应的电脉冲。在无线遥测遥控系统和无线电技术中调制就

9、是用基带信号控制高频载波的参数（振幅、频率和相位），使这些参数随基带信号变化。用来控制高频载波参数的基带信号称为调制信号。未调制的高频电振荡称为载波（可以是正弦波，也可以是非正弦波，如方波、脉冲序列等）。被调制信号调制过的高频电振荡称为已调波或已调信号。已调信号通过信道传送到接收端，在接收端经解调后恢复成原始基带信号。2 设计原理本课程设计采用matlab软件对正弦波和三角波信号进行FM仿真设计，并绘制相关波形及其频谱，对FM调制过程进行分析，同时通过输出信噪比和输入信噪比对解调系统的抗噪声性能进行分析。2.1MATLAB仿真软件MATLAB是美国MathWorks公司出品的商业数学软件，用于

10、算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。MATLAB是matrix和laboratory两个词的组合，意为矩阵工厂（矩阵实验室）。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际

11、科学计算软件的先进水平。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C,FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C+，JAVA的支持。MATLAB软件的优点在于

12、：（1）高效的数值计算及符号计算功能，能使用户从繁杂的数学运算分析中解脱出来；（2）具有完备的图形处理功能，实现计算结果和编程的可视化；（3）友好的用户界面及接近数学表达式的自然化语言，使学者易于学习和掌握；（4）功能丰富的应用工具箱（如信号处理工具箱、通信工具箱等），为用户提供了大量方便实用的处理工具。2.2FM调制原理连续波调制分为幅度调制，频率调制和相位调制。频率调整和相位调制都是使载波的相角发生变化，因此两者又统称为角度调制。调制在通信系统中具有十分重要的作用，通过调制，可对消息信号的频谱搬移，使已调信号适合信道传输的要求，同时也有利于实现信道复用。例如，将多路基带信号调制到不同的载频

13、上进行并行传输，实现信道的频分复用。调制方式往往对通信系统的性能有很大的影响。如果m(t)是连续信号，并且使某个参数连续地与m(t)相对应，称为模拟调制。大多数待传输的信号具有较低的频率成分。称之为基带信号，如果将基带信号直接传输，称为基带传输。但是，很多信道不适宜进行基带信号的传输，或者说，如果基带信号在其中传输，会产生很大的衰减和失真。因此，需要将基带信号进行调制，变换为适合信道传输的形式，调制是让基带信号m(t)去控制载波的某个(或某些)参数，使该参数按照信号m(t)的规律变化的过程。载波可以是正弦波，也可以是脉冲序列，以正弦信号作为载波的调制称连续波(CW)调制。对于连续波调制，已调信

14、号可表示为(221)S(t)=A(t)cosot+0(t)im它由振幅A(t)、角频率w和相位0(t)3个参数构成。控制3个参数中的任何一个都会实现调制，使之成为携带信息的信号。对于FM调制，图2.2.1FM调制模型其中，m(t)为基带调制信号，设调制信号为m(t)=Acos(2ft)(2.2.2)设正弦载波为c(t)=COSft)(2.2.3)信号传输信道为高斯白噪声信道，其功率为b2。在调制时，调制信号的频率去控制载波的频率的变化，载波的瞬时频偏随调制信号m(t)成正比例变化，即dt(2.2.4)式中，Kf为调频灵敏度(rad(sV)。这时相位偏移为申(t)=Kfim(t)dT则可得到调频

15、信号为sFJ()=Acosqt+Kpm(u)dr(226)2.3FM解调原理2.3.1两种解调方式调频信号的解调也分为相干解调和非相干解调。相干解调仅适用于NBFM信号，而非相干解调对NBFM信号和WBFM信号均适用。1)非相干解调2非相干解调器由限幅器、鉴频器和低通滤波器等组成，其方框图如图所示。限幅器输入为已调频信号和噪声，限幅器是为了消除接收信号在幅度上可能出现的畸变；带通滤波器的作用是用来限制带外噪声，使调频信号顺利通过。鉴频器中的微分器把调频信号变成调幅调频波，然后由包络检波器检出包络，最后通过低通滤波器取出调制信号。鉴频器图2.3.1.1非相干解调框图2)相干解调由于NBFM信号可

16、分解成同相分量与正交分量之和，因而可以采用相干解调法进行解调，如图sNBFM(t)si(t)sp(t)sd(t)c(t)图2.3.1.2相干解调框图相干解调可以恢复原始调制信号。但要求本地载波与调制载波同步，否则将使解调信号失真。2.3.2解调过程分析本课程设计采用同步解调方式。输入调频信号为相干载波为uFM(t)=Ucost+kJ'm(t)dt)f0(2.3.2.1)c(t)=cos(2*pi*fc*t)(2.3.2.2)乘法器的作用是把调频信号变成有多种频率的波的混合,s(t)=-sin2wt+KJm(t)dt丄(1-cos2®t)p2c2fc乘法器输出为(2.3.2.3

17、)经低通滤波器后取出器低频分量为s(t)=KJm(t)dtd2f(2.3.2.4)在经过微分器，即得出解调出的基带信号：m(t)=Km(t)(2.3.2.5)02f相干解调可以恢复出原来的基带信号，而且要求本地载波与调制载波同步，否则会使解调信号失真。2.4抗噪声性能分析设调频信号为3SFM(t)=Acos(3ct+KJtm(T)dT)g故其输入功率为输入噪声功率为A2(2.4.2)(2.4.3)Ni=noBFM因此输入信噪比为SiiNiA2(2.4.4)在大信噪比条件下，信号和噪声的相互作用可以忽略，这时可以把信号和噪声分开来算，这里，我们可以得到解调器的输出信噪比So=N_o3A2K2m2

18、(t)f8兀2nf3om(2.4.5)上式中，A为载波的振幅，Kf为调频器灵敏度，fm为调制信号m(t)的最高频率，no为噪声单边功率谱密度。我们如若考虑m(t)为单一频率余弦波时的情况，可得到解调器的制度增益为FM3=m22fnofm考虑在宽带调频时，信号带宽为BFM=2(mf+1)fm=2(纣+fm)(2.4.6)(2.4.7)则可以得到GFM=3m2f(mf+1)(2.4.8)可以看出，大信噪比时宽带调频系统的信噪比增益是很高的，它与调频指数的立方成正比。可见，加大调频指数/，可使调频系统的抗噪声性能迅速改善。3 仿真实现本课程设计采用的是MATLABR2016b版本。根据第2章的设计原

19、理，使用matlab编写相应的M文件实现FM的调制、解调以及抗噪声性能的分析。3.1 FM调制仿真打开matlab软件，新建脚本，编写M文件，根据设计要求，调制信号：分别为300Hz正弦信号和三角波信号；载波频率：30kHz；调制仿真的代码如下4：t0=0;t1=0.01;%调制信号的时域范围fs=100000;%抽样频率ts=1/fs;%采样率t=t0:ts:t1;fc=30000;%设定载波频率x=cos(2*pi*fc*t);%生成载波fm=300;%设定调制信号频率%mf=0.5;%mf=3;mf=1;m=cos(2*pi*fm*t);%生成调制信号%m=sawtooth(2*pi*3

20、00*t,0.5);s=modulate(m,fc,fs,'fm',(fc/fs)*2*pi/mf);%生成已调信号3.1.1 正弦波的情况原信号波形如图所示jFigure-1X文件(F)帚竭至着阿SfcACI工具E卓面Q)画口(W)帮助(H)a己£|矗鸟轨劭題原7-1|隹|口图3.1.1.1原调制信号已调信号波形如图所示&Figure1口X文件血'SSS(E)至看W插入工具CD宝®Dj国口CW)翻助0H)D厲Z®底饨飞r?暂*一x-罔口口廿tXOCMO0020.003-口0040.仃0忌口DOS0.007OnDBD口口号0.01图

21、3.1.1.2已调信号(mf=0.5)图3.1.1.3mf=1NFigure1口-1査看“描Ml!工具E墓面DI窗口CW)群助H1a1$|%气7從X7昼1国1B:O0.0010.0020.003-0.0040.0050.0060.007O.-0-0B0.0090.01图3.1.1.4mf=3已调信号的频谱如图所示Figure2-文件蕖垢(E)查春VJ插入工具E其面叨裔口的帮助旧】已巳©|恰|整凤紗渥摇彳凰|口固|口图3.1.1.5mf=0.5*Figure2圭伴(F)扁揖(E)苣吉M捆入H目E卓面0=0帮ESCH)nXU芒PIX|傀f<7倉跟晟”|良|目|OO123156709

22、1<Jx1O4图3.1.1.6mf=1匡Figure2IIX文件旧扁辑E章走（VIWACI工具E点JHTQ窗口CW辖取HD£3苟|AA汇/”|風|辰图3.1.1.7已调信号的频谱（mf=3）3.1.2三角波的情况原调制信号波形如图所示:養iFigure1立件(F)輪審=(巳描入工具点m歯口网哥則巴乂1Cl£3*3|%.C121口困|O图3.121原三角波信号1*|Figliu1口X文件阳為無旧i插入工具e至面cozi窗口cw）縉目HHDP4|4JH<7显+凰固na图3.1.2.2已调信号（mf=0.5）O0.0010.0O.OU3O.UCMt».OO5

23、O.U1360.0070.0DE?HFigure1立件(F)蕖垢EJ宜吉M折入H晨E点面esn<wD抹时H口X口心出&|冷I头铁r?®區虫十|型|口£=l|口图3.1.2.3mf=1060.2-o_e唧刖啊测呷-0.244-O-SM血山lit恤-1O0.0010.0020.0030.0040.0050.00600070.-0080.0090.01Figure1文件為旨Mm扬人T-am克面Q)裔口他耕助讯)u芒ta%企寳®电必*目a图3.124mf=3已调信号的频谱如图所示Figure?ZX立件舸領(E)至舌M病人(0工SE贞面Q)圜口(W)辛助(H)

24、U芒Hdl怪I代純區目，|曰罔I口DSBx1O4图3.1.2.5mf=0.5鼻jFigure-2丈?(牛(F)辐银1莊入工且E窗口CW)括国KH11dgre|【符|虫头O®膜必W丄口目1OO1J2345Ta91O图3126mf=1鼻jFiguire-2口X立件（円sb（e）se（.v）j®入工且Era（D）豳口wvjweo（hjU曰Ud|0|举H虫彰®归區|鈕|口罔|口图3.1.2.7mf=33.2FM解调仿真正弦波解调波形及其频谱如图所示Figure3口X|文件编唱童音M牺入工具E宾面1D裔口CW3群助H勺iU心材騒|險|九注凸®X”|岂|£

25、|旦*Figure4X文件旧輛E）查看（V）插入I）工具E宾面（D窗口刪冇旨甘国k%电紳勒從掐f反冋图3.2.1解调信号波形及频谱mf=0.5Figure3X袁:件扁揖査看M®A(I)工具E卓面(6誉口(W)稻助(H)ac3氐丨%紳汇家”|孱|直SFigure-4口X立件(F)漏35(曰®WtV5SACD工貝CD宾直(Q窗口闪和朋(H)勺U呂恪I汽虫紳®礙必十|凰|口目|型图3.2.2mf=1兰Figure3X交祥旧輪辑查看插入工具E臺面(D啬口WJ拜助(巴巴£Ife|：o®vx-1a|010*Figure4X宏件(F)需辑i：曰麦着M插入工具

26、CD壳面(巧窗口(W)腳助(H)x104口t±耳騎怎怯匚紳®也搖J毘T刃丁旦35030025020015010050O23图323mf=33.3抗噪声性能调制系数分别为0.5、1、3时解调器输入和输出信噪比的关系曲线如图所示莖Figure4X直件旧褊辑(E)查吾插入I)工具克面(叨窗口阿帮馭H吾吾fe:%図鋼L貳圍同25002D001500100050000510152025输入信噪比图3.3.1输出信噪比与输入信噪比曲线4 结果分析4.1仿真结果分析从第三章的仿真结果可以得出以下结论：从FM的调制仿真绘制的波形可以看出，无论是正弦波还是三角波，已调信号的频率会随着原信号的

27、幅度大小而变。而这种频率随原信号大小改变的程度会随着调制系数的不同而出现差异。所以，在调制系数不同的条件下，已调信号以及已调信号的频谱都会不同。从FM的解调仿真绘制的波形可以看出，调制系数为0.5的时候，解调输出的信号有很大程度的失真，原因是调制系数太小，解调器的抗噪声性能太差，导致解调输出的波形失真，随着调制系数的增大，调制系数为1时，解调输出的波形基本上与原调制信号相同。解调信号的频谱也很清晰。从不同调制系数下的解调器输出信噪比与输入信噪比的曲线可以看出，随着输入信噪比的增大，输出信噪比也会增大，同时调制系数增大时，输出信噪比随输入信噪比增大的程度会更明显，对相同的输入信噪比，输出信噪比会

28、更大，这时系统的抗噪声性能就越强。4.2仿真设计的完整代码t0=0;t1=0.01;%调制信号的时域范围fs=100000;%抽样频率ts=1/fs;%采样率t=t0:ts:t1;fc=30000;%设定载波频率x=cos(2*pi*fc*t);%生成载波fm=300;%设定调制信号频率%mf=0.5;%mf=3;mf=1;m=cos(2*pi*fm*t);%生成调制信号(正弦信号)%m=sawtooth(2*pi*300*t,0.5);%三角波信号s=modulate(m,fc,fs,'fm',(fc/fs)*2*pi/mf);%生成已调信号figure(1);plot(t,

29、s,t,m);df=fs/(length(t)-1);f=(0:length(t)-1)*df;figure(2);plot(f,abs(fftshift(fft(s);%已调信号的频谱y=demod(s,fc,fs,'fm');figure(3);plot(t,y);figure(4);plot(f,abs(fftshift(fft(y);%解调信号的频谱x3=(0.1:0.05:22.05)'mf=0.513;y3(:,1:3)=3*mf.A2.*(1+mf).*x3;figure(5);plot(x3,y3(:,1),'b');holdonplot(x3,y3(:,2),'r');plot(x3,y3(:,3),'g');%不同调制系数下的输出信噪比与输入信噪比的曲线legend('mf=0.5','mf=1','mf=3');xlabel(输入信噪比')；ylabel(输出信噪比');