7.8.3.m

来自「周克敏 鲁棒控制基础 （essential of robust control 」· M 代码 · 共 8 行

T3=tf(conv(conv(conv(conv([-0.05,1],[-0.05,1]),conv([-0.05,1],[-0.05,1])),conv(conv([-0.05,1],[-0.05,1]),conv([-0.05,1],[-0.05,1]))),conv([-0.05,1],[-0.05,1])),conv(conv(conv(conv(conv([0.05,1],[0.05,1]),conv([0.05,1],[0.05,1])),conv(conv([0.05,1],[0.05,1]),conv([0.05,1],[0.05,1]))),conv([0.05,1],[0.05,1])),[0.1,1]));
M3=ss(T3);
[Ab3,Bb3,Cb3,sig3,Tinv3]=balreal(M3.a,M3.b,M3.c);
sig3
N3=[Ab3,Bb3;Cb3,0];
Gb3=pss2sys(N3,11);
Gr3=strunc(Gb3,10)