潮流计算matlab-牛顿拉夫逊法法.zip

  • 杰拉德韦
    了解作者
  • matlab
    开发工具
  • 22KB
    文件大小
  • zip
    文件格式
  • 0
    收藏次数
  • 1 积分
    下载积分
  • 4
    下载次数
  • 2018-08-04 11:49
    上传日期
能计算任何节点的潮流分布,具有较强的实用性,程序易懂
潮流计算matlab-牛顿拉夫逊法法.zip
  • 潮流计算matlab-牛顿拉夫逊法法.doc
    98KB
内容介绍
<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/626a16a90990925c04778b0d/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/626a16a90990925c04778b0d/bg1.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x1 h3 y2 ff1 fs0 fc0 sc0 ls0 ws0">&#29992;<span class="_ _0"> </span><span class="ff2 sc1">matlab<span class="_ _0"> </span></span>&#28526;&#27969;&#35745;&#31639;<span class="ff2 sc1">(</span>&#29275;&#39039;&#25289;&#22827;&#36874;&#27861;<span class="ff2 sc1">)</span></div><div class="t m0 x2 h4 y3 ff2 fs1 fc0 sc1 ls0 ws0">%<span class="ff1 sc0">&#20027;&#31243;&#24207;&#65306;</span></div><div class="t m0 x2 h5 y4 ff3 fs1 fc0 sc1 ls0 ws0">[dfile,pathname]=uigetfile('ieee14.m','Select Data File');</div><div class="t m0 x2 h5 y5 ff3 fs1 fc0 sc1 ls0 ws0">if pathname==0</div><div class="t m0 x2 h5 y6 ff3 fs1 fc0 sc1 ls0 ws0"> error('you must select a valid data file')</div><div class="t m0 x2 h5 y7 ff3 fs1 fc0 sc1 ls0 ws0">else</div><div class="t m0 x2 h5 y8 ff3 fs1 fc0 sc1 ls0 ws0"> lfile=length(dfile);</div><div class="t m0 x2 h5 y9 ff3 fs1 fc0 sc1 ls0 ws0"> % strip off .m</div><div class="t m0 x2 h4 ya ff3 fs1 fc0 sc1 ls0 ws0"> eval(dfile(1:lfile-2)); %<span class="ff1">&#25171;&#24320;&#25968;&#25454;&#25991;&#20214;</span></div><div class="t m0 x2 h5 yb ff3 fs1 fc0 sc1 ls0 ws0">end</div><div class="t m0 x2 h5 yc ff3 fs1 fc0 sc1 ls0 ws0">global n;</div><div class="t m0 x2 h5 yd ff3 fs1 fc0 sc1 ls0 ws0">global m;</div><div class="t m0 x2 h4 ye ff3 fs1 fc0 sc1 ls0 ws0">[nb,mb]=size(bus); %<span class="ff1">&#33410;&#28857;&#37325;&#26032;&#32534;&#21495;</span></div><div class="t m0 x2 h5 yf ff3 fs1 fc0 sc1 ls0 ws0">[nl,ml]=size(line);</div><div class="t m0 x2 h4 y10 ff3 fs1 fc0 sc1 ls0 ws0"> nSW=0; % <span class="ff1">&#24179;&#34913;&#33410;&#28857;&#25968;&#30446;</span></div><div class="t m0 x2 h4 y11 ff3 fs1 fc0 sc1 ls0 ws0"> nPV=0; % PV<span class="_ _1"> </span><span class="ff1">&#33410;&#28857;&#25968;&#30446;</span></div><div class="t m0 x2 h4 y12 ff3 fs1 fc0 sc1 ls0 ws0"> nPQ=0; % PQ<span class="_ _1"> </span><span class="ff1">&#33410;&#28857;&#25968;&#30446;</span></div><div class="t m0 x2 h4 y13 ff3 fs1 fc0 sc1 ls0 ws0"> for i=1:nb, % nb<span class="_ _1"> </span><span class="ff1">&#20026;&#24635;&#33410;&#28857;&#25968;</span></div><div class="t m0 x2 h5 y14 ff3 fs1 fc0 sc1 ls0 ws0"> type=bus(i,6);</div><div class="t m0 x2 h5 y15 ff3 fs1 fc0 sc1 ls0 ws0"> if type==3,</div><div class="t m0 x2 h4 y16 ff3 fs1 fc0 sc1 ls0 ws0"> nSW=nSW+1; % <span class="ff1">&#32479;&#35745;&#24179;&#34913;&#33410;&#28857;&#25968;&#30446;</span></div><div class="t m0 x2 h5 y17 ff3 fs1 fc0 sc1 ls0 ws0"> SW(nSW<span class="_ _2"></span>,:)=bus(i,:);</div><div class="t m0 x2 h5 y18 ff3 fs1 fc0 sc1 ls0 ws0"> elseif type==2,</div><div class="t m0 x2 h4 y19 ff3 fs1 fc0 sc1 ls0 ws0"> nPV=nPV+1; % <span class="ff1">&#32479;&#35745;<span class="_ _1"> </span></span>PV<span class="_ _1"> </span><span class="ff1">&#33410;&#28857;&#25968;&#30446;</span></div><div class="t m0 x2 h5 y1a ff3 fs1 fc0 sc1 ls0 ws0"> PV(nPV<span class="_ _3"></span>,:)=bus(i,:);</div><div class="t m0 x2 h5 y1b ff3 fs1 fc0 sc1 ls0 ws0"> else</div><div class="t m0 x2 h4 y1c ff3 fs1 fc0 sc1 ls0 ws0"> nPQ=nPQ+1; % <span class="ff1">&#32479;&#35745;<span class="_ _1"> </span></span>PQ<span class="_ _1"> </span><span class="ff1">&#33410;&#28857;&#25968;&#30446;</span></div><div class="t m0 x2 h5 y1d ff3 fs1 fc0 sc1 ls0 ws0"> PQ(nPQ,:)=bus(i,:);</div><div class="t m0 x2 h5 y1e ff3 fs1 fc0 sc1 ls0 ws0"> end</div><div class="t m0 x2 h5 y1f ff3 fs1 fc0 sc1 ls0 ws0"> end;</div><div class="t m0 x2 h5 y20 ff3 fs1 fc0 sc1 ls0 ws0"> bus=[PQ;PV<span class="_ _2"></span>;SW];</div><div class="t m0 x2 h5 y21 ff3 fs1 fc0 sc1 ls0 ws0"> newbus=[1:nb]';</div><div class="t m0 x2 h5 y22 ff3 fs1 fc0 sc1 ls0 ws0"> f=bus(:,1);</div><div class="t m0 x2 h5 y23 ff3 fs1 fc0 sc1 ls0 ws0"> nodenum=[newbus bus(:,1)]; </div><div class="t m0 x2 h5 y24 ff3 fs1 fc0 sc1 ls0 ws0"> bus(:,1)=newbus;</div><div class="t m0 x2 h5 y25 ff3 fs1 fc0 sc1 ls0 ws0"> for i=1:nl</div><div class="t m0 x2 h5 y26 ff3 fs1 fc0 sc1 ls0 ws0"> for j=1:2</div><div class="t m0 x2 h5 y27 ff3 fs1 fc0 sc1 ls0 ws0"> for k=1:nb</div><div class="t m0 x2 h5 y28 ff3 fs1 fc0 sc1 ls0 ws0"> if line(i,j)==nodenum(k,2)</div><div class="t m0 x2 h5 y29 ff3 fs1 fc0 sc1 ls0 ws0"> line(i,j)=nodenum(k,1);</div><div class="t m0 x2 h5 y2a ff3 fs1 fc0 sc1 ls0 ws0"> break</div><div class="t m0 x2 h5 y2b ff3 fs1 fc0 sc1 ls0 ws0"> end</div><div class="t m0 x2 h5 y2c ff3 fs1 fc0 sc1 ls0 ws0"> end</div><div class="t m0 x2 h5 y2d ff3 fs1 fc0 sc1 ls0 ws0"> end</div><div class="t m0 x2 h5 y2e ff3 fs1 fc0 sc1 ls0 ws0"> end</div><div class="t m0 x2 h4 y2f ff3 fs1 fc0 sc1 ls0 ws0">Y=y(bus,line); %<span class="ff1">&#24418;&#25104;&#33410;&#28857;&#23548;&#32435;&#30697;&#38453;</span></div><div class="t m0 x2 h4 y30 ff3 fs1 fc0 sc1 ls0 ws0">K=0; %<span class="ff1">&#36845;&#20195;&#27425;&#25968;&#21021;&#20540;</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>
评论
    相关推荐
    • 潮流计算.zip
      30节点潮流计算,MATLAB编程,用于潮流分析,可计算节点过电压
    • 潮流计算.rar
      这是对电力系统进行潮流计算的一个例子,算法采用牛顿-拉夫逊法。
    • 潮流计算功率分布.zip
      给定电力系统网络拓扑、元件参数和发电、负荷参量条件下,计算有功功率、无功功率及电压在电力网中的分布。
    • 潮流计算.rar
      潮流计算一直以来都是大家关心的热点问题,本程序通过牛顿拉夫逊法以及PQ分解法进行了潮流计算的分析
    • IEEE14潮流计算.rar
      标准14节点的电力系统测试程序,用于算例验证等相关电力系统分析内容
    • 潮流计算.zip
      仿真和练习使用,用于小规模电力系统潮流计算
    • 潮流计算.zip
      能够实现任意给定参数的交流系统的潮流计算
    • 潮流计算系统仿真.zip
      潮流计算系统仿真 参数设置可改 系统仿真、波形图正常
    • 潮流计算
      IEEE30节点潮流计算,通过matlab实现
    • 潮流计算.rar
      潮流计算C++语言牛顿法 IEEE4、5、14、30、57、118、300节点系统