短路电流计算 .zip_电力系统_短路计算_能源行业(电力石油煤炭)

<html xmlns="http://www.w3.org/1999/xhtml"> <head> <meta charset="utf-8"> <meta name="generator" content="pdf2htmlEX"> <meta http-equiv="X-UA-Compatible" content="IE=edge,chrome=1"> <link rel="stylesheet" href="https://static.pudn.com/base/css/base.min.css"> <link rel="stylesheet" href="https://static.pudn.com/base/css/fancy.min.css"> <link rel="stylesheet" href="https://static.pudn.com/prod/directory_preview_static/6270e7407b370112144f78e0/raw.css"> <script src="https://static.pudn.com/base/js/compatibility.min.js"></script> <script src="https://static.pudn.com/base/js/pdf2htmlEX.min.js"></script> <script> try{ pdf2htmlEX.defaultViewer = new pdf2htmlEX.Viewer({}); }catch(e){} </script> <title></title> </head> <body> <div id="sidebar" style="display: none"> <div id="outline"> </div> </div> <div id="pf1" class="pf w0 h0" data-page-no="1"><div class="pc pc1 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/6270e7407b370112144f78e0/bg1.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x1 h3 y2 ff1 fs0 fc0 sc0 ls0 ws0">电力系统短路故障分析的<span class="_ _0"> </span><span class="ff2 sc1">MA<span class="_ _1"></span>TLAB<span class="_ _0"> </span><span class="ff1 sc0">辅助程序设计</span></span></div><div class="t m0 x2 h4 y3 ff1 fs1 fc0 sc1 ls0 ws0">电<span class="_ _2"></span>力<span class="_ _2"></span>系<span class="_ _2"></span>统<span class="_ _2"></span>短<span class="_ _2"></span>路<span class="_ _3"></span>故<span class="_ _2"></span>障<span class="_ _3"></span>可<span class="_ _2"> </span>分<span class="_ _2"></span>为<span class="_ _2"></span>三<span class="_ _2"></span>相<span class="_ _2"></span>对<span class="_ _3"></span>称<span class="_ _2"> </span>短<span class="_ _2"></span>路<span class="_ _2"></span>故<span class="_ _2"></span>障<span class="_ _4"> </span><span class="ff3">(three-pha<span class="_ _5"></span>se<span class="_ _6"></span> <span class="_ _6"></span>balanced</span></div><div class="t m0 x3 h4 y4 ff3 fs1 fc0 sc1 ls0 ws0">faults)<span class="ff1">和<span class="_ _7"></span>不对<span class="_ _7"></span>称短路<span class="_ _7"></span>故障<span class="_ _7"></span></span>(unbalance<span class="_ _7"></span>d faults )<span class="_ _7"></span><span class="ff1">。<span class="_ _7"></span>不对称<span class="_ _7"></span>短路故<span class="_ _7"></span>障又<span class="_ _7"></span>分为单<span class="_ _7"></span>相</span></div><div class="t m0 x3 h4 y5 ff1 fs1 fc0 sc1 ls0 ws0">接<span class="_ _3"></span>地<span class="_ _2"></span>短<span class="_ _3"></span>路<span class="_ _2"></span>故<span class="_ _3"></span>障<span class="_ _8"> </span><span class="ff3">(single<span class="_ _7"></span> <span class="_ _6"></span>line-to<span class="_ _7"></span>-gr<span class="_ _5"></span>ound<span class="_ _7"></span> <span class="_ _6"></span>fault)<span class="_ _2"> </span><span class="ff1">、<span class="_ _2"></span>两<span class="_ _3"></span>相<span class="_ _2"></span>短<span class="_ _2"></span>路<span class="_ _3"></span>故<span class="_ _2"></span>障<span class="_ _3"></span> <span class="_ _8"> </span></span>(line-to-<span class="_ _7"></span>line</span></div><div class="t m0 x3 h4 y6 ff3 fs1 fc0 sc1 ls0 ws0">fault)<span class="_ _7"></span><span class="ff1">以<span class="_ _7"></span>及<span class="_ _6"></span>两<span class="_ _7"></span>相<span class="_ _7"></span>接<span class="_ _7"></span>地<span class="_ _7"></span>短<span class="_ _7"></span>路<span class="_ _7"></span>故<span class="_ _7"></span>障<span class="_ _6"></span></span>(double<span class="_ _7"></span> li<span class="_ _7"></span>ne-to-ground <span class="_ _7"></span>fault)<span class="_ _7"></span><span class="ff1">。<span class="_ _6"></span>根<span class="_ _7"></span>据<span class="_ _7"></span>故<span class="_ _7"></span>障<span class="_ _7"></span>分<span class="_ _7"></span>析</span></div><div class="t m0 x3 h4 y7 ff1 fs1 fc0 sc1 ls0 ws0">结果<span class="_ _7"></span>可以<span class="_ _7"></span>对继<span class="_ _7"></span>电保<span class="_ _7"></span>护装<span class="_ _7"></span>置、<span class="_ _7"></span>自动<span class="_ _7"></span>装置<span class="_ _7"></span>进行<span class="_ _7"></span>整定<span class="_ _7"></span>计算<span class="_ _7"></span>，我<span class="_ _7"></span>们可<span class="_ _7"></span>以建<span class="_ _7"></span>立算<span class="_ _7"></span>法<span class="_ _7"></span>来形<span class="_ _7"></span>成</div><div class="t m0 x3 h4 y8 ff1 fs1 fc0 sc1 ls0 ws0">节点阻抗矩阵<span class="ff3">,</span>利用节点阻抗矩阵来计算短路故障情况下的节点电压和线路电流。</div><div class="t m0 x3 h5 y9 ff1 fs2 fc0 sc1 ls0 ws0">一、三相对称短路故障</div><div class="t m0 x2 h4 ya ff1 fs1 fc0 sc1 ls0 ws0">进<span class="_ _7"></span>行<span class="_ _6"></span>三<span class="_ _7"></span>相<span class="_ _6"></span>短<span class="_ _7"></span>路<span class="_ _7"></span>计<span class="_ _6"></span>算<span class="_ _7"></span>需<span class="_ _6"></span>要<span class="_ _7"></span>两<span class="_ _6"></span>个<span class="_ _7"></span>程<span class="_ _7"></span>序<span class="_ _6"></span>：<span class="_ _3"></span><span class="ff2 fs3">zbuild<span class="_ _7"></span> <span class="_ _7"></span>/<span class="fs1"> <span class="_ _7"></span>zbuildpi<span class="_ _0"> </span></span><span class="ff1 sc0">和<span class="_ _9"> </span></span>symfault<span class="_ _9"> </span><span class="ff1 sc0">程<span class="_ _6"></span>序<span class="_ _7"></span><span class="sc1"> <span class="_ _6"></span> <span class="_ _7"></span>，</span></span></span></div><div class="t m0 x3 h4 yb ff2 fs3 fc0 sc1 ls0 ws0">zbuild<span class="ff1 sc0">、<span class="_ _7"></span></span><span class="fs1">zbuildpi<span class="_ _9"> </span></span><span class="ff1">用来<span class="_ _7"></span>在<span class="_ _4"> </span><span class="ff3">MA<span class="_ _a"></span>TLAB<span class="_ _4"> </span><span class="ff1">中<span class="_ _7"></span>形成节<span class="_ _7"></span>点阻<span class="_ _7"></span>抗矩阵<span class="_ _7"></span>，<span class="_ _7"></span><span class="ff2">symfault<span class="_ _4"> </span></span><span class="fs1">用<span class="_ _7"></span>来计算<span class="_ _7"></span>三相</span></span></span></span></div><div class="t m0 x3 h4 yc ff1 fs1 fc0 sc1 ls0 ws0">对称故障<span class="fs3">。</span></div><div class="t m0 x2 h4 yd ff2 fs1 fc0 sc1 ls0 ws0">Zbus=zbuild(zdata)<span class="ff4"> <span class="ff1">这里的参数<span class="_ _9"> </span></span></span>zdata<span class="_ _4"> </span><span class="ff1">是一个<span class="ff3">(e×4)</span>阶矩阵，<span class="ff3">e<span class="_ _9"> </span></span>是拓</span></div><div class="t m0 x3 h4 ye ff1 fs1 fc0 sc1 ls0 ws0">扑图<span class="_ _7"></span>的总<span class="_ _7"></span>支路<span class="_ _7"></span>数目<span class="_ _7"></span>。第<span class="_ _7"></span>一列<span class="_ _7"></span>和第<span class="_ _7"></span>二列<span class="_ _7"></span>为元<span class="_ _7"></span>素两<span class="_ _7"></span>端的<span class="_ _7"></span>节点<span class="_ _7"></span>编号<span class="_ _7"></span>，第<span class="_ _7"></span>三列<span class="_ _7"></span>和<span class="_ _7"></span>第四<span class="_ _7"></span>列</div><div class="t m0 x3 h4 yf ff1 fs1 fc0 sc1 ls0 ws0">分别<span class="_ _7"></span>是线<span class="_ _7"></span>路的<span class="_ _7"></span>电阻<span class="_ _7"></span>、电<span class="_ _7"></span>抗的<span class="_ _7"></span>标幺<span class="_ _7"></span>值。<span class="_ _7"></span>连接<span class="_ _7"></span>在<span class="_ _0"> </span><span class="ff3">0<span class="_ _4"> </span></span>节<span class="_ _7"></span>点和<span class="_ _7"></span>发电<span class="_ _7"></span>机节<span class="_ _7"></span>点之<span class="_ _7"></span>间的<span class="_ _7"></span>发电<span class="_ _7"></span>机</div><div class="t m0 x3 h4 y10 ff1 fs1 fc0 sc1 ls0 ws0">阻抗<span class="_ _7"></span>可能<span class="_ _7"></span>是次<span class="_ _7"></span>暂态<span class="_ _7"></span>电抗<span class="_ _7"></span>、暂<span class="_ _7"></span>态电<span class="_ _7"></span>抗或<span class="_ _7"></span>同步<span class="_ _7"></span>电抗<span class="_ _7"></span>，而<span class="_ _7"></span>且这<span class="_ _7"></span>个矩<span class="_ _7"></span>阵中<span class="_ _7"></span>还包<span class="_ _7"></span>含<span class="_ _7"></span>并联<span class="_ _7"></span>电</div><div class="t m0 x3 h4 y11 ff1 fs1 fc0 sc1 ls0 ws0">抗器和负荷阻抗。</div><div class="t m0 x2 h4 y12 ff2 fs1 fc0 sc1 ls0 ws0">Zbus=zbuildpi(linedata,gendata,yload<span class="_ _7"></span><span class="ff4">)<span class="ff5"> <span class="_ _7"></span><span class="ff1">这<span class="_ _7"></span>个函<span class="_ _7"></span>数与<span class="_ _7"></span>潮流<span class="_ _7"></span>计算<span class="_ _7"></span>程</span></span></span></div><div class="t m0 x3 h4 y13 ff1 fs1 fc0 sc1 ls0 ws0">序是<span class="_ _7"></span>相容<span class="_ _7"></span>的，<span class="_ _7"></span>第一<span class="_ _7"></span>个参<span class="_ _7"></span>数<span class="_ _0"> </span><span class="ff2">linedat<span class="ff4">a<span class="_ _9"> </span></span></span>与潮<span class="_ _7"></span>流计<span class="_ _7"></span>算程<span class="_ _7"></span>序中<span class="_ _7"></span>的文<span class="_ _7"></span>件是<span class="_ _7"></span>一致<span class="_ _7"></span>的。<span class="_ _7"></span>第一</div><div class="t m0 x3 h4 y14 ff1 fs1 fc0 sc1 ls0 ws0">列和<span class="_ _7"></span>第二<span class="_ _7"></span>列为<span class="_ _7"></span>节点<span class="_ _7"></span>编号<span class="_ _7"></span>；第<span class="_ _7"></span>三列<span class="_ _7"></span>到第<span class="_ _7"></span>五列<span class="_ _7"></span>分别<span class="_ _7"></span>是线<span class="_ _7"></span>路的<span class="_ _7"></span>电阻<span class="_ _7"></span>、电<span class="_ _7"></span>抗以<span class="_ _7"></span>及<span class="_"> </span><span class="ff3">1/2<span class="_ _9"> </span></span>线</div><div class="t m0 x3 h4 y15 ff1 fs1 fc0 sc1 ls0 ws0">路电<span class="_ _7"></span>纳值<span class="_ _7"></span>，这<span class="_ _7"></span>三项<span class="_ _7"></span>都为<span class="_ _7"></span>在统<span class="_ _7"></span>一基<span class="_ _7"></span>准容<span class="_ _7"></span>量下<span class="_ _7"></span>的标<span class="_ _7"></span>幺值<span class="_ _7"></span>；最<span class="_ _7"></span>后一<span class="_ _7"></span>列是<span class="_ _7"></span>变压<span class="_ _7"></span>器<span class="_ _7"></span>分接<span class="_ _7"></span>头</div><div class="t m0 x3 h4 y16 ff1 fs1 fc0 sc1 ls0 ws0">位<span class="_ _3"></span>置<span class="_ _3"></span>，<span class="_ _3"></span>对<span class="_ _6"></span>线<span class="_ _3"></span>路<span class="_ _3"></span>来<span class="_ _3"></span>说<span class="_ _3"></span>，<span class="_ _3"></span>必<span class="_ _3"></span>须<span class="_ _3"></span>输<span class="_ _3"></span>入<span class="_ _b"> </span><span class="ff3">1<span class="_ _3"></span></span>；<span class="_ _6"></span>线<span class="_ _2"></span>路<span class="_ _6"></span>无<span class="_ _3"></span>输<span class="_ _3"></span>入<span class="_ _3"></span>顺<span class="_ _3"></span>序<span class="_ _3"></span>。<span class="_ _3"></span>发<span class="_ _3"></span>电<span class="_ _6"></span>机<span class="_ _2"></span>参<span class="_ _6"></span>数<span class="_ _3"></span>不<span class="_ _3"></span>包<span class="_ _3"></span>含<span class="_ _3"></span>在</div></div></div><div class="pi" data-data='{"ctm":[1.611850,0.000000,0.000000,1.611850,0.000000,0.000000]}'></div></div> </body> </html>

