IEEE 39节点系统使用matlab进行潮流计算,完美符合matpower数据结果,内含重新编号代码。
%数据
%负荷有功无功
%编号 类型 有功 无功 幅值初值 相角初值 基准电压
%bus_i type Pd Qd Vm Va baseKV
data.bus = [
1 1 97.6 44.2 1 0 345;
2 1 0 0 1 0 345;
3 1 322 2.4 1 0 345;
4 1 500 184 1 0 345;
5 1 0 0 1 0 345;
6 1 0 0 1 0 345;
7 1 233.8 84 1 0 345;
8 1 522 176.6 1 0 345;
9 1 6.5 -66.6 1 0 345;
10 1 0 0 1 0 345;
11 1 0 0 1 0 345;
12 1 8.53 88 1 0 345;
13 1 0 0 1 0 345;
14 1 0 0 1 0 345;
15 1 320 153 1 0 345;
16 1 329 32.3 1 0 345;
17 1 0 0 1 0 345;
18 1 158 30 1 0 345;
19 1 0 0 1 0 345;
20 1 680 103 1 0 345;
21 1 274 115 1 0 345;
22 1 0 0 1 0 345;
23 1 247.5 84.6 1 0 345;
24 1 308.6 -92.2 1 0 345;
25 1 224 47.2 1 0 345;
26 1 139 17 1 0 345;
27 1 281 75.5 1 0 345;
28 1 206 27.6 1 0 345;
29 1 283.5 26.9 1 0 345;
30 2 0 0 1.0499 0 345;
31 3 9.2 4.6 0.982 0 345;
32 2 0 0 0.9841 0 345;
33 2 0 0 0.9972 0 345;
34 2 0 0 1.0123 0 345;
35 2 0 0 1.0494 0 345;
36 2 0 0 1.0636 0 345;
37 2 0 0 1.0275 0 345;
38 2 0 0 1.0265 0 345;
39 2 1104 250 1.03 0 345;
];
%发电机有功无功
%编号 有功 无功 工作电压 功率基准 运行状态
%bus Pg Qg Vg mBase status
data.gen = [
30 250 0 1.0499 100 1;
% 31 677.871 0 0.982 100 1;
32 650 0 0.9841 100 1;
33 632 0 0.9972 100 1;
34 508 0 1.0123 100 1;
35 650 0 1.0494 100 1;
36 560 0 1.0636 100 1;
37 540 0 1.0275 100 1;
38 830 0 1.0265 100 1;
39 1000 0 1.03 100 1;
];
%线路参数
%起点编号 终点编号 r x b 支路变比 支路的相位角度 支路工作状态
%fbus tbus r x b ratio angle status
data.line = [
1 2 0.0035 0.0411 0.6987 0 0 1;
1 39 0.001 0.025 0.75 0 0 1;
2 3 0.0013 0.0151 0.2572 0 0 1;
2 25 0.007