bode_asymp

s = poly(0, "s");
h = syslin("c", (s^2+2*0.9*10*s+100)/(s^2+2*0.3*10.1*s+102.01));
clf(); bode(h, 0.01, 100); bode_asymp(h);

s = poly(0, "s");
h1 = syslin("c", (s^2+2*0.9*10*s+100)/(s^2+2*0.3*10.1*s+102.01));
h2 = syslin("c", (10*(s+3))/(s*(s+2)*(s^2+s+2)));
clf(); bode([h1; h2], 0.01, 100, ["h1"; "h2"]); bode_asymp([h1; h2]);

s = %s;
G = (s+10)/(s*(s+100)); // 有理行列
sys = syslin("c", G); // 伝達関数表現の連続時間線形系
f_min = .0001; f_max = 1600; // 単位: Hz
clf(); bode(sys, f_min, f_max, "rad"); // オプション引数 "rad" は Hzをrad/sに変換します
bode_asymp(sys, 10, 1000); // 指定した範囲(単位:rad/s)で漸近線をプロット

s = poly(0, "s");
h = syslin("c", (s^2+2*0.9*10*s+100)/(s^2+2*0.3*10.1*s+102.01));
h = tf2ss(h);
clf();
gainplot(h, 0.01, 100, "$\frac{s^2+18 s+100}{s^2+6.06 s+102.1}$"); // 単位: ヘルツ
bode_asymp(h)