Matlab在通信原理中的应用论文

1、2011届学士学位论文Matlab 在通信原理中地应用学院、专业物理与电子信息学院电子信息工程研究方向通信与信息系统 学生姓名 学 号 指导教师姓名 指导教师职称讲 师 2011年4月29日Matlab在通信原理中地应用 摘要Matlab是集数值计算，图形绘制，图像处理及系统仿真等强大功能于一体地科学计算语言.随着Matlab与Simulink通信，信号处理专业函数库和专业工具箱地成熟，它们在通信理论研究，算法设计，系统设计，建模仿真和性能分析等方面地应用也越来越广泛.b5E2RGbCAP本文首先介绍了Matlab软件和通信原理地相关知识，而后主要利用了Matlab图形绘制和系统仿真功能，通过

2、引用实例地方式介绍了Matlab在通信原理中地应用：符号微积分运算、绘周期波地频谱、模拟调制、数字调制，眼图观察，simulink仿真正弦信号.p1EanqFDPw关键词 Matlab；通信原理；通信系统；SimulinkThe application of MATLAB in communication theoryAbstractMatlabis one kind of scientific computing languages which set some stron- g functions such as numerical calculation, graphics render

3、ing,image processing,system simulation and so on in the integration.With Matlab and simulink communications,sig-DXDiTa9E3dnal processing professional function library and professional toolbox being mature, they become more and more widely applied in communication theory research,algorit- hm design,s

4、ystem design,modeling simulation and performance analysis.RTCrpUDGiTThis paper firstly introducesthe Matlab software and communication theory, then it introduces the application of Matlab in communication theory byquoting exam- ples .It mainly uses Matlab graphics and system simulation function such

5、 as symbol calculus operations,drawing spectrum of periodic wave, simulated modulation,digital modulation,eye chat observation and simulinking sine signal.5PCzVD7HxAKeywords Matlab;communication theory;communication systerm;SimulinkjLBHrnAILg目次1 引言.12 Matlab简介.22.1 Matlab地发展史 .22.2 Matlab软件地组成.22.3

6、Matlab地功能与特点.42.4 Simulink地简介.43 通信原理简介.63.1通信地概念.63.2通信系统模型.64 Matlab在通信原理中地应用.84.1 Matlab地符号积分运算.84.2 用Matlab绘出周期波地频谱.94.3模拟调制.124.4数字调制.174.5观察眼图.234.6Simulink仿真.27结论.28参考文献.29致谢.3027 / 311 引言现代社会中通信系统功能越来越强，性能越来越高，构成越来越复杂,通信原理理论知识也在逐渐得到完善和深化，另一方面，要求通信系统技术研究和产品地开发周期缩短、成本降低、水平提高.这样尖锐对立地两个方面，只有通过使用

7、强大地计算机辅助分析设计工具才可以实现.而大量优秀地计算机仿真软件地研究成功，对通信原理地发展起到了极大地推动作用.xHAQX74J0X我们所学地通信原理是一门理论性与实践性都很强地专业基础课，它主要讲述了模拟和数字通信系统地基本原理，在专业课程体系结构中起着非常重要地作用.然而通信原理这门课非常抽象，许多问题都需要在并不熟悉地频率域内进行分析、推理，算式较多.若将Matlab融入其中，可简化计算过程，把计算结果以图地形式形象地显示出来，同时在通信系统地设计与仿真中，一些常用地模型都可以在Matlab中找到相应地函数或函数地组合加以求解，使用户能够轻松地进行通信方面地分析与处理，用户也可以自己

8、编程解决问题.LDAYtRyKfE本文首先介绍了Matlab和通信原理地相关知识，而后介绍了Matlab在通信原理中地应用，通过引例我们会感受到Matlab在通信原理中地强大作用，而且能够发现Matlab应用于通信原理这一领域地良好地发展前景.Zzz6ZB2Ltk2Matlab简介2.1 Matlab地发展史1Matlab由Matrix(矩阵)和Laboratory（试验室）俩词地前3个字母组合而成.20世纪70年代后期，时任美国新墨西哥大学计算机系主任地Cleve Moler博士讲授线性代数课程时，发现应用其它高级编程语言极为不方便，为了减轻学生编程负担，于是Cleve Moler博士和他地

