gcq.m

这是一个研究生课程的matlab课件

%对应的观测器极点为-50 
a=[0 1 0 0;0 -5 -100 -9.14;95 0 -95 1;2025 0 -2025 -5];
 b=[0 100 0 100]';
 c=[1 0 0 0;0 0 1 0];
 d=0;
 t=0:0.1:10;
 t1=length(t);
 u=ones(t1,1)*0.1;
 [y,x]=lsim(a,b,c,d,u,t)
 keyboard
%  x0=[0 0 0 0];
% [y,x]=lsim(a,b,c,d,u,t,x0)
% plot(t,x(:,1),'-r',t,x(:,2),'-.g',t,x(:,3),'-b',t,x(:,4),'-.m')
x01=[0 0 0.15 0.15];
[y1,x1]=lsim(a,b,c,d,u,t,x01)
% a11=[0 1;-100,-14.14];b11=[0 100]';c11=[1 0];d=0;
% [y11,x11]=lsim(a11,b11,c11,d,u,t)
plot(t,y,'-r',t,y1,'-.g')
%对应的观测器极点为-10
%  a1=[0 1 0 0;0 -5 -100 -9.14;15 0 -15 1;25 0 -25 -5];
%  b=[0 100 0 100]';
%  c=[1 0 0 0];
%  d=0;
%  t=0:0.1:10;
%  t1=length(t);
%  u=ones(t1,1)*0.1;
%  [y1,x1]=lsim(a1,b,c,d,u,t)
% %  plot(t,x1)
%  plot(t,y,t,y1)
%  eig(a1)