COMPTE
RENDU
FAIT LE
01/05/2024
Fait Par :
Makhloufi Fares
Goupe 03 / Sous Group 5
Compte rendu 3 et 4 : lagrange de newton
1_ les commentaire et codes:
%effacer le contenu de la fenetre de commande (effacer I ecran):
clc
%Effacer les variables de I espace de travail et les liens des
fonctions
utilisees: clear all
% Fermer tous les fenetre de figures qui ont été ouvertes: Close all
X=input ('entrez les valeurs de xi')
Y=input('entrez les valeurs de yi')
Xx=input ('entrez les valeurs de xx')
%déterminer la longueur d'un vecteur x" et représentant le nombre
d'éléments dans le vecteur x: n=lenght(x)
% Initialisation d'une matrice b de taille n x n
avec des éléments initialisés à zéro
b = zeros(n, n);
% Remplissage de la première colonne de b avec des uns
b(:,1) = 1;
% Calcul des éléments de la matrice b en utilisant la méthode de
Newton
for j = 2:n
for i = j:n
% La formule suivante est utilisée pour calculer b(i,j)
en utilisant les éléments précédents de la matrice
b(i, j) = b(i, j-1)* (x(i) - x(j-i));
end
end
% Transposition du vecteur y
y1 = y';
% Calcul des coefficients a en utilisant la matrice b et le vecteur y1
a=b\y1;
% Initialisation d'un vecteur de longueur length(x) rempli de 1
Nones(length(x), 1);
% Calcul des éléments du polynôme de Newton en utilisant les
coefficients a et le vecteur N
for j = 2:length(x)
N(j) = N(j-1)*(xx-x(j-1));
end
% Calcul du polynôme de Newton en multipliant les coefficients a
avec
les éléments du vecteur N et en sommant les résultats
PN=A' * N;
% Tracé du graphique du polynôme de Newton
plot(xx, PN, 'r', x, y, 'o');
hold on;
% Définition de la fonction f(x) = x^2+1
f = @(z) z^2 + 1;
% Tracé du graphique de la fonction f(x) dans l'intervalle [0,3]
% Affichage des valeurs de PN et f(xx)
fprintf('p(%d)=%f\n', xx, PN);
fprintf('f(%d)=%f\n', xx, f(xx));
3_la foction:
f(x)-exp(x)+1
x=[0 134]
y=[3 4.72 22.09 56.6]
le code:
clc;clear all; close all;
x=[0 134];
y=[3 4.72 22.09 56.6];
xx=2;
n=length(x);
b=zeros(n,n);
b(:,1)=1;
for j = 2:n
for i = j:n
b(i, j)b(i, j-1) (x(i) - x(j-i));
end
end
y1 = y':
a=b\y1;
N = ones(length(x), 1);
for j = 2:length(x)
N(j) = N(j-1)*(xx - x(j-1));
end
PN=A' N;
plot(xx, PN, 'r, x, y, 'o');
hold on;
f= @(z) exp(z)+ 1;
fplot(f.[0 3], 'k')
fprintf('p(%d)=%f\n', xx, PN);
fprintf('f(%d) = %f\n', xx, f(xx));
les resultats et le graphique :
PN=7.9400
P(2)=7.940000
F(2)=8.389056
Tp N°4
1-Script Matlab
cle,clear all,close all;
%Saisir les valeur de xi et yi
X=[0123];
y=[101510]:
% Suisir la valeur de xx
Xx=25;
% le number des d'éléments dans matrice ou un
Vecteur
N=length(x);
%l'element dneutre de la somme
P=0
% for (les condution de la grange)
for j= 1:n
L=1;
For j=1:n
if(i~=j) L=(xx-x(j))/(x(i)-×(j))*L;
end
end
p=p+y(i)*L
end
%la fonction f(x)
F=@(x) x^2+1:
%afficher le graph de f(x) enter literval [03]
Fplot(f,[05],'+r')
%tracer touts les points dans le méme figur
hold on
%afficher les points dans la figur
Plot(x,y,'o',xx,p. 'k'xx.f(xx)+r)
%afficher la solution de p(x)
fprintf('p(%f)=%f\n',xx,p)
%afficher la soluti de f(x)
Fprintf('f(%f)=%f\n',xx, f(xx))
3-le code de f(x)=x^2+1 x=2.5
cle,clear all,close all;
X=[0 1 2 3];
y=[0 15 10];
% Suisir la valeur de xx
Xx=2.5;
N=length(x);
P=0
for j= 1:n
L=1;
For j=1:n
if(i~=j) L=(xx-x(j))/(x(i)-x(j))*L;
end
end
p=p+y(i)*L
end
F=@(x) x^2+1:
Fplot(f,[03],'+r')
hold on
Plot(x,y,'o',xx,p.*'k'xx.f(xx)+r)
fprintf('p(%f)=%f\n',xx,p)
Fprintf('f(%f)=%f\n',xx, f(xx))
3-le graphe :
4-les resultats:
P(2.500000)=7.500000
F(2.500000)=7.250000
5-les resultats pour la function f :exp(x)=+1
Code de programme
Clc:clear all;close all;
x=[0 134];
y=[3 4.72 22.09 56.6];
xx=2;
n=length(x);
p=0;
for i=1:n
L=1;
for j=1:n
if(i~=j)
L=(XX-X(j))/(x(i)-x(j))* L;
end
p=p+y(i) *L;
end
f=@(z) exp(z)+1:
fplot(f,[0 3], 'r')
hold on
plot(x,y, 'o',xx,p,'*k',xx, f(xx),'r+')
fprintf('p(%t) %f\n',xx,p)
fprintf('f(f)=%f\n',xx, f(xx))
les resultat et la figeur:
p(2.000000)=7.940000
f (2.000000)=8.389056
le graphe :