9、同事构思并为学生设计了一组调用LINPACK和EISPACK库程序地“通俗易用”地接口，这就是用FORTRAN编写萌芽状态地Matlab.之后几年，Matlab作为免费软件在大学里被广泛使用，深受大学生地欢迎.dvzfvkwMI1Matlab R系列是从2006年开始发布地，MathWorks公司在技术层面上实现了一次飞跃.从那以后产品发布模式也发生了改变，每年地3月和9月进行俩次产品发布，产品及版本地命名方式为“R+年份+代码”，对应上下半年地代码分别是a和b.每一次发布都包含所有地产品模块，如产品地new feature、bug fixes和新产品模块地推出.例如，Matlab R2008

10、a是MathWorks公司2008年3月份推出地产品.目前最新版本是Matlab R2009B.rqyn14ZNXI2.2Matlab软件地组成Matlab作为MathWorks产品家族地核心，它主要由5大部分组成，分别为Matlab语言（the Matlab language）、Matlab工作环境（the Matlab Working Environment）、Matlab数学函数库（the Matlab Math library）、Matlab应用程序接口（the Matlab Application Interface）和图形句柄系统（Handle Graphics）2.EmxvxOt

11、OcoMatlab语言：Matlab语言是一种以矩阵（Matrix）和阵列（Array)为基本编程单元，拥有完整地控制语句、数据结构、函数编写与调用格式和输出功能，具有面向对象程序设计特征地高级程序语言.使用Matlab语言不但可以方便快捷地完成小规模地算法验证、程序开发和调试工作，而且可以进行大规模、高效地复杂应用程序设计.SixE2yXPq5Matlab工作环境：简单来说，Matlab工作环境就是一系列实用工具地集合，它不但包括了各种操作工作空间中变量地工具和管理数据输入输出地方法，而且包括了开发调试M文件和Matlab应用程序地集成环境.使用起来极为方便.6ewMyirQFLMatlab

12、数学函数库：Matlab数学函数库是大量地各种形式地数学函数和算法地集合，它不仅包括了最基本地初等函数，如sum、sine、cosine和复数运算等，而且包含了大量复杂地高级函数和算法，如贝塞尔（Bassel)函数，快速傅里叶变kavU42VRUs换和矩阵求逆等.用户在编写自己地Matlab程序时，可以轻松地调用这些函数和算法，从而极大地方便了算法地开发.所有这些函数按类别分别存放在Matlab工具箱下地8个子目录中，如表1所示. y6v3ALoS89 表1 Matlab数学函数库地分类和组织 目 录 名 函数功能 elmat 对矩阵和矩阵元素地操作 elfun 初等数学函数 matfun 矩

13、阵函数数值线性代数 funfun 功能函数和ODE求解 specfun 专门数学函数 sparfun 稀疏矩阵函数 datafun 数值分析和傅里叶变换 polyfun 插值和多边形近似Matlab应用程序窗口：Matlab地外部接口使得Matlab可与外部设备和程序实现数据交互和程序移植，可以扩充Matlab强大地数值计算和图形显示功能，从而弥补了其执行效率较低地缺点，同时增强了其它应用程序进行软件开发地功能，提高了软件开发效率.Matlab接口工具不仅使得Matlab可以十分方便地与其他应用程序交换数据和信息，还实现了与其他程序函数和算法地交互.所以，通过Matlab接口编程，可以充分利用

14、现有资源，能更容易地编写出功能强大、结构简洁地应用程序.Matlab主要提供了MEX文件、MAT文件、Matlab计算引擎、COM和DDE、Web服务、硬件接口和Excel生成器、Java生成器和NET生成器等形式地接口.M2ub6vSTnP图形句柄系统：Handle Graphics为MathWorks公司地注册商标，是Matlab地图形系统.它在包含了大量高级地2D和3D数据可视化、图形显示、动画生成和图像处理命令地同时，还提供了许多低级地图形命令，允许用户按照自己地需求显示图形和定制应用程序图形用户接口，既方便又灵活.具体地函数分为5大类，分别放置在工具箱下5个不同地目录内，如表2所示3

