拉格朗日插值多项式

Question

我有分段线性函数

我尝试为2、6、14个主机构建Lagrange的多项式。 但是对于我的主机，我原来的功能不正确：

我认为我的多项式是对的，但是我很难计算和建立不准确度的度量

a = 1;
b = -1;
c = -1;
d = 1;
x = -1:0.33:1;
y = [];
for i = 1 : length(x)
 if (x(i) <= c/2)
    y(i) = x(i)+ 1;
elseif(x(i) >= d/2)
    y(i) = x(i) - 1;
 else
    y(i) = -x(i);
 end
end
plot(x, y);
hold on
nx = -1:0.01:1;
ny=lagrange(x,y,nx);
plot(x,y,'bo');
plot(nx,ny,'r');

function ny=lagrange(x,y,nx)
N = length(x);
ny=zeros(size(nx));
for k=1:N
  t=ones(size(nx));
   for j=[1:k-1, k+1:N]
      t=t.*(nx-x(j))/(x(k)-x(j));
  end
  ny = ny + y(k)*t;
end

Answer 1

好吧，我犯了一些错误。 首先，我将错误的参数传递给函数。 另外，我没有正确计算原始函数的值。 “原始”函数计算原始图形的值。 “ lagrange”函数计算值多项式。

x = -1:0.01:1; 
y = []; 
y = original(x); 
subplot (2, 1, 1); 
plot(x, y); 
xlabel('x'); 
ylabel('y'); 
hold on; 
grid on; 
%hosts
xx = -1:0.5:1;
%values of original func in hosts 
yy = original(xx); 
plot(xx,yy,'o'); 
nx = -1:0.01:1; 
ny=lagrange(xx,yy,nx);     
plot(nx,ny,'r'); 
subplot (2, 1, 2); 
%measure of inaccuracy 
plot(x,ny-y, ''); 
grid on;