Scribd
Importer
6 vues2 pages

Polynôme d'interpolation de Lagrange

Le document décrit le polynôme d'interpolation de Lagrange pour approximer une fonction à partir de quelques points échantillons. Le polynôme est exprimé comme une somme de termes, chacun associé à un point d'échantillonnage.

Transféré par

selmaaaab2703
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
6 vues2 pages

Polynôme d'interpolation de Lagrange

Le document décrit le polynôme d'interpolation de Lagrange pour approximer une fonction à partir de quelques points échantillons. Le polynôme est exprimé comme une somme de termes, chacun associé à un point d'échantillonnage.

Transféré par

selmaaaab2703
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

Le polynôme d'interpolation de Lagrange est :

(921216758193973*((4*X)/9 - 14/9)*(X - 8)*(X - 37/8)*(X - 55/8))/8349416423424 -


(3588903522234403*((8*X)/9 - 28/9)*(X - 8)*(X - 23/4)*(X - 55/8))/300578991243264 +
(4660585534417883*((8*X)/9 - 37/9)*(X - 8)*(X - 23/4)*(X - 55/8))/4809263859892224 +
(364177552073791*((2*X)/9 - 7/9)*(X - 23/4)*(X - 37/8)*(X - 55/8))/1043677052928 -
(4256321070142085*((8*X)/27 - 28/27)*(X - 8)*(X - 23/4)*(X - 37/8))/12524124635136
clc
clear all
f = @(x) exp(x);
a = 3.5;
b = 8;
n = 4; % Nombre de sous-intervalles
h = (b - a) / n; % Calcul du pas
x = a:h:b;
syms X;
P = 0;
for i = 1:length(x)
% Calcul du produit L_i(x)
L = 1;
for j = 1:length(x)
if j ~= i
L = L * (X - x(j)) / (x(i) - x(j));
end
end
% Ajout du terme au polynôme P(x)
P = P + f(x(i)) * L;
end

disp('Le polynôme d''interpolation de Lagrange est :');


disp(P);
X_plot = linspace(a, b, 1000);
Y_plot = subs(P, X, X_plot);
% Tracé du polynôme
figure;
plot(X_plot, Y_plot, 'LineWidth', 1.5);
hold on; % Pour conserver le tracé existant
plot(x, f(x), 'ro', 'MarkerSize', 8); % 'ro' pour des cercles rouges (marqueurs)
xlabel('x');
ylabel('f(x)');
title('Polynôme d''interpolation de Lagrange avec les points d''interpolation');
legend('Polynôme d''interpolation', 'Points d''interpolation');

Vous aimerez peut-être aussi