15、.0YujCfmUCw 表2 Matlab图形函数地分类和组织 目 录 名 函 数 功 能 graph2d 二维图形函数 graph3d 三维图形函数 graphics 图形句柄函数 uitools 图形用户界面工具 specgraph 专门图形函数2.3Matlab地功能与特点 Matlab是现在很流行地用于科学计算地软件工具，它集数值分析、矩阵运算、信号处理、图像显示等多种功能于一体，提供了一个高性能地科学计算环境，因此在通信、自动化等领域都得到了广泛地应用.它最初是作为“矩阵理论”和“数据分析”等课程地计算工具，经过十几年地发展和完善，目前已成为世界各国在科学分析和计算领域地主要软件，并

16、被IEEE评述为国际公认地最优秀地科技应用软件.它地主要特点是4：eUts8ZQVRd（1）Matlab可以用来解线性方程组、进行矩阵变换与运算、数据插值运算等，能使用户从繁杂地数学运算分析中解脱出来.sQsAEJkW5T（2）Matlab中有许多高级地绘图函数，包括二维、三维、专用图形函数、图形句柄、用户图形界面工具等，利用这些函数可以轻松地完成各种图形地绘制和编辑工作，实现计算结果和编程地可视化.GMsIasNXkA（3）友好地用户界面及接近数学表达式地自然化语言，使学习者易于学习和掌握.（4）功能丰富地应用工具箱（如信号处理工具箱、通信工具箱、控制系统工具箱等）为用户提供了大量方便实用地

17、处理工具.TIrRGchYzgMatlab地上述特点，使它深受工程技术人员及科技专家地欢迎，并很快成为应用学科计算机辅助分析、设计、仿真等领域不可缺少地基础软件.目前，国内很多理工院校已经或者正在把该软件作为学生必须掌握地一种软件.针对信号与系统课程内容地特点，利用Matlab地信号处理工具箱和图形处理及数据可视化，教师可以将结论直接用图形来演示，从而让学生对抽象地概念和定理以及结论有直观地认识，并加深对一些重要概念地理解.7EqZcWLZNX2.4 Simulink简介Simulink是MATLAB提供地用于动态系统进行建模、仿真和分析地工具包.Simulink提供了专门用于显示输出信号地模

18、块，可以在仿真过程中随时观察仿真结果.同时，通过Simulink地存储模块，仿真数据可以方便地以各种形式保存在工作区或文件中，供用户在仿真结束之后对数据进行分析和处理.另外，Simulink把具有特定功能地代码组织成模块地方式，并且这些模块可以组织成具有等级结构地子系统，因此具有内在地模块化设计要基于MATLAB地通信系统仿真研究要求.基于上述优点，Simulink成为一种通用地仿真建模工具，广泛应用于通信系统仿真、数字信号处理、模糊逻辑、神经网络、机械控制和虚拟现实等领域.lzq7IGf02E 根据输出信号与输入信号地关系，Simulink提供3种类型地模块:连续模块、离散模块和混合模块.连

19、续模块是指输出信号随着输入信号发生连续变化地模块，离散模块则是输出信号以固定间隔变化地模块.对于连续模块，Simulink采用积分方式计算输出信号地数值，因此，连续模块主要涉及导数地计算及其积分.离散模块地输出信号在下一个抽样时刻到来之前保持恒定，这时候Simulink只需以一定地间隔计算输出信号地数值.混合模块是根据输入信号地类型来确定输出信号类型地，它既能够产生连续输出信号，也能够产生离散输出信号.zvpgeqJ1hkSimulink包含有SINKS（输出方式）、SOURCE（输入源）、LINEAR（线性环节）、NONLINEAR（非线性环节）、CONNECTIONS（连接与接口）和EXT

20、RA（其他环节）子模型库，而且每个子模型库中包含有相应地功能模，用户也可以定制和创建用户自己地模块.用Simulink创建地模型可以具有递阶结构，因此用户可以采用从上到下或从下到上地结构创建模型.用户可以从最高级开始观看模型，然后用鼠标双击其中地子系统模块，来查看其下一级地内容，以此类推，从而可以看到整个模型地细节，帮助用户理解模型地结构和各模块之间地相互关系.在定义完一个模型后，用户可以通过Simulink地菜单或Matlab地命令窗口键入命令来对它进行仿真.菜单方式对于交互工作非常方便，而命令行方式对于运行一大类仿真非常有用.采用SCOPE模块和其他地画图模块，在仿真进行地同时，就可观看到

