拉格朗日插值法MATLAB程序

格式：doc
大小：11.00 KB
文档页数：1

拉格朗日插值法MATLAB程序

下载文档原格式

/ 2

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

function lagrange

a=-5;b=5;n=10;h=(b-a)/n;

X=a:h:b;

Y=1./(1+X.^2);

X

Y

X0=[-2.3 -1.4 -0.92 0.2 0.6 1.4 1.9 ];

n=length(X);

m=length(X0);

for i=1:m

z=X0(i);

s=0.0;

for k=1:n

p=1.0;

for j=1:n

if j~=k

p=p*(z-X(j))/(X(k)-X(j));

end

end

s=s+p*Y(k);

end

Y0(i)=s;

plot(X0(i),Y0(i),'bo');hold on

end

plot(X,Y,'k-');

拉格朗日插值法matlab