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

PAGE

PAGE1

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

组号 班级 学号 姓名 分数

一：实验目的

1、掌握计算多项式插值与拟合的方法。

2、掌握样条插值与拟合的方法。

3、通过实例学习插值方法与拟合方法的联系与区别。

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

插值的原理：

函数y?f?x? 在区间 [a,b] 上有定义，在点a?x

0

?x x

?1 n

?

?b 上的值为

y,y

1 2

, y

?n

?

，构造一个简单的函数P?x? ，使P?x

i

??y

i

(i?0,1,?,n)（1）成立，

?? ?? ? ?

则Px 就是插值函数，再利用 Px计算插值，即y*?Px*

L?x??axn

a xn?1

???a

x?a

l?x

拉格朗日多项式插值

????1,k?j

?

n n n?1

1 0 插值基函数

j k ?0,k?j (j,k?0,1,?,n) （2） 可得到 n 次插值基函数

l?x??

(x?x) (x?x

?0

?

k?1

)(x?x

k?1

) (x?x)

??nn

?

?

n

?k (x

?

k

x) (x

0 k

x

k?1

)(x

k

x

k?1

) (x

k

?x) （(k?0,1,?,n)） 它满足条件

（2），于是满足条件（1）的插值多项式格朗日插值多项式。

拟合的原理：

L?x?

n

可表示为L

????n

nxk?0

n

x

yl?x?

kk

就是拉

?已知一组二维数据，即平面上的n 个点(x,y)i?1,2, ,3 ,寻求一个函数

?

i i

y?f?x? ，使 f?x? 在某种准则下与所有数据点最为接近，即曲线拟合的最好。三次样条插值的基本思路为：

设在区间[a,b]上，已给n+1个互不相同的结点

a?x

0

?x???x ?b

1 n

而函数y?f?x?在这些结点的值f

?x??y

i i

，i?0,1,?,n.如果分段表示的函数S(x)满

?足下列条件，就称S(x)为f(x)的基点x,x, ,x

?

0 1 n

的三次样条插值函数，简称三次样条。

⑴S(x)在子区间[x,x

i i?1

]的表达式S

i

?x?都是次数不超过3的多项式；

⑵S?x??y ；

i i

⑶S?x?在整个区间[a,b]有连续的二阶导数。

i

三：实验问题及方法、步骤

拟合下列数据

x

0.1

0.2

0.15

0

-0.2

0.3

y

方法一

0.95

多项式插值与拟合

0.84

0.86

1.06

1.50

0.72

P=polyfit(x,y,k)用k次多项式拟合向量数据（x,y），返回多项式的降幂系数。当kn-1时，polyfit实现多项式插值。这里n是向量维数。

⑴二次拟合多项式p(1)x2+p(2)x+p(3)

clear;x=[0.1,0.2,0.15,0,-0.2,0.3]

y=[0.95,0.84,0.86,1.06,1.50,0.72]

p=polyfit(x,y,2)p=

1.7432 -1.6959 1.0850

⑵五次拟合多项式（等价于多项式插值）

xi=-0.2:0.01:0.3

yi=polyval(p,xi);subplot(2,2,1)

plot(x,y,o,xi,yi,k)

title(polyfit)

绘制出的图像为

（图1）

p=polyfit(x,y,5)p=

1.0e+003*

-1.8524 0.7560 0.0079 -0.0275 0.0010 0.0011

yi=polyval(p,xi);subplot(2,2,2)

plot(x,y,o,xi,yi,k)

title(ployinterp)

绘制出的图像为

方法二样条插值和拟合

（图2）

clear;x=[0.1,0.2,0.15,0,-0.2,0.3]

y=[0.95,0.84,0.86,1.06,1.50,0.72]

pp=spline(x,y)x=

y=

0.1000

0.9500

0.2000

0.8400

0.1500

0.8600

0

1.0600

-0.2000

1.5000

0.3000

0.7200

pp=

form:

pp

breaks:

[-0.200000.10000.15000.20000.3000]

coefs:

[5x4double]