Discussion :

[saad.pizzicato]

Bonjour,
En saisissant les valeurs à assigner à une matrice carrée K d'ordre n*m selon la boucle explicitée ci-dessous, je suis confronté au message d'erreur "subscript indices must be either positive integers less than 2^31 or logicals" qui résulterait du fait que pour certaine valeur de l'indice i, l'indice de la matrice K se trouve négatif.

Comment je peux, dans mon cas surmonter ce problème en assignant à tout élément de K ayant un indice négatif la valeur 0 ou en le sautant tout simplement??

Merci à l'avance pour votre aide.

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
for i=1:n*m
K(i,i)= (2*CBX(i)/(DELTAX^2))+(2*CBY(i)/(DELTAY^2))+(4*POISSON*(CBX(i)+CBY(i))/(DELTAX*DELTAY))+(CBX(i+1)/(DELTAX^2))+(CBX(i-1)/(DELTAX^2))+(CBY(i+n)/(DELTAY^2))
+(CBY(i-n)/(DELTAY^2))+((CTXY(i)+CTYX(i))/(DELTAX*DELTAY))+((CTXY(i-1)+CTYX(i-1))/(DELTAX*DELTAY))+((CTXY(i-n)+CTYX(i-n))/(DELTAX*DELTAY))+((CTXY(i-n-1)+CTYX(i-n-1))/(DELTAX*DELTAY));
K(i,i-1)= (CBX(i)/(DELTAX^2))+(POISSON*(CBX(i)+CBY(i))/(DELTAX*DELTAY))-(2*CBX(i-1)/(DELTAX^2))-(POISSON*(CBX(i-1)+CBY(i-1))/(DELTAX*DELTAY))+((CTXY(i-1)+CTYX(i-1))/(DELTAX*DELTAY))-((CTXY(i-n-1)+CTYX(i-n-1))/(DELTAX*DELTAY));
K(i,i+1)=-(CBX(i)/(DELTAX^2))-(POISSON*(CBX(i)+CBY(i))/(DELTAX*DELTAY))+(2*CBX(i+1)/(DELTAX^2))+(POISSON*(CBX(i+1)+CBY(i+1))/(DELTAX*DELTAY))+((CTXY(i)+CTYX(i))/(DELTAX*DELTAY))-((CTXY(i-n)+CTYX(i-n))/(DELTAX*DELTAY));
K(i,i+n)=-(CBY(i)/(DELTAY^2))-(POISSON*(CBX(i)+CBY(i))/(DELTAX*DELTAY))+(2*CBY(i+n)/(DELTAY^2))+(POISSON*(CBX(i+n)+CBY(i+n))/(DELTAX*DELTAY))+((CTXY(i)+CTYX(i))/(DELTAX*DELTAY))-((CTXY(i-1)+CTYX(i-1))/(DELTAX*DELTAY));
K(i,i-n)=(CBY(i)/(DELTAY^2))+(POISSON*(CBX(i)+CBY(i))/(DELTAX*DELTAY))-(CBY(i-n)/(DELTAY^2))-(POISSON*(CBX(i-n)+CBY(i-n))/(DELTAX*DELTAY))+((CTXY(i-n)+CTYX(i-n))/(DELTAX*DELTAY))-((CTXY(i-n-1)+CTYX(i-n-1))/(DELTAX*DELTAY));
K(i,i+2)=-(CBX(i+1)/(DELTAX^2));
K(i,i-2)=-(CBX(i-1)/(DELTAX^2));
K(i,i+n+1)=-(POISSON*(CBX(i+1)+CBY(i+1))/(2*DELTAX*DELTAY))-(POISSON*(CBX(i+n)+CBY(i+n))/(2*DELTAX*DELTAY))-((CTXY(i)+CTYX(i))/(DELTAX*DELTAY));
K(i,i+n-1)=(POISSON*(CBX(i-1)+CBY(i-1))/(2*DELTAX*DELTAY))+(POISSON*(CBX(i+n)+CBY(i+n))/(2*DELTAX*DELTAY))+((CTXY(i-1)+CTYX(i-1))/(DELTAX*DELTAY));
K(i,i-n+1)=(POISSON*(CBX(i+1)+CBY(i+1))/(2*DELTAX*DELTAY))+(POISSON*(CBX(i-n)+CBY(i-n))/(2*DELTAX*DELTAY))+((CTXY(i-n)+CTYX(i-n))/(DELTAX*DELTAY));
K(i,i-n-1)=-(POISSON*(CBX(i-1)+CBY(i-1))/(2*DELTAX*DELTAY))-(POISSON*(CBX(i-n)+CBY(i-n))/(2*DELTAX*DELTAY))+((CTXY(i-n-1)+CTYX(i-n-1))/(DELTAX*DELTAY));
K(i,i+(2*n))=-(CBY(i+n)/(DELTAY^2));
K(i,i-(2*n))=-(CBY(i-n)/(DELTAY^2));                 
end

[le fab]

salut

pour "sauter" le problème, rien de tel que les instructions try et catch
mais il vaut mieux bien maîtriser les indices avant de les utiliser !

Fabien

