DEtools - Maple Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Differential Equations : Lie Symmetry Method : Commands for ODEs : DEtools/power_equivalent

DEtools

  

power_equivalent

  

computes an equivalent ODE with invariants with minimal degrees with respect to power transformations

  

rational_equivalent

  

computes an equivalent ODE with invariants with minimal degrees with respect to rational transformations

 

Calling Sequence

Parameters

Description

Notes

Examples

Calling Sequence

power_equivalent(ODE, y(x), u(t), keep)

rational_equivalent(ODE, y(x), u(t), keep)

Parameters

ODE

-

linear ODE of any differential order or a first order nonlinear ODE of Abel or Riccati type

y(x)

-

dependent variable of the ODE; required only if the given ODE contains more than one function that is differentiated

u(t)

-

(optional); specify new variables t,ut used in the returned rational ODE

keep

-

(optional); specify that the variables used in the returned rational ODE be those of the given ODE.

Description

• 

The power_equivalent and rational_equivalent commands return an ODE (equivalent to the given ODE), with rational coefficients, and the transformation that maps the returned ODE back to the given ODE. The returned ODE has some additional relevant features.

1. 

It is in normal form.

2. 

All its invariants are rational.

3. 

The degree of each invariant with respect to the independent variable is the lowest one achievable using power (power_equivalent) or rational (rational_equivalent) transformations.

• 

The use of a NormalForm where in addition item 3 holds makes it possible to solve equivalence problems under compositions of Mobius and power transformations between a given ODE and ODEs admitting, for instance, hypergeometric solutions (see dsolve/hyper3).

• 

The transformation used to compute the equivalent ODE composes arbitrary powers (for power_equivalent) or rational functions (for rational_equivalent) of the independent variable with the most general linear transformation of the dependent variable, so they are of the form x=ta,yx=Ptut+Qt. In the case of linear homogeneous ODEs, Q=0. In terms of linear ODEs, the rational_equivalent command only handles those with order > 2. The transformations preserve the polynomial degree, the linear or nonlinear type and the homogeneous character in the case of linear equations. For example, second order linear equations are mapped into equations of the same type, the same way as Abel or Riccati type nonlinear equations are mapped into equations of the same type.

• 

By comparing the transformations used by power_equivalent and rational_equivalent with the general transformation x=Ft,yx=Ptut+Qt which defines the equivalence class associated to any polynomial ODE (linear or not), you can see that the restrictions imposed by power_equivalent and rational_equivalent respectively consist of taking Ft=ta and F(t)=rational function of t. Due to these restrictions, for some given ODE, the equivalence with a rational coefficients ODE may not exist. However, if in the given ODE, when written in normal form (see convert/NormalForm), all the nonrational objects are either fractional or abstract powers of the independent variable, then the equivalence problem has a solution, which the power_equivalent command should find. Caveat: Despite the theory behind the algorithm, in the case of abstract powers, the success relies on the ability of the simplification routines to rewrite these powers in a normal form. (This may be a difficult problem; see evala/Normal.)

• 

When the equivalence does not exist power_equivalent and rational_equivalent return NULL. Otherwise, the output sequence consists of a rational ODE and the transformation mapping it into the given ODE. New variables of global type are introduced in the result.

  

Optionally, you can specify that new dependent ut and independent t variables be used in the output by including the option ut in the calling sequence.

  

By using the keep option, you specify that the returned rational ODE should be expressed using the original variables x,yx. In that case, the returned transformation maps the returned rational ODE into the given ODE expressed using the new variables t,ut.

  

Note: The keep option supersedes the u(t) parameter. If superseded by the keep parameter, t,ut represents the original variables in the transformation back to the given ODE. (See the examples.)

Notes

• 

Parts of the code for the command rational_equivalent were written by Austin Roche.

Examples

withDEtools:

Linear ODEs

As the simplest example, consider a first order linear homogeneous ODE.

