Equations diferentielle du 2è degre

[Fermé]
Signaler
Messages postés
7
Date d'inscription
dimanche 14 janvier 2007
Statut
Membre
Dernière intervention
19 janvier 2007
-
 Ben Salem Manel -
Je voudrais savoir comment prg -t-on sous matlab pour resoudre les equation differentielles avec seconde membre avec la methode des elements finis ! qelq'un peut m'aider ?

7 réponses

Messages postés
110
Date d'inscription
mardi 12 avril 2005
Statut
Membre
Dernière intervention
17 avril 2012
28
bonjour,

ODE45 est une méthode pour équation différentielles de Matlab.
Une résolution Runge-Kutta ou approchant si je me souviens bien.
Il faut donner les différents coefficients et leux variations ensuite apres plusieurs itérations, tu peux ploter le résultats.

J'espère t'avoir aidé un peu.
Bonne chance ;)
Messages postés
110
Date d'inscription
mardi 12 avril 2005
Statut
Membre
Dernière intervention
17 avril 2012
28
Alors, je devrais avoir un bout de code qui traine quelque part.
Je vais essayer de te trouver ca.

Laisse-moi deux ou trois jours histoire de partir en quete du fichier perdu.
Et ca risque d'être une grande expédition :p

Je te tiens au courant
Messages postés
7
Date d'inscription
dimanche 14 janvier 2007
Statut
Membre
Dernière intervention
19 janvier 2007
>
Messages postés
110
Date d'inscription
mardi 12 avril 2005
Statut
Membre
Dernière intervention
17 avril 2012

haha merci Khulun bcp , j'ai aussi fait des recherche sur net et j'aimerais savoir comment on peut resoudre l'eq (1+x)Y''+Y'=-x en methodes des elements finis ? mais pas en methde Rung-Kutta !?
Messages postés
110
Date d'inscription
mardi 12 avril 2005
Statut
Membre
Dernière intervention
17 avril 2012
28 >
Messages postés
7
Date d'inscription
dimanche 14 janvier 2007
Statut
Membre
Dernière intervention
19 janvier 2007

Il n'y a pas de Y(x) dans ton équation.
A moins que je dises une grosse grosse bétise, tu peux faire un changement de variable, Y' = K

Ton équation te donne (1+x)K' + K = -x
Et tu te retrouve avec une equation différentielle de premier ordre.
Tu trouve la solution (facile ca ^^) et pis tu intègre la solution trouvée (car k=Y' et pas Y).

Pour la résolution, une méthode d'Euler toute bete et hop :)
Le secret de l'univers sera dans tes mains.
Messages postés
7
Date d'inscription
dimanche 14 janvier 2007
Statut
Membre
Dernière intervention
19 janvier 2007
>
Messages postés
110
Date d'inscription
mardi 12 avril 2005
Statut
Membre
Dernière intervention
17 avril 2012

j'ai l'ai vu le pb mais puisque je voudrais uiliser le met des elements finis pour le resoudre mais pas du methode d'Euler ou Runge-Kutta ! Est ce tu connais cette methode ou bien methode de Ritz ?
>
Messages postés
110
Date d'inscription
mardi 12 avril 2005
Statut
Membre
Dernière intervention
17 avril 2012

Si tu veux la méthode de Ritz en SCILAB:
http://www.arbredeslemuriens.com/Categorie.php?IDCategorie=AlgoScilab&IDTitre=146
bonjour;
j'ai un probleme pour resoudre sous matlab cette equation différentielle:

d²(d)/dt²+((4.7+0.75/(d(x)/dt+d(y)/dt)^0.5)*d(d)/dt=Fd(t)

Avec :
d : c’est x ou y
Fd(t) : c’est 1 ou bien -1
On donc 2 équations différentielle, l'une en fonction de ' x 'et l'autre en fonction de ' y ' a résoudre ?
aller au site : www.les-mathmathique.net
slt est merci pour votre aide j'ai besoin les étapes pour la création du programme d'une equation du 2 éme degré.
mon amie moi aussi je cherche queleque cours sur la mèthode des èlèment fins parceque mon projet de fin d'etude baser sur cette mèthode si je trouve quelque document sur cette mèthode je envoyer a toi merci et bonne chance.
salut , bonjour a tous, je veux résoudre une équation différentielle de premier ordre en matlab par l'utilisation de la fonction "ode45" , et malheureusement j'ai trouvé des problèmes sur cette fonction ,j'ai essayé de respecter tous les conditions de fonction de "ode45" et j'ai pas trouvé la faute ou exacte , alors sa me fais plaisir vraiment quant vous allez m'aider SVP et merci en tout cas.
voici le programme que je le construire:
%%%% LA FONCTION%%%%%%%%%
function dphi=rsprime(t,phi)
Qx=0.5;
phi=pi/12;
gamma1=0.08;
kphi=3;
dphi=gamma1+kphi*Qx
%%%%%%%%%%%%%%%
clc;
clear;
Qx=0.5;
phi0=pi/12;
gamma1=0.08;
t0=0;
tf=10;
kphi=3;
[t,phi]=ode45(rsprime,t0,tf,phi0);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
et l'erreur qui je le trouve c'est
" Error using ==> exist"
The first input to exist is a string.
et "error in ==> funfun\private\odearguments at 79"
if (exist(ode)==2) % M-file
et"Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ...

Error in ==> exemple2ode45 at 9
[t,phi]=ode45(rsprime,t0,tf,phi0);
aide moi SVP

Organigramme de runge kutta d ordre 4
Données Y0, n, f(y,t)
Dt=(t2-t1)/2n
i = 0
t0 = 0
((test)) i inférieur ou égal à n
si non : écrire y0,....., yn
stop
si oui : calcul de y chapeau (inconnue) i+1/2 =.....
calcul de y double chapeau =....
Dt = (t2-t1)/n
i = i+1/2
CALCUL DE y chapeau i+1 =.....
Dt/6 = (t2-t1)/6
calcul de y i+1+....