matlab与在信号与系统中的应用

matlab与在信号与系统中的应用

ID:25692236

大小:158.50 KB

页数:11页

时间:2018-11-22

matlab与在信号与系统中的应用_第1页
matlab与在信号与系统中的应用_第2页
matlab与在信号与系统中的应用_第3页
matlab与在信号与系统中的应用_第4页
matlab与在信号与系统中的应用_第5页
资源描述:

《matlab与在信号与系统中的应用》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、MATLAB与在信号与系统中的应用摘要:本论文通过MATLAB在信号与系统中的应用实例,探讨了MATLAB在信号与系统中的应用方法和技巧,对运用计算机辅助软件完成“信号与系统”课程的波形绘制,微分方程的求解,传递函数的分析、信号与系统分析具有较好的参考价值。关键字:MATLAB信号与系统微分方程传递函数0.引言“信号与系统”课程是一门实用性较强、涉及面较广的专业基础课,是电子信息类专业学生的必修课程。它是将电路分析的知识领域引入信号处理与传输领域的关键性课程,对后继专业课起着承上起下的作用。该课程的基本方法和理论大量运用于计算机信息处理的各个领域,特别是通信、数字语音处理、

2、图像处理、数字信号分析等领域,应用更为广泛。MATLAB最为一个辅助类的软件可以很好的完成数值计算、信号与系统分析的可视化建模及仿真调试。同时MATLAB中的SIMULINK也是非常好的建模工具,极大的简化了仿真和计算的工作量。1.用MATLAB绘制信号的波形本文运用MATLAB绘图,很快速和简便的得到响应函数的波形图,通过分析响应的波形可以对信号的分析有更深刻的认识。用MATLAB绘制f(t)=e(t+3)-2e(t),其中e(t)为step函数解:在MATLAB窗口中输入如下命令即可实现f(t)的绘制f=sym('heaviside(t+3)-2*heaviside(t

3、)')ezplot(f,[-5,4]),holdon,plot([0,0],[-1,1]),axis([-5,4,-1.1,1.1]),holdoff命令执行后,得出波形如图1。图1从图1可以清楚的看到f(t)的时域关系图,并且操作简单。2.利用MATLAB解微分方程在信号与系统中,连续时间系统处理连续时间信号,通常用微分方程来描述这类系统,也就是系统的输入与输出之间通过他们时间函数及其对时间t的各阶导数的线性组合联系起来。因此,在信号与系统中,求解微分方程对于研究连续时间系统的时域分析非常的重要。下面本文选择了一种简便的方法来处理这类问题,运用MATLAB的方法。求解线性

4、微分方程y’’’+5y’’+4y’+7y=3u’’+0.5u’+4u在输入u(t)为单位脉冲及单位阶跃信号时的解。解:两边进行拉普拉斯变换(脉冲输入U(S)=1;单位阶跃U(S)=1/S)求脉冲响应:A=[1547];B=[30.54];[r,p,k]=residue(B,A)时域解:t=0:0.2:10y=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t)+r(3)*exp(p(3)*t);plot(t,y)得到的波形图如图2-1图图2-1求阶跃响应:打开MATLAB中的simulink模块,建立一个模型文件,命名为”li.mdl”。如图2-2图2-2单击

5、仿真按钮,然后返回MATLAB中输入一下命令,即可得到如图2-3图2-3传统求系统阶跃响应和冲激响应的方法都是对传递函数进行拉氏变换,再和激励函数的拉氏形式相乘,左后求反变换。不仅工程量大,而且得出的结果不直观;本文采用MATLAB建模的方式解决问题,不但简单,而且非常直观的反映了响应函数的特性。3.根据传递函数求系统的零极点分布,单位冲激响应,单位阶跃响应以及幅频特性在分析系统的静态和动态特性时,常常可以通过对传递函数进行分析,从而解决系统是否稳定等问题,然而传统的分析方法,常常要借助于人工计算,不仅工作量比较大,而且效率也比较低。而MATLAB解决了这一难题,本文通过调

6、用MATLAB内部的函数进行计算,非常的高效的得到了零级点分布,单位冲激响应,单位阶跃响应以及幅频特性,非常的直观。例:传递函数是解:首先建立一个m文件,并且命名为“li.m”cleart=0:0.01:10;num=[1];den=[1221];sys=tf(num,den)poles=roots(den)figure(1),pzmap(sys);%零级点分布h1=impulse(num,den,t);figure(2),plot(t,h1);%单位冲激响应的波形title(‘ImpulseResponse’)%加标题h2=step(num,den,t);figure(3

7、),plot(t,h2);%单位阶跃响应波形[H,w]=freqs(num,den);figure(4),plot(w,abs(H));%幅频特性图xlabel(‘omega’)title(‘magnitudeResponse’)%加标题运行结果是Transferfunction:1---------------------s^3+2s^2+2s+1poles=-1.0000-0.5000+0.8660i-0.5000-0.8660i零极点分布如图3-1所示图3-1单位冲激响应的波形如图3-2所示图图3-2单位阶

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。