.
状态反馈设计的转置,可以用H=(acker(A’,C’,V’))’来确定一般系统的观测器矩阵,用命令eig(estim(sysold,H))来检验极点配置。
三、实验步骤及结果
step(A,B,C,D); num=[0 0 1]; den=[1 3 2];
[A,B,C,D]=tf2ss(num,den) A =-3 -2 1 0 B =1 0 C = 0 1 D = 0
2、配置后系统的时间响应曲线为:
Word文档
.
A=[-3 -2;1 0]; B=[1;0]; C=[0 1]; D=0;
P=[-1+sqrt(-1);-1-sqrt(-1)]; K=acker(A,B,P) K = -1 0
>> disp('极点配置后的闭环系统为') 极点配置后的闭环系统为 >> sysnew=ss(A-B*K,B,C,D) a = x1 x2 x1 -2 -2 x2 1 0 b = u1 x1 1 x2 0 c = x1 x2 y1 0 1 d = u1 y1 0
Continuous-time model. >> step(sysnew)
Word文档
.
所以:K=[-1 0] A=[-3 -2;1 0]; B=[1;0]; C=[0 1]; D=0; V=[-3;-3]; sysold=ss(A,B,C,D); p=eig(A) p =-2 -1 Q=obsv(A,C); m=rank(Q);
Word文档
.
n=length(A); if m==n H=acker(A',C',V')' else
disp('系统不是状态完全可观测') end H =-2 3
所以:H=[-2 3]
Word文档