本科生实验报告
课程名称： 电力系统分析实验
专业班级： 电力系统124 班 姓 名： 学 号：
所在学期： 2014-2015-2
2015年 6 月 20 日
实验四 配电网潮流计算实验
一、实验目的
本实验针对配电网具有辐射式开环结构的特性，基于前推回代的方法编制程序使系统潮流计算能够由计算机自行完成，即已知配电网首端节点电压和各末端节点的功率，由末端向首端计算网络功率分布，再由首端向末端计算节点电压分布，迭代直至收敛条件。通过实验教学加深学生对配电网潮流计算方法的理解，掌握树状辐射式配电网潮流计算的基本算法。熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。 二、实验器材
计算机、软件（已安装，包括各类编程软件C语言、C++、VB、VC等、应用软件MATLAB等）、移动存储设备（学生自备，软盘、U盘等） 三、实验内容
编制调试配电网潮流计算的计算机程序。程序要求根据已知的配电网参数，完成该配电系统的潮流计算，要求计算出节点电压、功率等参数。
1.先编制好的电力系统潮流计算的计算程序原代码由自备移动存储设备导入计算机。
2.应的编程环境下对程序进行组织调试。 3.应用计算例题验证程序的计算效果。 4.对调试正确的计算程序进行存储、打印。 5.完成本次实验的实验报告。 四、实验数据
31+j2? 11+j1?2S3M10kVS1S22+j2? 42+j1? 5S5
如图所示一个5节点的配电网系统，S1=S2=2+j2(MV·A)，S3=S4=S5=4+j4(MV·A)，假定所有负荷均为恒功率负荷，节点1为参考节点相
?6角为0o。计算电网个节点电压及支路功率分布，收敛条件为??10。
S4五、实验程序
% 配电网前推回代潮流计算程序
% 使用IEEE 33节点配电系统作为算例，可实现弱环网情况下的潮流计算
countnum=0; BranchData = [
1 2 0.0922 0.0470; 2 3 0.4930 0.2511; 3 4 0.3660 0.1864; 4 5 0.3811 0.1941; 5 6 0.8190 0.7070; 6 7 0.1872 0.6188; 7 8 0.7114 0.2351; 8 9 1.0300 0.7400; 9 10 1.0440 0.7400; 10 11 0.1966 0.0650; 11 12 0.3744 0.1238; 12 13 1.4680 1.1550; 13 14 0.5416 0.7129; 14 15 0.5910 0.5260; 15 16 0.7463 0.5450; 16 17 1.2890 1.7210; 17 18 0.7320 0.5740; 2 19 0.1640 0.1565; 19 20 1.5042 1.3554; 20 21 0.4095 0.4784; 21 22 0.7089 0.9373; 3 23 0.4512 0.3083; 23 24 0.8980 0.7091; 24 25 0.8960 0.7011; 6 26 0.2030 0.1034; 26 27 0.2842 0.1447; 27 28 1.0590 0.9337; 28 29 0.8042 0.7006; 29 30 0.5075 0.2585; 30 31 0.9744 0.9630; 31 32 0.3105 0.3619; 32 33 0.3410 0.5302;
];
NodeData = [
2 100.00 60.00; 3 90.00 40.00; 4 120.00 80.00;
% 支路，阻抗 5 60.00 30.00; 6 60.00 20.00; 7 200.00 100.00; 8 200.00 100.00; 9 60.00 20.00; 10 60.00 20.00; 11 45.00 30.00; 12 60.00 35.00; 13 60.00 35.00; 14 120.00 80.00; 15 60.00 10.00; 16 60.00 20.00; 17 60.00 20.00; 18 90.00 40.00; 19 90.00 40.00; 20 90.00 40.00; 21 90.00 40.00; 22 90.00 40.00; 23 90.00 50.00; 24 420.00 200.00; 25 420.00 200.00; 26 60.00 25.00; 27 60.00 25.00; 28 60.00 20.00; 29 120.00 70.00; 30 200.00 600.00; 31 150.00 70.00; 32 210.00 100.00; 33 60.00 40.00;
]; % 节点，负荷
UB = 12.66; % 电压基准 kV SB = 10; % 功率基准 MVA ZB = UB^2/SB; % 阻抗基准 ohm
BranchData(:,[3,4]) = BranchData(:,[3,4]) / ZB; % 阻抗标幺化 NodeData(:,[2,3]) = NodeData(:,[2,3]) / SB / 1000;% 功率标幺化
NN = 33; % 节点数
A0 = zeros(NN); for n = 1:NN-1
A0(BranchData(n,1),BranchData(n,2)) = 1;
end % 形成 A0
AssociatedMatrix=0;
for n=2:NN-1
AssociatedMatrix(n,n)=1; temp=BranchData(n-1,1);
AssociatedMatrix(n,1:n-1)=AssociatedMatrix(temp,1:n-1); end
A0T = A0'; % 形成 A0 的转置 S = [0;-NodeData(:,2) - i*NodeData(:,3)]; % 形成 S
ZL = [0;BranchData(:,3) + i*BranchData(:,4)]; % 形成 ZL
V = ones(NN,1);
V(1) = 1; % 各个节点电压赋初值
IL(NN,1) = -conj(S(NN) / V(NN)); % 最末支路电流赋初值
Delta = 1; % 收敛判据赋初值
TempV = V; % 赋初值，用于记忆上次迭代结果
while Delta > 1e-8
countnum=countnum+1;
IN = conj(S ./ V); % 节点注入电流
for n = 1:NN-1
IL(NN-n) = A0(NN-n,NN-n+1:end) * IL(NN-n+1:end) - IN(NN-n); end % 电流回代过程
for n = 2:NN
V(n) = A0T(n,1:n-1) * V(1:n-1) - ZL(n) * IL(n);
end % 电压前推过程
Delta = max(abs(V-TempV)); % 更新收敛判据 TempV = V; % 记忆迭代结果 end
Vangle(:,1)=abs(V);
Vangle(:,2)=angle(V)/3.1415*180;
for i=1:NN-1
st=BranchData(i,1); en=BranchData(i,2);