控制系统仿真应用.docx

1、控制系统仿真应用MATLAB在控制系统仿真中的应用从本章开始，正式进入到有关MATLAB在控制系统仿真的领域中去。我们知道，对控制系统进行仿真，首先要建立系统的数学模型，这是计算机仿真的基础。我们下面先介绍控制系统数学模型的描述方法。第一节 控制系统数学模型的基本描述方法在控制系统仿真中，主要用4种形式的数学模型：传递函数.、零极点模型、结构图形式和状态方程模型。这些模型之间存在着内在的等效关系。在不同的场合下可能使用的模型形式要求不同，需要了解模型之间的转换方法。这一节主要介绍它们的MATLAB实现。一控制系统的传递函数描述：（在MATLAB中称tf模型）对系统的微分方程在零初始条件下做拉氏

2、变换，则可得系统的传递函数(SISO系统)：G(s)= 对线性时不变（线性定常）系统（LTI）来说，a、b均为常数a10。前面讲过多项式的表示方法。这里分子分母都为多项式，可将分子分母分别表示出来。即用分子分母的系数构成两个向量，唯一的确定出来： num=b,b,bden=a,a,a注意：构成分子，分母向量按降幂排列的顺序。这只是一个简单的特例，很多时候，传递函数的分子、分母均为多项式相乘的形式，如：不能直接写出，可借助多项式乘法运算函数conv( )来处理，以便获得分子、分母多项式向量。conv( )函数的调用方式为：c=conv(a,b)a,b各表示一个多项式，c表示a和b的乘积多项式，此

