Discussion :

[mari0marie]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[Invité]

Bonjour,

De petites étourderies:

• poly=poly+Z(i).*Y(i); à remplacer par poly = poly + z*Y(i); avec Z minuscule (z de la taille de XX), et le .* non-obligatoire car Y(i) est un scalaire.
• L'oubli des parenthèses z=z.*(XX-X(j))/(X(i)-X(j));
• La première ligne z=zeros(size(XX)); est inutile.
• La variable poly à initialiser