基于MATLAB的信号仿真分析的设计___课设说明书

1、唐山学院信号仿真分析实践题目基于MATLAB的信号仿真分析的设计系（部）信息工程系班级姓名学号指导教师2014年6月16日至6月20日共1周如14年6月20日信号仿真分析实践任务书一、设计题目、内容及要求设计题目：1、基于MATLAB的信号仿真分析的设计设计内容：(1)MATLAB程序入门和基础应用。学习Matlab软件的基本使用方法，了解Matlab的数值计算、符号运算、可视化功能以及Matlab程序设计入门。(2)连续时间LTI系统的时域分析。利用MATLAB求解连续系统的零状态响应、冲激响应及阶跃响应，利用MATLAB实现连续信号卷积的方法。(3)连续时间LTI系统的频域分析。利用MAT

2、LAB实现连续时间信号傅里叶变换，分析LTI系统的频域特性和LTI系统的输出响应。(4)连续时间LTI系统的复频域分析。利用MATLAB进行部分分式展开,分析LTI系统的特性，利用MATLAB进行Laplace正、反变换。设计要求：1、根据题目要求进行基于MATLAB的信号仿真分析的总体设计。2 .完成算法具体部分的设计。(1)仿真原理图。(2)仿真分析说明。3 .仿真程序的设计。(1)连续时间系统时域和变换域相关特性的仿真分析。(2)完整源程序。4.书写设计说明书。、设计原始资料1郑君里.信号与系统(第三版).高等教育出版社，2011三、要求的设计成果(课程设计说明书、设计实物、图纸等)设计

3、结果能正确仿真演示设计说明书一份(包括总体设计、算法原理图及说明、系统仿真演示、源程序清单等)四、进程安排周周周周周资料收集利用MATLAB完成信号仿真分析的设计程序调试书写仿真分析实践说明书答辩五、主要参考资料1薛定宇，陈阳泉.基于MATLAB/Simulink的系统仿真技术与应用，清华大学出版社，20022钟麟，王峰.MATLAB仿真技术与应用教程，国防工业出版社，20053楼顺天，李博前.基于MATLAB的系统分析与设计一信号处理，西安电子科技大学出版社，19984吴湘淇.信号、系统与信号处坪的软硬件实现.电子工业出版社.2002指导教师(签名)：教研室主任(签名)：课程实践成绩评定表出

4、勤情况出勤天数缺勤天数成绩评定出勤情况及设计过程表现（20分）课程实践答辩（20分）设计成果（60分）总成绩（100分）提问（答辩）问题情况综合评定指导教师签名：年月日目录1引言11.1 MATLAB勺发展历史11.2 MATLAB勺系统结构11.3 MATLAB勺主要特点11.4 MATLAB勺影响22 .总体设计32.1 设计框图32.2 设计思想33 .详细设计53.1 基础运算53.2 连续时间LTI系统时域仿真分析73.2.1 设计方法与步骤73.2.2 连续时间系统零状态响应的数值求解83.2.3 卷积的计算83.2.4 连续时间系统冲激响应和阶跃响应的求解93.3 连续时间LTI

5、系统频域仿真分析133.3.1 实验目的133.3.2 设计框图143.3.3 系统的频率特性143.3.4 连续时间信号傅里叶变换的数值计算方法153.3.5 涉及的MATLAB数153.4 连续时间LTI系统复频域仿真分析173.4.1 实验目的173.4.2 利用MATLAEft行部分分式展开183.4.3 连续时间信号的拉普拉斯变换183.4.4 连续系统的复频域分析183.4.5 系统频率特性分析194 .总结24参考文献25唐山学院课程实践1引言MATLAB是一种面向科学与工程计算的高级语言,由于其集成了许多领域的工具箱，因此又被称之为巨人肩上的工具本文以该软件的MATLAB7.0

