Re: Matrix is singular to working precision.

Ray Zimmerman Mon, 21 Oct 2019 08:24:13 -0700

Did you check the output of case_info() for your case?

   Ray


> On Oct 21, 2019, at 10:04 AM, Sabhan Kanata <sabhankan...@yahoo.co.id> wrote:
> 
> 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?
> 
> "
OutmailID: 124038130, List: 'matpower-l', MemberID: 82861091
SCRIPT: "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 >  ]; >  >  > 
Regards > Sabhan K > Pada Senin, 21 Oktober 2019 20.55.34 WIB, Ray Zimmerman 
<r...@cornell.edu> 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 <sabhankan...@yahoo.co.id 
<mailto:sabhankan...@yahoo.co.id>> 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); >>  >> 
>> Regards >> Sabhan K >    --Apple-Mail=_A33335EE-D265-4CC7-AB1C-A53F457D7D13 
>> Content-Transfer-Encoding: quoted-printable Content-Type: text/html; 
>> charset="us-ascii"  <html><head><meta http-equiv="Content-Type" 
>> content="text/html; charset=us-ascii"></head><body style="word-wrap: 
>> break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" 
>> class="">Did you check the output of <font face="Courier" 
>> class="">case_info()</font> for your case?<div class=""><br 
>> class=""></div><div class="">&nbsp; &nbsp;Ray<br class=""><div><br 
>> class=""><blockquote type="cite" class=""><div class="">On Oct 21, 2019, at 
>> 10:04 AM, Sabhan Kanata &lt;<a href="mailto:sabhankan...@yahoo.co.id"; 
>> class="">sabhankan...@yahoo.co.id</a>&gt; wrote:</div><br 
>> class="Apple-interchange-newline"><div class=""><div class=""><div 
>> class="ydp26f789a1yahoo-style-wrap" style="font-family: Helvetica Neue, 
>> Helvetica, Arial, sans-serif; font-size: 13px;"><div class=""></div>         
>> <div dir="ltr" data-setdir="false" class="">Dear Prof Ray,</div><div 
>> dir="ltr" data-setdir="false" class=""><br class=""></div><div dir="ltr" 
>> data-setdir="false" class=""><span class="">Thank you for 
>> responding</span><br class=""></div><div dir="ltr" data-setdir="false" 
>> class="">My branch data, <span class="">There are no R and X parameters that 
>> are 0.&nbsp; <span class="">Where do you think the problem 
>> is?</span></span><br class=""></div><div dir="ltr" data-setdir="false" 
>> class=""><br class=""></div><div dir="ltr" data-setdir="false" class=""><div 
>> class="">"
TCL MERGE ERROR ( 10/21/2019 11:23:32 ): "extra characters after close-quote 
branch data<br class="">%&nbsp;&nbsp; &nbsp;fbus&nbsp;&nbsp; 
&nbsp;tbus&nbsp;&nbsp; &nbsp;r&nbsp;&nbsp; &nbsp;x&nbsp;&nbsp; 
&nbsp;b&nbsp;&nbsp; &nbsp;rateA&nbsp;&nbsp; &nbsp;rateB&nbsp;&nbsp; 
&nbsp;rateC&nbsp;&nbsp; &nbsp;ratio&nbsp;&nbsp; &nbsp;angle&nbsp;&nbsp; 
&nbsp;status<br class="">branch = [<br class="">&nbsp;&nbsp; 
&nbsp;1&nbsp;&nbsp; &nbsp; 2&nbsp;&nbsp; &nbsp;0.000062650&nbsp;&nbsp; 
&nbsp;0.0007008&nbsp;&nbsp; &nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; &nbsp;1; %Suralaya-Cilegon<br 
class="">&nbsp;&nbsp; &nbsp;1&nbsp;&nbsp; 24&nbsp;&nbsp; 
&nbsp;0.000367768&nbsp;&nbsp; &nbsp;0.0035333&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Suralaya-Balaraja<br class="">&nbsp;&nbsp; &nbsp;1&nbsp;&nbsp; 
31&nbsp;&nbsp; &nbsp;0.000006000&nbsp;&nbsp; &nbsp;0.0005600&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Suralaya-LBE<br class="">&nbsp;&nbsp; &nbsp;31&nbsp;&nbsp; 
&nbsp;24&nbsp;&nbsp; &nbsp;0.000072160&nbsp;&nbsp; &nbsp;0.0071762&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %LBE-Balaraja<br class="">&nbsp;&nbsp; &nbsp;2&nbsp;&nbsp; 
&nbsp;5&nbsp;&nbsp; &nbsp;0.001313332&nbsp;&nbsp; &nbsp;0.0146926&nbsp;&nbsp; 
&nbsp;0.2109700&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Cilegon-Cibinong<br class="">&nbsp;&nbsp; &nbsp;4&nbsp;&nbsp; 
&nbsp;3&nbsp;&nbsp; &nbsp;0.000151318&nbsp;&nbsp; &nbsp;0.0016982&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Gandul-Kembangan<br class="">&nbsp;&nbsp; &nbsp;18&nbsp;&nbsp; 
&nbsp;4&nbsp;&nbsp; &nbsp;0.000034700&nbsp;&nbsp; &nbsp;0.0003334&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Depok-Gandul<br class="">&nbsp;&nbsp; &nbsp;7&nbsp;&nbsp; 
&nbsp;5&nbsp;&nbsp; &nbsp;0.000444188&nbsp;&nbsp; &nbsp;0.0042675&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;1; %Bekasi-Cibinong<br class="">&nbsp;&nbsp; &nbsp;8&nbsp;&nbsp; 
&nbsp;5&nbsp;&nbsp; &nbsp;0.000621160&nbsp;&nbsp; &nbsp;0.0059678&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;1; %Muaratawar-Cibinong<br class="">&nbsp;&nbsp; &nbsp;11&nbsp;&nbsp; 
&nbsp;5&nbsp;&nbsp; &nbsp;0.000411138&nbsp;&nbsp; &nbsp;0.0045995&nbsp;&nbsp; 
&nbsp;0.2652900&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;1; %Saguling-Cibinong<br class="">&nbsp;&nbsp; &nbsp;6&nbsp;&nbsp; 
&nbsp;7&nbsp;&nbsp; &nbsp;0.000197365&nbsp;&nbsp; &nbsp;0.0018962&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Cawang-Bekasi<br class="">&nbsp;&nbsp; &nbsp;8&nbsp;&nbsp; 
&nbsp;6&nbsp;&nbsp; &nbsp;0.000562560&nbsp;&nbsp; &nbsp;0.0054048&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Muaratawar-Cawamg<br class="">&nbsp;&nbsp; &nbsp;8&nbsp;&nbsp; 
&nbsp;9&nbsp;&nbsp; &nbsp;0.000282206&nbsp;&nbsp; &nbsp;0.0027112&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Muaratawar-Cibatu<br class="">&nbsp;&nbsp; &nbsp;10&nbsp;&nbsp; 
&nbsp;9&nbsp;&nbsp; &nbsp;0.000273900 0.0026324&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Cirata-Cibatu<br class="">&nbsp;&nbsp; &nbsp;11&nbsp;&nbsp; 
&nbsp;10&nbsp;&nbsp; &nbsp;0.000147473 0.0014168&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Saguling-Cirata<br class="">&nbsp;&nbsp; &nbsp;12&nbsp;&nbsp; 
&nbsp;11&nbsp;&nbsp; &nbsp;0.000195780&nbsp;&nbsp; &nbsp;0.0021902&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Bandung Selatan-Saguling<br class="">&nbsp;&nbsp; 
&nbsp;13&nbsp;&nbsp; &nbsp;12&nbsp;&nbsp; &nbsp;0.013981960&nbsp;&nbsp; 
&nbsp;0.1343320&nbsp;&nbsp; &nbsp;0.1924070&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; &nbsp;1; %Mandiracan-Bandung Selatan<br 
class="">&nbsp;&nbsp; &nbsp;13&nbsp;&nbsp; &nbsp;26&nbsp;&nbsp; 
&nbsp;0.001012900&nbsp;&nbsp; &nbsp;0.0097320&nbsp;&nbsp; 
&nbsp;0.1393940&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Mandiracan-Ujung Berung<br class="">&nbsp;&nbsp; &nbsp;26&nbsp;&nbsp; 
&nbsp;12&nbsp;&nbsp; &nbsp;0.000385400&nbsp;&nbsp; &nbsp;0.0037034&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Ujung Berung-Bandung Selatan<br class="">&nbsp;&nbsp; 
&nbsp;30&nbsp;&nbsp; &nbsp;13&nbsp;&nbsp; &nbsp;0.000200000&nbsp;&nbsp; 
&nbsp;0.0022500&nbsp;&nbsp; &nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; &nbsp;1; %Pemalang-Mandiracan<br 
class="">&nbsp;&nbsp;&nbsp; 14&nbsp;&nbsp; &nbsp;30&nbsp;&nbsp; 
&nbsp;0.000320000&nbsp;&nbsp; &nbsp;0.0035300&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Ungaran-Pemalang<br class="">&nbsp;&nbsp; &nbsp;15&nbsp;&nbsp; 
&nbsp;14&nbsp;&nbsp; &nbsp;0.000067660&nbsp;&nbsp; &nbsp;0.0075703&nbsp;&nbsp; 
&nbsp;0.4366440&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Tanjung Djati-Ungaran<br class="">&nbsp;&nbsp; &nbsp;16&nbsp;&nbsp; 
&nbsp;14&nbsp;&nbsp; &nbsp;0.002979200&nbsp;&nbsp; &nbsp;0.0286229&nbsp;&nbsp; 
&nbsp;0.4099740&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Surabaya Barat-Ungaran<br class="">&nbsp;&nbsp; &nbsp;14&nbsp;&nbsp; 
&nbsp;20&nbsp;&nbsp; &nbsp;0.000903612&nbsp;&nbsp; &nbsp;0.0086815&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Ungaran-Pedan<br class="">&nbsp;&nbsp; &nbsp;16&nbsp;&nbsp; 
&nbsp;17&nbsp;&nbsp; &nbsp;0.000140000&nbsp;&nbsp; &nbsp;0.0013399&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Surabaya Barat-Gresik<br class="">&nbsp;&nbsp;&nbsp; 23&nbsp;&nbsp; 
&nbsp;16&nbsp;&nbsp; &nbsp;0.000398638&nbsp;&nbsp; &nbsp;0.0044597&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Grati-Surabaya Barat<br class="">&nbsp;&nbsp; &nbsp;18&nbsp;&nbsp; 
&nbsp; 5&nbsp;&nbsp; &nbsp;0.000091200&nbsp;&nbsp; &nbsp;0.0008765&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Depok-Cibinong<br class="">&nbsp;&nbsp; &nbsp;19&nbsp;&nbsp; 
&nbsp;18&nbsp;&nbsp; &nbsp;0.001403000&nbsp;&nbsp; &nbsp;0.0156967&nbsp;&nbsp; 
&nbsp;0.9015550&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Tasikmalaya-Depok<br class="">&nbsp;&nbsp; &nbsp;29&nbsp;&nbsp; 
&nbsp;19&nbsp;&nbsp; &nbsp;0.000544600&nbsp;&nbsp; &nbsp;0.0060933&nbsp;&nbsp; 
&nbsp;0.3499780&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Kesugihan-Tasikmalaya<br class="">&nbsp;&nbsp; &nbsp;27&nbsp;&nbsp; 
&nbsp;29&nbsp;&nbsp; &nbsp;0.000110000&nbsp;&nbsp; &nbsp;0.0012668&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Adipala-Kesugihan<br class="">&nbsp;&nbsp;&nbsp; 28&nbsp;&nbsp; 
&nbsp;27&nbsp;&nbsp; &nbsp;0.000038100&nbsp;&nbsp; &nbsp;0.0004387&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Cilacap-Adipala<br class="">&nbsp;&nbsp; &nbsp;20&nbsp;&nbsp; 
&nbsp;29&nbsp;&nbsp; &nbsp;0.000983900&nbsp;&nbsp; &nbsp;0.0110073&nbsp;&nbsp; 
&nbsp;0.6348830&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Pedan-Kesugihan<br class="">&nbsp;&nbsp; &nbsp;21&nbsp;&nbsp; 
&nbsp;20&nbsp;&nbsp; &nbsp;0.002053100&nbsp;&nbsp; &nbsp;0.0229694&nbsp;&nbsp; 
&nbsp;0.3298180&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Kediri-Pedan<br class="">&nbsp;&nbsp;&nbsp; 22&nbsp;&nbsp; 
&nbsp;21&nbsp;&nbsp; &nbsp;0.001026500&nbsp;&nbsp; &nbsp;0.0114847&nbsp;&nbsp; 
&nbsp;0.7239300&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Paiton-Kediri<br class="">&nbsp;&nbsp; &nbsp;22&nbsp;&nbsp; 
&nbsp;23&nbsp;&nbsp; &nbsp;0.000044100&nbsp;&nbsp; &nbsp;0.0049343&nbsp;&nbsp; 
&nbsp;0.3110350&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Paiton-Grati<br class="">&nbsp;&nbsp; &nbsp;24&nbsp;&nbsp; &nbsp; 
4&nbsp;&nbsp; &nbsp;0.000297922&nbsp;&nbsp; &nbsp;0.0028623&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Balaraja-Gandul<br class="">&nbsp;&nbsp; &nbsp;25&nbsp;&nbsp; 
&nbsp;14&nbsp;&nbsp; &nbsp;0.002347961&nbsp;&nbsp; &nbsp;0.0225581&nbsp;&nbsp; 
&nbsp;0.3244670&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Ngimbang-Ungaran<br class="">&nbsp;&nbsp; &nbsp;16&nbsp;&nbsp; 
&nbsp;25&nbsp;&nbsp; &nbsp;0.000596665&nbsp;&nbsp; &nbsp;0.0057324&nbsp;&nbsp; 
&nbsp;0.0000000&nbsp;&nbsp;&nbsp; 9900&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; 
&nbsp;0&nbsp;&nbsp; &nbsp;0&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 
&nbsp;1; %Surabaya Barat-Ngimbang<br class="">&nbsp;];<br class=""><br 
class=""></div><div class=""><br class=""></div></div><div dir="ltr" 
data-setdir="false" class="">Regards</div><div dir="ltr" data-setdir="false" 
class="">Sabhan K<br class=""></div>
        
        </div><div id="yahoo_quoted_1773503525" class="yahoo_quoted">
            <div style="font-family:'Helvetica Neue', Helvetica, Arial, 
