Discussion :

[TiMeWaLk]

Bonjour,

J'ai une équation de la forme du/dz=H(z)*u(z) où u est un vecteur ayant environ 4000 coefficients, H une énorme matrice et u(0)=u0 un vecteur connu.

J'ai vu dans l'aide de matlab que ce genre d'équation nécessite la création d'une fonction "handle" (par exemple pour ode45), le problème est le suivant: H est une matrice que j'ai crée coefficients par coefficients, pour différents z, en gros j'ai accès à H(z1), H(z2), .... H(z100).

Comment puis-je résoudre cette équation avec matlab et mes données?

[kmaniche]

http://www.developpez.net/forums/d68...fferentielles/

[TiMeWaLk]

Visiblement j'avais mal lu la documentation qui semble résoudre mon problème, même si le temps de calculation sera de plusieurs heures...

Merci

[kmaniche]

qu'as tu fais ?

[TiMeWaLk]

Pour le moment, j'ai pas encore lancer le programme mais je pense faire quelquchose comme ça:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
function dy = rigid(t,y)
dy = zeros(3,1);    % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);

Sauf que dy va de 1 à 4000, rigid sera définit comme ceci:

dy(1)=H*y(1), le problème vient du calcul de H à chaque évaluation de dy...

H est une matrice dont chaque coefficient est un produit de 3 facteurs, dont des coefficients de séries de fourier...