enrichir le TP avec un exemple en plus (binôme ADJOUDJ+RAHIL ST 22)
2 participants
Page 1 sur 1
enrichir le TP avec un exemple en plus (binôme ADJOUDJ+RAHIL ST 22)
fonction exp
résoudre l'équation suivante sous MATLAB par 2 méthodes (Newton, et Point fixe):
f(x)= exp(x-3)+4.
1- méthode de NEWTON
il faut que je transforme la fonction en suite. f(x) devient f(xn).
je dois résoudre le système suivant :
x0 = ce que tu cherches
x(n+1)=xn-f(xn)/f'(xn)
je devrais essayer le code suivant:
CODE
e=1e-10;
for i=m:n
x(i+1)=x(i)-f(x(i))/diff(f(x(i)));
if abs(x(i+1)-x(i))<=e %condition d'arret
xp=x(i);
fprintf('xp=%f\n',x(i));
return;
end
end
résoudre l'équation suivante sous MATLAB par 2 méthodes (Newton, et Point fixe):
f(x)= exp(x-3)+4.
1- méthode de NEWTON
il faut que je transforme la fonction en suite. f(x) devient f(xn).
je dois résoudre le système suivant :
x0 = ce que tu cherches
x(n+1)=xn-f(xn)/f'(xn)
je devrais essayer le code suivant:
CODE
e=1e-10;
for i=m:n
x(i+1)=x(i)-f(x(i))/diff(f(x(i)));
if abs(x(i+1)-x(i))<=e %condition d'arret
xp=x(i);
fprintf('xp=%f\n',x(i));
return;
end
end
adjoudjIbtissam- Messages : 18
Date d'inscription : 06/10/2015
Age : 30
Re: enrichir le TP avec un exemple en plus (binôme ADJOUDJ+RAHIL ST 22)
méthode de point fixe :
Pour finir je pourrais utiliser quelque chose de ce genre là pour la méthode du point fixe.
Code :
function x1 = FPI(f,x1,tol)
% FPI(f,x0) Iterations du point fixe
%
if nargin == 2
tol = 1e-8;
end
it = 0;
itmax = 100;
x0 = realmax ;
while ~converged (x0,x1,tol)
x0 = x1;
x1 = feval(f,x0);
it = it + 1;
if it > itmax
error('Maxit in FPI');
end
end
avec converged une fonction qui définie le convergence
Code
function res = converged(y0,y1,tol)
res = all(abs(y1-y0)./(abs(y0)+1) < tol);
end
Pour finir je pourrais utiliser quelque chose de ce genre là pour la méthode du point fixe.
Code :
function x1 = FPI(f,x1,tol)
% FPI(f,x0) Iterations du point fixe
%
if nargin == 2
tol = 1e-8;
end
it = 0;
itmax = 100;
x0 = realmax ;
while ~converged (x0,x1,tol)
x0 = x1;
x1 = feval(f,x0);
it = it + 1;
if it > itmax
error('Maxit in FPI');
end
end
avec converged une fonction qui définie le convergence
Code
function res = converged(y0,y1,tol)
res = all(abs(y1-y0)./(abs(y0)+1) < tol);
end
adjoudjIbtissam- Messages : 18
Date d'inscription : 06/10/2015
Age : 30
Re: enrichir le TP avec un exemple en plus (binôme ADJOUDJ+RAHIL ST 22)
je vous ai envoyés sous forme de fichier joint des images mais il ne marche pas il me signale se problème (( Le fichier envoyé n'est pas valide: la taille du fichier joint dépasse la limite autorisée (0 ko).))
adjoudjIbtissam- Messages : 18
Date d'inscription : 06/10/2015
Age : 30
mezrara khadidja- Messages : 3
Date d'inscription : 19/10/2015
Page 1 sur 1
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum
|
|