状态方程与输出方程的组合称为状态空间表达式，又称为动态方程，经典控制理论用传递函数将输入—输出关系表达出来，而现代控制理论则用状态方程和输出方程来表达输入—输出关系，揭示了系统内部状态对系统性能的影响。 自动控制网www.eadianqi.com版权所有

    在MATLAB中，系统状态空间用（A,B,C,D)矩阵组表示。 自动控制网www.eadianqi.com版权所有

    

    举例：

本文来自www.eadianqi.com

    系统为一个两输入两输出系统 本文来自www.eadianqi.com

    》A=[16910;31268;47911;5121314]; 本文来自www.eadianqi.com

    》B=[46;24;22;10]; 本文来自www.eadianqi.com

    》C=[0021;8022]; 自动控制网www.eadianqi.com版权所有

    》D=zeros(2,2); 本文来自www.eadianqi.com

 模型的转换与连接 自动控制网www.eadianqi.com版权所有

    在一些场合下需要用到某种模型，而在另外一些场合下可能需要另外的模型，这就需要进行模型的转换。

    模型转换的函数包括： 自动控制网www.eadianqi.com版权所有

    residue：传递函数模型与部分分式模型互换

自动控制网www.eadianqi.com版权所有

    ss2tf：状态空间模型转换为传递函数模型 自动控制网www.eadianqi.com版权所有

    ss2zp：状态空间模型转换为零极点增益模型

    tf2ss：传递函数模型转换为状态空间模型

    tf2zp：传递函数模型转换为零极点增益模型 自动控制网www.eadianqi.com版权所有

    zp2ss：零极点增益模型转换为状态空间模型 本文来自www.eadianqi.com

    zp2tf：零极点增益模型转换为传递函数模型 本文来自www.eadianqi.com

一、模型的转换

    用法举例： 自动控制网www.eadianqi.com版权所有

    1）已知系统状态空间模型为： 本文来自www.eadianqi.com

     自动控制网www.eadianqi.com版权所有

    》A=[01;-1-2];B=[0;1];

自动控制网www.eadianqi.com版权所有

    》C=[1,3];D=[1]; 自动控制网www.eadianqi.com版权所有

    》[num,den]=ss2tf(A,B,C,D,iu) 本文来自www.eadianqi.com

    ％iu用来指定第n个输入，当只有一个输入时可忽略。 自动控制网www.eadianqi.com版权所有

    》num=152;den=121;

自动控制网www.eadianqi.com版权所有

    》[z,p,k]=ss2zp(A,B,C,D,iu) 自动控制网www.eadianqi.com版权所有

    》z=-4.5616p=-1k=1 自动控制网www.eadianqi.com版权所有

    -0.4384-1

本文来自www.eadianqi.com

     自动控制网www.eadianqi.com版权所有

    2）已知一个单输入三输出系统的传递函数模型为：

    

    》num=[00-2;0-1-5;120];den=[16116];

    》[A,B,C,D]=tf2ss(num,den)

    》A=-6-11-6B=1C=00-2D=0

自动控制网www.eadianqi.com版权所有

    10000-1-50 自动控制网www.eadianqi.com版权所有

    01001200 自动控制网www.eadianqi.com版权所有

    

本文来自www.eadianqi.com

    3）系统的零极点增益模型： 自动控制网www.eadianqi.com版权所有

    

本文来自www.eadianqi.com

    》z=[-3];p=[-1,-2,-5];k=6;

    》[num,den]=zp2tf(z,p,k) 本文来自www.eadianqi.com

    》num=00618den=181710 自动控制网www.eadianqi.com版权所有

     自动控制网www.eadianqi.com版权所有

    》[a,b,c,d]=zp2ss(z,p,k) 自动控制网www.eadianqi.com版权所有

    》a=-1．000000b=1

    2．0000-7．0000-3.16231 自动控制网www.eadianqi.com版权所有

    03.162300 自动控制网www.eadianqi.com版权所有

    c=001.8974d=0

自动控制网www.eadianqi.com版权所有

    注意：零极点的输入可以写出行向量，也可以写出列向量。 本文来自www.eadianqi.com

    4）已知部分分式：

    

自动控制网www.eadianqi.com版权所有

    》r=[-0.25i,0.25i,-2]; 本文来自www.eadianqi.com

    》p=[2i,-2i,-1];k=2; 本文来自www.eadianqi.com

    》[num,den]=residue(r,p,k)

