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

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

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

3、e processing ,system simulation and so on in the integration.with matlab and simulink communications,sig-nal processing professional function library and professional toolbox being mature, they become more and more widely applied in communication theory research,algorit- hm design,system design,mode

4、ling simulation and performance analysis.this paper firstly introduces the matlab software and communication theory, then it introduces the application of matlab in communication theory by quoting exam- ples .it mainly uses matlab graphics and system simulation function such as symbol calculus opera

5、tions,drawing spectrum of periodic wave, simulated modulation,digital modulation,eye chat observation and simulinking sine signal.keywords matlab;communication theory;communication systerm;simulink目次1 引言.12 matlab简介.22.1 matlab的发展史 .22.2 matlab软件的组成.22.3 matlab的功能与特点.42.4 simulink的简介.43 通信原理简介.63.1通

6、信的概念.63.2通信系统模型.64 matlab在通信原理中的应用.84.1 matlab的符号积分运算.84.2 用matlab绘出周期波的频谱.94.3 模拟调制.124.4 数字调制.174.5 观察眼图.234.6 simulink仿真.27结论.28参考文献.29致谢.301 引言现代社会中通信系统功能越来越强，性能越来越高，构成越来越复杂,通信原理理论知识也在逐渐得到完善和深化，另一方面，要求通信系统技术研究和产品的开发周期缩短、成本降低、水平提高。这样尖锐对立的两个方面，只有通过使用强大的计算机辅助分析设计工具才可以实现。而大量优秀的计算机仿真软件的研究成功，对通信原理的发展起

7、到了极大的推动作用。我们所学的通信原理是一门理论性与实践性都很强的专业基础课，它主要讲述了模拟和数字通信系统的基本原理，在专业课程体系结构中起着非常重要的作用。然而通信原理这门课非常抽象，许多问题都需要在并不熟悉的频率域内进行分析、推理，算式较多。若将matlab融入其中，可简化计算过程，把计算结果以图的形式形象的显示出来，同时在通信系统的设计与仿真中，一些常用的模型都可以在matlab中找到相应的函数或函数的组合加以求解，使用户能够轻松地进行通信方面的分析与处理，用户也可以自己编程解决问题。本文首先介绍了matlab和通信原理的相关知识，而后介绍了matlab在通信原理中的应用，通过引例我们

8、会感受到matlab在通信原理中的强大作用，而且能够发现matlab应用于通信原理这一领域的良好的发展前景。2 matlab简介2.1 matlab的发展史1 matlab由matrix(矩阵)和laboratory（试验室）俩词的前3个字母组合而成。20世纪70年代后期，时任美国新墨西哥大学计算机系主任的cleve moler博士讲授线性代数课程时，发现应用其它高级编程语言极为不方便，为了减轻学生编程负担，于是cleve moler博士和他的同事构思并为学生设计了一组调用linpack和eispack库程序的“通俗易用”的接口，这就是用fortran编写萌芽状态的matlab.之后几年，ma

9、tlab作为免费软件在大学里被广泛使用，深受大学生的欢迎。matlab r系列是从2006年开始发布的，mathworks公司在技术层面上实现了一次飞跃。从那以后产品发布模式也发生了改变，每年的3月和9月进行俩次产品发布，产品及版本的命名方式为“r+年份+代码”，对应上下半年的代码分别是a和b.每一次发布都包含所有的产品模块，如产品的new feature、bug fixes和新产品模块的推出。例如，matlab r2008a是mathworks公司2008年3月份推出的产品。目前最新版本是matlab r2009b.2.2 matlab软件的组成matlab作为mathworks产品家族的核

10、心，它主要由5大部分组成，分别为matlab语言（the matlab language）、matlab工作环境（the matlab working environment）、matlab数学函数库（the matlab math library）、matlab应用程序接口（the matlab application interface）和图形句柄系统（handle graphics）2.matlab语言：matlab语言是一种以矩阵（matrix）和阵列（array)为基本编程单元，拥有完整的控制语句、数据结构、函数编写与调用格式和输出功能，具有面向对象程序设计特征的高级程序语言。使用m