[Jerome Briot]

Tu peux fixer la taille définitive de K avant la boucle. Ensuite, il faudra appliquer les conditions correspondantes à chaque ligne :

Code :

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

1
2
3
4
5
6
7
8
9
10
K = zeros(...,...);
 
for i=1:n*m
    if (i-1)>0 && (i-n)>0 && (i-n-1)>0 && (i+1)<?? && (i+n)<??
        K(i,i)= (2*CBX(i)/(DELTAX^2))+(2*CBY(i)/(DELTAY^2))+(4*POISSON*(CBX(i)+CBY(i))/(DELTAX*DELTAY))+(CBX(i+1)/(DELTAX^2))+(CBX(i-1)/(DELTAX^2))+(CBY(i+n)/(DELTAY^2))+(CBY(i-n)/(DELTAY^2))+((CTXY(i)+CTYX(i))/(DELTAX*DELTAY))+((CTXY(i-1)+CTYX(i-1))/(DELTAX*DELTAY))+((CTXY(i-n)+CTYX(i-n))/(DELTAX*DELTAY))+((CTXY(i-n-1)+CTYX(i-n-1))/(DELTAX*DELTAY));
    end
 
    ...
 
end

De manière plus générale, tu pourrais améliorer la lisibilité de ton code en supprimant les parenthèses inutiles et en raccourcissant le nom des variables.

Par exemple :

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
17
18
19
20
DXY = DELTAX*DELTAY;
DX2 = DELTAX^2;
DY2 = DELTAY^2;
P = POISSON;
 
for i = 1:n*m
    K(i,i)     = 2*CBX(i)/DX2 + 2*CBY(i)/DY2 + 4*P*(CBX(i)+CBY(i))/DXY + CBX(i+1)/DX2 + CBX(i-1)/DX2 + CBY(i+n)/DY2 + CBY(i-n)/DY2 + (CTXY(i)+CTYX(i))/DXY + (CTXY(i-1)+CTYX(i-1))/DXY + (CTXY(i-n)+CTYX(i-n))/DXY + (CTXY(i-n-1)+CTYX(i-n-1))/DXY;
    K(i,i-1)   =  CBX(i)/DX2 + P*(CBX(i)+CBY(i))/DXY - 2*CBX(i-1)/DX2 - P*(CBX(i-1)+CBY(i-1))/DXY + (CTXY(i-1)+CTYX(i-1))/DXY - (CTXY(i-n-1)+CTYX(i-n-1))/DXY;
    K(i,i+1)   = -CBX(i)/DX2 - P*(CBX(i)+CBY(i))/DXY + 2*CBX(i+1)/DX2 + P*(CBX(i+1)+CBY(i+1))/DXY + (CTXY(i)+CTYX(i))/DXY     - (CTXY(i-n)+CTYX(i-n))/DXY;
    K(i,i+n)   = -CBY(i)/DY2 - P*(CBX(i)+CBY(i))/DXY + 2*CBY(i+n)/DY2 + P*(CBX(i+n)+CBY(i+n))/DXY + (CTXY(i)+CTYX(i))/DXY     - (CTXY(i-1)+CTYX(i-1))/DXY;
    K(i,i-n)   =  CBY(i)/DY2 + P*(CBX(i)+CBY(i))/DXY - CBY(i-n)/DY2   - P*(CBX(i-n)+CBY(i-n))/DXY + (CTXY(i-n)+CTYX(i-n))/DXY - (CTXY(i-n-1)+CTYX(i-n-1))/DXY;
    K(i,i+2)   = -CBX(i+1)/DX2;
    K(i,i-2)   = -CBX(i-1)/DX2;
    K(i,i+n+1) = -P*(CBX(i+1)+CBY(i+1))/(2*DXY) - P*(CBX(i+n)+CBY(i+n))/(2*DXY) - (CTXY(i)+CTYX(i))/DXY;
    K(i,i+n-1) =  P*(CBX(i-1)+CBY(i-1))/(2*DXY) + P*(CBX(i+n)+CBY(i+n))/(2*DXY) + (CTXY(i-1)+CTYX(i-1))/DXY;
    K(i,i-n+1) =  P*(CBX(i+1)+CBY(i+1))/(2*DXY) + P*(CBX(i-n)+CBY(i-n))/(2*DXY) + (CTXY(i-n)+CTYX(i-n))/DXY;
    K(i,i-n-1) = -P*(CBX(i-1)+CBY(i-1))/(2*DXY) - P*(CBX(i-n)+CBY(i-n))/(2*DXY) + (CTXY(i-n-1)+CTYX(i-n-1))/DXY;
    K(i,i+2*n) = -CBY(i+n)/DY2;
    K(i,i-2*n) = -CBY(i-n)/DY2;
end

Il est encore possible de factoriser certains termes mais cela peut devenir nuisible si on veut clairement séparer leurs contributions propres dans l'équation.

[saad.pizzicato]

Bonsoir M. Fab
Vous pouvez me donner, s'il vous plait, un exemple pour son utilisation quant à mon problème?