ⅆⅆxyx+fxyx=0

ⅆⅆxyx+fxyx=0

(1)

power_equivalent

ⅆⅆ_a_b_a=0,_a=x,_b_a=yxⅇ∫fxⅆx

(2)

In this output we see the "rational normal form with minimal degrees" followed by the transformation mapping it into the given ODE. To perform this mapping and reobtain the given ODE, use dchange.

factorPDEtools[dchange]2,1,x,yx

ⅇ∫fxⅆxⅆⅆxyx+fxyx=0

(3)

A second order linear ODE with nonrational coefficients:

non_rational_odeⅆ2ⅆx2yx=1213x56+2x23yx4x+11+x16ⅆⅆxyx2x

non_rational_ode:=ⅆ2ⅆx2yx=1423x5/6+2x2/3yxx+121+x1/6ⅆⅆxyxx

(4)

rational_ode,trpower_equivalentnon_rational_ode,ut

rational_ode,tr:=ⅆ2ⅆt2ut=814+9t+2t2ut,t=x1/6,ut=yxⅇ∫141+x1/6xⅆxx5/12

(5)

Solve the equivalent rational ODE, and then change variables using the returned transformation to obtain the solution to the given ODE.

sol_rationaldsolverational_ode

sol_rational:=ut=_C1ⅇ92tt+_C2ⅇ92t81t218t+2t

(6)

sol_non_rationalevalsol_rational,tr

sol_non_rational:=yxⅇ∫141+x1/6xⅆxx5/12=_C1ⅇ92x1/6x1/6+_C2ⅇ92x1/681x1/318x1/6+2x1/6

(7)

Isolating yx, obtain the solution in explicit form.

simplifyisolate,yx

yx=18_C2ⅇ9x1/6x1/681_C2ⅇ9x1/6x1/32_C2ⅇ9x1/6_C1x1/4ⅇ14∫1+x1/6xⅆx92x1/6

(8)

odetest,non_rational_ode

0

(9)

The method works with abstract powers as well. (The 'keep' option specifies that the returned rational ODE be in yx.)

ⅆ2ⅆx2yx=1x121n2n+2xnyx4x+11+x1n2+12ⅆⅆxyx2x

ⅆ2ⅆx2yx=14x1212nn+2xnyxx+121+x12n+12ⅆⅆxyxx

(10)

power_equivalent,ut,keep

ⅆ2ⅆx2yx=94n12+32nn12x14n22nn12x2yx,x=t12n+12,yx=t1414nutⅇ∫141+t12n+12tⅆt

(11)

A higher order linear ODE example:

ⅆ5ⅆx5yx=5ⅆ4ⅆx4yxx151ⅆ3ⅆx3yx4x2+1ayx32x152

ⅆ5ⅆx5yx=5ⅆ4ⅆx4yxx154ⅆ3ⅆx3yxx2+132ayxx15/2

(12)

power_equivalent,ut,keep

ⅆ5ⅆx5yx=19008625x5+13125ax6yx18504625ⅆⅆxyxx4+725ⅆ2ⅆx2yxx3245ⅆ3ⅆx3yxx2,x=t5/2,yx=t3utⅇ∫1tⅆt

(13)

For the following example, already having rational coefficients, note the difference between its Normal form (see convert/NormalForm) and the Normal form with minimal degrees returned by power_equivalent.

odeⅆ5ⅆx5yx=xyx

ode:=ⅆ5ⅆx5yx=xyx

(14)

convertode,NormalForm,keep

ⅆ5ⅆx5yx=xyx,_ax=_C1yx

(15)

power_equivalentode,keep

ⅆ5ⅆx5yx=17776x4+5005162x5yx4865162ⅆⅆxyxx4+17512ⅆ2ⅆx2yxx317536ⅆ3ⅆx3yxx2,x=_a6,yx=_a10_b_a

(16)