11、atlab语言不但可以方便快捷地完成小规模的算法验证、程序开发和调试工作，而且可以进行大规模、高效的复杂应用程序设计。matlab工作环境：简单来说，matlab工作环境就是一系列实用工具的集合，它不但包括了各种操作工作空间中变量的工具和管理数据输入输出的方法，而且包括了开发调试m文件和matlab应用程序的集成环境。使用起来极为方便。matlab数学函数库：matlab数学函数库是大量的各种形式的数学函数和算法的集合，它不仅包括了最基本的初等函数，如sum、sine、cosine和复数运算等，而且包含了大量复杂的高级函数和算法，如贝塞尔（bassel)函数，快速傅里叶变换和矩阵求逆等。用户在

12、编写自己的matlab程序时，可以轻松地调用这些函数和算法，从而极大地方便了算法的开发。所有这些函数按类别分别存放在matlab工具箱下的8个子目录中，如表1所示。 表1 matlab数学函数库的分类和组织 目 录 名 函数功能 elmat 对矩阵和矩阵元素的操作 elfun 初等数学函数 matfun 矩阵函数数值线性代数 funfun 功能函数和ode求解 specfun 专门数学函数 sparfun 稀疏矩阵函数 datafun 数值分析和傅里叶变换 polyfun 插值和多边形近似matlab应用程序窗口：matlab的外部接口使得matlab可与外部设备和程序实现数据交互和程序移植，

13、可以扩充matlab强大的数值计算和图形显示功能，从而弥补了其执行效率较低的缺点，同时增强了其它应用程序进行软件开发的功能，提高了软件开发效率。matlab接口工具不仅使得matlab可以十分方便地与其他应用程序交换数据和信息，还实现了与其他程序函数和算法的交互。所以，通过matlab接口编程，可以充分利用现有资源，能更容易地编写出功能强大、结构简洁的应用程序。matlab主要提供了mex文件、mat文件、matlab计算引擎、com和dde、web服务、硬件接口和excel生成器、java生成器和net生成器等形式的接口。图形句柄系统：handle graphics为mathworks公司的

14、注册商标，是matlab的图形系统。它在包含了大量高级的2d和3d数据可视化、图形显示、动画生成和图像处理命令的同时，还提供了许多低级的图形命令，允许用户按照自己的需求显示图形和定制应用程序图形用户接口，既方便又灵活。具体的函数分为5大类，分别放置在工具箱下5个不同的目录内，如表2所示3。 表2 matlab图形函数的分类和组织 目 录 名 函 数 功 能 graph2d 二维图形函数 graph3d 三维图形函数 graphics 图形句柄函数 uitools 图形用户界面工具 specgraph 专门图形函数2.3 matlab的功能与特点 matlab是现在很流行的用于科学计算的软件工具

15、，它集数值分析、矩阵运算、信号处理、图像显示等多种功能于一体，提供了一个高性能的科学计算环境，因此在通信、自动化等领域都得到了广泛的应用。它最初是作为“矩阵理论”和“数据分析”等课程的计算工具，经过十几年的发展和完善，目前已成为世界各国在科学分析和计算领域的主要软件，并被ieee评述为国际公认的最优秀的科技应用软件。它的主要特点是4：（1）matlab可以用来解线性方程组、进行矩阵变换与运算、数据插值运算等，能使用户从繁杂的数学运算分析中解脱出来。（2）matlab中有许多高级的绘图函数，包括二维、三维、专用图形函数、图形句柄、用户图形界面工具等，利用这些函数可以轻松地完成各种图形的绘制和编辑

16、工作，实现计算结果和编程的可视化。（3）友好的用户界面及接近数学表达式的自然化语言，使学习者易于学习和掌握。（4）功能丰富的应用工具箱（如信号处理工具箱、通信工具箱、控制系统工具箱等）为用户提供了大量方便实用的处理工具。matlab的上述特点，使它深受工程技术人员及科技专家的欢迎，并很快成为应用学科计算机辅助分析、设计、仿真等领域不可缺少的基础软件。目前，国内很多理工院校已经或者正在把该软件作为学生必须掌握的一种软件。针对信号与系统课程内容的特点，利用matlab的信号处理工具箱和图形处理及数据可视化，教师可以将结论直接用图形来演示，从而让学生对抽象的概念和定理以及结论有直观的认识，并加深对一