<div id="pf2" class="pf w0 h0" data-page-no="2"><div class="pc pc2 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/6270e7407b370112144f78e0/bg2.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x3 h4 y17 ff2 fs1 fc0 sc1 ls0 ws0">Linedata<span class="_ _9"> </span><span class="ff1">参<span class="_ _7"></span>数<span class="_ _7"></span>中<span class="_ _7"></span>，而<span class="_ _7"></span>是<span class="_ _7"></span>包<span class="_ _7"></span>含在<span class="_ _7"></span>第<span class="_ _7"></span>二<span class="_ _7"></span>个参<span class="_ _7"></span>数<span class="_"> </span></span>gendata<span class="_ _9"> </span><span class="ff1">中<span class="_ _7"></span>，<span class="_ _7"></span></span>gendata<span class="_ _9"> </span><span class="ff1">是<span class="_ _7"></span>一<span class="_ _7"></span>个</span></div><div class="t m0 x3 h4 y18 ff3 fs1 fc0 sc1 ls0 ws0">g×4<span class="_ _4"> </span><span class="ff1">阶矩阵<span class="_ _7"></span></span>,g<span class="_ _4"> </span><span class="ff1">是发电机总数。第一列和第二列为<span class="_ _0"> </span></span>0<span class="_ _4"> </span><span class="ff1">节点、发电机节点编号，第</span></div><div class="t m0 x3 h4 y19 ff1 fs1 fc0 sc1 ls0 ws0">三列和第四<span class="_ _7"></span>列为发电机的暂态电阻<span class="_ _7"></span>和暂态电抗。最后一个<span class="_ _7"></span>参数<span class="_ _9"> </span><span class="ff2">yload<span class="_ _9"> </span></span>是可选择</div><div class="t m0 x3 h4 y1a ff1 fs1 fc0 sc1 ls0 ws0">的，这<span class="_ _7"></span>个矩阵有两<span class="_ _7"></span>列，第一列<span class="_ _7"></span>为节点编<span class="_ _7"></span>号，第二列<span class="_ _7"></span>为复数导<span class="_ _7"></span>纳值，<span class="_ _2"></span><span class="ff2">yload<span class="_ _9"> </span></span>可以</div><div class="t m0 x3 h4 y1b ff1 fs1 fc0 sc1 ls0 ws0">由潮流程序<span class="_ _4"> </span><span class="ff2">lfguss<span class="_ _7"></span></span>，<span class="ff2">lfnewton<span class="_ _4"> </span></span>或者<span class="_ _9"> </span><span class="ff2">decouple<span class="_ _4"> </span></span>自动生成。</div><div class="t m0 x2 h4 y1c ff2 fs1 fc0 sc1 ls0 ws0">Zbuild<span class="_ _c"> </span><span class="ff1">和<span class="_ _7"></span> <span class="_ _7"></span></span>zbuildpi<span class="_ _c"> </span><span class="ff1">两<span class="_ _7"></span>个<span class="_ _7"></span>函<span class="_ _7"></span>数<span class="_ _7"></span>可<span class="_ _7"></span>以<span class="_ _7"></span>通<span class="_ _7"></span>过<span class="_ _7"></span>建<span class="_ _7"></span>立<span class="_ _7"></span>算<span class="_ _7"></span>法<span class="_ _7"></span>求<span class="_ _7"></span>出<span class="_ _7"></span>节<span class="_ _7"></span>点<span class="_ _7"></span>阻<span class="_ _7"></span>抗<span class="_ _7"></span>矩<span class="_ _7"></span>阵<span class="_ _7"></span>。<span class="_ _7"></span>首</span></div><div class="t m0 x3 h4 y1d ff1 fs1 fc0 sc1 ls0 ws0">先添加所有与参考节点相连的树支，然后添加其余的树支，最后添加共轭连支。</div><div class="t m0 x2 h4 y1e ff1 fs1 fc0 sc1 ls0 ws0">程序<span class="_ _9"> </span><span class="ff2">symfault(zdata,Zbus,V<span class="ff4">)<span class="_ _7"></span></span></span>用来计算<span class="_ _7"></span>三相对称故<span class="_ _7"></span>障，程序要<span class="_ _7"></span>求输入</div><div class="t m0 x3 h4 y1f ff2 fs1 fc0 sc1 ls0 ws0">zdata<span class="_ _9"> </span><span class="ff1">和<span class="_ _9"> </span></span>Zbu<span class="_ _7"></span>s<span class="_ _9"> </span><span class="ff1">两个<span class="_ _7"></span>矩<span class="_ _7"></span>阵<span class="_ _7"></span>，<span class="_ _7"></span>第三<span class="_ _7"></span>个<span class="_ _7"></span>参<span class="_ _7"></span>数<span class="_ _c"> </span></span>V<span class="_ _c"> </span><span class="ff1">是<span class="_ _7"></span>可选<span class="_ _7"></span>的<span class="_ _7"></span>。<span class="_ _7"></span>如果<span class="_ _0"> </span></span>V<span class="_ _9"> </span><span class="ff1">不<span class="_ _7"></span>存<span class="_ _7"></span>在，<span class="_ _7"></span>程<span class="_ _7"></span>序<span class="_ _7"></span>将</span></div><div class="t m0 x3 h4 ya ff1 fs1 fc0 sc1 ls0 ws0">默认<span class="_ _7"></span>故障前<span class="_ _7"></span>所有<span class="_ _7"></span>的节点电<span class="_ _7"></span>压标<span class="_ _7"></span>幺值为<span class="_ _0"> </span><span class="ff3">1.0<span class="_ _7"></span></span>，如果<span class="_ _7"></span>变量<span class="_ _9"> </span><span class="ff2">V<span class="_ _9"> </span></span>存<span class="_ _7"></span>在，那<span class="_ _7"></span>么<span class="_ _9"> </span><span class="ff2">V<span class="_ _9"> </span></span>包括<span class="_ _7"></span>节点</div><div class="t m0 x3 h4 yb ff1 fs1 fc0 sc1 ls0 ws0">编<span class="_ _8"> </span>号<span class="_ _8"> </span>和<span class="_ _8"> </span>复<span class="_ _8"> </span>数<span class="_ _8"> </span>电<span class="_ _8"> </span>压<span class="_ _8"> </span>值<span class="_ _2"> </span>。<span class="_ _8"> </span>电<span class="_ _8"> </span>压<span class="_ _8"> </span>向<span class="_ _8"> </span>量<span class="_ _d"> </span><span class="ff2">V<span class="_"> </span></span>也<span class="_ _8"> </span>可<span class="_ _8"> </span>以<span class="_ _8"> </span>由<span class="_ _8"> </span>潮<span class="_ _8"> </span>流<span class="_ _8"> </span>计<span class="_ _8"> </span>算<span class="_ _8"> </span>程<span class="_ _8"> </span>序<span class="_ _2"> </span>自<span class="_ _8"> </span>动<span class="_ _8"> </span>生<span class="_ _8"> </span>成<span class="_ _8"> </span>。<span class="_ _8"> </span>当</div><div class="t m0 x3 h4 yc ff2 fs1 fc0 sc1 ls0 ws0">symfault<span class="_ _9"> </span><span class="ff1">程序<span class="_ _7"></span>运行时<span class="_ _7"></span>，用户<span class="_ _7"></span>要输<span class="_ _7"></span>入故障<span class="_ _7"></span>节点编<span class="_ _7"></span>号和<span class="_ _7"></span>故障阻抗<span class="_ _7"></span>，运<span class="_ _7"></span>行可得<span class="_ _7"></span>到总</span></div><div class="t m0 x3 h4 yd ff1 fs1 fc0 sc1 ls0 ws0">的故障电流，节点电压幅值以及故障情况下的线路电流。</div><div class="t m0 x2 h4 ye ff1 fs1 fc0 sc1 ls0 ws0">在三相<span class="_ _7"></span>短路计算中<span class="_ _7"></span>，<span class="ff2 fs3">zbuild<span class="_ _7"></span> <span class="ff1 sc0">和<span class="_ _4"> </span></span>symfault<span class="_ _4"> </span><span class="ff1 sc0">程<span class="_ _7"></span>序，</span><span class="fs1">zbuildpi<span class="_ _9"> </span></span><span class="ff1 sc0">和<span class="_ _4"> </span></span>symfault<span class="_ _4"> </span><span class="ff1 sc0">程序</span></span></div><div class="t m0 x3 h4 yf ff1 fs3 fc0 sc1 ls0 ws0">都可以进行计算，<span class="fs1">下面是三相短路计算使用的程序代码：</span></div><div class="t m0 x3 h4 y10 ff1 fs1 fc0 sc0 ls0 ws0">（<span class="ff2 sc1">1</span>）<span class="ff2 sc1">Zbuild.m<span class="_ _9"> </span></span>程序代码：</div><div class="t m0 x3 h6 y20 ff6 fs4 fc1 sc1 ls0 ws0">function<span class="fc0"> [Zbus] = zbuild(linedata)</span></div><div class="t m0 x3 h6 y21 ff6 fs4 fc0 sc1 ls0 ws0">nl = linedata(:,1); nr = linedata(:,2); R = linedata(:,3);</div><div class="t m0 x3 h6 y22 ff6 fs4 fc0 sc1 ls0 ws0">X = linedata(:,4);</div><div class="t m0 x3 h6 y23 ff6 fs4 fc0 sc1 ls0 ws0">nbr=length(linedata(:,1)); nbus = max(max(nl), max(nr));</div><div class="t m0 x3 h6 y24 ff6 fs4 fc1 sc1 ls0 ws0">for<span class="fc0"> k=1:nbr</span></div><div class="t m0 x3 h6 y25 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> R(k) == inf | X(k) ==inf</div><div class="t m0 x3 h7 y26 ff6 fs4 fc0 sc1 ls0 ws0"> R(k) = 999999999; X(k) = 999999999;<span class="fc2">%<span class="ff1">无穷</span></span></div><div class="t m0 x3 h6 y27 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">else</span>, <span class="fc1">end</span></div><div class="t m0 x3 h6 y28 ff6 fs4 fc1 sc1 ls0 ws0">end</div><div class="t m0 x3 h6 y29 ff6 fs4 fc1 sc1 ls0 ws0"> </div><div class="t m0 x3 h6 y2a ff6 fs4 fc0 sc1 ls0 ws0">ZB = R + j*X;</div></div></div><div class="pi" data-data='{"ctm":[1.611850,0.000000,0.000000,1.611850,0.000000,0.000000]}'></div></div>

