Dear Nitya I would do two things:
1. I would replace all commas "," by points "." in the data file. 2. I would try to draw the one line diagram of the network based on the data, to check whether the network is all connected, that is, that there are no isolated buses or group of buses, Good luck, Carlos Em seg, 8 de abr de 2019 às 14:42, Nitya Kirana <[email protected]> escreveu: > 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 > > --