17、些重要概念的理解。2.4 simulink简介 simulink是matlab提供的用于动态系统进行建模、仿真和分析的工具包。simulink提供了专门用于显示输出信号的模块，可以在仿真过程中随时观察仿真结果。同时，通过simulink的存储模块，仿真数据可以方便地以各种形式保存在工作区或文件中，供用户在仿真结束之后对数据进行分析和处理。另外，simulink把具有特定功能的代码组织成模块的方式，并且这些模块可以组织成具有等级结构的子系统，因此具有内在的模块化设计要基于matlab的通信系统仿真研究要求。基于上述优点，simulink成为一种通用的仿真建模工具，广泛应用于通信系统仿真、数字信号

18、处理、模糊逻辑、神经网络、机械控制和虚拟现实等领域。 根据输出信号与输入信号的关系，simulink提供3种类型的模块:连续模块、离散模块和混合模块。连续模块是指输出信号随着输入信号发生连续变化的模块，离散模块则是输出信号以固定间隔变化的模块。对于连续模块，simulink采用积分方式计算输出信号的数值，因此，连续模块主要涉及导数的计算及其积分。离散模块的输出信号在下一个抽样时刻到来之前保持恒定，这时候simulink只需以一定的间隔计算输出信号的数值。混合模块是根据输入信号的类型来确定输出信号类型的，它既能够产生连续输出信号，也能够产生离散输出信号。simulink包含有sinks（输出方式

19、）、source（输入源）、linear（线性环节）、nonlinear（非线性环节）、connections（连接与接口）和extra（其他环节）子模型库，而且每个子模型库中包含有相应的功能模，用户也可以定制和创建用户自己的模块。用simulink创建的模型可以具有递阶结构，因此用户可以采用从上到下或从下到上的结构创建模型。用户可以从最高级开始观看模型，然后用鼠标双击其中的子系统模块，来查看其下一级的内容，以此类推，从而可以看到整个模型的细节，帮助用户理解模型的结构和各模块之间的相互关系。在定义完一个模型后，用户可以通过simulink的菜单或matlab的命令窗口键入命令来对它进行仿真。菜

20、单方式对于交互工作非常方便，而命令行方式对于运行一大类仿真非常有用。采用scope模块和其他的画图模块，在仿真进行的同时，就可观看到仿真结果。除此之外，用户还可以在改变参数后来迅速观看系统中发生的变化情况。仿真的结果还可以存放到matlab的工作空间里做事后处理5。 模型分析工具包括线性化和平衡点分析工具、matlab的许多工具及matlab的应用工具箱。由于matlab和simulink的集成在一起的，因此用户可以在这两种环境下对自己的模型进行仿真、分析和修改。3 通信原理简介3.1 通信的概念 通信就是克服距离上的障碍，从一地向另一地传递和交换消息。消息是信息源所产生的，是信息的物理表现，

21、例如，语音、文字、数据、图形和图象等都是消息(message)。消息有模拟消息（如语音、图象等）以及数字消息（如数据、文字等）之分。所有消息必须在转换成电信号（通常简称为信号）后才能在通信系统中传输。所以，信号（signal）是传输消息的手段，信号是消息的物质载体。 相应的信号可分为模拟信号和数字信号，模拟信号的自变量可以是连续的或离散的，但幅度是连续的 ，如电话机、电视摄像机输出的信号就是模拟信号。数字信号的自变量可以是连续的或离散的，但幅度是离散的，如电船传机、计算机等各种数字终端设备输出的信号就是数字信号6。 通信的目的是传递消息，但对受信者有用的是消息中包含的有效内容，也即信息(inf

22、ormation) 。消息是具体的、表面的，而信息是抽象的、本质的，且消息中包含的信息的多少可以用信息量来度量。3.2 通信系统模型 通信系统就是传递信息所需要的一切技术设备和传输媒质的总和，包括信息源、发送设备、信道、接收设备和信宿(受信者) ，它的一般模型如图1所示。信息源发送端信道接受设备噪声源受信者图1 通信系统的一般模型各种通信系统中传输的消息可以分为离散消息和连续消息。离散消息也称为数字消息，它的消息状态是可数的或离散的。连续消息也称为模拟消息，其状态是连续变化的。按照信道中传输的消息是模拟消息还是数字消息可以把通信系统分为两类：模拟通信系统和数字通信系统，它们的模型如图2和图3所

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

