red_rock.m

来自「导弹控制系统的鲁棒控制设计实例」· M 代码 · 共 19 行

% Controller order reduction for the rocket stabilization system
%
omega = logspace(-6,4,100);
K_g = frsp(K,omega);
[Kb,hsig] = sysbal(K);
n_red = 8;
Kred = hankmr(Kb,hsig,n_red,'d');
Kred_g = frsp(Kred,omega);
%
% frequency responses of the controllers
figure(1)
vplot('liv,lm',vsvd(K_g),'r-',vsvd(Kred_g),'c--'), grid
tmp1 = 'Singular values of the';
tmp2 = ' controller transfer matrices';
title([tmp1 tmp2])
xlabel('Frequency (rad/sec)')
%legend('Full-order (n = 28) controller', ...
%       'Reduced-order (n = 8) controller',3)