Discussion :

[initil]

Bonjour,
Je cherche a résoudre un système de 6 équations à 6 inconnues du 1er ordre que j'ai définis a partir de la méthode de runge-kutta.
Voici le systeme:
---------
α1’= K2/J1*θ1- (K2/J1)* θ2 + (Cm-Cf1)/J1;
α2’= K2/J2* θ1 - (K2+K3)/J2* θ2 + (K3/J2)* θ3 - (Cf2/J2);
α3’ = K3/J3* θ1 - (-K3/J3)* θ2 - (Cf3+Cr)/J3;
θ1’= α1
θ2’= α2
θ3’= α3
---------
K1...K3, J1...J3, Cm,Cf1....Cf3 et Cr sont des constantes.(ne faite pas attention).
θ1,θ2,θ3 , α1,α2,α3 sont des variables que je recherche. Elles représentent en fait l'angle de position et les vitesses des arbres couplés a un moteur.

J'aimerai résoudre ce systeme grace a l'ode45 de matlab en passant par une matrice M(t, y) y' = F(t, y) dans mon odeset ('Mass',@mamatrice)
Comment défnir cette matrice ??
Je sais poser le systeme sous la forme [A].[y]+[B]=[y'], mon M(t,y)y' serait alors -[B]+[y'] ??

[christophe_halgand]

Si j'ai bien compris :

Code :

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

y=[a1 a2 a3 θ1 θ2 θ3]';

Dans ce cas, il est simple de trouver ta matrice A et b, non ?

Christophe

[initil]

Oui c'est bien ça y.
J'ai déja mes matrice A et B mais la fonction odeset me demande une matrice M(t). sachant que M(t, y) y' = F(t, y)