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);

## 参照

- frep2tf
- arl2
- time_id
- armax
- markp2ss
- ldiv