lyap リアプノフ方程式

## 呼び出し手順

    [X]=lyap(A,C,'c')

                [X]=lyap(A,C,'d')

## 引数

A, C

  実数正方行列, `C` は対称であることが必要

## 説明

`X= lyap(A,C,flag)` は連続時間または離散時間行列 リアプノフ方程式を解きます:

    A'*X + X*A = C       ( flag='c' )
    A'*X*A - X = C       ( flag='d' )

`A`の固有値が`-A`の 固有値でない場合(`flag='c'`) または`A`の固有値分の1の場合 (`flag='d'`)にのみユニークな 解が得られることに注意してください.

## 例

    A=rand(4,4);C=rand(A);C=C+C';
    X=lyap(A,C,'c');
    A'*X + X*A -C
    X=lyap(A,C,'d');
    A'*X*A - X -C

## 参照

- sylv
- ctr_gram
- obs_gram