Quadratization of ODEs: Monomial vs. Non-Monomial

Abstract

Quadratization is a transform of a system of ODEs with polynomial right-hand side into a system of ODEs with at most quadratic right-hand side via the introduction of new variables. It has been recently used as a pre-processing step for new model order reduction methods, so it is important to keep the number of new variables small. Several algorithms have been designed to search for a quadratization with the new variables being monomials in the original variables. To understand the limitations and potential ways of improving such algorithms, we study the following question: can quadratizations with not necessarily monomial new variables produce a model of substantially smaller dimension than quadratization with only monomial new variables? To do this, we restrict our attention to scalar polynomial ODEs. Our first result is that a scalar polynomial ODE x=p(x)=anxn+an-1xn-1+… + a0 with n≥slant 5 and an≠0 can be quadratized using exactly one new variable if and only if p(x-an-1n· an)=anxn+ax2+bx for some a, b ∈ C. In fact, the new variable can be taken z:=(x-an-1n· an)n-1. Our second result is that two non-monomial new variables are enough to quadratize all degree 6 scalar polynomial ODEs. Based on these results, we observe that a quadratization with not necessarily monomial new variables can be much smaller than a monomial quadratization even for scalar ODEs. The main results of the paper have been discovered using computational methods of applied nonlinear algebra (Gr\"obner bases), and we describe these computations.

0

Discussion (0)

Sign in to join the discussion.

Loading comments…