Equa diff pour Matlab
Réponses à toutes vos questions après le Bac (Fac, Prépa, etc.)
-
kimimaro
- Messages: 9
- Enregistré le: 24 Mar 2007, 08:07
-
par kimimaro » 24 Mar 2007, 08:20
Je dois transformer le système suivant en système d'équa diff d'ordre 1 afin d'effectuer une résolution numérique avec la fonction ode45 de Matlab :
}=y"(t) e^{t}+\left(y^{(3)}(t)\right)^3)
et y(1)=0
et y'(1)=2
Problème : les conditions initiales semblent insuffisantes
S'agirait-il d'une erreur de l'énoncé où y a-t-il un moyen de transformer l'équa diff pour qu'elle ne concerne plus que y, y' et y" par exemple ?
-
Flodelarab
- Membre Légendaire
- Messages: 6574
- Enregistré le: 29 Juil 2006, 14:04
-
par Flodelarab » 24 Mar 2007, 13:56
kimimaro a écrit:Je dois transformer le système suivant en système d'équa diff d'ordre 1 afin d'effectuer une résolution numérique avec la fonction ode45 de Matlab :
}=y"(t) e^{t}+\left(y^{(3)}(t)\right)^3)
et y(1)=0
et y'(1)=2
Problème : les conditions initiales semblent insuffisantes
S'agirait-il d'une erreur de l'énoncé où y a-t-il un moyen de transformer l'équa diff pour qu'elle ne concerne plus que y, y' et y" par exemple ?
Je vais ptet dire une bétise ... mais:
Et si tu posais Y(t)=y''(t) ???
n'aurais tu pas Y''(t)=Y(t)exp(t) + (Y'(t))³ ???
-
kimimaro
- Messages: 9
- Enregistré le: 24 Mar 2007, 08:07
-
par kimimaro » 24 Mar 2007, 16:15
Flodelarab a écrit:Je vais ptet dire une bétise ... mais:
Et si tu posais Y(t)=y''(t) ???
n'aurais tu pas Y''(t)=Y(t)exp(t) + (Y'(t))³ ???
Oui, mais c'est les conditions y(1)=0 et y'(1)=0 dont je ne sais pas quoi faire.
Si je suis ta proposition, la mise en forme pour Matlab devient :
\\<br />Y'(t)<br />\end{tabular}\right) '<br />=<br />\left(\begin{tabular}<br />Y'(t)\\<br />Y(t) e^{t} + \left(Y'(t)\right)^3<br />\end{tabular}\right)<br />\,<br />ie<br />\,<br />dZ<br />=<br />\left(\begin{tabular}<br />Z_2\\<br />Z_1 e^{t}+{Z_2}^3<br />\end{tabular}\right))
La fonction ode45 prend comme paramètres :
- la fonction calculant dZ en fonction de t et Z
- l'intervalle d'intégration, ici [1,b] où b est à déterminer visuellement
- Y(1) et Y'(1), ie la valeur initiale de Z, que je ne connais pas
-
nox
- Membre Complexe
- Messages: 2157
- Enregistré le: 14 Juin 2006, 09:32
-
par nox » 24 Mar 2007, 18:08
kimimaro a écrit:Y(0) et Y'(0), ie la valeur initiale de Z, que je ne connais pas
ba prends [1 B] comme intervalle d'intégration
évidemment si comme condition initiale on te donne y(1) ca sous-entend que 1 est le moment initial.
Donc si on te demande [0 B] je dirais que c'est une erreur d'énoncé
-
kimimaro
- Messages: 9
- Enregistré le: 24 Mar 2007, 08:07
-
par kimimaro » 24 Mar 2007, 18:18
En fait c'est moi qui ai fait une petite erreur dans mon message, j'aurais dû écrire :
La fonction ode45 prend comme paramètres :
- la fonction calculant dZ en fonction de t et Z
- l'intervalle d'intégration, ici [1,b] où b est à déterminer visuellement
- Y(1) et Y'(1), ie la valeur initiale de Z, que je ne connais pas
Désolé pour la confusion, nox.
Mon problème n'est pas sur les bornes de l'intervalle, mais sur le fait que les conditions limites font intervenir y' et y (je précise y 'minuscule', pas Y 'majuscule' qui est le changement proposé par FLodelarab) qui n'apparaissent pas dans l'équa diff.
-
fahr451
- Membre Transcendant
- Messages: 5142
- Enregistré le: 05 Déc 2006, 23:50
-
par fahr451 » 24 Mar 2007, 18:30
bonsoir
qui est ode (45)?
-
kimimaro
- Messages: 9
- Enregistré le: 24 Mar 2007, 08:07
-
par kimimaro » 24 Mar 2007, 18:38
fahr451 a écrit:bonsoir
qui est ode (45)?
ode45 est une fonction de Matlab permettant de résoudre numériquement une équa diff (ODE : Ordinary Differential Equation).
Voici ce que dit la doc sur l'algo utilisé :
ode45 is based on an explicit Runge-Kutta (4,5) formula, the Dormand-Prince pair.
-
fahr451
- Membre Transcendant
- Messages: 5142
- Enregistré le: 05 Déc 2006, 23:50
-
par fahr451 » 24 Mar 2007, 19:11
ah c'est runge kutta qui a changé de papiers
merci
Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 21 invités