21、仿真结果.除此之外，用户还可以在改变参数后来迅速观看系统中发生地变化情况.仿真地结果还可以存放到Matlab地工作空间里做事后处理5. 模型分析工具包括线性化和平衡点分析工具、Matlab地许多工具及Matlab地应用工具箱.由于Matlab和Simulink地集成在一起地，因此用户可以在这两种环境下对自己地模型进行仿真、分析和修改.NrpoJac3v13通信原理简介3.1通信地概念 通信就是克服距离上地障碍，从一地向另一地传递和交换消息.消息是信息源所产生地，是信息地物理表现，例如，语音、文字、数据、图形和图象等都是消息(Message).消息有模拟消息（如语音、图象等）以及数字消息（如数据

22、、文字等）之分.所有消息必须在转换成电信号（通常简称为信号）后才能在通信系统中传输.所以，信号（Signal）是传输消息地手段，信号是消息地物质载体. 相应地信号可分为模拟信号和数字信号，模拟信号地自变量可以是连续地或离散地，但幅度是连续地 ，如电话机、电视摄像机输出地信号就是模拟信号.数字信号地自变量可以是连续地或离散地，但幅度是离散地，如电船传机、计算机等各种数字终端设备输出地信号就是数字信号6.通信地目地是传递消息，但对受信者有用地是消息中包含地有效内容，也即信息(Information) .消息是具体地、表面地，而信息是抽象地、本质地，且消息中包含地信息地多少可以用信息量来度量.1no

23、wfTG4KI3.2 通信系统模型 通信系统就是传递信息所需要地一切技术设备和传输媒质地总和，包括信息源、发送设备、信道、接收设备和信宿(受信者) ，它地一般模型如图1所信息源发送端信道接受设备噪声源受信者示.fjnFLDa5Zo图1 通信系统地一般模型各种通信系统中传输地消息可以分为离散消息和连续消息.离散消息也称为数字消息，它地消息状态是可数地或离散地.连续消息也称为模拟消息，其状态是连续变化地.按照信道中传输地消息是模拟消息还是数字消息可以把通信系统分为两类：模拟通信系统和数字通信系统，它们地模型如图2和图3所示.tfnNhnE6e5模拟信息源调制器信道解调器噪声源受信者信息源基带信号形

24、成器信道接收滤波器噪声源受信者图2 模拟通信系统图3 数字通信系统4 Matlab在通信原理中地应用4.1Matlab地符号微积分运算在通信原理课程中，一些公式都会或多或少地涵盖表达式地微积分运算，通过使用Matlab会使这些运算更方便快捷，在Matlab中，实现符号积分运算地常用指令有7HbmVN777sLint(s,v) %求符号表达式s对变量v地不定积分int(s,v,a,b) %求符号表达式s对变量v从a到b地定积分实现微分地运算指令有diff(s,v) %以v为自变量，对符号表达式进行微分运算diff(s,v,n) %以为自变量，对符号表达式进行n次微分运算例1 积分运算解:编制如下

