Forum de Ti-Gen.org, le site francophone des TI

Bienvenue Invité ! Se connecter Mes sujets	Administrer
Utilisateurs actifs : Aucun membre + 1 visiteur

Avant de poster sur le forum, il y a des régles de bases à respecter pour une bonne entente et un respect de tous.
Veuillez lire la charte du forum.

:: Index » Forum Ti68K » Programmation TI-Basic » fonction de newton (20 réponse(s))

./POST DE DEPART (post n°0) Marquer comme non lu.

collinm

Ecrit le: Vendredi 15 juillet 2005 à 04:34

Membre: Visiteur

Messages: 19

Membre n°899

Inscrit le 01 juillet 2005

salut

jais le code suivant dans ma ti:

newton(a,to,trop)
func
local m,xn
0->m:0->xn
loop
    m+1->m
    if(d(f(x),x)|x=a)=0 then @@@@@
        exit
    else
        a-f(a)/d(f(x),x)|x=a)->xn
        if m>= trop or abs(xn-a)<tol
            exit
        endif
    xn->a
endloop
if (d(f(x),x|x=a)=0 then
    "la methode echoue"
else
    {m,xn}
endif
endfunc

si je fais:

newton(3*x,x,4)

j'obtient

a test did not resolve to true or false

et ça m'amène à la ligne @@@@@

une idée?

www.laboiteaprog.com

./Post n°1 Marquer comme non lu.

Kevin Kofler

Ecrit le: Vendredi 15 juillet 2005 à 04:47

Administrateur

Groupe: Staff de Ti-Gen

Messages: 4814

Membre n°40

Inscrit le 08 février 2004

f(x) n'est défini nulle part dans ton code.

Membre de l'équipe de TIGCC: http://tigcc.ticalc.org
Mainteneur du portage Linux/Unix de TIGCC: http://tigcc.ticalc.org/linux/
Membre de l'équipe de CalcForge: http://www.calcforge.org:70/

Participez à la reprise de Ti-Gen!

./Post n°2 Marquer comme non lu.

collinm

Ecrit le: Vendredi 15 juillet 2005 à 05:39

Membre: Visiteur

Messages: 19

Membre n°899

Inscrit le 01 juillet 2005

je remplace x par a

www.laboiteaprog.com

./Post n°3 Marquer comme non lu.

RHJPP

Ecrit le: Vendredi 15 juillet 2005 à 11:15

Administrateur

Groupe: Staff de Ti-Gen

Groupe: Omega staff

Messages: 2174

Membre n°776

Inscrit le 06 avril 2005

Non c'est la fonction f qui n'existe pas

./Post n°4 Marquer comme non lu.

collinm

Ecrit le: Samedi 16 juillet 2005 à 00:15

Membre: Visiteur

Messages: 19

Membre n°899

Inscrit le 01 juillet 2005

en gros je veux faire comme: http://lumimath.univ-mrs.fr/~jlm/cours/analnum/node3.html

comme paramètre f(x), x0 (approximation initiale), n n'ombre d'itération

newton(f(x),x0,n)
func
local m,xn
0->m:0->xn
loop
     m+1->m
     if d(f(x),x)=0 then
         exit
     else
         x-f(x)/d(f(x),x)->xn
         if m>=n trop or abs(xn-a)<x0
             exit
         endif
     xn->a
endloop
if  d(f(x),x)=0 then
     "la methode echoue"
else
     {m,xn}
endif
endfunc

j'ai l'erreur:

Error: Argument must be a variable name

www.laboiteaprog.com

./Post n°5 Marquer comme non lu.

Folco

Ecrit le: Samedi 16 juillet 2005 à 00:47

Administrateur

Groupe: Staff de Ti-Gen

Messages: 1286

Membre n°181

Inscrit le 05 avril 2004

à quelle ligne?
edit-> commence à ce moment là par détruire tes variables de a à z.
-Edité le Samedi 16 juillet 2005 à 00:47 par Martial Demolins-

<<< Kernel Extremist©®™ >>>
Pas la peine d'aller là plus d'une fois tous les six mois...

"Il faut apprendre pour savoir qu'il faut apprendre pour savoir."

./Post n°6 Marquer comme non lu.

collinm

Ecrit le: Samedi 16 juillet 2005 à 07:01

Membre: Visiteur

Messages: 19

Membre n°899

Inscrit le 01 juillet 2005

si on regarde sur: http://www.bibmath.net/dico/index.php3?action=affiche&quoi=./n/newton_meth.html

ça serait la première formule (on le voit par le dessin)

c'est déjà fait

www.laboiteaprog.com

./Post n°7 Marquer comme non lu.

RHJPP

Ecrit le: Samedi 16 juillet 2005 à 11:28

Administrateur

Groupe: Staff de Ti-Gen

Groupe: Omega staff

Messages: 2174

Membre n°776

Inscrit le 06 avril 2005

"f(x)" n'est pas un nom de variable en suite "d(f(x),x)=0" ne peut pas fonctionner.

./Post n°8 Marquer comme non lu.

collinm

Ecrit le: Samedi 16 juillet 2005 à 23:47

Membre: Visiteur

Messages: 19

Membre n°899

Inscrit le 01 juillet 2005

newton(3*x^2-x,1,3

newton(a,to,trop)
func
local f, m,xn
0->m:0->xn
a->f(x)
loop
     m+1->m
     if(d(a,x)=0 then
     exit
....

j'obtiens une erreur à la dérivé

a test did not resolve to true or false

www.laboiteaprog.com

./Post n°9 Marquer comme non lu.

RHJPP

Ecrit le: Dimanche 17 juillet 2005 à 00:00

Administrateur

Groupe: Staff de Ti-Gen

Groupe: Omega staff

Messages: 2174

Membre n°776

Inscrit le 06 avril 2005

Normal, la dérivée de ta fonction dépend de x donc la calculatrice ne sait pas si ça fait 0 ou pas de plus il doit manquer ")".
Dans ta boucle rien ne dépend de m donc elle risque de ne jamais se terminer.

-Edité le Dimanche 17 juillet 2005 à 00:04 par Thepro-

./Post n°10 Marquer comme non lu.

collinm

Ecrit le: Dimanche 17 juillet 2005 à 00:21

Membre: Visiteur

Messages: 19

Membre n°899

Inscrit le 01 juillet 2005

newton(a,x0,n)
  func
  local f, m,xn
  0->m:0->xn
  a->f(x)
  loop
      m+1->m
      if(d(a,x)=0 then
      exit
...

idem

www.laboiteaprog.com

./Post n°11 Marquer comme non lu.

RHJPP

Ecrit le: Dimanche 17 juillet 2005 à 00:51

Administrateur

Groupe: Staff de Ti-Gen

Groupe: Omega staff

Messages: 2174

Membre n°776

Inscrit le 06 avril 2005

Regarde

"if(d(a,x)=0 then" ne peut pas fonctionner d(a,x)=a'(x) il y a toujours x qui est en trop
dans ton problème on utilise x=abscisse de la tangente (à l'ancien x) lorsque cette tangente coupe l'axe.
Il faut donc calculer cette tangent et son point d'intersection avec l'axe des abscisses.

./Post n°12 Marquer comme non lu.

collinm

Ecrit le: Dimanche 17 juillet 2005 à 02:39

Membre: Visiteur

Messages: 19

Membre n°899

Inscrit le 01 juillet 2005

pour le moment j'ai enlever le calcul de la dérivé 0... étant donné que ça cause problème...

j'ai fait un programme et la méthode de newton fonctionne j'ai fait la même chose avec une fonction.... et ça passe pu

le programme

ttt()
prgm
local f,m,xn,uu,n,x0
clrio
input "valeur ini", x0
x0->uu
1->
3->n
x^2-2->f(x)
loop
    uu-f(uu)/d(f(uu),uu))->xn
    if m>=n then
        exit
    endif
    xn->uu
    m+1->m
endloop
disp "nb ite: "&string(m),"
u(n)="&string(uu)
endprgm

j'obtient si je met 1 pour la valeur initiale:

nb ite: 3
u(n)=17/12

avec le programme

newton(x^2-2,1,3)

newton(a,x0,n)
func
clrio
local f,m,xn,uu
1->m
a->f(x)
x0->uu
loop
    uu-f(uu)/d(f(uu),uu))->xn
    if m>=n then
        exit
    endif
    xn->uu
    m+1->m
endloop
return uu
endfunc

j'obtient underf

une idée

www.laboiteaprog.com

./Post n°13 Marquer comme non lu.

Kevin Kofler

Ecrit le: Dimanche 17 juillet 2005 à 04:37

Administrateur

Groupe: Staff de Ti-Gen

Messages: 4814

Membre n°40

Inscrit le 08 février 2004

a->f(x) ne fonctionnera pas, tu dois écrire expr(string(a)&"->f(x)").

./Post n°14 Marquer comme non lu.

RHJPP

Ecrit le: Dimanche 17 juillet 2005 à 13:51

Administrateur

Groupe: Staff de Ti-Gen

Groupe: Omega staff

Messages: 2174

Membre n°776

Inscrit le 06 avril 2005

Voilà ce que je ferais :

newton(a,xn,n)
Func
Local d,m
d(a,x)->d
0->m
While m<n
xn-a/d|x=xn->xn
m+1->m
EndWhile
Return xn
EndFunc

./Post n°15 Marquer comme non lu.

limmt

Ecrit le: Dimanche 17 juillet 2005 à 19:19

Administrateur

Groupe: Staff de Ti-Gen

Groupe: Orage Studio

Groupe: Omega staff

Messages: 4712

Membre n°20

Inscrit le 01 février 2004

pour ton problème de if qui marche pas il faut mettre:
if d(a,x)=0 then
et non pas:
if(d(a,x)=0 then

http://www.falco-fr.com/ - http://www.jump67.com/ - http://www.msf-league.com/

./Post n°16 Marquer comme non lu.

RHJPP

Ecrit le: Dimanche 17 juillet 2005 à 19:24

Administrateur

Groupe: Staff de Ti-Gen

Groupe: Omega staff

Messages: 2174

Membre n°776

Inscrit le 06 avril 2005

En fait limmt ya pas que ça comme problème
d(a,x) varie en fonction de x (généralement) donc on ne peut pas savoir si ça vaut 0

./Post n°17 Marquer comme non lu.

limmt

Ecrit le: Dimanche 17 juillet 2005 à 22:32

Administrateur

Groupe: Staff de Ti-Gen

Groupe: Orage Studio

Groupe: Omega staff

Messages: 4712

Membre n°20

Inscrit le 01 février 2004

d(a,x) renverra une fonction qui peut etre nulle mais tu peut pas savoir si ca ser aégal a 0 apsque ca accepte aps la condition, tu peut tester avec des valeurs discrete smasi bon c'est aps bien comme méthode

http://www.falco-fr.com/ - http://www.jump67.com/ - http://www.msf-league.com/

./Post n°18 Marquer comme non lu.

Link

Ecrit le: Dimanche 17 juillet 2005 à 23:16

Membre: Magister

Messages: 302

Membre n°896

Inscrit le 25 juin 2005

enfin, il peut toujours tester avec un when: si la fonction dépend de x, elle ne sera donc pas égale à zéro, donc:

if when(d(a, x)=0, true, false, false) : exit

pourrait marcher.
Attention par contre, je ne vois pas trop si c'est la seule condition pour quitter la boucle, mais j'espère bien que non, car la dérivée d'un sinus ou d'un cos ne sera jamais nulle...