hankelsv ハンケル特異値

## 呼び出し手順

    [nk2,W]=hankelsv(sl [,tol])
                [nk2]=hankelsv(sl [,tol])

## パラメータ

sl

  線形システム(状態空間)を表す `syslin` リスト.

tol

  虚軸モードを検出するための許容誤差パラメータ (デフォルト値: `1000*%eps`).

## 説明

`sl`の ハンケル特異値の二乗 `nk2`と, `W = P*Q` = 可制御性グラミアンと 可観測性グラミアンの積を返します.

`nk2` は`W`の 固有値のベクトルです.

## 例

    A=diag([-1,-2,-3]);
    sl=syslin('c',A,rand(3,2),rand(2,3));[nk2,W]=hankelsv(sl)
    [Q,M]=pbig(W,nk2(2)-%eps,'c');
    slr=projsl(sl,Q,M);hankelsv(slr)

## 参照

- balreal
- equil
- equil1