Font Size: a A A

Some Research On Rational Interpolation

Posted on:2010-08-23Degree:DoctorType:Dissertation
Country:ChinaCandidate:S T ChenFull Text:PDF
GTID:1100360302965945Subject:Computational Mathematics
Abstract/Summary:PDF Full Text Request
As a generalization of polynomial interpolation, rational interpolation has many advantages. In this paper we study matrix valued rational interpolation problem and Birkhoff-type rational interpolation.Ⅰ. Matrix Valued Rational InterpolationIn this paper we convert the task of seeking the weak solution (P(X), q(X)) of matrix valued osculatory rational interpolation into computing the Groebner bases of R-submodule M of the free module over the polynomial ring. This Groebner bases can be computed recursively. Then we can obtain a parametric rational interpolation function with the Groebner bases. By choosing these parameters properly, we may get the desired matrix valued rational interpolation function.A subsetΑ(?) Nn is called delta set, if it closed under the divison order; that is, ifα∈Αthenβ∈Αfor allα= (α1,... ,αn),β= (β1,...,βn) componentwise.Let D a differential operator.The problem of osculatory rational interpolation can be stated as follows:Definition 1. Given a set of L distinct points in space Rn {Y1,...,YL}. Each point have multiplicity defined by the sets {Α1,... ,ΑL}. And corresponding function values (?). Construct a rational interpolation functions.tFrom the definition of matrix valued oscillatory rational interpolation we know that the equivalent definition is: the Taylor series expansion of R(X) about the points X = Yi satisfiesDefinition 2.For each point Yi and the lower setΑi, define polynomial matrix Hi(X)that isDefinition 3.(weak interpolations) (P,q)∈(Rd1×d2[X],R[X]) is called a weak interpolation. ifWhere (?) is the vanishing ideal ofΑi.Let M = {(P,q)|P≡qH modⅠ} (?) (Rd1×d2[X],R[X]) where H∈Rd1×d2[X] is a polynomial matrix, I (?) R[X] is an ideal. We defined (P1,q1) + (P2,q2) = (P1 + P2,q1+q2),g·(P1,q1) = (gP1,gq1). Thus M is a R-submoduleof (Rd1×d2[X],R[X]).Definition 4. (order (?)) 1. We say (?) if |α| > |β| or |α| = |β| and k2 < k1 , or |α| = |β| , k1 = k2 and l2 < l1 or |α| = |β| , k1= k2, l2 = l1 and (?);2. We say (?) if |α| > |β| +ξ;3. (?) if |α| > |β| or |α| = |β| and (?).where (?) is the lexicographic order on R[X], andξis a given integer.We can check that the order (?) is a monomial order on Rd+1.Definition 5.Let (?) = {{Y1,Α1,H1),..., (YL,ΑL,HL)}, We defined(?) is called matrix valued interpolation set, M(?) is called the matrix valued interpolation modules of (?).Proposition 1. Fix an order (?). Let M (?) (Rd1×d2[X],R[X]) be the submodule. We will denote by〈LT(M)〉the monomial submodule generated by the leading terms of all (?)∈M with respect to (?), then1. for any (?)∈(Rd1×d2[X],R[X]), there exists a unique element (?)∈(Rd1×d2[X],R[X]), s.t. (?) -(?)∈M, and (?) is a R-linear combination of monomials in (?).2. if (?) is finite-dimensional, thentherefore, dim (?). Where (?). Theorem 2. Fox an order on (Rd1×d2[X],R[X). Let W ={(Y1,Α1,H1), ...,(YL,ΑL,HL)} be the matrix valued interpolation set. Let V ={(Y1,Α1), ...,(YL,ΑL)}, I(V) be the vanishing ideal, thenProposition 3. M(?) is the matrix valued interpolation modules of (?). if (?), and (?), then (?) is a Groebner bases for M(?).Let M0=(Rd1×d2[X],R[X]),(?),m=1.....N.Where (?) is the vanishing ideal ofΑi,Hi(X) is polynomial matrix.Obviously M = MN and (?).We compute the Groebner bases of MN recursively. It is easy to see that (?)is the Groebner bases of M0. Using M0, if we can compute the Groebner bases of M1, the Groebner bases of M2 recursively. We can compute the Groebner bases of MN.If we have computed the Groebner bases of (?).Thuswhere (?) are the parameters.Ⅱ. Birkhoff-type rational interpolationThe problem of Hermite-Birkhoff rational interpolation.
Keywords/Search Tags:Interpolation
PDF Full Text Request
Related items