技術士(総合技術監理・機械部門)のブログ

技術士がお届けする,技術,工学,技術者倫理などの話題

SCILABによる制御系設計 その7

2017-03-17 20:21:52 | 工学
Scilab7
3 過渡応答解析 (その4)

インパルス応答をやってみよう

plot2d([t',t'],[(csim('impulse',t,des2ss(A,B,C,D,E)))',0*t'])

でいけるはず。

-->A=[0 1;-1 -1];

-->B=[0;1];

-->C=[1 0];

-->D=[0];

-->E=[1 0;0 1];

-->t=0:0.01:10;

-->plot2d([t',t'],[(csim('impuls',t,des2ss(A,B,C,D,E)))',0*t'])
警告: csim: 入力引数 #1 は連続時間系で指定してください.




任意の入力の場合、


s=poly(0,'s');

-->sy=5/(s^2+s+5);

として

-→u=2+t;

つまりランプ入力

plot2d(csim(u,t,tf2ss(sy)));

とすれば任意の入力uに対する応答が得られる。


SCILABによる制御系設計 その7

2017-03-17 20:21:52 | 工学
Scilab7
3 過渡応答解析 (その4)

インパルス応答をやってみよう

plot2d([t',t'],[(csim('impulse',t,des2ss(A,B,C,D,E)))',0*t'])

でいけるはず。

-->A=[0 1;-1 -1];

-->B=[0;1];

-->C=[1 0];

-->D=[0];

-->E=[1 0;0 1];

-->t=0:0.01:10;

-->plot2d([t',t'],[(csim('impuls',t,des2ss(A,B,C,D,E)))',0*t'])
警告: csim: 入力引数 #1 は連続時間系で指定してください.




任意の入力の場合、


s=poly(0,'s');

-->sy=5/(s^2+s+5);

として

-→u=2+t;

つまりランプ入力

plot2d(csim(u,t,tf2ss(sy)));

とすれば任意の入力uに対する応答が得られる。