东北大学 自动化 控制系统计算机辅助设计实验.pdfVIP

RTO

RTO

RTO

11

控制系统计算机辅助设计

第一部分

functiondx=lorenzeq(t,x)

dx=[-x(2)-x(3);

x(l)+0.2*x(2);

0.2+(x(l)-5.7)*x(3);]

»x0=[0;0;0];

»|t,y|=odc45(lorcnzcq,[0,1001,x0);

»plot(t,y)

»figure;plot(y(:,I),y(:,2),y(:,)),grid,

5

-100

-5

•20-10

2

function[y,yeq]=f2(x)

yeq=[];

AA

y=4*x(l)2+x(2)2-4;

»Aeq=[];Beq=[];A=[];B=[];

»xm=[0;0];xM=[];xO=[0;0];

»fl=inline(x(1)A2-2*x(l)+x(2));

»[x,f]=fmincon(fl,xO,A,B,A,eq,Beq,xin,xM,f2a);x,f

ans=

1.00000

f=

-I

3

(a)

AAAAA

s=tf(s);G=(s3+4*s+2)/s3/(s2+2)/[(s2+l)3+2*s+5]

Transferfunction:

sA+4s+2

AAAAAAA

sll+5s9+9s7+2s6+12s5+4s4+!2s3

(b)

»z=tf(z,0.1);

»H=(zA2+0.568)/(z-1)/(zA2-0.2*z+0.99)

Transferfunction:

zA2+0.568

AA

z-1.2z2+1.19z-0.99

Samplingtime:0.1

4

»A=fO10;00l;-5-4-13];

»B=[0;0;2];

»C=n00;000;000);

»D=[0];

»G=ss(A,B,C,D);

»G

a=

xlx2x3

xl010

x2001

x3-5-4-13

b=

ul

xl0

x20

x32

xlx2x3

yi।00

y2000

y30