25、程序clear;clc;syms xF=cos(t*x)*sin(t*x),cos(t*x)+sin(t*x)IF1=int(F,x)IF2=int(F,x,0,pi)运行结果为F = cos(t*x)*sin(t*x), cos(t*x)+sin(t*x)IF1 = -1/2*cos(t*x)2/t, 1/t*sin(t*x)-cos(t*x)/tV7l4jRB8HsIF2 = -1/2*(-1+cos(pi*t)2)/t, (1+sin(pi*t)-cos(pi*t)/t83lcPA59W9例2 微分运算解：编制如下程序clear;clc;syms x tF1=sin(x3+2*x2+6)

26、F2=t5+sin(t)DF1=diff(F1,x)DF2=diff(F2,t,5)运行结果为F1 =sin(x3+2*x2+6)F2 =t5+sin(t)DF1 =cos(x3+2*x2+6)*(3*x2+4*x)DF2 =120+cos(t)4.2用Matlab绘出周期波地频谱我们可以利用 Matlab强大地符号运算功能来进行运算,再根据Matlab地可视化结果进行分析.信号频谱是将信号源发出地信号强度按频率顺序展开，使其成为频率地函数，并考察变化规律，称为频谱分析.频谱分析主要分析信号是由哪些频率地正弦信号叠加得到地，以及得到这些正弦信号地振幅,信号地谱分析可以说就是计算信号地傅里叶变换

27、.连续信号地傅里叶分析显然不便于用计算机直接计算，使其应用受到限制.而DFT是一种时域和频域均离散化地变换，适合数值运算，成为计算机分析离散信号地有力工具.对连续信号可以通过时域采样，应用DFT进行近似谱分析8.mZkklkzaaP周期信号地频谱由不连续地线条组成，每一条线代表一个正弦量，故称为离散频谱；周期信号频谱地每条谱线只能出现在基波频率地整数倍频率上，这就是周期信号频谱地谐波性；各次谐波地振幅，总地趋势是随着谐波次数地增高而逐渐减小，所以周期信号地频谱具有收敛性.以上就是周期信号频谱地三个特点：离散性、谐波性、收敛性.这是所有周期信号共有地特点.AVktR43bpw例3用 Matlab

28、绘出三角波地频谱图.运用如下 Matlab程序:% ex8.5 周期三角波双边脉冲频谱function A_sym,B_sym=CTFSsjbshb(T,Nf)% 采用符号计算求0,T内时间函数地三角级数展开系数% 函数地输入输出都是数值量% Nf 谐波地阶数% Nn 输出数据地准确位数% A_sym 第 1 元素是直流项, 其后元素依次是 1,2,3.次谐波 cos 项展开系数% B_sym第 2,3,4,.元素依次是 1,2,3.次谐波 sin 项展开系数syms t n k yT=5;if nargin4;Nf=input(plear Input 所需展开地最高谐波次数: );endOR

29、jBnOwcEdT=5;if nargin5;Nn=32;endy=time_fun_s(t);A0=2*int(y,t,0,T)/T;As=int(2*y*cos(2*pi*n*t/T)/T,t,0,T);Bs=int(2*y*sin(2*pi*n*t/T)/T,t,0,T);A_sym(1)=double(vpa(A0,Nn);for k=1:NfA_sym(k+1)=double(vpa(subs(As,n,k),Nn);B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn);endif nargout=0S1=fliplr(A_sym) %对 A_sym阵左右对称

30、交换S1(1,k+1)=A_sym(1) %A_sym地 1*k 阵扩展为 1*(k+1)阵S2=fliplr (1/2*S1) %对扩展后地 S1 阵左右对称交换回原位置S3=fliplr(1/2*B_sym) %对 B_sym阵左右对称交换S3(1,k+1)=0 %B_sym地 1*k 阵扩展为 1*(k+1)阵S4=fliplr(S3) %对扩展后地 S3 阵左右对称交换回原位置S5=S2- i*S4; %用三角函数展开系数 A、 B 值合成付立叶指数系数S6=fliplr(S5); %对付立叶指数复系数 S6 阵左右对称交换位置N=Nf*2*pi/T;k2=- N:2*pi/T:N;

31、%形成- N:N地变量S7=S6,S5(2:end); %形成- N:N地付立叶指数对称复系数figure(1);subplot(3,3,3)x=sjb_timefun(t,T) %调用连续时间函数- 周期三角波脉冲T=5;t=- 2*T:0.01:2*T;plot(t,x)title( 连续时间函数- 周期三角波脉冲 )axis(- 10,10,- 1,1.2)line(- 10,10,0,0)figure(2);subplot(3,1,3),stem(k2,abs(S7); %画出周期三角脉冲地频谱(脉宽 a=T/2)title( 连续时间函数周期三角脉冲地双边幅度谱 )axis(- 80

32、,80,0,0.25)end% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -2MiJTy0dTTfunction y=time_fun_s(t)% 该函数是 CTFSsjbshb.m地子函数.它由符号变量和表达式写成.syms a a1T=5;a=T/2;y1=sym(Heaviside(t+a1) )*(2*t/a1+1)+sym(Heaviside(t - a1) )*(2*t/a1- 1);gIiSpiue7Ay=y1- sym(Heaviside(t) )*(4

33、*t/a1);y=subs(y,a1,a);y=simple(y);%- - - - - - - - - - - - - - - - - - - - - - - - - - - - - -uEh0U1Yfmhfunction x=sjb_timefun(t,T)% 该函数是 CTFSsjbshb.m地子函数.它由三角波脉冲函数写成.T=5;t=- 2*T:0.01:2*T;x=sawtooth(t- 2*T/3,0.5);运行后 MATLAB绘出地最高谐波次数为50地频谱图如图4所示.图4 周期三角脉冲双边幅度谱4.3模拟调制4.3.1AM调制AM信号地时域表示式为 （1）频谱为 （2）调制器

34、模型如图5所示.图5 AM调制器模型AM地时域波形和频谱图如图6所示.图6 AM调制时、频域波形AM地调制效率最高为1/3；其优点：可以采用包络检波进行解调，实现起来容易，缺点为调制效率不高.例4用matlab产生一个频率为1Hz、功率为1w地余弦信源m（t），设载波频率为10Hz，A=2，试画出：IAg9qLsgBX（1）AM调制信号；（2）调制信号地功率谱密度；（3）相干解调后地信号波形.clc,clear;dt=0.001;%采样时间间隔fm=1;%信源最高频率fc=10;%载波中心频率T=5;%信号时长t=0:dt:T;mt=sqrt(2)*cos(2*pi*t);%信源%AM调制A=

35、2;s_am=(A+mt).*cos(2*pi*fc*t);B=2*fm;figure(1);subplot(311);plot(t,s_am);%绘制AM信号波形hold on;plot(t,mt,r-);%标示AM地包络title(AM调制信号);xlabel(t);%AM解调rt=s_am.*cos(2*pi*fc*t);rt=rt-mean(rt);f,rf=T2F(t,rt);t,rt=lpf(f,rf,2*fm);subplot(312);plot(t,rt);hold on;plot(t,mt/2,r-);title(相干解调后地波形与输入信号地比较);xlabel(t);%AM

36、调制信号频谱subplot(313);f,sf=T2F(t,s_am);psf=(abs(sf).2)/T;%调制信号地功率谱密度plot(f,psf);axis(-2*fc 2*fc 0 max(psf);title(AM信号功率);xlabel(f);仿真结果如图7所示.图7 AM调制信号4.3.2DSB-AM调制在DSB-AM系统中，已调信号地幅度正比与消息信号.这种调制通过使用乘法器完成，将消息信号吗m(t)与载波，如图8所示，表示为：WwghWvVhPE(3)图8 DSB-AM调制原理结构框图其中(4)是载波，而m(t)是消息信号.若以单频正弦信号调制为例，则调制信号u(t)经傅立叶

37、变换，可以得到DSB-AM信号地频域表示为：asfpsfpi4k(5)其中M(f)是m(t)地傅立叶变换.很明显可以看出，这种调制方式将消息信号地频谱进行了搬移，并在幅度上乘以Ac/2，传输带宽Br是消息信号带宽地两倍，即是：ooeyYZTjj1Br=2W (6)图9显示了一个典型地消息信号地频谱及其相对应地DSB-AM已调信号地频谱.图9 消息信号与DSB-AM已调信号地频谱已调信号地功率为(7)其中Pm是消息信号地功率.在DSB-AM通信系统中，信噪比SNR等于：(8)其中PR是接收到地功率（在接收端已调信号地功率），N0是噪声功率谱密度（假定为白噪声），W是信号噪声地带宽9.例5 对频率

38、为30Hz地余弦信号进行双边带幅度调制，载波频率为300Hz,并采用相干解调法实现解调.clear;clc;close all;fm=30;fc=300;T=1;t=0:0.001:T;m=2*cos(2*pi*fm*t);dsb= m.*cos(2*pi*fc*t);subplot(211);plot(t,dsb);title(DSB-AM调制信号);xlabel(t);r=dsb.*cos(2*pi*fc*t);r= r-mean(r);b=firl(40,0.01);rt=filter(b,1,r);subplot(212);plot(t,rt);title(相干解调后地信号);xlab

39、el(t);仿真结果如图10所示.图10 DSB_AM调制信号 4.4数字调制数字基带信号是低通型信号，其功率谱集中在零频附近，它可以直接在低通型信道中传输.实际信道很多是带通型地，数字基带信号无法直接通过带通型信道.因此，在发送端需要把数字基带信号地频谱搬移到带通信道地通带范围内，以便信号在带通型信道中传输，这个频谱地搬移过程称为数字调制.下面主要介绍了基于MATLAB地两种数字调制：二进制频移键控调制、正交幅度调制.BkeGuInkxI4.4.1二进制频移键控调制（2FSK）二进制频率调制是用二进制数字信息控制正弦波地频率，使正弦波地频率随二进制数字信息地变化而变化.由于二进制数字信息只有

40、两个不同地符号，所以调制后地已调信号有两个不同地频率f1和f2，对应数字信息”1”，对应数字信息”0”.图11为2FSK信号地产生方法及波形示例10.PgdO0sRlMo 图 11 2FSK信号地产生方法及波形相位不连续地2FSK信号，可看作是两个不同载频、时间交替地2ASK信号之和.其数学表示式可写为 (9)或 (10)相位不连续地2FSK信号可看作两个2ASK信号地叠加，因此其功率谱是两个2ASK信号功率谱之和. (11)2FSK信号地功率谱既有连续谱又有离散谱，离散谱位于两个载波频率处，连续谱分布在载波频率附近，若取功率谱第一个零点以内地成分计算带宽，显然2FSK信号地带宽为3cdXwc

41、km15(12)功率谱以 fc为中心对称分布,在h较小时功率谱呈现单峰，随着h地增大,f1和f2之间地距离增大，功率谱呈现双峰.h8c52WOngM为了节约频带，同时也能区分f1和f2，通常取此时2FSK信号带宽为(13)当时，2FSK地功率谱由双峰变成单峰，此时带宽为(14)频带利用率为(15)例6用MATLAB产生独立等概地二进制信源，对其进行2FSK调制，画出2FSK信号波形及功率谱图.程序如下： clc;clear;close all;M=2;N=200;nsample=8;fc=2;dt=1/fc/nsample;t=0:dt:N-dt;s=sign(randn(1,N);d=zer

42、os(fc*nsample,length(s+1)/2);v4bdyGiousd(1,:)=s;d=reshape(d,1,fc*nsample*length(s+1)/2);J0bm4qMpJ9g=ones(1,fc*nsample);dd=conv(d,g);sfsk=2*dd-1;fsk=cos(2*pi*fc*t+2*pi*sfsk(1:length(t).*t);sfft=abs(fft(fsk);sfft=sfft.2/length(sfft);subplot(2,1,1);plot(1:200,fsk(1:200);title(2FSK时域波形);subplot(2,1,2);p

43、lot(sfft);title(2FSK功率谱图);XVauA9grYP程序运行结果如图12所示. 图12 2FSK时域波形及其功率谱4.4.2 正交幅度调制正交幅度调制（QAM）信号可以看做是把信息序列an分离成俩路独立地基带数字波形，然后分别调制俩个正交载波地幅度和相位，最后把它们合并起来进行传输.这种信号地一个码元可以表示为11bR9C6TJscw(16)式中，k = 整数；和分别可以取多个离散值在信号表示式中，若qk值仅可以取p/4和-p/4，值仅可以取+A和-A，则此QAM信号就成为QPSK信号，如图13所示.pN9LBDdtrd 图13 QPSK信号矢量图所以QPSK信号是一种简单

44、地QAM信号，有代表性地QAM信号是16进制地，记为16QAM，它地矢量图示于图14中所示.DJ8T7nHuGTAk 图14 16QAM矢量图下面是用MATLAB实现16QAM信号及其功率谱图.例7 产生一个每码元4个样点地16QAM信号，采用升余弦脉冲成型，滚降系数为0.35，画出其功率谱图.QF81D7bvUA程序如下：clear;close all;clc;M=16;L=512;P=4;ini_phase=0;roll_off=0.35;a=2*randint(1,L,sqrt(M)-(sqrt(M)-1);b=2*randint(1,L,sqrt(M)-(sqrt(M)-1);x=a+

45、j*b;N=L*P;y=zeros(1,N);for n=1:Ny(n)=0;for k=1:Lt=(n-1)/P-(k-1);y(n)=y(n)+x(k)*(sin(pi*t+eps)/(pi*t+eps)*.(cos(roll_off*pi*t+eps)/(1-(2*roll_off*t)2)+eps);4B7a9QFw9hendendsfft=abs(fft(y);sfft=sfft.2/length(sfft);subplot(311);plot(real(x),imag(x),.);axis equal;title(16QAM信号星座图);ix6iFA8xoXsubplot(312)

46、;plot(1:length(sfft),sfft);title(16QAM基带信号功率谱图);for n=1:Nz(n)=y(n)*exp(j*2*pi*1*n/P);endsfft=abs(fft(z);sfft=sfft.2/length(sfft);subplot(313);plot(1:length(sfft),sfft);title(16QAM调制信号功率谱图);仿真结果如图15所示.图 15 16QAM信号地星座图及其功率谱通过使用MATLAB仿真可使16QAM信号及其功率谱图直观地反映出来.4.5 观察眼图码间干扰和噪声是影响数字基带通信系统性能地两个重要因素.码间干扰问题与系

47、统地发送滤波器、信道特性、接收滤波器特性等因素有关, 当系统总地特性为理想低通时, 可以完全消除码间干扰, 但是理想低通滤波器在现实地通信系统中是无法实现地.因此在现实地通信系统中码间干扰是一定存在地, 设计者只能让系统函数逼近理想低通来提高系统地性能, 降低误码率.为了对由码间干扰所引起地误率有一个直观地认识, 观察眼图是一个很好地方法.眼图可以借助于通信原理实验箱来观察, 也可以借助于MATLAB 地系统仿真功能来实现. 例8 下面就是运用matlab来仿真通信系统，模拟眼图地程序.wt6qbkCyDEFc=10; %载频Fs=40; %系统采样频率Fd=1; %码速率N=Fs/Fd;df

48、=10;numSymb=25;%进行仿真地信息代码个数M=2; %进制数SNRpBit=60;%信噪比SNR=SNRpBit/log2(M);seed=12345 54321;numPlot=15;%产生25个二进制随机码x=randsrc(numSymb,1,0:M-1);%产生25个二进制随机码figure(2)subplot(2,1,1)stem(0:numPlot-1,x(1:numPlot),bx);title(二进制随机序列)xlabel(Time);ylabel(Amplitude);%调制y=dmod(x,Fc,Fd,Fs,fsk,M,df);numModPlot=numPlo

49、t*Fs;t=0:numModPlot-1./Fs;subplot(2,1,2)plot(t,y(1:length(t),b-);axis(min(t) max(t) -1.5 1.5);title(调制后地信号)xlabel(Time);ylabel(Amplitude);%在已调信号中加入高斯白噪声randn(state,seed(2);y=awgn(y,SNR-10*log10(0.5)-10*log10(N),measured,dB);%在已调信号中加入高斯白噪声Kp5zH46zRkfigure(3)plot(t,y(1:length(t),b-);%画出经过信道地实际信号axis(m

50、in(t) max(t) -1.5 1.5);title(加入高斯白噪声后地已调信号)xlabel(Time);%相干解调figure(4)z1=ddemod(y,Fc,Fd,Fs,fsk/eye,M,df);title(相干解调后地信号地眼图)%带输出波形地相干M元频移键控解调figure(5)stem(0:numPlot-1,x(1:numPlot),bx);hold on;stem(0:numPlot-1,z1(1:numPlot),ro);hold off;axis(0 numPlot -0.5 1.5);title(相干解调后地信号原序列比较)legend(原输入二进制随机序列,相干

51、解调后地信号)xlabel(Time);ylabel(Amplitude); %非相干解调figure(6)z2=ddemod(y,Fc,Fd,Fs,fsk/eye/noncoh,M,df);title(非相干解调后地信号地眼图)%带输出波形地非相干M元频移键控解调figure(7)stem(0:numPlot-1,x(1:numPlot),bx);hold on;stem(0:numPlot-1,z2(1:numPlot),ro);hold off;axis(0 numPlot -0.5 1.5);title(非相干解调后地信号)legend(原输入二进制随机序列,非相干解调后地信号)xlabel(Time);ylabel(Amplitude);%误码率统计errorSym ratioSym=symerr(x,z1);figure(8)subplot(2,1,1)simbasebandex(0:1:5);title(相干解调后误码率统计)errorSym ratioSym=symerr(x,z2); subplot(2,1,2)simbasebandex(0:1:5);title(非相干解调后误码率统计)%滤除高斯白噪声Delay=3;R=0.5;PropD=0