牛顿插值流程图
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
I=zeros(1,N);
last_I=I;
NIout=0;
t=1;
fori=1:N
I(1)=ym(i);
forj=2:i
I(j)=(I(j-1)-last_I(j-1))./(xm(i)-xm(i-j+1));
end
last_I=I;
NIout=NIout+I(i)*t;
t=t*(x-xm(i));
1.流程图
2.算法
functionnewtonInterpolation
xm=-4:0.01:4;
ym=1./(1+xm.^2+xm.^4);
plot(xm,ym,'k');
set(gca,'color',[0.7,1,1]);
holdon;
forN=[3,5,7,9,11,13,15]
n=8/(N-1);
end
end
end
3.runge现象
xm=-4:n:4;
ym=1./(1+xm.^2+xm.^4);
Y=zeros(1,801);
forX=-4:0.01:4;
Y(round((X+4)./0.01+1))=NI(xm,ym,X);
endБайду номын сангаас
X=-4:0.01:4;
plot(X,Y);
holdon;
end
functionNIout=NI(xm,ym,x)
last_I=I;
NIout=0;
t=1;
fori=1:N
I(1)=ym(i);
forj=2:i
I(j)=(I(j-1)-last_I(j-1))./(xm(i)-xm(i-j+1));
end
last_I=I;
NIout=NIout+I(i)*t;
t=t*(x-xm(i));
1.流程图
2.算法
functionnewtonInterpolation
xm=-4:0.01:4;
ym=1./(1+xm.^2+xm.^4);
plot(xm,ym,'k');
set(gca,'color',[0.7,1,1]);
holdon;
forN=[3,5,7,9,11,13,15]
n=8/(N-1);
end
end
end
3.runge现象
xm=-4:n:4;
ym=1./(1+xm.^2+xm.^4);
Y=zeros(1,801);
forX=-4:0.01:4;
Y(round((X+4)./0.01+1))=NI(xm,ym,X);
endБайду номын сангаас
X=-4:0.01:4;
plot(X,Y);
holdon;
end
functionNIout=NI(xm,ym,x)