Hi,
I'm actually stuck with trying to solve the load flow for a radial
distribution, and i need your help because i really don't know what to
do.
The system is composed of 91 nodes and 1 generator, i attached to this
message the test case with null and non null values.
This case don't converge to a result, and matlab tells me:
"Warning: Matrix is singular to working precision.
In newtonpf at 186
In runpf at 226"
So i fisrtly thought that null values in bus Pd, Qd was the problem but
even when i put small value to it, it don't solve the problem.
I tried to see where the values become null or infinite and it comes
from the computation of "dx=-(J/F)"
J and F have non null value but are big (~1e+6)
I also think that BR_R and BR_X are to small, so to be sure here is my
per unit calculation:
S_ref=10Mva, V_ref=20Kv
Z_ref=(20000^2)/(10000000/3)
And for P an Q values are in MW and MVar
And i also tried to run a continuation PF but it can't handle "nan"
values.
So what can you advice to me to try to find a solution ?
Thank you very much
Best regards
Abdelkrim Ali zazou
function mpc = case_01_1_dep
%CASE30 Power flow data for 30 bus, 6 generator case.
% Please see CASEFORMAT for details on the case file format.
%
% Based on data from ...
% Alsac, O. & Stott, B., "Optimal Load Flow with Steady State Security",
% IEEE Transactions on Power Apparatus and Systems, Vol. PAS 93, No. 3,
% 1974, pp. 745-751.
% ... with branch parameters rounded to nearest 0.01, shunt values divided
% by 100 and shunt on bus 10 moved to bus 5, load at bus 5 zeroed out.
% Generator locations, costs and limits and bus areas were taken from ...
% Ferrero, R.W., Shahidehpour, S.M., Ramesh, V.C., "Transaction analysis
% in deregulated power systems using game theory", IEEE Transactions on
% Power Systems, Vol. 12, No. 3, Aug 1997, pp. 1340-1347.
% Generator Q limits were derived from Alsac & Stott, using their Pmax
% capacities. V limits and line |S| limits taken from Alsac & Stott.
% MATPOWER
% $Id: case30.m,v 1.12 2010/03/10 18:08:13 ray Exp $
%% MATPOWER Case Format : Version 2
mpc.version = '2';
mpc.baseMVA = 10;
%% Bus Data
mpc.bus= [
1 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
2 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0 ;
3 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0 ;
4 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
5 1 0.03375 0.01634587 0 0 1 1
25.84193 20 1 1.05 0 0;
6 1 0.0918 0.04446078 0 0 1 1
25.84193 20 1 1.05 0 0;
7 1 0.0594 0.02876874 0 0 1 1
25.84193 20 1 1.05 0 0;
8 1 0.0828 0.04010188 0 0 1 1
25.84193 20 1 1.05 0 0;
9 1 0.1044 0.05056324 0 0 1 1
25.84193 20 1 1.05 0 0;
10 1 0.0351 0.01699971 0 0 1 1
25.84193 20 1 1.05 0 0;
11 1 0.0189 0.00915369 0 0 1 1
25.84193 20 1 1.05 0 0;
12 1 0.0567 0.02746107 0 0 1 1
25.84193 20 1 1.05 0 0;
13 1 0.0981 0.04751201 0 0 1 1
25.84193 20 1 1.05 0 0;
14 1 0.0243 0.01176903 0 0 1 1
25.84193 20 1 1.05 0 0;
15 1 0.1602 0.07758842 0 0 1 1
25.84193 20 1 1.05 0 0;
16 1 0.0162 0.00784602 0 0 1 1
25.84193 20 1 1.05 0 0;
17 1 0.0666 0.03225586 0 0 1 1
25.84193 20 1 1.05 0 0;
18 1 0.027 0.0130767 0 0 1 1
25.84193 20 1 1.05 0 0;
19 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
20 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
21 1 0.0828 0.04010188 0 0 1 1
25.84193 20 1 1.05 0 0;
22 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
23 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
24 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0 ;
25 1 0.0108 0.00523068 0 0 1 1
25.84193 20 1 1.05 0 0;
26 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
27 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
28 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
29 1 0.0513 0.02484573 0 0 1 1
25.84193 20 1 1.05 0 0;
30 1 0.0216 0.01046136 0 0 1 1
25.84193 20 1 1.05 0 0;
31 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
32 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0 ;
33 1 0.1053 0.05099913 0 0 1 1
25.84193 20 1 1.05 0 0;
34 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
35 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
36 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
37 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
38 1 0.0711 0.03443531 0 0 1 1
25.84193 20 1 1.05 0 0;
39 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
40 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
41 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
42 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
43 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
44 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
45 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
46 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
47 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
48 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
49 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
50 1 0.0162 0.00784602 0 0 1 1
25.84193 20 1 1.05 0 0;
51 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
52 3 0.00001 0.00001 0 0 1 1.03 0 20
1 1.05 0 0 ;
53 1 0.1008 0.04881968 0 0 1 1
25.84193 20 1 1.05 0 0;
54 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0 ;
55 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
56 1 0.0324 0.01569204 0 0 1 1
25.84193 20 1 1.05 0 0;
57 1 0.0999 0.04838379 0 0 1 1
25.84193 20 1 1.05 0 0;
58 1 0.1323 0.06407583 0 0 1 1
25.84193 20 1 1.05 0 0;
59 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
60 1 0.0027 0.00130767 0 0 1 1
25.84193 20 1 1.05 0 0;
61 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0 ;
62 1 0.16875 0.08172938 0 0 1 1
25.84193 20 1 1.05 0 0;
63 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0 ;
64 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
65 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
66 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
67 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
68 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
69 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
70 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
71 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
72 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
73 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
74 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
75 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
76 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
77 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
78 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
79 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
80 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
81 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
82 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
83 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
84 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
85 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
86 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
87 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
88 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
89 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
90 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
91 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
92 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
93 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
];
%% generator data
mpc.gen= [
52 3 3 999999 0 10 1 999999 999999 0 0 0 0 0 0 0 0 0;
];
%% Branch data
mpc.branch= [
1 2 0.00079091481 0.0005140946265000001 0 8 8 8 0 0 1 -360 360 ;
4 1 2.1957460133333335e-5 9.033132333333336e-6 0 6 6 6 0 0 1 -360 360 ;
2 3 0.0005720009866666666 0.0003718006413333333 0 8 8 8 0 0 1 -360 360 ;
3 63 1.085185e-6 7.0537025e-7 0 8 8 8 0 0 1 -360 360 ;
65 19 0.0004964366608499999 0.00018061624875 0 2 2 2 0 0 1 -360 360 ;
18 53 0.0004413146533333334 0.0002868545246666667 0 8 8 8 0 0 1 -360 360 ;
51 18 5.060353333333334e-6 3.289229666666667e-6 0 8 8 8 0 0 1 -360 360 ;
26 68 0.0034558943249500002 0.0027615593920833335 0 4 4 4 0 0 1 -360 360 ;
35 23 0.0007813979533333334 0.00032146118333333337 0 6 6 6 0 0 1 -360 360 ;
69 70 0.0057047391906 0.004558581545 0 4 4 4 0 0 1 -360 360 ;
71 26 0.0035307791776500006 0.0028213988862499997 0 4 4 4 0 0 1 -360 360 ;
35 72 0.0018416064926000003 0.0014716033616666666 0 4 4 4 0 0 1 -360 360 ;
73 35 0.0023422466793499995 0.0018716583054166664 0 4 4 4 0 0 1 -360 360 ;
74 6 0.0017254742886333334 0.0006277713108333334 0 2 2 2 0 0 1 -360 360 ;
56 34 0.0016937272865500002 0.0006162209462499999 0 2 2 2 0 0 1 -360 360 ;
12 75 0.0031468124616499996 0.0011448901887499998 0 2 2 2 0 0 1 -360 360 ;
36 30 4.3935405000000013e-5 2.855801325e-5 0 8 8 8 0 0 1 -360 360 ;
53 31 0.0005246659250000001 0.00034103285125 0 8 8 8 0 0 1 -360 360 ;
48 27 5.1925578200000006e-5 2.13617885e-5 0 6 6 6 0 0 1 -360 360 ;
72 33 0.0044105512957000005 0.0035244131358333334 0 4 4 4 0 0 1 -360 360 ;
73 50 0.00036866939974999994 0.0002945988354166666 0 4 4 4 0 0 1 -360 360 ;
76 47 0.0010931357150500002 0.0008735102745833333 0 4 4 4 0 0 1 -360 360 ;
76 73 0.0009360721964 0.0007480028966666667 0 4 4 4 0 0 1 -360 360 ;
19 5 0.0008050901360666667 0.00029291221166666666 0 2 2 2 0 0 1 -360 360 ;
37 38 0.0000001 0.0000001 0 2 2 2 0 0 1 -360 360 ;
38 77 0.0000001 0.0000001 0 2 2 2 0 0 1 -360 360 ;
39 37 0.0000001 0.0000001 0 2 2 2 0 0 1 -360 360 ;
77 33 0.00032999706500000007 0.00021449809225000004 0 8 8 8 0 0 1 -360 360 ;
34 32 0.00022064490123333333 9.077163658333333e-5 0 6 6 6 0 0 1 -360 360 ;
34 40 0.0008979513868 0.00032669748999999997 0 2 2 2 0 0 1 -360 360 ;
40 72 0.004124282800516667 0.001500518690416667 0 2 2 2 0 0 1 -360 360 ;
78 7 0.0011785434983833333 0.0004287840170833333 0 2 2 2 0 0 1 -360 360 ;
78 41 0.0009989464838166668 0.00036344206791666666 0 2 2 2 0 0 1 -360 360 ;
64 9 0.00060720191805 0.00022091545874999998 0 2 2 2 0 0 1 -360 360 ;
64 42 0.00018438561006666666 6.708416166666666e-5 0 2 2 2 0 0 1 -360 360 ;
80 52 3.309269183333333e-5 1.3614082083333334e-5 0 6 6 6 0 0 1 -360 360 ;
81 13 0.0004699604888333333 0.0001709835458333333 0 2 2 2 0 0 1 -360 360 ;
81 44 0.0034228753730166663 0.0012453288779166664 0 2 2 2 0 0 1 -360 360 ;
82 21 0.0018645889964000002 0.00067838477 0 2 2 2 0 0 1 -360 360 ;
82 81 0.004525359672166667 0.0016464406291666663 0 2 2 2 0 0 1 -360 360 ;
14 82 0.0017651154466833334 0.0006421937695833334 0 2 2 2 0 0 1 -360 360 ;
83 46 0.00026861798535 9.773003625e-5 0 2 2 2 0 0 1 -360 360 ;
16 83 0.0033262158496666664 0.0012101616916666666 0 2 2 2 0 0 1 -360 360 ;
66 55 0.0026072950037333334 0.0009486000533333333 0 2 2 2 0 0 1 -360 360 ;
66 47 0.0012992917181833335 0.0004727152820833334 0 2 2 2 0 0 1 -360 360 ;
84 85 0.0044402999367999995 0.00161549374 0 2 2 2 0 0 1 -360 360 ;
50 84 0.0010149622526333332 0.0003692690108333333 0 2 2 2 0 0 1 -360 360 ;
86 28 3.911652333333333e-6 1.4231583333333332e-6 0 2 2 2 0 0 1 -360 360 ;
70 46 0.0004901945154166666 0.00017834519791666664 0 2 2 2 0 0 1 -360 360 ;
69 44 0.001162468045483333 0.00042293535958333326 0 2 2 2 0 0 1 -360 360 ;
69 45 0.0010746510497666667 0.00039098530916666666 0 2 2 2 0 0 1 -360 360 ;
68 42 0.0010267242739 0.00037354833249999994 0 2 2 2 0 0 1 -360 360 ;
87 71 0.0007138366438666667 0.00025971187666666665 0 2 2 2 0 0 1 -360 360 ;
86 87 0.0028728646676166665 0.0010452210329166664 0 2 2 2 0 0 1 -360 360 ;
29 86 0.0023903865145333335 0.0008696832433333333 0 2 2 2 0 0 1 -360 360 ;
88 10 1.6498781e-6 6.002674999999998e-7 0 2 2 2 0 0 1 -360 360 ;
23 24 0.0007645104867333334 0.00031451380783333336 0 6 6 6 0 0 1 -360 360 ;
67 88 0.0036057164840666667 0.0013118511116666667 0 2 2 2 0 0 1 -360 360 ;
58 87 2.6781052666666663e-6 1.1017521666666667e-6 0 6 6 6 0 0 1 -360 360 ;
59 58 2.3557649899999998e-5 9.69143825e-6 0 6 6 6 0 0 1 -360 360 ;
89 41 0.0006974911335166666 0.00025376496541666666 0 2 2 2 0 0 1 -360 360 ;
89 67 0.0029338324597833333 0.0010674026620833333 0 2 2 2 0 0 1 -360 360 ;
43 89 0.0010194760047333331 0.0003709112283333333 0 2 2 2 0 0 1 -360 360 ;
60 61 8.525732666666668e-6 3.507421666666667e-6 0 6 6 6 0 0 1 -360 360 ;
60 62 0.00017986328333333335 0.00011691113416666666 0 8 8 8 0 0 1 -360 360 ;
54 62 0.0008368968083333333 0.0005439829254166667 0 8 8 8 0 0 1 -360 360 ;
54 90 3.958368333333333e-6 2.5729394166666664e-6 0 8 8 8 0 0 1 -360 360 ;
17 90 0.004442779884616666 0.0016163960079166664 0 2 2 2 0 0 1 -360 360 ;
90 83 0.0005929933864833333 0.0002157460345833333 0 2 2 2 0 0 1 -360 360 ;
74 43 0.0031021882478166666 0.0011286547679166667 0 2 2 2 0 0 1 -360 360 ;
65 74 0.0017193581812333333 0.0006255461158333333 0 2 2 2 0 0 1 -360 360 ;
91 57 0.0011836329674000002 0.00043063569499999995 0 2 2 2 0 0 1 -360 360 ;
91 65 0.0005486792172333333 0.00019962341583333333 0 2 2 2 0 0 1 -360 360 ;
11 91 0.0016498279156666666 0.0006002492416666667 0 2 2 2 0 0 1 -360 360 ;
75 21 0.0007215457030166667 0.00026251662791666665 0 2 2 2 0 0 1 -360 360 ;
75 48 0.0029021376381833334 0.0010558712820833331 0 2 2 2 0 0 1 -360 360 ;
15 63 0.0009936733289166667 0.0003615235604166667 0 2 2 2 0 0 1 -360 360 ;
63 45 0.0010615525539666667 0.00038621974416666664 0 2 2 2 0 0 1 -360 360 ;
19 20 0.0007126708825 0.0004632360925 0 8.00000000 8.00000000 8.00000000 0 0 1
-360 360 ;
64 88 0.0086973985 0.0031643341 0 2.00000000 2.00000000 2.00000000 0 0 1 -360
360 ;
66 49 0.004808182516666666 0.0017493387083333332 0 2.00000000 2.00000000
2.00000000 0 0 1 -360 360 ;
67 8 0.003954792516666667 0.0014388538916666666 0 2.00000000 2.00000000
2.00000000 0 0 1 -360 360 ;
68 25 0.00583548645 0.004663059608333334 0 4.00000000 4.00000000 4.00000000 0 0
1 -360 360 ;
69 25 0.0017030360833333335 0.0013608735083333332 0 4.00000000 4.00000000
4.00000000 0 0 1 -360 360 ;
70 22 0.000738570715 0.0005901820958333334 0 4.00000000 4.00000000 4.00000000 0
0 1 -360 360 ;
71 80 0.005592975516666667 0.004469272991666666 0 4.00000000 4.00000000
4.00000000 0 0 1 -360 360 ;
76 22 0.0008360097833333333 0.0006680443883333334 0 4.00000000 4.00000000
4.00000000 0 0 1 -360 360 ;
78 79 0.001037387175 0.00037742778666666665 0 2.00000000 2.00000000 2.00000000
0 0 1 -360 360 ;
80 31 0.0011012909333333334 0.0008019279691666666 0 4.00000000 4.00000000
4.00000000 0 0 1 -360 360 ;
];
function mpc = radial_case_null
%CASE30 Power flow data for 30 bus, 6 generator case.
% Please see CASEFORMAT for details on the case file format.
%
% Based on data from ...
% Alsac, O. & Stott, B., "Optimal Load Flow with Steady State Security",
% IEEE Transactions on Power Apparatus and Systems, Vol. PAS 93, No. 3,
% 1974, pp. 745-751.
% ... with branch parameters rounded to nearest 0.01, shunt values divided
% by 100 and shunt on bus 10 moved to bus 5, load at bus 5 zeroed out.
% Generator locations, costs and limits and bus areas were taken from ...
% Ferrero, R.W., Shahidehpour, S.M., Ramesh, V.C., "Transaction analysis
% in deregulated power systems using game theory", IEEE Transactions on
% Power Systems, Vol. 12, No. 3, Aug 1997, pp. 1340-1347.
% Generator Q limits were derived from Alsac & Stott, using their Pmax
% capacities. V limits and line |S| limits taken from Alsac & Stott.
% MATPOWER
% $Id: case30.m,v 1.12 2010/03/10 18:08:13 ray Exp $
%% MATPOWER Case Format : Version 2
mpc.version = '2';
mpc.baseMVA = 10;
%% Bus Data
mpc.bus= [
1 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
2 1 0 0 0 0 1 1 0 20 1 1.05 0 0 ;
3 1 0 0 0 0 1 1 0 20 1 1.05 0 0 ;
4 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
5 1 0.03375000 0.01634587 0 0 1 1.0 25.84193 20 1 1.05 0 0;
6 1 0.09180000 0.04446078 0 0 1 1.0 25.84193 20 1 1.05 0 0;
7 1 0.05940000 0.02876874 0 0 1 1.0 25.84193 20 1 1.05 0 0;
8 1 0.08280000 0.04010188 0 0 1 1.0 25.84193 20 1 1.05 0 0;
9 1 0.10440000 0.05056324 0 0 1 1.0 25.84193 20 1 1.05 0 0;
10 1 0.03510000 0.01699971 0 0 1 1.0 25.84193 20 1 1.05 0 0;
11 1 0.01890000 0.00915369 0 0 1 1.0 25.84193 20 1 1.05 0 0;
12 1 0.05670000 0.02746107 0 0 1 1.0 25.84193 20 1 1.05 0 0;
13 1 0.09810000 0.04751201 0 0 1 1.0 25.84193 20 1 1.05 0 0;
14 1 0.02430000 0.01176903 0 0 1 1.0 25.84193 20 1 1.05 0 0;
15 1 0.16020000 0.07758842 0 0 1 1.0 25.84193 20 1 1.05 0 0;
16 1 0.01620000 0.00784602 0 0 1 1.0 25.84193 20 1 1.05 0 0;
17 1 0.06660000 0.03225586 0 0 1 1.0 25.84193 20 1 1.05 0 0;
18 1 0.02700000 0.01307670 0 0 1 1.0 25.84193 20 1 1.05 0 0;
19 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
20 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
21 1 0.08280000 0.04010188 0 0 1 1.0 25.84193 20 1 1.05 0 0;
22 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
23 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
24 1 0 0 0 0 1 1 0 20 1 1.05 0 0 ;
25 1 0.01080000 0.00523068 0 0 1 1.0 25.84193 20 1 1.05 0 0;
26 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
27 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
28 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
29 1 0.05130000 0.02484573 0 0 1 1.0 25.84193 20 1 1.05 0 0;
30 1 0.02160000 0.01046136 0 0 1 1.0 25.84193 20 1 1.05 0 0;
31 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
32 1 0 0 0 0 1 1 0 20 1 1.05 0 0 ;
33 1 0.10530000 0.05099913 0 0 1 1.0 25.84193 20 1 1.05 0 0;
34 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
35 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
36 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
37 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
38 1 0.07110000 0.03443531 0 0 1 1.0 25.84193 20 1 1.05 0 0;
39 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
40 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
41 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
42 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
43 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
44 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
45 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
46 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
47 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
48 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
49 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
50 1 0.01620000 0.00784602 0 0 1 1.0 25.84193 20 1 1.05 0 0;
51 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
52 3 0 0 0 0 1 1.03 0 20 1 1.05 0 0 ;
53 1 0.10080000 0.04881968 0 0 1 1.0 25.84193 20 1 1.05 0 0;
54 1 0 0 0 0 1 1 0 20 1 1.05 0 0 ;
55 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
56 1 0.03240000 0.01569204 0 0 1 1.0 25.84193 20 1 1.05 0 0;
57 1 0.09990000 0.04838379 0 0 1 1.0 25.84193 20 1 1.05 0 0;
58 1 0.13230000 0.06407583 0 0 1 1.0 25.84193 20 1 1.05 0 0;
59 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
60 1 0.00270000 0.00130767 0 0 1 1.0 25.84193 20 1 1.05 0 0;
61 1 0 0 0 0 1 1 0 20 1 1.05 0 0 ;
62 1 0.16875000 0.08172938 0 0 1 1.0 25.84193 20 1 1.05 0 0;
63 1 0 0 0 0 1 1 0 20 1 1.05 0 0 ;
64 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
65 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
66 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
67 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
68 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
69 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
70 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
71 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
72 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
73 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
74 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
75 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
76 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
77 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
78 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
79 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
80 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
81 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
82 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
83 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
84 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
85 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
86 1 0 0 0 0 1 1 25.8 20 1 1.05 0 0 ;
87 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
88 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
89 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
90 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
91 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
92 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
93 1 0 0 0 0 1 1 0 20 1 1.05 0 0;
];
%% generator data
mpc.gen= [
52 0 0 999999 0 10 1 999999 999999 0 0 0 0 0 0 0 0 0;
];
%% Branch data
mpc.branch= [
1 2 0.00079091481 0.0005140946265000001 0 8 8 8 0 0 1 -360 360 ;
4 1 2.1957460133333335e-5 9.033132333333336e-6 0 6 6 6 0 0 1 -360 360 ;
2 3 0.0005720009866666666 0.0003718006413333333 0 8 8 8 0 0 1 -360 360 ;
3 63 1.085185e-6 7.0537025e-7 0 8 8 8 0 0 1 -360 360 ;
65 19 0.0004964366608499999 0.00018061624875 0 2 2 2 0 0 1 -360 360 ;
18 53 0.0004413146533333334 0.0002868545246666667 0 8 8 8 0 0 1 -360 360 ;
51 18 5.060353333333334e-6 3.289229666666667e-6 0 8 8 8 0 0 1 -360 360 ;
26 68 0.0034558943249500002 0.0027615593920833335 0 4 4 4 0 0 1 -360 360 ;
35 23 0.0007813979533333334 0.00032146118333333337 0 6 6 6 0 0 1 -360 360 ;
69 70 0.0057047391906 0.004558581545 0 4 4 4 0 0 1 -360 360 ;
71 26 0.0035307791776500006 0.0028213988862499997 0 4 4 4 0 0 1 -360 360 ;
35 72 0.0018416064926000003 0.0014716033616666666 0 4 4 4 0 0 1 -360 360 ;
73 35 0.0023422466793499995 0.0018716583054166664 0 4 4 4 0 0 1 -360 360 ;
74 6 0.0017254742886333334 0.0006277713108333334 0 2 2 2 0 0 1 -360 360 ;
56 34 0.0016937272865500002 0.0006162209462499999 0 2 2 2 0 0 1 -360 360 ;
12 75 0.0031468124616499996 0.0011448901887499998 0 2 2 2 0 0 1 -360 360 ;
36 30 4.3935405000000013e-5 2.855801325e-5 0 8 8 8 0 0 1 -360 360 ;
53 31 0.0005246659250000001 0.00034103285125 0 8 8 8 0 0 1 -360 360 ;
48 27 5.1925578200000006e-5 2.13617885e-5 0 6 6 6 0 0 1 -360 360 ;
72 33 0.0044105512957000005 0.0035244131358333334 0 4 4 4 0 0 1 -360 360 ;
73 50 0.00036866939974999994 0.0002945988354166666 0 4 4 4 0 0 1 -360 360 ;
76 47 0.0010931357150500002 0.0008735102745833333 0 4 4 4 0 0 1 -360 360 ;
76 73 0.0009360721964 0.0007480028966666667 0 4 4 4 0 0 1 -360 360 ;
19 5 0.0008050901360666667 0.00029291221166666666 0 2 2 2 0 0 1 -360 360 ;
37 38 0.0000001 0.0000001 0 2 2 2 0 0 1 -360 360 ;
38 77 0.0000001 0.0000001 0 2 2 2 0 0 1 -360 360 ;
39 37 0.0000001 0.0000001 0 2 2 2 0 0 1 -360 360 ;
77 33 0.00032999706500000007 0.00021449809225000004 0 8 8 8 0 0 1 -360 360 ;
34 32 0.00022064490123333333 9.077163658333333e-5 0 6 6 6 0 0 1 -360 360 ;
34 40 0.0008979513868 0.00032669748999999997 0 2 2 2 0 0 1 -360 360 ;
40 72 0.004124282800516667 0.001500518690416667 0 2 2 2 0 0 1 -360 360 ;
78 7 0.0011785434983833333 0.0004287840170833333 0 2 2 2 0 0 1 -360 360 ;
78 41 0.0009989464838166668 0.00036344206791666666 0 2 2 2 0 0 1 -360 360 ;
64 9 0.00060720191805 0.00022091545874999998 0 2 2 2 0 0 1 -360 360 ;
64 42 0.00018438561006666666 6.708416166666666e-5 0 2 2 2 0 0 1 -360 360 ;
80 52 3.309269183333333e-5 1.3614082083333334e-5 0 6 6 6 0 0 1 -360 360 ;
81 13 0.0004699604888333333 0.0001709835458333333 0 2 2 2 0 0 1 -360 360 ;
81 44 0.0034228753730166663 0.0012453288779166664 0 2 2 2 0 0 1 -360 360 ;
82 21 0.0018645889964000002 0.00067838477 0 2 2 2 0 0 1 -360 360 ;
82 81 0.004525359672166667 0.0016464406291666663 0 2 2 2 0 0 1 -360 360 ;
14 82 0.0017651154466833334 0.0006421937695833334 0 2 2 2 0 0 1 -360 360 ;
83 46 0.00026861798535 9.773003625e-5 0 2 2 2 0 0 1 -360 360 ;
16 83 0.0033262158496666664 0.0012101616916666666 0 2 2 2 0 0 1 -360 360 ;
66 55 0.0026072950037333334 0.0009486000533333333 0 2 2 2 0 0 1 -360 360 ;
66 47 0.0012992917181833335 0.0004727152820833334 0 2 2 2 0 0 1 -360 360 ;
84 85 0.0044402999367999995 0.00161549374 0 2 2 2 0 0 1 -360 360 ;
50 84 0.0010149622526333332 0.0003692690108333333 0 2 2 2 0 0 1 -360 360 ;
86 28 3.911652333333333e-6 1.4231583333333332e-6 0 2 2 2 0 0 1 -360 360 ;
70 46 0.0004901945154166666 0.00017834519791666664 0 2 2 2 0 0 1 -360 360 ;
69 44 0.001162468045483333 0.00042293535958333326 0 2 2 2 0 0 1 -360 360 ;
69 45 0.0010746510497666667 0.00039098530916666666 0 2 2 2 0 0 1 -360 360 ;
68 42 0.0010267242739 0.00037354833249999994 0 2 2 2 0 0 1 -360 360 ;
87 71 0.0007138366438666667 0.00025971187666666665 0 2 2 2 0 0 1 -360 360 ;
86 87 0.0028728646676166665 0.0010452210329166664 0 2 2 2 0 0 1 -360 360 ;
29 86 0.0023903865145333335 0.0008696832433333333 0 2 2 2 0 0 1 -360 360 ;
88 10 1.6498781e-6 6.002674999999998e-7 0 2 2 2 0 0 1 -360 360 ;
23 24 0.0007645104867333334 0.00031451380783333336 0 6 6 6 0 0 1 -360 360 ;
67 88 0.0036057164840666667 0.0013118511116666667 0 2 2 2 0 0 1 -360 360 ;
58 87 2.6781052666666663e-6 1.1017521666666667e-6 0 6 6 6 0 0 1 -360 360 ;
59 58 2.3557649899999998e-5 9.69143825e-6 0 6 6 6 0 0 1 -360 360 ;
89 41 0.0006974911335166666 0.00025376496541666666 0 2 2 2 0 0 1 -360 360 ;
89 67 0.0029338324597833333 0.0010674026620833333 0 2 2 2 0 0 1 -360 360 ;
43 89 0.0010194760047333331 0.0003709112283333333 0 2 2 2 0 0 1 -360 360 ;
60 61 8.525732666666668e-6 3.507421666666667e-6 0 6 6 6 0 0 1 -360 360 ;
60 62 0.00017986328333333335 0.00011691113416666666 0 8 8 8 0 0 1 -360 360 ;
54 62 0.0008368968083333333 0.0005439829254166667 0 8 8 8 0 0 1 -360 360 ;
54 90 3.958368333333333e-6 2.5729394166666664e-6 0 8 8 8 0 0 1 -360 360 ;
17 90 0.004442779884616666 0.0016163960079166664 0 2 2 2 0 0 1 -360 360 ;
90 83 0.0005929933864833333 0.0002157460345833333 0 2 2 2 0 0 1 -360 360 ;
74 43 0.0031021882478166666 0.0011286547679166667 0 2 2 2 0 0 1 -360 360 ;
65 74 0.0017193581812333333 0.0006255461158333333 0 2 2 2 0 0 1 -360 360 ;
91 57 0.0011836329674000002 0.00043063569499999995 0 2 2 2 0 0 1 -360 360 ;
91 65 0.0005486792172333333 0.00019962341583333333 0 2 2 2 0 0 1 -360 360 ;
11 91 0.0016498279156666666 0.0006002492416666667 0 2 2 2 0 0 1 -360 360 ;
75 21 0.0007215457030166667 0.00026251662791666665 0 2 2 2 0 0 1 -360 360 ;
75 48 0.0029021376381833334 0.0010558712820833331 0 2 2 2 0 0 1 -360 360 ;
15 63 0.0009936733289166667 0.0003615235604166667 0 2 2 2 0 0 1 -360 360 ;
63 45 0.0010615525539666667 0.00038621974416666664 0 2 2 2 0 0 1 -360 360 ;
19 20 0.0007126708825 0.0004632360925 0 8.00000000 8.00000000 8.00000000 0 0 1
-360 360 ;
64 88 0.0086973985 0.0031643341 0 2.00000000 2.00000000 2.00000000 0 0 1 -360
360 ;
66 49 0.004808182516666666 0.0017493387083333332 0 2.00000000 2.00000000
2.00000000 0 0 1 -360 360 ;
67 8 0.003954792516666667 0.0014388538916666666 0 2.00000000 2.00000000
2.00000000 0 0 1 -360 360 ;
68 25 0.00583548645 0.004663059608333334 0 4.00000000 4.00000000 4.00000000 0 0
1 -360 360 ;
69 25 0.0017030360833333335 0.0013608735083333332 0 4.00000000 4.00000000
4.00000000 0 0 1 -360 360 ;
70 22 0.000738570715 0.0005901820958333334 0 4.00000000 4.00000000 4.00000000 0
0 1 -360 360 ;
71 80 0.005592975516666667 0.004469272991666666 0 4.00000000 4.00000000
4.00000000 0 0 1 -360 360 ;
76 22 0.0008360097833333333 0.0006680443883333334 0 4.00000000 4.00000000
4.00000000 0 0 1 -360 360 ;
78 79 0.001037387175 0.00037742778666666665 0 2.00000000 2.00000000 2.00000000
0 0 1 -360 360 ;
80 31 0.0011012909333333334 0.0008019279691666666 0 4.00000000 4.00000000
4.00000000 0 0 1 -360 360 ;
];
