潮流计算的计算机算法实验报告
姓名：学号：班级：
一、实验目的
掌握潮流计算的计算机算法。
熟悉MATLAB，并掌握MATLAB程序的基本调试方法。
二、实验准备
根据课程内容，熟悉MATLAB软件的使用方法，自行学习MATLAB程序的基础语法，并根据所学知识编写潮流计算牛顿拉夫逊法(或PQ分解法) 的计算程序，用相应的算例在MATLAB上进行计算、调试和验证。
三、实验要求
每人一组，在实验课时内，用MATLAB调试和修改运行程序，用算例计算输出潮流结果。
四、实验程序
clear; %清空内存
n=input('请输入节点数:n=');
n1=input('请输入支路数:n1=');
isb=input('请输入平衡节点号:isb=');
pr=input('请输入误差精度:pr=');
B1=input('请输入支路参数:B1=');
B2=input('请输入节点参数:B2=');
X=input('节点号和对地参数:X=');
Y=zeros(n);
Times=1;
%一：创建节点导纳矩阵
for i=1:n1
if B1(i,6)==0 %不含变压器的支路
p=B1(i,1);
q=B1(i,2);
Y(p,q)=Y(p,q)-1/B1(i,3);
Y(q,p)=Y(p,q);
Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4);
Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4);
else %含有变压器的支路
p=B1(i,1);
q=B1(i,2);
Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5));
Y(q,p)=Y(p,q);
Y(p,p)=Y(p,p)+1/B1(i,3);
Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3));
end
end
Y;
%将OrgS、DetaS初始化
OrgS=zeros(2*n-2,1);
DetaS=zeros(2*n-2,1);
%二：创建OrgS，用于存储初始功率参数
h=0;
j=0;
for i=1:n %对PQ节点的处理
if i~=isb&B2(i,6)==2
h=h+1;
for j=1:n
OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j ,3))-imag(Y(i,j))*I
mag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i, j))*real(B2(j,3)));
OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3)) -imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2( j,3))+imag(Y(i,j))*real(B2(j,3)));
end
end
end
%三：对PV节点的处理，注意这时不可再将h初始化为0
for i=1:n
if i~=isb&B2(i,6)==3
h=h+1;
for j=1:n
OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j ,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag (B2(j,3))+imag(Y(i,j))*real(B2(j,3)));
OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3)) -imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2( j,3))+imag(Y(i,j))*real(B2(j,3)));
end
end
end
OrgS;
%四：创建PVU 用于存储PV节点的初始电压
PVU=zeros(n-h-1,1);
t=0;
for i=1:n
if B2(i,6)==3
t=t+1;
PVU(t,1)=B2(i,3);
end
end
PVU;
%五：创建DetaS,用于存储有功功率、无功功率和电压幅值的不平衡量h=0;
for i=1:n %对PQ节点的处理
if i~=isb&B2(i,6)==2
h=h+1;
DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);
DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);
end
end
t=0;
for i=1:n
%六：对PV节点的处理，注意这时不可再将h初始化为0
if i~=isb&B2(i,6)==3
h=h+1;
t=t+1;
DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);
DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,
3))^2-imag( B2(i,3))^2;
end
end
DetaS;
%七：创建I，用于存储节点电流参数
i=zeros(n-1,1);
h=0;
for i=1:n
if i~=isb
h=h+1;
I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));
end
end
I;
%八：创建Jacbi(雅可比矩阵)
Jacbi=zeros(2*n-2);
h=0;
k=0;
for i=1:n %对PQ节点的处理
if B2(i,6)==2
h=h+1;
for j=1:n
if j~=isb
k=k+1;
if i==j %对角元素的处理
Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag( B2(i,3))+imag(I(h,1));
Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2( i,3))+real(I(h,1));
Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));
Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));
else %非对角元素的处理
Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag( B2(i,3));
Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2( i,3));
Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);
Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);
end
if k==(n-1) %将用于内循环的指针置于初始值，以确保雅可比矩阵换行
k=0;
end
end
end
end
end
k=0;
for i=1:n %对PV节点的处理
if B2(i,6)==3
h=h+1;
for j=1:n
if j~=isb
k=k+1;
if i==j %对角元素的处理
Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag( B2(i,3))+imag(I(h,1));
Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2( i,3))+real(I(h,1));
Jacbi(2*h,2*k-1)=2*imag(B2(i,3));
Jacbi(2*h,2*k)=2*real(B2(i,3));
else %非对角元素的处理
Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag( B2(i,3));
Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2( i,3));
Jacbi(2*h,2*k-1)=0;
Jacbi(2*h,2*k)=0;
end
if k==(n-1) %将用于内循环的指针置于初始值，以确保雅可比矩阵换行
k=0;
end
end
end
end
end
Jacbi;
%九：求解修正方程，获取节点电压的不平衡量
DetaU=zeros(2*n-2,1);
DetaU=inv(Jacbi)*DetaS;
DetaU;
%修正节点电压
j=0;
for i=1:n %对PQ节点处理
if B2(i,6)==2
j=j+1;
B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);
end
end
for i=1:n %对PV节点的处理
if B2(i,6)==3
j=j+1;
B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);
end
end
B2;
%十：开始循环
********************************************************************* *
while abs(max(DetaU))>pr
OrgS=zeros(2*n-2,1); %!!!初始功率参数在迭代过程中是不累加的，所以在这里必须将其初始化为零矩阵
h=0;
j=0;
for i=1:n
if i~=isb&B2(i,6)==2
h=h+1;
for j=1:n
OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j ,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag (B2(j,3))+imag(Y(i,j))*real(B2(j,3)));
OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3)) -imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2( j,3))+imag(Y(i,j))*real(B2(j,3)));
end
end
end
for i=1:n
if i~=isb&B2(i,6)==3
h=h+1;
for j=1:n
OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j ,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag (B2(j,3))+imag(Y(i,j))*real(B2(j,3)));
OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3)) -imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2( j,3))+imag(Y(i,j))*real(B2(j,3)));
end
end
end
OrgS;
%创建DetaS
h=0;
for i=1:n
if i~=isb&B2(i,6)==2
h=h+1;
DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);
DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);
end
end
t=0;
for i=1:n
if i~=isb&B2(i,6)==3
h=h+1;
t=t+1;
DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1); DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-im ag(B2(i,3))^2;
end
end
DetaS;
%创建I
i=zeros(n-1,1);
h=0;
for i=1:n
if i~=isb
h=h+1;
I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));
end
end
I;
%创建Jacbi
Jacbi=zeros(2*n-2);
h=0;
k=0;
for i=1:n
if B2(i,6)==2
h=h+1;
for j=1:n
if j~=isb
k=k+1;
if i==j
Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag( B2(i,3))+imag(I(h,1));
Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2( i,3))+real(I(h,1));
Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));
Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));
else
Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag( B2(i,3));
Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2( i,3));
Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);
Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);
end
if k==(n-1)
k=0;
end
end
end
end
end
k=0;
for i=1:n
if B2(i,6)==3
h=h+1;
for j=1:n
if j~=isb
k=k+1;
if i==j
Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag( B2(i,3))+imag(I(h,1));
Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2( i,3))+real(I(h,1));
Jacbi(2*h,2*k-1)=2*imag(B2(i,3));
Jacbi(2*h,2*k)=2*real(B2(i,3));
else
Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag( B2(i,3));
Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2( i,3));
Jacbi(2*h,2*k-1)=0;
Jacbi(2*h,2*k)=0;
end
if k==(n-1)
k=0;
end
end
end
end
end
Jacbi;
DetaU=zeros(2*n-2,1);
DetaU=inv(Jacbi)*DetaS;
DetaU;
%修正节点电压
j=0;
for i=1:n
if B2(i,6)==2
j=j+1;
B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);
end
end
for i=1:n
if B2(i,6)==3
j=j+1;
B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);
end
end
B2;
Times=Times+1; %迭代次数加1
end
Times;
五、实验流程
六、实验结果
参数输入：
运行结果：
七、实验体会
通过这次实验，让我第一次接触到了MATLAB,并深切体会到了它的强大之处；潮流计算的计算机算法的实现不仅巩固了我的学过的知识，还让我学到一些MATLAB的编程，虽然在实验的过程中出现了很多的错误，但在老师的细心指导下，问题都解决啦；计算机为我们省去了大量的人工计算，希望在以后的学习中能接触到更多的软件，学习到更多的知识。
第四章潮流计算的计算机算法 第一节概述 潮流计算是电力系统最基本、最常用的计算。根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线的电压（幅值及相角），各元件中流过的功率、整个系统的功率损耗等。潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。因此潮流计算在电力系统的规划设计、生产运行、调度管理及科学研究中都有着广泛的应用。 电力系统潮流计算分为离线潮流计算和在线潮流计算。前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。本章主要讨论离线潮流计算问题，它的基本算法同样适用于在线潮流计算。 潮流计算在数学上是多元非线性方程组的求解问题，求解的方法有很多种。自从五十年代计算机应用于电力系统以来，当时求解潮流的方法是以节点导纳矩阵为基础的逐次代入法（导纳法），后来为解决导纳法的收敛性较差的问题，出现了以阻抗矩阵为基础的逐次代入法（阻抗法）。到六十年代，针对阻抗法占用计算机内存大的问题又出现了分块阻抗法及牛顿－拉夫逊（Newton－Raphson）法。Newton —Raphson法是数学上解非线形方程式的有效方法，有较好的收敛性。将N－R法用于潮流计算是以导纳矩阵为基础的，由于利用了导纳矩阵的对称性、稀疏性及节点编号顺序优化等技巧，使N－R法在收敛性、占用内存、计算速度方面的优点都超过了阻抗法，成为六十年代末期以后普遍采用的方法。同时国内外广泛研究了诸如非线形规划法、直流法、交流法等各种不同的潮流计算方法。七十年代以来，又涌现出了更新的潮流计算方法。其中有1974年由B、Stott、O、Alsac 提出的快速分解法以及1978年由岩本伸一等提出的保留非线性的高 129
哈尔滨理工大学毕业设计（论文）任务书 学生姓名：孙聪学号：0903010909 学院：电气与电子工程学院专业：电气工程及其自动化 任务起止时间：2013 年2月25 日至2013年 6 月20 日 毕业设计（论文）题目： 基于前推回代法的配电网潮流计算 毕业设计工作内容： 1、查阅国内外相关参考文献，要求阅读20篇以上文献，了解当今电力 系统的发展状况，及目前研究的热点问题； 2、复习并熟练掌握电力系统潮流计算步骤及计算过程； 3、自学前推回代法潮流计算的基本原理及过程； 4、熟悉C语言，编写配电网潮流计算程序； 5、通过实际算例验证所编写程序的可靠性和准确性； 6、撰写论文，准备答辩。 资料： 1、王守相，王成山．现代配电系统分析[M]．北京：高等教育出版社， 2007． 2、刘健，毕鹏翔，董海鹏．复杂配电网简化分析与优化[M]．北京：中 国电力出版社，2002． 3、何仰赞，温增银．电力系统分析(上册)(第三版)[M] ．武汉：华中科 技大学出版社，2002． 4、李光琦．电力系统暂态分析[M]．北京：中国电力出版社，1998． 指导教师意见： 签名： 年月日系主任意见： 签名： 年月日 教务处制表
基于前推回带法的配电网潮流计算的研究 摘要 电力系统的潮流计算在电力系统稳态分析和电力系统设计中有很重要的作用，潮流计算也是电力系统暂态分析的基础。潮流计算是根据给定的系统运行条件来计算系统各个部分的运行状况，主要包括电压和功率的计算。 配电网潮流计算是配电管理系统高级应用软件功能组成之一。本课题在分析配电网元件模型的基础上，建立了配电网潮流计算的数学模型。由于配电网的结构和参数与输电网有很大的区别，因此配电网的潮流计算必须采用相适应的算法。配电网的结构特点呈辐射状，在正常运行时是开环的；配电网的另一个特点是配电线路的总长度较输电线路要长且分支比较多，配电线路的线径比输电网细导致配电网的R/X较大，且线路的充电电容可以忽略。配电网的潮流计算采用的方法是前推回代法，文中对前推回代法的基本原理、收敛性及计算速度等进行了理论分析比较。经过C语言编程，运行算例表明，前推回代法具有编程简单、计算速度快、收敛性好的特点，此方法是配电网潮流计算的有效算法，具有很强的实用性。 关键词：电力系统；配电网；潮流计算；前推回代法
目录 摘要 (1) 1.任务及题目要求 (2) 2.计算原理 (3) 2.1牛顿—拉夫逊法简介 (3) 2.2牛顿—拉夫逊法的几何意义 (7) 3计算步骤 (7) 4.结果分析 (9) 小结 (11) 参考文献 (12) 附录：源程序 (13) 本科生课程设计成绩评定表 (32)
摘要 电力系统的出现，使高效，无污染，使用方便，易于调控的电能得到广泛应用，推动了社会生产各个领域的变化，开创了电力时代，发生率第二次技术革命。电力系统的规模和技术水准已经成为一个国家经济发展水平的标志之一。 电力系统稳态分析包括潮流计算和静态安全分析。电力系统潮流计算是电力系统最基本的计算，也是最重要的计算。所谓潮流计算，就是已知电网的接线方式与参数及运行条件，计算电力系统稳态运行各母线电压、个支路电流与功率及网损。对于正在运行的电力系统，通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。 在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。 关键词：电力系统潮流计算牛顿-拉夫逊法
潮流计算的基本算法及使用方法 一、 潮流计算的基本算法 1. 牛顿－拉夫逊法 1．1 概述 牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。这种方法的特点就是把对非线 性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。 牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏 导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。 1．2 一般概念 对于非线性代数方程组 ()0=x f 即 ()0,,,21=n i x x x f ()n i ,2,1= (1－1) 在待求量x 的某一个初始计算值() 0x 附件，将上式展开泰勒级数并略去二阶及以上的高 阶项，得到如下的线性化的方程组 ()()()() ()0000=?'+x x f x f (1－2) 上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量 ()() ()[]()()0 1 00x f x f x -'-=? (1－3) 将() 0x ?和() 0x 相加，得到变量的第一次改进值()1x 。接着再从() 1x 出发，重复上述计算 过程。因此从一定的初值() 0x 出发，应用牛顿法求解的迭代格式为 ()()()()() k k k x f x x f -=?' (1－4) ()()()k k k x x x ?+=+1 (1－5) 上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代
题目：潮流计算与matlab 教学单位电气信息学院姓名 学号 年级 专业电气工程及其自动化指导教师 职称副教授
摘要 电力系统稳态分析包括潮流计算和静态安全分析。本文主要运用的事潮流计算，潮流计算是电力网络设计与运行中最基本的运算，对电力网络的各种设计方案及各种运行方式进行潮流计算，可以得到各种电网各节点的电压，并求得网络的潮流及网络中的各元件的电力损耗，进而求得电能损耗。本位就是运用潮流计算具体分析，并有MATLAB仿真。 关键词：电力系统潮流计算 MATLAB Abstract Electric power system steady flow calculation and analysis of the static safety analysis. This paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, various design scheme and the operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the components of the power loss, and getting electric power. The standard is to use the power flow calculation and analysis, the specific have MATLAB simulation. Key words: Power system; Flow calculation; MATLAB simulation
目录 摘要................................................. - 1 - 1.设计意义与要求..................................... - 2 - 1.1设计意义 ...................................... - 2 - 1.2设计要求（具体题目）........................... - 2 - 2.题目解析........................................... - 3 - 2.1设计思路 ...................................... - 3 - 2.2详细设计 ...................................... - 4 - 2.2.1节点类型.................................. - 4 - 2.2.2待求量 ................................... - 4 - 2.2.3导纳矩阵.................................. - 4 - 2.2.4潮流方程.................................. - 5 - 2.2.5牛顿—拉夫逊算法.......................... - 6 - 2.2.5.1牛顿算法数学原理：................... - 6 - 2.2.5.2修正方程............................. - 7 - 2.2.5.3收敛条件............................. - 9 - 3.结果分析.......................................... - 10 - 4.小结.............................................. - 11 - 参考文献............................................ - 12 -
由于本人参加我们电气学院的电气小课堂，主讲的是计算机算法计算潮流这章，所以潜心玩了一个星期，下面整理给大家分享下。 本人一个星期以来的汗水，弄清楚了计算机算法计算潮流的基础，如果有什么不懂的可以发信息到邮箱：zenghao616@https://www.doczj.com/doc/9410687360.html, 接下来开始弄潮流的优化问题，吼吼！ 电力系统的潮流计算的计算机算法：以MATLAB为环境 这里理论不做过多介绍，推荐一本专门讲解电力系统分析的计算机算法的书籍---------《电力系统分析的计算机算法》—邱晓燕、刘天琪编著。 这里以这本书上的例题【2-1】说明计算机算法计算的过程，分别是牛顿拉弗逊算法的直角坐标和极坐标算法、P-Q分解算法。主要是简单的网络的潮流计算，其实简单网络计算和大型网络计算并无本质区别，代码里面只需要修改循环迭代的N即可，这里旨在弄清计算机算法计算潮流的本质。代码均有详细的注释. 其中简单的高斯赛德尔迭代法是以我们的电稳教材为例子讲，其实都差不多，只要把导纳矩阵Y给你，节点的编号和分类给你，就可以进行计算了，不必要找到原始的电气接线图。 理论不多说，直接上代码： 简单的高斯赛德尔迭代法： 这里我们只是迭代算出各个节点的电压值，支路功率并没有计算。 S_ij=P_ij+Q_ij=V_i(V_i* - V_j*) * y_ij* 可以计算出各个线路的功率 在显示最终电压幅角的时候注意在MATLAB里面默认的是弧度的形式，需要转化成角度显示。 clear;clc; %电稳书Page 102 例题3-5 %计算网络的潮流分布 --- 高斯-赛德尔算法 %其中节点1是平衡节点 %节点2、3是PV节点，其余是PQ节点 % 如果节点有对地导纳支路 %需将对地导纳支路算到自导纳里面 %------------------------------------------------% %输入原始数据，每条支路的导纳数值，包括自导和互导纳； y=zeros(5,5); y(1,2)=1/(0.0194+0.0592*1i); y(1,5)=1/(0.054+0.223*1i); y(2,3)=1/(0.04699+0.198*1i); y(2,4)=1/(0.0581+0.1763*1i); %由于电路网络的互易性，导纳矩阵为对称的矩阵
摘要 配电网潮流计算是配电管理系统应用软件功能组成之一。本设计在分析配电网元件模型的基础上，建立了配电网潮流计算的数学模型。由于配电网的结构参数与输电网有很大的区别，因此配电网的潮流计算采用相适应的算法。配电网的结构特点呈辐射状，在正常运行时是开环的；配电网的另一个特点是配电线路的总长度较输电线路要长并且分支较多，配电线路的线径比输电网的细以至于配电网的R/X较大，且线路的充电电容可以忽略。配电网的潮流计算采用的方法是前推回代法，文中对前推回代法的基本原理，收敛性及计算速度等进行了理论分析比较仿真和算例表明，前推回代法具有编程简单、计算速度快、收敛性好的特点，这个方法是配电网潮流计算的有效算法，具有很强的实用性。 关键词配电网，潮流计算，前推回代法
Abstract Flow solution of distribution networks is one of software in DMS. Because of the different structures between transmission networks and distribution networks, the corresponding methods in flow solution of distribution networks must be applied. Distributions network is radial shape and in the condition of regular is annular. Another characteristic of distribution networks is cabinet minister of distribution long than transmission networks. The line diameter of distribution networks is thin than transmission networks, it cause R/X is large of distribution networks and the line’s capacitance can neglect. Load flow calculation of distributions network use back/ forward sweep. It has some peculiarities such as simple procedures and good restrain and so on. This method of distribution network is an effective method of calculating the trend, with some practicality. Key words :distribution network，load flow calculation，back/ forward sweep
电力系统分析潮流计算报告
目录 一．配电网概述 (3) 1.1 配电网的分类 (3) 1.2 配电网运行的特点及要求 (3) 1.3 配电网潮流计算的意义 (4) 二．计算原理及计算流程 (4) 2.1 前推回代法计算原理 (4) 2.2 前推回代法计算流程 (7) 2.3主程序清单： (9) 2.4 输入文件清单： (11) 2.5计算结果清单： (12) 三．前推回代法计算流程图 (13) 参考文献 (14)
一．配电网概述 1.1 配电网的分类 在电力网中重要起分配电能作用的网络就称为配电网； 配电网按电压等级来分类，可分为高压配电网（35—110KV），中压配电网（6—10KV，苏州有20KV的），低压配电网（220/380V）； 在负载率较大的特大型城市，220KV电网也有配电功能。 按供电区的功能来分类，可分为城市配电网，农村配电网和工厂配电网等。 在城市电网系统中,主网是指110KV及其以上电压等级的电网，主要起连接区域高压（220KV及以上）电网的作用。 配电网是指35KV及其以下电压等级的电网，作用是给城市里各个配电站和各类用电负荷供给电源。 从投资角度看，我国与国外先进国家的发电、输电、配电投资比率差异很大，国外基本上是电网投资大于电厂投资，输电投资小于配电投资。我国刚从重发电轻供电状态中转变过来，而在供电投资中，输电投资大于配电投资。从我国城网改造之后，将逐渐从输电投资转入配电建设为主。 本文是基于前推回代法的配电网潮流分析计算的研究，研究是是以根节点为10kV的电压等级的配电网。 1.2 配电网运行的特点及要求 配电系统相对于输电系统来说，由于电压等级低、供电范围小，但与用户直接相连，是供电部门对用户服务的窗口，因而决定了配电网运行有如下特点和基本要求：
毕业设计（论文）题目配电网潮流计算与程序设计 完成日期年月日
目录 1绪论 (4) 1.1潮流计算的目的和意义 (4) 1.2潮流计算的发展 (4) 1.3本文的主要工作 (6) 2 配电网络模型 (6) 2.1线路模型 (6) 2.2变压器的模型 (9) 2.3 负荷模型 (13) 3基于前推回代法的潮流计算 (14) 4基于matlab的程序设计 (15) 4.1 MATLAB简介 (15) 4.2 程序设计 (16) 4.3 算例 (20) 5 总结 (24)
配电网潮流计算与程序设计 摘要：在电力系统的正常运行中，随着用电负荷的变化和系统运行方式的改变，网络中的损耗也将发生变化。要严格保证所有的用户在任何时刻都有额定的电压是不可能的，因此系统运行中个节点出现电压的偏移是不可避免的。为了保证电力系统的稳定运行，要进行潮流调节。 Abstract:In the normal operation of the power system, with the change of power load and the change of system operation mode, the network losses will also change. To strictly ensure that all users at any time have a rated voltage is impossible, so the node appears in the operation of the system voltage deviation is inevitable. In order to guarantee the stable operation of power system, must carry on the trend of the regulation. 关键词: 电网潮流计算前推后代法MATLAB Keywords:Power flow calculation Push the offspring before MATLAB 前言 电力是以电能作为动力的能源。发明于19世纪70 年代，电力的发明和应用掀起了第二次工业化高潮。成为人类历史18世纪以来，世界发生的三次科技革命之一，从此科技改变了人们的生活。既是是当今的互联网时代我们仍然对电力有着持续增长的需求，因为我们发明了电脑、家电等更多使用电力的产品。不可否认新技术的不断出现使得电力成为人们的必需品。 20世纪出现的大规模电力系统是人类工程科学史上最重要的成就之一，是由发电、输电、变电、配电和用电等环节组成的电力生产与消费系统。它的功能是将自然界的一次能源通过发电动力装置转化成电能，再经输电、变电和配电将电能供应到各用户。为实现这一功能，电力系统在各个环节和不同层次还具有相应的信息与控制系统，对电能的生产过程进行测量、调节、控制、保护、通信和调度，以保证用户获得安全、经济、优质的电能。 电力工业发展初期，电能是直接在用户附近的发电站（或称发电厂）中生产的，各发电站孤立运行。随着工农业生产和城市的发展，电能的需要量迅速增加，而热能资源和水能资源丰富的地区又往往远离用电比较集中的城市和工矿区，为了解决这个矛盾，就需要在动力资源丰富的地区建立大型发电站，然后将电能远距离输送给电力用户。同时，为了提高供电的可靠性以及资源利用的综合经济性，又把许多分散的各种形式的发电站，通过送电线路和变电所联系起来。这种由发电机、升压和降压变电所，送电线路以及用电设备有机连接起来的整体，即称为电力系统。 随着电力系统及在线应用的发展，计算机网络已经形成，为电力系统的潮流计算提供了物质基础。电力系统潮流计算是电力系统分析计算中最基本的内容，也是电力系统运行及设计中必不可少的工具。根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线电压的幅值及相角、各元件中流过
《电力系统潮流上机》课程设计报告 院系：电气工程学院 班级：电088班 学号： 0812002221 学生姓名：刘东昇 指导教师：张新松 设计周数：两周 日期：2010年 12 月 25 日
一、课程设计的目的与要求 目的：培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识 要求：基本要求： 1.编写潮流计算程序； 2.在计算机上调试通过； 3.运行程序并计算出正确结果； 4.写出课程设计报告 二、设计步骤： 1.根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 三、设计原理 1．牛顿-拉夫逊原理 牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为
额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤： （1）形成各节点导纳矩阵Y。 （2）设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。 （3）计算各个节点的功率不平衡量。 （4）根据收敛条件判断是否满足，若不满足则向下进行。 （5）计算雅可比矩阵中的各元素。 （6）修正方程式个节点电压 （7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。 （8）计算平衡节点输出功率和各线路功率 2．网络节点的优化 １）静态地按最少出线支路数编号 这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 四、设计内容
潮流计算数学模型与数值方法 1. 什么是潮流计算潮流计算的主要作用有哪些 潮流计算，电力学名词，指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下，计算有功功率、无功功率及电压在电力网中的分布。 潮流计算是电力系统非常重要的分析计算，用以研究系统规划和运行中提出的各种问题。对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求；对运行中的电力系统，通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件(线路、变压器等)是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。 2. 潮流计算有哪些待求量、已知量 （已知量：1、电力系统网络结构、参数 2、决定系统运行状态的边界条件 待求量：系统稳态运行状态 例如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等） 3. 潮流计算节点分成哪几类分类根据是什么 （分成三类：PQ 节点、PV 节点和平衡节点，分类依据是给定变量的不同） 4. 教材牛顿-拉夫逊法及有功-无功分解法是基于何种电路方程可否采用其它类型方程 答:基于节点电压方程，还可以采用回路电流方程和割集电压方程等。但是后两者不常用。 5. 教材牛顿-拉夫逊法是基于节点阻抗方程、还是基于节点导纳方程进行迭代计算的试阐述这两种方程的优点与缺点。 1.不能由等值电路直接求出 2.满秩矩阵内存量大 3.对角占优矩阵。。 节点导纳矩阵的特点：1.直观容易形成2.对称阵3.稀疏矩阵（零元素多）：每一行的零元素个数=该节点直接连出的支路数。 6. 说出至少两种建立节点导纳矩阵的方法，阐述其中一种方法的原理与过程。 方法：1.根据自导纳和互导纳的定义直接求取2.运用一节点关联矩阵计算3.阻抗矩阵的逆矩阵 节点导纳矩阵的形成：1.对角线元素ii Y 的求解)1,,0(=≠==i j I i ii U i j U U I Y 【除i 外的其他节点接地，0=j U ，只在i 节点加单位电压值】解析ii Y 等于与i 节点直接相连的的所有支路导纳和2.互导纳),0,1(j k U U U I Y k j j i ij ≠===，ji ij Y Y =（无源网络导纳之间是对称的）解析：ij Y 等于j i ,节点之间直接相连的支路导纳的负值。 7. 潮流计算需要考虑哪些约束条件 答: 为了保证系统的正常运行必须满足以下的约束条件： 对控制变量
电力系统分析课程设计 设计题目9节点电力网络潮流计算 指导教师 院（系、部）电气与控制工程学院 专业班级 学号 姓名 日期
电气工程系课程设计标准评分模板
目录 1 PSASP软件简介 (1) 1.1 PSASP平台的主要功能和特点 (1) 1.2 PSASP的平台组成 (2) 2 牛顿拉夫逊潮流计算简介 (3) 2.1 牛顿—拉夫逊法概要 (3) 2.2 直角坐标下的牛顿—拉夫逊潮流计算 (5) 2.3 牛顿—拉夫逊潮流计算的方法 (6) 3 九节点系统单线图及元件数据 (7) 3.1 九节点系统单线图 (7) 3.2 系统各项元件的数据 (8) 4 潮流计算的结果 (10) 4.1 潮流计算后的单线图 (10) 4.2 潮流计算结果输出表格 (10) 5 结论 (14)
电力系统分析课程设计任务书9节点系统单线图如下： 基本数据如下：
表3 两绕组变压器数据 负荷数据
1 PSASP软件简介 “电力系统分析综合程序”（Power System Analysis Software Package,PSASP）是一套历史悠久、功能强大、使用方便的电力系统分析程序，是高度集成和开发具有我国自主知识产权的大型软件包。 基于电网基础数据库、固定模型库以及用户自定义模型库的支持，PSASP可进行电力系统（输电、供电和配电系统）的各种计算分析，目前包括十多个计算机模块，PSASP的计算功能还在不断发展、完善和扩充。 为了便于用户使用以及程序功能扩充，在PSASP7.0中设计和开发了图模一体化支持平台，应用该平台可以方便地建立电网分析的各种数据，绘制所需要的各种电网图形（单线图、地理位置接线图、厂站主接线图等）；该平台服务于PSASP 的各种计算，在此之外可以进行各种分析计算，并输出各种计算结果。 1.1PSASP平台的主要功能和特点 PSASP图模一体化支持平台的主要功能和特点可概括为： 1. 图模支持平台具备MDI多文档操作界面，是一个单线图图形绘制、元件数据录入编辑、各种计算功能、结果显示、报表和曲线输出的集成环境。用户可以方便地建立电网数据、绘制电网图形、惊醒各种分析计算。人机交互界面全部汉化，界面良好，操作方便。 2. 真正的实现了图模一体化。可边绘图边建数据，也可以在数据已知的情况下进行图形自动快速绘制；图形、数据自动对应，所见即所得。 3. 应用该平台可以绘制各种电网图形，包括单线图、地理位置接线图、厂站主接线图等。 ●所有图形独立于各种分析计算，并为各计算模块所共享； ●可在图形上进行各种计算操作，并在图上显示各种计算结果； ●同一系统可对应多套单线图，多层子图嵌套； ●单线图上可细化到厂站主接线结构；
3.3复杂电力网潮流计算的计算机解法 3.3.1 导纳矩阵的形成 1．自导纳 节点i的自导纳，亦称输入导纳，在数值上等于在节点i施加单位电压，其他节点全部接地时，经节点i注入网络的电流。 主对角线元素，更具体地说，就等于与节点连接的所有支路导纳的和。 2．互导纳 节点i、j间的互导纳，在数值上等于在节点i施加单位电压，其他节点全部接地时，经节点j注入网络的电流。非对角线元素。 更具体地说，是连接节点j和节点i支路的导纳之和再加上负号而得。 3．导钠矩阵的特点： （1）因为，导纳矩阵Y是对称矩阵； （2）导纳矩阵是稀疏矩阵，每一非对角元素是节点i和j间支路导纳的负值，当i和j间没有直接 相连的支路时，即为零，根据一般电力系统的特点，每一节点平均与3-5个相邻节点有直接联系，所以导纳矩阵是一高度稀疏的矩阵； （3）导纳矩阵能从系统网络接线图直观地求出。 4．节点导纳矩阵的修改 （1）从原有网络引出一支路，同时增加一节点，设i为原有网络结点，j为新增节点，新增支路ij的导纳为y ij。如图3-17（a）所示。 因新增一节点，新的节点导纳阵需增加一阶。且新增对角元Y jj=y ij，新增非对角元Y ij=Y ji=－y ij，同时对原阵中的对角元Y ii进行修改，增加ΔY ii＝y ij。 （2）在原有网络节点i、j间增加一支路。如图3-17（b）所示。 设在节点i增加一条支路，由于没有增加节点数，节点导纳矩阵Y阶次不变，节点的自导纳Y ii、Y jj和互导纳Y ij分别变化量为 （3-57） 图 3-17 网络接线的变化图 （a）网络引出一支路，（b）节点间增加一支路，（c）节点间切除一支路，（d）节点间导纳改变 （3）在原有网络节点i、j间切除一支路。如图3-17（c）所示。 设在节点i切除一条支路，由于没有增加节点数，节点导纳矩阵Y阶次不变，节点的自导纳Y ii、Y jj和互导纳Y ij分别发生变化，其变化量为
3 简单电力系统潮流计算 3．1 思考题、习题 1）电力线路阻抗中的功率损耗表达式是什么？电力线路始、末端的电容功率表达式是什 么？ 2）电力线路阻抗中电压降落的纵分量和横分量的表达式是什么？ 3）什么叫电压降落、电压损耗、电压偏移、电压调整及输电效率？ 5）对简单开式网络、变电所较多的开式网络和环形网络潮流计算的内容及步骤是什么？ 6）变压器在额定状况下，其功率损耗的简单表达式是什么？ 9）为什么要对电力网络的潮流进行调整控制？调整控制潮流的手段主要有哪些？ 10）欲改变电力网络的有功功率和无功功率分布，分别需要调整网络的什么参数？ 16）110kV 双回架空线路，长度为150kM ，导线型号为LGJ-120，导线计算外径为15.2mm ， 三相导线几何平均距离为5m 。已知电力线路末端负荷为30+j15MVA ，末端电压为106kV ，求 始端电压、功率，并作出电压向量图。 17）220kV 单回架空线路，长度为200kM ，导线型号为LGJ-300，导线计算外径为24.2mm ， 三相导线几何平均距离为7.5m 。已知电力线路始端输入功率为120+j50MVA ，始端电压为 240kV ，求末端电压、功率，并作出电压向量图。 18）110kV 单回架空线路，长度为80kM ，导线型号为LGJ-95，导线计算外径为13.7mm ，三 相导线几何平均距离为5m 。已知电力线路末端负荷为15+j10MVA ，始端电压为116kV ，求末 端电压和始端功率。 19）220kV 单回架空线路，长度为220kM ，电力线路每公里的参数分别为： kM S b kM x kM r /1066.2,/42.0,/108.06111-?=Ω=Ω=、 线路空载运行，当线路末端电压为205kV ，求线路始端的电压。 20）有一台三绕组变压器，其归算至高压侧的等值电路如图3-1所示，其中 ,68~,45~,8.3747.2,5.147.2,6547.232321MVA j S MVA j S j Z j Z j Z T T T +=+=Ω+=Ω-=Ω+=当变压器变比 为110/38.5（1+5%）/6.6kV ，U 3=6kV 时，试计算高压、中压侧的实际电压。 图3- 1 图3-2
第三章 电力系统的潮流计算 3-1 电力系统潮流计算就是对给定的系统运行条件确定系统的运行状态。系 统运行条件是指发电机组发出的有功功率和无功功率（或极端电压），负荷的有 功功率和无功功率等。运行状态是指系统中所有母线（或称节点）电压的幅值和 相位，所有线路的功率分布和功率损耗等。 3-2 电压降落是指元件首末端两点电压的相量差。 电压损耗是两点间电压绝对值之差。当两点电压之间的相角差不大时， 可以近似地认为电压损耗等于电压降落的纵分量。 电压偏移是指网络中某点的实际电压同网络该处的额定电压之差。电压 偏移可以用kV 表示，也可以用额定电压的百分数表示。 电压偏移= %100?-N N V V V 功率损耗包括电流通过元件的电阻和等值电抗时产生的功率损耗和电压 施加于元件的对地等值导纳时产生的损耗。 输电效率是是线路末端输出的有功功率2P 与线路首端输入的有功功率 1P 之比。 输电效率= %1001 2 ?P P 3-3 网络元件的电压降落可以表示为 ()? ? ? ? ? +=+=-2221V V I jX R V V δ? 式中，?2V ?和? 2V δ分别称为电压降落的纵分量和横分量。 从电压降落的公式可见，不论从元件的哪一端计算，电压降落的纵、横分量计算公式的结构都是一样的，元件两端的电压幅值差主要有电压降落的纵分量决定，电压的相角差则由横分量决定。在高压输电线路中，电抗要远远大于电阻，即R X ??，作为极端的情况，令0=R ，便得 V QX V /=?，V PX V /=δ 上式说明，在纯电抗元件中，电压降落的纵分量是因传送无功功率而产生的，而电压降落的横分量则是因为传送有功功率产生的。换句话说，元件两端存在电压幅值差是传送无功功率的条件，存在电压相角差则是传送有功功率的条件。 3-4 求解已知首端电压和末端功率潮流计算问题的思路是，将该问题转化成 已知同侧电压和功率的潮流计算问题。
毕业设计（论文）题，目配电网潮流计算与程序设计 学生姓名石昊晨学号2010151107 专业发电厂及电力系统班级20109091 指导教师刘会家 评阅教师 完成日期年月日
目录 摘要 一．配电网潮流概述 (5) 1.1配电网潮流计算的目的与意义 (5) 1.2潮流计算方法概述 (5) 1.2.1 牛顿——拉夫逊法 (6) 1.2.2 快速解耦法 (6) 1.2.3 回路阻抗法 (9) 1.2.4 前推回代法 (11) 1.3 本文工作 (11) 二．配电网网络模型 (11) 2.1元件模型 (11) 2.1.1 电力线路的数学模型 (11) 2.1.2 变压器的等值电路 (13) 2.2网络模型 (15) 三：基于matlab的配电网潮流计算算法 (16) 3.1配电网潮流计算算法原理 (16) 3.2 matlab的概述 (19) 3.3程序设计 (21) 3.3.1 牛顿--拉夫逊法潮流求解过程 (21) 3.3.2牛顿—拉夫逊法的程序框图 (25) 四：算例 (27) 参考文献 (28) 致谢 (29)
配电网潮流计算与程序设计 学生：石昊晨 指导教师：刘会家 （三峡大学国际文化交流学院） 摘要：本文首先分析了配电网的特点及对算法的要求，然后建立配电网潮流计算模型。针对配电网潮流计算的现状进行了全面分析，深入讨论了目前各方法的特点，并从收敛性及其他性能指标进行了比较分析；详细研究用的比较广泛的牛顿——拉夫逊法，并以广度优先顺序搜索策略作为理论基础。针对某地区配电网的具体情况，选取IOKV的配电网子系统进行潮流计算。利用MATLAB 2009a 进行了基于牛顿——拉夫逊法的配电网的潮流计算程序。由计算结果可知，该算法具有一定的优越性，软件的开发具有一定的实用性。 关键词：电力系统，配电网潮流，牛顿——拉夫逊法，MATLAB程序设计
信息工程学系 2011-2012学年度下学期电力系统分析课程设计 题目:电力系统潮流计算 专业：电气工程及其自动化 班级： 学号： 学生姓名： 指导教师：钟建伟 2012年3月10日
信息工程学院课程设计任务书
目录 1 任务提出与方案论证 (4) 1.1潮流计算的定义、用途和意义 (4) 1.2 运用软件仿真计算 (5) 2 总体设计 (7) 2.1潮流计算设计原始数据 (7) 2.2总体电路设计 (8) 3 详细设计 (10)
3.1数据计算 (10) 3.2 软件仿真 (14) 4 总结 (24) 5参考文献 (25)
1任务提出与方案论证 1.1潮流计算的定义、用途和意义 1.1.1潮流计算的定义 潮流计算，指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下，计算有功功率、无功功率及电压在电力网中的分布。潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。 1.1.2潮流计算的用途 电力系统潮流计算是电力系统最基本的计算，也是最重要的计算。所谓潮流计算，就是已知电网的接线方式与参数及运行条件，计算电力系统稳态运行各母线电压、个支路电流与功率及网损。对于正在运行的电力系统，通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。
配电网络的拓扑分析及潮流计算 李晨 在当前经济迅猛发展、供电日趋紧张的情况下，通过配电网络重构，充分发挥现有配电网的潜力，提高系统的安全性和经济性，具有很大的经济效益和社会效益。本文对配电网拓扑分析、对配电网络潮流计算作分析研究，应用MATLAB编程来验证并分析配电网结构特点。配电网的拓扑分析用树搜索法，并采用前推回代法进行潮流计算分析，通过树搜索形成网络拓扑表，然后利用前推回代法计算潮流分布。 1 配电网的接线分析 配电网是指电力系统中二次降压侧直接或降压后向用户供电的网络。配电网由馈线、降压变压器、断路器、各种开关构成。就我国电力系统而言，配电网是指110kV及以下的电网。在配电网中，通常把110kV，35kV级称为高压，10kV级称为中压，0.4kV级称为低压。从体系结构上，配电网可以分作辐射状网、树状网和环状网，如图2.3所示。我国配电网大部分是呈树状结构。 辐射网树状网环状网 图1-1配电网的体系结构 1.1 配电网的支路节点编号 通过简化可把一个复杂的配电网络简化成一个节点一边关系的树状网络，于是就可以运行图论的知识进行网络拓扑分析。按照这种简化模型，易知：节点数目比支路数目和开关数目多1，所以节点从0开始编号，而支路数和开关数从1开始编号，这样编号三者在序号上就可以完全一致，为后面的网损计算打下良好的基础。联络线支路和上面的联络开关编号放在最后处理。 图1-2节点支路编号示意图 图中①为节点号，1为支路号，其它节点、支路编号的含义相同。 节点、支路编号原则：将根节点编为0，并按父节点小于子节点号的原则由根节点向下顺序编号，规定去路正方向为父节点指向子节点，且支路编号与其子节点同号，则网络结构