将拉格朗日插值法推广到高维情况

转载已于 2022-09-28 23:11:54 修改 · 502 阅读

CC 4.0 BY-SA版权

原文链接：https://kns.cnki.net/KXReader/Detail?invoice=PahdFD%2FXpAIkHcQCk0UGqcuS0pOxQMo5NJ7b3NQ6WN1J59dgnOxqJb7RHBBORfJQUqdNstOoErAkas2WNIIw%2FNIOjOf5SmKH0dwLMIX8H%2BiHz0Ormr8rSO72WFlrxzgdX4LwaoHhbLAfxBwIceE0T8dBIGz20xMF4RIcAU0%2FeJE%3D&DBCODE=CJFD&FileName=GJJS20

文章标签：

#matlab #算法

于 2022-09-28 23:09:56 首次发布

分享一篇关于高维拉格朗日插值法的论文，通过MATLAB代码展示了如何在二维和三维空间进行复杂数据的插值操作，适用于数值计算和数据拟合的实际应用。

无意中在知网发现了一篇关于高维拉格朗日插值法研究的论文，故将之分享给大家~

中国知网-登录

以下是我根据该篇论文所编写的MATLAB代码：

%realize Lagrange polynomial-2
function [x,v,y,u,z,w] = lagpol2(a,b,n)
s = [a:(b-a)/(3*(n-1)):b];
t = [a:(b-a)/(4*(n-1)):b];
[x,y] = meshgrid(s);
[v,u] = meshgrid(t);
z = cos(x).*cos(y);
lens = length(s);
lent = length(t);
w = zeros(lent);
for i0 = 1:lent
for j0 = 1:lent
%二维插值开始
g = linspace(0,0,lens);
for k = 1:lens
%对x的一维插值开始
for i1 = 1:lens
temp1 = 1;
for j1 = 1:lens
if j1~= i1
temp1 = temp1 * (v(i0,j0) - x(1,j1))/(x(1,i1) - x(1,j1));
end
end
g(k) = g(k) + temp1 * z(k,i1);
end
%对x的一维插值结束
end
%对y的一维插值开始
for i2 = 1:lens
temp2 = 1;
for j2 = 1:lens
if j2 ~= i2
temp2 = temp2 * (u(i0,j0) - y(j2,1))/(y(i2,1) - y(j2,1));
end
end
w(i0,j0) = w(i0,j0) + temp2 * g(i2);
end
%对y的一维插值结束
%二维插值结束
end
end