Why am I getting the error "Array indices must be positive integers or logical values" in MATPOWER?

6 views (last 30 days)
Lucía
Lucía on 19 Dec 2020
Edited: VBBV on 11 Aug 2023
I'm trying to run an optimal power flow in MATPOWER and I'm getting this error:
Array indices must be positive integers or logical values.
Error in ext2int (line 162)
bs(n2i(mpc.branch(:, F_BUS))) & ...
Error in opf (line 222)
mpc = ext2int(mpc, mpopt);
Error in runopf (line 75)
[r, success] = opf(casedata, mpopt);
Error in Iberia_runopf (line 9)
[results,success]=runopf('Iberia'); % llama a runopf
When I click on line 162 it takes me to this part of the code:
brs = ( mpc.branch(:, BR_STATUS) & ... %% branch status
bs(n2i(mpc.branch(:, F_BUS))) & ...
bs(n2i(mpc.branch(:, T_BUS))) );
I've checked the matrix mpc.branch and the columns BR_STATUS and F_BUS but all the data seems to be fine... I don't know what's wrong. Here is the whole matrix:
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
mpc.branch = [
108 119 0.007594 0.07594 13.16829076 0 0 0 0 0 1 -360 360;
108 104 0.004456 0.04456 22.44165171 0 0 0 0 0 1 -360 360;
109 116 0.003121 0.03121 32.0410125 0 0 0 0 0 1 -360 360;
109 142 0.012049 0.12049 8.299443937 0 0 0 0 0 1 -360 360;
109 105 0.004398 0.04398 22.737608 0 0 0 0 0 1 -360 360;
102 110 0.005339 0.05339 18.73009927 0 0 0 0 0 1 -360 360;
103 128 0.003647 0.03647 27.41979709 0 0 0 0 0 1 -360 360;
103 90 0.006485 0.06485 15.42020046 0 0 0 0 0 1 -360 360;
100 101 0.008493 0.08493 11.77440245 0 0 0 0 0 1 -360 360;
101 111 0.005131 0.05131 19.48937829 0 0 0 0 0 1 -360 360;
101 122 0.00295 0.0295 33.89830508 0 0 0 0 0 1 -360 360;
101 532 0.00424 0.0424 23.58490566 0 0 0 0 0 1 -360 360;
106 115 0.00319 0.0319 31.34796238 0 0 0 0 0 1 -360 360;
106 96 0.006157 0.06157 16.24167614 0 0 0 0 0 1 -360 360;
107 75 0.007719 0.07719 12.95504599 0 0 0 0 0 1 -360 360;
107 119 0.001344753 0.013447529 74.3631019 0 0 0 0 0 1 -360 360;
107 131 0.004107 0.04107 24.348673 0 0 0 0 0 1 -360 360;
107 97 0.003857 0.03857 25.92688618 0 0 0 0 0 1 -360 360;
107 136 0.003916 0.03916 25.53626149 0 0 0 0 0 1 -360 360;
104 120 0.004252 0.04252 23.51834431 0 0 0 0 0 1 -360 360;
104 144 0.002618 0.02618 38.19709702 0 0 0 0 0 1 -360 360;
104 97 0.003486 0.03486 28.68617326 0 0 0 0 0 1 -360 360;
104 110 0.003123 0.03123 32.02049312 0 0 0 0 0 1 -360 360;
105 146 0.008242 0.08242 12.13297743 0 0 0 0 0 1 -360 360;
105 123 0.010867 0.10867 9.202171713 0 0 0 0 0 1 -360 360;
105 140 0.004299 0.04299 23.26122354 0 0 0 0 0 1 -360 360;
105 51 0.002751 0.02751 36.35041803 0 0 0 0 0 1 -360 360;
105 99 0.003099749 0.030997486 32.26068102 0 0 0 0 0 1 -360 360;
105 114 0.00701 0.0701 14.26533524 0 0 0 0 0 1 -360 360;
105 63 0.023205 0.23205 4.309416074 0 0 0 0 0 1 -360 360;
105 116 0.001375 0.01375 72.72727273 0 0 0 0 0 1 -360 360;
105 69 0.001969491 0.019694913 50.77453199 0 0 0 0 0 1 -360 360;
105 92 0.021646 0.21646 4.619791185 0 0 0 0 0 1 -360 360;
105 95 0.001991 0.01991 50.22601708 0 0 0 0 0 1 -360 360;
99 88 0.006002151 0.060021508 16.6606944 0 0 0 0 0 1 -360 360;
98 80 0.014331 0.14331 6.97788012 0 0 0 0 0 1 -360 360;
98 116 0.016496 0.16496 6.062075655 0 0 0 0 0 1 -360 360;
91 80 0.004809 0.04809 20.79434394 0 0 0 0 0 1 -360 360;
90 86 0.00734 0.0734 13.6239782 0 0 0 0 0 1 -360 360;
93 83 0.005018 0.05018 19.92825827 0 0 0 0 0 1 -360 360;
93 122 0.008512 0.08512 11.7481203 0 0 0 0 0 1 -360 360;
93 86 0.003013 0.03013 33.18951211 0 0 0 0 0 1 -360 360;
93 118 0.013996 0.13996 7.144898542 0 0 0 0 0 1 -360 360;
93 84 0.010716 0.10716 9.331840239 0 0 0 0 0 1 -360 360;
92 87 0.004796 0.04796 20.85070892 0 0 0 0 0 1 -360 360;
95 146 0.008139 0.08139 12.28652169 0 0 0 0 0 1 -360 360;
95 59 0.009639 0.09639 10.37452018 0 0 0 0 0 1 -360 360;
95 70 0.002719 0.02719 36.77822729 0 0 0 0 0 1 -360 360;
95 64 0.007112 0.07112 14.06074241 0 0 0 0 0 1 -360 360;
95 127 0.004652657 0.046526574 21.49309356 0 0 0 0 0 1 -360 360;
94 89 0.003254 0.03254 30.7314075 0 0 0 0 0 1 -360 360;
94 97 0.002965 0.02965 33.72681282 0 0 0 0 0 1 -360 360;
97 82 0.006533 0.06533 15.30690341 0 0 0 0 0 1 -360 360;
97 64 0.007379 0.07379 13.55197181 0 0 0 0 0 1 -360 360;
96 83 0.010203 0.10203 9.80103891 0 0 0 0 0 1 -360 360;
179 182 0.009069 0.09069 11.02657404 0 0 0 0 0 1 -360 360;
179 171 0.00866 0.0866 11.54734411 0 0 0 0 0 1 -360 360;
178 166 0.013029 0.13029 7.675186123 0 0 0 0 0 1 -360 360;
178 199 0.019119 0.19119 5.230399079 0 0 0 0 0 1 -360 360;
178 190 0.01345 0.1345 7.434944238 0 0 0 0 0 1 -360 360;
177 175 0.003048 0.03048 32.80839895 0 0 0 0 0 1 -360 360;
176 183 0.002686 0.02686 37.23008191 0 0 0 0 0 1 -360 360;
176 175 0.001261 0.01261 79.30214116 0 0 0 0 0 1 -360 360;
176 172 0.001541 0.01541 64.89292667 0 0 0 0 0 1 -360 360;
175 183 0.002729 0.02729 36.64345914 0 0 0 0 0 1 -360 360;
175 168 0.006836 0.06836 14.62843768 0 0 0 0 0 1 -360 360;
175 189 0.013195 0.13195 7.578628268 0 0 0 0 0 1 -360 360;
175 172 0.000862279 0.008622792 115.9717121 0 0 0 0 0 1 -360 360;
174 199 0.004227 0.04227 23.65744027 0 0 0 0 0 1 -360 360;
174 193 0.007017 0.07017 14.25110446 0 0 0 0 0 1 -360 360;
174 157 0.002157 0.02157 46.36068614 0 0 0 0 0 1 -360 360;
173 181 0.003684 0.03684 27.14440825 0 0 0 0 0 1 -360 360;
173 167 0.003568 0.03568 28.02690583 0 0 0 0 0 1 -360 360;
172 159 0.001484134 0.014841337 67.37937458 0 0 0 0 0 1 -360 360;
172 154 0.013266 0.13266 7.53806724 0 0 0 0 0 1 -360 360;
172 170 0.004223 0.04223 23.67984845 0 0 0 0 0 1 -360 360;
171 160 0.006985 0.06985 14.31639227 0 0 0 0 0 1 -360 360;
170 159 0.007019 0.07019 14.24704374 0 0 0 0 0 1 -360 360;
168 158 0.011281 0.11281 8.86446237 0 0 0 0 0 1 -360 360;
168 144 0.002627 0.02627 38.06623525 0 0 0 0 0 1 -360 360;
168 180 0.002762 0.02762 36.20564808 0 0 0 0 0 1 -360 360;
168 154 0.006981 0.06981 14.32459533 0 0 0 0 0 1 -360 360;
169 161 0.007987 0.07987 12.52034556 0 0 0 0 0 1 -360 360;
169 163 0.00625 0.0625 16 0 0 0 0 0 1 -360 360;
164 133 0.003534 0.03534 28.29654782 0 0 0 0 0 1 -360 360;
164 188 0.007042 0.07042 14.20051122 0 0 0 0 0 1 -360 360;
164 154 0.006071 0.06071 16.47175095 0 0 0 0 0 1 -360 360;
164 193 0.010527 0.10527 9.49938254 0 0 0 0 0 1 -360 360;
164 149 0.005181 0.05181 19.30129319 0 0 0 0 0 1 -360 360;
165 18 0.000575 0.00575 173.9130435 0 0 0 0 0 1 -360 360;
165 163 0.002603 0.02603 38.41721091 0 0 0 0 0 1 -360 360;
166 151 0.032483 0.32483 3.078533387 0 0 0 0 0 1 -360 360;
166 153 0.008044 0.08044 12.43162606 0 0 0 0 0 1 -360 360;
166 155 0.013015 0.13015 7.683442182 0 0 0 0 0 1 -360 360;
166 161 0.004552 0.04552 21.96836555 0 0 0 0 0 1 -360 360;
167 160 0.002117 0.02117 47.23665564 0 0 0 0 0 1 -360 360;
160 151 0.003235 0.03235 30.91190108 0 0 0 0 0 1 -360 360;
160 156 0.004445 0.04445 22.49718785 0 0 0 0 0 1 -360 360;
162 151 0.023442 0.23442 4.265847624 0 0 0 0 0 1 -360 360;
163 157 0.002528 0.02528 39.55696203 0 0 0 0 0 1 -360 360;
59 58 0.003497 0.03497 28.59593938 0 0 0 0 0 1 -360 360;
59 48 0.003903 0.03903 25.62131694 0 0 0 0 0 1 -360 360;
59 49 0.001855 0.01855 53.9083558 0 0 0 0 0 1 -360 360;
58 49 0.004732 0.04732 21.13271344 0 0 0 0 0 1 -360 360;
55 46 0.006678 0.06678 14.97454328 0 0 0 0 0 1 -360 360;
55 74 0.004521 0.04521 22.11900022 0 0 0 0 0 1 -360 360;
54 51 0.003914 0.03914 25.54931017 0 0 0 0 0 1 -360 360;
54 42 0.003793 0.03793 26.36435539 0 0 0 0 0 1 -360 360;
54 70 0.002292 0.02292 43.63001745 0 0 0 0 0 1 -360 360;
57 56 0.009463 0.09463 10.56747332 0 0 0 0 0 1 -360 360;
57 42 0.004225 0.04225 23.66863905 0 0 0 0 0 1 -360 360;
56 88 0.007783 0.07783 12.848516 0 0 0 0 0 1 -360 360;
56 52 0.002997 0.02997 33.36670003 0 0 0 0 0 1 -360 360;
51 32 0.003249 0.03249 30.77870114 0 0 0 0 0 1 -360 360;
51 31 0.002187 0.02187 45.72473708 0 0 0 0 0 1 -360 360;
51 63 0.004515 0.04515 22.14839424 0 0 0 0 0 1 -360 360;
50 63 0.006031 0.06031 16.58099818 0 0 0 0 0 1 -360 360;
53 69 0.001599 0.01599 62.53908693 0 0 0 0 0 1 -360 360;
531 75 0.004977 0.04977 20.09242516 0 0 0 0 0 1 -360 360;
151 146 0.0017 0.017 58.82352941 0 0 0 0 0 1 -360 360;
151 154 0.000701 0.00701 142.6533524 0 0 0 0 0 1 -360 360;
151 158 0.011528 0.11528 8.674531575 0 0 0 0 0 1 -360 360;
151 149 0.001796 0.01796 55.67928731 0 0 0 0 0 1 -360 360;
150 159 0.003026 0.03026 33.04692664 0 0 0 0 0 1 -360 360;
150 141 0.003501 0.03501 28.56326764 0 0 0 0 0 1 -360 360;
152 154 0.003349 0.03349 29.8596596 0 0 0 0 0 1 -360 360;
155 149 0.017768 0.17768 5.628095453 0 0 0 0 0 1 -360 360;
154 146 0.000740174 0.007401737 135.10343 0 0 0 0 0 1 -360 360;
154 144 0.00219 0.0219 45.66210046 0 0 0 0 0 1 -360 360;
154 157 0.005558 0.05558 17.99208348 0 0 0 0 0 1 -360 360;
154 156 0.004328 0.04328 23.10536044 0 0 0 0 0 1 -360 360;
154 149 0.001204 0.01204 83.05647841 0 0 0 0 0 1 -360 360;
159 136 0.00481 0.0481 20.79002079 0 0 0 0 0 1 -360 360;
1 2 0.00117325 0.011732495 85.23336106 0 0 0 0 0 1 -360 360;
1 79 0.002073 0.02073 48.23926676 0 0 0 0 0 1 -360 360;
146 145 0.001077 0.01077 92.85051068 0 0 0 0 0 1 -360 360;
146 148 0.001932258 0.019322582 51.75291763 0 0 0 0 0 1 -360 360;
146 149 0.001961342 0.019613415 50.98551144 0 0 0 0 0 1 -360 360;
147 145 0.001845 0.01845 54.20054201 0 0 0 0 0 1 -360 360;
147 116 0.006756 0.06756 14.80165779 0 0 0 0 0 1 -360 360;
147 140 0.003363 0.03363 29.73535534 0 0 0 0 0 1 -360 360;
147 149 0.002037 0.02037 49.09180167 0 0 0 0 0 1 -360 360;
142 149 0.009905 0.09905 10.09591116 0 0 0 0 0 1 -360 360;
143 131 0.008289 0.08289 12.06418145 0 0 0 0 0 1 -360 360;
141 138 0.000582109 0.005821089 171.7891653 0 0 0 0 0 1 -360 360;
141 118 0.002375 0.02375 42.10526316 0 0 0 0 0 1 -360 360;
141 134 0.002836 0.02836 35.26093089 0 0 0 0 0 1 -360 360;
515 39 0.008438 0.08438 11.85114956 0 0 0 0 0 1 -360 360;
2 6 0.000783609 0.007836094 127.6146039 0 0 0 0 0 1 -360 360;
11 9 0.005734 0.05734 17.43983258 0 0 0 0 0 1 -360 360;
11 13 0.003911 0.03911 25.56890821 0 0 0 0 0 1 -360 360;
10 12 0.004151 0.04151 24.09058058 0 0 0 0 0 1 -360 360;
10 6 0.006738 0.06738 14.84119917 0 0 0 0 0 1 -360 360;
13 9 0.006538 0.06538 15.29519731 0 0 0 0 0 1 -360 360;
13 15 0.005141626 0.051416263 19.44909916 0 0 0 0 0 1 -360 360;
12 15 0.004595 0.04595 21.76278564 0 0 0 0 0 1 -360 360;
15 14 0.00239 0.0239 41.84100418 0 0 0 0 0 1 -360 360;
15 17 0.007751 0.07751 12.90156109 0 0 0 0 0 1 -360 360;
15 16 0.005326 0.05326 18.77581675 0 0 0 0 0 1 -360 360;
17 20 0.002653354 0.026533541 37.6881477 0 0 0 0 0 1 -360 360;
17 16 0.005277 0.05277 18.95016108 0 0 0 0 0 1 -360 360;
19 18 0.002887 0.02887 34.63803256 0 0 0 0 0 1 -360 360;
19 20 0.003769 0.03769 26.53223667 0 0 0 0 0 1 -360 360;
201 183 0.005611 0.05611 17.82213509 0 0 0 0 0 1 -360 360;
201 196 0.0068 0.068 14.70588235 0 0 0 0 0 1 -360 360;
201 198 0.003227 0.03227 30.98853424 0 0 0 0 0 1 -360 360;
201 200 0.001162 0.01162 86.05851979 0 0 0 0 0 1 -360 360;
200 207 0.00553 0.0553 18.08318264 0 0 0 0 0 1 -360 360;
203 212 0.005509 0.05509 18.15211472 0 0 0 0 0 1 -360 360;
203 214 0.010994 0.10994 9.095870475 0 0 0 0 0 1 -360 360;
203 215 0.002369 0.02369 42.21190376 0 0 0 0 0 1 -360 360;
203 199 0.000805594 0.008055939 124.1320242 0 0 0 0 0 1 -360 360;
202 199 0.00696426 0.069642597 14.35902804 0 0 0 0 0 1 -360 360;
205 208 0.002217698 0.022176975 45.09181157 0 0 0 0 0 1 -360 360;
205 204 0.003279 0.03279 30.49710278 0 0 0 0 0 1 -360 360;
205 207 0.009817 0.09817 10.18641133 0 0 0 0 0 1 -360 360;
204 194 0.012498 0.12498 8.001280205 0 0 0 0 0 1 -360 360;
204 206 0.00417 0.0417 23.98081535 0 0 0 0 0 1 -360 360;
207 192 0.007307 0.07307 13.68550705 0 0 0 0 0 1 -360 360;
206 211 0.005144 0.05144 19.44012442 0 0 0 0 0 1 -360 360;
209 211 0.010143 0.10143 9.85901607 0 0 0 0 0 1 -360 360;
208 215 0.001845786 0.018457858 54.17746636 0 0 0 0 0 1 -360 360;
208 210 0.001596 0.01596 62.6566416 0 0 0 0 0 1 -360 360;
208 193 0.005739 0.05739 17.42463844 0 0 0 0 0 1 -360 360;
3 5 0.008701 0.08701 11.49293185 0 0 0 0 0 1 -360 360;
3 4 0.005409 0.05409 18.48770568 0 0 0 0 0 1 -360 360;
60 88 0.008339 0.08339 11.99184555 0 0 0 0 0 1 -360 360;
60 45 0.007579 0.07579 13.19435282 0 0 0 0 0 1 -360 360;
61 80 0.006031 0.06031 16.58099818 0 0 0 0 0 1 -360 360;
61 63 0.012575 0.12575 7.952286282 0 0 0 0 0 1 -360 360;
61 29 0.011974 0.11974 8.351428094 0 0 0 0 0 1 -360 360;
61 79 0.006433 0.06433 15.54484688 0 0 0 0 0 1 -360 360;
62 67 0.002955 0.02955 33.84094755 0 0 0 0 0 1 -360 360;
63 76 0.002757 0.02757 36.27130939 0 0 0 0 0 1 -360 360;
63 47 0.005378 0.05378 18.59427296 0 0 0 0 0 1 -360 360;
63 73 0.00435 0.0435 22.98850575 0 0 0 0 0 1 -360 360;
63 69 0.000605 0.00605 165.2892562 0 0 0 0 0 1 -360 360;
63 80 0.002861 0.02861 34.9528137 0 0 0 0 0 1 -360 360;
64 38 0.010463 0.10463 9.557488292 0 0 0 0 0 1 -360 360;
64 49 0.004336 0.04336 23.06273063 0 0 0 0 0 1 -360 360;
64 45 0.001624 0.01624 61.57635468 0 0 0 0 0 1 -360 360;
64 70 0.00479 0.0479 20.87682672 0 0 0 0 0 1 -360 360;
64 67 0.003905 0.03905 25.60819462 0 0 0 0 0 1 -360 360;
65 77 0.005136 0.05136 19.47040498 0 0 0 0 0 1 -360 360;
65 82 0.004502 0.04502 22.21235007 0 0 0 0 0 1 -360 360;
65 49 0.005542191 0.055421906 18.04340684 0 0 0 0 0 1 -360 360;
66 86 0.001924 0.01924 51.97505198 0 0 0 0 0 1 -360 360;
67 88 0.011507 0.11507 8.690362388 0 0 0 0 0 1 -360 360;
67 85 0.007118 0.07118 14.04889014 0 0 0 0 0 1 -360 360;
68 69 0.004837 0.04837 20.67397147 0 0 0 0 0 1 -360 360;
69 27 0.00358 0.0358 27.93296089 0 0 0 0 0 1 -360 360;
69 73 0.00397 0.0397 25.18891688 0 0 0 0 0 1 -360 360;
216 215 0.001487 0.01487 67.24949563 0 0 0 0 0 1 -360 360;
216 1251 0.000001 1.00E-05 100000 0 0 0 0 0 1 -360 360;
214 215 0.010238 0.10238 9.767532721 0 0 0 0 0 1 -360 360;
215 212 0.005251 0.05251 19.04399162 0 0 0 0 0 1 -360 360;
215 213 0.011106 0.11106 9.004141905 0 0 0 0 0 1 -360 360;
213 210 0.006172 0.06172 16.2022035 0 0 0 0 0 1 -360 360;
4 8 0.005136427 0.051364267 19.4687874 0 0 0 0 0 1 -360 360;
5 7 0.003086 0.03086 32.404407 0 0 0 0 0 1 -360 360;
199 195 0.020385 0.20385 4.905567819 0 0 0 0 0 1 -360 360;
199 188 0.00398 0.0398 25.12562814 0 0 0 0 0 1 -360 360;
198 197 0.002869 0.02869 34.8553503 0 0 0 0 0 1 -360 360;
195 191 0.003323 0.03323 30.0932892 0 0 0 0 0 1 -360 360;
194 193 0.002788 0.02788 35.86800574 0 0 0 0 0 1 -360 360;
194 184 0.008127 0.08127 12.30466347 0 0 0 0 0 1 -360 360;
197 189 0.009314 0.09314 10.73652566 0 0 0 0 0 1 -360 360;
191 187 0.002009 0.02009 49.77600796 0 0 0 0 0 1 -360 360;
191 193 0.002691 0.02691 37.16090673 0 0 0 0 0 1 -360 360;
190 24 0.008605 0.08605 11.62115049 0 0 0 0 0 1 -360 360;
190 188 0.00175 0.0175 57.14285714 0 0 0 0 0 1 -360 360;
192 185 0.001462 0.01462 68.3994528 0 0 0 0 0 1 -360 360;
6 20 0.004165355 0.041653547 24.00755958 0 0 0 0 0 1 -360 360;
6 7 0.001149283 0.011492825 87.0108075 0 0 0 0 0 1 -360 360;
24 23 0.00243685 0.024368496 41.0365903 0 0 0 0 0 1 -360 360;
25 27 0.000992386 0.009923859 100.7672489 0 0 0 0 0 1 -360 360;
26 31 0.004173 0.04173 23.96357537 0 0 0 0 0 1 -360 360;
26 28 0.001544 0.01544 64.76683938 0 0 0 0 0 1 -360 360;
27 29 0.002338 0.02338 42.77159966 0 0 0 0 0 1 -360 360;
20 21 0.001142942 0.011429418 87.49352114 0 0 0 0 0 1 -360 360;
20 22 0.01024 0.1024 9.765625 0 0 0 0 0 1 -360 360;
20 23 0.00365 0.0365 27.39726027 0 0 0 0 0 1 -360 360;
21 22 0.004487 0.04487 22.28660575 0 0 0 0 0 1 -360 360;
21 23 0.002735 0.02735 36.5630713 0 0 0 0 0 1 -360 360;
28 31 0.003084 0.03084 32.42542153 0 0 0 0 0 1 -360 360;
29 46 0.005976 0.05976 16.73360107 0 0 0 0 0 1 -360 360;
29 36 0.009907 0.09907 10.09387302 0 0 0 0 0 1 -360 360;
29 40 0.005288 0.05288 18.9107413 0 0 0 0 0 1 -360 360;
29 34 0.005167 0.05167 19.35359009 0 0 0 0 0 1 -360 360;
29 79 0.012431 0.12431 8.044405116 0 0 0 0 0 1 -360 360;
7 8 0.002553622 0.025536215 39.16007121 0 0 0 0 0 1 -360 360;
182 181 0.003099 0.03099 32.2684737 0 0 0 0 0 1 -360 360;
183 189 0.00198 0.0198 50.50505051 0 0 0 0 0 1 -360 360;
183 186 0.00151 0.0151 66.22516556 0 0 0 0 0 1 -360 360;
180 185 0.001563 0.01563 63.97952655 0 0 0 0 0 1 -360 360;
181 184 0.002944 0.02944 33.9673913 0 0 0 0 0 1 -360 360;
187 184 0.002453 0.02453 40.76640848 0 0 0 0 0 1 -360 360;
77 75 0.012732 0.12732 7.854225573 0 0 0 0 0 1 -360 360;
76 80 0.001579 0.01579 63.33122229 0 0 0 0 0 1 -360 360;
75 84 0.007673 0.07673 13.03271211 0 0 0 0 0 1 -360 360;
74 79 0.002004206 0.020042061 49.8950679 0 0 0 0 0 1 -360 360;
73 70 0.002963388 0.029633881 33.74515805 0 0 0 0 0 1 -360 360;
72 83 0.002887 0.02887 34.63803256 0 0 0 0 0 1 -360 360;
71 80 0.002548494 0.025484937 39.23886457 0 0 0 0 0 1 -360 360;
71 79 0.00549 0.0549 18.21493625 0 0 0 0 0 1 -360 360;
70 45 0.00555 0.0555 18.01801802 0 0 0 0 0 1 -360 360;
79 80 0.001766 0.01766 56.62514156 0 0 0 0 0 1 -360 360;
78 115 0.004348 0.04348 22.99908004 0 0 0 0 0 1 -360 360;
9 8 0.003506463 0.035064625 28.5187704 0 0 0 0 0 1 -360 360;
39 43 0.002113 0.02113 47.32607667 0 0 0 0 0 1 -360 360;
39 49 0.001273075 0.012730746 78.54999097 0 0 0 0 0 1 -360 360;
38 42 0.002651 0.02651 37.72161449 0 0 0 0 0 1 -360 360;
38 35 0.003555 0.03555 28.12939522 0 0 0 0 0 1 -360 360;
33 31 0.001922 0.01922 52.02913632 0 0 0 0 0 1 -360 360;
32 31 0.001114 0.01114 89.76660682 0 0 0 0 0 1 -360 360;
31 30 0.002725 0.02725 36.69724771 0 0 0 0 0 1 -360 360;
31 47 0.013426 0.13426 7.448234768 0 0 0 0 0 1 -360 360;
30 35 0.009233954 0.092339545 10.82959637 0 0 0 0 0 1 -360 360;
37 45 0.001664895 0.016648953 60.06383779 0 0 0 0 0 1 -360 360;
37 43 0.001839 0.01839 54.37737901 0 0 0 0 0 1 -360 360;
37 49 0.002353 0.02353 42.49893753 0 0 0 0 0 1 -360 360;
35 42 0.002622 0.02622 38.13882532 0 0 0 0 0 1 -360 360;
8 133 0.004468 0.04468 22.38137869 0 0 0 0 0 1 -360 360;
8 127 0.00596 0.0596 16.77852349 0 0 0 0 0 1 -360 360;
133 127 0.000488527 0.004885275 204.6967768 0 0 0 0 0 1 -360 360;
132 115 0.00537 0.0537 18.62197393 0 0 0 0 0 1 -360 360;
132 122 0.006811 0.06811 14.68213185 0 0 0 0 0 1 -360 360;
132 139 0.004847 0.04847 20.63131834 0 0 0 0 0 1 -360 360;
131 121 0.005178 0.05178 19.31247586 0 0 0 0 0 1 -360 360;
131 123 0.009032 0.09032 11.07174491 0 0 0 0 0 1 -360 360;
131 128 0.002364459 0.023644594 42.29296579 0 0 0 0 0 1 -360 360;
131 139 0.021881 0.21881 4.570175038 0 0 0 0 0 1 -360 360;
130 122 0.001936 0.01936 51.65289256 0 0 0 0 0 1 -360 360;
130 137 0.010249 0.10249 9.757049468 0 0 0 0 0 1 -360 360;
137 124 0.004968 0.04968 20.12882448 0 0 0 0 0 1 -360 360;
136 138 0.001808 0.01808 55.30973451 0 0 0 0 0 1 -360 360;
135 114 0.007034 0.07034 14.21666193 0 0 0 0 0 1 -360 360;
135 127 0.011191 0.11191 8.935751944 0 0 0 0 0 1 -360 360;
134 129 0.0021 0.021 47.61904762 0 0 0 0 0 1 -360 360;
138 115 0.003535 0.03535 28.28854314 0 0 0 0 0 1 -360 360;
138 128 0.002568 0.02568 38.94080997 0 0 0 0 0 1 -360 360;
138 129 0.002527 0.02527 39.57261575 0 0 0 0 0 1 -360 360;
88 45 0.01284 0.1284 7.788161994 0 0 0 0 0 1 -360 360;
89 85 0.004601 0.04601 21.73440556 0 0 0 0 0 1 -360 360;
80 87 0.005234 0.05234 19.10584639 0 0 0 0 0 1 -360 360;
81 84 0.004475 0.04475 22.34636872 0 0 0 0 0 1 -360 360;
121 119 0.002835 0.02835 35.27336861 0 0 0 0 0 1 -360 360;
122 115 0.005224 0.05224 19.1424196 0 0 0 0 0 1 -360 360;
122 126 0.002666 0.02666 37.50937734 0 0 0 0 0 1 -360 360;
122 128 0.006966 0.06966 14.35544071 0 0 0 0 0 1 -360 360;
122 129 0.001128 0.01128 88.65248227 0 0 0 0 0 1 -360 360;
122 118 0.007012 0.07012 14.2612664 0 0 0 0 0 1 -360 360;
124 117 0.008545 0.08545 11.70275015 0 0 0 0 0 1 -360 360;
125 116 0.00145 0.0145 68.96551724 0 0 0 0 0 1 -360 360;
127 114 0.002525274 0.025252736 39.59966885 0 0 0 0 0 1 -360 360;
127 116 0.006064 0.06064 16.49076517 0 0 0 0 0 1 -360 360;
128 117 0.006601 0.06601 15.14921982 0 0 0 0 0 1 -360 360;
129 118 0.000514 0.00514 194.5525292 0 0 0 0 0 1 -360 360;
115 111 0.005426 0.05426 18.42978253 0 0 0 0 0 1 -360 360;
115 118 0.001711 0.01711 58.44535359 0 0 0 0 0 1 -360 360;
114 113 0.002787 0.02787 35.88087549 0 0 0 0 0 1 -360 360;
113 112 0.004113555 0.041135546 24.3098757 0 0 0 0 0 1 -360 360;
48 49 0.003258 0.03258 30.6936771 0 0 0 0 0 1 -360 360;
48 41 0.004416 0.04416 22.64492754 0 0 0 0 0 1 -360 360;
44 40 0.004046 0.04046 24.71576866 0 0 0 0 0 1 -360 360;
45 41 0.002078 0.02078 48.12319538 0 0 0 0 0 1 -360 360;
];