Recalling the definition, a linear ODE of order n is in normal form when the coefficient of the (n-1)th derivative is zero. For a given ODE, there are many normal forms while the normal form with minimal degrees (one of these possible normal forms) is unique in that all the invariants (the coefficients of yx and its derivative) have the lowest degrees achievable through power transformations. Here, the term degree of a rational expression means the degree of its numerator minus the degree of its denominator.

First order ODEs of Abel type

A first order ODE of Abel type is one that, when written in First Kind form, is of degree 3 in the dependent variable (see odeadvisor,Abel). For example,

ⅆⅆxyx=3x12+1xx12+1yx32x11215yx24x134

ⅆⅆxyx=32x+1xx+1yx3x11/2154yx2x13/4

(17)

As in the linear ODE case, the Normal form returned by power_equivalent is one of the possible Rational Normal Forms, that is, where the coefficient of yx2 is zero and the independent term is equal to 1 (see convert/Abel_RNF). The Normal form by power_equivalent however always has rational coefficients and is of lowest degree (with respect to power transformations) in the independent variable.

power_equivalent

ⅆⅆ_a_b_a=125116642+9_a2_b_a3_a3_a+13+54_a2+99_a+38_b_a108_a3+132_a2+24_a,_a=x3/2,_b_a=1856x3yx+5x15/4+12x3/2yx+5x9/4+6yxx3/49x3/2+2

(18)

An example where the input Abel ODE is in second kind form (see convert/FirstKind):

ⅆⅆxyx=1+12x72+ax52yx

ⅆⅆxyx=1+12x7/2+ax5/2yx

(19)

power_equivalent

ⅆⅆ_a_b_a=5025004744_a9a2_b_a3_a212_a+a3+2217172_a22931_aa27a2_b_a_a44_a9a12_a+a+1,_a=x7/2,_b_a=635x12yxx6+432x7+yxx5/2a+72x7/2a+3a2yx44x7/29a

(20)

In the following example, the Abel ODE has rational coefficients, but its Rational Normal Form (see convert/Abel_RNF) is different from the normal form returned from power_equivalent in that the degrees in x of the invariants in the canonical form are significantly lower.

ⅆⅆxyx=6x5+Ayx3x4yx2

ⅆⅆxyx=6x5+Ayx3x4yx2

(21)

convert,Abel_RNF,keep

ⅆⅆxyx=116729x214x5+9A2yx36x5+A3+127518x10+1479Ax581A2yx6x5+A149x5+Ax,_ax=127x356x5+36Ayx54x69Ax6x5+A2

(22)

power_equivalent,keep

ⅆⅆxyx=169112514x+9A2yx3x26x+A3+11527A2+1959Ax1526x2yxx14x+9A6x+A+1,x=_a5,yx=454_a108_a10_b_a+6_a9+36A_a5_b_a+_a4A+3A2_b_a14_a5+9A

(23)

First order ODEs of Riccati type

As an example of Riccati type, consider

ⅆⅆxyx=yx2+x5x5

ⅆⅆxyx=yx2+x5x5

(24)

An equivalent Riccati equation in u(t), in normal form with minimal degrees, and the corresponding transformation are

power_equivalent,ut

ⅆⅆtut=ut21367t4t3,t=x3,ut=13yxx776x3

(25)

Changing variables back we reobtain the starting equation

PDEtools:-dchange2,1,normal

11821x4+2ⅆⅆxyxx14yxx10=11821x814yxx4+2x5+2yx2x14

(26)

normalisolate,ⅆⅆxyx

ⅆⅆxyx=yx2+x5x5

(27)

See Also

convert/Abel_RNF

convert/FirstKind

convert/NormalForm

convert/ODE

DEtools

DEtools[Lie]

dsolve

dsolve/hyper3

eval

factor

PDEtools[dchange]

simplify

 


Download Help Document

Was this information helpful?



Please add your Comment (Optional)
E-mail Address (Optional)
What is ? This question helps us to combat spam