我添加了输出函数称为z。
clc;明确
x = (2, 1, 0, 1, 2, 3];
y =(1, 4, 11日,16日,13日,4];
d0 = difdiv (x, y);
d = d0 (1);
xx = linspace (2, 2, 10);
z = 0(1、长度(xx));
为i = 1:长度(xx)
z (i) =拉格朗日(xx (i), x, y);
结束
z
函数z =拉格朗日(a, x, y)
d_before = difdiv (x, y);
d = d_before (1);
z = deflatedpolylagrange (a, x, d);
结束
函数z = deflatedpolylagrange (x,间,d)
n_d =长度(d);
n_x =长度(间);
如果n_d ~ = n_x
disp (“不可能的”。)
结束
v (1) = d (n_x);
为我= 2:n_x
v (i) =(张)* (x -间(n_x - i + 1)) + d (n_x - i + 1);
结束
z = v(结束);
结束
函数一个= difdiv (x, y)
n =长度(x);
一个= 0 (n, n);
(:1)= y ';
为j = 2: n
为i = 1: (n - j + 1)
(i, j) = ((i + 1, j - 1) - (i, j - 1)) / (x (i + j - 1) - x (i));
结束
结束
结束