<div id="pf3" class="pf w0 h0" data-page-no="3"><div class="pc pc3 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/6270e7407b370112144f78e0/bg3.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x3 h6 y2b ff6 fs4 fc0 sc1 ls0 ws0">Zbus = zeros(nbus, nbus);</div><div class="t m0 x3 h6 y2c ff6 fs4 fc0 sc1 ls0 ws0">tree=0; </div><div class="t m0 x3 h7 y2d ff6 fs4 fc2 sc1 ls0 ws0">%<span class="ff1">从参考总线<span class="_ _8"> </span></span>0<span class="_ _4"> </span><span class="ff1">上添加一个分支</span></div><div class="t m0 x3 h6 y2e ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">for</span> I = 1:nbr</div><div class="t m0 x3 h6 y2f ff6 fs4 fc0 sc1 ls0 ws0"> ntree(I) = 1;</div><div class="t m0 x3 h6 y30 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> nl(I) == 0 | nr(I) == 0</div><div class="t m0 x3 h6 y31 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> nl(I) == 0 n = nr(I);</div><div class="t m0 x3 h6 y32 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">elseif</span> nr(I) == 0 n = nl(I);</div><div class="t m0 x3 h6 y33 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">end</span></div><div class="t m0 x3 h6 y34 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> abs(Zbus(n, n)) == 0 Zbus(n,n) = </div><div class="t m0 x3 h6 y35 ff6 fs4 fc0 sc1 ls0 ws0">ZB(I);tree=tree+1; <span class="fc2">%%new</span></div><div class="t m0 x3 h6 y36 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">else</span> Zbus(n,n) = Zbus(n,n)*ZB(I)/(Zbus(n,n) + ZB(I));</div><div class="t m0 x3 h6 y37 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">end</span></div><div class="t m0 x3 h6 y38 ff6 fs4 fc0 sc1 ls0 ws0"> ntree(I) = 2;</div><div class="t m0 x3 h6 y39 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">else</span>,<span class="fc1">end</span></div><div class="t m0 x3 h6 y3a ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">end</span></div><div class="t m0 x3 h7 y3b ff6 fs4 fc2 sc1 ls0 ws0">% <span class="ff1">添加一个新总线分支到现有总线上</span></div><div class="t m0 x3 h6 y3c ff6 fs4 fc1 sc1 ls0 ws0">while<span class="fc0"> tree < nbus </span></div><div class="t m0 x3 h6 y3d ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">for</span> n = 1:nbus</div><div class="t m0 x3 h6 y3e ff6 fs4 fc0 sc1 ls0 ws0"> nadd = 1;</div><div class="t m0 x3 h6 y3f ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> abs(Zbus(n,n)) == 0</div><div class="t m0 x3 h6 y40 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">for</span> I = 1:nbr</div><div class="t m0 x3 h6 y41 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> nadd == 1;</div><div class="t m0 x3 h6 y42 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> nl(I) == n | nr(I) == n</div><div class="t m0 x3 h6 y43 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> nl(I) == n k = nr(I);</div><div class="t m0 x3 h6 y44 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">elseif</span> nr(I) == n k = nl(I);</div><div class="t m0 x3 h6 y45 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">end</span></div><div class="t m0 x3 h6 y46 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> abs(Zbus(k,k)) ~= 0</div><div class="t m0 x3 h6 y47 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">for</span> m = 1:nbus</div><div class="t m0 x3 h6 y48 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> m ~= n</div><div class="t m0 x3 h6 y49 ff6 fs4 fc0 sc1 ls0 ws0"> Zbus(m,n) = Zbus(m,k);</div><div class="t m0 x3 h6 y4a ff6 fs4 fc0 sc1 ls0 ws0"> Zbus(n,m) = Zbus(m,k);</div><div class="t m0 x3 h6 y4b ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">else</span>, <span class="fc1">end</span></div><div class="t m0 x3 h6 y20 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">end</span></div><div class="t m0 x3 h6 y21 ff6 fs4 fc0 sc1 ls0 ws0"> Zbus(n,n) = Zbus(k,k) + ZB(I); tree=tree+1;</div><div class="t m0 x3 h6 y22 ff6 fs4 fc0 sc1 ls0 ws0"> nadd = 2; ntree(I) = 2;</div><div class="t m0 x3 h6 y23 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">else</span>, <span class="fc1">end</span></div><div class="t m0 x3 h6 y24 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">else</span>, <span class="fc1">end</span></div><div class="t m0 x3 h6 y25 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">else</span>, <span class="fc1">end</span></div><div class="t m0 x3 h6 y4c ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">end</span></div><div class="t m0 x3 h6 y27 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">else</span>, <span class="fc1">end</span></div><div class="t m0 x3 h6 y28 ff6 fs4 fc1 sc1 ls0 ws0">end</div><div class="t m0 x3 h6 y29 ff6 fs4 fc1 sc1 ls0 ws0"> </div><div class="t m0 x3 h6 y2a ff6 fs4 fc1 sc1 ls0 ws0">end<span class="fc0"> </span></div></div></div><div class="pi" data-data='{"ctm":[1.611850,0.000000,0.000000,1.611850,0.000000,0.000000]}'></div></div>