自动控制网www.eadianqi.com版权所有

    》num=

    2091

本文来自www.eadianqi.com

    》den=

本文来自www.eadianqi.com

    1144

    注意余式一定要与极点相对应。 自动控制网www.eadianqi.com版权所有

二、模型的连接 本文来自www.eadianqi.com

    （一）并联：parallel

本文来自www.eadianqi.com

    格式：

本文来自www.eadianqi.com

    [a,b,c,d]=parallel(a1,b1,c1,d1,a2,b2,c2,d2) 本文来自www.eadianqi.com

    ％并联连接两个状态空间系统。

    [a,b,c,d]=parallel(a1,b1,c1,d1,a2,b2,c2,d2,inp1,inp2,out1,out2)

自动控制网www.eadianqi.com版权所有

    ％inp1和inp2分别指定两系统中要连接在一起的输入端编号，从u1,u2,…,un依次编号为1,2,…,n；out1和out2分别指定要作相加的输出端编号，编号方式与输入类似。inp1和inp2既可以是标量也可以是向量。out1和out2用法与之相同。如inp1=1,inp2=3表示系统1的第一个输入端与系统2的第三个输入端相连接。 自动控制网www.eadianqi.com版权所有

    若inp1=[13],inp2=[21]则表示系统1的第一个输入与系统2的第二个输入连接，以及系统1的第三个输入与系统2的第一个输入连接。

    [num,den]=parallel(num1,den1,num2,den2) 自动控制网www.eadianqi.com版权所有

    ％将并联连接的传递函数进行相加。

    （二）串联：series 本文来自www.eadianqi.com

    格式：

本文来自www.eadianqi.com

    [a,b,c,d]=series(a1,b1,c1,d1,a2,b2,c2,d2) 自动控制网www.eadianqi.com版权所有

    ％串联连接两个状态空间系统。

自动控制网www.eadianqi.com版权所有

    [a,b,c,d]=series(a1,b1,c1,d1,a2,b2,c2,d2,out1,in2) 本文来自www.eadianqi.com

    ％out1和in2分别指定系统1的部分输出和系统2的部分输入进行连接。 本文来自www.eadianqi.com

    [num,den]=series(num1,den1,num2,den2)

自动控制网www.eadianqi.com版权所有

    ％将串联连接的传递函数进行相乘。 本文来自www.eadianqi.com

    （三）反馈：feedback

自动控制网www.eadianqi.com版权所有

    格式： 本文来自www.eadianqi.com

    [a,b,c,d]=feedback(a1,b1,c1,d1,a2,b2,c2,d2)

本文来自www.eadianqi.com

    ％将两个系统按反馈方式连接，一般而言，系统1为对象，系统2为反馈控制器。 本文来自www.eadianqi.com

    [a,b,c,d]=feedback(a1,b1,c1,d1,a2,b2,c2,d2,sign)

本文来自www.eadianqi.com

    ％系统1的所有输出连接到系统2的输入，系统2的所有输出连接到系统1的输入，sign用来指示系统2输出到系统1输入的连接符号，sign缺省时，默认为负，即sign=-1。总系统的输入/输出数等同于系统1。

    [a,b,c,d]=feedback(a1,b1,c1,d1,a2,b2,c2,d2,inp1,out1) 自动控制网www.eadianqi.com版权所有

    ％部分反馈连接，将系统1的指定输出out1连接到系统2的输入，系统2的输出连接到系统1的指定输入inp1，以此构成闭环系统。 自动控制网www.eadianqi.com版权所有

    [num,den]=feedback(num1,den1,num2,den2,sign) 自动控制网www.eadianqi.com版权所有

    ％可以得到类似的连接，只是子系统和闭环系统均以传递函数的形式表示。sign的含义与前述相同。

本文来自www.eadianqi.com

    （四）闭环：cloop（单位反馈） 自动控制网www.eadianqi.com版权所有

    格式： 自动控制网www.eadianqi.com版权所有

    [ac,bc,cc,dc]=cloop(a,b,c,d,sign)

自动控制网www.eadianqi.com版权所有

    ％通过将所有的输出反馈到输入，从而产生闭环系统的状态空间模型。当sign=1时采用正反馈；当sign=-1时采用负反馈；sign缺省时，默认为负反馈。 自动控制网www.eadianqi.com版权所有

    [ac,bc,cc,dc]=cloop(a,b,c,d,outputs,inputs) 本文来自www.eadianqi.com

    ％表示将指定的输出outputs反馈到指定的输入inputs，以此构成闭环系统的状态空间模型。一般为正反馈，形成负反馈时应在inputs中采用负值。

    [numc,denc]=cloop(num,den,sign)

    ％表示由传递函数表示的开环系统构成闭环系统，sign意义与上述相同。