3、函数允许嵌套使用：c=conv(a,conv(b,e)，则上例的G(s)可用下面的语句来输入：num=4*conv(1,2,conv(1, 6, 6,1, 6, 6);den=conv(1,0,conv(1,1,conv(1,1,conv(1,1,1,3,2,5);还可以进一步地编写一个convs( )函数来一次性地求出若干个（十个）多项式的连乘积：function a=convs(a1,a2,a3,a4,a5,a6,a7,a,8,a9,a10)a=a1;for i=2:nargin %在MATLAB中它也是一个固定变量，表示在函数调用时实际输入变量的个数。 eval(a=conv(a,a i

4、nt2str(i) );) %eval( )函数执行MATLAB语句构成的字符串。end这里定义的函数，可以一次执行十个多项式的乘积。如上例可以写成：convs(1,0,1,1,1,1,1,1,1,3,2,5)注意：写完函数后一定要按函数名来保存该文件。相应地，对离散时间系统的动态模型，用脉冲传递函数描述，其输入方法与此类似，不多举例。如果是MIMO系统，则用传递矩阵描述，如： 则可表示为：num= 1,1;1den=1 2 2 ;1 1 二控制系统的零极点模型：（zpk模型）是传递函数的另一种表现形式，对SISO系统：G(s)= 将零点、极点及K值输入即可建立零极点模型。 z=-z,-z ,

5、-z p=-p,-p ,-p k=k对于给出的传函来说，分别对分子分母作因式分解，则可以得出系统的零极点模型。这可以通过求出分子，分母多项式的根来实现。MATLAB中提供了多项式求根的函数，roots ( )调用格式： z=roots(a) 其中：z 各个根所构成的向量 a 多项式系数向量系统增益k即为原传递函数分子的最高项系数与分母最高项系数的比值。分别求出分子分母多项式的根，即可得到系统的零，极点模型。对于多输入多输出系统，应分别对每个输入求出系统的零极点模型最后才可以获得整个系统的零极点模型（为矩阵形式）。两种模型之间的转换函数：z ,p ,k=tf2zp(num , den);num

6、, den =zp2tf(z ,p ,k)三. 控制系统状态方程模型：( ss 模型)LTI系统的状态方程： = Ax + Bu Y = Cx + Du只要将A，B，C，D几个矩阵输入进去即可。对于离散系统来说，也与上面类似。MATLAB还提供了由系统状态方程转换为系统零极点模型及传递函数模型的函数。num , den=ss2tf(A ,B ,C ,D , iu ) %iu表示输入的序号（对多输入系统）z, p ,k=ss2zp(A ,B ,C ,D , iu ) %z ,p ,k 表示对第iu个输入信号的传递函数的零极点。A , B ,C ,D= tf2ss (num , den)A , B

7、 ,C ,D=zp2ss(z ,p ,k)四 在MATLAB控制系统工具箱中的LTI对象：为了避免对一个系统采用多个分离变量进行描述，新版本的控制系统工具箱，将LTI系统的各种描述封装成一个对象，即用一个变量来描述。 在控制系统工具箱中，有以上讲述的三种对象，即ss 对象，tf对象和zpk对象。每种系统模型的生成和模型间的转换均可以通过一个函数来实现。下面介绍这些函数：1 dss( )函数：生成系统的状态空间模型。格式： sys= dss (a ,b ,c ,d ,e ) 该调用生成连续系统的状态空间模型：= = E为非奇异阵（E阵奇异，为奇异系统）。sys = (a ,b , c, d ,e

8、 ,Ts)生成离散系统的状态空间模型。= Ts为采样周期2filt()函数：生成DSP（数字信号处理）形式的离散传递函数： 格式：sys=filt(num ,den) - 生成离散传递函数模型 sys=filt(num ,den ，Ts)-定义模型的采样周期注：DSP（数字信号处理）形式的离散传递函数如下面的形式：3ss()函数：生成状态空间模型，或者将传递函数及零极点模型转换成状态空间模型。 格式：1）sys=ss(a,b,c,d)-生成连续系统的状态空间模型。形式为： y =Cx + Du 2）sys=ss(a,b,c,d,Ts)-生成离散状态空间模型x(n+1)=Ax(n)+Bu(n)y

9、(n)=Cx(n)+Du(n) 3）sys_ss=ss(sys)-将任意的LTI对象sys转换成状态空间模型。4tf()函数：生成传递函数模型，或将零极点模型及状态空间模型转换成传递函数模型。 格式为：sys=tf(num,den)-生成连续时间系统传函模型: G(s)= sys=tf(num ,den, Ts)-生成离散时间系统传函。tfsys=tf(sys)-将任意的LTI对象转换成传递函数模型。5zpk()函数：生成零极点模型或者将其他模型转化成零极点模型。 格式：sys=zpk(z ,p ,k) - 连续系统的零极点增益模型。sys=zpk(z ,p ,k,Ts)- 离散时间系统的零极

10、点增益模型。zsys=zpk(sys) - 将任意LTI对象转换成零极点增益模型。五. 控制系统的结构图描述及转换： 由于复杂的结构图建模可以用SIMUKINK。在这里我们仅介绍简单的结构图描述，也可称为系统建模。 典型连接：1串联连接： =子系统1： =Ax+B u 子系统2： =Ax+Bu y= cx+Du y= Cx+Du串联后系统模型为： =+ y= +D2D1u1 2并联连接： +并联后系统模型： =+u y=+(+)u3负反馈连接： - 为了简化,设，则反馈后：=+=+ (r-) =+r- (+) =+r- = 所以，状态方程为： 在MATLAB工具箱中，提供了子系统的连接处理函数

11、：1） series()函数：系统串联实现格式：sys=series(,)A,B,C,D=series(A,B,C,D, A,B,C,D)2） parallel()函数：系统并联实现。格式：sys=parallel(,)3） feedback()函数：系统反馈连接格式：sys=feedback(,) sys=feedback(,sign) sign定义反馈形式：正反馈，sign=+1；负反馈，sign=-1。 4） append函数：多个LTI系统的组合格式:sys=append(,)对N个子系统组合： 设各个子系统的模型为： (s), (s), (s)则合成后子系统为分块对角阵： 如子系统为

12、状态空间模型，则组合后： 第二节 控制系统的稳定性分析在计算机仿真中，控制系统稳定性的分析方法有：1 求闭环特征方程的根；2 化为零极点模型，看极点是否在s右半平面；3 对状态空间形式（闭环），求A阵的特征值eig(A)；4 用Lyapunov方程求解。 在自控原理中，我们均采用间接方法劳斯判据，奈氏判据等，但由于在MATLAB中很容易地求解多项式方程。因此，我们可以直接求出特征方程的根，可用多项式求根函数roots()。 例1：判断如下系统的稳定性： -可编程如下：numg=1; deng=1 1 2 23; numf=1; denf=1; (num,den)= feedback(numg,

13、deng,numf,denf,-1); roots(den)例2：判断状态空间形式描述的系统的稳定性，两个输入。A=0 1 1；-6 11 6；-6 11 5；B=0 0 1；C=1 0 0； D=0；,=ss2zp(A,B,C,D,1),=ss2zp(A,B,C,D,2)ii=find(real()0), =length(ii);%find（ ）：找出数组中满足实部值大于0的所有元素的下标，并将结果返回到数组中。ii1= find(real()0), =length(ii1);if (+0),disp(system is Unstable)else disp(system is Stable

14、);end一般当系统中含有非线性环节时，即不能用上述方法，一般需求解李雅普诺夫方程： P+PA=-Q 求出P为正定对称矩阵，则系统是稳定的。李雅普诺夫（Lyapunov）方程求解函数Lyap()（采用Schur方法编写的），其调用格式：x=lyap(A,B,C)该函数可求解：AX+XB=-C方程如要求解： P+PA=-Q 则可写成：p=lyap(A,Q)注：判断稳定性的方法还可确定是否为最小相位系统（没有右半平面零极点的系统）。第三节 控制系统的根轨迹在控制系统工具箱中提供了系统根轨迹绘制及分析函数，还有一个GUI分析工具。下面介绍一下这些根轨迹的绘制及分析函数：1 rlocfind():计算给定根的根轨迹增益 格式：k,poles=rlocfind(sys) k,poles= rlocfind(sys,p)1 在LTI对象的根轨迹图中显示出十字光标，当用户选择其中一点时，其相应的增益由k记录，与增益相关的所有极点记录poles中。若要使用该函数，必须首先在当前窗口上绘制系统的根轨迹。2 定义要得到增益的根矢量P，即事先给出极点。除了显示出该根对应的增益以外，还显示出