基于matlab的FM通信系统仿真设计与实现课程设计报告

1、成绩沮南呻兔学院课程设计报告目：基于matlab的FM通信系统仿真设计与实现学生姓名：杨丽君学生学号：1008030317系别：电气信息工程学院专业：电子信息工程届别：14-指导教师：马立宪电气信息工程学院制2021年5月基于matlab的FM通信系统仿真设计与实现学生：杨丽君指导教师：马立宪电气信息工程学院电子信息工程i课程设计的任务与要求i.i课程设计的任务(1)熟悉MATLAB文件中M文件的使用方法,包括函数、原理和方法的应用.(2)加深对FM信号调制原理的理解.(3)画出基于MATLAB的FM通信系统仿真设计与实现设计的原理图.1.2 课程设计的要求(1)学会MATLAB软件的安装.(

2、2)在做完FM调制仿真之后,在今后遇到类似的问题,学会对所面对的问题进行系统的分析,并能从多个层面进行比拟.(3)熟练并且掌握对MATLAB软件的使用,学会输入程序并且加以运行.1.3 课程设计的研究根底通信的目的是传输信号.通信系统的作用是将信息从信息源发送到一个或者多个目的地.模拟信号是时间和幅值上都连续的信号.调制是用原始信号即调制信号去限制高频载波信号的某一参数,是指随着原始信号幅度的变化而变化.而FM频率调制是高频载波信号的频率随着原始信号幅度变化而变化.解调是将已调制的信号恢复成原始信号即基带调制信号.以下是通信系统的一般模型：图1通信系统的一般模型(1)信息源信息源(简称信源)的

3、作用是把各种消息转换成电信号.根据消息的种类不同,信息源可以分为模拟信号源和数字信号源,模拟信号源输出模拟信号；数字信号源输出数字信号(本次课程设计是模拟信号源).(2)发送设备发送设备的作用是产生适合在信道中传输的信号即使发送信号的特性与信道特性相匹配,具有抗信道干扰水平,并且具有足够的功率以满足远距离传输的需要.因此,发送设备涵盖的内容很多,包含变换、放大、滤波、编码、调制的过程.(3)信道信道是一种物理媒介,用于将来自发送设备的信号发送到接收端.信道分为无线信道和有线信道.在有线信道中可以是明线、电缆、光纤.在无线信道中,信道可以是自由空间.信道的固有特性及引入的干扰与噪声直接关系到通信

4、的质量.(4)噪声源噪声源是信道中及分布在系统中的其他各处噪声集中表示,噪声是随机的、形式是多样的,它的出现直接干扰信号的传输.(5)接收设备接收设备的功能是将信号放大及反变换(如译码、解调等).目的是从受到减损的信号中正确恢复出原始信号,减少在传输过程中噪声与干扰所带来的影响.2 FM通信系统方案制定2.1 方案提出下列图为大体模拟通信系统模型：图2模拟通信系统模型本次课程设计FM模拟通信系统模型中对于调制信号通过调制器产生的调频波有两种方法其一为直接调制即宽带调制,其二为间接调制及窄带调制.解调器对应的解调方法也有两种其一为相干解调另外一种为非相干解调.所以据此我提出了四种方案：方案一：模

5、型中调制器中的调制方法为直接调制即宽带调制,解调器对应的解调方法为相干解调.方案二：模型中调制器中的调制方法为直接调制即宽带调制,解调器对应的解调方法为非相干解调.方案三：模型中调制器中的调制方法为间接调制即窄带调制,解调器对应的解调方法为非相干解调方案四：模型中调制器中的调制方法为间接调制即宽带调制,解调器对应的解调方法为相干解调.2.2 方案论证窄带调频的应用更广泛与宽带调频,我们对此也更为熟悉,技术也更为成熟.止匕外,它的最大频率偏移较小,占据带宽较窄、抗干扰性能更好等.所以本次课程设计选择窄带调频.一般情况下,相干解调法较适用于窄带调频.所以在以上选择的前提下,本次课程设计的解调方法选

