矩陣課件-第5章-5.2.4hermite與分段低次2017春

ixoxPnxi

f

xi

Pnxi

f

xi

L

L

L

Li

若有相同的切線若彎曲方向相同近似程度越來越好Pn

xi

f

xi

,

i

0

,

1

,

L

,

nyf

(

x)

Pn

(

x)Pi

y

f

(

x)i在x

處相交Lagrange插值多項滿in次多項式pn(x)，使其滿足插值條件：p

(

i

)

(x

)

f(i

)

(x

)

y(i

)

,其中1

,

2

,L

,

s(5-18)5.2.4

Hermite插值理論和應(yīng)用中

某些插值問題，要求插值函數(shù)p(x)具有一定的光滑度，即在插值節(jié)點處滿足一定的導(dǎo)數(shù)條件，這類插值問題稱為Hermite插值問題。Hermite插值問題的一般提法是：設(shè)已知函數(shù)

f(x)在s

個互異點

x1,x2

,…,xs

處的函數(shù)值和導(dǎo)數(shù)值：為正整數(shù)，記1

2

L

s

n

1,

構(gòu)造一個i

1,

2,

L

,

s

;i

0,1,L,i

1。1

1f

x

,

f

x

,

L

,f

(1

1)

x

;2

2

f

x

,

f

x

,

L

,

fx

1;22(

1)

ss

sf

x

,

f

x

,

L

,

fsx

,(

1)L

LL這樣將得到如下形式的n次插值多項可以采用類似于構(gòu)造Lagrange插值基函數(shù)li(x)的方法來解Hermite插值問

構(gòu)造一批

n

次多項式作為插值基函數(shù)，Li,k

(x),

i

1,2,

L,

s

;

k

0,1,L,i

1,li

xixi

1

,

2

,

L

,

sif

x

iif

x

i

1sx

f

xis

lip

x

i1if

x

f

x

Li

1

,

2

,

L

,

six回p

xn？si1LLi

i,0

i

i,1y

x

y

xyiii,1

L

x

1

L

i

1

,

2

,

L,

ssk

x

yinp

x

i1

k

0且滿足插值條件（5-18）L

i,ki

1

yi

i

1

,

2

,

L,

s

,i

pi

i

0,1,L,i

1。滿問二點三次Hermite插值多項設(shè)已知函數(shù)f(x)在2

個互異點x1,x2

處的函數(shù)值和導(dǎo)數(shù)值：f

x1

,

f

x1

;

f

x2

,

f

x2

;構(gòu)造一個三次多項式H3

x

ax

bx32

cx

d使其滿足插值條件H3

x1

f

x1

,

H3

x1

f

x1

;H3

x2

f

x2

,

H3

x2

f

x2

;1

1,0

2

12,0

1,12(x)

f(x

)

L

(x)

f

(x

)

L

(x)

f

(x

)

L2,1(x)L1,0

(x1

)

1L1,1(x1

)

1L1,0

(x2

)

L1,0(x1

)

L1,0(x2

)

0L1,1(x1

)

L1,1(x2

)

L1,1(x2

)

0L2,0

(x2)

1

L2,0

(x1

)

L2,0

(x1

)

L2,0

(x2

)

0L2,1(x2

)

1

L2,1(x1

)

L2,1(x1

)

L2,1(x2

)

0以L1,0(x)為例計算之,L2,0(x),L1,1(x),L2,1(x)同例1用基函數(shù)法來構(gòu)造三次多項式H3(x)。H3

(x)

f

(x

)

L其中L1,0(x),L2,0(x),L1,1(x),L2,1(x),為插值基函它們滿足設(shè)21,0

2x

ax

bL

(x)

x

21,0

22L

(x)

2

x

xax

b

a

x

x由于L1,0(x)為三次多項式,又L1,0

(x2

)L1,0(x2

)

0,

故應(yīng)有又由L1,0

(x1

)

1,L1,0(x1

)0,

進一步21

2

1x

x

ax

b

12121

ax

b

1

2112

ax

b

a

x

x

0

31

2x

x

2x

x

a

231x1

x2

3x

xb

1,0L2

1 2

1

2

(x)

2

1

1

x

x

2x

x

x

xx

x

代入上式那么，有

x

x

2L1,1(x)

2

x

x1

,2,0Lx2

x1

(x)

1

1

2

2

x2

x1

x

x x

x2

x2

x1

x1

x2

x

x

2L2,1

(x)

1

x

x2

1212x

x

x

xx

xx

x1 2

1 2

2

x2

x1

x

x

f

(x1

)(x

x1

)22122x

x

x

x

x

x

x

x

f

(x

2

)1

2

1

2 1

12

2

x

x

2

f

(x

)(x

x

)

x

x

2 1

從而得3次插值多項式:3H

(x)

f

(x

)

1

2這就是二點三次Hermite插值多項式,其滿足插值條H3

xk

f

xk

H3

xk

f

xk

k

1,

24 2

x

,

,構(gòu)造三次Hermite多項式練習(xí)已知f(x)=sinx,2

,2f

sin

4

4

2

,2f

cos

4

4

f

sin

1

,

2

2f

cos

0

,

2

2

則

2

x

x

4

2

f1

2

2

x

f

x

2

4

4

4 2

4

4 2

4 2

2x

x

2

4

f

21

2

解

已知3H

x

2228x

4x

2

3

22

2 4

2 4

2

x

4x

2

2

4

2

3

4x

4x

32

3f

(x)

p

(x)

min(x1

,x2

)

max(x1

,x2

)。其中設(shè)f(x)∈C3[a,b]，在（a,b）內(nèi)4階可

22124!(4)f

()x

x x

x

,

x

