hin_rock.m

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

% H_infinity design for the rocket stabilization system
%
nmeas = 2;
ncon = 1;
gmin = 0.15;
gmax = 1000000;
hin_ic = sel(sys_ic,[8:11],[8:11]);
tol = 0.001;
[K_hin,clp] = hinfsyn(hin_ic,nmeas,ncon,gmin,gmax,tol);
disp(' ')
minfo(K_hin)
disp(' ')
disp('Closed-loop poles')
sp = spoles(clp)
omega = logspace(-3,2,500);
clp_g = frsp(clp,omega);
figure(1)
vplot('liv,lm',sel(vsvd(clp_g),1,1),'m-'), grid
title('Singular Value Plot of clp')
xlabel('Frequency (rad/sec)')
ylabel('Magnitude')
K = K_hin;