17,324 research outputs found
Roots of bivariate polynomial systems via determinantal representations
We give two determinantal representations for a bivariate polynomial. They
may be used to compute the zeros of a system of two of these polynomials via
the eigenvalues of a two-parameter eigenvalue problem. The first determinantal
representation is suitable for polynomials with scalar or matrix coefficients,
and consists of matrices with asymptotic order , where is the degree
of the polynomial. The second representation is useful for scalar polynomials
and has asymptotic order . The resulting method to compute the roots of
a system of two bivariate polynomials is competitive with some existing methods
for polynomials up to degree 10, as well as for polynomials with a small number
of terms.Comment: 22 pages, 9 figure
Multiparameter spectral analysis for aeroelastic instability problems
This paper presents a novel application of multiparameter spectral theory to
the study of structural stability, with particular emphasis on aeroelastic
flutter. Methods of multiparameter analysis allow the development of new
solution algorithms for aeroelastic flutter problems; most significantly, a
direct solver for polynomial problems of arbitrary order and size, something
which has not before been achieved. Two major variants of this direct solver
are presented, and their computational characteristics are compared. Both are
effective for smaller problems arising in reduced-order modelling and
preliminary design optimization. Extensions and improvements to this new
conceptual framework and solution method are then discussed.Comment: 20 pages, 8 figure
The linearization problem of a binary quadratic problem and its applications
We provide several applications of the linearization problem of a binary
quadratic problem. We propose a new lower bounding strategy, called the
linearization-based scheme, that is based on a simple certificate for a
quadratic function to be non-negative on the feasible set. Each
linearization-based bound requires a set of linearizable matrices as an input.
We prove that the Generalized Gilmore-Lawler bounding scheme for binary
quadratic problems provides linearization-based bounds. Moreover, we show that
the bound obtained from the first level reformulation linearization technique
is also a type of linearization-based bound, which enables us to provide a
comparison among mentioned bounds. However, the strongest linearization-based
bound is the one that uses the full characterization of the set of linearizable
matrices. Finally, we present a polynomial-time algorithm for the linearization
problem of the quadratic shortest path problem on directed acyclic graphs. Our
algorithm gives a complete characterization of the set of linearizable matrices
for the quadratic shortest path problem
Perturbation, extraction and refinement of invariant pairs for matrix polynomials
Generalizing the notion of an eigenvector, invariant subspaces are frequently used in the context of linear eigenvalue problems, leading to conceptually elegant and numerically stable formulations in applications that require the computation of several eigenvalues and/or eigenvectors. Similar benefits can be expected for polynomial eigenvalue problems, for which the concept of an invariant subspace needs to be replaced by the concept of an invariant pair. Little has been known so far about numerical aspects of such invariant pairs. The aim of this paper is to fill this gap. The behavior of invariant pairs under
perturbations of the matrix polynomial is studied and a first-order perturbation expansion is given. From a computational point of view, we investigate how to best extract invariant pairs from a linearization of the matrix polynomial. Moreover, we describe efficient refinement procedures directly based on the polynomial formulation. Numerical experiments
with matrix polynomials from a number of applications demonstrate the effectiveness of our extraction and refinement procedures
A framework for structured linearizations of matrix polynomials in various bases
We present a framework for the construction of linearizations for scalar and
matrix polynomials based on dual bases which, in the case of orthogonal
polynomials, can be described by the associated recurrence relations. The
framework provides an extension of the classical linearization theory for
polynomials expressed in non-monomial bases and allows to represent polynomials
expressed in product families, that is as a linear combination of elements of
the form , where and
can either be polynomial bases or polynomial families
which satisfy some mild assumptions. We show that this general construction can
be used for many different purposes. Among them, we show how to linearize sums
of polynomials and rational functions expressed in different bases. As an
example, this allows to look for intersections of functions interpolated on
different nodes without converting them to the same basis. We then provide some
constructions for structured linearizations for -even and
-palindromic matrix polynomials. The extensions of these constructions
to -odd and -antipalindromic of odd degree is discussed and
follows immediately from the previous results
- …