实验一函数插值方法
一、问题提出
nge公式求其插值多项式或分段二次Lagrange插值多项式。
数据如下
(1)
0.65 04 0.55 0.8 0.9 1.05
O 5
0.6967 0.9 1.0 0.4107 0.5781 Y 1.2538
5 5 5 2
(x) 求五次Lagrange多项式 和分段三次插值多项式,计算 f(0.596)f(0.99)
的值。(提示:结果为
f(0.596)0625732f(0.99)1.05423
21
3 2 A 5 大
0.135 0.05 0.01 0.00 0.36 0.000.00
8 7 8 21
试构造Lagrange多项式 L.(x) 计算的
f(1.8)0.164762f(6.15)0.001266
二、要求
1、利用Lagrange插值公式
L(x)= X-x
0 编写出插值多项式程序
2 给出插值多项式或分段三次插值多项式的表达式:
3 根据节点选取原则,对问题(2)用三点插值或二点插值,其结果如何:
4、对此插值问题用Newton插值多项式其结果如何。Newton插值多项式如下
(x-X;)
O.jk
f(x)
其中: 107
三、自的和意义
1、学会常用的插值方法,求函数的近似表达式,以解决其它实际问题:
2、 明确插值多项式和分段插值多项式各自的优缺点:
3. 熟悉插值方法的程序编制:
4、如果绘出插值函数的曲线,观察其光滑性。
四、实验学时:2学时
五、实验步骤:
1.进入C或matlab开发环境
2.根据实验内容和要求编写程序:
function yy=lagrange(x1,y1,xx);
x1=(0.4,0.55,0.65,0.8,0.9,1.05);
y1=(0.41075,0.57815,0.69675,0.90,1.00,1.25382);
xx=(0.596,0.99);
n=length(x1)'
l=ones(1,n);
for k=1:n
for j=1:n
if j~=k
l(k)=l(k).*(xx-x1(j))./(x1(k)-x1(j));
end
end
end
sums(y1.*l);
这是我的代码,但是运行不出来