资源描述：
《牛顿-拉夫逊迭代法电力网潮流计算方法与程序》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、6牛顿-拉夫逊迭代法电力网潮流计算方法与程序编写佘名寰牛顿-拉夫逊迭代法在电力网潮流计算中因其收敛性较好获得广泛运用，该算法的难点是需反复计算功率方程中雅可比矩阵各个元素表达式。本文简叙了牛顿-拉夫逊迭代法潮流计算的基本公式，通过例题介绍了用牛顿-拉夫逊法计算电力网潮流电压的MATLAB程序。程序采用MATLAB语言的符号矩阵简化了雅可比矩阵系数的计算。本文可供电力系统电气技术人员和大专院校电力类专业师生参考。2．牛顿-拉夫逊迭代法潮流计算2.1极坐标型式节点功率方程式由电源注入或从负载流出节点的电流统称节点电流，由节点电压和节点电流所求得的节点
2、注入功率为：Sti=ViIiIi=YijVj(I,j=1,2,….n)Pti=Vi1nVjGijcosθij+BijsinθijQti=Vi1nVjGijsinθij-Bijcosθij(i,j=1,2…..n)（2-01）式中Pti,QtI——节点注入的有功和无功功率Vi，Vj——节点i,j电压幅值θij=θi-θj节点i,j电压的相角差Gij,Bij节点导纳矩阵的元素，YIJ=GIJ+jBIJ节点功率平衡关系为：Pgi-PlI=PtiQgi-Qli=QtiΔPi=Pgi-PlI-Pti=0ΔQi=Qgi-Qli-Qti=0（2-02）Pgi,
3、Qgi ——节点i发电机输入有功和无功功率Pli,Qli——节点i负荷有功和无功功率ΔPi，ΔQi--节点i不平衡功率不平衡功率的微分d(ΔPi),d(ΔQi)d∆pi=-(∂pti∂v1∆v1+∂pti∂v2∆v2…+∂pti∂ϑ1∆ϑ1+∂pti∂ϑ2∆ϑ2….)（i=1,2,…,n）d∆qi=-(∂qti∂v1∆v1+∂q∂v2∆v2…+∂qti∂ϑ1∆ϑ1+∂qti∂ϑ2∆ϑ2….)（2-03）对于n个节点系统可得如下矩阵形式修正方程式；6∆P1∆P2:∆Pn∆Q1∆Q2:∆Qn=∂∆P1∂ϑ1∂∆P2∂ϑ1:∂∆Pn∂ϑ1∂∆Q1∂ϑ1
4、∂∆Q2∂ϑ1:∂∆Qn∂ϑ1∂∆P1∂ϑ2∂∆P2∂ϑ2:∂∆Pn∂ϑ2∂∆Q1∂ϑ2∂∆Q2∂ϑ2:∂∆Qn∂ϑ2⋯⋯⋯⋯⋯⋯⋯⋯∂∆P1∂ϑn∂∆P2∂ϑn:∂∆Pn∂ϑn∂∆Q1∂ϑn∂∆Q2∂ϑn:∂∆Qn∂ϑn∂∆P1∂v1∂∆P2∂v1:∂∆Pn∂v1∂∆Q1∂v1∂∆Q2∂v1:∂∆Qn∂v1∂∆P1∂v2∂∆P2∂v2:∂∆Pn∂v2∂∆Q1∂v2∂∆Q2∂v2:∂∆Qn∂v2⋯⋯⋯⋯⋯⋯⋯⋯∂∆P1∂vn∂∆P2∂vn:∂∆Pn∂vn∂∆Q1∂vn∂∆Q2∂vn:∂∆Qn∂vnΔϑ1∆ϑ2:∆ϑn∆v1∆v2:∆vn（2-04
5、）式中偏微分矩阵为雅可比矩阵。若系统总共有n个节点，其中r个为PV节点，去掉平衡节点则电压角度θ变量有n-1-个，电压幅值v变量有n-1-r个，方程式共有2n-r-2个。2.2牛顿-拉夫逊迭代法解算过程；⑴指定各节点电压的初值；v1,v2,…vn;ϑ1,ϑ2,…ϑn.一般设V=1,θ=0；⑵将所给定的电压值代入式（2-01），（2-02）计算各节点不平衡功率；⑶将电压初值代入求雅可比矩阵系数；⑷解方程式（2-04），求节点电压幅值和相角的修正值∆vi,∆ϑi;⑸用所求得的电压修正值修正各节点电压：v(1)=v(0)+∆v(0)ϑ(1)=ϑ(0)+
6、∆ϑ(0)⑹用新的节点电压重复以上五个步骤计算直到计算精度合符要求。牛顿-拉夫逊迭代法计算的难点是求雅可比矩阵系数，本文采用MATLAB语言符号矩阵运算jacobian(f,x)命令求解函数f(x)对x的雅可比矩阵。2.3牛顿-拉夫逊迭代法计算程序我们通过例题说明牛顿-拉夫逊迭代法计算潮流程序。【例2.3.1】网络接线如图2-1所示。支路阻抗分别为Z12=j0.1,Z23=j0.1,Z31=j0.1;三条支路两端的对地电纳皆是j0.01。给定的注入功率分别为：S1=-1.7192-j0.7346;S2=0.6661。节点1是负荷节点即PQ节点；节
7、点2给定的电压大小为V2=1.05,其属于PV节点；节点3是平衡节点，其电压保持定值，大小为V3=1.0。Z13Z12Z23YC3YC1YC1V3=1.0312YC2P1=1.7192Q1=0.7346P2=0.6661V2=1.05Z21=j0.1Z23=j0.1Z13=j0.1YC1=J0.02YC2=J0.02YC3=J0.026图2-1[例2.3.1]网络接线图2.3.1.1程序参变量说明Y=G+j*B节点导纳矩阵，G为电导，B为电纳；P节点有功功率给定值，正为发电机注入，负为负荷功率；q节点无功功率给定值，正为发电机注入，负为负荷功率；
8、u节点电压幅值的符号矩阵；delt节点电压相角的符号矩阵；pt节点有功功率符号表达式；qt节点无功功率符号表达式；pp节点有功功率不平衡
