ztgcq.m

来自「这是一个研究生课程的matlab课件」· M 代码 · 共 25 行

A=[0 1 0;980 0 -2.8;0 0 -100];
B=[0;0 ;100];
C=[1 0 0];
D=[0];
op=[-100;-102;-103];
sysold=ss(A,B,C,D);
disp('原系统的闭还极点为');
p1=eig(A)
n=length(A);
Q=zeros(n);
Q(1,:)=C;
for i=2:n
     Q(i,:)=Q(i-1,:)*A;
 end
m=rank(Q);
 if m==n
     G=place(A',C',op')'
 else
     disp('系统不是状态完全可观测')
 end
 disp('状态观测器模型');
est=estim(sysold,G)
 disp('配置后观测器的极点为');
p2=eig(est)