85,689 research outputs found
Fraction-free algorithm for the computation of diagonal forms matrices over Ore domains using Gr{\"o}bner bases
This paper is a sequel to "Computing diagonal form and Jacobson normal form
of a matrix using Groebner bases", J. of Symb. Computation, 46 (5), 2011. We
present a new fraction-free algorithm for the computation of a diagonal form of
a matrix over a certain non-commutative Euclidean domain over a computable
field with the help of Gr\"obner bases. This algorithm is formulated in a
general constructive framework of non-commutative Ore localizations of
-algebras (OLGAs). We split the computation of a normal form of a matrix
into the diagonalization and the normalization processes. Both of them can be
made fraction-free. For a matrix over an OLGA we provide a diagonalization
algorithm to compute and with fraction-free entries such that
holds and is diagonal. The fraction-free approach gives us more information
on the system of linear functional equations and its solutions, than the
classical setup of an operator algebra with rational functions coefficients. In
particular, one can handle distributional solutions together with, say,
meromorphic ones. We investigate Ore localizations of common operator algebras
over and use them in the unimodularity analysis of transformation
matrices . In turn, this allows to lift the isomorphism of modules over an
OLGA Euclidean domain to a polynomial subring of it. We discuss the relation of
this lifting with the solutions of the original system of equations. Moreover,
we prove some new results concerning normal forms of matrices over non-simple
domains. Our implementation in the computer algebra system {\sc
Singular:Plural} follows the fraction-free strategy and shows impressive
performance, compared with methods which directly use fractions. Since we
experience moderate swell of coefficients and obtain simple transformation
matrices, the method we propose is well suited for solving nontrivial practical
problems.Comment: 25 pages, to appear in Journal of Symbolic Computatio
Computing diagonal form and Jacobson normal form of a matrix using Gr\"obner bases
In this paper we present two algorithms for the computation of a diagonal
form of a matrix over non-commutative Euclidean domain over a field with the
help of Gr\"obner bases. This can be viewed as the pre-processing for the
computation of Jacobson normal form and also used for the computation of Smith
normal form in the commutative case. We propose a general framework for
handling, among other, operator algebras with rational coefficients. We employ
special "polynomial" strategy in Ore localizations of non-commutative
-algebras and show its merits. In particular, for a given matrix we
provide an algorithm to compute and with fraction-free entries such
that holds. The polynomial approach allows one to obtain more precise
information, than the rational one e. g. about singularities of the system.
Our implementation of polynomial strategy shows very impressive performance,
compared with methods, which directly use fractions. In particular, we
experience quite moderate swell of coefficients and obtain uncomplicated
transformation matrices. This shows that this method is well suitable for
solving nontrivial practical problems. We present an implementation of
algorithms in SINGULAR:PLURAL and compare it with other available systems. We
leave questions on the algorithmic complexity of this algorithm open, but we
stress the practical applicability of the proposed method to a bigger class of
non-commutative algebras
Fast, deterministic computation of the Hermite normal form and determinant of a polynomial matrix
Given a nonsingular matrix of univariate polynomials over a
field , we give fast and deterministic algorithms to compute its
determinant and its Hermite normal form. Our algorithms use
operations in ,
where is bounded from above by both the average of the degrees of the rows
and that of the columns of the matrix and is the exponent of matrix
multiplication. The soft- notation indicates that logarithmic factors in the
big- are omitted while the ceiling function indicates that the cost is
when . Our algorithms are based
on a fast and deterministic triangularization method for computing the diagonal
entries of the Hermite form of a nonsingular matrix.Comment: 34 pages, 3 algorithm
Denominator Bounds and Polynomial Solutions for Systems of q-Recurrences over K(t) for Constant K
We consider systems A_\ell(t) y(q^\ell t) + ... + A_0(t) y(t) = b(t) of
higher order q-recurrence equations with rational coefficients. We extend a
method for finding a bound on the maximal power of t in the denominator of
arbitrary rational solutions y(t) as well as a method for bounding the degree
of polynomial solutions from the scalar case to the systems case. The approach
is direct and does not rely on uncoupling or reduction to a first order system.
Unlike in the scalar case this usually requires an initial transformation of
the system.Comment: 8 page
- âŠ