6、为例，在简单介绍MATLAB的基础上，详细阐述了其仿真功能在信号处理中的应用。1.1 MATLAB的发展历史MATLAB语言是由美国的CleverMoler博士于1980年开发的；设计者的初衷是为解决线性代数”课程的矩阵运算问题；取名MATLAB即MatrixLaboratory矩阵实验室的意思；1984年，MathWorks公司成立，MATLAB以商品形式出现，随着市场的推广，以其良好的开放性和运行的可靠性，成功地应用于各工程学科的研究领域。1.2 MATLAB的系统结构MATLAB系统结构有三个层次、八个部分构成，每一个部分完成不同的功能，主要包括：(1) MATLAB主包(2)工具箱(T

7、oolboxes)(3)编译器(Compiler)(4)建模仿真(Simulink)(5)模块集(Blockset)(6)实时仿真(Real-TimeWorksho)(7)状态流程(Stateflow)(8)状态代码生成(StateflowCoder)1.3MATLAB的主要特点(1)简单易学：与C语言几乎一致。(2)编程简洁、效率高；语言规则不强，无需编译。(3)强大而简单的绘图功能；二、三维绘图，直观展示结果唐山学院课程实践(4)扩展性强一一丰富的工具箱：MATLAB主工具箱、符号数学工具箱、SIMULINK仿真工具箱、控制系统工具箱、信号处理工具箱、图象处理工具箱、通信工具箱、系统辨识工

8、具箱、神经元网络工具箱、金融工具箱等。1.4MATLAB的影响(1) matlab在教学中的应用在大学里，诸如应用代数、数理统计、自动控制、数字信号处理、动态系统仿真等课程的教科书都把MATLAB作为内容。在欧美是九十年代教科书与旧版书籍的区别性标志。越来越多国内教科书也加入了Matlab相关内容。MATLAB是攻读学位的大学生、硕士生、博士生必须掌握的基本工具。(2) MATLAB行业运用相关图书：MATLAB与选煤/选矿数据处理、MATLAB在工程数学上的应用、电磁场数值计算法与MATLAB实现、数学物理方程的MATLAB解法与可视化、高等应用数学问题的MATLAB求解、MATLAB有限元

9、分析与应用、MATLAB7及工程问题解决方案、MATLAB及在电子信息类课程中的应用、MATLAB在电气工程中的应用等等。(3) MATLAB在本专业中的应用不管世界如何改变，MATLAB在通信专业的地位基本不变。(通信专业计算器”)几乎所有通信与信号处理类课程：高等数学赠通物理峨性代数既率论与数理统计算方法(数值分析)号与系统t路！信原理砌机信号处理数字信号处理峨字图像处理毗纤通信移动通信！信网唐山学院课程实践2 .总体设计总体设计的目标是：(1)MATLAB程序入门和基石运用。学习MATLAB软件的基本使用方法，了解MATLAB的数值计算、符号运算、可视化功能以及MATLAB程序入门设计。

10、(2)连续时间LTI系统的时域分析。利用MATLAB求解连续系数的零状态响应、冲激响应及阶跃响应，利用MATLAB实现连续信号卷积的方法。(3)连续时间LTI系统的频域分析。利用MATLAB实现连续时间信号傅里叶变换，分析LTI系统的频域特性和LTI系统的输出响应。(4)连续时间LTI系统的复频域分析。利用MATLAB进行Laplace正、反变换。2.1设计框图图212.2设计思想唐山学院课程实践LTI系统各种分析方法的理论基础是信号的分解特性与系统的线性、时不变特性，其出发点是：激励信号可以分解为若干基础信号单元的线性组合；系统对激励所产生的零状态响应是系统对各基本信号单元分别激励下响应的叠

11、加。利用MATLAB进行MATLAB程序入门和基础运用、连续时间LTI系统的时域分析、连续时间LTI系统的频域分析和连续时间LTI系统的复频域分析。唐山学院课程实践3 .详细设计3.1 基础运算设计步骤：在MATLAB命令窗口输入运算命令进行运算下面一道题是用MATLAB进行加减乘除”运算1 1利用MATLAB计算A=-2-5i解在MATLAB命令窗口键入如下语句：A=-2-i*5运行A=-2.0000-5.0000i运行图如图31所小图31下面一道题是利用MATLAB表示一个复常量极坐标12一个复常量极坐标表示解在MATLAB命令窗口键入如下语句：B=3*exp(i*pi/4)运行B=唐山学

12、院课程实践2.1213+2.1213i运行图如图32所小图32下面一道题是利用MATLAB进行向量和逗号运算13向量运算和逗号解在MATLAB命令窗口键入如下语句：A=4,5,8,9运行A=4运行图如图33所小图33唐山学院课程实践下面一道题是利用MATLAB进行分号运算14分号运算解在MATLAB命令窗口键入如下语句：B=8;6;8;1B=8681运行图如图34所小图343.2连续时间LTI系统时域仿真分析实验目的：(1)学会运用MATLAB符号求解连续系统的零输入响应；(2)学会运用MATLAB数值求解连续系统的零状态响应；(3)学会运用MATLAB求解连续系统的冲击响应和阶跃响应;(4)