<div id="pf4" class="pf w0 h0" data-page-no="4"><div class="pc pc4 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/6270e7407b370112144f78e0/bg4.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x3 h6 y4d ff6 fs4 fc0 sc1 ls0 ws0"> </div><div class="t m0 x3 h7 y4e ff6 fs4 fc2 sc1 ls0 ws0">%<span class="ff1">增加两个原有总线间的支路阻抗</span></div><div class="t m0 x3 h6 y4f ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">for</span> n = 1:nbus</div><div class="t m0 x3 h6 y2e ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">for</span> I = 1:nbr</div><div class="t m0 x3 h6 y2f ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> ntree(I) == 1</div><div class="t m0 x3 h6 y30 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> nl(I) == n | nr(I) == n</div><div class="t m0 x3 h6 y31 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> nl(I) == n k = nr(I);</div><div class="t m0 x3 h6 y32 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">elseif</span> nr(I) == n k = nl(I);</div><div class="t m0 x3 h6 y33 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">end</span></div><div class="t m0 x3 h6 y34 ff6 fs4 fc0 sc1 ls0 ws0"> DM = Zbus(n,n) + Zbus(k,k) + ZB(I) - 2*Zbus(n,k);</div><div class="t m0 x3 h6 y35 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">for</span> jj = 1:nbus</div><div class="t m0 x3 h6 y36 ff6 fs4 fc0 sc1 ls0 ws0"> AP = Zbus(jj,n) - Zbus(jj,k);</div><div class="t m0 x3 h6 y37 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">for</span> kk = 1:nbus</div><div class="t m0 x3 h6 y38 ff6 fs4 fc0 sc1 ls0 ws0"> AT = Zbus(n,kk) - Zbus(k, kk);</div><div class="t m0 x3 h6 y39 ff6 fs4 fc0 sc1 ls0 ws0"> DELZ(jj,kk) = AP*AT/DM;</div><div class="t m0 x3 h6 y3a ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">end</span></div><div class="t m0 x3 h6 y50 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">end</span></div><div class="t m0 x3 h6 y3c ff6 fs4 fc0 sc1 ls0 ws0"> Zbus = Zbus - DELZ;</div><div class="t m0 x3 h6 y3d ff6 fs4 fc0 sc1 ls0 ws0"> ntree(I) = 2;</div><div class="t m0 x3 h6 y3e ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">else</span>,<span class="fc1">end</span></div><div class="t m0 x3 h6 y3f ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">else</span>,<span class="fc1">end</span></div><div class="t m0 x3 h6 y40 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">end</span></div><div class="t m0 x3 h6 y41 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">end</span></div><div class="t m0 x3 h4 y51 ff1 fs1 fc0 sc0 ls0 ws0">（<span class="ff2 sc1">2</span>）<span class="ff2 sc1">Zbuildpi.m<span class="_ _9"> </span></span>程序代码：</div><div class="t m0 x3 h6 y46 ff6 fs4 fc2 sc1 ls0 ws0">% This program forms the complex bus impedance matrix by the method</div><div class="t m0 x3 h6 y47 ff6 fs4 fc2 sc1 ls0 ws0">% of building algorithm. Bus zero is taken as reference.</div><div class="t m0 x3 h6 y48 ff6 fs4 fc2 sc1 ls0 ws0">% This program is compatible with power flow data.</div><div class="t m0 x3 h6 y52 ff6 fs4 fc2 sc1 ls0 ws0"> </div><div class="t m0 x3 h6 y4a ff6 fs4 fc1 sc1 ls0 ws0">function<span class="fc0"> [Zbus, linedata] = zbuildpi(linedata, gendata, yload)</span></div><div class="t m0 x3 h6 y4b ff6 fs4 fc0 sc1 ls0 ws0">ng=length(gendata(:,1));</div><div class="t m0 x3 h6 y20 ff6 fs4 fc0 sc1 ls0 ws0">nlg=gendata(:,1);</div><div class="t m0 x3 h6 y21 ff6 fs4 fc0 sc1 ls0 ws0">nrg=gendata(:,2);</div><div class="t m0 x3 h6 y53 ff6 fs4 fc0 sc1 ls0 ws0"> </div><div class="t m0 x3 h6 y23 ff6 fs4 fc0 sc1 ls0 ws0">zg= gendata(:,2) + j*gendata(:,3);</div><div class="t m0 x3 h6 y24 ff6 fs4 fc0 sc1 ls0 ws0">nl = linedata(:,1); nr = linedata(:,2); R = linedata(:,3);</div><div class="t m0 x3 h6 y25 ff6 fs4 fc0 sc1 ls0 ws0">X = linedata(:,4);</div><div class="t m0 x3 h6 y4c ff6 fs4 fc0 sc1 ls0 ws0">nbr=length(linedata(:,1)); nbus = max(max(nl), max(nr));</div><div class="t m0 x3 h6 y27 ff6 fs4 fc0 sc1 ls0 ws0">nc = length(linedata(1,:));</div><div class="t m0 x3 h6 y28 ff6 fs4 fc1 sc1 ls0 ws0">for<span class="fc0"> k=1:nbr</span></div><div class="t m0 x3 h6 y54 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> R(k) == inf | X(k) == inf</div><div class="t m0 x3 h6 y2a ff6 fs4 fc0 sc1 ls0 ws0"> R(k) = 99999999; X(k) = 99999999;</div></div></div><div class="pi" data-data='{"ctm":[1.611850,0.000000,0.000000,1.611850,0.000000,0.000000]}'></div></div>