sans-serif;font-size:13px;color:#26282a;" class="">
                
                <div class="">
                    Pada Senin, 21 Oktober 2019 20.55.34 WIB, Ray Zimmerman 
&lt;<a href="mailto:r...@cornell.edu"; class="">r...@cornell.edu</a>&gt; menulis:
                </div>
                <div class=""><br class=""></div>
                <div class=""><br class=""></div>
                <div class=""><div id="yiv6822524658" class=""><div class="">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 <font class="yiv6822524658" 
face="Courier">case_info()</font> for your case to see if you see anything 
unusual.<div class="yiv6822524658"><br class="yiv6822524658" 
clear="none"></div><div class="yiv6822524658">&nbsp; &nbsp; Ray</div><div 
class="yiv6822524658yqt3153578238" id="yiv6822524658yqt88331"><div 
class="yiv6822524658"><br class="yiv6822524658" clear="none"><div class=""><br 
class="yiv6822524658" clear="none"><blockquote class="yiv6822524658" 
type="cite"><div class="yiv6822524658">On Oct 20, 2019, at 11:59 PM, Sabhan 
Kanata &lt;<a rel="nofollow" shape="rect" class="yiv6822524658" 
ymailto="mailto:sabhankan...@yahoo.co.id"; target="_blank" 
href="mailto:sabhankan...@yahoo.co.id";>sabhankan...@yahoo.co.id</a>&gt; 
wrote:</div><br class="yiv6822524658Apple-interchange-newline" 
clear="none"><div class="yiv6822524658"><div class="yiv6822524658"><div 
class="yiv6822524658yahoo-style-wrap" style="font-family:Helvetica Neue, 
Helvetica, Arial, sans-serif;font-size:13px;"><div class="yiv6822524658" 
dir="ltr"><div class="yiv6822524658"><div class="yiv6822524658" dir="ltr"><div 
class="yiv6822524658"><div class="yiv6822524658" dir="ltr">Dear Prof 
Ray,</div><div class="yiv6822524658" dir="ltr"><br class="yiv6822524658" 
clear="none"></div><div class="yiv6822524658">I get a warning when running the 
power flow as below:</div><br class="yiv6822524658" clear="none">What do you 
think is the solution?</div><div class="yiv6822524658"><br 
class="yiv6822524658" clear="none"></div></div><div class="yiv6822524658"><br 
class="yiv6822524658" clear="none"></div><div class="yiv6822524658">Error 
using&nbsp; \ </div>Matrix is singular to working precision.<br 
class="yiv6822524658" clear="none"><br class="yiv6822524658" clear="none">Error 
in newtonpf (line 100)<br class="yiv6822524658" clear="none">&nbsp;&nbsp; dx = 
-(J \ F);<br class="yiv6822524658" clear="none"><br class="yiv6822524658" 
clear="none">Error in runpf (line 148)<br class="yiv6822524658" 
clear="none">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
[V, success, iterations] = newtonpf(Ybus, Sbus, V0, ref, pv, pq, 
mpopt);</div><div class="yiv6822524658"><br class="yiv6822524658" 
clear="none"></div><div class="yiv6822524658" dir="ltr">Regards</div><div 
class="yiv6822524658" dir="ltr">Sabhan K<br class="yiv6822524658" 
clear="none"></div></div></div></div></div></blockquote></div><br 
class="yiv6822524658" clear="none"></div></div></div></div></div>
            </div>
        </div></div></div></blockquote></div><br class=""></div></body></html>

Re: Matrix is singular to working precision.

Reply via email to