13、学会运用MATLAB卷积积分法求解系统的零状态响应。3.2.1设计方法与步骤一般的连续时间系统分析有以下几个步骤：唐山学院课程实践求解系统的零输入响应；求解系统的零状态响应；求解系统的全响应；分析系统的卷积；画出它们的图形。下面以具体的微分方程为例说明利用MATLAB软件分析系统的具体方法3.2.2连续时间系统零状态响应的数值求解我们知道，LTI连续系统可用如下所示的线性常系数微分方程来描述，NM、ay(i)(t)八bf(t)i=0jJ3例如，对于以下方程：-一一一一_一一a3y(t)a2y(t)a1y(t)a0y(t)=b3f(t)b2f(t)bj(t)bf(t)可用a=a3,82,81,8

14、0,b=b3,b2,bI,b3,输入函数u=f(t)，得出它的冲击响应h,再根据LTI系统的零状态响应y(t)是激励u(t)与冲击响应h(t)的卷积积分。注意，如果微分方程的左端或右端表达式中有缺项，则其向量a或b中的对应元素应为零，不能省略不写，否则出错。求函数的零状态响应y(t)5y(t)4y(t)=2f(t)-4f(t)及初始状态yzs(0_)=yZs(0_)=0o输入函数f(t)=sin(3*t)+cos(2*t)。建模先求出系统的冲击响应，写出其特征方程2_一254:0求出其特征根为p和p,及相应的留数r,r；则冲击响应为Pitp2th(t)-rer2e输入y(t)可用输入u(t)与

