《通信系统课程设计之基于MATLABFM通信系统.doc》由会员分享，可在线阅读，更多相关《通信系统课程设计之基于MATLABFM通信系统.doc（17页珍藏版）》请在三一办公上搜索。
1、西南科技大学课 程 设 计 报 告课程名称： 通信系统课程设计 设计名称： 基于MATLAB的FM通信系统设计 姓 名： 容晓庆 学 号: 20096025 班 级： 通信0901班 指导教师： 侯宝林 起止日期： 2012.6.17-2012.6.25 西南科技大学信息工程学院课 程 设 计 任 务 书学生班级： 通信0901班 学生姓名： 容晓庆 学号： 20096025 设计名称： 基于MATLAB的FM通信系统仿真 起止日期： 2012.06.10-2012.06.25 指导教师： 侯宝林 设计要求：1)对FM通信系统有整体深入的理解，画出对应的通信系统的原理框图； 2）提出仿真方案，
2、用FM对模拟输入信号m(t)进行调制与解调；3）完成仿真软件的设计并演示；4）提交详细的设计报告。课 程 设 计 学 生 日 志时间设计内容2012.06.11-2012.06.13领取课程设计题目，学习题目要求，预览相关资料2012.06.14-2012.06.16根据设计要求，确定相应方案，准备课程设计要用到的软件2012.06.17-2012.06.19查阅资料，在软件上实现功能设计要求 2012.06.20-2012.06.22调试程序，发现不足之处，完善功能要求2012.06.23-2012.06.24按要求完成设计报告书写2012.06.25答辩课 程 设 计 评 语 表指导教师评
3、语： 成绩： 指导教师： 年 月 日基于MATLAB的FM通信系统仿真一、 设计目的和意义（1）熟悉MATLAB文件中M文件的使用方法，包括函数、原理和方法的应用。（2）加深对FM信号调制原理的理解。（3）增强在通信原理仿真方面的动手能力与自学能力。（4）在做完FM调制仿真之后，在今后遇到类似的问题，学会对所面对的问题进行系统的分析，并能从多个层面进行比较。二、设计原理 图1 模拟通信系统模型调制器: 使信号与信道相匹配, 便于频分复用等。发滤波器: 滤除调制器输出的无用信号。收滤波器: 滤除信号频带以外的噪声，一般设N(t)为高斯白噪声,则Ni(t)为窄带白噪声。在通信系统中一般需要将信号进
4、行相应调制,以利于信号在信道上的传输，调制是将用原始信号去控制高频振荡信号的某一参数，使之随原始信号的变化而成规律变化。调制可分为线性调制和非线性调制。线性调制有AM、DSB等，非线性调制有FM、PM等，这里主要讨论FM调制通信系统1.FM调制原理角调制不是线性调制，角调制中已调信号和调制信号频谱之间不是线性关系而是产生出新的与频谱搬移不同的新的频率分量，呈现非线性特性，故又成为非线性调制。FM调制中瞬时角频率是关于调制信号的线性函数，瞬时角频率偏移量 ,则， 瞬时角频率为 为频偏指数则， 调频信号为 当调制信号是单频余弦时，调制信号为 ，为调频指数，;调制信号的信号的产生 图 2 图 3图2
5、为调频信号的直接产生，图3为调频信号的间接产生，在间接产生方法中，因为实际调相的相位调节器的范围在（-，），故而间接调频只能用于相位偏移和频率偏移不大的窄调制情况，而直接调频则常用于宽带调制情况。窄带角频率调制时，最大瞬时相位偏移远小于30即 调频信号的时域表达式 ，当为窄带调制时，假设调制线号的频谱为F(W),而且假设飞f（t）的平均值为零，有傅氏变换可得 此时调频信号的带宽为调制信号的两倍。其相应图形如下： 调制信号频谱 调频信号频谱2.解调原理1）非相干解调由于调频信号的瞬时频率正比于调制信号的幅度，因而调制信号的调节器必须产生正比于输入频率的输出电压，即输入为调制信号为，则 解调器的输
6、出应为 最简单的解调器是鉴频器 原理图如下：低通滤波器包络检波微分器限幅器及带通Si(t)So(t)包络检波输出为 So(t)= ， 为鉴频器的灵敏度，微分器和包络检波构成鉴频器。2）相干解调由于窄带调频信号可分解为同相分量和正交分量，因而可以用线性调制中的相干解调来进行解调。原理图如下：相乘器低通滤波器微分器带通滤波器Si(t) Sp(t) So(t)带通滤波器的作用是抑制信号带宽以外的噪声，低通滤波器输出 ；上图所示的相干解调法只适用与窄带调频3.调频系统的抗噪声性能主要讨论非相干解调的抗噪声性能，已知输入信号为，输入功率 ，输入噪声,信噪比为,非相干解调的解调器的输入端是调频信号与噪声的
7、叠加，为，在大信噪比情况下，解调器的输出端的信噪比为，考虑单频余弦信号调制，故可得大信噪比情况下的信噪比增益为 ，单频时，带宽,所以增益可化解为 .由此可看出，性噪比增益和调频指数的三次方成正比。加大调频指数，可使调频系统的抗噪声性能迅速改善。二 详细设计步骤1. 产生调制信号和载波信号,画出相应的时域图和频域图2. 利用modulate函数产生调制信号3. 对调制信号进行解调4. 对信道中加入加性高斯白噪声，绘制调频信号的时域图和频域图5.绘制加噪声的调频信号的解调时域图和频域图二、 设计结果及分析1. 对已调制未加噪声的调频信号而言，在时域图中出现疏密交错的变化规律，有理论和实际都可以看到
8、，在瞬时最大正频偏处，波形最密；在瞬时最低处，波形最疏.2. 由图可知，输入的调制信号通过调制之后，波形发生了明显的变化，原本规则的正弦信号变成了不规则的上下起伏波动的图形，而且调制后的图形也没有原本正弦信号般圆滑，出现了十分尖锐的突起。说明正弦信号通过FM调制之后波形发生了明显的改变。通过频谱图的对照比较我们可以看出FM调制并不是使原正弦信号的频谱在原来位置上通过移动得到调制波形，调制后的波形与调制前的完全不同，这证明FM调制并不是线性的，而是非线性的.3. 改变抽样频率，可以看到当抽样频率低于调制信号的时，无法得出信号波形，由奈奎斯特定义可以知道，当抽样频率低于调制信号时，信号会混叠失真.
9、4. 加入噪声后时域波形与原来的时域波形相比，波形明显失真，波形不仅不如原本波形般规则，而且曲线之间还出现了未链接在一起的断裂，但随着信噪比的增大，与原有的波形的相似度也增大了，说信噪比越大，噪声对信号的影响也变小了。从加入噪声的图形与未加入噪声的对比中我们还可以看出噪声对时域图的变化明显比频域图的变化更为突出，白噪声在整个频谱内每个频点的能量为常数，且基本恒定，所以他对于时域的影响更大。三、 体会此次课程设计让自己对FM调制及MATLAB软件有了较深了解，设计时遇到许多难题，但通过自己努力，以及与同学一起探讨，以及在老师的帮助下，逐步对自己的课题有了更加深入的了解，刚开始对自己所要做的FM调
10、制原理没有搞清楚，书籍知识也只理解到皮毛而且对MATLAB软件的使用和其语言的使用也不是很会，虽然在数字信号处理过程中学过，但仅仅停留在表面上，只会简单的调试，对其未进行深入学习，这也增加了此次设计的对自己的难度，但通过相关书籍的阅读和与老师的交流，在老师，同学的帮助下，完成了课程设计的基本要求。在课程设计的这段时间内主要遇到了以下问题：1、 对MATLAB中的函数的理解不够，在信号的调制时，一开始并不知道要用到modulate函数来对未调制信号进行调制，还很辛苦的在在网上找了一大堆程序来完成调制过程，后来经过同学的提醒，简化了程序，简单的完成了对信号的调制。2、 在做时域到频域的变化步骤时，
11、由于要利用傅立叶函数进行变化，刚开始完全不知如何下手，很多书籍上有DFT，IFFT函数，不能明白其意思，后来通过老师的指点，使用FFT可以使时域直接到频域进行变化，从而完成时域到频域的转换。3、 加噪声是对系统抗噪声性能的一种测试，用awgn这个函数时，不知是什么原因，系统总是报错，后来经过对软件的重新启动又可以执行相关功能了，软件的性能也很重要呀。4、 在此次课程设计中遇到的棘手问题只要在MATLAB软件语言上的使用，虽然对FM调制有了一定理解，但却不能将其用相应的语言来表示出FM调制过程，在程序方面主要还是通过网络方面收集，并对其进行整理和分析，综合成相关的程序要求，简单的完成了这次的课程
12、设计。四、 参考文献1曹志刚 钱亚生 现代通信原理. 清华大学出版社2程培青 数字信号处理 清华大学出版社 3聂翔飞、王海宝 MATLAB程序设计及其在信号处理中的应用 西南交通大学出版社 4沈伟慈 通信电路 西安电子科技大学出版社 5 百度文库等网上资源.echo off close allclear allclc.%调制信号的产生.dt=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(
13、调制信号时域图);fc=50; %设定载波频率ct=cos(2*pi*fc*t); %生成载波figure(2)plot(t,ct) ; %设定调频指数xlabel(t)title(载波信号时域图);fs=1000; %调制，产生已调信号sfm=modulate(mt,fc,fs,FM);ts=0.001; %抽样间隔fs=1/ts; %抽样频率df=0.25; %所需的频率分辨率，用在求傅里叶变换%时，它表示FFT的最小频率间隔%*对调制信号m(t)求傅里叶变换*mt=am*cos(2*pi*fm*t); %原调信号fs=1/ts;if nargin=2 n1=0;else n1=fs/df
14、;endn2=length(mt);n=2(max(nextpow2(n1),nextpow2(n2);M=fft(mt,n);mt=mt,zeros(1,n-n2);df1=fs/n; %以上程序是对调制后的信号u求傅里变换M=M/fs; %缩放，便于在频铺图上整体观察f=0:df1:df1*(length(mt)-1)-fs/2; %时间向量对应的频率向量figure(11)plot(f,abs(fftshift(M) %fftshift:将FFT中的DC分量移到频谱中心xlabel(频率f)title(调制信号的频谱图);figure(3)plot(t,sfm);xlabel(t)tit
15、le(调频信号的时域图);axis(0 0.15 -1 1);b=fft(sfm,1024); f=(0:length(b)-1)*fs/length(b) -fs/2;figure(4)plot(f,abs(fftshift(b); %FM信号频谱图xlabel(Frequence(Hz); %频率分量ylabel(Power Spectrum(y); %频率谱能量title(调频信号的频谱图);.%调频信号的解调.nsfm=sfm; for i=1:length(t)-1 %接受信号通过微分器处理 diff_nsfm(i)=(nsfm(i+1)-nsfm(i)./dt;enddiff_ns
16、fmn = abs(hilbert(diff_nsfm); %hilbert变换，求绝对值得到瞬时幅度（包络检波）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); title(无噪声条件下解调信号的时域图);%对已调信号进行傅氏变换x=fft(diff_nsfmn1,1024);f=(0:length(x)-1)*fs/l
17、ength(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(sfm,10,0); % FM信号加入噪声 figure(7)plot(t,y1); %FM信号时域图xlabel(t(s);
18、ylabel(sfm);title(加噪声之后的调频信号时域图,性噪比为10)axis(0 0.15 -1 1);b=fft(y1,1024); f=(0:length(b)-1)*fs/length(b) -fs/2;figure(8)plot(f,abs(fftshift(b); %FM信号频谱图xlabel(Frequence(Hz); %频率分量ylabel(Power Spectrum(y); %频率谱能量title(加噪声后的调频信号频谱图);%对加性噪声调频信号进行解调ny1=y1; for i=1:length(t)-1 %接受信号通过微分器处理 diff_ny1(i)=(ny
19、1(i+1)-ny1(i)./dt;enddiff_ny1n = abs(hilbert(diff_ny1); %hilbert变换，求绝对值得到瞬时幅度（包络检波）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(加性噪声条件下解调信号的时域图);%对已调信号进行傅氏变换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(加噪声的已解调信号的频谱图);
