郭佳欣

（北京建筑大學(xué)，北京 102616）

實(shí)際問題中經(jīng)常有這樣的函數(shù)y=f(x)，其在某個(gè)區(qū)間[a，b]上有有限個(gè)離散點(diǎn)x0,x1,…xn，且這些點(diǎn)對應(yīng)函數(shù)值為yi=f(xi)(i=0,1,2...)，若想得到其它點(diǎn)的值就必須找一個(gè)滿足上述條件的函數(shù)表達(dá)式。

1 Lagrange 插值函數(shù)

1.求作n 次多項(xiàng)式pn(x)，使?jié)M足條件：

這就是所謂的拉格朗日（Lagrange）插值。點(diǎn)xi（它們互不相同）稱為插值節(jié)點(diǎn)。用幾何語言來表達(dá)這類差值，就是通過曲線y=f(x)上給定的n+1 個(gè)點(diǎn)(xi,yi)(i=0,1,...,n)，求作一條n 次代數(shù)曲線y=pn(x)作為y=f(x)的近似。

2.拉格朗日插值公式。

（1）首先考察線性插值的簡單情形。若y=f(x)表示過兩點(diǎn)(x0,y0),(x1,y1)的直線，這個(gè)問題是我們所熟悉的，它的解可表示為對稱式：

此類為一次插值，稱為線性插值。若令：

由此可得：l0(x0)=1,l1=(x1)=1,l1(x0)=0，則有：p1(x)=l0(x)y0+l1(x)y1，這里的l0(x),l1(x)可以看作是滿足條件l0(x0)=1,l0(x1)=0l1(x1)=1,l1(x0)=0的插值多項(xiàng)式，這兩個(gè)特殊的插值多項(xiàng)

式稱作上述問題的插值基函數(shù)[1-2]。

（2）拋物插值。線性插值僅僅利用了兩個(gè)節(jié)點(diǎn)的信息，精度自然很低，為了提高精度，進(jìn)一步考察二次插值。若y=f(x)表示過三點(diǎn)(x0,y0),(x1,y1),(x2,y2)的曲線，求它的解。二次插值的幾何解釋是，用通過三點(diǎn)(x0,y0),(x1,y1),(x2,y2)的拋物線y=p2(x)來近似所考察的曲線y=f(x)，因此這類插值亦稱為拋物插值。

根據(jù)上述解決方法，同理可以得出：

容易看出，這樣構(gòu)造的p2(x)為問題的解。

（3）推廣到一般：已知函數(shù)在n+1 個(gè)不同點(diǎn)x0,x1,...xn上的函數(shù)值分別為y0,y1,...yn求一個(gè)次數(shù)不超過n 的多項(xiàng)式pn(x)，使其滿足：pn(xi)=yi(i=0,1,...n)即n+1 個(gè)不同的點(diǎn)可以決定的一個(gè)n 次多項(xiàng)式。過n+1 個(gè)不同的點(diǎn)分別決定n+1 個(gè)n 次插值基函數(shù)。l0(x),l1(x),...ln(x)每個(gè)插值基多項(xiàng)式滿足：a.li(x)是n 次多項(xiàng)式；b.li(x)=1，而在其它n 個(gè)點(diǎn)li(xk)=0,(k≠i)。由于li(xk)=0,(k≠i)，故有因子(x-x0)...(x-xi-1)(x-xi+1)...(x-xn)因其已經(jīng)是n 次多項(xiàng)式，故而僅相差一個(gè)常數(shù)因子。令li(x)=a(x-x0)...(x-xi-1)(x-xi+1)...(x-xn)由li(xi)=1，可以定出a,進(jìn)而得到：

則n 次拉格朗日型插值多項(xiàng)式pn(x)=y0l0(x)+y1l1(x)+...ynln(x)，從而pn(x)是一個(gè)次數(shù)不超過n 的多項(xiàng)式，且滿足pn(xi)=yi(i=0,1,...n)[3-5]。

2 Newton 插值函數(shù)

2.1 Newton 法的簡述

設(shè)(x-xk)是f(x)的一個(gè)近似根，把f(x)在xk處泰勒展開：

若取前兩項(xiàng)來近似代替f(x)，則f(x)=0 的近似方程為f(x)=f(xk)+f'(x)(x-xk)=0。設(shè)f'(x)≠0，設(shè)其根為xk+1，則xk+1的計(jì)算公式為：

這即為牛頓法，其迭代函數(shù)為：

牛頓法其實(shí)是一種逐步的線性化的方法，它的本質(zhì)是將f(x)=0 非線性方程的求根問題歸結(jié)為類似的一系列的線性方程f(xk)+f'(x)(x-xk)=0 的根。

