目录：
一、软件需求说明书 (3)
二、概要设计说明书 (4)
1、编写潮流计算程序 (4)
2、数据的输入测试 (4)
3、运行得出结果 (4)
4、进行实验结果验证 (4)
三、详细设计说明书 (5)
1、数据导入模块 (5)
2、节点导纳矩阵模块 (5)
3、编号判断模块 (5)
4、收敛条件判定模块 (5)
5、雅可比矩阵模块 (5)
6、迭代计算模块 (5)
7、计算输出参数模块 (5)
四、程序代码 (6)
五、最测试例 (15)
1、输入结果 (15)
2、输出结果 (15)
3、结果验证 (15)
一、软件需求说明书
本次设计利用MATLAB/C++/C（使用MATLAB）编程工具编写潮流计算，实现对节点电压和功率分布的求取。
潮流方程的求解基本方法是迭代，包括牛顿-拉夫逊法，以及P-Q分解法，本次设计采用牛顿迭代法。
牛顿迭代法（Newton's method）又称为牛顿-拉夫逊方法，它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。多数方程不存在求根公式，因此求精确根非常困难，甚至不可能，从而寻找方程的近似根就显得特别重要。方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的根。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。
牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。
二、概要设计说明书
1、编写潮流计算程序
本程序主要分为七个模块：数据导入模块、节点导纳矩阵模块、编号判断模块、收敛条件判定模块、雅可比矩阵模块、迭代计算模块、计算输出参数模块。
下图为潮流迭代框图。
是
2、数据的输入测试
本次设计是将《电力系统分析》课本P88例题进行潮流计算。
3、运行得出结果
得出电压、有功功率、无功功率、角度的数据。
4、进行实验结果验证
三、详细设计说明书
1、数据导入模块：
利用Excel输入已知节点、支路数据，通过“读取”将数据导入MATLAB中。
2、节点导纳矩阵模块：
利用已知的电阻电抗及导纳的值，根据导纳的计算公式，计算出节点的自导纳及节点间的互导纳的值，按照节点编号组成导纳矩阵，利用MATLAB“real”和“imag”调用导纳矩阵中的实部和虚部，分别形成电导和电纳的矩阵。
3、编号判断模块：
当首节点在变压器左侧，设为1，位于变压器右侧，设为2，既非1也非0为不含变压器；
节点类型为PQ时，为1，节点类型为PV时，为2，节点类型为V?时，为3。
4、收敛条件判定模块：
根据节点的类型赋初值，并进行失配功率的初步计算，判断是否符合收敛条件max{ΔPi，ΔQi}<104-。如不符合，则进行后续的计算。
5、雅可比矩阵模块：
根据节点类型确定雅克比矩阵的阶数，然后根据n维非线性方程组的修正方程求出雅克比矩阵。
6、迭代计算模块：
解修正方程，并进行收敛条件判断，如不符合条件则进行下一次迭代，以一直到符合条件为止。
7、计算输出参数模块：
当满足收敛条件max{ΔPi，ΔQi}<104-
时，结束迭代计算，计算电压、有功功率、无功功率、电压角度等参数。
四、程序代码
clear %清除变量；
filename='E:\MATLAB\zhilu.xlsx';
a=xlsread(filename); %读取支路信息；
filename='E:\MATLAB\jiedian.xlsx';
c=xlsread(filename); %读取节点信息
b=zeros(a(1,7)); %定义节点导纳矩阵
G=zeros(a(1,7));
B=zeros(a(1,7));
for i1=1:(a(1,7))
if (a(i1,5)==0) %%%首节点在变压器左侧
b(a(i1,1),a(i1,1))=b(a(i1,1),a(i1,1))+1/(a(i1,3)+1j*a(i1,4));
b(a(i1,2),a(i1,2))=b(a(i1,2),a(i1,2))+1/((a(i1,3)+1j*a(i1,4))*a(i1,6)
^2);
b(a(i1,1),a(i1,2))=b(a(i1,1),a(i1,2))+(-1/((a(i1,3)+1j*a(i1,4))*a(i1,
6)));
b(a(i1,2),a(i1,1))=b(a(i1,2),a(i1,1))+(-1/((a(i1,3)+1j*a(i1,4))*a(i1, 6)));%%%进行导纳计算
G(a(i1,1),a(i1,1))=real(b(a(i1,1),a(i1,1)));
G(a(i1,2),a(i1,2))=real(b(a(i1,2),a(i1,2)));
G(a(i1,1),a(i1,2))=real(b(a(i1,1),a(i1,2)));
G(a(i1,2),a(i1,1))=real(b(a(i1,2),a(i1,1)));
B(a(i1,1),a(i1,1))=imag(b(a(i1,1),a(i1,1)));
B(a(i1,2),a(i1,2))=imag(b(a(i1,2),a(i1,2)));
B(a(i1,1),a(i1,2))=imag(b(a(i1,1),a(i1,2)));
B(a(i1,2),a(i1,1))=imag(b(a(i1,2),a(i1,1)));
elseif (a(i1,5)==1) %%%首节点在变压器右侧
b(a(i1,1),a(i1,1))=b(a(i1,1),a(i1,1))+1/((a(i1,3)+1j*a(i1,4))*a(i1,6)
^2);
b(a(i1,2),a(i1,2))=b(a(i1,2),a(i1,2))+1/(a(i1,3)+1j*a(i1,4));
b(a(i1,1),a(i1,2))=b(a(i1,1),a(i1,2))+(-1/((a(i1,3)+1j*a(i1,4))*a(i1,
6)));
b(a(i1,2),a(i1,1))=b(a(i1,2),a(i1,1))+(-1/((a(i1,3)+1j*a(i1,4))*a(i1,
6))); %%%进行导纳计算
G(a(i1,1),a(i1,1))=real(b(a(i1,1),a(i1,1)));
G(a(i1,2),a(i1,2))=real(b(a(i1,2),a(i1,2)));
G(a(i1,1),a(i1,2))=real(b(a(i1,1),a(i1,2)));
G(a(i1,2),a(i1,1))=real(b(a(i1,2),a(i1,1)));
B(a(i1,1),a(i1,1))=imag(b(a(i1,1),a(i1,1)));
B(a(i1,2),a(i1,2))=imag(b(a(i1,2),a(i1,2)));
B(a(i1,1),a(i1,2))=imag(b(a(i1,1),a(i1,2)));
B(a(i1,2),a(i1,1))=imag(b(a(i1,2),a(i1,1)));
else %%%不含变压器支路
b(a(i1,1),a(i1,1))=b(a(i1,1),a(i1,1))+1/(a(i1,3)+1j*a(i1,4))+1j*a(i1, 5);
b(a(i1,2),a(i1,2))=b(a(i1,2),a(i1,2))+1/(a(i1,3)+1j*a(i1,4))+1j*a(i1, 5);
b(a(i1,1),a(i1,2))=b(a(i1,1),a(i1,2))+(-1/((a(i1,3)+1j*a(i1,4))));
b(a(i1,2),a(i1,1))=b(a(i1,2),a(i1,1))+(-1/((a(i1,3)+1j*a(i1,4)))); %%%进行导纳计算
G(a(i1,1),a(i1,1))=real(b(a(i1,1),a(i1,1)));
G(a(i1,2),a(i1,2))=real(b(a(i1,2),a(i1,2)));
G(a(i1,1),a(i1,2))=real(b(a(i1,1),a(i1,2)));
G(a(i1,2),a(i1,1))=real(b(a(i1,2),a(i1,1)));
B(a(i1,1),a(i1,1))=imag(b(a(i1,1),a(i1,1)));
B(a(i1,2),a(i1,2))=imag(b(a(i1,2),a(i1,2)));
B(a(i1,1),a(i1,2))=imag(b(a(i1,1),a(i1,2)));
B(a(i1,2),a(i1,1))=imag(b(a(i1,2),a(i1,1)));
end
end
%%%%计数各个节点个数
PQjd=0;
PVjd=0;
VOjd=0;
for i2=1:a(1,7)
if c(i2,6)==1
PQjd=PQjd+1;
end
if c(i2,6)==2
PVjd=PVjd+1;
end
if c(i2,6)==3
VOjd=VOjd+1;
end
end
%注入功率的计算
zhuru_PQ=zeros(PQjd*2+PVjd,1);
k=1;
for i3=1:a(1,7)
if c(i3,6)==1
zhuru_PQ(k,1)=c(i3,2)-c(i3,4);
zhuru_PQ(k+a(1,7)-1,1)=0;
for i10=1:a(1,7)
zhuru_PQ(k+a(1,7)-1,1)=zhuru_PQ(k+a(1,7)-1,1)+B(i3,i10);
end
zhuru_PQ(k+a(1,7)-1,1)=zhuru_PQ(k+a(1,7)-1,1)-c(i3,5)
k=k+1;
end
if c(i3,6)==2
zhuru_PQ(k,1)=c(i3,2)-c(i3,4);
k=k+1;
end
end
n=a(1,7);
%%%%%后续矩阵变量定义
angle_u=zeros(PQjd*2+PVjd,1); %迭代角度、电压矩阵
angle_u(n:PQjd*2+PVjd,1)=1; %电压初值设为1，角度初值设为0 delta_au=zeros(PQjd*2+PVjd,1);
delta_au(n:PQjd*2+PVjd,1)=0;
%角度、电压修正量向量
%%%%%失配功率的计算
k=1;
k_=PQjd+1;
delta_PQ= zeros(PQjd*2+PVjd,1);
for i4=1:n
%%%%%%%%%%%%%%%%%%%PQ节点的计算
if c(i4,6)==1
delta_PQ(k,1)=zhuru_PQ(k,1);
delta_PQ(k+n-1,1)=zhuru_PQ(k+n-1,1);
i6=1;
for i5=1:n
if c(i5,6)==1
delta_PQ(k,1)=delta_PQ(k,1)-angle_u(n-1+k)*angle_u(n-1+i6)*(G(i4,i5)*
cos(angle_u(k,1)-angle_u(i6,1))+B(i4,i5)*sin(angle_u(k,1)-angle_u(i6, 1)));
delta_PQ(k+n-1,1)=delta_PQ(k+n-1,1)-angle_u(n-1+k)*angle_u(n-1+i6)*(G (i4,i5)*sin(angle_u(k,1)-angle_u(i6,1))-B(i4,i5)*cos(angle_u(k,1)-ang le_u(i6,1)));
i6=i6+1;
end
if c(i5,6)==2
delta_PQ(k,1)=delta_PQ(k,1)-angle_u(n-1+k)*c(i5,1)*(G(i4,i5)*cos(angl e_u(k,1)-angle_u(i6,1))+B(i4,i5)*sin(angle_u(k,1)-angle_u(i6,1)));
delta_PQ(k+n-1,1)=delta_PQ(k+n-1,1)-angle_u(n-1+k)*c(i5,1)*(G(i4,i5)* sin(angle_u(k,1)-angle_u(i6,1))-B(i4,i5)*cos(angle_u(k,1)-angle_u(i6, 1)));
i6=i6+1;
end
if c(i5,6)==3
delta_PQ(k,1)=delta_PQ(k,1)-angle_u(n-1+k)*c(i5,1)*(G(i4,i5)*cos(angl e_u(k,1))+B(i4,i5)*sin(angle_u(k,1)));
delta_PQ(k+n-1,1)=delta_PQ(k+n-1,1)-angle_u(n-1+k)*c(i5,1)*(G(i4,i5)* sin(angle_u(k,1))-B(i4,i5)*cos(angle_u(k,1)));
end
end
k=k+1;
end
%%%%%%%%%%%%%%%%%%%%%PV节点的计算
if c(i4,6)==2
delta_PQ(k_,1)=zhuru_PQ(k_,1);
i6=1;
for i5=1:n
if c(i5,6)==1
delta_PQ(k_,1)=delta_PQ(k_,1)-c(i4,1)*angle_u(n-1+i6)*(G(i4,i5)*cos(a ngle_u(k_,1)-angle_u(i6,1))+B(i4,i5)*sin(angle_u(k_,1)-angle_u(i6,1)) );
i6=i6+1;
end
if c(i5,6)==2
delta_PQ(k_,1)=delta_PQ(k_,1)-c(i4,1)*c(i5,1)*(G(i4,i5)*cos(angle_u(k _,1)-angle_u(i6,1))+B(i4,i5)*sin(angle_u(k_,1)-angle_u(i6,1)));
i6=i6+1;
end
if c(i5,6)==3
delta_PQ(k_,1)=delta_PQ(k_,1)-c(i4,1)*c(i5,1)*(G(i4,i5)*cos(angle_u(k _,1))+B(i4,i5)*sin(angle_u(k_,1)));
end
end
k_=k_+1;
end
end
zhongjian_delta_PQ=abs(delta_PQ); %%%%取绝对值
jingdu=max(zhongjian_delta_PQ); %%%%为迭代精度的判断做赋值准备
while (jingdu>0.0001)
%%%%%%%雅各比矩阵的形成
J=zeros(2*PQjd+PVjd); %%%%矩阵定义
U=ones(n,1); %%%%为方便运算将电压取出放入一个新的矩阵
for i=1:PQjd
U(i,1)=angle_u(i+n-1,1);
end
angle=zeros(n,1); %%%%为方便运算将角度取出放入一个新的矩阵
for i=1:n-1
angle(i,1)=angle_u(i,1);
end
%%%%%H%%%%% %%%%%H部分生成
for i1=1:n-1
for i2=1:n-1
if i1~=i2
J(i1,i2)=-U(i1,1)*U(i2,1)*(G(i1,i2)*sin(angle(i1,1)-angle(i2,1))-B(i1 ,i2)*cos(angle(i1,1)-angle(i2,1)));
end
if i1==i2
for i3=1:n
if i3~=i1
J(i1,i2)=J(i1,i2)+U(i1,1)*U(i3,1)*(G(i1,i3)*sin(angle(i1,1)-angle(i3, 1))-B(i1,i3)*cos(angle(i1,1)-angle(i3,1)));
end
end
end
end
end
%%%%%N%%%%% %%%%%N部分生成
for i1=1:n-1
for i2=1:PQjd
if i1~=i2
J(i1,i2+n-1)=-U(i1,1)*(G(i1,i2)*cos(angle(i1,1)-angle(i2,1))-B(i1,i2) *sin(angle(i1,1)-angle(i2,1)));
end
if i1==i2
J(i1,i2+n-1)=(-2)*U(i1,1)*U(i1,1)*G(i1,i1);
for i3=1:n
if i3~=i1
J(i1,i2+n-1)=J(i1,i2+n-1)-U(i3,1)*(G(i1,i3)*cos(angle(i1,1)-angle(i3, 1))+B(i1,i3)*sin(angle(i1,1)-angle(i3,1)));
end
end
end
end
end
%%%%%K%%%%% %%%%%K部分生成
for i1=1:PQjd
for i2=1:n-1
if i1~=i2
J(i1+n-1,i2)=U(i1,1)*U(i2,1)*(G(i1,i2)*cos(angle(i1,1)-angle(i2,1))-B (i1,i2)*sin(angle(i1,1)-angle(i2,1)));
end
if i1==i2
for i3=1:n
if i3~=i1
J(i1+n-1,i2)=J(i1+n-1,i2)-U(i1,1)*U(i3,1)*(G(i1,i3)*cos(angle(i1,1)-a
ngle(i3,1))+B(i1,i3)*sin(angle(i1,1)-angle(i3,1)));
end
end
end
end
end
%%%%%L%%%%% %%%%%L部分生成
for i1=1:PQjd
for i2=1:PQjd
if i1~=i2
J(i1+n-1,i2+n-1)=-U(i1,1)*(G(i1,i2)*sin(angle(i1,1)-angle(i2,1))-B(i1 ,i2)*cos(angle(i1,1)-angle(i2,1)));
end
if i1==i2
J(i1+n-1,i2+n-1)=2*U(i1,1)*U(i1,1)*B(i1,i1);
for i3=1:n
if i3~=i1
J(i1+n-1,i2+n-1)=J(i1+n-1,i2+n-1)-U(i3,1)*(G(i1,i3)*sin(angle(i1,1)-a ngle(i3,1))-B(i1,i3)*cos(angle(i1,1)-angle(i3,1)));
end
end
end
end
end
delta_au=-inv(J)*delta_PQ; %%%解修正方程
angle_u=angle_u+delta_au; %%%形成新的初始值%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%再次计算失配功率
k=1;
k_=PQjd+1;
delta_PQ= zeros(PQjd*2+PVjd,1);
for i4=1:n
%%%%%%%%%%%%%%%%%%%
if c(i4,6)==1
delta_PQ(k,1)=zhuru_PQ(k,1);
delta_PQ(k+n-1,1)=zhuru_PQ(k+n-1,1);
i6=1;
for i5=1:n
if c(i5,6)==1
delta_PQ(k,1)=delta_PQ(k,1)-angle_u(n-1+k)*angle_u(n-1+i6)*(G(i4,i5)* cos(angle_u(k,1)-angle_u(i6,1))+B(i4,i5)*sin(angle_u(k,1)-angle_u(i6, 1)));
delta_PQ(k+n-1,1)=delta_PQ(k+n-1,1)-angle_u(n-1+k)*angle_u(n-1+i6)*(G (i4,i5)*sin(angle_u(k,1)-angle_u(i6,1))-B(i4,i5)*cos(angle_u(k,1)-ang le_u(i6,1)));
i6=i6+1;
end
if c(i5,6)==2
delta_PQ(k,1)=delta_PQ(k,1)-angle_u(n-1+k)*c(i5,1)*(G(i4,i5)*cos(angl e_u(k,1)-angle_u(i6,1))+B(i4,i5)*sin(angle_u(k,1)-angle_u(i6,1)));
delta_PQ(k+n-1,1)=delta_PQ(k+n-1,1)-angle_u(n-1+k)*c(i5,1)*(G(i4,i5)* sin(angle_u(k,1)-angle_u(i6,1))-B(i4,i5)*cos(angle_u(k,1)-angle_u(i6, 1)));
i6=i6+1;
end
if c(i5,6)==3
delta_PQ(k,1)=delta_PQ(k,1)-angle_u(n-1+k)*c(i5,1)*(G(i4,i5)*cos(angl e_u(k,1))+B(i4,i5)*sin(angle_u(k,1)));
delta_PQ(k+n-1,1)=delta_PQ(k+n-1,1)-angle_u(n-1+k)*c(i5,1)*(G(i4,i5)* sin(angle_u(k,1))-B(i4,i5)*cos(angle_u(k,1)));
end
end
k=k+1;
end
%%%%%%%%%%%%%%%%%%%%%
if c(i4,6)==2
delta_PQ(k_,1)=zhuru_PQ(k_,1);
i6=1;
for i5=1:n
if c(i5,6)==1
delta_PQ(k_,1)=delta_PQ(k_,1)-c(i4,1)*angle_u(n-1+i6)*(G(i4,i5)*cos(a ngle_u(k_,1)-angle_u(i6,1))+B(i4,i5)*sin(angle_u(k_,1)-angle_u(i6,1))
);
i6=i6+1;
end
if c(i5,6)==2
delta_PQ(k_,1)=delta_PQ(k_,1)-c(i4,1)*c(i5,1)*(G(i4,i5)*cos(angle_u(k _,1)-angle_u(i6,1))+B(i4,i5)*sin(angle_u(k_,1)-angle_u(i6,1)));
i6=i6+1;
end
if c(i5,6)==3
delta_PQ(k_,1)=delta_PQ(k_,1)-c(i4,1)*c(i5,1)*(G(i4,i5)*cos(angle_u(k _,1))+B(i4,i5)*sin(angle_u(k_,1)));
end
end
k_=k_+1;
end
end
%%%%再次进行精度判断前的准备
m=0;
zhongjian_delta_PQ=delta_PQ;
zhongjian_delta_PQ=abs(delta_PQ);
jingdu=max(zhongjian_delta_PQ);
m=m+1;
end
%%%%%%迭代完毕
%%%%%%电压与角度换算结果输出
U=ones(n,1);
for i=1:PQjd
U(i,1)=angle_u(i+n-1,1);
end
angle=zeros(n,1);
for i=1:n-1
angle(i,1)=angle_u(i,1);
end
angle_jiaodu=angle/(2*pi)*360;
%%%%%%%%%计算各个节点功率
P=zeros(n,1);
Q=zeros(n,1);
for i1=1:n
for i2=1:n
P(i1,1)=P(i1,1)+U(i1,1)*U(i2,1)*(G(i1,i2)*cos(angle(i1,1)-angle(i2,1) )+B(i1,i2)*sin(angle(i1,1)-angle(i2,1)));
Q(i1,1)=Q(i1,1)+U(i1,1)*U(i2,1)*(G(i1,i2)*sin(angle(i1,1)-angle(i2,1) )-B(i1,i2)*cos(angle(i1,1)-angle(i2,1)));
end
end
五、最测试例
本次设计试例为《电力系统分析》课本P88例题3-3。
1、输入数据
支路数据如下：
节电数据如下：
2、输出结果
J矩阵
k 电压U 角度?有功功率P 无功功率Q
3、结果验证
输出结果与课本P90结果一致，潮流计算程序正确。
电气2013级卓班电力电子技术与电力系统分析 课程实训报告 专业：电气工程及其自动化 班级： 姓名： 学号： 指导教师：
兰州交通大学自动化与电气工程学院 2016 年 1 月日
电力电子技术与电力系统分析课程实训报告 1 电力电子技术实训报告 1.1 实训题目 1.1.1电力电子技术实训题目一 一．单相半波整流 参考电力电子技术指导书中实验三负载，建立MATLAB/Simulink环境下三相半波整流电路和三相半波有源逆变电路的仿真模型。仿真参数设置如下： (1)交流电压源的参数设置和以前实验相关的参数一样。 (2)晶闸管的参数设置如下： R=0.001Ω，L =0H，V f=0.8V，R s=500Ω，C s=250e-9F on (3)负载的参数设置 RLC串联环节中的R对应R d，L对应L d，其负载根据类型不同做不同的调整。 (4)完成以下任务： ①仿真绘出电阻性负载（RLC串联负载环节中的R d= Ω，电感L d=0，C=inf，反电动势为0）下α=30°，60°，90°，120°，150°时整流电压U d，负载电流L 和晶闸管两端电压U vt1的波形。 d ②仿真绘出阻感性负载下（负载R d=Ω，电感L d为，反电动势E=0）α=30°，60°，90°，120°，150°时整流电压U d，负载电流L d和晶闸管两端电压U vt1的波形。 ③仿真绘出阻感性反电动势负载下α=90°，120°，150°时整流电压U d，负载电流L d和晶闸管两端电压U vt1的波形，注意反电动势E的极性。 (5)结合仿真结果回答以下问题： ①该三项半波可控整流电路在β=60°，90°时输出的电压有何差异?
三、计算分析题 1、图1.5.1所示电路，已知U =3V ，求R 。（2k Ω） 2、图1.5.2所示电路，已知U S ＝3V ，I S ＝2A ，求U AB 和I 。（1V 、5A ） 3、电路如图1.5.5所示，求10V 电压源发出的 功率。 （－35W ） 4、分别计算S 打开与闭合时图1.5.6电路中A 、B 两点的电位。（S 打开：A －10.5V,B －7.5V S 闭合：A 0V ，B 1.6V ） 5、试求图1.5.7所示电路的入端电阻R AB 。（150Ω） U － 图1.5.1 1Ω 图1.5.2 6V 图1.5.5 B －图1.5.6 Ω 图1.5.7
6、试求图2.4.1所示电路的电压U 。 7、已知图2.5.1电路中电压U =4.5V ，试应用已经学过的电路求解法求电阻R 。 （18Ω） 8、求解图2.5.2所示电路的戴维南等效电路。 （U ab =0V ，R 0=8.8Ω） 9、列出图2.5.4所示电路的结点电压方程。 解：画出图2.5.4等效电路图如下： 图2.5.1 9V 图2.5.2 2A Ω U 图2.4.1题电路
对结点A 对结点B 10、应用等效变换求图示电路中的I的值。（10分） 解：等效电路如下： 11、应用等效变换求图示电路中的I的值。
12、应用戴维南定理求解图示电路中的电流I 13、如下图所示，RL等于何值时，能得到最大传输功率P0max？并计算P0max 。
16、图示电路中，开关闭合之前电路已处于稳定状态，已知R1=R2=2Ω 解开关闭合后电感电流iL的全响应表达式。 17、图示电路中，t=0时开关闭合，闭合之前电路已处于稳定状态，请用三要素法求解开关闭合后电容电压uc的全响应表达式。
目录 摘要 (1) 1.任务及题目要求 (2) 2.计算原理 (3) 2.1牛顿—拉夫逊法简介 (3) 2.2牛顿—拉夫逊法的几何意义 (7) 3计算步骤 (7) 4.结果分析 (9) 小结 (11) 参考文献 (12) 附录：源程序 (13) 本科生课程设计成绩评定表 (32)
摘要 电力系统的出现，使高效，无污染，使用方便，易于调控的电能得到广泛应用，推动了社会生产各个领域的变化，开创了电力时代，发生率第二次技术革命。电力系统的规模和技术水准已经成为一个国家经济发展水平的标志之一。 电力系统稳态分析包括潮流计算和静态安全分析。电力系统潮流计算是电力系统最基本的计算，也是最重要的计算。所谓潮流计算，就是已知电网的接线方式与参数及运行条件，计算电力系统稳态运行各母线电压、个支路电流与功率及网损。对于正在运行的电力系统，通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。 在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。 关键词：电力系统潮流计算牛顿-拉夫逊法
《电力系统潮流上机》课程设计报告 院系：电气工程学院 班级：电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 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 四、设计内容
电力系统的潮流计算 西安交通大学自动化学院 2012.10 3．1 电网结构如图3—11所示，其额定电压为10KV 。已知各节点的负荷功率及参数： MVA j S )2.03.0(2 +=， MVA j S )3.05.0(3+=， MVA j S )15.02.0(4+= Ω+=)4.22.1(12j Z ，Ω+=)0.20.1(23j Z ，Ω+=)0.35.1(24j Z 试求电压和功率分布。 解：（1）先假设各节点电压均为额定电压，求线路始端功率。 0068.00034.0)21(103.05.0)(2 2223232232323j j jX R V Q P S N +=++=++=?0019.00009.0)35.1(10 15.02.0)(2 2 224242242424j j jX R V Q P S N +=++=++=?
则： 3068.05034.023323j S S S +=?+= 1519.02009.024424j S S S +=?+= 6587.00043.122423' 12 j S S S S +=++= 又 0346 .00173.0)4.22.1(106587.00043.1)(2 2 212122'12'1212j j jX R V Q P S N +=++=++=? 故： 6933.00216.112'1212 j S S S +=?+= （2） 再用已知的线路始端电压kV V 5.101 =及上述求得的线路始端功率 12 S ，求出线 路 各 点 电 压 。
kV V X Q R P V 2752.05 .104.26933.02.10216.1)(11212121212=?+?=+=? kV V V V 2248.101212=?-≈ kV V V V kV V X Q R P V 1508.100740.0) (24242 2424242424=?-≈?=+=? kV V V V kV V X Q R P V 1156.101092.0) (23232 2323232323=?-≈?=+=? （3）根据上述求得的线路各点电压，重新计算各线路的功率损耗和线路始端功率。 0066.00033.0)21(12.103.05.02 2 223j j S +=++=? 0018.00009.0)35.1(15 .1015.02.02 2 224j j S +=++=? 故 3066.05033.023323j S S S +=?+= 1518.02009.024424j S S S +=?+= 则 6584.00042.122423' 12 j S S S S +=++= 又 0331.00166.0)4.22.1(22 .106584.00042.12 2 212j j S +=++=? 从而可得线路始端功率 6915.00208.112 j S +=
题目：潮流计算与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
吉林大学网络教育学院 2019-2020学年第二学期期末考试《电力系统分析》大作业 学生姓名专业 层次年级学号 学习中心成绩 年月日
作业要求：大作业要求学生手写完成，提供手写文档的清晰扫描图片，并将图片添加到word文档内，最终wod文档上传平台，不允许学生提交其他格式文件（如JPG，RAR等非word文档格式），如有雷同、抄袭成绩按不及格处理。 一计算题 (共9题，总分值90分 ) 1. 有一台型10kv网络供电的降压变压器，铭牌给出的试验数据为：。 试求（1）计算折算到一次（二次）侧的变压器参数，并作其Г型Π型等值电路 变压器不含励磁之路时的Π型等值电路。（10 分）
2. 降压变压器及等效电路示于图5-7a、b。折算至一次侧的阻抗为Ω。已知在最大负荷和最小负荷时通过变压器的功率分别为，一次侧的电压分别为=110KV和113KV。要求二次侧母线的变化不超过6.0—6.6KV的范围，试选择分接头。 图5-19 习题5-8a 5-8b （10 分）
3. 简单电力系统如图7-52习题7-7所示，已知元件参数如下：发电机：，=0.16， =0.19；变压器：，=10。5，k点分别发生单相接地、两相短路、两相接地和三相短路时，试计算短路点短路电流的有名值，并进行比较分析。 图7-52 习题7-7（10 分）
4.已知一200km长的输电线，R=0.1Ω/km，L=2.0mH/km，C=0.01μF/km，系统频率为50Hz。使用（1）短线路，（2）中程线路，（3）长线路模型求其π形等效电路。（10 分） 解： （1）短线路一字型等值电路参数： （2）中程线路∏形等值电路参数（不需修正）： （3）长线路：
2、停电有可能导致人员伤亡或主要生产设备损坏的用户的用电设备属于（ ）。 A 、一级负荷； B 、二级负荷； C 、三级负荷； D 、特级负荷。 4、衡量电能质量的技术指标是（ ）。 A 、电压偏移、频率偏移、网损率； B 、电压偏移、频率偏移、电压畸变率； C 、厂用电率、燃料消耗率、网损率； D 、厂用电率、网损率、电压畸变率 5、用于电能远距离输送的线路称为（ ）。 A 、配电线路； B 、直配线路； C 、输电线路； D 、输配电线路。 7、衡量电力系统运行经济性的主要指标是（ ）。 A 、燃料消耗率、厂用电率、网损率； B 、燃料消耗率、建设投资、网损率； C 、网损率、建设投资、电压畸变率； D 、网损率、占地面积、建设投资。 8、关于联合电力系统，下述说法中错误的是（ ）。 A 、联合电力系统可以更好地合理利用能源； B 、在满足负荷要求的情况下，联合电力系统的装机容量可以减少； C 、联合电力系统可以提高供电可靠性和电能质量； D 、联合电力系统不利于装设效率较高的大容量机组。 9、我国目前电力系统的最高电压等级是（ ）。 A 、交流500kv ，直流kv 500±； B 、交流750kv ，直流kv 500±； C 、交流500kv ，直流kv 800±；； D 、交流1000kv ，直流kv 800±。 10、用于连接220kv 和110kv 两个电压等级的降压变压器，其两侧绕组的额定电压应为（ ）。 A 、220kv 、110kv ； B 、220kv 、115kv ； C 、242Kv 、121Kv ； D 、220kv 、121kv 。 11、对于一级负荷比例比较大的电力用户，应采用的电力系统接线方式为（ ）。 A 、单电源双回路放射式； B 、双电源供电方式； C 、单回路放射式接线； D 、单回路放射式或单电源双回路放射式。 12、关于单电源环形供电网络，下述说法中正确的是（ ）。 A 、供电可靠性差、正常运行方式下电压质量好； B 、供电可靠性高、正常运行及线路检修（开环运行）情况下都有好的电压质量； C 、供电可靠性高、正常运行情况下具有较好的电压质量，但在线路检修时可能出现电压质量较差的情况； D 、供电可靠性高，但电压质量较差。 13、关于各种电压等级在输配电网络中的应用，下述说法中错误的是（ ）。 A 、交流500kv 通常用于区域电力系统的输电网络； B 、交流220kv 通常用于地方电力系统的输电网络； C 、交流35kv 及以下电压等级通常用于配电网络； D 、除10kv 电压等级用于配电网络外，10kv 以上的电压等级都只能用于输电网络。 14、110kv 及以上电力系统应采用的中性点运行方式为（ ）。 A 、直接接地； B 、不接地； C 、经消弧线圈接地； D 、不接地或经消弧线圈接地。 16、110kv 及以上电力系统中，架空输电线路全线架设避雷线的目的是（ ）。
电力系统分析潮流计算报告
目录 一．配电网概述 (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 配电网运行的特点及要求 配电系统相对于输电系统来说，由于电压等级低、供电范围小，但与用户直接相连，是供电部门对用户服务的窗口，因而决定了配电网运行有如下特点和基本要求：
电力系统分析大作业
一、设计题目 本次设计题目选自课本第五章例5-8,美国西部联合电网WSCC系统的简化三机九节点系统，例题中已经给出了潮流结果，计算结果可以与之对照。取ε=0.00001 。
二、计算步骤 第一步，为了方便编程，修改节点的序号，将平衡节点放在最后。如下图： 第二步，这样得出的系统参数如下表所示： 第三步，形成节点导纳矩阵。 9 2 1 3 2 7 4 5 6 8 3
第四步，设定初值： ο 01)0(6)0(5)0(4)0(3)0(2)0(1∠======??????U U U U U U ； 0)0(8)0(7==Q Q ，0)0(8)0(7==θθ。 第五步，计算失配功率 )0(1P ?=0，)0(2P ?=-1.25，)0(3P ?=-0.9，) 0(4P ?=0，)0(5P ?=-1，)0(6P ?=0，)0(7P ?=1.63， )0(8P ?=0.85； )0(1Q ?=0.8614，)0(2Q ?=-0.2590，)0(3Q ?=-0.0420，) 0(4Q ?=0.6275，)0(5Q ?=-0.1710， )0(6Q ?=0.7101。 显然，5108614.0|},max {|-=>=??εi i Q P 。 第六步，形成雅克比矩阵（阶数为14×14） 第七步，解修正方程，得到： =?)0(1θ-0.0371，=?)0(2θ-0.0668，=?)0(3θ-0.0628，=?)0(4θ0.0732，=?)0(5θ0.0191，=?)0(6θ0.0422，=?)0(7θ0.1726，=?)0(8θ0.0908； =?)0(1U 0.0334，=?)0(2U 0.0084，=?)0(3U 0.0223，=?)0(4U 0.0372，=?)0(5U 0.0266，
第一章 电力系统的基本概念 1-1 什么叫电力系统、电力网及动力系统 1-2 电力线、发电机、变压器和用电设备的额定电压是如何确定的 1-3 我国电网的电压等级有哪些 1-4 标出图1-4电力系统中各元件的额定电压。 1-5 请回答如图1－5所示电力系统中的二个问题： ⑴ 发电机G 、变压器1T 2T 3T 4T 、三相电动机D 、单相电灯L 等各元件的额定电压。 ⑵ 当变压器1T 在+%抽头处工作，2T 在主抽头处工作，3T 在%抽头处工作时，求这些变压器的实际变比。 1-6 图1-6中已标明各级电网的电压等级。试标出图中发电机和电动机的额定电压及变压器的额定变比。 1-7 电力系统结线如图1-7所示，电网各级电压示于图中。试求： ⑴发电机G 和变压器1T 、2T 、3T 高低压侧的额定电压。 习题1-4图
⑵设变压器1T 工作于+%抽头， 2T 工作于主抽头，3T 工作于-5%抽头，求这些变压器的实际变比。 1-8 比较两种接地方式的优缺点，分析其适用范围。 1-9 什么叫三相系统中性点位移它在什么情况下发生中性点不接地系统发生单相接地时，非故障相电压为什么增加3倍 1-10 若在变压器中性点经消弧线圈接地，消弧线圈的作用是什么 第二章 电力系统各元件的参数及等值网络 2-1 一条110kV 、80km 的单回输电线路，导线型号为LGJ —150，水平排列，其线间距离为4m ，求此输电线路在40℃时的参数，并画出等值电路。 2-2 三相双绕组变压器的型号为SSPL —63000/220，额定容量为63000kVA ，额定电压为242/，短路损耗404=k P kW ，短路电压45.14%=k U ，空载损耗93=o P kW ，空载电流 41.2%=o I 。求该变压器归算到高压侧的参数，并作出等值电路。 2-3 已知电力网如图2-3所示： 各元件参数如下： 变压器：1T ：S =400MVA ，12%=k U ， 242/ kV 2T ：S =400MVA ，12%=k U ， 220/121 kV 线路：2001=l km, /4.01Ω=x km (每回路) 602=l km, /4.01Ω=x km 115kV T 1 T 2 l 1 l 2 习题2-3图
电力系统分析练习题 一、单项选择题 1.对电力系统的基本要求就是(A) A.保证供电可靠性、保证良好的电能质量、保证系统运行的经济性 B.保证供电可靠性、保证良好的电能质量、保证继电保护动作的选择性 C.保证供电可靠性与良好的电能质量 D.保证供电可靠性与系统运行的经济性 2.下图所示的电力系统中,变压器T1的额定变比应为(B) A.242/121/11kV B.220/121/11kV C.220/110/10kV D.242/121/10kV 3.连接220kV 电力系统与110kV 电力系统的降压变压器,其额定变比应为(D) A.220/110kV B.220/115、5kV C.242/121kV D.220/121kV 4.我国110kV 以上电力系统中性点通常采用的运行方式就是(B) A 、不接地 B 、直接接地 C 、经消弧线圈接地 D 、经电容器接地 5.电力系统中性点经消弧线圈接地时,应采用的补偿方式为(A) A.过补偿 B.欠补偿 C.全补偿 D.全补偿或欠补偿 6.采用同一型号导线的三相架空输电线路,相间距离增大时,其电容(B) A.增大 B.减小 C.不变 D.无法确定 7.架空输电线路采用分裂导线的目的就是(A) A 、减小线路电抗 B 、减小线路电阻 C 、减小线路电容 D 、增大线路电抗 8.三相架空输电线路导线全换位的目的就是(A) A.减小三相参数的不平衡 B.减小线路电抗 C.减小线路电容 D.提高输电线路的输送容量 9.变压器参数B T 由试验数据 确定。(2章)(B) A.%K U B.%0I C.0P ? D.K P ? 10.变压器的电抗参数 T X ,由实验数据 确定。(A) A 、%K U B 、%0I C 、0P ? D 、K P ? 11.变压器的电纳参数T G ,由实验数据 确定。(C) A.%K U B.%0I C.0P ? D.K P ? 12.变压器的电阻参数T R ,由实验数据 确定。(D) A.%K U B.%0I C.0P ? D.K P ? 13.如果三相功率基准值为b S 、线电压基准值为b U ,则阻抗基准值为(D) A.b b U S / B. b b U S 3/ C.2/b b U S D.b b S U /2 14.输电线路运行时,其始端电压电压与末端电压的大小关系就是(C) A.始端电压一定高于末端电压 B. 始端电压始终等于末端电压 C.在不计电压降落横分量影响情况下,PR+QX>0时,始端电压高于末端电压 D.在不计电压降落横分量影响情况下,PR+QX>0时,始端电压低于末端电压 15.输电线路的电压降落就是指(A) A 、线路始端电压与末端电压的相量差 B 、线路始端电压与末端电压的数值差 C 、线路始端电压与额定电压的数值差 D 、线路末端电压与额定电压的数值差
课程设计 课程名称：电力系统分析 设计题目：基于Matlab计算程序地电力系统运行分析学院：电力工程学院 专业：电气工程自动化 年级： 学生姓名： 指导教师： 日期： 教务处制
目录 前言 (1) 第一章参数计算 (2) 一、目标电网接线图 (2) 二、电网模型地建立 (3) 第二章潮流计算 (6) 一．系统参数地设置 (6) 二．程序地调试 (7) 三、对运行结果地分析 (13) 第三章短路故障地分析计算 (15) 一、三相短路 (15) 二、不对称短路 (16) 三、由上面表对运行结果地分析及在短路中地一些问题 (21) 心得体会 (26) 参考文献 (27)
前言 电力系统潮流计算是电力系统分析中地一种最基本地计算，是对复杂电力系统正常和故障条件下稳态运行状态地计算.潮流计算地目标是求取电力系统在给定运行状态地计算.即节点电压和功率分布，用以检查系统各元件是否过负荷.各点电压是否满足要求，功率地分布和分配是否合理以及功率损耗等.对现有电力系统地运行和扩建，对新地电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础.潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算地模型和方法有直接影响. 在电力系统中可能发生地各种故障中，危害最大且发生概率较高地首推短路故障.产生短路故障地主要原因是电力设备绝缘损坏.短路故障分为三相短路、两相短路、单相接地短路及两相接地短路.其中三相短路时三相电流仍然对称，其余三类短路统成为不对称短路.短路故障大多数发生在架空输电线路.电力系统设计与运行时，要采取适当地措施降低短路故障地发生概率.短路计算可以为设备地选择提供原始数据.
电力系统分析课程设计 设计题目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. 应用该平台可以绘制各种电网图形，包括单线图、地理位置接线图、厂站主接线图等。 ●所有图形独立于各种分析计算，并为各计算模块所共享； ●可在图形上进行各种计算操作，并在图上显示各种计算结果； ●同一系统可对应多套单线图，多层子图嵌套； ●单线图上可细化到厂站主接线结构；
《电力系统潮流上机》课程设计报告院系：电气工程学院 班级：电088班 学号： 学生姓名：刘东昇 指导教师：张新松 设计周数：两周 日期：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 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也
电力系统分析计算题 1.假设电压U 1固定不变，试就图示系统分析为什么投入电容器C 后可以降低线损和提高电压U 2。 2. 试就图示系统分析为什么对中枢点U 1进行逆调压可以对负荷点电压U 2进行控制的原理。 3. 试就图示系统分析（a ）、（b ）两种情况下线路的电能损耗ΔA ，你的结论是什么？ 4.某负荷由发电厂经电压为110kV 的输电线路供电。 Ω, X=41Ω，C Q 2 1 =1.7Mvar ， 发 线路的参数为：R=17 电厂高压母线电压U 1=116kV, 线路末端负荷为20+j10 MVA,求输电线路的功率损耗和末端电压U 2（计及电压降落的横分量)。(12分) 5、某降压变电所装有两台并联工作的有载调压变压器，电压为110±5×2.5/11 kV ，容量为31.5MVA, 已知最大负荷时:高压母线电压为103 kV,两台变压器并列运行时的电压损耗为5.849 kV,最小负荷时: 高压母线电压为108.5 kV, 两台变压器并列运行时的电压损耗为2.631 kV 。变电所低压母线要求逆调压，试选择有载调压变压器分接头。(12分) 6.一双电源电力系统如下图所示，如在f 点发生b 、 P+jQ P+jQ S ~ （h ） （a ） （b ） 20+j10 MVA U 2 2 T 1 1 MVA X X j E d G 3013.013.01"2., ,. 1. === U K %=10.5 60MV A U K %=10.5 31.5MV A 80km X 1=0.4Ω/km X 0=3.5 X 1 MVA X X j E d G 60125.0125.01".2, ,. 2. ===
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
电力系统分析课程报告姓名 ******* 学院自动化与电气工程学院 专业控制科学与工程 班级 ******* 指导老师 ******* 二〇一六年五月十三
一、同步发电机三相短路仿真 1、仿真模型的建立 选取三相同步发电机模型，以三相视图表示。励磁电压和原动机输入转矩Ef 与Tm均为定常值，且发电机空载。当运行至时，发电机发生三相短路故障。同步发电机三相短路实验仿真模型如图1所示。 图1 同步发电机三相短路实验仿真模型 2、发电机参数对仿真结果的影响及分析 衰减时间常数Ta对于直流分量的影响 三相短路电流的直流分量大小不等，但衰减规律相同，均按指数规律衰减，衰减时间常数为Ta，由定子回路的电阻和等值电感决定（大约）。pscad同步发电机模型衰减时间常数Ta对应位置如图3所示（当前Ta=）。 图3 同步发电机模型参数Ta对应位置
1）Ta=时，直流分量的衰减过程（以励磁电流作为分析）如图4所示。 图4 Ta=发生短路If波形 2）Ta=时，直流分量的衰减过程（以励磁电流作为分析）如图5所示。 图5 Ta=发生短路If波形 短路时刻的不同对短路电流的影响 由于短路电流的直流分量起始值越大，短路电流瞬时值就越大，而直流分量的起始值于短路时刻的电流相位有关，即直流分量是由于短路后电流不能突变而产生的。 Pscad模型中对短路时刻的设置如图6所示 图6 Pscad对于短路时刻的设置 1）当在t=时发生三相短路，三相短路电流波形如图7所示。 图7 t=时三相短路电流波形 2）当在t=时发生三相短路，三相短路电流波形如图8所示。 图8 t=6时三相短路电流波形 Xd、Xd`、Xd``对短路电流的影响 1) Xd的影响 Pscad中对于Xd的设置如图9所示： 图9 Pscad对于D轴同步电抗Xd的设置 下面验证不同Xd时A相短路电流的稳定值。 i.Xd=（标幺制，下同）时，仿真波形如图10所示 图10 Xd=时A相短路电流波形 ii.Xd=10时，仿真波形如图11所示 图11 Xd=时A相短路电流波形 2）Xd`的影响 在Pscad中暂态电抗Xd`的设置如图13所示： 图13 Pscad对于暂态电抗Xd的设置 下面验证不同Xd`时A相短路电流的暂态过程。 i.Xd`=时A相短路电流的波形如图14所示： 图14 Xd`=时A相短路电流波形 ii.Xd`=1时A相短路电流的波形如图15所示： 图15 Xd``=1时A相短路电流波形 3）Xd``的影响 这里次暂态电抗Xd``与暂态电抗Xd`相似，Xd``影响的是短路后的次暂态过程。
电力系统分析潮流计算程序设计报告题目：13节点配电网潮流计算 学院电气工程学院 专业班级 学生姓名 学号 班内序号 指导教师房大中 提交日期 2015年05月04日
目录 一、程序设计目的 (1) 二、程序设计要求 (3) 三、13节点配网潮流计算 (3) 3.1主要流程................................................................................................... 错误!未定义书签。 3.1.1第一步的前推公式如下（1-1）-（1-5）： ....................................... 错误!未定义书签。 3.1.2第二步的回代公式如下（1-6）—（1-9）： ..................................... 错误!未定义书签。 3.2配网前推后代潮流计算的原理 (6) 3.3配网前推后代潮流计算迭代过程 (7) 3.3计算原理 (8) 四、计算框图流程 (9) 五、确定前推回代支路次序.......................................................................................... 错误!未定义书签。 六、前推回代计算输入文件 (10) 主程序： (10) 输入文件清单： (11) 计算结果： (12) 数据分析： (12) 七、配电网潮流计算的要点 (13) 八、自我总结 (13) 九、参考文献 (14) 附录一 MATLAB的简介 (14)