OrderBasis - Maple Programming Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Linear Algebra : OrderBasis

OrderBasis

compute an order basis

 

Calling Sequence

Parameters

Description

Examples

References

Calling Sequence

OrderBasis([f1, f2, ..., fn], x, N, [d1, d2, ..., dn])

Parameters

f1, ..., fn

-

expressions; represent the functions to be approximated

x

-

variable appearing in the fis

N

-

(optional) non-negative integer; specify the order of approximation. You must specify at least one of N or d1, ..., dn.

d1, ..., dn

-

(optional) non-negative integers; specify the degree bounds. You must specify at least one of N or d1, ..., dn.

Description

• 

The OrderBasis([f1, ..., fn], x, N, [d1, ..., dn]) command computes an order basis for the functions f1, ..., fn with respect to the variable x, the degrees d1, ..., dn, and the order N. It finds all polynomial coefficients that provide an identity of the form

f1xv1x++fnxvnx=0

  

up to a certain number of terms and with the degree of each vi bounded. This is similar to what the IntegerRelations[PSLQ] algorithm does for finding integer relations for floating-point numbers.

• 

More precisely, given functions fi assumed to have a series expansion about 0, the OrderBasis function returns a matrix whose columns provide a basis for the (mathematical) module defined by

L={v1x,...,vnx|f1xv1x+...+fnxvnx=OxN,degreev1d1,...,degreevndn}

• 

That is, for every vector v of polynomials in L there exist n polynomials a1,a2,,an 

v=a1xcolumnM,1+a2xcolumnM,2++anxcolumnM,n.

  

Here the degree of ai is bounded by didegreeMi,i.

• 

If there are three arguments with the third argument as a positive integer N, that is, OrderBasis([f1, ..., fn], x, N), then the degree bounds d1, ..., dn are assumed to be N, ..., N.

• 

If there are three arguments with the third argument as a list, that is, OrderBasis([f1, ..., fn], x, [d1, ..., dn]), then the order is determined by dn1++dnn+n1.

Examples

F1+x2x7+x12,sinx,ⅇx:

MOrderBasisF,x,8,3,1,2

M:=x4+682383604251353x3+5119423444251353x2+9717291364251353x26699764251353201169912754059x3231284251353x242526544251353x+346924251353553379912754059x3+170494524251353x2+337088984251353x6606042513538298194044251353+12936111604251353xx21416598502706x+84214698502706687851458502706+937995118502706x2669976425135318042185164251353x346924251353+1532238502706xx21363350618502706x+660604251353

(1)

Each column of M has order 8.

mapseries,mapexpand,Matrix1,3,F.M,x,8

Ox8Ox8Ox8

(2)

mapdegree,M,x

433121112

(3)

Implies that a basis for all v1,v2,v3 satisfying v1F1+v2F2+v3F3 with degreev13,degreev21,degreev32 is a constant * column(M,3).

In the next example, OrderBasis(F,x,8) is the same as OrderBasis(F,x,8,[8,8,8]).

MOrderBasisF,x,8

M:=x3693842011699x2127579622011699x+1040762011699405437882011699x21052661122011699x+464712201169980977402011699x2+214862162011699x764162011699127540592011699x24249774023398x+252644074023398x3+324253202011699x2954986402011699x30079248201169955337992011699x2+222847052011699x+1079330420116991040762011699+4596694023398x4647122011699+1358100722011699xx2323559362011699x+764162011699

(4)

Each column of M has order 8.

mapseries,mapexpand,Matrix1,3,F.M,x,8

Ox8Ox8Ox8

(5)

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.

  

Beckermann, B. and Labahn, G. "A Uniform Approach for the Fast Computation of Matrix-Type Pade Approximants." SIAM Journal on Matrix Analysis and Applications. Vol. 15 No. 3. (1994): 804-823.

See Also

IntegerRelations[PSLQ]

numapprox[pade]

series

 


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