例
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 |  |  |
最小公倍数対角分解
[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 | | |