《基于MATLAB的FM通信系统仿真(共17页).doc》由会员分享，可在线阅读，更多相关《基于MATLAB的FM通信系统仿真(共17页).doc（17页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上西南科技大学课 程 设 计 报 告课程名称： 通信系统课程设计 设计名称： 基于MATLAB的FM通信系统仿真姓 名： 学 号: 班 级： 通信0901 指导教师： 起止日期： 2012.6.4-2011.6.24 西南科技大学信息工程学院制专心-专注-专业课 程 设 计 任 务 书学生班级： 通信0901 学生姓名： 学号： 设计名称： 基于MATLAB的FM通信系统仿真 起止日期： 2012.06.042012.06.24 指导教师： 杨俊平 设计要求：1) 对FM通信系统有整体深入的理解，画出对应的通信系统的原理框图；2) 提出仿真方案，用FM对模拟输入信号m(
2、t)进行调制与解调；3) 完成仿真软件的设计并演示；4) 提交详细的设计报告。课 程 设 计 学 生 日 志时间设计内容2012.06.09查阅资料，确定方案2012.06.12设计总体方案 2012.06.132012.06.20开始整体设计，用FM对模拟输入信号m(t)进行调制与解调，完成仿真软件的设计2012.06.21检查设计内容是否符合要求，并对不足的地方进行修改2012.06.23提交详细设计报告2012.06.25答辩课 程 设 计 评 语 表指导教师评语： 成绩： 指导教师： 年 月 日基于MATLAB的FM通信系统仿真一、设计目的和意义1、熟悉MATLAB的使用方法，了解简单
3、函数、原理和掌握操作方法；2、加深对FM信号调制及解调原理的理解；3、了解FM调制、解调基本结论、规律；4、增加自己在通信原理方面的动手能力；5、通过查资料等提高自学能力，对通信原理在最后有更深入的了解。二、设计原理 1、FM调制原理 如果幅度不变，起始相位为0，而瞬时角频率是调制信号的线性函数，则这种调制方式称为频率调制。 调频信号可表达为 （1）调制信号为单余弦时 （2）当对载波进行频率调制时，可得调频信号表达式为 （3） （4）上式中称为调频指数，也可记作 （5） 图1 图2 图1为直接调频，图2为间接调频。由于实际相位调制器的调节范围不可能超出，因而间接调频的方法仅适合于相位偏移和频率
4、偏移不大的窄调制情形，而直接调频则常用于宽带调制情形。根据调制后载波瞬时相位偏移的大小，可以将角频率分为宽带与窄带两种。宽带与窄带调制的区分并无严格的界限，但通常认为由调频所引起的最大瞬时相位偏移远小于30时，即 （6）称为窄带调频。否则，称为宽带调频。 为方便起见，假设正弦载波的振幅A1，则由式（1）调频信号的一般表达式，有 = （7）通过化解，利用傅立叶变化公式可得NBFM信号的频域表达式： （8） 在NBFM中，由于下边频为负，因而合成矢量不与载波同相，而是存在相位偏移，当最大相位偏移满足式（2-2）时，合成矢量的幅度基本不变，这样就形成了FM信号。 NBFM信号频谱 有两类产生调频信号
5、的方法：直接法和倍频法、在直接法中采用压控振荡器作为产生调频信号的调制器，压控振荡器的输出频率正比于所加的控制电压。直接调频法的优点是可以得到很大的频偏，其主要缺点是载频会发生飘移，因而需要附加的稳频电路。在倍频法中吗，首先用类似于线性调制的的方法产生窄带调频信号，然后用倍频的方法变换为宽带调频信号。2、 FM解调原理（1） 非相干解调由于调频信号的瞬时频率正比于调制信号的幅度，因而调频信号的解调器必须能产生正比于输入频率的输出电压，也就是当输入调频信号为 （9）则解调输出应为：限幅器及带通 （10）微分器包络检波低通滤波器 调频信号的非相干解调 最简单的解调器是鉴频器，理想鉴频器可看成由微分
6、器和包络检波器的级联，限幅器的作用是消除信道中的噪声和其他原因引起的调频波的幅度起伏，带通滤波器（BPF）是让调频信号顺利通过，同时滤除带外噪声及高次谐波分量。微分器的作用是把幅度恒定的调频波变成幅度和频率都随原始信号变化的调幅调频，即 （11）微分器的输出是一个调幅调频信号，如果只取其包络信息，则正比于调制信号，因而滤去直流后，包络检波器的输出为 (12) 上述解调方法称为包络检波，其缺点之一是包络检波器对于由信道中噪声和其他原因引起的幅度起伏也有反应，因而在使用中常在微分器之前加一个限幅器和带通滤波器。（2） 相干解调带通 相乘器低通微分 窄带调频信号的相干解调由于窄带调频信号可分解成同相
7、分量和正交分量之和，因而可以采用线性调制中的相干解调法来进行解调。图中带通滤波器用来限制信道所引入的噪声，但调频信号应能正常通过。图中所示相干解调器的输出正比于调制信号f(t)，相干解调法只适用于窄带调频。3、 模拟通信系统简介 传递信息所需的一切技术设备的总和称为通信系统。通信系统由信息源和收信者、发送设备、传输媒介、接收设备几部分组成。通信系统的一般模型如下图所示：信息源收信者接收设备传输媒介传输设备干扰通信系统是为了有效可靠的传输信息，信息由信源发出，以语言、图像、数据为媒体,通过电(光)信号将信息传输，由信宿接收。通信系统又可分为数字通信与模拟通信。角调制属于模拟通信。信源是模拟信号，
8、信道中传输的也是模拟信号的系统为模拟通信。模拟通信系统的模型如下图所示：调制器: 使信号与信道相匹配, 便于频分复用等。发滤波器: 滤除调制器输出的无用信号。收滤波器: 滤除信号频带以外的噪声，一般设N(t)为高斯白噪声,则Ni(t)为窄带白噪声。4、高斯白噪声信道特性设正弦波通过加性高斯白噪声信道后的信号为 （13）其中，白噪声的取值的概率分布服从高斯分布。调用MATLAB本身自带的标准高斯白噪声函数awgn，再加上调频信号就可得到加有噪声的信号。在实验中选择10db和30db两种不同信噪比的高斯白噪声进行比较。三、 详细设计步骤 在程序设计时，采用直接调频的调制方法，利用matlab中的函
9、数modulate，直接得到调频信号表达式。要求出调制信号的时域和频域曲线，可以根据表达式利用plot函数画出时域波形，再将进行傅里叶变换得到，根据表达式利用plot函数画出频域波形。要求出调频信号sfm的时域和频域曲线。根据调频信号的一般表达式 （14）可知，要想求出调频信号sfm的表达式，须得先求出的积分，然后再求出调频信号的表达式，根据求出的表达式sfm利用plot函数画出时域波形，再将sfm进行傅里叶变换，再利用plot函数画出频域波形。要求出解调信号diff_sfm的时域和频域曲线。由于非相干解调对NBFM信号和WBFM信号均适用，而相干解调只适用于NBFM信号，根据题目要求，采用非
10、相干解调方法。此时解调器的输出应为：。根据表达式利用plot函数画出时域波形，再将进行傅里叶变换，再利用plot函数画出频域波形。1、求已调信号（无噪声条件下）的时域图如下：2、求无噪声条件下解调信号的时域曲线 由于非相干解调对NBFM信号及WBFM信号均适用，所以采用非相干的解调方法。3、 加入小信噪比高斯白噪声之后的调频、解调时域图4、 加入大信噪比高斯白噪声之后的调频、解调时域图5、 求原调制信号、调频信号、解调信号的频谱图对FM调制前后的时域波形用FFT函数进行傅利叶变换，进行傅立叶变化便分别得到调制信号与调制之后的FM信号以及解调信号的频谱图。利用plot函数画出频域波形。如下图所示
11、：6、 求原调制信号、加入小信噪比高斯白噪声后的调频信号、解调信号的频谱图7、求原调制信号、加入大信噪比高斯白噪声后的调频信号、解调信号的频谱图四、 设计结果及分析 从调频信号的时域波形图可以看出，调频信号的频率随调制信号而变化，产生频偏，成为疏密波形。调制信号电平最高处对应的瞬时正频偏最大，波形最密，而在调制信号电平最低处对应的瞬时频偏最小，波形最疏。在无噪声干扰的情况下，解调出的波形和原调制信号波形一样，只是幅度有点不同。当加入小信噪比高斯白噪声之后，调频信号不是很光滑，有很多毛刺，解调出的信号波形也出现失真，有毛刺；当加入大信噪比高斯白噪声后，噪声干扰相对较小，所以调频信号波形还是比较光
12、滑，解调出的波形虽然存在一定的失真，但是波形还是和原调制信号波形差不多。大信噪比时宽带调频系统的信噪比增益是很高的，它与调频指数的立方成正比。大信噪比时宽带调频系统的信噪比增益是很高的，它与调频指数的立方成正比。因此，加大调频指数，可使调频系统的抗噪声性能迅速改善。从它们的频谱图可以看出，调频信号的频谱是由载频和无穷多组上、下边频组成，调频信号是一种非线性频率变换的过程，其解调也是一种非线性频率变换过程。理论上调频信号的频带应为无限宽，但是从频谱图可以看出，其具有较大振幅的频率分量还是集中在载频附近，且上下边频在振幅上是对称的。 FM调制与幅度调制相比，最突出的优势在于其较高的抗噪声性能，但获
13、得这种优势的代价是角度调制占用比幅度调制信号更宽的带宽。调频是幅度恒定的已调信号，发射时，可采用高效率的丙类谐振功率放大器将它放大到所需的发射功率，而在接收这些已调信号时将呈现出很强的抗干扰能力。FM波的幅度恒定不变，这使得它对非线性器件不甚敏感，给FM带来了抗快衰落能力。利用自动增益控制和带通限幅还可以消除快衰落造成的幅度变化效应。这些特点使得NBFM对微波中继系统颇具吸引力。WBFM的抗干扰能力强，可以实现带宽与信噪比的互换，因而WBFM广泛应用于长距离高质量的通信系统中，如空间和卫星通信、调频立体声广播、短波电台等。WBFM的缺点是频带利用率低，存在门限效应，因此在接收信号弱、干扰大的情
14、况下宜采用NBFM，这就是小型通信机常采用NBFM的原因。五、体会在做课程设计时，由于对FM和MATLAB学得不是很好，也没有更深入的了解过FM，所以刚开始时遇到了很多问题。由于我对MATLAB也不是很了解，使得我的编程过程更加艰难，之后通过借来的书籍及老师的辅导后，自己逐步对该课题有了一定的了解，逐步完成了课程设计。还有就是在编程时，由于对MATLAB中函数的理解不够，而造成我在利用modulate这个函数时不知如何运用，之后通过MATLAB中的帮助选项后，再通过自己不断调试逐步学会了这个函数的运用。解调时，本来想调用matlab里的解调函数的，但是用了之后解调出来的图形不对，修改了好几次都
15、不对，所以后来就在网上搜了一个用于解调的函数，完成了解调这一过程。在做时域到频域的变化步骤时，由于要利用傅立叶函数进行变化，刚开始完全不知该如何用代码来表示傅立叶变换。最后查资料发现可以直接用matlab函数里的fft实现。通过这次课程设计，我了解了很多相关专业的知识，之前在学数字信号处理的时候也有用到matlab，但那时没有认真学过这个软件，在做课程设计时，查了很多关于matlab的资料，发现这个软件真的很强大，应该好好研究下这个软件。虽然在课程设计开始时遇到许多难题，但通过自己努力以及与同学一起探讨，逐步加深了理解，在老师的辅导下逐步完成了此次的课程设计。六、参考文献【1】现代通信原理，曹
16、志刚、钱亚生编著，清华大学出版社【2】MATLAB程序设计，聂祥飞、王海宝、谭泽富主编，西南交通大学出版社【3】数字信号处理教程，程佩青编著，清华大学出版社【4】信号与系统（第二版），刘树棠，西安交通大学出版社【5】通信电路（第二版），沈伟慈编著，西安电子科技大学出版社实验程序：dt=0.001; %设定时间步长t=0:dt:1.5; %产生时间向量am=5; %设定调制信号幅度fm=5; %设定调制信号频率mt=am*cos(2*pi*fm*t); %生成调制信号fc=50; %设定载波频率ct=cos(2*pi*fc*t); %生成载波fs=1000; %抽样频率sfm=modulate(
17、mt,fc,fs,FM); %调频信号sfm1=sfm+awgn(sfm,10,0); %加入小信噪比后的调频信号sfm2=sfm+awgn(sfm,30,0); %加入大信噪比后的调频信号%*FM解调*for i=1:length(t)-1 %接受信号通过微分器处理 diff_sfm(i)=(sfm(i+1)-sfm(i)./dt;enddiff_sfm1 = abs(hilbert(diff_sfm); %hilbert变换，求绝对值得到瞬时幅度（包%络检波）zero=(max(diff_sfm1)-min(diff_sfm1)/2;diff_nsfm1=diff_sfm1-zero; %
18、无噪声干扰下的解调信号for i=1:length(t)-1 %接受信号通过微分器处理 diff_sfm1(i)=(sfm1(i+1)-sfm1(i)./dt;enddiff_sfm2 = abs(hilbert(diff_sfm1); %hilbert变换，求绝对值得到瞬时幅度（包%络检波）zero=(max(diff_sfm2)-min(diff_sfm2)/2;diff_nsfm2=diff_sfm2-zero;%小信噪比高斯白噪声下的解调信号for i=1:length(t)-1 %接受信号通过微分器处理 diff_sfm2(i)=(sfm2(i+1)-sfm2(i)./dt;endd
19、iff_sfm3 = abs(hilbert(diff_sfm2); %hilbert变换，求绝对值得到瞬时幅度（包%络检波）zero=(max(diff_sfm3)-min(diff_sfm3)/2;diff_nsfm3=diff_sfm3-zero;%大信噪比高斯白噪声下的解调信号%*%*对调制信号m(t)求傅里叶变换*x=fft(mt,1024);f=(0:length(x)-1)*fs/length(x)-fs/2;x=x/1000; %缩放，便于在频谱图上整体观察%*对已调信号求傅里叶变换*y=fft(sfm,1024);f=(0:length(y)-1)*fs/length(y)-
20、fs/2;y=y/1000;y1=fft(sfm1,1024);f=(0:length(y1)-1)*fs/length(y1)-fs/2;y1=y1/1000;y2=fft(sfm2,1024);f=(0:length(y2)-1)*fs/length(y2)-fs/2;y2=y2/1000;%*对解调信号求傅里叶变换*z=fft(diff_nsfm1,1024);f=(0:length(z)-1)*fs/length(z)-fs/2;z=z/1000;z2=fft(diff_nsfm2,1024);f=(0:length(z2)-1)*fs/length(z2)-fs/2;z2=z2/10
21、00;z3=fft(diff_nsfm3,1024);f=(0:length(z3)-1)*fs/length(z3)-fs/2;z3=z3/1000;%原调制信号、载波信号和已调信号波形%*figure(1)*figure(1)subplot(3,1,1);plot(t,mt); %绘制调制信号的时域图xlabel(时间t);title(调制信号的时域图);subplot(3,1,2);plot(t,ct); %绘制载波的时域图xlabel(时间t);title(载波的时域图);subplot(3,1,3);plot(t,sfm); %绘制已调信号的时域图xlabel(时间t);title(
22、无噪声干扰下的已调信号的时域图);%原调制信号、无噪声条件下已调信号和解调信号波形%*figure(2)*figure(2)subplot(3,1,1);plot(t,mt); %绘制调制信号的时域图xlabel(时间t);title(调制信号的时域图);subplot(3,1,2);plot(t,sfm); %绘制已调信号的时域图xlabel(时间t);title(无噪声条件下已调信号的时域图); subplot(3,1,3); %绘制无噪声条件下解调信号的时域图plot(1:length(diff_nsfm1)./1000,diff_nsfm1./400,r);xlabel(时间t); t
23、itle(无噪声条件下解调信号的时域图);%原调制信号、加入小信噪比高斯白噪声后的载波信号和已调信号的曲线%*figure(3)*figure(3)subplot(3,1,1);plot(t,mt); %绘制调制信号的时域图xlabel(时间t);title(调制信号的时域图);subplot(3,1,2);plot( t,sfm1); %绘制含小信噪比高斯白噪声已调信号的时域图xlabel(时间t);title(含小信噪比高斯白噪声已调信号的时域图);subplot(3,1,3); %绘制含小信噪比高斯白噪声解调信号的时域图plot(1:length(diff_nsfm2)./1000,di
24、ff_nsfm2./400,r);xlabel(时间t); title(含小信噪比高斯白噪声解调信号的时域图);%原调制信号、加入大信噪比高斯白噪声后的载波信号和已调信号的曲线%*figure(4)*figure(4)subplot(3,1,1);plot(t,mt); %绘制调制信号的时域图xlabel(时间t);title(调制信号的时域图);subplot(3,1,2);plot( t,sfm2); %绘制含大信噪比高斯白噪声已调信号的时域图xlabel(时间t);title(含大信噪比高斯白噪声已调信号的时域图);subplot(3,1,3); %绘制含大信噪比高斯白噪声解调信号的时域
25、图plot(1:length(diff_nsfm3)./1000,diff_nsfm3./400,r);xlabel(时间t); title(含大信噪比高斯白噪声解调信号的时域图);%原调制信号和已调信号、解调信号在频域内的图形%*figure(5)*figure(5)subplot(3,1,1);plot(f,abs(fftshift(x) ; %fftshift:将FFT中的DC分量移到频谱中心xlabel(频率f);title(原调制信号的频谱图);subplot(3,1,2);plot(f,abs(fftshift(y);xlabel(频率f);title(已调信号的频谱图);subp
26、lot(3,1,3);plot(f,abs(fftshift(z);xlabel(频率f);title(解调信号的频谱图);%原调制信号和加入小信噪比噪声后的已调信号、解调信号在频域内的图形%*figure(6)*figure(6)subplot(3,1,1);plot(f,abs(fftshift(x) ; %fftshift:将FFT中的DC分量移到频谱中心xlabel(频率f);title(原调制信号的频谱图);subplot(3,1,2);plot(f,abs(fftshift(y1);xlabel(频率f);title(小信噪比下已调信号的频谱图);subplot(3,1,3);pl
27、ot(f,abs(fftshift(z2);xlabel(频率f);title(小信噪比下解调信号的频谱图);%原调制信号和加入大信噪比噪声号后的已调信号、解调信号在频域内的图形%*figure(7)*figure(7)subplot(3,1,1);plot(f,abs(fftshift(x) ; %fftshift:将FFT中的DC分量移到频谱中心xlabel(频率f);title(原调制信号的频谱图);subplot(3,1,2);plot(f,abs(fftshift(y2);xlabel(频率f);title(大信噪比下已调信号的频谱图);subplot(3,1,3);plot(f,abs(fftshift(z3);xlabel(频率f);title(大信噪比下解调信号的频谱图);