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

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

Newsletters