Question

已知 利用lagrange插值法求 的值

Answer 1

楼主的程序基本没错，就是在t = t .* (xx - x(1,i))./(x(1,j) - x(1,i));这一行里xx前面多了个（
以下程序正确：（要是只做2次插值的话令n=3就行了）
function yy = Nlagrange(x,y,xx)
yy = 0;
j = 1;
n = 3;
while j <= n
t = 1;
for i = 1:n
if i ~= j
t = t .* (xx - x(1,i))./(x(1,j) - x(1,i));
end
end
yy = yy + t .* y(1,j);
j = j + 1;
end
end
另外输入的变量y多个逗号就不吐槽了，但采用大括号很令人费解，一般向量都是用中括号的。
像这样在Command Window输入：
x = [4 9 16];
y = [2 3 4];
xx = 7;
yy = Nlagrange(x,y,xx)
结果：
yy =
2.6286