Bonjour,

Je dois utiliser le script ci-dessous permettant de calculer les coefficients de polynôme interpolé mais je n'arrive pas à en comprendre le fonctionnement.

Quelqu'un peut-il me commenter le script ?

Merci

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
function [l,L] = lagranp(x,y) 
N=length(x)-1; 
l=0; 
for m=1:N+1 
P=1; 
for k=1:N+1 
if k~=m, P=conv(P,[1-x(k)])/(x(m)-x(k)); 
end 
end 
L(m,:) = P; 
l=1+y(m)*P; 
end