大学Matlab课程第3讲控制系统的数学描述与建模.ppt

第二章控制系统的数学描述与建模,引言,控制系统研究中的三大问题建模分析设计,CH2、控制系统的数学描述与建模,研究控制系统的数学模型的重要性：要对系统进行仿真处理，首先应当知道系统的数学模型，然后才可以对系统进行模拟。知道了系统的模型，才可以在此基础上设计一个合适的控制器，使得系统响应达到预期的效果，从而符合工程实际的需要。,在线性系统理论中，一般常用的数学模型形式有：传递函数模型（系统的外部模型）、状态方程模型（系统的内部模型）、零极点增益模型、部分分式模型等。这些模型之间都有着内在的联系，可以相互进行转换。,2.1常见模型建立,微分方程是控制系统模型的基础；一般来讲，利用机械学、电学、力学等物理规律，便可以得到控制系统的动态方程，这些方程就是系统的微分方程。控制的分析问题就是在给定控制输入时，求解微分方程的解并分析其动态特性；通过拉氏变换和反变换，可以得到线性定常系统的解析解，解析解是精确的，然而通常寻找解析解是困难的。MATLAB提供了ode45等微分方程的数值解法函数。控制系统的综合问题就是为使微分方程的解达到期望的轨迹，设计微分方程的控制输入。,2.1.1微分方程模型,例exp3_1.m,电路图如下，R=1.4欧，L=2亨，C=0.32法，初始状态：电感电流为零，电容电压为0.5V，t=0时刻接入1V的电压，求0t15s时，i(t)，vo(t)的值，并且画出电流与电容电压的关系曲线。,传递函数是经典控制论描述系统的数学模型之一，它表达了输入量和输出量之间的关系。在Matlab中，可以利用分别定义的传递函数分子、分母多项式对其加以描述。num=c1,c2,cn-1,cnden=1,a1,a2,an-1,an注意：它们都是按s的降幂进行排列的。,一、连续系统的传递函数模型,2.1.2传递函数模型,一、连续系统的传递函数模型举例：,2.1.2传递函数模型,已知系统传递函数为：,试在Matlab中将上述传递函数模型表示出来。,num=2,9;den=13246;sys1=tf(num,den),一、连续系统的传递函数模型举例：,2.1.2传递函数模型,已知系统传递函数为：,试在Matlab中将上述传递函数模型表示出来。,num=7*23;den=conv(conv(conv(100,31),conv(12,12),5038);sys1=tf(num,den),在Matlab中，对于离散系统，同样可以建立相应的系统模型。num=cm,cm-1,c1,c0den=an,an-1,a1,a0sys=tf(num,den,T)其中，T为系统采样周期。,离散系统的脉冲传递函数模型,2.1.2传递函数模型,零极点模型实际上是传递函数模型的另一种表现形式，其原理是分别对原系统传递函数的分子、分母进行分解因式处理，以获得系统的零点和极点的表示形式。,二、零极点增益模型,式中：K为系统增益，zi为零点，pj为极点,在MATLAB中，零极点增益模型用z,p,K矢量组表示。即：z=z1,z2,zmp=p1,p2,.,pnK=ksys=zpk(z,p,K),零极点增益模型的Matlab表示,举例：num=1,11,30,0;den=1,9,45,87,50;z,p,k=tf2zp(num,den),结果表达式：,控制系统常用到并联系统，这时就要对系统函数进行分解，使其表现为一些基本控制单元的和的形式，也就是部分分式表示，如下式所示：,三、部分分式模型,MATLAB提供了函数r,p,k=residue(b,a)，它的功能是对两个多项式的比进行部分展开，以及把传递函数分解为微分单元的形式。其中向量b和a是按s的降幂排列的多项式系数。部分分式展开后，余数返回到向量r，极点返回到列向量p，常数项返回到k.,部分分式模型的MATLAB表示：,举例：num=1,11,30,0;den=1,9,45,87,50;r,p,k=residue(num,den),结果表达式为：,经典控制理论用传递函数描述输入输出关系现代控制理论则用状态空间描述法（状态方程和输出方程）来表达输入输出关系，它揭示了系统内部状态对系统性能的影响。,2.1.3状态空间描述,在MATLAB中，系统状态空间用（A,B,C,D)矩阵组表示。表达式如下：sys=ss(A,B,C,D)对于离散系统，同样可建立系统的状态空间模型：dsys=ss(A,B,C,D,T),T为采样周期,2.1.3状态空间描述,举例：,系统为两输入、两输出系统A=16910;31268;47911;5121314;B=46;24;22;10;C=0021;8022;D=zeros(2,2);sys=ss(A,B,C,D),2.2模型的相互转换,描述系统的主要模型有：传递函数模型零极点模型状态空间模型等在进行系统分析时，往往需要根据不同的要求选择不同形式的数学模型。,2.2系统模型的相互转换,形式不同，实质内容等价,ss2tf：状态空间模型转换为传递函数模型ss2zp：状态空间模型转换为零极点增益模型tf2ss：传递函数模型转换为状态空间模型tf2zp：传递函数模型转换为零极点增益模型zp2ss：零极点增益模型转换为状态空间模型zp2tf：零极点增益模型转换为传递函数模型,模型转换函数,ss2tf,State-spacetotransferfunctionconversion.NUM,DEN=SS2TF(A,B,C,D,iu)calculatesthetransferfunction:ofthesystem:fromtheiuthinput.,iu用来指定第n个输入，当只有一个输入时可忽略。,用法举例：,例：已知连续系统的系数矩阵如下：,求取该系统相应的传递函数模型,A=200;041;004;B=1;0;1;C=110;D=0;num,den=ss2tf(A,B,C,D)G=tf(num,den),A=200;041;004;B=1;0;1;C=110;D=0;sys=ss(A,B,C,D)G2=tf(sys),例exp3_5.m,已知连续系统的系数矩阵为：求系统相应的传递函数。,ss2zp,tf2zp,z,p,k=ss2zp(A,B,C,D,iu)z,p,k=tf2zp(num,den)Gzp=zpk(sys),将非零极点形式的模型转换成零极点模型,用法举例：已知系统状态空间模型为：求其零极点模型,A=01;-1-2;B=0;1;C=1,3;D=1;iu=1;z,p,k=ss2zp(A,B,C,D,iu)sys=zpk(z,p,k）,tf2ss,zp2ss,A,B,C,D=tf2ss(num,den)A,B,C,D=zp2ss(z,p,k)Syss=ss(sys),将非状态空间形式的模型转换成状态空间模型,例：系统的零极点增益模型如下：求系统的传递函数模型和状态空间模型,z=-3;p=-1,-2,-5;k=6;num,den=zp2tf(z,p,k)a,b,c,d=zp2ss(z,p,k)注意：零极点的输入可以写成行向量，也可以写成列向量。,已知一个单输入三输出系统的传递函数模型为：,num=00-2;0-1-5;120;den=16116;A,B,C,D=tf2ss(num,den);sys=ss(A,B,C,D),例：,连续和离散系统之间的转换,sysd=tf(num,den,T)sysd=ss(a,b,c,d,T)sysd=zpk(z,p,k,T)sysd=c2d(sysc,T,method),T:采样周期method:表示指定转换方式zoh：零阶保持器foh：三角形近似tustin：双线性变换,连续和离散系统之间的转换,例exp3_6.m已知系统传递函数为采样周期T=0.2秒，试将其进行离散化处理,2.3状态方程的变换,2.3状态方程的变换及实现,传递函数,状态空间模型,唯一,不唯一,？,2.3.1状态方程的相似变换,S2=ss2ss(S1,T)S2：变换后的系统状态方程模型S1：原系统的状态方程模型T：非奇异变换阵,2.3.2规范型状态方程的实现,基本格式：G1=canon(sys,type)G1,T=canon(sys,type)sys表示原系统状态方程模型type确定规范形式的类型modal:约当标准型companion:伴随矩阵形式G1:规范后的状态方程模型T:状态变换阵,必须为状态空间模型,例：,已知系统的系数阵为:对其进行规范型变换(约当变换),并给出变换阵.,A=5210;0460;0-3-50;0-3-6-1;B=1;2;3;4;C=1252;D=0;sys=ss(A,B,C,D);G,T=canon(sys,modal),2.3.3系统的均衡实现,基本格式：Ab,Bb,Cb,G,T=balreal(A,B,C)T:均衡变换阵G:均衡系统的Gram阵满足下列变换关系功能：实现系统的均衡变换，克服计算误差,例：,已知系统的系数阵为:对其进行均衡变换,并给出变换阵.,A=-100;0-25;B=10(-5);105;C=105,10(-5);D=0;Ab,Bb,Cb,G,T=balreal(A,B,C),2.3.4系统的降阶实现,基本格式：RSYS=modred(sys,ELIM)RSYS=modred(sys,ELIM,mdc)RSYS=modred(sys,ELIM,del)ELIM：待消去的状态mdc：表示在降阶中保证增益的匹配del：表示降阶中不能保证增益的匹配,例：,已知系统的系数矩阵为：在尽可能保持系统基本特征的情况下进行降阶处理。,A=-310-1;-0.5-11-1;-1.51-20;-1.521-4;B=1;0;0;0;C=10-10;Ab,Bb,Cb,G,T=balreal(A,B,C)sys=ss(Ab,Bb,Cb,0)sysr=modred(sys,3:4,mdc),例：,2.4常用模型的建立与模型连接,2.4.1基本模型系统的建立,二阶系统的生成：A,B,C,D=ord2(Wn,z)num,den=ord2(Wn,z)Wn：自然角频率，z：阻尼因子,2.4.1基本模型系统的建立,随机n阶系统的建立：num,den=rmodel(n):随机生成n阶稳定传递函数模型num,den=rmodel(n,p):随机生成单入p出的n阶稳定传递函数模型A,B,C,D=rmodel(n):随机生成n阶稳定SISO状态方程模型A,B,C,D=rmodel(n,p,m):随机生成n阶稳定p出m入状态方程模型drmodel():生成相应的离散模型,程序如下:wn=2.2;z=0.8;A,B,C,D=ord2(wn,z);sys1=ss(A,B,C,D)num,den=ord2(wn,z);sys2=tf(num,den)n=4;p=4;m=2;a,b,c,d=rmodel(n,p,m);sys3=ss(a,b,c,d)num1,den1=drmodel(n,2);sys4=tf(num1(1,:);num1(2,:),den1),例:基本模型的生成,2.4.1基本模型系统的建立,系统模型的重构：子系统的选取与删除Ae,Be,Ce,De=ssselect(A,B,C,D,inputs,outputs)Ae,Be,Ce,De=ssselect(A,B,C,D,inputs,outputs,states)Ar,Br,Cr,Dr=ssdelete(A,B,C,D,inputs,outputs)Ar,Br,Cr,Dr=ssdelete(A,B,C,D,inputs,outputs,states),例：,现利用ssselect函数在原系统的基础上构造新系统，保留1、3输入信号，1、2输出以及1、2、4状态变量。,状态的增广在对系统进行分析研究时，往往需要对状态在系统（输出）中加以增广。如对系统进行全状态反馈研究时，考虑到u=Kx，因此一般就需要在输出方程中增广状态。MATLAB中提供了一种状态增广函数augstate,即：asys=augstate(sys)构成新系统为：,例：将所有的状态增广到系统输出中。,2.4.2系统组合与连接,系统组合的方式:串联:series并联:parallel反馈:feedback,格式：sys=series(sys1,sys2);sys=sys1*sys2对于多输入多输出系统,串联表达式为:sys=series(sys1,sys2,outputs1,inputs2)上述函数实现由outputs1指定的sys1输出端连接到由inputs2指定的sys2输入端.,1、模型串联：series,格式：sys=parallel(sys1,sys2)对于多输入多输出系统,串联表达式为:sys=parallel(sys1,sys2,in1,in2,out1,out2)in1,in2指定了相连接的输入端.out1,out2指定了相连接的输出端,2、模型并联：parallel,格式：sys=feedback(sys1,sys2,sign)sign缺省时为负反馈,sign=1时为正反馈对于多输入多输出系统,表达式为:sys=feedback(sys1,sys2,feedin,feedout,sign)Feedin:sys1的输入向量,指定哪些sys1的输入与反馈环相连Feedout:sys1的输出向量,指定sys1的哪些输出端用于反馈,3、反馈连接：feedback,举例应用：exp3_2.m1）系统1为：系