24、达式进行n次微分运算例1 积分运算解:编制如下程序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)/tif2 = -1/2*(-1+cos(pi*t)2)/t, (1+sin(pi*t)-cos(pi*t)/t例2 微分运算解：编制如下程序clear;clc;syms x tf1=sin(x3+2*x2

25、+6)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的可视化结果进行分析。信号频谱是将信号源发出的信号强度按频率顺序展开，使其成为频率的函数，并考察变化规律，称为频谱分析。频谱分析主要分析信号是由哪些频率的正弦信号叠加得到的，以及得到这些正弦信号的振幅,信号的谱分析可以说就是计算信号的

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

27、 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 所需展开的最高谐波次数: );endt=5;if nargin5;n

28、n=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阵左右对称交换s1(1,k+1)=a_sym(1) %a_

29、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; %形成- n:n的变量s7=s6,s5(2:en

30、d); %形成- 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,80,0,0.25)end% - - - -

31、- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -function 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);y=y1- sym(heaviside(t) )*(4*t/a1);y=subs(y,a1,a);y=simple(y);%- - - - -

32、 - - - - - - - - - - - - - - - - - - - - - - - - -function 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.1 am调制am信号的时域表示式为 （1）频谱为 （2） 调制器模型如图5所示。图5 am调制器模型am的时域波形和频谱图如图6所示。 图6 am调制时、频域波形

33、am的调制效率最高为1/3；其优点：可以采用包络检波进行解调，实现起来容易，缺点为调制效率不高。例4 用matlab产生一个频率为1hz、功率为1w的余弦信源m（t），设载波频率为10hz，a=2，试画出：（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=2;s_am=(a+mt).*cos(2*pi*fc*t);b=2*fm;figure(1);subplot(31

34、1);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调制信号频谱subplot(313);f,sf=t2f(t,s_am);psf=(abs(sf).2)/t;%调制

