Código polinomios de Lagrange ;
function[c]=lagrangeinterpol(x,y)
%la forma de inicializar a esta varible va a ser para que el resultado se
%almacene en la matriz c, y esta matriz c va a representar los
coeficientes
%del polinomio de lagrange es decir si tenemos tres puntos nuestra matriz
%ce nos dara los coeficientes de ese polinomio
%generar procedimeinto
n1=length(x);% tamaño del vector x
n=n1-1; %grado del polinomio
L=zeros(n1,n1); %matriz de ceros de tamaño n1,n1. matriz cuadrada donde
voy a almacenar todos los coeficientes de lagrange; es decir los valores
que se vana ir actualizando, matrix de los cioesficientes de lagrange
for k=1:n+1
v=1; %parametro de ayuda para multiplicar todos los factores de los
coeficientes de lagrange
for j=1:n+1
if k~=j
v=conv(v,poly(x(j)))/(x(k)-x(j));
end
end
L(k,:)=v
end
c=y*L
plot(x,y)
RESULTADOS OBTENIDOS AL CORRER EL PROGRAMA:
Ítem a
>> x=[3 7]
x =
3 7
>> y=[5 -1]
y =
5 -1
>> lagrangeinterpol(x,y)
L =
-0.2500 1.7500
0 0