2.2 具有承襲性的插值公式

先考察線性插值的插值公式：

由于p0(x)=f(x0)可看作是零次插值多項(xiàng)式，上式表明p1(x)=p0(x)+c1(x-x0)其中，修正項(xiàng)的系數(shù)：

再修正p1(x)以進(jìn)一步得到拋物插值公式的解p2(x),為此，令：

顯然，不管系數(shù)c2 如何取值，p2(x)均能滿足兩個(gè)條件：p2(x0)=f(x0),p2(x1)=f(x1)，再用剩下的一個(gè)條件p2(x2)=f(x2)來確定c2，結(jié)果有：

記c0=f(x0)，從而有：p2(x)=c0+c1(x-x0)+c2(x-x0)(x-x1)

以上論述表明，為了建立具有承襲性的插值公式，需要引進(jìn)差商并研究其性質(zhì)。

2.3 差商及其性質(zhì)

對于給定的函數(shù)f(x),記f(x0,x1,...,xn)表示關(guān)于節(jié)點(diǎn)x0,x1,...,xn的n 階差商。一階差商定義為：

二階差商定義為一階差商的差商：

一般地，n 階差商遞推定義為：

為統(tǒng)一起見，補(bǔ)充定義函數(shù)f(x1)為零階差商。

一般地，用數(shù)學(xué)歸納法易證：

2.4 差商形式的插值形式

這種差商形式的插值公式稱作牛頓插值公式，牛頓公式其實(shí)只是拉格朗日公式的一種變形[6-7]。

3 Hermite 插值法

1.在某些問題中，為了保證插值函數(shù)能更好與原來的函數(shù)重合，使插值函數(shù)更接近原來函數(shù)軌跡，不但要求“過點(diǎn)”，即兩者在節(jié)點(diǎn)上具有相同的函數(shù)值，而且要求“相切”，即在節(jié)點(diǎn)上還具有相同的導(dǎo)數(shù)值，這類插值稱作切觸插值，或稱埃爾米特（Hermite）插值。顯然，埃爾米特插值是泰勒插值和拉格朗日插值的綜合和推廣。

2.Hermite 插值基本原理。通常如上條件的Hermite型插值是通過構(gòu)造相應(yīng)的插值基函數(shù)來完成的，為方便起見，以n=1 為例，說明傳統(tǒng)的求解方法，設(shè)給定的x0,x1和相應(yīng)的函數(shù)值f(x0),f(x1)及微商值f'(x0),f'(x1)構(gòu)造插值函數(shù)H3(x)。由構(gòu)造函數(shù)的辦法可知：對應(yīng)于x0和x1點(diǎn)函數(shù)值的插值函數(shù)分別為：

而對應(yīng)的x0和x1點(diǎn)導(dǎo)數(shù)值的插值基函數(shù)分別為：

因此所要求的插值函數(shù)H3(x)=f(x0)h0(x)+f(x1)h1(x)+f'(x0)H0(x)+f'(x1)H1(x)。

4 插值法的應(yīng)用

加工零件的輪廓曲線：該問題是有一個(gè)待加工零件，它的外形根據(jù)工藝要求由一組數(shù)據(jù)（x,y）給出（在平面情況下），要求要用數(shù)控機(jī)床加工時(shí)刀具必須沿這些數(shù)據(jù)點(diǎn)前進(jìn)，并且由于刀具每次只能沿x 方向或y方向走非常小的一步，所以需要將已知數(shù)據(jù)加密，得到加工所要求的步長很小的（x,y）坐標(biāo)。

表1 給出了輪廓線上x 每間隔0.2（長度單位）的加工坐標(biāo)x,y（順時(shí)針方向?yàn)樾?，由輪廓線的左右對稱性，表中只給出右半部的數(shù)據(jù)），假設(shè)需要得到x 或y坐標(biāo)每改變0.05 時(shí)的坐標(biāo)，試完成加工所需的加密數(shù)據(jù)，畫出該零件的輪廓曲線。運(yùn)用matlab 繪畫出的加工零件輪廓線，用不同的插值方法，畫出最好的效果圖。最后得出運(yùn)用三次樣條插值，畫出的輪廓圖最好。

表1 x 間隔0.2 的加工坐標(biāo)（x,y）（右半部的數(shù)據(jù)）

5 結(jié)語

我們通過對函數(shù)插值的分析及研究，用解決實(shí)際問題的方法，對不同的插值方法都有了深刻的了解，并掌握了運(yùn)用matlab 解決插值問題，明白了各個(gè)插值方法的含義。這讓我們不僅僅局限于理論知識(shí)的理解，更上升到了實(shí)際問題中，加深了我們對于插值函數(shù)的認(rèn)識(shí)。