Scilabでベクトル軌跡を書くワンライナー

Gのベクトル軌跡
:sをjωとおいてωを0→∞に飛ばした時のGの複素平面上での軌跡

微分要素

plot2d(real(horner(1/(1/%s),(0:0.01:100)*%i)),imag(horner(1/(1/%s),(0:0.01:100)*%i)),axesflag=1)

積分要素

plot2d(real(horner(1/%s,(0:0.01:100)*%i)),imag(horner(1/%s,(0:0.01:100)*%i)),axesflag=1)

1次遅れ系

plot2d(real(horner(1/(1+%s),(0:0.01:100)*%i)),imag(horner(1/(1+%s),(0:0.01:100)*%i)),axesflag=1)

1次遅れの微分要素

plot2d(real(horner(%s/(1+%s),(0:0.01:100)*%i)),imag(horner(%s/(1+%s),(0:0.01:100)*%i)),axesflag=1)

2次遅れ系

plot2d(real(horner(1/(%s^2+2*%s+1),(0:0.01:100)*%i)),imag(horner(1/(%s^2+2*%s+1),(0:0.01:100)*%i)),axesflag=1)

2次遅れ系の微分要素

plot2d(real(horner(%s/(%s^2+2*%s+1),(0:0.01:100)*%i)),imag(horner(%s/(%s^2+2*%s+1),(0:0.01:100)*%i)),axesflag=1)

位相進み要素

plot2d(real(horner((1+2*%s)/(1+%s),(0:0.01:100)*%i)),imag(horner((1+2*%s)/(1+%s),(0:0.01:100)*%i)),axesflag=1)

位相遅れ要素

plot2d(real(horner((1+%s)/(1+2*%s),(0:0.01:100)*%i)),imag(horner((1+%s)/(1+2*%s),(0:0.01:100)*%i)),axesflag=1)

一覧

clf()
subplot(4,2,1)
plot2d(real(horner(1/(1/%s),(0:0.01:100)*%i)),imag(horner(1/(1/%s),(0:0.01:100)*%i)),axesflag=1)
subplot(4,2,2)
plot2d(real(horner(1/%s,(0:0.01:100)*%i)),imag(horner(1/%s,(0:0.01:100)*%i)),axesflag=1)
subplot(4,2,3)
plot2d(real(horner(1/(1+%s),(0:0.01:100)*%i)),imag(horner(1/(1+%s),(0:0.01:100)*%i)),axesflag=1)
subplot(4,2,4)
plot2d(real(horner(%s/(1+%s),(0:0.01:100)*%i)),imag(horner(%s/(1+%s),(0:0.01:100)*%i)),axesflag=1)
subplot(4,2,5)
plot2d(real(horner(1/(%s^2+2*%s+1),(0:0.01:100)*%i)),imag(horner(1/(%s^2+2*%s+1),(0:0.01:100)*%i)),axesflag=1)
subplot(4,2,6)
plot2d(real(horner(%s/(%s^2+2*%s+1),(0:0.01:100)*%i)),imag(horner(%s/(%s^2+2*%s+1),(0:0.01:100)*%i)),axesflag=1)
subplot(4,2,7)
plot2d(real(horner((1+2*%s)/(1+%s),(0:0.01:100)*%i)),imag(horner((1+2*%s)/(1+%s),(0:0.01:100)*%i)),axesflag=1)
subplot(4,2,8)
plot2d(real(horner((1+%s)/(1+2*%s),(0:0.01:100)*%i)),imag(horner((1+%s)/(1+2*%s),(0:0.01:100)*%i)),axesflag=1)