趙前進，侯中麗

（安徽理工大學(xué)理學(xué)院，安徽 淮南232001）

當插值節(jié)點數(shù)較大時，Thiele型連分式有理插值可能比多項式插值的逼近效果更好。然而，有理插值函數(shù)難以避免在插值區(qū)間內(nèi)出現(xiàn)極點，也難以控制極點的位置，另外還可能有不可達點。重心有理插值比Thiele型連分式有理插值計算量小，數(shù)值穩(wěn)定性好，選擇適當?shù)臋?quán)可以不出現(xiàn)極點和不可達點。Berrut，Schneider，Nguyen等對重心有理插值進行了深入的研究［4-13］。在文獻［1］中，F(xiàn)loater和 Hormann通過在子節(jié)點集上構(gòu)造插值多項式，然后用特定的權(quán)函數(shù)對這些插值多項式進行重心型的混合，構(gòu)造了一種無極點、高精度的復(fù)合重心有理插值。在文獻［2］中，Klein就等距節(jié)點的情形又對此進行了改進，通過修正子節(jié)點集上的插值多項式構(gòu)造了一種新的復(fù)合重心有理插值。

本文將文獻［1］中的方法推廣到矩形域上的二元復(fù)合重心型混合有理插值。首先在小矩形域上構(gòu)造二元Newton插值多項式［3］，然后基于特定的權(quán)函數(shù)進行重心型的復(fù)合，構(gòu)造了二元復(fù)合重心型混合有理插值，并進一步分析了新的二元復(fù)合重心型混合有理插值的一些性質(zhì)，如無極點和不可達點等，最后由給出的數(shù)值例子驗證了新方法的有效性。

1 二元復(fù)合重心型混合有理插值

設(shè)f（x，y）在D 上有定義，且記f（xi，yj）＝fij，（i＝0，1，2，…，m；j＝0，1，2，…，n）。

對任意整數(shù)d1和d2（0≤d1≤m，0≤d2≤n），對于每個i＝0，1，2，…，m-d1，j＝0，1，2，…，n-d2，設(shè)Pij（x，y）為｛（xk，yq）｜k＝i，i＋1，…，i＋d1；q＝j(luò)，j＋1，…，j＋d2｝上的二元 Newton插值多項式，基于重心型復(fù)合，構(gòu)造二元有理函數(shù)

2 插值性質(zhì)

性質(zhì)2．1 二元有理函數(shù)R（x，y）在矩形域D內(nèi)沒有極點。

證明 只需證明R（x，y）的分母大于零?，F(xiàn)在對（3）式的分子、分母同時乘以（-1）n-d2（y-y0）（y-y1）…（y-yn）有

3 數(shù)值例子

例1 取被插值函數(shù)f（x，y）＝（x＋y）ln（x2＋y2＋1），在［-5，5］×［-5，5］上取等距節(jié)點，將m＝10，n＝10，d1＝5，d2＝5及m＝20，n＝20，d1＝5，d2＝5的插值函數(shù)和被插值函數(shù)用MATLAB繪制如下圖，并求出了最大絕對誤差（f（x，y）＝（x＋y）ln（x2＋y2＋1）與插值函數(shù)在插值區(qū)間上的最大誤差）。

圖1 被插值函數(shù)

圖2 m＝10，n＝10插值函數(shù)

圖3 m＝10，n＝10誤差函數(shù)

圖4 m＝20，n＝20插值函數(shù)

圖5 m＝20，n＝20誤差函數(shù)

表1 誤差比較（d1＝5，d2＝5）

由上表可見，m、n越大，插值誤差越小。

4 結(jié)論

筆者給出矩形域上的二元復(fù)合重心型混合有理插值新方法，首先在小矩形域上構(gòu)造二元Newton插值多項式，然后基于重心型復(fù)合，構(gòu)造出了二元復(fù)合重心型混合有理插值，證明了二元復(fù)合重心型混合有理插值無極點和不可達點，最后用數(shù)值例子驗證了新方法的有效性。

［1］Michael S．Floater，KaiHormann．Barycentric Rational Interpolation with No Poles and High Rates of Approximation．Numer．Math．，2007，107：315-331．

［2］Georges Klein，An Extension of the Floater-Hormann Family of Barycentric Rational Interpolants．Mathematics of Computation，2013，82（284），2273-2292．

［3］Gasca M，Sauer T．On the History of Multivariate Polynomial Interpolation［J］．Comput．Appl．Math．，2000，122：23-35．

［4］Berrut J．-P．，Trefethen L N．，Barycentric Lagrange Interpolation．SIAM．Rev．，2004，46：501-517．

［5］Schneider C．，Werner W．，Some New Aspects of Rational Interpolation．Math．Comp．，1986，175（47）：285-299．

［6］Berrut J．-P．Mittelmann H．Matrices for the Direct Determination of the Barycentric Weights of Rational Interpolation．J．Comput．Appl．Math．，1997，78：355-370．

［7］Jesus M．Carnicer，Weighted Interpolation for Equidistant Nodes，Numerical Algorithms，2010，55（2-3）：223-232．

［8］Schneider C．，Werner W．Hermite Interpolation：The Barycentric Approach．Computing．，1991，46：35-51．

［9］Berrut J．P．，Rational Functions for Guaranteed and Experimentally Well-conditioned Global Interpolation．Computers ＆ Mathematics with Applications，1988，15（1）：1-16．

［10］L Knockaert．A Simple and Accurate Algorithm for Barycentric Rational Interpolation．Signal Processing Letters［DB／OL］．IEEE，2008（15）：154-157．

［11］Berrut J．P．A Matrix for Determining Lower Complexity Barycentric Representations of Rational Interpolants．Numerical Algorithms，2000，24（1-2）：17-29．

［12］Higham N．J．The Numerical Stability of Barycentric Lagrange Interpolation ［J］．IMA J．Numer．Anal．，2004，24（4）：547-556．

［13］H．T．Nguyen，A．Cuyt，O．S．Celis．Shape Control in Multivariate Barycentric Rational Interpolation ［J］．Proc．ICNAAM，2010，1281：543-548．