Discussion :

[redui]

Bonjour,

Dans le cadre dans mon stage sur les écoulements en milieu poreux, nous utilisons MATLAB comme outil. Je dois déterminer un champ de vitesse sous la forme :



K et p sont une matrice. Nous avons des opérations de vecteurs à chaque boucle, mon code est le suivant :

Code :

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

1
2
3
4
5
6
for i=2:nx
           A=((xs(i+1,1)-xs(i-1,1))\(((xs(i+1)-xs(i))\KXX(i,:))+((xs(i)-xs(i-1))\KXX(i-1,:))));
           B=(2*(p(i,:)-p(i-1,:)))/((xs(i+1)+xs(i))-(xs(i)+xs(i-1)));
 
           u(i,:)=A*B;
end

J'ai une erreur de type :

Error using  * Inner matrix dimensions must agree

J'ai bien regardé et mes matrices sont bien à la même taille. Je doute aussi que mes expressions sur MATLAB soient correctes mathématiquement

[Jerome Briot]

Que renvoi ceci ?

Code :

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

whos A B

L'erreur dit que les deux matrices doivent avoir la même dimension "interne"

Tu peux multiplier une matrice 4x1 par une matrice 1x4 mais pas par une autre matrice 4x1 (règle du calcul matriciel)

Ne souhaites-tu pas plutôt faire une multiplication terme à terme ?

Code :

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

u(i,:)= A.*B;

Voir la FAQ : Inner matrix dimensions must agree.

[redui]

Bonjour,

cela me renvoie :

Code :

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

1
2
3
4
5
 
Name      Size             Bytes  Class     Attributes
 
  A         1x100              800  double              
  B         1x100              800  double

[Jerome Briot]

Il te suffit donc de bien relire mon précédent message