fig11_27.m

现代控制系统（matlab程序）

% Chapter 11: Figure 11.27, p. 670
%
%   Satellite controllability with a radial thruster only 
%   (i.e. failed tangential thruster).
%
w=0.0011;
A=[0 1 0  0;3*w^2 0 0 2*w;0 0 0 1;0 -2*w 0 0];
b1=[0;1;0;0];
Pc=ctrb(A,b1);
n=det(Pc);  
if abs(n) < eps
 disp('Satellite is uncontrollable with radial thruster only!')
else
 disp('Satellite is controllable with radial thruster only!')
end