fonction non linéaire par raphson

**membreComplexe12** · 10/03/2010, 14h32

Envoyé par crippella

bonjour ,
j'ai un problème de calculer la dérivée en MATLAB je ne sais pas quelle est la commende convenable ou bien l'expression si elle y'en a pas il ya pas
svp aidez moi je suis coinsé

moi je suis comme toi je ne sais pas, pour cela il faudrai que matlab est une partie calcul formel, et je ne sais pas si c'est le cas.....

**phryte** · 10/03/2010, 14h37

Bonjour.

matlab est une partie calcul formel

C'est ce qui est dans cette discussion.

http://www.developpez.net/forums/d87...ewton-raphson/

**membreComplexe12** · 10/03/2010, 14h39

Envoyé par phryte

Bonjour.

C'est ce qui est dans cette discussion.

http://www.developpez.net/forums/d87...ewton-raphson/

merci

**crippella** · 13/03/2010, 21h25

merci pour votre réponses
finalement le programme est bien executer mais il ya une toute petite erreur c'est que par exemple si je fais entrer un fonction qui s'annule deux fois ou plus le programme calcule un point dont lequel la fonction il s'annule et pas plus,
j'espere que vous avez compris

voila le code MATLAB

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
function NewtonRaphson(x0,f)
 
clc
format long
x(1)=x0; 
e=1e-10;
n=5000;
for i=2:n 
df=(f(x(i-1)+e)-f(x(i-1)))/e; 
%on calcule une valeur approchée de la dérivée de f en x(i-1).
x(i)=x(i-1)-f(x(i-1))/df;
if abs(x(i)-x(i-1))<=e
disp('La valeur finale est:')
xp=x(i)
break 
end
end
disp('les valeurs des itérations successives sont:')
x'
figure(1)
clf
hold on 
plot(x) 
plot(x,'*r')
title('Approximations successives de la methode de New.-Raph.');
xlabel('numero de l''iteration')
ylabel('valeur de l''approximation')
 
end

j'attends votre réponses
et merci d'avance

**phryte** · 14/03/2010, 11h28

Bonjour.
Une proposition (à améliorer) :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
clear
clc
xi=[0 10];
pas=0.1;
f=inline('x^3-6*x^2+11*x-6');
e=1e-12;
xp=0;
hold on 
for x0=xi(1):pas:xi(2)
    clear x
    x(1)=x0;
i=2;
 while abs(f(x(i-1)))>=e   
df=(f(x(i-1)+e)-f(x(i-1)))/e;
x(i)=x(i-1)-f(x(i-1))/df;
i=i+1;
end
xp=[xp x(end)+e];
plot(1:length(x),x) 
end
disp('Les valeurs finales sont :')
unique(fix(1000*xp(2:end))/1000)
title('Approximations successives de la methode de New.-Raph.');
xlabel('numero de l''iteration')
ylabel('valeur de l''approximation')
grid

**crippella** · 14/03/2010, 21h33

merci beaucoup pour la réponse,
je sais que je pose beaucoup de questions ,mais est ce que vous pouvez m'expliquer la boucle for et ces instructions

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
xp=[xp x(end)+e];
plot(1:length(x),x)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

unique(fix(1000*xp(2:end))/1000)

merci d'avance

**phryte** · 15/03/2010, 08h17

Bonjour.

xi=[0 10];
pas=0.1;

Plage de recherche en x et pas de recherche.
La boucle for effectue tous les calculs de recherche en initialisant le x0 à chaque itération.

xp=[xp x(end)+e];

Résultats en enlevant le premier terme d'initialisation et en calculant par excès.

plot(1:length(x),x)

Tracé de tous les calculs en ordonnant les valeurs.

unique(fix(1000*xp(2:end))/1000)

Extrait des valeurs uniques des racines (car il y a plusieurs valeurs communes avec arrondi).

**crippella** · 15/03/2010, 13h55

bonjour;
est ce que vous pouvez m'expliquer de plus les instructions de la boucle for et pourquoi vous avez supprimer le x parce que j'ai pas bien comprendre
et merci d'avance