35、信号的功率谱密度plot(f,psf);axis(-2*fc 2*fc 0 max(psf);title(am信号功率);xlabel(f);仿真结果如图7所示。 图7 am调制信号4.3.2 dsb-am调制在dsb-am系统中，已调信号的幅度正比与消息信号。这种调制通过使用乘法器完成，将消息信号吗m(t)与载波，如图8所示，表示为： (3) 图8 dsb-am调制原理结构框图其中 (4)是载波，而m(t)是消息信号。若以单频正弦信号调制为例，则调制信号u(t)经傅立叶变换，可以得到dsb-am信号的频域表示为：(5)其中m(f)是m(t)的傅立叶变换。很明显可以看出，这种调制方式将消息信号

36、的频谱进行了搬移，并在幅度上乘以ac/2，传输带宽br是消息信号带宽的两倍，即是：br=2w (6)图9显示了一个典型的消息信号的频谱及其相对应的dsb-am已调信号的频谱。图9 消息信号与dsb-am已调信号的频谱已调信号的功率为 (7)其中pm是消息信号的功率。在dsb-am通信系统中，信噪比snr等于： (8)其中pr是接收到的功率（在接收端已调信号的功率），n0是噪声功率谱密度（假定为白噪声），w是信号噪声的带宽9。例5 对频率为30hz的余弦信号进行双边带幅度调制，载波频率为300hz,并采用相干解调法实现解调。clear;clc;close all;fm=30;fc=300;t=1

37、;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(相干解调后的信号);xlabel(t);仿真结果如图10所示。 图10 dsb_am调制信号 4.4数字调制 数字基带信号是低通型信号，其功率谱集中在零频附近，它可以直接在低通型信道中传输

38、。实际信道很多是带通型的，数字基带信号无法直接通过带通型信道。因此，在发送端需要把数字基带信号的频谱搬移到带通信道的通带范围内，以便信号在带通型信道中传输，这个频谱的搬移过程称为数字调制。下面主要介绍了基于matlab的两种数字调制：二进制频移键控调制、正交幅度调制。4.4.1二进制频移键控调制（2fsk） 二进制频率调制是用二进制数字信息控制正弦波的频率，使正弦波的频率随二进制数字信息的变化而变化。由于二进制数字信息只有两个不同的符号，所以调制后的已调信号有两个不同的频率f1和f2，对应数字信息”1”，对应数字信息”0”。图11为2fsk信号的产生方法及波形示例10。 图 11 2fsk信号

39、的产生方法及波形相位不连续的2fsk信号，可看作是两个不同载频、时间交替的2ask信号之和。其数学表示式可写为 (9)或 (10) 相位不连续的2fsk信号可看作两个2ask信号的叠加，因此其功率谱是两个2ask信号功率谱之和。 (11)2fsk信号的功率谱既有连续谱又有离散谱，离散谱位于两个载波频率处，连续谱分布在载波频率附近，若取功率谱第一个零点以内的成分计算带宽，显然2fsk信号的带宽为 (12)功率谱以 fc为中心对称分布,在h较小时功率谱呈现单峰，随着h的增大,f1和f2之间的距离增大，功率谱呈现双峰。为了节约频带，同时也能区分f1和f2，通常取此时2fsk信号带宽为 (13) 当

40、时，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=zeros(fc*nsample,length(s+1)/2);d(1,:)=s;d=reshape(d,1,fc*nsample*length(s+1)/2);g=ones(1,fc*nsample);dd=conv(d,

41、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);plot(sfft);title(2fsk功率谱图);程序运行结果如图12所示。 图12 2fsk时域波形及其功率谱4.4.2 正交幅度调制 正交幅度调制（qam）信号可以看做是把信息序列an分离成俩路独立的基带数字波形，然后分别调制俩个正交载波的幅度和相位

42、，最后把它们合并起来进行传输。这种信号的一个码元可以表示为11 (16)式中，k = 整数；和分别可以取多个离散值在信号表示式中，若qk值仅可以取p/4和-p/4，值仅可以取+a和-a，则此qam信号就成为qpsk信号，如图13所示。 图13 qpsk信号矢量图所以qpsk信号是一种简单的qam信号，有代表性的qam信号是16进制的，记为16qam，它的矢量图示于图14中所示。 ak 图14 16qam矢量图下面是用matlab实现16qam信号及其功率谱图。例7 产生一个每码元4个样点的16qam信号，采用升余弦脉冲成型，滚降系数为0.35，画出其功率谱图。程序如下：clear;close

43、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+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);endendsfft=abs(fft(y

44、);sfft=sfft.2/length(sfft);subplot(311);plot(real(x),imag(x), . );axis equal;title(16qam信号星座图);subplot(312);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调制信号功率谱图);仿真结果

45、如图15所示。 图 15 16qam信号的星座图及其功率谱通过使用matlab仿真可使16qam信号及其功率谱图直观的反映出来。4.5 观察眼图码间干扰和噪声是影响数字基带通信系统性能的两个重要因素。码间干扰问题与系统的发送滤波器、信道特性、接收滤波器特性等因素有关, 当系统总的特性为理想低通时, 可以完全消除码间干扰, 但是理想低通滤波器在现实的通信系统中是无法实现的。因此在现实的通信系统中码间干扰是一定存在的, 设计者只能让系统函数逼近理想低通来提高系统的性能, 降低误码率。为了对由码间干扰所引起的误率有一个直观的认识, 观察眼图是一个很好的方法。眼图可以借助于通信原理实验箱来观察, 也可

46、以借助于matlab 的系统仿真功能来实现。 例8 下面就是运用matlab 来仿真通信系统，模拟眼图的程序。fc=10; %载频fs=40; %系统采样频率fd=1; %码速率n=fs/fd;df=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

47、),bx);title(二进制随机序列)xlabel(time);ylabel(amplitude);%调制y=dmod(x,fc,fd,fs,fsk,m,df);nummodplot=numplot*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);%在已调信号中加入高斯白噪声figure(3)plot(t,y(1:length(t),b-);%画出经过信道的实际信号axis(min(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);h