双線形変換
[sl1]=cls2dls(sl,T [,fp])
線形システム (syslin リスト)
syslin
実数, サンプリング周期
プリワープ周波数 (単位: Hz)
sl=[A,B,C,D] (syslin リスト), 連続時間システム cls2dls は双線形変換s=(2/T)*(z-1)/(z+1) により得られたサンプル値システムを返します.
sl=[A,B,C,D]
cls2dls
s=(2/T)*(z-1)/(z+1)
s=poly(0,'s');z=poly(0,'z'); sl=syslin('c',(s+1)/(s^2-5*s+2)); //連続時間系伝達関数表現 slss=tf2ss(sl); //状態空間表現に変換 sl1=cls2dls(slss,0.2); //sl1= cls2dlsの出力 sl1t=ss2tf(sl1) // 伝達関数表現に変換 sl2=horner(sl,(2/0.2)*(z-1)/(z+1)) //sl2 と sl1を比較