<div id="pf5" class="pf w0 h0" data-page-no="5"><div class="pc pc5 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/6270e7407b370112144f78e0/bg5.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x3 h6 y2b ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">else</span>, <span class="fc1">end</span></div><div class="t m0 x3 h6 y2c ff6 fs4 fc1 sc1 ls0 ws0">end</div><div class="t m0 x3 h6 y2d ff6 fs4 fc1 sc1 ls0 ws0"> </div><div class="t m0 x3 h6 y2e ff6 fs4 fc1 sc1 ls0 ws0">if<span class="fc0"> nc > 4</span></div><div class="t m0 x3 h6 y2f ff6 fs4 fc0 sc1 ls0 ws0">BC = linedata(:,5);</div><div class="t m0 x3 h6 y30 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">for</span> n = 1:nbus</div><div class="t m0 x3 h6 y31 ff6 fs4 fc0 sc1 ls0 ws0"> yc(n) = 0;</div><div class="t m0 x3 h6 y32 ff6 fs4 fc0 sc1 ls0 ws0"> nlc(n) = 0; nrc(n) = n;</div><div class="t m0 x3 h6 y33 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">for</span> k = 1:nbr</div><div class="t m0 x3 h6 y34 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> nl(k) == n | nr(k) == n</div><div class="t m0 x3 h6 y35 ff6 fs4 fc0 sc1 ls0 ws0"> yc(n) = yc(n) + j*BC(k);</div><div class="t m0 x3 h6 y36 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">else</span>, <span class="fc1">end</span></div><div class="t m0 x3 h6 y37 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">end</span></div><div class="t m0 x3 h6 y38 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">end</span></div><div class="t m0 x3 h6 y39 ff6 fs4 fc1 sc1 ls0 ws0">elseif<span class="fc0"> nc==4 yc= zeros(1, nbr); </span>end</div><div class="t m0 x3 h6 y3a ff6 fs4 fc0 sc1 ls0 ws0">nlc=nlc'; nrc=nrc'; yc=yc.';</div><div class="t m0 x3 h6 y50 ff6 fs4 fc0 sc1 ls0 ws0">ZB = R + j*X;</div><div class="t m0 x3 h6 y3c ff6 fs4 fc1 sc1 ls0 ws0">if<span class="fc0"> exist(<span class="fc3">'yload'</span>) == 1</span></div><div class="t m0 x3 h6 y3d ff6 fs4 fc0 sc1 ls0 ws0"> yload = yload.';</div><div class="t m0 x3 h6 y3e ff6 fs4 fc0 sc1 ls0 ws0"> yc =yc + yload;</div><div class="t m0 x3 h6 y3f ff6 fs4 fc1 sc1 ls0 ws0">else<span class="fc0">, </span>end</div><div class="t m0 x3 h6 y40 ff6 fs4 fc0 sc1 ls0 ws0">m = 0;</div><div class="t m0 x3 h6 y41 ff6 fs4 fc1 sc1 ls0 ws0">for<span class="fc0"> n = 1:nbus</span></div><div class="t m0 x3 h6 y42 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> abs(yc(n)) ~=0</div><div class="t m0 x3 h6 y43 ff6 fs4 fc0 sc1 ls0 ws0"> m=m+1;</div><div class="t m0 x3 h6 y44 ff6 fs4 fc0 sc1 ls0 ws0"> nlcc(m) = nlc(n);</div><div class="t m0 x3 h6 y45 ff6 fs4 fc0 sc1 ls0 ws0"> nrcc(m) = nrc(n);</div><div class="t m0 x3 h6 y46 ff6 fs4 fc0 sc1 ls0 ws0"> zc(m) = 1/yc(n);</div><div class="t m0 x3 h6 y47 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">else</span>, <span class="fc1">end</span></div><div class="t m0 x3 h6 y48 ff6 fs4 fc1 sc1 ls0 ws0">end</div><div class="t m0 x3 h6 y49 ff6 fs4 fc0 sc1 ls0 ws0">nlcc=nlcc'; nrcc=nrcc'; zc=zc.';</div><div class="t m0 x3 h6 y4a ff6 fs4 fc0 sc1 ls0 ws0">nl=[nlg; nlcc; nl]; nr = [nrg; nrcc; nr]; ZB = [zg; zc; ZB];</div><div class="t m0 x3 h6 y4b ff6 fs4 fc0 sc1 ls0 ws0">linedata=[nl nr real(ZB) imag(ZB)];</div><div class="t m0 x3 h6 y20 ff6 fs4 fc0 sc1 ls0 ws0">nbr= length(nl);</div><div class="t m0 x3 h6 y21 ff6 fs4 fc0 sc1 ls0 ws0">Zbus = zeros(nbus, nbus);</div><div class="t m0 x3 h6 y22 ff6 fs4 fc0 sc1 ls0 ws0">tree=0; </div><div class="t m0 x3 h7 y55 ff6 fs4 fc2 sc1 ls0 ws0">%%<span class="ff1">从参考总线<span class="_ _8"> </span></span>0<span class="_ _4"> </span><span class="ff1">上添加一个分支</span></div><div class="t m0 x3 h6 y24 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">for</span> I = 1:nbr</div><div class="t m0 x3 h6 y25 ff6 fs4 fc0 sc1 ls0 ws0"> ntree(I) = 1;</div><div class="t m0 x3 h6 y4c ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> nl(I) == 0 | nr(I) == 0</div><div class="t m0 x3 h6 y27 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> nl(I) == 0 n = nr(I);</div><div class="t m0 x3 h6 y28 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">elseif</span> nr(I) == 0 n = nl(I);</div><div class="t m0 x3 h6 y54 ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">end</span></div><div class="t m0 x3 h6 y2a ff6 fs4 fc0 sc1 ls0 ws0"> <span class="fc1">if</span> abs(Zbus(n, n)) == 0 Zbus(n,n) = ZB(I); </div></div></div><div class="pi" data-data='{"ctm":[1.611850,0.000000,0.000000,1.611850,0.000000,0.000000]}'></div></div>