freq 周波数応答

## 呼び出し手順

    [x]=freq(A,B,C [,D],f)
                [x]=freq(NUM,DEN,f)

## パラメータ

A, B, C, D

  適当な次元`nxn, nxp, mxn, mxp`の実数行列 .

NUM,DEN

  次元 `mxp`の多項式行列

x

  実数または複素数の行列

## 説明

`x=freq(A,B,C [,D],f)` は次の関係を満たす実数または複素数の `mxp*t` 行列を返します.:

`x(:,k*p:(k+1)*p)= C*inv(f(k)*eye()-A)*B + D`.

ここで, `f` が虚軸上または単位円上に 値を有する場合, `x`は連続系または離散系における `(A,B,C,D)`の周波数応答となります.

`x=freq(NUM,DEN,f)` は, `x` の `k*(p-1)+1`列から `k*p`列が行列`NUM(f(k))./DEN(f(k))`を 含むような 実数または複素数の 行列`x`を返します.

## 例

    s=poly(0,'s');
    sys=(s+1)/(s^3-5*s+4)
    rep=freq(sys("num"),sys("den"),[0,0.9,1.1,2,3,10,20])
    [horner(sys,0),horner(sys,20)]
    //
    Sys=tf2ss(sys);
    [A,B,C,D]=abcd(Sys);
    freq(A,B,C,[0,0.9,1.1,2,3,10,20])

## 説明

- repfreq
- horner