实验八matlab状态空间分析.pdfVIP

-

实验八线性系统的状态空间分析

§8.1用MATLAB分析状态空间模型

1、状态空间模型的输入

线性定常系统状态空间模型

将各系数矩阵按常规矩阵形式描述。

在MATLAB里，用函数SS()来建立状态空间模型

例8.1*系统微分方程

求该系统的状态空间模型。

解：将上述微分方程写成状态空间形式

010

A，B

731

C50，D0

调用MATLAB函数SS()，执行如下程序

%MATLABPrograme*ample6.1.m

A=[01;-7-3];

B=[0;1];

C=[50];

D=0;

sys=ss(A,B,C,D)

运行后得到如下结果

a=

*1*2

*101

*2-7-3

b=

u1

*10

*21

c=

*1*2

y150

d=

u1

y10

Continuous-timemodel.

2、状态空间模型与传递函数模型转换

z.

-

状态空间模型用sys表示，传递函数模型用G表示。

G=tf(sys)

sys=ss(G)

状态空间表达式向传递函数形式的转换

G=tf(sys)

Or[num,den]=ss2tf(A,B,C,D)多项式模型参数

[num,den]=ss2tf(A,B,C,D,iu)

[z,p,k]=ss2zp(A,B,C,D,iu)零、极点模型参数

iu用于指定变换所需的输入量，iu默认为单输入情况。

传递函数向状态空间表达式形式的转换

sys=ss(G)

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

[A,B,C,D]=zp2ss(z,p,k)

例8.2

试用矩阵组[a，b，c，d]表示系统，并求出传递函数。

%MATLABPrograme*ample6.2.m

a=[-0.560.05;-0.250];

b=[0.031.14;0.110];

c=[10;01];

d=zeros(2,2);

sys=ss(a,b,c,d)

G1=tf(sys)

G2=zpk(sys)

运行后得到如下结果

a=

*1*2

*1-0.560.05

*2-0.250

b=

u1u2

*10.031.14

*20.110

c=

*1*2

y110

y201

d=

z.

-

u1u2

y100

y200

Continuous-timemodel.

Transf