|
MatrixPolynomialAlgebra[MatrixLCLM] - compute a least-common left multiple of 2 matrices of polynomials
MatrixPolynomialAlgebra[MatrixLCRM] - compute a least-common right multiple of 2 matrices of polynomials
|
|
Calling Sequence
|
|
MatrixLCLM(A, B, x)
MatrixLCRM(A, B, x)
|
|
|
Parameters
|
|
|
A
|
-
|
Matrix of polynomials
|
|
B
|
-
|
Matrix of polynomials
|
|
x
|
-
|
variable name of the polynomial domain
|
|
|
|
|
|
Description
|
|
|
•
|
The MatrixLCLM(A, B, x) command computes a least-common left multiple of two matrices of polynomials. Both input matrices of polynomials can be square or rectangular but must have the same number of columns. The entries are either univariate polynomials in x over the field of rational numbers Q, or rational expressions over Q, that is, univariate polynomials in x with coefficients in Q(a1,...,an).
|
|
|
The matrix of polynomials ![[A^T*B^T]^T](/support/helpjp/helpview.aspx?si=7300/file02411/math85.png) must have more rows than columns and full column rank.
|
|
•
|
The MatrixLCRM(A, B, x) command computes a least-common right multiple of two matrices of polynomials. The matrix of polynomials ![[A*B]](/support/helpjp/helpview.aspx?si=7300/file02411/math91.png) must have more columns than rows and full row rank.
|
|
|
The method is a fraction-free algorithm by Beckermann and Labahn that computes a matrix GCD using Mahler systems.
|
|
|
|
Examples
|
|
|
>
|
|
|
>
|
![A := Matrix(2, 2, [[-9*z^2-3*z+1, 12*z^2+10*z], [-3*z^3+2*z^2-z, 4*z^3+2*z-2*z^2]])](/support/helpjp/helpview.aspx?si=7300/file02411/math106.png)
|
|
>
|
![B := Matrix(2, 2, [[-3*z^3+6*z^2+5*z+1, -12*z^2-13*z], [z^4+z^3+z^2, -4*z^3-3*z+3*z^2]])](/support/helpjp/helpview.aspx?si=7300/file02411/math110.png)
|
Left matrix LCMs:
|
>
|
|
|
![C := Matrix([[66*z^2+108*z^5-228*z^4-92*z^3+27*z+2, -128*z^2-144*z^5+504*z^4+28*z^3-26*z], [36*z^5-84*z^4-24*z^3+26*z^2+7*z, -48*z^5+184*z^4-12*z^3-40*z^2-6*z]])](/support/helpjp/helpview.aspx?si=7300/file02411/math119.png)
| (1) |
Right matrix LCMs:
|
>
|
|
|
![C := Matrix([[12*z^2+60*z^4+59*z^3+2+10*z, 12*z^3+z^2-13*z], [20*z^5-13*z^4+2*z^2+17*z^3, 4*z^4-7*z^3+6*z^2-3*z]])](/support/helpjp/helpview.aspx?si=7300/file02411/math128.png)
| (2) |
|
|
|
References
|
|
|
|
Beckermann, B., and Labahn, G. "Fraction-free Computation of Matrix Rational Interpolants and Matrix GCDs." SIAM Journal on Matrix Analysis and Applications, Vol. 22 No.1, (2000): 114-144.
|
|
|
Download Help Document
Was this information helpful?
