ctr_gram 可制御グラミアン

## 呼び出し手順

    [Gc]=ctr_gram(A,B [,dom])
                [Gc]=ctr_gram(sl)

## パラメータ

A,B

  適当な次元の実数行列

dom

  文字 (`'c'` (デフォルト値) または `'d'`)

sl

  線形システム, `syslin` リスト

## 説明

`(A,B)` または `sl`の 可制御グラミアン (`syslin` 線形システム).

`dom` 文字は時間領域を次のように指定します : 離散時間システムの場合 `"d"`,連続時間システムの場合 `"c"` (デフォルト).

## 例

    A=diag([-1,-2,-3]);B=rand(3,2);
    Wc=ctr_gram(A,B)
    U=rand(3,3);A1=U*A/U;B1=U*B;
    Wc1=ctr_gram(A1,B1)    //Not invariant!

## 参照

- equil1
- obs_gram
- contr
- cont_mat
- cont_frm
- contrss