6、择相干解调针对以上的分析选择的最正确方案为方案四.3 FM通信系统方案设计3.1 FM通信系统模型设计m(t)模拟信源一>调制器>发滤波器一信道I收滤波器I1|ni(t)m0(t)(s/n)o噪声n(t)解调器信宿no(t)(s/n)i图3模拟通信系统模型设计3.2 FM通信系统各局部的功能调制器：使信号与信道相匹配,有利于信号在信道中传输.发滤波器：滤除调制器输出的无用信号.收滤波器：滤除信号频带以外的噪声.一般设N(t)为高斯白噪声,那么Ni(t)为窄带白噪声.在通信系统中一般需要将信号进行相应调制,以利于信号在信道上的传输,调制是将用原始信号去限制高频振荡信号的某一参数,使之

7、随原始信号的变化而成规律变化.调制可分为线性调制和非线性调制.线性调制有AM、DSB等,非线性调制有FM、PM等,这里主要讨论FM调制通信系统3.3 FM通信系统参数的计算及原理(1) FM调制原理角调制不是线性调制,角调制中已调信号和调制信号频谱之间不是线性关系而是产生出新的与频谱搬移不同的新的频率分量,呈现非线性特性,故又成为非线性调制.FM调制中瞬时角频率是关于调制信号的线性函数,瞬时角频率偏移量Aw=KFMf那么,瞬时角频率为：w=wc+KFMf(t);KFM为频偏指数,那么调频信号为:Sfm(t)=cos(wt+KFMff(t)dt;当调制信号是单频余弦时,调制信号为:Sfm(t)=

8、Acos(wctKFMAmcoswmtfdt)=Acos(wct：FMsinwmt);PFM为调频指数,PFM=KFMAm/Wm;调制信号的信号的产生m(t)FMSfm(t)图4调频信号的直接产生m(t)(0dtkPMSfm(t)图5调频信号的间接产生如上两图,在间接产生方法中,由于实际调相的相位调节器的范围在(-冗),故而间接调频只能用于相位偏移和频率偏移不大的窄带调制情况,而直接调频那么常用于宽带调制情况.窄带角频率调制时,最大瞬时相位偏移远小于30°即KFMff(tdt)max«n/6,调频信号的时域表达式Sfm(t)=cos(wj+KFM/出),当为窄带调制时,Sf

9、m(t)=Acoswct-AFMJf(t)dtsinwd,假设调制信号的频谱为F(w),而且假设f(t)的平均值为零,由傅氏变换可得：Snbfm(w)=A二(w-wc)二(wwc)AKFM/2F(w-wc)/(w-wc)-F(wwc)/(wwc),此时调频信号的带宽为调制信号的两倍.其相应图形如下：wm“m(w)wm图6调制信号频谱非相干解调由于调频信号的瞬时频率正比于调制信号的幅度,因而调制信号的调节器必须产生正比于输入频率的输出电压,即输入为调制信号为：S(t)=AcosWct+KFMf(t)dt,那么解调器的输出应为：SoOAKfmf最简单的解调器是鉴频器,原理图如下：图8非相干解调原理

10、图包络检波输出为S0(t)=KFMkd,kd为鉴频器的灵敏度,微分器和包络检波构成鉴频器相干解调由于窄带调频信号可分解为同相分量和正交分量,因而可以用线性调制中的相干解调来进行解调.原理图如下:图9相干解调原理图带通滤波器的作用是抑制信号带宽以外的噪声,低通滤波器输出So(t)=AKFMf(t)/2；相干解调不仅要恢复基带调制信号,还要恢复相干载波,上图所示的相干解调法只适用于窄带调频.(3)调频系统的抗噪声性能主要讨论非相干解调的抗噪声性能,输入信号Sfm=cos(wct+KFMfff(t)dt),输入功率S=A2/2,输入噪声Ni=n°BFM,信噪比为S/Ni=A*A/2n

11、76;BFM,非相干解调的解调器的输入端是调频信号与噪声的叠加,为S(t)+n(t)=SFM(t)+n(t),在大信噪比情况下,解调器的输出端的信噪比为So/NO=3A2Kf2m2(t)/8n2nofm3,考虑单频余弦信号调制,故可得大信噪比情况下的信噪比增益为Gfm=(&/N°)/(Si+Ni)=(3/2)(mf2)(A2/2)/(nofo),单频时,带宽Bfm=2(饵+*)所以增益可化解为Gfm=3PFM(PFM+1).由此可看出,性噪比增益和调频指数的三次方成正比.加大调频指数Bfm,可使调频系统的抗噪声性能迅速改善.4 FM通信系统的仿真4.1 FM通信系统仿真软件M

12、ATLAB的介绍MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境.它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学计算软件的先进水平.MATLAB和mathematicamaple并称为三大数学软件.它以矩阵为根本数据单位,在数学类科技应用软件中在数值计算方面首屈一指.MATLAB可以进行矩阵

13、运算、绘制函数和数据、实现算法、创立用户界面、连接其他编程语言的程序等,主要应用于工程计算、限制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域4.2 FM通信系统仿真实现图10程序流程图4.3 MATLAB源代码程序如下：echooffcloseallclearallclcdt=0.001;t=0:dt:1.5;am=5;fm=5;mt=am*cos(2*pi*fm*t);figure(1)plot(t,mt);xlabel('t');title('调制信号时域图');fc=50;ct=cos(2*pi*fc*t);figure(2)plot

14、(t,ct);xlabel('t')title('载波信号时域图');fs=1000;sfm=modulate(mt,fc,fs,'FM');ts=0.001;fs=1/ts;df=0.25;mt=am*cos(2*pi*fm*t);fs=1/ts;ifnargin=2n1=0;elsen1=fs/df;endn2=length(mt);n=2A(max(nextpow2(n1),nextpow2(n2);M=fft(mt,n);mt=mt,zeros(1,n-n2);df1=fs/n;M=M/fs;f=0:df1:df1*(length(mt)

15、-1)-fs/2;figure(11)plot(f,abs(fftshift(M)xlabel('频率f)title('调制信号的频谱图');figure(3)plot(t,sfm);xlabel('t')title('调频信号的时域图');axis(00.15-11);b=fft(sfm,1024);f=(0:length(b)-1)*fs/length(b)-fs/2;figure(4)plot(f,abs(fftshift(b);xlabel('Frequence(Hz)');ylabel('PowerSpe

16、ctrum(y)');title('调频信号的频谱图');nsfm=sfm;fori=1:length(t)-1diff_nsfm(i)=(nsfm(i+1)-nsfm(i)./dt;enddiff_nsfmn=abs(hilbert(diff_nsfm);zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn-zero;figure(5) plot(1:length(diff_nsfmn1)./1000,diff_nsfmn1./400,'r');xlabel('时间t

17、9;);title('无噪声条件下解调信号的时域图');x=fft(diff_nsfmn1,1024);f=(0:length(x)-1)*fs/length(x)-fs/2;x=x/1000;figure(6)plot(f,abs(fftshift(x);xlabel('频率f)title('已解调信号的频谱图')%有噪声情况下的调频信号fc=50;fs=1000;t=(0:0.001:0.15);am=5;fm=5;mt=am*cos(2*pi*fm*t)sfm=modulate(mt,fc,fs,'FM');y1=sfm+awgn(

18、sfm,10,0);figure(7)plot(t,y1);xlabel('t(s)');ylabel('sfm');title('加噪声之后的调频信号时域图,性噪比为10')axis(00.15-11);b=fft(y1,1024);f=(0:length(b)-1)*fs/length(b)-fs/2;figure(8)plot(f,abs(fftshift(b);xlabel('Frequence(Hz)');ylabel('PowerSpectrum(y)');title('加噪声后的调频信号频谱图

19、');ny1=y1;fori=1:length(t)-1diff_ny1(i)=(ny1(i+1)-ny1(i)./dt;enddiff_ny1n=abs(hilbert(diff_ny1);zero=(max(diff_ny1n)-min(diff_ny1n)/2;diff_ny1n1=diff_ny1n-zero;figure(9)plot(1:length(diff_ny1n1)./1000,diff_ny1n1./400,'r');xlabel('时间t');title('加性噪声条件下解调信号的时域图');%对已调信号进行傅氏变

20、换x=fft(diff_ny1n1,1024);f=(0:length(x)-1)*fs/length(x)-fs/2;x=x/1000;figure(10)plot(f,abs(fftshift(x);xlabel('频率f);title('加噪声的已解调信号的频谱图');4.4仿真结果(1)产生调制信号和载波信号,画出相应的时域图和频域图图11调制信号时域图图12调制信号的频谱图图13载波信号时域图(2)利用modulate函数产生调制信号FileEditViewjnsertToolsDesktopWindowHlelp口用P昌住ARC：目BQFigure4图14调

21、频信号的时域图FileEditViewInsertTealsDesktopWindowHelp口今Q41,0争黝吃口目口回调频信号的频谱图5004aoFrequence(Hz)'-4OI(AEEoCDdmLc.-5(3)对调制信号进行解调图16无噪声条件下解调信号的时域图(4)对信道中参加加性高斯白噪声,绘制调频信号的时域图和频域图FileEditViewInsertToolsDesktopWindowHelp口齿口聋fe敏Rf?®)要SaE图18加噪声之后的调频信号时域图,信噪比为10加噪声之后的调频信号时域图性噪比为1口tMJFigure8FileEditViewInse

22、rtToolsDesktopWindov*Help膏诒包a琳1)要口目sQ加噪声后的调频信号频谱图Frequence(Hz)(IEnJoOJdsEWJd(5)绘制加噪声的调频信号的解调时域图和频域图图20加性噪声条件下解调信号的时域图5总结5.1 设计小结(1)对已调制未加噪声的调频信号而言,在时域图中出现疏密交错的变化规律,有理论和实际都可以看到,在瞬时最大正频偏处,波形最密；在瞬时最低处,波形最疏(2)由图可知,输入的调制信号通过调制之后,波形发生了明显的变化,原本规那么的正弦信号变成了不规那么的上下起伏波动的图形,而且调制后的图形也没有原本正弦信号般圆滑,出现了十分锋利的突起.说明正弦信

23、号通过FM调制之后波形发生了明显的改变.通过频谱图的对照比拟我们可以看出FM调制并不是使原正弦信号的频谱在原来位置上通过移动得到调制波形,调制后的波形与调制前的完全不同.通过频谱图的对照比拟我们也可以看出FM调制并不是使原正弦信号的频谱在原来位置上通过移动得到调制波形,调制后的波形与调制前的完全不同,这证实FM调制并不是线性的,而是非线性的.(3)改变抽样频率,可以看到当抽样频率低于调制信号的时,无法得出信号波形,由奈奎斯特定义可以知道,当抽样频率低于调制信号时,信号会混叠失真.(4)在参加噪声后时域波形与原来的时域波形比照中,参加噪声后时域波形明显失真,波形不仅不如原本波形般规那么,而且曲线

24、之间还出现了未链接在一起的断裂,但随着信噪比的增大,与原有的波形的相似度也增大了,说明信噪比越大,噪声对信号的影响也变小了.从参加噪声的时域图与频域图的比照中,参加噪声的时域图的变化更为突出,白噪声在整个频谱内每个频点的能量为常数,且根本恒定,所以他对于时域的影响更大.5.1 收获体会本次课程设计让自己对FM调制及MATLAB软件有了较深了解,设计时遇到许多难题,但通过自己努力,以及在老师的帮助下,逐步对自己的课题有了更加深入的了解,刚开始对自己所要做的FM调制原理没有搞清楚,书籍知识也只理解到皮毛而且对MATLAB软件的使用和其语言的使用也不是很会,虽然在数字信号处理过程中学过,但仅仅停留在外表上,只会简单的调试,这也增加了此次设计的对自己的难度,但通过相关书籍的阅读和与老师的交流,在老师,同学的帮助下,完成了课程设计的根本要求.在课程设计的这段时间内主要遇到了以下问题：1 对MATLAB中的函数的理解不够,在信号的调制时,一开始并不知道要用到modulate函数来对未调制信号进行调制,后来经过同学的提醒,简