Université Sultane Moulay
Slimane Faculté des Sciences
et Techniques
TP d’analyse numérique
Compte rendu :
TP n° 2 :
Objectifs : maîtriser les commandes de base pour
crée un programme.
Encadré d’après la Prof [Link] , et l’étudiants :
ABDELWAHED NOURABI & MOHAMMED NAJI & SAIDA KHAFIFE .
Rédigé d’après : ABDELWAHED NOURABI & MOHAMMED NAJI & SAIDA KAFIFE
2/ Méthode de HÖRNER :
Algorithme 1 Evaluation :
function s = evaluation(a,z)
tic;
n = length(a);
s = 0;
for i = 1 : n
s = s + a(i) * z ^(i-1);
end
s;
toc;
end
Algorithm 2 HÖRNER:
function s = HORNER(a,x)
tic;
n = length(a);
s = a(n);
for i = n : -1 : 2
s = s * x + a(i-1);
end
s;
toc;
end
Rédigé d’après : ABDELWAHED NOURABI & MOHAMMED NAJI & SAIDA KAFIFE
3/ Matrice de Vandermonde :
function [Z,a] = vandermonde(X,Y)
n = length(X)
Z = [];
for i = 1 : n
for j = 1 : n
Z(i,j) = (X(i))^(j-1);
end
end
Z;
a = Z\Y'
plot([-2 -1 1 2],[10 4 6 3],'O')
hold on
fplot('4.5*x.^0 + 1.9197*x.^1 + 0.5*x.^2 + -0.9197*x.^3', [-2 2],'r')
end
Rédigé d’après : ABDELWAHED NOURABI & MOHAMMED NAJI & SAIDA KAFIFE
4/ Interpolation de Lagrange :
function [p] = lagrange(X,Y,s)
n = length(X)
p = 0;
for i = 1 : n
L = 1;
for j = 1 : n
if(i~= j)
L = L.*(s-X(j))./(X(i) - X(j));
end
end
p = p + L * Y(i);
end
end
Rédigé d’après : ABDELWAHED NOURABI & MOHAMMED NAJI & SAIDA KAFIFE
ALGORITHM 5 :
function P = interpole1(X,Y,x)
P = 0;
L = [];
L = Lagrange(X,x);
P = L.*Y;
end
ALGORITHM 6 :
function Pe = interpole2(a,b,n1,n2)
h1 = (a-b)/ (n1 + 1);
h2 = (a-b)/ (n2 + 1);
X1 = a : h1 : b;
X2 = a : h2 : b;
Y = exp(X1);
for i = 1 : length(X2)
Pe(i) = interpole1(X1,Y,X2(i));
end
%plot(X2,Pe,'.-',X2,exp(X2),'o-');
%legend('interpole','exp(x)');
%grid on
%xlabel('x')
%ylabel('y')
%title('intrepolation de exp(x) sur [1 2]')
end
Rédigé d’après : ABDELWAHED NOURABI & MOHAMMED NAJI & SAIDA KAFIFE
Rédigé d’après : ABDELWAHED NOURABI & MOHAMMED NAJI & SAIDA KAFIFE
