Polynôme d'interpolation de Lagrange MATLAB Résolu/Fermé

Question

Bonjour,
je veux faire un programme qui calcule le polynome d'interpolation avec la methode d'agrange dans MATLAB 
mon programme d'essaye me donne toujours une erreur ;voici le code

function YY=lagrange(X,Y,XX)
YY=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
YY=YY+Z(i).*Y(i);
end 

comment faire marcher cette fonction merci de vouloire m'aider?

Abdou · Accepted Answer

normalement voici le programme
 
function []=fonc_lagrange(a,x,y)
dif= fopen('lagrange.res','w');
n=length(x)
p=0;
l=1;
for i=1:n;
for j=1:n;
   if j~=i
        l=((a-x(j))/(x(i)-x(j))*l);
    end 
end 
fprintf(dif,'%d
',l);
fprintf(dif,'%d
',a);
 disp(['lag(',num2str(i),')=',num2str(l)]) 
t=(y(i)*l);
end
p=p+t;
disp(['Pn(',num2str(a),')=',num2str(p)]) 
fprintf(dif,'%d
',p);
fclose(dif);

gumuslive · Answer

merci beaucoup je vais la essayer

Sacabouffe · Answer

Salut

Comme ce sujet est remonté, j'en profite pour signaler que sous Matlab, l'utilisation de polyfit suffit pour calculer le polynôme d'interpolation de Lagrange associé aux valeurs x et y.
P=polyfit(x,y,numel(x)-1) Ciao

Mourad. H · Answer

%%% le programme fonctionne correctement
 
function p=lagrange(xd,x,y)
a=xd;
p=0;
n=length(x); 

 
for i=1:n; 
    l=1;
for j=1:n; 
if j~=i 
l=((a-x(j))/(x(i)-x(j))*l); 
end 
end 
t=(y(i)*l);
p=p+t;
end

Polynôme d'interpolation de Lagrange MATLAB

4 réponses

Discussions similaires