资源描述：
《MATLAB在控制系统仿真中的应用》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、第四章MATLAB在控制系统仿真中的应用控制系统的数学模型对于控制系统的研究具有重要意义。要对系统进行仿真研究，首先应建立系统的数学模型，在此模型的基础上建立系统的仿真模型，然后进行仿真，分析研究系统，并设计出相应的控制器对系统进行控制，使系统响应达到预期目标。本章重点：控制系统数学模型的基本描述方法控制系统的稳定性分析控制系统的根轨迹系统的频率响应线性系统时间响应分析一、控制系统数学模型的基本描述方法G(s)=G(s)=控制系统仿真中，主要用4种形式的数学模型：传递函数、零极点模型、结构图形式和状态方程模型。传递函数零极点状态方程结构图形式X=Ax+B
2、uY=Cx+Du4.1.1控制系统的传递函数描述用分子分母的系数构成两个向量，唯一的确定出来：num=[b1,b2,…,bn]den=[a1,a2,…,an]注意：构成分子，分母向量按降幂排列的顺序。传递函数的分子、分母均为多项式相乘的形式时，不能直接写出，可借助多项式乘法运算函数conv()来处理，以便获得分子、分母多项式向量。举例4.1.2控制系统的零极点模型将零点、极点及K值输入即可建立零极点模型。z=[-z,-z…,-z]p=[-p,-p…,-p]k=k对于给出的传函来说，分别对分子分母作因式分解，则可以得出系统的零极点模型。这可以通过求出分子，
3、分母多项式的根来实现。MATLAB中提供了多项式求根的函数，roots()调用格式：z=roots(a)其中：z—各个根所构成的向量a—多项式系数向量系统增益k即为原传递函数分子的最高项系数与分母最高项系数的比值。分别求出分子分母多项式的根，即可得到系统的零，极点模型。4.1.3控制系统状态方程模型只要将A，B，C，D几个矩阵输入进去即可。对于离散系统来说，也与上面类似。MATLAB还提供了由系统状态方程转换为系统零极点模型及传递函数模型的函数。[num,den]=ss2tf(A,B,C,D,iu)%iu—表示输入的序号（对多输入系统）[z,p,k]=s
4、s2zp(A,B,C,D,iu)%[z,p,k]—表示对第iu个输入信号的传递函数的零极点。[A,B,C,D]=tf2ss(num,den)[A,B,C,D]=zp2ss(z,p,k)4.1.4控制系统工具箱中的LTI对象为了避免对一个系统采用多个分离变量进行描述，新版本的控制系统工具箱，将LTI系统的各种描述封装成一个对象，即用一个变量来描述。在控制系统工具箱中，有以上讲述的三种对象，即ss对象，tf对象和zpk对象。每种系统模型的生成和模型间的转换均可以通过一个函数来实现dss()函数：生成系统的状态空间模型。filt()函数：生成DSP（数字信号处
5、理）形式的离散传递函数ss()函数：生成状态空间模型，或者将传递函数及零极点模型转换成状态空间模型tf()函数：生成传递函数模型，或将零极点模型及状态空间模型转换成传递函数模型zpk()函数：生成零极点模型或者将其他模型转化成零极点模型4.1.5控制系统的结构图描述及转换注意：复杂的结构图建模可以用SIMUKINK简单的结构图描述及典型连接串联并联负反馈series()Parallel()Feedback()二、控制系统的稳定性分析计算机仿真中，控制系统稳定性的分析方法有：①求闭环特征方程的根；②化为零极点模型，看极点是否在s右半平面；③对状态空间形式（
6、闭环），求A阵的特征值eig(A)；④用Lyapunov方程求解。在自控原理中，均采用间接方法劳斯判据，奈氏判据等，但由于在MATLAB中很容易地求解多项式方程。因此，我们可以直接求出特征方程的根，可用多项式求根函数roots()。控制系统的稳定性分析numg=1;deng=[11223];numf=1;denf=1;(num,den)=feedback(numg,deng,numf,denf,-1);roots(den)例：三、控制系统的根轨迹控制系统工具箱中提供了系统根轨迹绘制及分析函数，还有一个GUI分析工具。根轨迹的绘制及分析函数：rlocfin
7、d():计算给定根的根轨迹增益格式：①[k,poles]=rlocfind(sys)②[k,poles]=rlocfind(sys,p)rlocus()函数：功能为求系统根轨迹。格式：rlocus(sys)---计算SISO开环LTI对象的根轨迹，增益自动选取。rlocus(sys,k)---显式设置增益。[r,k]=rlocus(sys)---返回系统的增益k（向量）和闭环极点r（向量），即对应于增益的闭环极点为r(i)sgrid()函数：绘制连续时间系统根轨迹和零极点图中的阻尼系数和自然频率栅格线控制系统的根轨迹num=[251];den=[123]
8、;sys=tf(num,den);rlocus(sys);[r,k]=rlocu
