Scilab Home page | Wiki | Bug tracker | Forge | Mailing list archives | ATOMS | File exchange
Please login or create an account
Change language to: English - Français - Português -

Please note that the recommended version of Scilab is 6.0.1. This page might be outdated.
See the recommended documentation of this function

Scilab help >> CACSD > csim

csim

線形システムのシミュレーション (時間応答)

呼び出しの手順

[y [,x]]=csim(u,t,sl,[x0 [,tol]])

パラメータ

u

関数, リストまたは文字列 (制御入力)

t

時間を指定するための実数ベクトルで、t(1) は 初期時間 (x0=x(t(1)))を表す.

sl

リスト (syslin)

y

y=[y(t(i)], i=1,..,n となる行列

x

x=[x(t(i)], i=1,..,n となる行列

tol

2つの要素 [atol rtol] からなるベクトルであり、それぞれ ODEソルバ(ode参照)の絶対許容誤差および相対許容誤差を定義する

説明

線形制御系 slのシミュレーションを行う. slは、 syslinリストで表された連続時間システムの 入力を前提とする.

u は制御入力、x0 は状態量初期値である.

y は出力、x は状態量である.

制御入力は以下のいずれかとすることができる:

1. 関数 : [inputs]=u(t)

2. リスト : list(ut,parameter1,....,parametern). ただし、inputs=ut(t,parameter1,....,parametern) (ut は関数)

3. インパルス応答の計算を表す文字列 "impuls" (ここで、sl は直達項を有さず、x0=0である SISO系であると仮定)

4. ステップ応答の計算を表す文字列 "step" (ここで、sl は直達項を有さず、x0=0である SISO系であると仮定)

5. t の各値に対応する u の値を指定するベクトル.

s=poly(0,'s');rand('seed',0);w=ssrand(1,1,3);w('A')=w('A')-2*eye();
t=0:0.05:5;
//impulse(w) = step (s * w)
clf(0);xset("window",0);show_window();
plot2d([t',t'],[(csim('step',t,tf2ss(s)*w))',0*t'])
clf(1);xset("window",1);show_window();
plot2d([t',t'],[(csim('impulse',t,w))',0*t'])
//step(w) = impulse (s^-1 * w)
clf(3);xset("window",3);show_window();
plot2d([t',t'],[(csim('step',t,w))',0*t'])
clf(4);xset("window",4);show_window();
plot2d([t',t'],[(csim('impulse',t,tf2ss(1/s)*w))',0*t'])

//input defined by a time function
deff('u=input(t)','u=abs(sin(t))')
clf();plot2d([t',t'],[(csim(input,t,w))',0*t'])

参照

  • syslin — 線形システムを定義する
  • dsimul — 離散時間状態空間シミュレーション
  • flts — 時間応答 (離散時間, 離散化システム)
  • ltitr — 離散時間応答 (状態空間)
  • rtitr — 離散時間応答 (伝達行列)
  • ode — 常微分方程式ソルバ
  • impl — 微分代数方程式
Scilab Enterprises
Copyright (c) 2011-2017 (Scilab Enterprises)
Copyright (c) 1989-2012 (INRIA)
Copyright (c) 1989-2007 (ENPC)
with contributors
Last updated:
Thu Mar 03 11:00:53 CET 2011