Discussion :

[lematlabeur]

salut,
j'ai un problème.
J'ai créé des sous fonctions avec des boucles for. j'ai fait des tests pour voir si elles fonctionnent.
Mais quand je fait appel à l'une d'elles dans le programme principal,elle ne fonctionne plus.
j'ai le message suivant:

Subscripted assignment dimension mismatch

j'ai fait des vérifications et j'ai trouvé que ça bloque sur cette sous-fonction.
qq peut m'aider svp?

[Jerome Briot]

Tu peux déjà jeter un oeil à la : Subscripted assignment dimension mismatch.

Ensuite, montre nous le bout de code où le problème se produit.
Voir encore la (Comment exploiter les messages d'erreurs ?) pour déterminer la zone du code qui pose problème.

[lematlabeur]

voila la sous fonction:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
function T=R1(Te,C,Vmax,Ko,Km,Ro)
    [nbca] = size(Te);
    [nblb,nbcb] = size(C); 
      B= [0 0 0 0 0 ]; 
 
for i=1:nblb-3
    for j=1:nbcb-1
  B(i,j)=-Vmax*((C(i,j)/((C(i,j)+Km))*((C(i+2,j)/((C(i+2,j)+Ko))))))*2.78*10^(-10)*(Ro/(1+Te(j)));
    end
end
 
   T=[1;-0.8;1.8;0]*B;
end

voila le bout de programme où je fait appel à la sous-fonction

Code :

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

1
2
3
4
5
6
7
 while t{k}<30
       k=k+1
            t{k}=t{k-1}+h;
    for i=1:4
        for j=1:5
 
               C{k}(i,j)= C{k-1}(i,j)+h*(((f1(C{k-1}(i,j),ka,u,deltaX,D,Tp)-f2(C{k-1}(i,j),ka,u,deltaX,D,Tp))/(deltaX*a)))+R1(Te,C{k-1}(i,j),Vmax,Ko,Km,Ro)/a;

[Jerome Briot]

Quel est le message d'erreur complet ?
(voir le lien vers la dans mon précédent message)

[lematlabeur]

voila le message d'erreur complet:

k =

2

??? Subscripted assignment dimension mismatch.

Error in ==> fihamdoullilah22 at 32
C{k}(i,j)= C{k-1}(i,j)+h*(((flux(C{k-1}(i,j),ka,u,deltaX,D,Tp)-flux2(C{k-1}(i,j),ka,u,deltaX,D,Tp))/(deltaX*a)))+R1(Teau,C{k-1}(i,j),Vmax,Koxygene,Kmethane,Ro)/a;

[Jerome Briot]

A première vue, C{k}(i,j) contient un scalaire (un nombre).

Donc l'autre membre de l'affectation, C{k-1}(i,j)+h*(((flux(C{k-1}(i,j),ka,u,deltaX,D,Tp)-flux2(C{k-1}(i,j),ka,u,deltaX,D,Tp))/(deltaX*a)))+R1(Teau,C{k-1}(i,j),Vmax,Koxygene,Kmethane,Ro)/a;, devrait aussi être un nombre.

Il faudrait vérifier ses dimensions en stockant sa valeur dans une variable temporaire.

Essaie ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
while t{k}<30
       k=k+1
            t{k}=t{k-1}+h;
    for i=1:4
        for j=1:5
 
            temp=h*(((f1(C{k-1}(i,j),ka,u,deltaX,D,Tp)-f2(C{k-1}(i,j),ka,u,deltaX,D,Tp))/(deltaX*a)))+R1(Te,C{k-1}(i,j),Vmax,Ko,Km,Ro)/a;
 
            try
               C{k}(i,j)= C{k-1}(i,j)+temp;
            catch
               size(temp)
               size(C{k}(i,j))
               size(C{k-1}(i,j))
               error('Probleme')
            end

Et montre nous ce que MATLAB renvoie

Sinon pourquoi utilises-tu des tableaux de cellules pour C ou t ?

[lematlabeur]

voila la réponse obtenue

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
k =
 
     2
 
 
ans =
 
     4     4
 
??? Index exceeds matrix dimensions.
 
Error in ==> fihamdoullilah22 at 45
                size(C{k}(i,j))

[Jerome Briot]

Mouais... mais non..; j'ai fait une petite boulette

On reprend avec ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
while t{k}<30
       k=k+1
            t{k}=t{k-1}+h;
    for i=1:4
        for j=1:5
 
            temp=h*(((f1(C{k-1}(i,j),ka,u,deltaX,D,Tp)-f2(C{k-1}(i,j),ka,u,deltaX,D,Tp))/(deltaX*a)))+R1(Te,C{k-1}(i,j),Vmax,Ko,Km,Ro)/a;
 
            try
               C{k}(i,j)= C{k-1}(i,j)+temp;
            catch
               size(temp)
               size(C{k-1}(i,j))
               error('Probleme')
            end

Parce que C{k}(i,j) n'est pas affecté dans le TRY-CATCH

[lematlabeur]

salut,
quand j'ai maintenant le message d'erreur suivant:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
k =
 
     2
 
 
ans =
 
     4     4
 
 
ans =
 
     1     1
 
??? Error using ==> filhamdoullilah22 at 46
Probleme

[Jerome Briot]

Donc, en simplifiant, tu essaies de stocker temp de dimension 4x4 dans C{k}(i,j) de dimension 1x1... d'où le message d'erreur Subscripted assignment dimension mismatch.

[lematlabeur]

C'est vrai. merci la fonction ne bloque plus.
mais il ya encore un problème.
la fonction "temp" calcule toujours une matrice nulle alors que si je fais le calcul de la première matrice de temp avec (C{k-1}=C{1})manuellement, elle différente de 0.

[lematlabeur]

non je me suis rendu compte que tout va bien.
merci pour votre aide.
bonne soirée.