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
参照
Report an issue | ||
<< invr | Polynomials | ldiv >> |