B =
0 1 0 1 C =
1 3 1 D =
1
单位负反馈连接 A =
0 1 0 -1 -2 -1 0 0 0 1 3 -2 B =
0 1 0 1 C =
1 3 -1 D =
1
4 0 -4 1 -7 -4 单位正反馈连接 A =
0 1 0 0 -1 -2 1 4 0 0 0 1 1 3 0 1 B =
0 1 0 1 C =
1 3 1 4 D =
1
源程序:
A=[0,-2;1,-3]; t=.2;
F=expm(A*t) %转移矩阵 B=[2;0]; C=[0,3]; D=[0]; x0=[1,1]; t=[0,.2]; u=0*t;
[y,x]=lsim(A,B,C,D,u,t,x0)
输出:
F =
0.9671 -0.2968 0.1484 0.5219 y =
3.0000 2.0110 x =
1.0000 1.0000
0.6703 0.6703
结论:
t=0.2时,系统响应为x1(0)?x2(0)?0.6703,y(0.2)=2.0110
源程序:
A=[-3,1;1,-3]; B=[1,1;1,1]; C=[1,1;1,-1]; D=[0];
Qc=ctrb(A,B) Qo=obsv(A,C) Rc=rank(Qc) Ro=rank(Qo)
输出:
Qc =
1 1 -2 -2 1 1 -2 -2
Qo =
1 1 1 -1 -2 -2 -4 4 Rc =
1 Ro =
2
结论:
能控性矩阵和能观性矩阵的秩分别为1,2,又系统阶次是2,故系统是不可控的,是可观测的。
源程序:
a=[0 1;-2 -3]; b=[0;1];
p=[b a*b a^2*b] c=[2 0];
q=[c;c*a;c*a*a] rank(q)
a=[0 1;-2 -3];
b=[0;1]; c=[2 0]; pe=[-3;-3];
G=(acker(a',c',pe))'
输出:
p =
0 1 -3 1 -3 7 q =
2 0 0 2 -4 -6
ans =
2 G =
1.5000 -1.0000
结论:
反馈增益矩阵为[1.5 -1],观测器的状态方程为
.?^??^??x1???31??x1??0??1.5??.?? ??????u???y??^??0?3??^??1???1???x2???x2?T