%本程序的功能是用牛拉法进行潮流计算
%原理介绍详见鞠平著《电气工程》
%默认数据为鞠平著《电气工程》例8.4所示数据
%B1是支路参数矩阵
%第一列和第二列是节点编号。节点编号由小到大编写
%对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点编号
%第三列为支路的串列阻抗参数，含变压器支路此值为变压器短路电抗
%第四列为支路的对地导纳参数，含变压器支路此值不代入计算
%第五烈为含变压器支路的变压器的变比，变压器非标准电压比
%第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，“0”为不含有变压器
%B2为节点参数矩阵
%第一列为节点注入发电功率参数
%第二列为节点负荷功率参数
%第三列为节点电压参数
%第四列
%第五列
%第六列为节点类型参数，“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数
%X为节点号和对地参数矩阵
%第一列为节点编号
%第二列为节点对地参数
clear;
clc;
num=input('是否采用默认数据？(1-默认数据;2-手动输入)');
if num==1
n=4;
n1=4;
isb=4;
pr=0.00001;
B1=[1 2 0.1667i 0 0.8864 1;1 3 0.1302+0.2479i 0.0258i 1 0;1 4 0.1736+0.3306i 0.0344i 1 0;3 4 0.2603+0.4959i 0.0518i 1 0];
B2=[0 0 1 0 0 2;0 -0.5-0.3i 1 0 0 2;0.2 0 1.05 0 0 3;0 -0.15-0.1i 1.05 0 0 1];
X=[1 0;2 0.05i;3 0;4 0];
else
n=input('请输入节点数:n=');
n1=input('请输入支路数:n1=');
isb=input('请输入平衡节点号:isb=');
pr=input('请输入误差精度:pr=');
B1=input('请输入支路参数:B1=');
B2=input('请输入节点参数:B2=');
X=input('节点号和对地参数:X=');
end
Times=1; %迭代次数
%创建节点导纳矩阵
Y=zeros(n);
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)-B1(i,5)/B1(i,3);
Y(q,p)=Y(p,q);
Y(p,p)=Y(p,p)+B1(i,5)/B1(i,3)+(1-B1(i,5))/B1(i,3);
Y(q,q)=Y(q,q)+B1(i,5)/B1(i,3)+(B1(i,5)*(B1(i,5)-1))/B1(i,3);
end
end
for i=1:n1
Y(i,i)=Y(i,i)+X(i,2); %计及补偿电容电纳
end
disp('导纳矩阵为：');
disp(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 %不是平衡点&是PQ点
h=h+1;
for j=1:n
%公式8-74
%Pi=ei*(Gij*ej-Bij*fj)+fi*(Gij*fj+Bij*ej)
%Qi=fi*(Gij*ej-Bij*fj)-ei*(Gij*fj+Bij*ej)
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 %对PV节点的处理，注意这时不可再将h初始化为0
if i~=isb&B2(i,6)==3 %不是平衡点&是PV点
h=h+1;
for j=1:n
%公式8-75-a
%Pi=ei*(Gij*ej-Bij*fj)+fi*(Gij*fj+Bij*ej)
%Qi=fi*(Gij*ej-Bij*fj)-ei*(Gij*fj+Bij*ej)
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
%创建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
%创建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); %delPi
DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1); %delQi 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,1))-OrgS(2*h-1,1); %delPi
DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2; %delUi
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));%conj求共轭end
end
%创建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
disp('初始雅可比矩阵为：');
disp(Jacbi);
%求解修正方程，获取节点电压的不平衡量
DetaU=zeros(2*n-2,1);
DetaU=inv(Jacbi)*DetaS; %inv矩阵求逆
% 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,1)=real(B2(i,1))-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
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));
附录1 使用牛顿拉夫逊法进行潮流计算的Matlab程序代码 % 牛拉法计算潮流程序 %----------------------------------------------------------------------- % B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 % 5、支路的变比；6、支路首端处于K侧为1，1侧为0 % B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值 % 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量 % 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点；3为PV节点； %------------------------------------------------------------------------ clear all; format long; n=input('请输入节点数:nodes='); nl=input('请输入支路数:lines='); isb=input('请输入平衡母线节点号:balance='); pr=input('请输入误差精度:precision='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=input('请输入各节点参数形成的矩阵:B2='); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); %------------------------------------------------------------------ for i=1:nl %支路数 if B1(i,6)==0 %左节点处于1侧 p=B1(i,1);q=B1(i,2); else %左节点处于K侧 p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元 Y(q,p)=Y(p,q); %非对角元 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4); %对角元K侧 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4); %对角元1侧 end %求导纳矩阵 disp('导纳矩阵 Y='); disp(Y) %------------------------------------------------------------------- G=real(Y);B=imag(Y); %分解出导纳阵的实部和虚部 for i=1:n %给定各节点初始电压的实部和虚部
内蒙古科技大学 电力系统稳态分析课程设计 题目：两机五节点网络潮流计算 —牛拉法 姓名：朱润民 学号：1167130230 学院：信息工程学院 专业：电气工程及其自动化 班级：11级电气2班 指导教师：刘景霞
目录 目录 ........................................................................... - 1 - 摘要 ........................................................................ - 2 - ABSTRACT ....................................................................... - 3 - 内蒙古科技大学课程设计任务书.................................................... - 3 - 第一章电力系统潮流计算简述...................................................... - 7 - 1.1 潮流计算简介............................................................ - 7 - 1.2潮流计算的意义及发展史.................................................. - 7 - 第二章潮流计算的数学模型....................................................... - 9 - 2.1 导纳矩阵的原理及计算方法............................................... - 9 - 2.2 潮流计算的基本方程.................................................... - 11 - 2.2 电力系统节点分类...................................................... - 11 - 2.4 潮流计算的约束条件.................................................... - 11 - 第三章牛顿－拉夫逊法概述...................................................... - 16 - 3.1 牛顿-拉夫逊法基本原理................................................. - 16 - 3.2 牛顿-拉夫逊法潮流求解过程............................................. - 16 - 3.3 牛顿—拉夫逊法的程序框图.............................................. - 22 - 第四章关于电力系统潮流计算手工计算........................................... - 23 - 4.1.节点导纳矩阵........................................................... - 23 - 4.2简化雅可比矩阵......................................................... - 24 - 4.3 修正、迭代 ............................................................ - 24 - 第五章牛顿—拉夫逊法潮流具体计算............................................. - 25 - 5.1 牛顿—拉夫逊直角坐标潮流计算Matlab程序及运行结果..................... - 25 - 5.1.1 Matlab程序...................................................... - 25 - 5.1.2 Matlab程序运行结果.............................................. - 25 - 5.1.3本程序的符号说明................................................. - 49 - 总结及感想 .................................................................... - 50 - 参考文献及资料； .............................................................. - 51 -
%本程序的功能是用牛拉法进行潮流计算 %原理介绍详见鞠平著《电气工程》 %默认数据为鞠平著《电气工程》例8.4所示数据 %B1是支路参数矩阵 %第一列和第二列是节点编号。节点编号由小到大编写 %对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点编号 %第三列为支路的串列阻抗参数，含变压器支路此值为变压器短路电抗 %第四列为支路的对地导纳参数，含变压器支路此值不代入计算 %第五烈为含变压器支路的变压器的变比，变压器非标准电压比 %第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，“0”为不含有变压器 %B2为节点参数矩阵 %第一列为节点注入发电功率参数 %第二列为节点负荷功率参数 %第三列为节点电压参数 %第四列 %第五列 %第六列为节点类型参数，“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数 %X为节点号和对地参数矩阵 %第一列为节点编号 %第二列为节点对地参数 clear; clc; num=input('是否采用默认数据？(1-默认数据;2-手动输入)'); if num==1 n=4; n1=4; isb=4; pr=0.00001; B1=[1 2 0.1667i 0 0.8864 1;1 3 0.1302+0.2479i 0.0258i 1 0;1 4 0.1736+0.3306i 0.0344i 1 0;3 4 0.2603+0.4959i 0.0518i 1 0]; B2=[0 0 1 0 0 2;0 -0.5-0.3i 1 0 0 2;0.2 0 1.05 0 0 3;0 -0.15-0.1i 1.05 0 0 1]; X=[1 0;2 0.05i;3 0;4 0];
自动化07-1班段佳 07051101 function nl; %------------------------------------------------------------------------ %=================================================================== %======================牛顿——拉夫逊法============================== %===========================潮流计算================================= %=================================================================== %----------------------------------------------------------------------- % % %－－－－－－－－－－－－－－－使用说明部分－－－－－－－－－－－－－－－－－－－－－－－－－－－ display('% %本程序的功能是用牛顿——拉夫逊法进行潮流计算'); display('% %本程序要求用户按照一定的格式将电力系统的参数制成excel表格，系统运行时将从excel中加载这些参数，随后后即可进行潮流计算'); display('% %为了方便运算，用户再给系统节点进行编号时，请按照先PQ节点，再PV节点，最后平衡节点的顺序从小到大编号'); display('% %电力系统潮流计算excel格式——支路参数：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳;5、支路的变比K:1；6、支路首端处于K侧为1，1侧为0'); display('% %电力系统潮流计算excel格式——节点参数：1、节点号；2、电压大小；3、相位角；4、发电机有功；5、发电机无功；6、负载有功；7、负载无功；8、节点类型'); %=================================================================== %==============================数据准备============================== %=================================================================== % %---------------------电力系统数据加载部分----------------------------------------------- clear x=0; Branch=0;%支路参数 Note=0;%节点参数 [filename, pathname] = uigetfile('*.xls', 'please choose the excel file with your powersystem parameters ');%从外部excel导入电力系统潮流计算相关参数 try if filename ~= 0 x=xlsread([pathname,filename],'sheet1', 'A3:F3'); Branch=xlsread([pathname,filename],'sheet1', 'A5:G10');%读支路参数 Note=xlsread([pathname,filename],'sheet1', 'A15:H19');%读节点参数 end catch %进行出错处理 errmsg = lasterr; errordlg(errmsg,'Save as Error'); rethrow(lasterror); end % %---------------------支路参数初始化部分-----------------------------------------------
电力系统分析潮流计算程序设计报告题目：13节点配电网潮流计算 学院电气工程学院 专业班级 学生姓名 学号 班内序号 指导教师房大中 提交日期 2015年05月04日
目录 一、程序设计目的 (2) 二、程序设计要求 (4) 三、13节点配网潮流计算 (4) 3.1主要流程................................................................................................ 错误！未定义书签。 3.1.1第一步的前推公式如下（1-1）-（1-5）： .................................. 错误！未定义书签。 3.1.2第二步的回代公式如下（1-6）—（1-9）： ................................ 错误！未定义书签。 3.2配网前推后代潮流计算的原理 (7) 3.3配网前推后代潮流计算迭代过程 (8) 3.3计算原理 (9) 四、计算框图流程 (10) 五、确定前推回代支路次序....................................................................................... 错误！未定义书签。 六、前推回代计算输入文件 (11) 主程序： (11) 输入文件清单： (12) 计算结果： (13) 数据分析： (13) 七、配电网潮流计算的要点 (14) 八、自我总结 (14) 九、参考文献 (15) 附录一 MATLAB的简介 (15)
1.短路时有哪些异常和危害？ 2.比较直轴同步电抗，暂态电抗，次暂态电抗的大小空载电势，暂态电势Q 轴分量，次暂态电势Q轴分量的大小？ 3.三相对称短路，单相接地短路，两相短路，两相接地短路情况下正序电流如何计算比较他们的大小？ 4.小干扰法分析电力系统静态稳定？ 5.发生短路时，发电机定，转子上有那些电流分量它们如何变化？ 6.高压断路器有那些如何选择？ 7.电气主接线有那些基本要求？ 8.提高静态稳定的措施有那些你认为其中那些措施有比较好的发展前景 保定校区 专业：电自面试题：英语：电力系统的组成，电压·电流·潮流·用功英文怎么说 专业：程序设计的方法，潮流计算的高级软件有哪些（一问这蒙了），电力系统调度自动化组（远动方面的），限制短路电流的方法（发电厂电气方面的），为什么要确定功率分点，怎么确定，复杂系统怎么确定（电分方面的），逻辑元件有哪些，进而问异或的功能表述（电路方面的）。， 保定校区----电力系统及其自动化专业 1.先是中文自我介绍 2.专业英语:'s the power system .母线 3.断路器 4.继电器 5.潮流计算 3.电力系统的强行励磁的原因 4.引起电力系统暂态的原因 5.功角是什么的夹角（有两种情况:1.发电机并网的时候---发电机的空载电动势和系统的电压；2.发电机未并网的时候---转子磁极与旋转磁场的夹角） 6.发电机的特殊运行状态有哪些（我只知道迟相和进相） 7.编写程序的方法有哪些？ 语言和C++有何区别？ 和ISO是分别是什么意思？ 10.你所知道的电力电子的器件有哪些？ 11.雅克比矩阵的组成及其作用 12.并列操作的含义及其条件 1.频率下降、电压下降的影响 2.限制短路电流的方法（《发电气》－－－三种） 3.远动和电网调度的区别 4.潮流的含义 5.电力系统的组成 保定校区电力系统及其自动化（电自） 面试：1'单侧电气量的速断保护为什么不能保护线路全长 2、谈谈对电力系统的理解3、超高压是多少千伏4、交直流输电的区别和各自的优缺点,为什么,直流电怎样传输和转换5、为什么电压和无功相关(这个问题没有回答好,到现在我心里还耿耿于怀,大家要注意)6、电力电子:交流电怎样变换成直流电。7什么是运放 英语面试问题：先自我介绍，然后家乡介绍(我家是秦皇岛的所以我估计他会问家乡,所以有准备,时间用的比较长就不再问别的问题了)
%本程序的功能是用牛顿——拉夫逊法进行潮流计算 % B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 % 5、支路的变比；6、支路首端处于K侧为1，1侧为0 % B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值 % 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量 % 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点； % 3为PV节点； clear; n=input('请输入节点数:n='); nl=input('请输入支路数:nl='); isb=input('请输入平衡母线节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=input('请输入各节点参数形成的矩阵:B2='); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); % % %--------------------------------------------------- for i=1:nl %支路数 if B1(i,6)==0 %左节点处于1侧 p=B1(i,1);q=B1(i,2); else %左节点处于K侧 p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元 Y(q,p)=Y(p,q); %非对角元 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2; %对角元K侧 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; %对角元1侧 end %求导纳矩阵 disp('导纳矩阵Y='); disp(Y) %---------------------------------------------------------- G=real(Y);B=imag(Y); %分解出导纳阵的实部和虚部 for i=1:n %给定各节点初始电压的实部和虚部 e(i)=real(B2(i,3)); f(i)=imag(B2(i,3)); V(i)=B2(i,4); %PV节点电压给定模值 end for i=1:n %给定各节点注入功率 S(i)=B2(i,1)-B2(i,2); %i节点注入功率SG-SL B(i,i)=B(i,i)+B2(i,5); %i节点无功补偿量 end %=================================================================== P=real(S);Q=imag(S); %分解出各节点注入的有功和无功功率 ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; %迭代次数ICT1、a；不满足收敛要求的节点数IT2
《电力系统潮流上机》课程设计报告 院系：电气与电子工程学院 班级：电气1405 学号： 1141180505 学生姓名： 指导教师：孙英云 设计周数：两周 成绩： 日期：2017年7月5日
一、课程设计的目的与要求 培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识 二、设计正文 1．掌握计算机潮流计算的原理： a)复习电力系统分析基础中潮流的计算机算法一章，重点掌握节点分类、潮流算法介 绍 b)详细阅读牛拉法部分，掌握潮流方程（极坐标、直角坐标）的写法，掌握雅可比矩 阵的公式及排列顺序和潮流方程、变量顺序的关系，掌握迭代法收敛条件及迭代法 的基本原理 c)设计程序框图，划分功能模块、并对每个模块的输入输出量进行细化。 2．编写计算机潮流计算程序 a)学习了解IEEE标准格式数据，学习掌握C/C++读取数据的方法 b)设计计算机数据存储母线、支路数据的结构，并将所读取的数据存放于所设计的结 构当中 c)学习节点排序、节点导纳阵计算方法，编写节点导纳阵生成模块 d)编写潮流方程不平衡量计算模块 e)编写雅可比矩阵生成子模块 f)利用给定的pfMatrix类，编写修正量计算模块 g)实现潮流计算主程序，并利用IEEE标准节点数据进行校验，要求能够输出计算结 果、支路潮流等必要信息 3．思考题 1．潮流计算的方法有哪些？各有何特点？ 答：潮流计算分为简单电力网络的手算和复杂电力网络的机算两大类，其中机算又有高斯-赛德尔法、牛顿-拉夫逊法和P-Q分解法。 各方法特点如下所示： 手算求解潮流一般只用于简单的网络中，计算量大，对于多节点的网络用手算一般难以解决问题。但是通过手算可以对物理概念的理解，还可以在运用计算机计算前由手算的形式求取某些原始数据。
如图所示系统，试计算潮流分布，相关数据见《版潮流计算用户手册》P121。 #include <> #include <> float divRe(float b1,float b2,float b3,float b4) { float a1r; a1r=(b1*b3+b2*b4)/(b3*b3+b4*b4); return(a1r); } float divIm(float b1,float b2,float b3,float b4) { float a1i; a1i=(b2*b3-b1*b4)/(b3*b3+b4*b4); return(a1i); } float mulRe(float b1,float b2,float b3,float b4) { float a2r; a2r=b1*b3-b2*b4; return(a2r); } float mulIm(float b1,float b2,float b3,float b4) { float a2i; a2i=b2*b3+b1*b4; return(a2i);
} float Max(float a[],int n) {int i; float max; max= fabs(a[0]); for(i=1;i
最新文件---------------- 仅供参考--------------------已改成-----------word文本 --------------------- 方便更改 %本程序的功能是用牛顿——拉夫逊法进行潮流计算 % B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 % 5、支路的变比；6、支路首端处于K侧为1，1侧为0 % B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值 % 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量 % 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点； % 3为PV节点； clear; n=input('请输入节点数:n='); nl=input('请输入支路数:nl='); isb=input('请输入平衡母线节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=input('请输入各节点参数形成的矩阵:B2='); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); % % %--------------------------------------------------- for i=1:nl %支路数 if B1(i,6)==0 %左节点处于1侧 p=B1(i,1);q=B1(i,2); else %左节点处于K侧 p=B1(i,2);q=B1(i,1);
目录 一、任务书…………………………………….. 二、潮流计算及其要求………………………………. 三、MATLAB简介…………………………………………. 四、实验内容及设计要求…………………………………………. 五、电力系统等值电路与牛拉法程序设计流图………………………. 六、实验设计步骤及程序……………………….. 七、实验结果…………………………………………. 八、心得体会……………………………………………. 九、.参考文献……………………………….
《电力系统分析》课程设计任务书
一.潮流计算及其要求 1.潮流计算 潮流计算是研究电力系统稳态运行情况的一种基本电气计算，常规潮流计算的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。潮流计算的结果是电力系统稳定计算和故障分析的基础。通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。 在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要和基础的计算。 2潮流计算的要求 2.1.节点电压应满足 min max (1,2,)i i i U U U i n ≤≤= 从保证电能质量和供电安全的要求来看，电力系统的所有电气设备都必须运行在额定电压附近。PU 节点电压幅值必须按上述条件给定。因此，这一约束条件对PQ 节点而言。 2.2.节点的有功功率和无功功率应满足 min max min max Gi Gi Gi Gi Gi Gi P P P Q Q Q ≤≤??≤≤? PQ 节点的有功功率和无功功率，以及PU 节点的有功功率，在给定是就必须满足上述条件，因此，对平衡节点的P 和Q 以及PU 节点的Q 应按上述条件进行检验。 2.3．节点之间电压的相位差应满足 max ||||||ij i j i j θθθθθ=-<- 为了保证系统运行的稳定性，要求某些输电线路两端的电压相位不超过一定的数值。这一约束的主要意义就在于此。
1. 电力网，电力系统和动力系统的定义是什么？（p2）答: 电力系统：由发电机、发电厂、输电、变电、配电以及负荷组成的系统。 电力网：由变压器、电力线路、等变换、输送、分配电能的设备组成的部分。 动力系统：电力系统和动力部分的总和。 2. 电力系统的电气接线图和地理接线图有何区别？ （p4-5） 答：电力系统的地理接线图主要显示该系统中发电厂、变电所的地理位置，电力线路的路径以及它们相互间的连接。但难以表示各主要电机电器间的联系。 电力系统的电气接线图主要显示该系统中发电机、变压器、母线、断路器、电力线路等主要电机电器、线路之间的电气结线。但难以反映各发电厂、变电所、电力线路的相对位置。 3. 电力系统运行的特点和要求是什么？（p5）(考要求)答：特点：（1）电能与国民经济各部门联系密切。（2）电能不能大量储存。(3)生产、输送、消费电能各环节所组成的统一整体不可分割。（4）电能生产、输送、消费工况的改变十分迅速。（5）对电能质量的要求颇为严格。 要求：（1）保证可靠的持续供电。（2）保证良好的电能质量。（3）保证系统运行的经济性。 4. 电网互联的优缺点是什么？（p7） 答：可大大提高供电的可靠性，减少为防止设备事故引起供电中断而设置的备用容量；可更合理的调配用电，降低联合系统的最大负荷，提高发电设备的利用率，减少联合系统中发电设备的总容量；可更合理的利用系统中各类发电厂提高运行经济性。同时，由于个别负荷在系统中所占比重减小，其波动对系统电能质量影响也减小。联合电力系统容量很大，个别机组的开停甚至故障，对系统的影响将减小，从而可采用大容高效率的机组。 5. 我国电力网的额定电压等级有哪些？与之对应的平均额定电压是多少？系统各元件的额定电压如何确定？（p8-9）(必考看图标电压) 答：额定电压等级有（kV）：3、6、10、20、35、60、110、154、220、330、500、750、1000 平均额定电压有（kV）：3.15、6.3、10.5、37、115、230、345、525 系统各元件的额定电压如何确定：发电机母线比额定电压高5%。变压器一次接发电机比额定电压高5%，接线路为额定电压；二次接线路比额定电压高10%，变压器如果直接接负荷，则这一侧比额定电压高5%。6. 电力系统为什么不采用一个统一的电压等级，而要设置多级电压？（p8） 答：三相功率S和线电压U、线电流I之间的固定关系为。当功率一定时电压越高电流越小，导线的载流面积越小，投资越小；但电压越高对绝缘要求越高，杆塔、变压器、断路器等绝缘设备投资越大。综合考虑，对应一定的输送功率和输送距离应有一最合理的线路电压。但从设备制造角度考虑，又不应任意确定线路电压。考虑到现有的实际情况和进一步发展，我国国家标准规定了标准电压等级。 7. 导线型号LGJ-300/40中各字母和数字代表什么? （p27） 答：L表示铝，G表示钢，J表示多股导线绞合。300表示铝线额定截面积为300 ， 40表示钢线额定截面积为40 。 8. 什么是电晕现象,它和输电线路的哪个参数有关? 答：电晕指在强电磁场作用下导线周围空气的电离现象。它和输电线路的电导G有关。 9. 我国中性点接地方式有几种?为什么110kv以上电网采用中性点直接接地?110kv以下电网采用中性点不接地方式?(p10-11) 答:有不接地、直接接地、经消弧线圈接地三种接地方式。110kv以上电网采用中性点直接接地防止单相故障时某一相的电压过高。110kv以下电网采用中性点不接地方式可提高供电可靠性。 10. 架空输电线路为什么要进行换位?(p28-29) 答：为了平衡线路的三相参数。 11. 中性点不接地的电力系统发生单相接地故障时,各相对地电压有什么变化?单相接地电流的性质如何?怎样计算? （p11） 答：故障相电压等于0，非故障相电压升高倍。单相接地电流为容性。（计算见书p11.） 12. 电力系统的接线方式有哪些？各自的优缺点有哪些? （p7） 答：接线方式:有备用接线和无备用接线。 有备用接线优点：提高供电可靠性，电压质量高。缺点：不够经济。 无备用接线优点：简单，经济，运行方便。缺点：供电可靠性差。
基于Matpower2机5节点系统的潮流计算方法 引言 Matpower是基于Matlab M文件的组建包，主要用来解决电力潮流和优化潮流的问题。Matpower的设计理念是尽可能简单易懂，它可以执行电力常规潮流运算，如牛顿拉夫逊法，P-Q分解法等，也可以执行最优潮流程序。本文主要对执行常规的潮流计算进行分析。 1Matpower的简介 Matpower是基于Matlab M文件的组建包，主要用来解决电力潮流和优化潮流的问题。Matpower的设计理念是尽可能简单易懂，它可以执行电力常规潮流运算，如牛顿拉夫逊法，P-Q分解法等，也可以执行最优潮流程序。本文主要对执行常规的潮流计算进行分析。Matpower简介Matpower所用的所有数据文件均为Matlab的M文件或者MAT文件，可用来定义和返回变量BaseMV A，bus，branch，gen等[4]。其中，BaseMV A变量是一个标量，用来设置基准容量。bus 变量是一个矩阵，用来设置电网中各母线参数，其格式为bus-i，type，Pd，Qd，Gs，Bs，area，Vm，Va，baseKV，zone，Vmax，Vmin。格式中的bus-i 用来设置母线编号，范围为1~299 970；type用来设置母线类型，1为PQ节点母线，2为PV节点母线，3为平衡(参考)节点母线；Pd和Qd用来设置母线注入负荷的有功和无功功率；Gs，Bs 用来设置与母线并联电导和电纳；baseKV用来设置该母线的基准电压；Vm和Va用来设置母线电压的幅值和相位初值；Vmax 和Vmin用来设置工作时母线的最高与最低电压幅值；area和zone用来设置电网断面号和省耗分区号，一般都设置为1，设置范围分别为1~100和1~9 990[5]。 branch变量也是一个矩阵，用来设置电网中各支路参数，其格式为fbus，tbus，r，x，b，rateA，rateB，rateC，ratio，angle，status。fbus和tbus 用来设置该支路由起始节点(母线)编号和终点节点(母线)编号；r，x和b用来设置该支路的电阻、电抗和充电电纳；rateA，rateB和rateC分别用来设置该支路的长期、短期和紧急允许功率；ratio用来设置该支路的变比，如果支路元件仅仅
一、牛顿拉夫逊法进行潮流计算的MATLAB程序 %本程序的功能是用牛顿拉夫逊法进行潮流计算 %n=input('请输入节点数:n='); %nl=input('请输入支路数:nl='); %isb=input('请输入平衡母线节点号:isb='); %pr=input('请输入误差精度:pr='); %B1=input('请输入由各支路参数形成的矩阵:B1='); %B2=input('请输入各节点参数形成的矩阵:B2='); n=5; %输入节点数 nl=5; %输入支路数 isb=1; %输入平衡母线节点号 pr=0.00001; %输入误差精度 B1=[1 2 0.03i 0 1.05 0;2 3 0.08+0.3i 0.5i 1 0;2 4 0.1+0.35i 0 1 0;3 4 0.04+0.25i 0.5i 1 0;3 5 0.015i 0 1.05 1];%输入由各支路参数形成的矩阵 B2=[0 0 1.05 1.05 0 1;0 3.7+1.3i 1.05 0 0 2;0 2+1i 1.05 0 0 2;0 1.6+0.8i 1.05 0 0 2;5 0 1.05 1.05 0 3];%输入各节点参数形成的矩阵 Y=zeros(n); e=zeros(1,n);f=zeros(1,n);V=zeros(1,n); O=zeros(1,n);S1=zeros(nl); for i=1:nl %形成节点导纳矩阵 if B1(i,6)==0 %判断变压器变比是否在低压侧 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %节点互导纳 Y(q,p)=Y(p,q); %节点互导纳 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;%节点导纳=自导纳+X修正量 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; end %求导纳矩阵 disp('导纳矩阵Y='); disp(Y); G=real(Y);B=imag(Y);%分解出导纳阵的实部和虚部 for i=1:n %给定各节点初始电压的实部和虚部 e(i)=real(B2(i,3)); %取出实部 f(i)=imag(B2(i,3)); %取出虚部 V(i)=B2(i,4); %PV节点电压给定模值 end for i=1:n %给定各节点注入功率 S(i)=B2(i,1)-B2(i,2);%i节点注入功率SG-SL B(i,i)=B(i,i)+B2(i,5);%i节点无功补偿量 end P=real(S);Q=imag(S); %分解出各节点注入的有功和无功功率 ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;
辽宁工程技术大学 电力系统分析综合训练一 设计题目5节点电力网络潮流计算 指导教师 院（系、部）电气与控制工程学院 专业班级智能电网信息工程 学号1305080116 姓名 日期2016/05/04
智能电网系综合训练标准评分模板 综合设计成绩评定表 学期2015/2016第2学期姓名 专业班级 课程名称电力系统分析 设计题目 5节点电力网络潮流计算 成绩 评分项目 合格评定不合格评定 设计表现1.独立工作能力独立完成设计不能独立完成设计 2.上交设计时间按时迟交 设计说明书3.设计内容 设计思路清晰，结 构方案良好，设计 参数选择正确，条 理清楚，内容完整， 结果正确 设计思路不清晰，结 构方案不合理，关键 设计参数选择有错 误，调理清楚，内容 不完整，有明显错误4.设计书写、字 体、排版 规范、整洁、有条 理，排版很好 不规范、不整洁、无 条理，排版有问题很 大 5.封面、目录、 参考文献 完整不完整 图纸6.绘图效果满足要求很差 7.布局合理布局混乱 8.绘图工程标准符合标准不符合标准 答辩9.回答问题 回答基本正确或正 确 回答不正确总评定 评定说明： （1）不合格标准 1）设计说明书不合格否决制，即3、4两项达不到要求，不予合格； 2）9项评分标准中，有6项达不到要求，不予合格。 （2）合格标准 除设计说明书的3、4、5项必须满足要求外，其余6项，至少有4项满足要求，给予合格。 （3）请在评定栏里打“√”评定，若全部满足要求，不必分项评定，只需在总评定中打“√”即可，最后给出最终成绩，并签字。 最终成绩：评定教师签字：
电网13-1班数据（mile） 序 姓名线路1 线路2 线路3 号 1 200 100 50 2 180 125 45 3 150 250 100 4 250 80 130 5 225 100 70 6 200 80 150 7 150 150 120 8 120 150 80 9 125 175 40 10 100 150 60 11 110 200 50 12 120 180 55 13 140 160 80 14 苏小平100 200 50 15 125 180 45 16 250 150 100 17 80 250 130 18 100 225 70 19 80 200 150 20 150 150 120 21 150 120 80 22 175 125 40 23 150 100 60 24 200 110 50 25 180 120 55 说明：1. 线路1为双回线路，双回线路参数完全相同，其余线路为单回线路；2. 变压器T2为两台并联，并联运行的变压器参数相同，其他变压器为单台运行。
牛顿-拉夫逊法潮流计算的Matlab实现(2009-09-26 13:23:26)转载 标签：电力系统潮流计算牛顿-拉夫逊牛顿拉夫逊matlab newton-raphson 写了近一个礼拜的程序，藏着可惜，拿出来给需要的人 转载请注明出处 有不少朋友反映程序报错： ??? Output argument "bus" (and maybe others) not assigned during call to "E:\MAT\best\OpDF_.m>OpDF_". 是因为输入数据的格式问题输入数据中是包含bus,line的，并非只有两个矩阵，如：bus = [ 1 1.00 0.00 -0.30 -0.18 1; 2 1.00 0.00 -0.55 -0.1 3 1; 3 1.10 0.00 0.50 0.00 2; 4 1.0 5 0.00 0.00 0.00 3]; line = [ 1 2 0.10 0.40 0.0 0.01528 0; 4 2 0.08 0.40 0.0 0.01413 0; 1 4 0.1 2 0.50 0.0 0.0192 0; 3 1 0.00 0.3 0.0 0.0 -1.1]; 所有蓝色文字都为数据文件中的内容 ---------------------分---------------------割---------------------线--------------------- 【注意】 一、此程序只适用于求解节点电压以极坐标形式表示的潮流方程，没有考虑节点优化编号 二、程序在Matlab 6.5上测试通过，应该适用于目前其所有后续版本 三、程序主要通过文件方式输入输出（同时程序也返回结果向量），对输入文件格式有严格要求，具体如下： 1.输入文件可以直接在Matlab中新建m文件编写，也内容可以以文本方式编写，但最后必须改后缀名成为“.m”文件。文件名第一个字符必须是字母，后面可以跟字母、数字和下划线的任何组合，但不能和已有文件和函数冲突，不能含中文。 2.输入文件内容格式： 内容中包括"bus"（节点数据）和"line"（线路数据）两个矩阵 其中"bus"的格式为： bus = [节点编号节点电压节点相角（弧度制）有功注入无功注入节点类型]; "line"的格式为： line = [节点i 节点j 线路电阻电抗电导电纳变压器变比（普通线路为零）]; 四、程序分主程序和各子程序，需全部放入Matlab当前目录下调用，计算开始后会在当前目录下生成文件"Result.m"保存运算结果。 ---------------------分---------------------割---------------------线---------------------