Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter
Méthode de Runge-Kutta
Voici deux posts avec deux implémentations (C++ et Java) de la méthode de Runge-Kutta :
C++ :
http://www.developpez.net/forums/sho...77&postcount=2
Java :
http://www.developpez.net/forums/sho....php?p=2528269
Ces posts datent, je viens seulement de penser à les poster ici.
A noter que j'ai proposé une solution pour résoudre une équation de la forme :
dy/dx = f(x,y) (je vous passe les détails sur les hypothèses à faire sur les fonctions)
cool...
On va pouvoir poster des filtres a base de PDE...
ALGORITHME (n.m.): Méthode complexe de résolution d'un problème simple.
Méthode de Runge-Kutta
Pour tous ceux qui cherchent des algorithmes performants pour le calcul numérique, je ne peux que recommander "Numerical Recipes" déjà souvent cité sur ce forum, et plus encore le site www.netlib.org .
En ce qui concerne l'intégration des équations différentielles ordinaires, il n'y a pas que Runge-Kutta, et je vous recommande de faire un tour sur www.unige.ch/~hairer/software.html
Bonne chance à tous.
Jean-Marc Blanc
Bonjour,
je viens d'essayer le lien java qui ne marche pas
Consignes aux jeunes padawans : une image vaut 1000 mots !
- Dans ton message respecter tu dois : les règles de rédaction et du forum, prévisualiser, relire et corriger TOUTES les FAUTES (frappes, sms, d'aurteaugrafe, mettre les ACCENTS et les BALISES) => ECRIRE clairement et en Français tu DOIS.
- Le côté obscur je sens dans le MP => Tous tes MPs je détruirai et la réponse tu n'auras si en privé tu veux que je t'enseigne.(Lis donc ceci)
-ton poste tu dois marquer quand la bonne réponse tu as obtenu.
GrrrrrrrrrEnvoyé par ToTo13
Bonjour,
je viens d'essayer le lien java qui ne marche pas
Je sais pourquoi... J'ai une sauvegarde quelque part, je remettrais quand je pourrais.
Méthode de Runge-Kutta
Si quelqu'un est intéressé à l'utiliser en Fortran, j'ai retrouvé ça dans mon répertoire à vieilleries:
Subroutine D006(SM,N,T,DT,Y)
C
C***********************************************************************
C
C Bibliothèque JMBPACK
C
C Série D: Equations différentielles
C
C***********************************************************************
C
C Méthode de Runge-Kutta d'ordre 4 (1/6 + 2/6 + 2/6 + 1/6)
C
C***********************************************************************
C
C Version 3.0: Jean-Marc Blanc, juillet 2004
C
C **********************************************************************
C
C SM Nom d'un sous-programme fourni par l'utilisateur pour le
C calcul des dérivées à partir des valeurs de la variable
C indépendante et des variables d'état. Ce sous-programme
C doit être de la forme
C Subroutine SM(N,T,Y,DY)
C DY(1)= ...
C DY(2)= ...
C .
C .
C DY(N)= ...
C Return
C End
C et avoir été déclaré dans une instruction External.
C
C N Ordre du système différentiel.
C
C T Variable indépendante.
C
C DT Pas d'intégration.
C
C Y Vecteur des variables d'état.
C
C **********************************************************************
C
Implicit None
C
Integer N
Real*8 DT,T,Y(N)
C
Integer I
Real*8 DY(N),K1(N),K2(N),K3(N),K4(N),Z(N)
C
Call SM(N,T,Y,DY)
Do I=1,N
K1(I)=DT*DY(I)
Z(I)=Y(I)+K1(I)/2.d0
End Do
C
Call SM(N,T+DT/2.d0,Z,DY)
Do I=1,N
K2(I)=DT*DY(I)
Z(I)=Y(I)+K2(I)/2.d0
End Do
C
Call SM(N,T+DT/2.d0,Z,DY)
Do I=1,N
K3(I)=DT*DY(I)
Z(I)=Y(I)+K3(I)
End Do
C
Call SM(N,T+DT,Z,DY)
Do I=1,N
K4(I)=DT*DY(I)
Y(I)=Y(I)+(K1(I)+2.d0*K2(I)+2.d0*K3(I)+K4(I))/6.d0
End Do
C
T=T+DT
Return
C
End
Salut à tous
Jean-Marc Blanc
On a réussi à recupérer les données :Bonjour,
je viens d'essayer le lien java qui ne marche pas
http://www.developpez.net/forums/sho....php?p=2528269
Répondre avec citation
Partager