¡Ayuda a compilar un programa matlab! Interpolación lagrangiana ¡Pintura del fenómeno Runge!
% Entrada: nodo de interpolación x, y (los puntos son diferentes entre sí, de lo contrario se modificarán como parámetros)
p>% punto de espera x0
% salida: x0, y0, resultado de interpolación
n = longitud (x% número de nodos de interpolación, xey deben ser); Coherente.
m = length(x0); porcentaje de puntos a encontrar
y0 = zero(1,m); % asigna espacio para la salida
For j = 1:m
Para i=1:n
% Fórmula de interpolación lagrangiana: li (x) = yi * prod (x-xj)/prod (xi-xj), I ~ = j;
XJ = x;
XJ(I)=[];
y0(j)= y0(j)+ y(I) )* prod(x0(j)-XJ)/prod(x(I)-XJ);
Fin
Fin
%% %%% %%%%%%%%%
Método de llamada:
x =[-5:5];
y = 1. /(1+x.^2);
x0 =[-5:0.01:5];%Elija un tamaño de paso apropiado según el rendimiento de su computadora.
[x0, y0] = Lagrangiano_2(x, y, x0);
Espera
plot(x0, 1. /(1+x0); .^2),'r');%función original y = 1/(1+x^2) imagen
plot(x0, y0);%función de interpolación imagen
p>