Dear Prof Ray,
Thank you for responding
My branch data, There are no R and X parameters that are 0. Where do you think
the problem is?
%% branch data
% fbus tbus r x b rateA rateB rateC ratio angle
status
branch = [
1 2 0.000062650 0.0007008 0.0000000 9900 0 0 0
0 1; %Suralaya-Cilegon
1 24 0.000367768 0.0035333 0.0000000 9900 0 0 0
0 1; %Suralaya-Balaraja
1 31 0.000006000 0.0005600 0.0000000 9900 0 0 0
0 1; %Suralaya-LBE
31 24 0.000072160 0.0071762 0.0000000 9900 0 0 0
0 1; %LBE-Balaraja
2 5 0.001313332 0.0146926 0.2109700 9900 0 0 0
0 1; %Cilegon-Cibinong
4 3 0.000151318 0.0016982 0.0000000 9900 0 0 0
0 1; %Gandul-Kembangan
18 4 0.000034700 0.0003334 0.0000000 9900 0 0 0
0 1; %Depok-Gandul
7 5 0.000444188 0.0042675 0.0000000 9900 0 0 0
0 1; %Bekasi-Cibinong
8 5 0.000621160 0.0059678 0.0000000 9900 0 0 0
0 1; %Muaratawar-Cibinong
11 5 0.000411138 0.0045995 0.2652900 9900 0 0 0
0 1; %Saguling-Cibinong
6 7 0.000197365 0.0018962 0.0000000 9900 0 0 0
0 1; %Cawang-Bekasi
8 6 0.000562560 0.0054048 0.0000000 9900 0 0 0
0 1; %Muaratawar-Cawamg
8 9 0.000282206 0.0027112 0.0000000 9900 0 0 0
0 1; %Muaratawar-Cibatu
10 9 0.000273900 0.0026324 0.0000000 9900 0 0 0
0 1; %Cirata-Cibatu
11 10 0.000147473 0.0014168 0.0000000 9900 0 0 0
0 1; %Saguling-Cirata
12 11 0.000195780 0.0021902 0.0000000 9900 0 0 0
0 1; %Bandung Selatan-Saguling
13 12 0.013981960 0.1343320 0.1924070 9900 0 0 0
0 1; %Mandiracan-Bandung Selatan
13 26 0.001012900 0.0097320 0.1393940 9900 0 0 0
0 1; %Mandiracan-Ujung Berung
26 12 0.000385400 0.0037034 0.0000000 9900 0 0 0
0 1; %Ujung Berung-Bandung Selatan
30 13 0.000200000 0.0022500 0.0000000 9900 0 0 0
0 1; %Pemalang-Mandiracan
14 30 0.000320000 0.0035300 0.0000000 9900 0 0 0
0 1; %Ungaran-Pemalang
15 14 0.000067660 0.0075703 0.4366440 9900 0 0 0
0 1; %Tanjung Djati-Ungaran
16 14 0.002979200 0.0286229 0.4099740 9900 0 0 0
0 1; %Surabaya Barat-Ungaran
14 20 0.000903612 0.0086815 0.0000000 9900 0 0 0
0 1; %Ungaran-Pedan
16 17 0.000140000 0.0013399 0.0000000 9900 0 0 0
0 1; %Surabaya Barat-Gresik
23 16 0.000398638 0.0044597 0.0000000 9900 0 0 0
0 1; %Grati-Surabaya Barat
18 5 0.000091200 0.0008765 0.0000000 9900 0 0 0
0 1; %Depok-Cibinong
19 18 0.001403000 0.0156967 0.9015550 9900 0 0 0
0 1; %Tasikmalaya-Depok
29 19 0.000544600 0.0060933 0.3499780 9900 0 0 0
0 1; %Kesugihan-Tasikmalaya
27 29 0.000110000 0.0012668 0.0000000 9900 0 0 0
0 1; %Adipala-Kesugihan
28 27 0.000038100 0.0004387 0.0000000 9900 0 0 0
0 1; %Cilacap-Adipala
20 29 0.000983900 0.0110073 0.6348830 9900 0 0 0
0 1; %Pedan-Kesugihan
21 20 0.002053100 0.0229694 0.3298180 9900 0 0 0
0 1; %Kediri-Pedan
22 21 0.001026500 0.0114847 0.7239300 9900 0 0 0
0 1; %Paiton-Kediri
22 23 0.000044100 0.0049343 0.3110350 9900 0 0 0
0 1; %Paiton-Grati
24 4 0.000297922 0.0028623 0.0000000 9900 0 0 0
0 1; %Balaraja-Gandul
25 14 0.002347961 0.0225581 0.3244670 9900 0 0 0
0 1; %Ngimbang-Ungaran
16 25 0.000596665 0.0057324 0.0000000 9900 0 0 0
0 1; %Surabaya Barat-Ngimbang
];
RegardsSabhan K
Pada Senin, 21 Oktober 2019 20.55.34 WIB, Ray Zimmerman <[email protected]>
menulis:
It is possible that you have a disconnected network with an island that does
not have a reference bus, or maybe you have a zero impedance branch. I would
begin by looking at the output of case_info() for your case to see if you see
anything unusual.
Ray
On Oct 20, 2019, at 11:59 PM, Sabhan Kanata <[email protected]> wrote:
Dear Prof Ray,
I get a warning when running the power flow as below:
What do you think is the solution?
Error using \ Matrix is singular to working precision.
Error in newtonpf (line 100)
dx = -(J \ F);
Error in runpf (line 148)
[V, success, iterations] = newtonpf(Ybus, Sbus, V0, ref, pv, pq,
mpopt);
RegardsSabhan K
