construction de la matrice RR
Bonsoir,
je veut construire une matrice qui contient les éléments des résistance suivante:
Code:
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
| SE RE Resistance
[0 1 0.117
1 2 0.10725
2 3 0.16445
3 4 0.1495
4 5 0.1495
5 6 0.3144
6 7 0.2096
7 8 0.3144
8 9 0.2096
9 10 0.131
10 11 0.1048
2 12 0.1572
12 13 0.2096
13 14 0.1048
14 15 0.0524
5 16 0.1794
16 17 0.16445
17 18 0.2079
18 19 0.189
19 20 0.189
20 21 0.262
21 22 0.262
22 23 0.3144
23 24 0.2096
24 25 0.131
25 26 0.1048
6 27 0.1572
27 28 0.1572
28 29 0.1572
9 30 0.1572
30 31 0.2096
31 32 0.1572
32 33 0.1048] |
je veut faire une matrice RR(32,33) telle que si il y a une relation entre SE et RE on met la valeur de la résistance, si il n'y a pas on met 0, par exemple
dans le cas de SE=0 et RE=1 il y a une résistance R=0.117 , mais dans le cas de SE=0 et RE=0 sa n'existe pas alors on met 0, un autre exemple dans le cas ou SE=1 et RE=1 n'existe pas alors on met 0.
Code:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
| RR= zeros(nnb,nnb+1)
%des éléments non-diagonaux Yij
for k=1:nbranch
RR(SE(k),RE(k)) = R(SE(k),RE(k));
RR(SE(k),RE(k)) = RR(SE(k),RE(k));
end
%des éléments diagonaux Yii
for m=1:nnb
for n=1:nbranch
if SE(n) == m & RE(n) == m
RR(m,m) = 0;
end
end
end
RR; |