imrep2ss
インパルス応答の状態空間実現
呼び出し手順
[sl]=imrep2ss(v [,deg])
パラメータ
- v
インパルス応答のベクトル係数,
v(:,k)はk番目のサンプル- deg
整数 (必要な次数)
- sl
syslinリスト
説明
インパルス応答の線形システムへの変換 (1入力).
v の列数は偶数である必要があります.
例
s=poly(0,'s'); H=[1/(s+0.5);2/(s-0.4)] //厳密にプロパー np=20;w=ldiv(H('num'),H('den'),np); rep=[w(1:np)';w(np+1:2*np)']; //インパルス応答 H1=ss2tf(imrep2ss(rep)) z=poly(0,'z'); H=(2*z^2-3.4*z+1.5)/(z^2-1.6*z+0.8) //プロパーな伝達関数 u=zeros(1,20);u(1)=1; rep=rtitr(H('num'),H('den'),u); //インパルス応答 // <=> rep=ldiv(H('num'),H('den'),20) w=z*imrep2ss(rep) //シフトされたインパルス応答を有する実現 // すなわち厳密にプロパーからプロパーへ H2=ss2tf(w);
参照
| Report an issue | ||
| << findx0BD | Identification | inistate >> |