含分布式电源的配网三相不平衡潮流计算程序（前推回代法）
目录
1 主要内容
2 部分代码
3 程序结果
4 下载链接
1 主要内容
该程序主要采用前推回代法的三相不平衡潮流计算模型，采用对称分量法，并提出考虑相间互阻抗的精确网络损耗计算方法，采用IEEE33节点模型，通过改变branch矩阵中线路不同相阻抗参数，得到线路参数不平衡条件下的潮流结果，通过不同相增加分布式电源，也得到系统三相不平衡结果，该程序通用性强，学习三相不平衡潮流计算的必备程序。
2 部分代码
al=exp(j*2*pi/3);%
SS=[1,1,1;1,al^2,al;1,al,al^2];%对称分量变换矩阵
%% 数据初始化：branch、bus为matpower格式，A、S、U为本程序格式
% fbus tbus r g x b rateA rateB rateC ratio angle status angmin angmax
branch = [
1 2 0.0922 0.0922 0.0922 0 0 0 0.047 0.047 0.047 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 1 1
2 3 0.4930 0.3930 0.4930 0 0 0 0.2511 0.2511 0.2511 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 2 1
3 4 0.3660 0.3660 0.3660 0 0 0 0.1864 0.1864 0.1864 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 3 1
4 5 0.3811 0.3811 0.3811 0 0 0 0.1941 0.1941 0.1941 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 4 1
5 6 0.8190 0.8190 0.8190 0 0 0 0.7070 0.7070 0.7070 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 4 1
6 7 0.5872 0.1872 0.1872 0 0 0 0.6188 0.6188 0.6188 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 5 1
7 8 0.7114 0.7114 0.7114 0 0 0 0.2351 0.2351 0.2351 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 6 1
8 9 1.0300 1.0300 1.0300 0 0 0 0.7400 0.7400 0.7400 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 7 1
9 10 1.0400 1.0400 1.0400 0 0 0 0.7400 0.7400 0.7400 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 8 1
10 11 0.1966 0.1966 0.1966 0 0 0 0.0650 0.0650 0.0650 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 9 1
11 12 0.3744 0.3744 0.3744 0 0 0 0.1238 0.1238 0.1238 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 10 1
12 13 1.4680 1.4680 1.4680 0 0 0 0.1550 0.1550 0.1550 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 11 1
13 14 0.5416 0.5416 0.5416 0 0 0 0.7192 0.7192 0.7192 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 12 1
14 15 0.5910 0.5910 0.4910 0 0 0 0.5260 0.5260 0.5260 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 13 1
15 16 0.7463 0.7463 0.7463 0 0 0 0.5450 0.5450 0.5450 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 14 1
16 17 1.2890 1.2890 1.2890 0 0 0 1.7210 1.7210 1.7210 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 15 1
17 18 0.3720 0.3720 0.3720 0 0 0 0.5740 0.5740 0.5740 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 16 1
%17 18 0.7429 0.7407 0.7393 0.0073 0.0109 0.0087 0.5826 0.5808 0.5797 0.0057 0.0086 0.0068 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 17 1
2 19 0.1640 0.1640 0.1640 0 0 0 0.1565 0.1565 0.1565 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 18 1
19 20 1.5042 1.5042 1.5042 0 0 0 1.3554 1.3554 1.3554 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 19 1
20 21 0.4095 0.4095 0.4095 0 0 0 0.4784 0.4784 0.4784 0 0 0 0 0 0 2.356 5.0593 250 250 250 5.9859 0 1 1 1 -360 360 20 1
21 22 0.7089 0.7089 0.7089 0 0 0 0.9373 0.9373 0.9373 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 21 1
3 23 0.4512 0.4512 0.4512 0 0 0 0.3083 0.3083 0.3083 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 22 1
23 24 0.8980 0.8980 0.8980 0 0 0 0.7091 0.7091 0.7091 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 23 1
24 25 0.8960 0.8960 0.8960 0 0 0 0.7011 0.7011 0.7011 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 23 1
6 26 0.2030 0.1030 0.2030 0 0 0 0.1034 0.1034 0.1034 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 25 1
26 27 0.2842 0.2842 0.3842 0 0 0 0.1477 0.1477 0.1477 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 26 1
27 28 1.0590 1.0590 1.0590 0 0 0 0.9377 0.9377 0.9377 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 27 1
28 29 2.3942 0.8042 0.8042 0 0 0 0.7006 0.7006 0.7006 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 28 1
29 30 0.5075 0.5075 0.5075 0 0 0 0.2585 0.2585 0.2585 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 29 1
30 31 0.9704 0.9704 0.9704 0 0 0 0.9630 0.9630 0.9630 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 30 1
31 32 0.3105 0.2105 0.3105 0 0 0 0.3619 0.3619 0.3619 0 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 30 1
32 33 0.3410 0.3410 0.3410 0 0 0 0.5381 0.5362 0.5362 0.5362 0 0 0 0 0 0 0 250 250 250 0 0 1 1 1 -360 360 32 1
];
% % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
branch(:,6:8)=0.1.*branch(:,3:5);%更新g
branch(:,12:14)=0.1.*branch(:,9:11);%更新b3 程序结果
4 下载链接
点击直达！