interpolation de Lagrange
Bonjour ,
J'ai besoin d'un coup de main
Je dois programmer une fonction qui prend en entrée deux vecteurs (un vecteur pour les abcisses des points considérés et un vecteur pour les cordonnées ) et donne en sortie le polynome interpolateur en ces points calculé par la méthode d Lagrange .
donc je doit appliquer ça a la fonction f(x)=exp(x) sur l'intervalle [-1,1] avec points d'interpolation équidistants puis tracer la fonction considérée et le polynome sur le meme graphique:et utiliser linspace pour les définir
voila mon essai
Code:
1
2
3
4
5
6
7
8
9
10
11
12
| function poly=lagrange(X,Y,XX)
z=zeros(size(XX));
n=length(Y);
for i=1:n
z=ones(size(XX));
for j=1:n
if i~=j
z=z.*(XX-X(j)/X(i)-X(j));
end
end
poly=poly+Z(i).*Y(i);
end |
