Re: Matrix is singular to working precision.

Sabhan Kanata Mon, 21 Oct 2019 08:36:40 -0700

 Dear Prof Ray..
How do I check with the case_info command, sir? Please the instructions
myfiledata : casejawabali.m


RegardsSabhan

    Pada Senin, 21 Oktober 2019 22.23.34 WIB, Ray Zimmerman <[email protected]> 
menulis:  
 
 Did you check the output of case_info() for your case?

  Ray

> On Oct 21, 2019, at 10:04 AM, Sabhan Kanata <[email protected]> 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: 28632318
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 <[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] <mailto:[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); >>  >> 
>> 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:[email protected]"; 
>> class="">[email protected]</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:31 ): "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:[email protected]"; class="">[email protected]</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:[email protected]"; target="_blank" 
href="mailto:[email protected]";>[email protected]</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