拉格朗日插值多项式MATLAB实现
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1、Lagrange 插值多项式源代码：
%未增加所谓的x,y维数判别，基本上y与x一一对应，维数都是一样的。
functionya=lag(x,y,xa)
%x 所有已知插值点
%y 插值点对应函数值
%xa所求点，自变量
%ya所求点插值估计量
ya=0;
mu=1;
%初始化
%循环方式求L系数，并求和：
fori = 1:length(y)
for j = 1:length(x)
if i ~= j
mu = mu * (xa - x(j) ) / ( x(i) - x(j) );
else
continue
end
end
ya = ya + y(i) * mu ;
mu = 1;
end