大家好,又见面了,我是你们的朋友全栈君。
lagrange插值法:求拉格朗日插值多项式matlab实现(内附代码及例题)
关于拉格朗日插值法相关理论知识,在这里小编不在赘述,请不明白的小伙伴自行百度。小编只负责给出matlab源码。
**例题:**看下面例题(如图):
matlab代码:
代码语言:javascript复制%%%% 求拉格朗日多项式及基函数 %%%%
%%%% Liu Deping %%%%
%%%% 2020.06.14 %%%%
%输入的量:n 1个节点(x_i,y_i)(i = 1,2, ... , n 1)横坐标向量X,纵坐标向量Y
%输出的量:n次拉格朗日插值多项式L和基函数l
X=input('请输入横坐标向量X:nX='); %输入的数据为一维数组,例如:[1,3,4,5](下同);
Y=input('请输入纵坐标向量Y:nY=');
m = length(X);
L = ones(m,m);
for k = 1 : m
V = 1;
for i = 1 : m
if k ~= i
V = conv(V,poly(X(i))) / (X(k) - X(i));
end
end
L1(k, :) = V;
l(k, :) = poly2sym(V);
end
fprintf('基函数为:n');
for k=1:m
fprintf('q%d(x)=%sn',k,l(k));
end
L = Y * l;
fprintf('拉格朗日多项式为:nP(x)=%sn',L);
操作步骤: step1:复制代码至matlab编辑器,并点击运行; step2:按照提示,输入坐标向量(输入完坐标向量后敲回车键); 上面例题的运行结果如下:
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。