ex3_21.m

来自「反馈控制系统分析与设计,虽然有点抽象」· M 代码 · 共 14 行

G=tf([1,7,24,24],[1,10,35,50,24]); G1=ss(G);
ghinf=normhinf(G1.a,G1.b,G1.c,G1.d)

gamma1=0.8, gamma2=1.2;
A=G1.a; B=G1.b; C=G1.c; D=G1.d;
H_g1=[A, -gamma1^(-2)*B*B'; C'*C, -A'];
H_g2=[A, -gamma2^(-2)*B*B'; C'*C, -A'];
[eig(H_g1), eig(H_g2)]

[m,p,w]=bode(G); % or [m,p,w]=bode(G1);
max(m)

semilogx(w,m(:))