高压输电网潮流的计算机算法程序设计课程设计.docx

贵州大学电气工程学院电力系统分析课程设计汇报书题目：高压输电网时尚的计算机算法程序设计毕业设计（论文）原创性申明和使用授权阐明原创性申明本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及获得的成果。尽我所知，除文中尤其加以标注和道谢日勺地方外，不包括其他人或组织已经刊登或公布过时研究成果，也不包括我为获得及其他教育机构的学位或学历而使用过的材料。对本研究提供过协助和做出过奉献的个人或集体，均已在文中作了明确日勺阐明并表达了谢意。作者签名：日期：指导教师签名：B期：使用授权阐明本人完全理解大学有关搜集、保留、使用毕业设计（论文）的规定，即：按照学校规定提交毕业设计（论文）的印刷本和电子版本；学校有权保留毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其他复制手段保留论文；在不以获利为目日勺前提下，学校可以公布论文的部分或所有内容。日期:作者签名:学位论文原创性申明本人郑重申明：所呈交的论文是本人在导师的指导下独立进行研究所获得的研究成果。除了文中尤其加以标注引用的内容外，本论文不包括任何其他个人或集体已经刊登或撰写日勺成果作品。对本文的研究做出重要奉献的个人和集体，均已在文中以明确方式标明。本人完全意识到本申明的法律后果由本人承担。作者签名：日期：年月日学位论文版权使用授权书本学位论文作者完全理解学校有关保留、使用学位论文日勺规定，同意学校保留并向国家有关部门或机构送交论文日勺复印件和电子版，容许论文被查阅和借阅。本人授权大学可以将本学位论文日勺所有或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保留和汇编本学位论文。涉密论文按学校规定处理。作者签名：日期：年月日导师签名:日期：年月指导教师评阅书指导教师评价：一、撰写（设计）过程1、学生在论文（设计）过程中的治学态度、工作精神口优口良口中口及格口不及格2、学生掌握专业知识、技能的扎实程度口优口良口中口及格口不及格3、学生综合运用所学知识和专业技能分析和处理问题的能力口优口良口中口及格口不及格4、研究措施日勺科学性；技术线路的可行性；设计方案的合理性口优口良口中口及格口不及格5、完毕毕业论文（设计）期间的出勤状况口优口良口中口及格口不及格二、论文（设计）质量1、论文（设计）的整体构造与否符合撰写规范？ 优口良口中口及格口不及格2、与否完毕指定的论文（设计）任务（包括装订及附件）？ 优口良口中口及格口不及格三、论文（设计）水平1、论文（设计）日勺理论意义或对处理实际问题的指导意义口优口良口中口及格口不及格2、论文日勺观念与否有新意？设计与否有创意？口优口良口中口及格口不及格3、论文（设计阐明书）所体现的整体水平口优口良口中口及格口不及格提议成绩：口优口良口中口及格口不及格（在所选等级前的口内画“）指导教师：（签名）单位：（盖章）年月曰评阅教师评阅书评阅教师评价：一、论文（设计）质量1、论文（设计）日勺整体构造与否符合撰写规范？ 优口良口中口及格口不及格2、与否完毕指定的论文（设计）任务（包括装订及附件）？优口良口中口及格口不及格二、论文（设计）水平1、论文（设计）的理论意义或对处理实际问题的指导意义口优口良口中口及格口不及格2、论文的观念与否有新意？设计与否有创意？口优口良口中口及格口不及格3、论文（设计阐明书）所体现的整体水平口优口良口中口及格口不及格提议成绩：口优口良口中口及格口不及格（在所选等级前的口内画“J”）评阅教师:教研室（或答辩小组）及教学系意见教研室（或答辩小组）评价：一、答辩过程1、毕业论文（设计）的基本要点和见解的论述状况口优口良中口及格口不及格2、对答辩问题的反应、理解、体现实状况况口优口良中口及格口不及格3、学生答辩过程中的精神状态口优口良中口及格口不及格二、论文（设计）质量1、论文（设计）的整体构造与否符合撰写规范？口优口良中口及格口不及格2、与否完毕指定的论文（设计）任务（包括装订及附件）？口优口良中口及格口不及格三、论文（设计）水平1、论文（设计）的理论意义或对处理实际问题的指导意义口优口良中口及格口不及格2、论文的观念与否有新意？设计与否有创意？口优口良中口及格口不及格3、论文（设计阐明书）所体现的整体水平优口良中及格口不及格评估成绩：优口良中口及格口不及格（在所选等级前的口内画“J”）教研室主任（或答辩小组组长）:（签名）年月日教学系意见：系主任：（签名）年月日目录第一章课程设计概述-4-1.1 设计目的-4-13设计题目51.4 设计内容-5-1.5 设计时间-5-第二章时尚计算设计题目及思绪-6-2.1 时尚计算题目-6-2.2 对题目的分析及设计思绪-6-第三章电力系统时尚计算概述-7-3.1 电力系统概述-7-3.2 时尚计算的意义-7-3.3 牛顿一拉夫逊法时尚计算基本环节-8-第四章程序流程图及程序代码-9-4.1 程序流程图-9-4.2 运用MAT1.AB编程程序代码-10-4.3 例题数据代码-25-第五章运行成果分析-26-5.1 运用MAT1.AB运行成果-26-5.2 运用powerworld软件仿真-38-5.2.1 使用powerworld仿真时尚分布-38-第六章总结-40-参照文献-40-第一章课程设计概述1.1 设计目的1 .掌握电力系统时尚计算的基本原理和电力系统运行方式B¾变化；2 .掌握并能纯熟运用一门计算机语言（MAT1.AB语言或C语言或C+语言）；3 .采用计算机语言对时尚计算进行计算机编程计算。1.2 设计规定1 .程序源代码；2 .选定算例0输入，输出文献；3 .程序阐明；4 .选定算例的程序计算过程；5 .选定算例的手算过程（至少迭代2次）（可选）。1.3 设计题目高压输电网时尚的计算机算法程序设计（PQ分解法、牛顿-拉夫逊法）或中压配电网时尚的计算机算法程序设计（前推后裔法、同伦延拓法等）或电力系统短路故障的计算机算法程序设计（规定不限）1.4 设计内容1 .根据电力系统网络推导电力网络数学模型，写出节点导纳矩阵；2 .赋予各节点电压变量（直角坐标系形式）初值后，求解不平衡量；3 .形成雅可比矩阵；4 .求解修正量后，重新修改初值，从2开始重新循环计算；5 .求解的电压变量到达所规定的精度时，再计算各支路功率分布、功率损耗和平衡节点功率；6 .上机编程调试；连调；7 .计算分析给定系统时尚分析并与手工计算成果作比较分析。8 .准备计算机演示答辩，书写该课程设计阐明书（必须计算机打印）。1.5设计时间春季第17周至第18周第二章时尚计算设计题目及思绪2.1 时尚计算题目题目：设计原始资料网络接线如下图所示：其中已知：各支路阻抗表幺值参数分别为Z12=0.04+j00.25,Z13=0.10+j.35,Z23=O.O8+j.3O,ZT1=j.15,ZT2=j.O3,K=1.05,1.1、1.2电纳为0.5。运用以极坐标表达的高斯一塞德尔法、牛顿一拉夫逊法、P-Q分解法计算该系统的时尚分布。计算精度规定各节点电压修正量不不小于IT，°2.2 对题目的分析及设计思绪此电力系统是一种5节点，3支路的电力网络。然后确定节点类型，由于一般无发电设备的变电所、功率固定的发电厂为PQ节点；有可调无功设备的变电所、有励磁储备的发电厂为PV节点；主调频电厂或出线多的发电厂为平衡节点。通过对此电力系统分析，我们很轻易得以节点1、2、3为PQ节点，节点4为PV节点，节点5为平衡节点。综合比较牛顿拉夫逊法（直角坐标、极坐标）、PQ分解法等多种求解措施的特点，最终确定采用牛顿拉夫逊法（极坐标）。由于此措施所需解的方程组至少。第三章电力系统时尚计算概述3.1 电力系统概述电力工业发展初期，电能是直接在顾客附近的发电站（或称发电厂）中生产的，各发电站孤立运行。伴随工农业生产和都市aJ发展，电能B需要量迅速增长，而热能资源和水能资源丰富的地区又往往远离用电比较集中的都市和工矿区，为了处理这个矛盾，就需要在动力资源丰富B¾地区建立大型发电站，然后将电能远距离输送给电力顾客。同步，为了提高供电B¾可靠性以及资源运用B¾综合经济性，又把许多分散0¾多种形式0发电站，通过送电线路和变电所联络起来。这种由发电机、升压和降压变电所，送电线路以及用电设备有机连接起来的整体，即称为电力系统。3.2 时尚计算的意义（D在电网规划阶段，通过时尚计算，合理规划电源容量及接入点，合理规划网架，选择无功赔偿方案，满足规划水平的大、小方式下时尚互换控制、调峰、调相、调压的规定。（2）在编制年运行方式时，在估计负荷增长及新设备投运基础上，选择经典方式进行时尚计算，发现电网中微弱环节，供调度员平常调度控制参照,并对规划、基建部门提出改善网架构造，加紧基建进度0¾提议。(3)正常检修及特殊运行方式下B¾时尚计算，用于日运行方式B编制，指导发电厂开机方式，有功、无功调整方案及负荷调整方案，满足线路、变压器热稳定规定及电压质量规定。(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。总结为在电力系统运行方式和规划方案的研究中，都需要进行时尚计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同步，为了实时监控电力系统的运行状态，也需要进行大量而迅速的时尚计算。因此，时尚计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线时尚计算；在电力系统运行状态的实时监控中，则采用在线时尚计算。U(O)皿&。)=Pix-22UiUj(Gi.cosSiy÷Bijsin%。)3&牛酶层逊翅睛计算基本坏节。，)牛顿一拉夫逊法时尚计算的基本环节不外乎有如下儿步。1、形成节点导纳矩阵。2、静曲率电压甲慨S值H和求出功率偏差qJ,(w3、根据公式算出雅可比矩阵。4、解如下线性方程组，得到根据如下公式计算修正量，并形成迭代公式:这样反更迭代计算，直至所有节点IAUIV£和A<为止。5、计算线路功率分布，PV节点无功功率和平衡节点注入功率。第四章程序流程图及程序代码4.1 程序流程图根据牛顿一拉夫逊法时尚计算的基本环节可以得到如下程序流程图4.2 运用MAT1.AB编程程序代码铲*主程序开始*clc;clearallformatlong;告*3*globalnSWnPQnPV;,平衡节点，PQ节点，PV节点个数globalnb;,节点个数globalnl;3支路数globalbus;告bus:节点数据globalline;%line:支路数据globalY;告Y:节点导纳矩阵globalnodenum;考节点编号对应表global1PQ；%线路损耗矩阵globalmyf;%输出文献句柄¾*1.*子周月*openfile;考打开busIine数据文献change;务节点重新编号ybus;先建立节点导纳矩阵NR;名牛顿拉夫逊法解非线性方程序组PQ；%计算节点注入功率flow;告计算线路功率及损耗ret;告恢复原节点编号output;%计算成果输出*主程序结束*才丁pbusX*functionopenfileglobalnbnl;globalline;globalbus;globalmyf;(dfiIezpathname=Uigetfile(,x.m,z,SelectDataFile');3打开数据文献ifpathname=0error(,youmustselectavaliddatafile,)elseIfile=Iength(dfile);%stripoff.meval(dfile(1:Ifile-2);先执行数据文献中语句endnl,ml=size(line);告求线路数nlnb,mb=size(bus);官求节点数nbOutfile=,;forI=I:length(dfile)ifdfile(I)break;endOutfile=Strcat(outfile,dfile(I);endOutfile=Strcat(outfile,.txt,);为输出文献名形成myf=foen(outfilez,w,);%打开输出文献莒*给节点编号*functionchangeglobalnb;globalnl;globalnPQ;globalbus;globalline;globalnodenum;nSW=O;nPV=O;nPQ=O;forI=1:nb,%nb为总节点数type=bus(I,6);iftype3,务求平衡节点数nSW=nSW+1;SW(nSW,z)=bus(I,:);elseiftype=2,告求PV节点数nPV=nPV+1;PV(nPVz:)=bus(I,:);else行求PQ节点数nPQ=nPQ+1;PQ(nPQ,:)=bus(I,:);endendbus=PQ;PV;SW;按PQPVSW次序排列H勺节点数据矩阵newbus=1:nb,;nodenum=newbusbus(:,1);当形成节点号对应表bus(:,1)=newbus;forI=l:nl%根据节点号对应表重新排列线路矩阵节点编号forJ=I:2fork=l:nbifline(I,J)=nodenum(kz2)line(I,J)=nodenum(k,1);breakendendendend%*建立节点导纳矩阵*functionY=ybus(bus,line)Y=zeros(nb,nb);fork=l:nl3对导纳矩阵赋初值0globalnl;globalnb;globalbus;globalline;globalY;globalmyf;I=Iine(k.,1);当读入线路参数J=Iine(k,2);Zt=Iine(k,3)+j*line(k,4);ifJ-=O名接地支路分母为零Yt=lZt;彩非接地支路为阻抗的倒数endYm=Iine(k,5)+j*line(k,6);告求YmK=Iine(k,7);落求Kif(K=O)&(J=0)%i般线路：=0;Y(IzI)=Y(I,I)+Yt+Ym;Y(J,J)=Y(JzJ)+Yt+Ym;Y(IrJ)=Y(IzJ)-Yt;Y(JzI)=Y(IzJ);endif(K=O)&(J=O)%对地支路：K=O,J=OtR=X=O;Y(I,I)=Y(IzI)+Ym;endifK>0%变压器线路：Zt和Ym为折算到i侧H勺值，K在j侧Y(I,I)=Y(IzI)+Yt+Ym;Y(JzJ)=Y(JzJ)+YtKK;Y(IzJ)=Y(IzJ)-YtK;Y(JzI)=Y(IzJ);endifK<0%变压器线路：Zt和Ym为折算到K所在侧的值，K在i侧Y(I,I)=Y(IzI)+Yt+Ym;Y(J,J)=Y(J,J)+K*K*Yt;Y(I,J)=Y(IzJ)+K*Yt;Y(JzI)=Y(IzJ);end%*functionNRglobalnb;globalnPQ;globalbus;globalline;globalY;globalmyf;maxl=100;%最大迭代次数*牛顿拉夫逊法解非线性方程序组*%收敛精度epsl=l.0e-10;eps2=l.0e-10;fprintf(myf,节点导纳矩阵Yn,);告输出节点导纳矩阵forI=I:nbforJ=I:nbfprintf(myfz,%10f+j*(%10f)',real(Y(I,J)zimag(Y(I,J);endfprintf(myf,n,);endfori=l:maxlangl(:,1)=bus(1:nb-l,3);u(:,1)=bus(1:nPQ,2);x=angl;u;电压形成的列向量作为XJac=form_jac(bus,Y);del=dPQ(Y,bus);dx=Jacdel;告从bus矩阵中得到PQPV节点的相角、PQ节点的也得到雅可比矩阵告得到功率偏差列向量当求得dxfprintf(myf,第Wd次迭代成果n,i);fprintf(myf,第%d次迭代的雅比矩阵Jn,zi);forI=I:nb+nPQ-1forJ=I:nb+nPQ-1fprintf(myfz,%10f',Jac(I,J);当输由每次迭代肚J雅可比矩阵endfprintf(myf,n');endfprintf(myf,第冬d次迭代的功率偏差dP和dQn,i);forI=I:nb+nPQ-1fprintf(myf,%10e'zdel(I,l);当输出每次迭代的功率偏差fprintf(myf,'n')fprintf(myfz第*d次迭代的节点相角和电压的偏差dxn,i);forI=I:nb÷nPQ-lfprintf(myf,%10e,zdx(1,1);原输出每次迭代U勺节点相角和电压的偏差fprintf(myf,n');endforI=nb:nb+nPQ-1dx(Izl)=dx(Izl)*x(I,l);理%由于求得的是delUU,故还需作对应处end书修正x=x-dx;fprintf(myf,第当d次迭代I向节点相角delta和电压Un,i);angl=x(1:nb-l,1);为将合一起Bx,拆成相角、电压两个列向量u=x(nb:nb+nPQ-1,1);forI=I:nb-1%输出修正的JbUS矩阵中PQ、PV结点的相角bus(1,3)=angl(1,1);fprintf(myf,ang%d%10fn,I,angl(Iz1);endforI=IznPQ%输出修正H勺bus矩阵中PQ结点出)电压bus(I,2)=u(1,1);fprintf(myf,U%d%10fn,zI,u(1,1);endif(max(abs(dx)<epsl)&(max(abs(dPQ(Y,bus)<eps2)当判断与否到达所指定的精度规定breakendendifi=maxl%超过最大迭代次数。出错指示error超过最大迭代次数，不收敛停机！,)；end&*计算节点注入功率*functionPQglobalnb;globalnPQ;globalbus;globalline;globalY;forI=nPQ+l:nb为求PV结点的Q,平衡节点的P,Qifbus(I,6)=33对于平衡节点，求其Psum=0;forJ=I:nbang=bus(I,3)-bus(J,3);A=real(Y(I,J)*cos(ang)+imag(Y(IzJ)*sin(ang);sum=sum+bus(1,2)*bus(J,2)*（八）;endbus(I,4)=sum;3求取口勺平衡节点的P存回bus矩阵endsum=0;forJ=I:nbang=bus(I,3)-bus(J,3);B=real(Y(IzJ)*sin(ang)-imag(Y(I,J)*cos(ang);sum=sum+bus(I,2)*bus(Jz2)*B;endbus(I,5)=sum;务求取卧JPV、平衡节点型JQ存回bus矩阵end告*计算线路功率及损耗*functionflowglobalnl;global1PQ；globalbus;globalline;fork=l:nlI=Iine(k,1);告读入线路参数J=Iine(k,2);lPQ(kzl)=I;当用IPQ矩阵第一、二列保留线路两端节点编号lPQ(k,2)=J;ifJ-=OZt=Iine(k,3)+j*line(k,4);Yt=lZt;endYm=Iine(k,5)+j*line(k,6);K=Iine(kz7);Ui=bus(I,2)*(cos(bus(I,3)+j*sin(bus(Iz3);电用复数表达UiifJ=0Uj=bus(J,2)*(cos(bus(Jz3)+j*sin(bus(Jz3);%非接地节点,则表达Ujif(K=O)&(J-=O)当一般线路非接地支路;Iij=Ui*(Yt+Ym)-Uj*Yt;Iji=Uj*(Yt+Ym)-Ui*Yt;endif(K=O)&(J=O)与接地支路Iij=Ui*Ym;Iji=O;endifK>0%变压器线路：Zt和Ym为折算到i侧时值，K在j侧Iij=(Ui-Uj)*YtK+Ui*(Ym+Yt*(K-I)/K);Iji=(Uj-Ui)*YtK+Uj*Yt*(I-K)/K2;endifK<0舍变压器线路：Zt和Ym为折算到长侧的值，K在i侧K=-K;Iij=(Ui-Uj)*Yt*K+Ui*(Ym+Yt*(I-K);Iji=(Uj-Ui)*Yt*K+Uj*Yt*K*(K-I);endSij=Ui*conj(Iij);省求取节点注入功率Sji=Uj*conj(Iji);delS=Sij+Sji;当求线路损耗lPQ(k,3:5)=SijSjidels;告保留到IPQ矩阵end节点functionretglobalnl;globalnb;globalbus;globalline;globalY;globalnodenum;global1PQ;forI=I:nb钻US矩阵中的节点号恢复本来的编号forJ=I:nb%查找其原编号ifnodenum(Jz2)=Ibreak;endendtem(I,:)=bus(nodenum(J,1),:);当bus矩阵中R节点号恢复本来的编号tem(Iz1)=1;endbus=tem;forI=I:nl告Iine矩阵中的1节点号恢复本来的）编号forJ=I:2fork=l:nbifIPQ(I,J)=nodenum(k,1)IPQ(I,J)=nodenum(k,2);breakendendendend*将计算成果输出*functionoutputglobalnb;globalnl;globalbus;global1PQ;globalmyf;fprintf(myf,n牛顿-拉夫逊法时尚计算成果n');fprintf(myfz，节点计算成果：r);偷出节点计算成果到文献fprintf(myf,JA1Ai工点电压节点相角(角度)节点注入功率r);forI=I:nb,fprintf(myf,'%2d%10f%10f%10f+j%10fbus(1,1),bus(Iz2),bus(I,3)*180pi,bus(I,4)fbus(工，5);endfprintf（myf,，n线路计算成果：r）;吉输出线路计算成果到文献fprintf(myf,点工节点J线路功率SQ,J)线路功率S(J,工)线路损耗dS(工，J)r);forI=I:nl,fprintf(myf,%2d%2d%10f+j%10f%10f+j%10f%10f+j%10fn,lPQ(Izl)zlPQ(Iz2)zreal(1PQ(I,3)zimag(1PQ(1,3)zreal(lPQ(Iz4),imag(lPQ(I,4)zreal(lPQ(I,5),imag(1PQ(Iz5);endfclose(myf);4.3例题数据代码%(bus#)(volt)(ang)(P)(q)(bustype)bus：=11.000.00-1.60-0.801；21.000.00-2.00-1.001;31.000.00-3.70-1.301;41.050.005.000.002;51.050.000.000.003;%b#lb#2(R)(X)(G)(B)(K)120.040.250.00.250;130.100.350.00.00;230.080.300.00.250;530.000.030.00.01.05;420.000.0150.00.01.05;line第五章运行成果分析5.1运用MAT1.AB运行成果节点导纳矩阵Y-1.378742+j*(-6.291665)-0.624025+j*(3.900156)-0.754717+j*(2.641509)0.000000+j*(0.000000)0.000000+j*(0.000000)-0.624025+j*(3.900156)1.453900+j*(-66.980821)-0.829876+j*(3.112033)0.000000+j*(63.492063)0.000000+j*(0.000000)-0.754717+j*(2.641509)-0.829876+j*(3.112033)1.584592+j*(-35.737859)0.000000+j*(0.000000)0.000000+j*(31.746032)0.000000+j*(0.000000)0.000000+j*(63.492063)0.000000+j*(0.000000)0.000000+j*(-66.666667)0.000000+j*(0.000000)0.000000+j*(O.OOOOOO)0.000000+j*(0.000000)0.000000+j*(31.746032)0.000000+j*(0.000000)O.OOOOOO+j*(-33.333333)第1次迭代成果第1次迭代的雅比矩阵J-6.5416653.9001562.641509-0.000000-1.3787420.6240250.7547173.900156-73.6788563.11203366.6666670.624025-1.4539000.8298762.6415093.112033-39.086876-0.0000000.7547170.829876-1.584592-0.00000066.666667-0.000000-66.666667-0.000000-0.000000-0.0000001.378742-0.624025-0.7547170.000000-6.0416653.9001562.641509-0.6240251.453900-0.8298760.0000003.900156-60.2827863.112033-0.754717-0.8298761.5845920.0000002.6415093.112033-32.388841第1次迭代的功率偏差dP和dQ-1.600000e÷000-2.000000e+000-3.700000e+0005.000000e+000-5.500000e-0015.698O35e+OOO2.049017e+OOO第1次迭代的节点相角和电压的偏差dx3.34820le-002-3.607046e-0016.900000e-002-4357046e-0013.356932e-002-1.053820e-001-5.88131le-002第1次迭代的节点相角delta和电压Uangl-0.033482ang20.360705ang3-0.069000ang40.435705Ul0.966431U21.105382U31.058813第2次迭代成果-第2次迭代的雅比矩阵J2.215658-6.3195833.5908882.7286950.0000000.3160084.102936-81.3037223.71581473.484972-0.9846460.349720-0.6344112.6738472.906538-40.7901720.0000000.8677762.4003722.148559-0.00000073.484972-0.000000-73.484972-0.000000-5.521730-0.0000002.891465-2.215658-0.6758080.000000-5.4331003.5908882.7286950.9846463.9026730.634411-5.5217304.102936-82.3799083.715814-0.867776-2.4003725.7014860.0000002.6738472.906538-39.340197第2次迭代的功率偏差dP和dQ3.736505e-0031.261967e-0012.25O225e-OOl-5.217300e-001-3.567586e-001-1.538O93e+OOO-5.750124e-001-第2次迭代时节点相角和电压的偏差dx4.313772e-0025.062770e-0035.188733e-OO29.447229e-0022.354765e-OO21.971867e-002第2次迭代时节点相角delta和电压Uangl-0.076620ang20.314148ang3-0.074063ang40.383817Ul0.875130U21.079353U31.037935第3次迭代成果第3次迭代的雅比矩阵J-5.5793553.1817602.3975950.0000000.5280651.9482270.6916633.630791-78.9919833.57889871.782294-0.8582190.304143-0.4591922.4011012.875049-39.7791270.0000000.6793932.1802442.005516-0.00000071.782294-0.000000-71.782294-0.000000-5.009149-0.0000002.639890-1.948227-0.6916630.000000-4.0576163.18176023975950.8582193.6917380.459192-5.0091493.630791-77.0736773.578898-0.679393-2.1802445.4197060.0000002.4011012.875049-37.222270第3次迭代的功率偏差dP和dQ-1.602219e-002-2.059526e-0031.261084e-002-9.148976e-003-3.913O38e-OO2-4.084703e-002-2.1571