自动控制网www.eadianqi.com版权所有

    举例应用：

自动控制网www.eadianqi.com版权所有

    1）系统1为： 自动控制网www.eadianqi.com版权所有

    系统2为： 自动控制网www.eadianqi.com版权所有

    求按串联并联正反馈负反馈连接时的系统状态方程及系统1按单位负反馈连接时的状态方程。 自动控制网www.eadianqi.com版权所有

    exp3_2．m 自动控制网www.eadianqi.com版权所有

    clc

    clear 自动控制网www.eadianqi.com版权所有

    moreon 自动控制网www.eadianqi.com版权所有

    a1=[01;-1-2]; 本文来自www.eadianqi.com

    b1=[0;1];

    c1=[13];d1=[1]; 本文来自www.eadianqi.com

    a2=[01;-1-3];

    b2=[0;1];

本文来自www.eadianqi.com

    c2=[14];d2=[0];

自动控制网www.eadianqi.com版权所有

    %串联连接

    disp('串联连接')

    [a,b,c,d]=series(a1,b1,c1,d1,a2,b2,c2,d2) 自动控制网www.eadianqi.com版权所有

    %并联连接 本文来自www.eadianqi.com

    disp('并联连接') 自动控制网www.eadianqi.com版权所有

    [a,b,c,d]=parallel(a1,b1,c1,d1,a2,b2,c2,d2) 自动控制网www.eadianqi.com版权所有

    %正反馈

本文来自www.eadianqi.com

    disp('正反馈连接')

本文来自www.eadianqi.com

    [a,b,c,d]=feedback(a1,b1,c1,d1,a2,b2,c2,d2,+1) 自动控制网www.eadianqi.com版权所有

    %负反馈 本文来自www.eadianqi.com

    disp('负反馈连接') 本文来自www.eadianqi.com

    [a,b,c,d]=feedback(a1,b1,c1,d1,a2,b2,c2,d2)

自动控制网www.eadianqi.com版权所有

    %单位负反馈

    disp('单位负反馈连接')

    [a,b,c,d]=cloop(a1,b1,c1,d1) 本文来自www.eadianqi.com

三、模型的属性 自动控制网www.eadianqi.com版权所有

    ctrb和obsv函数可以求出状态空间系统的可控性和可观性矩阵。 本文来自www.eadianqi.com

    格式：co=ctrb(a,b)ob=obsv(a,c) 本文来自www.eadianqi.com

    对于n×n矩阵a，n×m矩阵b和p×n矩阵c

自动控制网www.eadianqi.com版权所有

    ctrb(a,b)可以得到n×nm的可控性矩阵 自动控制网www.eadianqi.com版权所有

    co=[baba2b…an-1b]

本文来自www.eadianqi.com

    obsv(a,c)可以得到nm×n的可观性矩阵

    ob=[ccaca2…can-1]’

    当co的秩为n时，系统可控；当ob的秩为n时，系统可观。

    exp3_3．m 自动控制网www.eadianqi.com版权所有

    线性系统H(S)=(s+alph)/(s^3+10s^2+27s+18)，当alph分别取－1，0，1时，判别系统的可控性和可观性，并求出相应的状态方程。

自动控制网www.eadianqi.com版权所有

    clc

本文来自www.eadianqi.com

    clear 本文来自www.eadianqi.com

    moreon

    foralph=[-1:1] 本文来自www.eadianqi.com

    alph 本文来自www.eadianqi.com

    num=[1,alph]; 本文来自www.eadianqi.com

    den=[1102718]; 自动控制网www.eadianqi.com版权所有

    [a,b,c,d]=tf2ss(num,den) 本文来自www.eadianqi.com

    cam=ctrb(a,b) 自动控制网www.eadianqi.com版权所有

    rcam=rank(cam) 自动控制网www.eadianqi.com版权所有

    oam=obsv(a,c) 自动控制网www.eadianqi.com版权所有

    roam=rank(oam) 自动控制网www.eadianqi.com版权所有

    end 本文来自www.eadianqi.com

    moreoff 本文来自www.eadianqi.com