数值分析实验一 插值与拟合3.docxVIP

PAGE

PAGE1

数值分析实验一 插值与拟合

组号 班级 学号 姓名 分数

一：实验目的

1、掌握三次样条插值的计算方法；

2、会用分段线性插值法求插值；

3、会用最小二乘法求拟合曲线；

4、掌握计算方法中插值与拟合的区别与联系；

5、熟悉Matlab软件中有关的命令，用Matlab做插值计算。

二：实验内容及基本知识介绍

1、三次样条插值函数：

）、若函数 S（x）∈c2[a,b],且在每个小区间[

x ,x ]上是三次多项式，其中

j j?1

a?x

0

?x? ? x

…1 n

…

?b是给定节点，则称S(x)是节点x

0

,x ,..., x

1 n

上的三次样条

函数。若在节点 x

j

上给定函数值 y

j

?f(x

j

)(?j

0,1n,并.成立

s(x

j

)?y

j

j(? 0 ,n1则,称S(x)为三次样条插值函数。

、三次样条插值的计算方法：

①因为在每个小区间上S

i

以写出它的表达式

(x)是三次多项式，所以S(x)在每个小区间上是直线，可

i

x ? x x? x

S(x)? m i?1 ?m i ,

i i x

x

i?1 i

i?1

x ? x

i?1 i

?1其中m ,

?1

i i

②把它积分两次，得到

(x ?x)3 (x?x)3

S (x)?m

i i

i?1 ?m

6h

i

i?1

i ?cx?d,6h

i

这里的c和d是积分常数，h

i

?x ?x

。i?1 i

。

利用S(x)

i i

?y和S(x )?y 可以确定c,d,于是有

i i i?1 i?1

(x ?x)3 (x?x)3

S(x)?m

i i

i?1

6h

i

m i i?1 6h

i

mh2 x

x m h2 x?x

(y

???i i) i?1

i 6 h

i

?(y ???i?1 i ) i,

i?1 6 h

i

将其求导数得到

(x

?x)2 (x?x)2

S(x)??m

i i

i?1 ?m

2h

i

i i?1 2h

i

y ?y m ?m

???i?1 i

h

i

???i?1 ih.

6 i

至此，我们把S (x)以及它的一、二阶导函数都用两个参数表示出来。

i

?1 ?1 ?1③我们令 S(x )?S (x ),i?0,1,...n?2, 得到一

?1 ?1 ?1

i i i i

0 1m ,m ,...,m

0 1

n

?1 ?1?m ?2m ??m ?d,i?1,2,..., n?

?1 ?1

i i i i i

y ?y y ?y

i?1 i ???i i?1

其中，? ?

h

i?1

,? ?1??,d ?6

h h

i i?1 .

i h ?h

i

i?1

i i i

h ?h

i

i?1

该方程有n?1个未知数，n?1个方程。针对不同的边界条件可以有相应的附加方程，

最常用到的是m ??,m ? ?.解出（1.1）及其附加方程得到m 再代进S (x)的

0 n i i

表达式，就得到了全部解。

2、分段线性插值：

所谓分段线性插值就是通过插值点用折线段连接起来逼近f(x).设已知节点

a?x ?x???x

?b上的函数值f

,f,?,f

,记h ?x ?x,h?maxh

,求一折

0

线函数I

h

1 n

(x)满足：

0 1 n

k k?1 k kk

10I

h

(x)?[a,b],

?20I(x)?f （k?0,1, ,n），

?

k k k

30I (x)在每个小区间[x,x ]上是线性函数。

h k k?1

则称I

h

(x)为分段线性插值函数。

3、插值理论：

设函数y=f(x)在区间[a,b]上连续，在[a,b]上有互异点x0,x1,…,xn处取值y0,y1,…,yn。如果函数φ(x)在点xi上满足φ(xi)=yi(i=0,1,2,…,n)，则称φ(x)是