% ex8_16
% 確定例8-15中原系統狀態和降階觀測器狀態的軌線
eval('ex8_15');
ao=[a  zeros(2,1); L*c  F];
bo=[b; J];
co=eye(3);
do=[0; 0; 0];
x0=[1;0;0];
t=[0:0.01:4];
u=4*ones(length(t),1); 
stp=100;
u([stp:length(t)],1)=8*ones((length(t)+1-stp),1); 
[y,x]=lsim(ao,bo,co,do,u,t,x0);
est=zeros(length(t),2);
figure('pos',[50 100,200,200],'color','w');
for k = 1:length(t),
    est(k,:)=[Iphe*[y(k,1); y(k,3)]]';
end;
y=[y est];
plot(t,y(:,[1 2 4 5])),
axis([0 4 -1 10]);
set(gca,'fonts',8,'fontw','b')