Discussion :

[dptmt]

Bonjour,

J'essaye de résoudre ce système d'équations différentielles avec la fonction ODE de matlab :

dH1(omega,X0)/dX0 = Y1/(a1*X0)-H1(omega,X0)/X0-((omega*b1)/(X0^2*a1))*H2(omega,X0)

dH1(omega,X0)/dX0 = Y1/(conj(a1)*X0)-H1(omega,X0)/X0+((omega*conj(b1))/(X0^2*conj(a1)))*H2(omega,X0)

les inconnues sont H1 et H2 (matrices de tailles M*N)
Y1, a1 , b1 , omega et X0 sont tous connues.
le problème c'est que Y1 , a1, b1 sont des matrices de tailles M*N
omega et X0 sont des vecteurs de tailles 1*M et 1*N respectivement (sachant qu'on peut les dupliquer pour en former des matrices).
on intègre de X0min à X0max.

J'arrive pas à résoudre ce système à cause la dimension matricielle du problème.
Si vous pouvez me donner quelques pistes ça serait gentille

[FR119492]

Salut!
Comme cela arrive très souvent sur ce site, tu es confronté à un problème (probablement de simulation d'un phénomène physique) et tu as essayé de le résoudre en commençant par le modéliser, c'est-à-dire le mettre en équations. Malheureusement, tu t'y es mal pris, ce qui t'a conduit à des équations qui n'ont pas de sens. En plus du fait que ta seconde équation devrait vraisemblablement commencer par dH2(... , j'y vois plusieurs objections:

• dH1/dX0: que représente la dérivée d'une matrice par un vecteur?
• H1(omega,X0): quelle est l'expression de cette fonction?
• H1(omega,X0)/X0: que représente le quotient d'une matrice par un vecteur?
• X0^2: comment calcules-tu le produit d'un vecteur par lui-même?
• conj(a1): qu'est-ce que des coefficients complexes viennent faire dans un système différentiel?

Comme tu le vois, il y a encore beaucoup de points d'interrogation. Pour que nous puissions t'aider, il me semble indispensable que tu nous expose ton problème depuis le début (le phénomène physique), et, comme cela n'a rien à voir avec Matlab, je transfère cette discussion dans le forum algo/maths.
Jean-Marc Blanc