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
>
>

--

Reply via email to