Dear all,
I have run the power flow of 64 radial test bus system using Newton Raphson
in Matpower 4.1. It shows the result of power losses but it doesn't
converge, whenever I run the matlab, the result is always different, and I
don't know what's the problem. Can anyone help me? because this is my final
project. It shows "Warning: Matrix is singular to working precision."
So, here is my data :
bus data
bus type Pd Qd Gs Bs area Vm
Va baseKV zone Vmax Vmin
1 3 0 0 0 0 1 1 0 20 1 1 0,9
2 3 0 0 0 0 1 1 0 20 1 1 0,9
3 3 0 0 0 0 1 1 0 20 1 1 0,9
4 1 0,0136 0,016256 0 0 1 1 0 20 1 1 0,9
5 1 0,0204 0,024384 0 0 1 1 0 20 1 1 0,9
6 1 0 0 0 0 1 1 0 20 1 1 0,9
7 1 0,07684 0,091846 0 0 1 1 0 20 1 1 0,9
8 1 0,02856 0,034138 0 0 1 1 0 20 1 1 0,9
9 1 0,03808 0,045517 0 0 1 1 0 20 1 1 0,9
10 1 0,0238 0,028448 0 0 1 1 0 20 1 1 0,9
11 1 0,01768 0,021133 0 0 1 1 0 20 1 1 0,9
12 1 0,017 0,02032 0 0 1 1 0 20 1 1 0,9
13 1 0,022304 0,02666 0 0 1 1 0 20 1 1 0,9
14 1 0,0204 0,002438 0 0 1 1 0 20 1 1 0,9
15 1 0,05848 0,069901 0 0 1 1 0 20 1 1 0,9
16 1 0,26384 0,315366 0 0 1 1 0 20 1 1 0,9
17 1 0,03774 0,04511 0 0 1 1 0 20 1 1 0,9
18 1 0,0204 0,002438 0 0 1 1 0 20 1 1 0,9
19 1 0,05576 0,06665 0 0 1 1 0 20 1 1 0,9
20 1 0,0136 0,016256 0 0 1 1 0 20 1 1 0,9
21 1 0,0306 0,016256 0 0 1 1 0 20 1 1 0,9
22 1 0,0204 0,024384 0 0 1 1 0 20 1 1 0,9
23 1 0,0272 0,032512 0 0 1 1 0 20 1 1 0,9
24 1 0,11764 0,140614 0 0 1 1 0 20 1 1 0,9
25 1 0,0068 0,008128 0 0 1 1 0 20 1 1 0,9
26 1 0,0136 0,016256 0 0 1 1 0 20 1 1 0,9
27 1 0,0238 0,028448 0 0 1 1 0 20 1 1 0,9
28 1 0,112948 0,135006 0 0 1 1 0 20 1 1 0,9
29 1 0,02176 0,02601 0 0 1 1 0 20 1 1 0,9
30 1 0 0 0 0 1 1 0 20 1 1 0,9
31 1 0,01768 0,021133 0 0 1 1 0 20 1 1 0,9
32 1 0,05168 0,061773 0 0 1 1 0 20 1 1 0,9
33 1 0,03468 0,041453 0 0 1 1 0 20 1 1 0,9
34 1 0,02856 0,034138 0 0 1 1 0 20 1 1 0,9
35 1 0,0272 0,032512 0 0 1 1 0 20 1 1 0,9
36 1 0,01768 0,021133 0 0 1 1 0 20 1 1 0,9
37 1 0,07616 0,091034 0 0 1 1 0 20 1 1 0,9
38 1 0,0816 0,097536 0 0 1 1 0 20 1 1 0,9
39 1 0,06868 0,082093 0 0 1 1 0 20 1 1 0,9
40 1 0,02448 0,029261 0 0 1 1 0 20 1 1 0,9
41 1 0,0408 0,048768 0 0 1 1 0 20 1 1 0,9
42 1 0,03536 0,042266 0 0 1 1 0 20 1 1 0,9
43 1 0,03536 0,042266 0 0 1 1 0 20 1 1 0,9
44 1 0,03536 0,042266 0 0 1 1 0 20 1 1 0,9
45 1 0,02856 0,034138 0 0 1 1 0 20 1 1 0,9
46 1 0,0272 0,032512 0 0 1 1 0 20 1 1 0,9
47 1 0,03672 0,043891 0 0 1 1 0 20 1 1 0,9
48 1 0,0136 0,016256 0 0 1 1 0 20 1 1 0,9
49 1 0,0306 0,036576 0 0 1 1 0 20 1 1 0,9
50 1 0,01088 0,013005 0 0 1 1 0 20 1 1 0,9
51 1 0,06324 0,07559 0 0 1 1 0 20 1 1 0,9
52 1 0,04488 0,053645 0 0 1 1 0 20 1 1 0,9
53 1 0,01088 0,013005 0 0 1 1 0 20 1 1 0,9
54 1 0,06528 0,013005 0 0 1 1 0 20 1 1 0,9
55 1 0,06528 0,078029 0 0 1 1 0 20 1 1 0,9
56 1 0,0272 0,032512 0 0 1 1 0 20 1 1 0,9
57 1 0,01768 0,021133 0 0 1 1 0 20 1 1 0,9
58 1 0,12886 0,154026 0 0 1 1 0 20 1 1 0,9
59 1 0,18836 0,225146 0 0 1 1 0 20 1 1 0,9
60 1 0,0068 0,008128 0 0 1 1 0 20 1 1 0,9
61 1 0,0068 0,008128 0 0 1 1 0 20 1 1 0,9
62 1 0,03128 0,037389 0 0 1 1 0 20 1 1 0,9
63 1 0 0 0 0 1 1 0 20 1 1 0,9
64 1 0,0612 0,073152 0 0 1 1 0 20 1 1 0,9
*generator databus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2
Qc1min Qc1max Qc2min Qc2max ramp_agc ramp_10 ramp_30 ramp_q apf1 0 0 0 0
1.00 100 1 0 0 0 0 0 0 0 0 0 0 0 0 0;2 0 0 0 0 1.00 100 1
0 0 0 0 0 0 0 0 0 0 0 0 0;3 0 0 0 0 1.00 100 1 0 0 0 0 0
0 0 0 0 0 0 0 0;*
*branch data fbus tbus r x b rateA rateB rateC
ratio angle status angmin angmax*
1 4 0,242144 0,3752 0 0 0 0 0 0 1 -360 360
4 5 0,041078 0,06365 0 0 0 0 0 0 1 -360 360
5 6 0,047564 0,0737 0 0 0 0 0 0 1 -360 360
6 7 0,045402 0,07035 0 0 0 0 0 0 1 -360 360
7 8 0,04324 0,067 0 0 0 0 0 0 1 -360 360
8 9 0,05405 0,08375 0 0 0 0 0 0 1 -360 360
6 10 0,082156 0,1273 0 0 0 0 0 0 1 -360 360
10 11 0,051888 0,0804 0 0 0 0 0 0 1 -360 360
11 12 0,036754 0,05695 0 0 0 0 0 0 1 -360 360
12 13 0,041078 0,06365 0 0 0 0 0 0 1 -360 360
13 14 0,062698 0,09715 0 0 0 0 0 0 1 -360 360
14 15 0,034592 0,0536 0 0 0 0 0 0 1 -360 360
15 16 0,069184 0,1072 0 0 0 0 0 0 1 -360 360
16 17 0,051888 0,0804 0 0 0 0 0 0 1 -360 360
2 18 0,25944 0,402 0 0 0 0 0 0 1 -360 360
18 19 0,058374 0,09045 0 0 0 0 0 0 1 -360 360
19 20 0,071346 0,11055 0 0 0 0 0 0 1 -360 360
20 21 0,036754 0,05695 0 0 0 0 0 0 1 -360 360
21 22 0,025944 0,0402 0 0 0 0 0 0 1 -360 360
22 23 0,025944 0,0402 0 0 0 0 0 0 1 -360 360
23 24 0,06486 0,1005 0 0 0 0 0 0 1 -360 360
24 60 0,051888 0,0804 0 0 0 0 0 0 1 -360 360
24 25 0,036754 0,05695 0 0 0 0 0 0 1 -360 360
25 26 0,034592 0,402 0 0 0 0 0 0 1 -360 360
26 27 0,019458 0,03015 0 0 0 0 0 0 1 -360 360
27 28 0,05405 0,08375 0 0 0 0 0 0 1 -360 360
28 29 0,060536 0,0938 0 0 0 0 0 0 1 -360 360
28 62 0,051888 0,0804 0 0 0 0 0 0 1 -360 360
29 30 0,034592 0,0536 0 0 0 0 0 0 1 -360 360
30 31 0,02162 0,0335 0 0 0 0 0 0 1 -360 360
31 32 0,036754 0,05695 0 0 0 0 0 0 1 -360 360
32 33 0,034592 0,402 0 0 0 0 0 0 1 -360 360
33 34 0,067022 0,10385 0 0 0 0 0 0 1 -360 360
34 35 0,036754 0,05695 0 0 0 0 0 0 1 -360 360
35 36 0,041078 0,06365 0 0 0 0 0 0 1 -360 360
36 37 0,015134 0,02345 0 0 0 0 0 0 1 -360 360
37 38 0,02162 0,0335 0 0 0 0 0 0 1 -360 360
38 39 0,05405 0,08375 0 0 0 0 0 0 1 -360 360
39 40 0,03243 0,05025 0 0 0 0 0 0 1 -360 360
3 41 0,025944 0,0402 0 0 0 0 0 0 1 -360 360
41 42 0,045402 0,07035 0 0 0 0 0 0 1 -360 360
42 43 0,045402 0,07035 0 0 0 0 0 0 1 -360 360
43 44 0,05405 0,08375 0 0 0 0 0 0 1 -360 360
44 45 0,060536 0,0938 0 0 0 0 0 0 1 -360 360
45 46 0,049726 0,07705 0 0 0 0 0 0 1 -360 360
46 47 0,045402 0,07035 0 0 0 0 0 0 1 -360 360
47 48 0,049726 0,07705 0 0 0 0 0 0 1 -360 360
48 49 0,034592 0,0402 0 0 0 0 0 0 1 -360 360
49 50 0,069184 0,1072 0 0 0 0 0 0 1 -360 360
50 51 0,095128 0,1474 0 0 0 0 0 0 1 -360 360
51 52 0,088642 0,13735 0 0 0 0 0 0 1 -360 360
52 61 0,034592 0,0402 0 0 0 0 0 0 1 -360 360
52 53 0,101614 0,15745 0 0 0 0 0 0 1 -360 360
53 54 0,047564 0,0737 0 0 0 0 0 0 1 -360 360
54 55 0,017296 0,0268 0 0 0 0 0 0 1 -360 360
55 56 0,006486 0,01005 0 0 0 0 0 0 1 -360 360
56 57 0,019458 0,03015 0 0 0 0 0 0 1 -360 360
57 63 0,11891 0,18425 0 0 0 0 0 0 1 -360 360
57 58 0,155664 0,2412 0 0 0 0 0 0 1 -360 360
58 59 0,038916 0,0603 0 0 0 0 0 0 1 -360 360
59 64 0,025944 0,0402 0 0 0 0 0 0 1 -360 360
6 60 0,03243 0,05025 0 0 0 0 0 0 1 -360 360
24 61 0,062698 0,09715 0 0 0 0 0 0 1 -360 360
11 62 0,030268 0,0469 0 0 0 0 0 0 1 -360 360
30 63 0,002162 0,00335 0 0 0 0 0 0 1 -360 360
39 64 0,006486 0,01005 0 0 0 0 0 0 1 -360 360
*Best Regards,*
*Putu Diah Nitya Kirana,*
Assistant and Researcher at Electric Energy Conversion Laboratory B101
Electrical Engineering Department
Faculty of Electrical Technology
Institut Teknologi Sepuluh Nopember
Surabaya, Indonesia
+62 8170 3178 248
