(完整版)几种插值法比较与应用

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

多种插值法比较与应用

(一)Lagrange 插值 1. Lagrange 插值基函数 n+1个n 次多项式

∏≠=--=

n

k

j j j k

j

k x x

x x x l 0)( n k ,,1,0 =

称为Lagrange 插值基函数 2. Lagrange 插值多项式

设给定n+1个互异点))(,(k k x f x ,n k ,,1,0 =,j i x x ≠,j i ≠,满足插值条件

)()(k k n x f x L =,n k ,,1,0 =

的n 次多项式

∏∏

∏=≠==--==n

k n

k

j j j

k j k k n

k k n x x x x x f x l x f x L 0

00

))(()()()(

为Lagrange 插值多项式,称

∏=+-+=-=n

j j x n n x x n f x L x f x E 0

)1()()!1()()()()(ξ 为插值余项,其中),()(b a x x ∈=ξξ (二)Newton 插值 1.差商的定义 )(x f 关于i x 的零阶差商

)(][i i x f x f = )(x f 关于i x ,j x 的一阶差商

i

j i j j i x x x f x f x x f --=

][][],[

依次类推,)(x f 关于i x ,1+i x ,……,k i x +的k 阶差商

i

k i k i i k i i k i i i x x x x f x x f x x x f --=

+-+++++]

,,[],,[],,,[111

2. Newton 插值多项式

设给定的n+1个互异点))(,(k k x f x ,n k ,,1,0 =,j i x x ≠,j i ≠, 称满足条件

)()(k k n x f x N =,n k ,,1,0 =

的n 次多项式

)()](,,,[)](,[][)(10100100---++-+=n n n x x x x x x x f x x x x f x f x N

为Newton 插值多项式,称

],[,)(],,,[)()()(010b a x x x x x x f x N x f x E n

j j n n ∈-=-=∏=

为插值余项。 (三)Hermite 插值

设],[)(1b a C x f ∈,已知互异点0x ,1x ,…,],[b a x n ∈及所对应的函数值为0f ,1f ,…,n f ,导数值为'0f ,'1f ,…,'n f ,则满足条件

n i f x H f x H i i n i i n ,,1,0,)(,)(''1212 ===++ 的12+n 次Hermite 插值多项式为

)()()(0

'12x f x f x H j n

j j j n

j

i n βα∏∏=++=

其中

)())((,)]()(21[)(2

2'x l x x x l x l x x x j j j j j j j j ---=βα

称为Hermite 插值基函数,)(x l j 是Lagrange 插值基函数,若

],[22b a C f n +∈,插值误差为

220)

22(12)()()!

22()

()()(n x n n x x x x n f

x H x f --+=

-++ ξ,),()(b a x x ∈=ξξ

(四)分段插值

设在区间],[b a 上给定n+1个插值节点 b x x x a n =<<<= 10

和相应的函数值0y ,1y ,…,n y ,求作一个插值函数)(x ϕ,具有性质

①i i y x =)(ϕ (n i ,,2,1,0 =)。

②)(x ϕ在每个小区间内],[1+i i x x (n i ,,2,1,0 =)上是线性函数。 (五)样条插值

设在区间],[b a 上取n+1个节点

b x x x a n =<<<= 10 给定这些点的函数值)(i i x f y =。

若函数)(x s 满足条件: ①i i y x s =)(,n i ,,2,1,0 =;

②在每个区间],[1+i i x x (n i ,,2,1,0 =)上是3次多项式; ③],[)(2b a C x s i ∈; ④取下列边界条件之一:

(ⅰ)第一边界条件:)()(0'0'x f x s =,)()(''n n x f x s =,

(ⅱ)第二边界条件:)()(0''0''x f x s =,)()(''''n n x f x s =或

0)()(''0''==n x s x s

(ⅲ)周期边界条件:)()(0n k k x s x s =,2,1=k

称)(x s 为3次样条插值函数。 (六)有理插值

设在区间],[b a 上给定n+m+1个互异节点

0x ,1x ,2x ,……,1-+m n x ,m n x + 上的函数值)(i i x f y =,m n i +=,,2,1,0 ,构造一个有理插值

m

m m m n

n n n m n mn b x b x b x b a x a x a x a x Q x p x R ++++++++=

=----11101110)()()( , 满足条件:

)()(i i mn x f x R =,m n i +=,,2,1,0

则称)(x R mn 为点集{0x ,1x ,2x ,……,1-+m n x ,m n x +}上的有理插值函数。

例1.设0x ,1x ,…,n x 为n+1个互异的插值节点,)(0x l ,)(1x l ,…,

)(x l n 为Lagrange 插值基函数,证明 ∏=≡n

j j x l 01)(

证 考虑1)(≡x f ,利用Lagrange 插值余项定理

)())(()!

1()

()()(101n n n x x x x x x n f x L x f ---+=-+ ξ

显然 1)()(≡=x f x L n 。

利用Lagrange 基函数插值公式,有

k j n

j k j j n

j j n x x l x x l x f x L =⋅==∏∏==)()()()(0

例2 给出下列表格:

相关文档
最新文档