Discussion :

[crippella]

Bonjour,
svp est ce que vous pouvez m'indiquer comment faire un programme sur MATLAB ce programme doit demander à l'utilisateur d'entrer une fonction non linéaire quelconque et le programme va calculer ou elle s'annule en appliquant la methode de newton raphson qui permet d’approcher par itérations la valeur au moyen de la
relation suivante :
x(n+1)=x(n)-[f(x(n))/f '(x(n))].
j'attends vos reponses avec patient

[FR119492]

Salut!
Le problème est que la réponse de l'utilisateur sera du texte, c'est-à-dire une chaîne de caractères, ou plutôt deux, car il doit saisir l'expression de ta fonction et de sa dérivée.
J'ai une idée, mais je ne sais pas si ça marche parce que je n'ai jamais essayé: il faudrait que ton programme écrive la fonction dans un fichier fonc.m et sa dérivée dans un fichier deri.m puis que Newton-Raphson fasse appel à ces deux fonctions.
Essaie et dis-nous en des nouvelles.
Jean-Marc Blanc

[phryte]

Bonjour.
Débat sur ce forum :
http://www.developpez.net/forums/d87...ewton-raphson/

[crippella]

merci pour ta réponce mais quelle est la fonction que je vais utiliser pour que l'utilisateur entre une équation est ce que la fonction "input" ou bien une autre??

[crippella]

voila mon essai que je l'ai fait dans un fichier mais lorsque je l'appelais dans MATLAB il ne fonctionne pas portant il me demande d'entrer une équation.
svp vous pouvez me trouver ou est le problème exact.

voila le code:

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
x=-10:0.1:10;
f(x)=input('taper une équation :');% taper l'equation par l'utilisateur
x(1)=input('taper la valeur initiale x(1):\n');
e=1e-10;
n=1000; 
for i=2: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));
         break;
    end
end
j=1:i;
plot(j,x(j),'*r',j,x(j));
xlabel('Nombre d''itérations');
disp('Les valeurs successives de x(i) sont :');

[Jerome Briot]

Code :

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

1
2
 
f(x)=input('taper une équation :');% taper l'equation par l'utilisateur

Essaie plutôt ceci :

Code :

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

1
2
str = input('taper une équation :','s');
f = inline(str);

[crippella]

merci pour la réponse ;
mais quand je fais entrer la fonction et la valeur initiale x(1) il me donne "matix dimensions must agree"et une erreur au niveau de la formule .
svp aider moi pour résoudre ce programme

[phryte]

Bonjour.
Il faut au moins deux valeurs pour diff.

[crippella]

j'ai pas compris !!!!!

[Jean Dumoncel]

Bonjour,

comme te l'a précisé FR119492 dans ta précédente discussion, et il te le redit ici, pour appliquer la méthode de Newton Raphson, il faut connaître la fonction et sa dérivée. Donc tu peux demander à l'utilisateur de renseigner ces 2 fonctions.

Si tu veux utiliser une approximation de la dérivée, comme te l'a précisé phryte, il faut 2 arguments en entrées de diff, et là il faut que tu précises la méthode que tu souhaites utiliser si tu veux qu'on t'aide.

Et sinon le message d'erreur complet que tu obtenais, c'était quoi?

[crippella]

pardon ,
mais mon miniprojet nécessite que l'utilisateur entre une fonction non linéaire seulement et le programme va calculer la dérivée et ou il s'annule par cette methode si vous connaissez une autre fonction qui calcule la dérivée veuillez m'informer car j'ai pas compris comment je vais utiliser deux argument pour calculer la dérivée? vous pouvez m'expliquer?
merci d'avance

[phryte]

Bonjour.
Une piste pour le début (parmi d'autres):

Code :

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

1
2
3
4
5
clear
xi=-10:0.1:10;
str = input('taper une équation :','s');
f = inline(str);
der=diff(f(xi))./diff(xi);

[FR119492]

Salut!

mais mon miniprojet nécessite que l'utilisateur entre une fonction non linéaire seulement

Alors, il ne s'agit pas de la véritable méthode de Newton-Raphson, mais de la méthode appelée "quasi-Newton".
Jean-Marc Blanc

[membreComplexe12]

Salut!

Alors, il ne s'agit pas de la véritable méthode de Newton-Raphson, mais de la méthode appelée "quasi-Newton".
Jean-Marc Blanc

non par forcement,

il dit qu'en entrée il doit rentrer seulement f(x) et que son programme pour pouvoir calculer la formule ci dessous génère automatiquement f'(x) à partir de ce f(x).

x(n+1)=x(n)-[f(x(n))/f '(x(n))].

ps: desolé mais je peux pas t'aider pour l'edition du programme (je suis trop nul en matlab)

[FR119492]

il dit qu'en entrée il doit rentrer seulement f(x) et que son programme pour pouvoir calculer la formule ci dessous génère automatiquement f'(x) à partir de ce f(x).

Alors il a tort: Dans la méthode de Newton-Raphson originale, un utilise l'expression analytique de f'(x) tandis que dans Quasi-Newton, on fait un calcul approché de f'(x) à partir de valeurs de f(x).
Jean-Marc Blanc

[membreComplexe12]

Alors il a tort: Dans la méthode de Newton-Raphson originale, un utilise l'expression analytique de f'(x) tandis que dans Quasi-Newton, on fait un calcul approché de f'(x) à partir de valeurs de f(x).
Jean-Marc Blanc

je suis d'accord avec toi mais à partir de la seule connaissance de f(x) on peut calculer f'(x) et donc apres appliquer Newton-Raphson, c'est se qu'il veut faire si j'ai bien compris....

le seul probleme avec matlab je pense et de lui dire de calculer la derivée d'une fonction f(x) quelconque l'on rentre manuellement...

[crippella]

boujour ,
merci pour vos reponses
j'ai essaié de modifier mon programme d'aprés votre réponces
mais il ya encore des erreurs .............
voici le code:

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
x=-10:0.1:10;
str=input('taper une équation :','s');%taper l'equation par l'utilisateur
f=inline(str);
x(1)=input('entrere la premiere valeur :');
e=1e-10;
n=1000; 
for i=2: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));
         break;
    end
end
j=1:i;
plot(j,x(j),'*r',j,x(j));
xlabel('Nombre d''itérations');
disp('Les valeurs successives de x(i) sont :',x);

merci d'avance

[FR119492]

Bonjour à tous.

à partir de la seule connaissance de f(x) on peut calculer f'(x)

De deux chose l'une: ou bien on connait une expression de f(x) dérivable analytiquement, donc on connait g(x)=f'(x), et on applique Newton-Raphson; ou bien on essaie de calculer approximativement des valeurs discrètes de f'(x) à partir de valeurs discrètes de f(x). Alors, j'attends avec impatience celui qui me proposera un algorithme et/ou un sous-programme accomplissant cette tâche et que je ne sois pas capable de faire dérailler.
Jean-Marc Blanc

[membreComplexe12]

il n'y a pas sous matlab des fonctions de calcul formel qui à partir de f(x) nous donnerai f'(x)?

[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é