matlab 关于内维尔插值问题

来源:百度知道 编辑:UC知道 时间:2024/06/02 06:31:26
利用下表所给数据用内维尔插值法求正弦积分对sint/t的定积分的负数,上限是无穷大,下线时X。在X=0.462的值
0.3 0.4 0.5 0.6 0.7
0.29850 0.39646 0.49311 0.58813 0.68122

你的问题看得不是很清楚。
我只做了一半,就是内维尔插值法。

%%%算法参:http://class.htu.cn/SZBJ/2/2_10.htm
x=[0.3 0.4 0.5 0.6 0.7];
y=[0.29850 0.39646 0.49311 0.58813 0.68122];
xx=0.462;

n=length(x);
f=y;
for i=1:n-1
for j=1:n-i
f(j)=((xx-x(j))*f(j+1)-(xx-x(j+i))*f(j))/(x(i+j)-x(j));
end
end

yy=f(1)