[a,

b]定理5.3’又設(shè)a≤x1＜x2≤b則兩點三次Hermite插值多項式p3(x)有如下的誤差估計

3f

x

p

x

2x

,

x

K對上述給定的

x

,引進輔助函數(shù):

2212x

xx

x

K

x

3

t

f

t

p

t

2212K

x t

xt

x

,證

若x為

x1,x2中的某一個,

則誤差估計式顯然成以下假設(shè)

x≠xi

(i=1,2),

由插值條件,可3f

(x)

p

(x)

min(x1

,x2

)

max(x1

,x2

)。其中設(shè)f(x)∈C3[a,b]，在（a,b）內(nèi)4階可

221

24!(4)f

()x

x x

x,

x

[a,

b]定理5.3’又設(shè)a≤x1＜x2≤b則兩點三次Hermite插值多項式p3(x)有如下的誤差估計4!f

(4)

K

x

于是，代入估計式即知結(jié)論成立。顯然

xi

xi

0,

(x)

0

。i

1,

2

,使得(

4

)

()

f

(4)

()

0

K

x

4!

0

,2231

2(t)

f

(t)

p

(t)

K

(x)(t

x

) (t

x

)

t

有2個二重零點

x1,x2

和一個單重零點x

。反復(fù)運用Rolle定理可證，至少有一個ξ,且min

x,

x1,

x2

max

x,

x1,

x2

5.2.5

分段低次插值利用插值法構(gòu)造近似函數(shù)時，為了提高近精度，經(jīng)常需要增加插值節(jié)點，加密插值節(jié)點會使插值函數(shù)與

值函數(shù)在

節(jié)點上的取值相同，那么誤差是否會隨之減小呢？答案是否定的。原因在于插值節(jié)點增多導(dǎo)致插值多項式的次數(shù)增高，而高次多項式的振蕩次數(shù)增多有可能使插值多項式在非節(jié)點處的誤差變得很大。51

x2nkx

5

10

k在[-5,5]上構(gòu)造等距節(jié)k

0,1,L

,n。例如，對于函數(shù)

f

x

分別取

n=6、n=8

和

n=10作出插值多項式pn(x)

近

5

4

3

2

10

5

1

x2f

x

p8

xp6

x等距節(jié)點高次插值多項式的Rung現(xiàn)象yp10

xaxben

max

In

f

,

x

In

f

,

xni0nni0,

I

f

,

xli

x

fi

。nnil

x

,

I

maxa

xb

i0

插值函數(shù)的穩(wěn)定性的分析，得到插值函數(shù)的舍入誤差項為：其顯然，對等距節(jié)點的高次的Largrange多項式插值ηn是隨著n增長Runge現(xiàn)象對等距節(jié)點的高次插值多項式的是典型nmaxaxb

i0

1in

li

xmax

fi

fini1in

max

f

fi故得出結(jié)論h(5-26)為了克服高次插值多項式的上述弊端，通常采用分段低次

插值的方法，即以插值節(jié)點為分點，將[a,b]分成若干個小區(qū)間，并在每個小區(qū)間上進行低次的多項式插值。一、分段線性Lagrange插設(shè)插值節(jié)點

x0,x1,…,xn滿足a≤x0<x1<…<xn≤b,在每一個間[xk,xk+1](k=0,1,…,n-1)上做線性插值多項式x

[xk

,xk

1]。Lh

(x)

hL(0)

(x),hL(1)

(x),x

[x0

,

x1],x

[x1

,

x2

],Mx

[xn1

,

xn

],hML(n1)

(x),令(5-27)L(k

)

(x)

ykkk

1x

xk

yk

1

xk

1x

xk

1

x

xk

x

,y=Lh(x)的圖形是平面上連接點1(x)(k

)hR

(x)

f

(x)

L2k

1k (x

x

)(x

x

),212Mk

k1axbmax

|

R

(x)

|M

2

max

|

f

(x)

|,

h

max

hk

,a

xb

0k

n1從而其中顯然

Lh(xi)=yi

(i=0,1,…,n),

Lh(x)稱為

f(x)在[a,b]上的分線性插值多項式(5-29)(5-28)8hk

xk

1

xk|

(x

x

)(x

x

)

|

M

2

h2

,

5

4

3

2

10(x1,y1)、…、(xn,yn)的一條折線（如圖）。y由插值余項定理，當(dāng)f(x)在[a,b]上二次可微時，對任意x∈[xk

,xk+1]，余項

f

()

5

1

x2f

x

hL(k

)

(x)hL(n1)

(x)h0n作為f(x)的近似值。k

k+1對x∈[a,b],若x∈[x

,x

],則h0易證,當(dāng)f(x)∈C[a,b]時,lim

Lh

(x)

f

(x)

在[a,b]上一致成立L(k

)(x)

作為f(x)的近似值h,則以

L(0)

(x)

作為f(x)的近似值

若x≥x

,則若x≤x二、分段二次Lagrange插值1當(dāng)給定的函數(shù)表

點的個數(shù)遠多于3的時候,為了提高計精度,或根據(jù)實際問題需要,有時采取分段二次插值法對于x∈[a,b],應(yīng)選擇靠近x的三個節(jié)點做二次插值多項式當(dāng)x∈[xk,xk+1],且x偏向xk時,選擇xk-1,xk,xk+1作為插值節(jié)點；當(dāng)x∈[xk,xk+1],且x偏向xk+1時,選擇xk,xk+1,xk+2作為插值節(jié)點；234當(dāng)x∈[x0,x1),或x＜x0時,選擇x0,x1,x2作為插值節(jié)點當(dāng)x∈(xn-1,xn],或x＞xn時,選擇xn-2,xn-1,xn作為插值節(jié)點根據(jù)實際問題的需要,還可采用分段