资源描述：
《电力系统潮流计算完整程序和详细理论说明》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、电力系统潮流计算完整程序及详细理论说明——秦羽风在我刚开始学习潮流程序时，总是找不到一个正确的程序开始模仿学习。后来经过多方努力，终于自己写岀了一个结构清晰、完整的潮流程序。此程序是一个通用的程序，只需要修改输入数据的了函数(PowerFlowsData_K)里面的母线、支路、发电机、负荷，就能算任意一个网络结构的交流系统潮流。很适合初学者学习。为了帮助电力系统的同学一起学习，我将我编写的潮流计算程序分享下来给大家；此程序是在基于牛顿拉夫逊算法的基础上，编写的快速解耦算法。每一个子程序我都有备注说明。如果有不对的地方,希望大家指正！下文中呈现的顺序为：网络结构、
2、子程序、主程序、运算结果、程序设计理论说明。一、网络结构：5节点网络如下图。夫逊法要少得多，而且并不需豪毎次迭代邯对眾可比龙畔環也图4・6五母线测试网络以及棉说廿算结果摘自G.W.StaggandA.H.El-Abiad,ComputerMnlhodsinPowerSywtcui4ma1v«u.I90MM€<•«h.,1二、子程序（共有9个子程序）子程序1:（其他系统，只需要修改Bus、Branch、Generator、Load,这四个矩阵就行了）function[Bus,Branch,Generator,Load]=PowerFlowsData_K%%节点数据
3、%类型：1-平衡节点；2■发电机PV节点；3■负荷PQ节点；4-发电机PQ节点;Bus=[%类型电压相角11.060;210;310;310;310%%线路数据Branch=[%发送接收电阻电感（电导电容）并联00.06;130.080.2400.05;120.020.06230.060.1800.04;240.060.1800.04;250.040.1200.03;340.010.0300.02;450・080.2400・05%%发电机数据Generator=[节点定有功10020.40下限）无功定无功（上限5-5;3-3］；%%负载数据Load=[%节点定有
4、功定无功10.20.1;20.450.15;30.40.05;40.60.1子程序2：%求解网络的导纳矩阵；function[YR,YI]=YBus_K(nbb,ntl,tlsend,tlrecztlresisztlreac,tlcondztlsuscep)YR=zeros(nbbznbb);YI=zeros(nbb,nbb);forkk=l:ntlii=tlsend(kk);jj=tlrec(kk);demon=tlresis(kk)A2+tlreac(kk)A2;YR(ii,ii)=YR(iifii)+tlresis(kk)/demon+0.5*tlcond
5、(kk);YI(iifii)=YI(iifii)-tlreac(kk)/demon+0.5*tlsuscep(kk);YR(ii,jj)=YR(ii,jj)-tlresis(kk)/demon;YI(ii,jj)=YI(ii,jj)+tlreac(kk)/demon;YR(jj,ii)=YR(jj,ii)-tiresis(k.k)/demon;YI(jjzii)=YI(jjzii)+tlreac(kk)/demon;YR(jj,jj)=YR(jj,jj)+tlresis(kk)/demon+0.5*tlcond(kk);YI(jj,jj)=YI(jj,jj)-t
6、lreac(kk)/demon+0.5*tlsuscep(kk);end子程序M%求节点的功率净值function[PNET,QNET]=NetPowers_K(nbb,ngn,nld,genbus,PGEN,QGEN,PLOAD,QLOAD,loadbus)PNET=zeros(1^nbb);QNET=zeros(1znbb);forkk=l:ngnii=genbus(kk);PNET(ii)=PNET(ii)+PGEN(kk);QNET(ii)=QNET(ii)+QGEN(kk);endforkk=l:nldii=loadbus(kk);PNET(ii)=P
7、NET(ii)-PLOAD(kk);QNET(ii)=QNET(ii)-QLOAD(kk);end子程序4：%求网络节点的注入功率function[PCAL,QCAL]=NetPowers_in_K(nbb,VM,VA,YR,YI)PCAL=zeros(1Anbb);QCAL=zeros(1,nbb);forii=l:nbbforjj=l:nbbPCAL(ii)=PCAL(ii)+VM(ii)*VM(jj)*(YR(ii,jj)*cos(VA(ii)-VA(jj))+YI(iizjj)*sin(VA(ii)-VA(jj)))；QCAL(ii)=QCAL(ii)+
8、VM(ii)*VM(jj
