陳艷秋，張臘娥

(湖南有色金屬職業(yè)技術(shù)學(xué)院，湖南 株洲 412006)

?

三角網(wǎng)格上Lagrange-Thiele型有理插值

陳艷秋，張臘娥

(湖南有色金屬職業(yè)技術(shù)學(xué)院，湖南 株洲 412006)

從Lagrange插值多項式出發(fā),結(jié)合Thiele型連分式,構(gòu)造了三角網(wǎng)格上Lagrange—Thiele型二元有理插值函數(shù),通過定義偏逆差商,建立遞推算法,構(gòu)造的插值函數(shù)滿足有理插值問題中所給的插值條件,并給出了插值的特征定理,最后給出的數(shù)值例子,驗證了所給算法的有效性。

三角網(wǎng)格；有理插值；遞推算法；特征定理

眾所周知，多項式插值結(jié)構(gòu)緊湊，思路清晰，運算簡單，在整個數(shù)軸上都有任意階導(dǎo)數(shù)，能夠進行函數(shù)值和微積分的運算。但是,利用多項式插值所求的插值多項式即使通過了所有給定的插值節(jié)點，但在其它點上的誤差卻可能很大，并且插值多項式的次數(shù)在七次以上時會出現(xiàn)Runge現(xiàn)象，導(dǎo)致數(shù)值不穩(wěn)定。為了克服多項式插值的缺點，引入了有理分式函數(shù)插值，即有理插值，有理插值比多項式插值復(fù)雜很多，因為有理插值不但計算更加復(fù)雜，而且可能不適定，但是有理插值卻為某些非線性問題提供了良好的解決途徑。連分式是構(gòu)造有理插值的常用方法之一，常見的有理插值的構(gòu)造都是假定有理插值問題有解的條件下給出的。對于矩形網(wǎng)格上的有理插值的研究，已經(jīng)吸引了許多作者的興趣，取得了豐碩的成果[1-6]，但考慮到某些給定的數(shù)據(jù)，如矩形網(wǎng)格是病態(tài)的，可能缺失部分?jǐn)?shù)據(jù)，此時將不可能在矩形網(wǎng)格上構(gòu)造一個有理插值函數(shù)[7-9]，這也是我們研究三角網(wǎng)格上的有理插值的原因。并且本文巧妙地將Lagrange多項式插值與Thiele型連分式插值結(jié)合起來，通過定義偏逆差商，建立遞推算法，構(gòu)造了二元Lagrange—Thiele型理插值函數(shù)，給出了插值的特征性質(zhì)和數(shù)值例子。還可以利用Samelson廣義逆將本文和插值方法推廣到向量值有理插值和矩陣值有理插值的情形。

所謂Lagrange插值，構(gòu)造多項式：

(1.1)

其中

(1.2)

設(shè)平面上的點集由下表給出：

(1.3)

其中xi≠xj,yi≠yj(i≠j)，稱之為左下三角網(wǎng)格，記為

(1.4)

1 三角網(wǎng)格上插值公式的構(gòu)造

(2.1)

其中

(2.2)

遞推算法定義如下：

(2.3)

(2.4)

則由(2.1)～(2.4)可得如下插值定理：

定理1 令aij=φ[xi;y0,y1,…yj]如(2.3)(2.4)式所定義，若所有的aij都存在且不為零，f(x,y)是定義在包含LB的區(qū)域D上的二元函數(shù)，則由(2.1)、(2.2)定義的R(x,y)為f(x,y)在LB上的Lagrange—Thiele型二元有理插值函數(shù)，且滿足：R(xi,yj)=f(xi,yj)，(xi,yj)∈LB。

R(xi,yj)=A0(yj)L0(xi)+A1(yj)L1(xi)+…+An(yj)Ln(xi)=Ai(yj)Li(xi)=Ai(yj)

=φ[xi;yj]=f(xi,yj)

故：R(xi,yj)=f(xi,yj)，(xi,yj)∈LB。定理得證。

2 特征定理

3 數(shù)值例子

例1 設(shè)f(x,y)在LB上的初始數(shù)據(jù)如下表：

表1

解 利用混合逆差商的遞推算法(2.3)～(2.4)式，列表計算如下：

表2 中間結(jié)果1

表3 中間結(jié)果2

表4 中間結(jié)果3

可驗證滿足插值條件，并且滿足特征定理的結(jié)論。

5 結(jié)束語

本文從Lagrange插值基函數(shù)出發(fā)，構(gòu)造了三角網(wǎng)格上的Lagrange—Thiele型二元有理插值函數(shù)，對于另外三種三角網(wǎng)格，我們也可以通過定義相應(yīng)的遞推算法，分別構(gòu)造其上的Lagrange—Thiele型二元有理插值函數(shù)。本文插值算法結(jié)構(gòu)簡單，使用方便，并且利用Samelson廣義逆可以將這種插值算法推廣到向量值有理插值和矩陣值有理插值情況。

[1] W.Semaszko.Thiele-type branched continued fractions for two variable functions[J].J.Comput.Appl.Math.9(1983)：137-153.

[2] A.Cuyt,Verdonk B.Multivaite rational interpolation[J].Computing,1985,34：141-161.

[3] Tan Jieqing，Tang Shuo. Composite schemes for multivariate blending rational interpolation[J].J ComputAppl Math,2002,144(1/2)：263-275

[4] 王家正.Stieltjes-Newton型有理插值[J].應(yīng)用數(shù)學(xué)與計算數(shù)學(xué)學(xué)報，2006,20(2):15-22

[5] Tan.J ,Q.Bivariate blending rational interpolants.Approx Theory & Its Appl,15(2) (1999) 74-83

[6] Jieqing Tan and Yi Fang.Newton-Thiele’s rational inter Polants，Numevical Algorithms.24(2000)：141-157.

[7] 王家正,梁艷.三角網(wǎng)格上的對稱型混合有理插值[J].西安工程大學(xué)學(xué)報,2008, 22(3): 354-357.

[8] Tan J Q,Song B R,Zhu G Q.Vector valued rational interpolants over triangular grid[J]Computers and Mathematics with Applications,2002，44(10211):1357-1367.

[9] 余小磊，唐爍.三角網(wǎng)格上新的插值公式的構(gòu)造[J].中國科學(xué)技術(shù)大學(xué)學(xué)報，2011，41(6)：504-511.

[10] 檀結(jié)慶.連分式理論及其應(yīng)用[M].北京：科學(xué)出版社，2007.

Lagrange-Thiele Type Rational Interpolation Over Triangular Grid

CHEN Yanqiu, ZHANG Lae

(HunanVocationalandTechnicalCollegeofNonferrousMetals,Zhuzhou412006,China)

In this paper, Lagrange-Thiele type bivariate rational interpolation has been constructed over triangular grid, which is based on Lagrange interpolating polynomial and combined with Thiele' continued fractions, by defining partial inverse differences, the recursive algorithms is given, The Lagrange- Thiele rational interpolating function is satisfied with the given interpolating conditions, the characterization theorems of the rational functions is obtained. At last, an example is given to illustrate the effectiveness of the interpolating algorithms.

Triangular grids; rational interpolation; recursive algorithms; characterization theorems

2017-02-06

湖南有色金屬職業(yè)技術(shù)學(xué)院院級項目

陳艷秋(1983-)，女，河南商水人，碩士，講師，研究方向:數(shù)值逼近。

O241.3

A

1674-2273(2017)03-0006-03