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 » conversion de programme ti83 à ti92 (5 réponse(s))

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

collinm

Ecrit le: Jeudi 4 août 2005 à 01:01

Membre: Visiteur

Messages: 19

Membre n°899

Inscrit le 01 juillet 2005

salut

j'ai trouvé une fonction que j'ai besoin: jacobi (les itérations...)
le hic c'est qu'elle est faite pour les ti83

voici le code:

Instruction#1_____If [A](1,1) = 0..........(See note at.end of program.)
Instruction#2_____Go to 1
Instruction#3_____dim([A]) -> L1..........(Note: We use the symbol "->" to
Instruction#4_____L1 -> dim([D])...........represent the operation "Store".)
Instruction#5_____Fill(0,[D])
Instruction#6_____L1(1) -> N
Instruction#7_____For(I,1,N)
Instruction#8_____[A](I,I) -> [D](I,I)
Instruction#9_____End
Instruction#10____[D] - [A] -> [E]
Instruction#11____[D]-1[E] -> [F]
Instruction#12____[D]-1[B] -> [G]
Instruction#13____[A](1,1) -> S
Instruction#14____0 -> [A](1,1)
Instruction#15____Lbl 1
Instruction#16____[F][C] + [G] -> [H]
Instruction#17____[H] - [C] -> [I]
Instruction#18____[H] -> [C]
Instruction#19____augment([H],[I]) -> [J]

comment convertir cela pour une ti-92

j'ai trouvé la fonction sur la page: http://qcpages.qc.cuny.edu/~roskes/NumAnalTI83.htm#Program #11

merci

www.laboiteaprog.com

./Post n°1 Marquer comme non lu.

RHJPP

Ecrit le: Jeudi 4 août 2005 à 01:20

Administrateur

Groupe: Staff de Ti-Gen

Groupe: Omega staff

Messages: 2174

Membre n°776

Inscrit le 06 avril 2005

Bizarre !? [c] vient d'où ?

./Post n°2 Marquer comme non lu.

supprimé000

Ecrit le: Jeudi 4 août 2005 à 03:29

Membre: Newbie

Messages: 0

Membre n°901

Inscrit le 01 janvier 1970

Soit il n'est pas complet ton code soit tu n'as pas tout copié(présence de [c]) et en plus ton lien est mort.
Normalement si tu comprends ce code ça doit pas être difficile de le reproduire pour TI-92.
Poste le code complet et ça sera plus claire.

./Post n°3 Marquer comme non lu.

Lionel Debroux

Ecrit le: Jeudi 4 août 2005 à 17:21

Envoyer un message privé à Lionel Debroux

Membre: Maître Ultime

Messages: 1721

Membre n°10

Inscrit le 31 janvier 2004

N'y a-t-il pas de Jacobi dans MathTools ?

Lionel Debroux - membre de TICT.

./Post n°4 Marquer comme non lu.

collinm

Ecrit le: Vendredi 5 août 2005 à 05:56

Membre: Visiteur

Messages: 19

Membre n°899

Inscrit le 01 juillet 2005

http://qcpages.qc.cuny.edu/~roskes/NumAnalTI83.htm

le code est bien complet... mais j'ai peut-être trouver quelque chose de plus intéressant...

salut

je suis en train de porter cette fonction provenant de mathlab... sur la ti-92

et j'aurais un peu besoin d'aide...

// Résolution de système du type A.X=B
// par la méthode itérative de JACOBI
// Initialisation
A=[-16 6 -2 -5;
      3 10 -5 1;
     -4 1 18 2;
      1 2 2 -14];
B=[-19;
      1;
     12;
      1];
X0=[0.1;
     0.1;
     0.1;
     0.1];
Xk=X0;
iter=0;
max_it=500;
tol = 0.0000000000001;

// Conditionnement des matrices D, L et U
n=4;
for i=1:n,
          for j=1:n,
                    I(i,i)=1;
                    D(i,i)=A(i,i);
                    if i>j then, L(i,j)=-A(i,j);,
                    else L(i,j)=0;, end,
                    if i<j then, U(i,j)=-A(i,j);,
                    else U(i,j)=0;, end ;
end ; end;

// première itération
// JACOBI
Xkplus1 = ( inv(D)*(L+U) * Xk + inv(D)* B);
// itérations suivantes
while norm(Xkplus1 - Xk)/norm(Xk) > tol,
      iter=iter+1;
      Xk = Xkplus1;
      // JACOBI
      Xkplus1 = ( inv(D)*(L+U) * Xk + inv(D)* B);
end;
// Solution obtenue
iter
Xkplus1

sur la ti ça me donne jusqu'à maintenant

jacobi(a,b,x0)
func
local xk,item,max_it,tol,i,j,n,d,u,l
coldim(a)->n
x0->xk
0->ite
500->max_it
.000000000001->tol

for i,1,n,1
  for j,1,n,1
    a[i,j]->d[i,i]
    if i>j then
      -a[i,j]->l[i,j]
    else
      0->l[i,j]
     endif

     if i<j then
      -a[i,j]->u[i,j]
    else
      0->u[i,j]
     endif
   endfor
endfor
xkplus1=d^(-1)*(l+u)*xk+d^(-1)*b
while norm(xkplus1-xk)/norm(xk) > tol
  iter+1>iter
  xkplus1->xk
  xkplus1=d^(-1)*(l+u)*xk+d^(-1)*b
endwhile

endfunc

lorsque je fais exécuter ce code
on me ramène a:
a[i,j]->d[i,i]

et on dit:
error: dimension

merci

www.laboiteaprog.com

./Post n°5 Marquer comme non lu.

Managaia

Ecrit le: Dimanche 28 août 2005 à 23:01