最小公倍数対角分解
[N,D]=lcmdiag(H) [N,D]=lcmdiag(H,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 | ![]() | ![]() |