lcmdiag

パラメータ

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 — 多項式またｈ整数のべズー方程式