Sign in to answer this question.

Answers (3)

Image Analyst
Image Analyst on 19 Dec 2020
  1 Comment
Lucía
Lucía on 19 Dec 2020
Thank you for your answer. I have checked the affected columns in my matrix (F_BUS & BR_STATUS) and they're all positive integers.

Sign in to comment.

Lechuan Piao
Lechuan Piao on 23 Mar 2023
I have the same issues here with MATPOWER. Wonder if you fixed that? @Lucía
VBBV
VBBV on 11 Aug 2023
Edited: VBBV on 11 Aug 2023
Ran in:
The problem is with nested matrices bs or n2i matrices.
% brs = ( mpc.branch(:, BR_STATUS) & ... %% branch status
% bs(n2i(mpc.branch(:, F_BUS))) & ...
% bs(n2i(mpc.branch(:, T_BUS))) );
You seem to access elements of both matrices using the logical index method, You can do it as follows
A = rand(100,45); % assume this as mpc.branch matrix
F_BUS = 4; % column index
T_BUS = 10;
BS = rand(100,20); % the problem is here ... with this nested matrix
idx = A(:,F_BUS) > 0.3; % find the logical indices
BS(idx) % extract from the matrix
ans = 78×1
0.3528 0.6659 0.2760 0.7496 0.5803 0.7379 0.9666 0.6620 0.3568 0.2353
Instead you were trying to plug in the elements (some +ve or -ve decima or some otherl values) of bs and n2i matrices as indices of the matrix BS, which will throw same error. Even though F_BUS & T_BUS are positive integers, yet Matlab throws error, because the elements of bs and n2i are not positive integers which will now act as array indices. Since you have enclosed them in parenthesis all together.
A = rand(100,45); % assume this as mpc.branch matrix
F_BUS = 4; % column index
T_BUS = 10;
BS = rand(100,20); % the problem is here ... with this nested matrix
A(BS(:,F_BUS))
Array indices must be positive integers or logical values.

Categories

Find more on MATLAB in Help Center and File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!