MATLAB与在信号与系统中的应用.docVIP

. . MATLAB与在信号与系统中的应用 摘要：本论文通过MATLAB在信号与系统中的应用实例，探讨了MATLAB在信号与系统中的应用方法和技巧，对运用计算机辅助软件完成“信号与系统”课程的波形绘制，微分方程的求解，传递函数的分析、信号与系统分析具有较好的参考价值。 关键字：MATLAB 信号与系统 微分方程 传递函数 0.引言 “信号与系统”课程是一门实用性较强、涉及面较广的专业基础课，是电子信息类专业学生的必修课程。它是将电路分析的知识领域引入信号处理与传输领域的关键性课程，对后继专业课起着承上起下的作用。该课程的基本方法和理论大量运用于计算机信息处理的各个领域，特别是通信、数字语音处理、图像处理、数字信号分析等领域，应用更为广泛。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)) ezplot(f,[-5,4]), hold on,plot([0,0],[-1,1]), axis([-5,4,-1.1,1.1]),hold off 命令执行后，得出波形如图1。 图1 从图1可以清楚的看到f(t)的时域关系图，并且操作简单。 2.利用MATLAB解微分方程 在信号与系统中，连续时间系统处理连续时间信号，通常用微分方程来描述这类系统，也就是系统的输入与输出之间通过他们时间函数及其对时间t的各阶导数的线性组合联系起来。因此，在信号与系统中，求解微分方程对于研究连续时间系统的时域分析非常的重要。下面本文选择了一种简便的方法来处理这类问题，运用MATLAB的方法。 求解线性微分方程 y’’’+5y’’+4y’+7y=3u’’+0.5u’+4u 在输入u(t)为单位脉冲及单位阶跃信号时的解。 解：两边进行拉普拉斯变换（脉冲输入U(S)=1;单位阶跃U(S)=1/S） 求脉冲响应:A=[1 5 4 7];B=[3 0.5 4];[r ,p,k]=residue(B,A) 时域解:t=0:0.2:10 y=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 单击仿真按钮，然后返回MATLAB中输入一下命令，即可得到如图2-3 图2-3 传统求系统阶跃响应和冲激响应的方法都是对传递函数进行拉氏变换，再和激励函数的拉氏形式相乘，左后求反变换。不仅工程量大，而且得出的结果不直观；本文采用MATLAB建模的方式解决问题，不但简单，而且非常直观的反映了响应函数的特性。 3.根据传递函数求系统的零极点分布，单位冲激响应，单位阶跃响应以及幅频特性 在分析系统的静态和动态特性时，常常可以通过对传递函数进行分析，从而解决系统是否稳定等问题，然而传统的分析方法，常常要借助于人工计算，不仅工作量比较大，而且效率也比较低。而MATLAB解决了这一难题，本文通过调用MATLAB内部的函数进行计算，非常的高效的得到了零级点分布，单位冲激响应，单位阶跃响应以及幅频特性，非常的直观。 例：传递函数是 解：首先建立一个m文件，并且命名为“li.m” clear t=0:0.01:10; num=[1]; den=[1 2 2 1]; sys=tf(num,den) poles=roots(den) figure(1),pzmap(sys); %零级点分布 h1=impulse(num,den,t); figure(2),plot(t,h1); %单位冲激响应的波形 title(‘Impulse Response’) %加标题 h2=step(num,den,t); fig