Dear Carlos A. Castro,

Thank you before. I've done the replacement and the buses are all correct,
because I draw the single line diagram first and then input the data in
Matlab. Is there any other solution?
*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



On Tue, 9 Apr 2019 at 00:47, Carlos A. Castro <[email protected]> wrote:

> 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