Bonsoir M.DUT
J'ai essayé cette commande de "if, else if, else" mais dans mon cas cela va générer un très grand nombre de conditions (environ 16 pour chaque expression).

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
%K(i,i)%
if (i+1)<=(n*m) && (i-1)>=1 && (i+n)<=(n*m) && (i-n)>=1 && (i-n-1)>=1
    K(i,i)     = 2*CBX(i)/DX2 + 2*CBY(i)/DY2 + 4*P*(CBX(i)+CBY(i))/DXY + CBX(i+1)/DX2 + CBX(i-1)/DX2 + CBY(i+n)/DY2 + CBY(i-n)/DY2 + (CTXY(i)+CTYX(i))/DXY + (CTXY(i-1)+CTYX(i-1))/DXY + (CTXY(i-n)+CTYX(i-n))/DXY + (CTXY(i-n-1)+CTYX(i-n-1))/DXY;
else if (i+1)<=(n*m) && (i-1)>=1 && (i+n)<=(n*m) && (i-n)>=1 && (i-n-1)<1
    K(i,i)     = 2*CBX(i)/DX2 + 2*CBY(i)/DY2 + 4*P*(CBX(i)+CBY(i))/DXY + CBX(i+1)/DX2 + CBX(i-1)/DX2 + CBY(i+n)/DY2 + CBY(i-n)/DY2 + (CTXY(i)+CTYX(i))/DXY + (CTXY(i-1)+CTYX(i-1))/DXY + (CTXY(i-n)+CTYX(i-n))/DXY;
else if (i+1)<=(n*m) && (i-1)>=1 && (i+n)<=(n*m) && (i-n)<1 && (i-n-1)<1 
    K(i,i)     = 2*CBX(i)/DX2 + 2*CBY(i)/DY2 + 4*P*(CBX(i)+CBY(i))/DXY + CBX(i+1)/DX2 + CBX(i-1)/DX2 + CBY(i+n)/DY2 + (CTXY(i)+CTYX(i))/DXY + (CTXY(i-1)+CTYX(i-1))/DXY;
else if (i+1)<=(n*m) && (i-1)>=1 && (i+n)>(n*m) && (i-n)<1 && (i-n-1)<1
 
....etc

Il n'y a pas une façon plus simple de faire de sorte que les expressions avec indice négatif se trouvent annulée??

Merci beaucoup pour votre aide.

[thecrazydonut]

Bonsoir M.DUT
J'ai essayé cette commande de "if, else if, else" mais dans mon cas cela va générer un très grand nombre de conditions (environ 16 pour chaque expression).

Les if else sont facilement lisible et extrêmement rapides à exécuter. Je pense pas que cela nuise à ton code. Ensuite, tu n'es pas obligé de vérifier la condition avant chaque ligne de code, ça fonctionne par bloc. Tout ce qui se trouve entre le if et le end correspondant sera exécute si jamais la condition est remplie. Donc si jamais tu n'as qu'une condition à vérifier, vérifie là une seule fois.

[le fab]

Bonsoir M. Fab
Vous pouvez me donner, s'il vous plait, un exemple pour son utilisation quant à mon problème?

un truc du style :

Code :

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

1
2
3
4
5
6
7
8
for ii = 1: ....
  try 
    % mon code qui risque de foirer
  catch
    % si ca foire j'atteri ici, et je continue ma boucle for à l'iteration suivante
    continue
  end
end

mais ca, c'est seulement si on est pas en mesure de garantir le bon déroulement de son code
si c'est juste un question d'indice pouvant éventuellement être négatif, un test via if avant utilisation est bien plus propre

Fab

[saad.pizzicato]

Merci pour vos aides précieuses,
Finalement j'ai pensé à cette formulation,
J'aimerai bien avoir vos avis:

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
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
 
 
K = zeros(n*m,n*m);
 
DXY = DELTAX*DELTAY;
DX2 = DELTAX^2;
DY2 = DELTAY^2;
P = POISSON;
 
for i = 1:n*m
 
%Calcul de K(i,i)%
K1(i,i)=2*CBX(i)/DX2 + 2*CBY(i)/DY2 + 4*P*(CBX(i)+CBY(i))/DXY + (CTXY(i)+CTYX(i))/DXY;
if i+1<=n*m
K2(i,i)= K1(i,i)+ CBX(i+1)/DX2;
else 
K2(i,i)= K1(i,i);
if i-1>=1
K3(i,i)=K2(i,i)+ CBX(i-1)/DX2 + (CTXY(i-1)+CTYX(i-1))/DXY;
else
K3(i,i)= K2(i,i);
if i+n<=n*m
K4(i,i)=K3(i,i)+ CBY(i+n)/DY2;
else
K4(i,i)=K3(i,i);
if i-n>=1
K5(i,i)= K4(i,i) + CBY(i-n)/DY2   + (CTXY(i-n)+CTYX(i-n))/DXY;
else
K5(i,i)=K4(i,i);
if i-n-1>=1
K6(i,i)=K5(i,i) + (CTXY(i-n-1)+CTYX(i-n-1))/DXY;
else
K6(i,i)=K5(i,i);
 
K(i,i) = K6(i,i);

Idem pour les autres termes.