# lcmdiag

最小公倍数対角分解

## 呼び出し手順

    [N,D]=lcmdiag(H)
                [N,D]=lcmdiag(H,flag)

## パラメータ

H

  有理行列

N

  多項式行列

D

  対角多項式行列

flag

  文字列: `'row'` または `'col'` (デフォルト)

## 説明

`[N,D]=lcmdiag(H,'row')` は, 分解`D*H=N`, すなわち,`H=D^(-1)*N`を計算します ただし, D は D(k,k)= H('den')のk行目の最小公倍数となる 対角行列です.

`[N,D]=lcmdiag(H)` または `[N,D]=lcmdiag(H,'col)` は,D(k,k)=H('den')のk列目の最小公倍数となる`H=N*D^(-1)`を返します.

## 例

    s=poly(0,'s');
    H=[1/s,(s+2)/s/(s+1)^2;1/(s^2*(s+2)),2/(s+2)];
    [N,D]=lcmdiag(H);
    N/D-H

## 参照

- lcm
- gcd
- bezout