Discussion :

[galaadb]

Bonjour,

je cherche à résoudre le système suivant: dG=MG.
G est le vecteur d'état et contient deux sous-fonctions.
M est une matrice 2x2 dont les éléments dépendent de x. M est définie via une autre fonction, et dépend également de w qu'on prend constant.
x est une liste dont les éléments correspondent à ceux de xspan.

J'essaie d'utiliser la fonction ode45 pour y parvenir:

Code :

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

1
2
 
[x,G]=ode45(@green,xspan,[0 1],[],M);

avec la fonction green.m qui définit le système:

Code :

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

1
2
3
4
5
6
7
 
function dG = green(x,G,w,M)
 
M=mat_sys(x,w);
 
dG = zeros(2,1);
dG=M*G;

Pour l'instant, Matlab me signale une erreur de dimension:

Error using ==> mtimes
Inner matrix dimensions must agree.

Est-ce que quelqu'un aurait une idée de la manière de résoudre ce genre de système?
Merci d'avance.

[captain_usopp]

Je ne sais pas résoudre ton problème mais tu peux pour t'y retrouver utiliser la fonction size qui te renvoie la taille d'une matrice.
Cela peut t'aider à vérifier ce que tu fais et si tes matrices sont bien de dimension voulue.
Désolé si tu connais cette fonction et si je ne te suis d'aucune utilité mais on ne sait jamais. D'autant que personne ne t'a encore répondu.

[galaadb]

Merci d'avoir jeté un coup d'oeil à ma demande.
Le problème c'est que dans le workspace ma matrice est de taille 2x2, donc ça devrait coller... sauf que ça ne prend pas en compte la dépendance en x...
Je suis un peu coincée là...

[Caro-Line]

As-tu bien regardé la taille de M et de G ?
M est 2x2
et G ?
Ne serait-il pas 1x2 ?
Si oui, ne faudrait-il pas qu'il soit plutôt 2x1 ?

[galaadb]

G est de dimension 2x1.

En fait, comme x est une liste, M est bien de dimension 2x2, mais il existe une version de M pour chaque valeur x(i), c'est ça le souci...

C'est justement cette dépendance en x que je cherche à obtenir;
je pense que j'arriverais à résoudre le système si M était constante en x, mais là je ne vois pas comment intégrer la dépendance dans le système pour le résoudre.

Des idées?