15、冲击响应h(t)的卷积求得。3.2.3卷积的计算连续时间信号fi(t)和22。)的卷积运算可用信为的分段求和来实现，即:f(t)=fi(t)*f2(t)=.fi(t)f2(t-)d=l呵%fi(k：)f2(t-k)：k一：J如果只求当t=n(n为整数)时f(t)的值f(n一),则上式可得：唐山学院课程实践oOOCf(n.：)=、f1(k；)f2(t-k.).:=.:、f1(k.)f2(n-k).-kk=_二二QO式中的Zfi(k&f2(nA-k)网实际上就是连续时间信号fi(t)和f2(t)经等时k二问间隔均匀抽样的离散序列(2和f2(Ak)的卷积和。当足够小时，f(nA)就是卷积积分的结果一

16、一连续时间信号f(t)的较好数值近似。下面是利用MATLAB实现连续信号卷积的通用程序conv(),该程序在计算出卷积积分的数值近似的同时，还绘制出f(t)的时域波形图。应注意，程序中是如何设定f(t)的时间长度。MATLAB程序：f1=input(输入函数f1=);f2=input(输入函数f2=);dt=input(dt=);y=conv(f1,f2);plot(dt*(1:length(y)-1),y);gridon;title(卷积);xlabel(t);ylabel(f1*f2)程序运行结果：输入以下数据：f1=sin(3*t)f2=cos(3*t+2)dt=0.01得出图形如图35

17、所示：图353.2.4连续时间系统冲激响应和阶跃响应的求解唐山学院课程实践在MATLAB中，对于连续LTI系统的冲激响应和阶跃响应的数值解，可分别用控制系统工具箱提供的函数impulse和step来求解。其语句格式分别为：y=impulse(sys,t);y=step(sys,t);其中，t表示计算机系统响应的时间抽样点向量，sys表示LTI系统模型。下面我们用几个例子来验证以上知识点(1)求零输入相应：求齐次微分方程在给定初始条件下的零输入相应，y(t)+4y(t)=0,y(0_)=1,y,(0_)=1。解在MATLAB命令窗口键入如下语句：eq=D2y+4*y=0;con=y(0)=1,D

18、y(0)=1;y=dsolve(eq,con);y=simplify(y)y=1/2*sin(2*t)+cos(2*t)ezplot(y),grid零输入响应图如图36所示图36(2)求零状态响应：已知系统y(t)+4y(t)=x(t),输入x(t)为cos(t)u(t),求系统的零状态响应。解在MATLAB命令窗口键入如下语句eq=D2y+4*y=cos(t)*heaviside(t);con=y(-0.01)=0,Dy(-0.01)=0;y=dsolve(eq,con);y=simplify(y)10唐山学院课程实践1/3*heaviside(t)*(cos(t)-cos(2*t)ezpl

19、ot(y,0,10)零状态响应图如图37所示图37(3)求冲击响应和阶跃响应：已知某LTI系统的微分方程为V,，(t)+2y(t)+32y(t)=f国珊1MATLA(t),命令绘出0三t三4范围内系统的冲激响应h(t)和阶跃响应g(t)0解在MATLAB命令窗口键入如下语句t=0:0,001:4;sys=tf(1,16,1,2,32); h=impulse(sys,t); g=step(sys,t); subplot(211)plot(t,h),gridonxlabel(Time(sec),),ylabel(,h(t)title(冲激响应,)冲激响应图如图38(1)所示11唐山学院课程实践图3

20、8(1)subplot(212)plot(t,g),gridonxlabel(Time(sec),ylabel(g(t) title(阶跃响应)阶跃响应图如图3-8(2)所示图38(2)12唐山学院课程实践(4)求卷积：求信号fl=t和f2(t)=sin的卷积的解析解解在MATLAB命令窗口键入如下语句 symssttao; f1=t;f2=sin(t);g=subs(f1,t,tao)*subs(f2,t,(t-tao); f=int(g,tao,0,t)t-sin(t)运行图如图39所小图393.3连续时间LTI系统频域仿真分析3.3.1 实验目的(1) .掌握连续时间信号傅立叶变换和傅立

21、叶逆变换的实现方法，以及傅立叶变换的时移特性，傅立叶变换的频移特性的实现方法；(2) .了解傅立叶变换的频移特性及其应用；(3) .掌握函数fourier和函数ifourier的调用格式及作用；(4) .掌握傅立叶变换的数值计算方法，以及绘制信号频谱图的方法。13唐山学院课程实践3.3.2设计框图设计框图如图310所示主窗口门信号尺度变换邳号信彳单边指数信号频移特性号信nJ门信号频域拉氏变换拉氏变换幅度图图3103.3.2系统的频率特性14时移特性即唐山学院课程实践连续LTI系统的频率特性又称频率响应特性，是指系统在正弦信号激励下的稳态响应随频率变化的情况，又称系统函数H侬对于一个零状态的线性

22、系统,如图311所小。Xj)Hj)Y(W)图311LTI系统框图Y(j)其系统函数H(o)定义为式H(j)=-X(j)中，X)为系统激励信号的傅里叶变换，Y侬)为系统在零状态条件下输出响应的傅里叶变换。系统函数H侬)反映了系统内在的固有特性，它取决于系统自身的结构及组成系统元器件的参数，与外部激励无关，是描述系统特性的一个重要参数。H(s)是缶的复函数，可以表示为H侬=H(ge我曾。其中，|H(o)随口变化的规律称为系统的幅频特性：中侬)随与变化的规律称为系统的相频特性。频率特性不仅可用函数表达式表示，还可以随频率f变化的曲线来表示。当频率特性曲线采用对数坐标时，又称为波特图。算法理论依据:3

23、.3.3 连续时间信号傅里叶变换的数值计算方法-bo8F(j)=_ftejtdt=lim%f(n)e-jn当f为时限信号时，或和近似的看做时限信号时，式中的n取值可认作是有限的，k=N；k设为N,则可得,编程中需要注意的是：要正确生成信号j”e0N1F(k)=fZf(ntJe-jwknT,0=k symst; f=fourier(exp(-2*abs(t); ezplot(f);傅里叶变换图如图3-12所示口0。刍玲10&号京口目r图312FileEditViewInsertToolsDesktopWindow3t(2)利用MATLAB绘出特性曲线：试回出f(t)=2AU(t)的波形及其频录3

24、e特性曲线。16唐山学院课程实践解在MATLAB命令窗口键入如下语句: symstvwf f=2/3*exp(-3*t)*sym(Heaviside(t); F=fourier(f); subplot(2,1,1); ezplot(f); subplot(2,1,2); ezplot(abs(F); ，、一22t信号f(t)=2eU(t)的波形及其幅频特性曲线如图3-13所小3e图3133.4连续时间LTI系统复频域仿真分析3.4.1 实验目的17唐山学院课程实践(1)利用MATLAB进行部分分式展开；(2)利用MATLAB进行Laplace正、反变换；(3)分析LTI系统的特性。3.4.2

25、利用MATLAB进行部分分式展开为了方便LAPLACE反变换，先对F(s)进行部分分式展开。根据F(s)分为具有不同极点的部分分式展开和具有多重极点的部分分式展开。分别讨论。(1)不同极点的部分分式展开：F(s)=B(s)/A(s)=num/den=(b0*sAn+b1*sA(n-1)+.+bn)/(sAn+a1*sA(n-1)+.an)在matlab行向量中，num和den分别表示传递函数分子和分母的系数num=b0b1.bnden=a0a1.an将求出F(s)部分展开式的留数，级点，余项r,p,kB(s)/A(s)=r(1)/(s-p(1)+r(2)/(s-p(2)+.+r(n)/(s-p

26、(n)+k(s)(2)多重极点部分分式展开：r(j)/(s-p(j)+r(j+1)/(s-p(j+1)A2+.r(j+m-1)/(s-p(j)Am3.4.3 连续时间信号的拉普拉斯变换连续时间信号的拉普拉斯正变换和逆变换分别为：二1j二7F(s)二一f(t)e,tf(t)=一F(s)estds_2j-j-Matlab的符号数学工具箱(SymbolicMathToolbox)提供了能直接求解拉普拉斯变换和逆变换的符号运算函数laplace()和ilaplace()。下面举例说明两函数的调用方法。3.4.4 连续系统的复频域分析若描述系统的微分方程为：NM、aiy(i)(t)4bjf(t)i=0j

27、=0则系统函数为:18唐山学院课程实践H(s)=Y(s)F(s)MVbjSj二-N-idsizzPB(s)A(s)其中，。通过分析系统函数H(s)的零、极点分布，可以掌握系统的特性。若H(s)的所有极点都分布于左半开s平面，则该因果系统是稳定的；否则，该因果系统是不稳定的。在MATLAB中，可利用多项式求根函数root()来确定系统函数H(s)的零、极点。(1)系统稳定性分析NMA(s)=、aisi,B(s)=、bjsji0j=0A(s)、B(s)是关于s的多项式，由其系数分别构成向量：a=aN,aN,ai,a0,b=/凤,.也息若执行如下命令，即可完成H(s)的零、极点的计算：p=roots

28、(a)%求人缶)的根(H(s)的极点)，并将结果返回给向量pq=roots(b)%求B(s)的根(H(s)的零点)，并将结果返回给向量q3.4.5系统频率特性分析对于因果系统来说，若系统函数H(s)的收敛域包含s平面的虚轴(j轴)，则其频率响应H(j/在，且H(j(o)=H(s)smy因此，分析系统频率特性之前，要先对系统的稳定性作分析。卜面让我们用几个例子来验证以上知识点:(1)利用MATLAB求解拉普拉斯变换(1) fi(t)=et；(t)(2)f2(t)=te：解在MATLAB命令窗口键入如下语句：symstf1=sym(exp(-2*t)*Heaviside(t);%求f1(t)的拉普

29、拉斯变换F1=laplace(f1)F1=1/(s+2)f2=sym(t*exp(-t)*Heaviside(t);19唐山学院课程实践F2=laplace(f2)F2=1/(s+1)A2运行图如图314所小图314(2)利用MATLAB拉普拉斯逆变换求解：若系统的系统函数为。H(s)=1s23s2Res-1求冲激响应h(t)。(拉普拉斯逆变换)解在MATLAB命令窗口键入如下语句H=sym(1/(sA2+3*s+2);h=ilaplace(H)20唐山学院课程实践2*exp(-3/2*t)*sinh(1/2*t)运行图如图315所示图315利用MATLAB进行系统稳定的判定:已知一因果系统的

30、系统函数为:H(s)=432s2s-3s2s1s2-4画出H(s)的零、极点分布图，并判断系统的稳定性。解在MATLAB命令窗口键入如下语句a=12-321;b=10-4;p=roots(a)21唐山学院课程实践-3.13000.9口后0.4口n0-0.2-0.4nh-0.9iI_用L一1寓1-3-2-1012实部H的等极点图图3160.7247+0.6890i0.7247-0.6890i-0.3195q=roots(b)q=2.0000-2.0000holdonplot(real(p),imag(p),x);plot(real(q),imag(q),o);title(H(s)的零极点图);gridon;ylabel(虚部);xlabel(实部)如图316所示由运行结果如图3-15所示。由图可知，H(s)的零点为q1=2.0,q2=-2.0;极点为Pi=-3.13,p2=0.7247+0.6890i,p3=0.7247-0.6890i,p4=-0.3195。H(s)的极点P2、P3处于右半开s平面，故系统