基于粒子群算法的配电网无功优化.rar

可在matlab中直接运行出结果。通过matlab程序用粒子群算法对无功补偿容量进行调整，得到一组最优无功补偿装置容量，将装置投入IEEE33节点得到优化后的节点电压和系统网损，可验证粒子群算法在配电网中无功优化的可行性。

function result=fitness11（xk）
global B1;
global B2;




%相关原始数据格式说明如下：
%n——节点个数；n1——支路条数；isb——平衡节点号；H——PQ节点个数（为后面形成PVU存储PV节点初始电压用）；pr——误差精度。
%B1——支路参数矩阵，其中第一列和第二列是起始节点编号和终点节点编号第三列、第四列、第五列、第六列分别为：支路电阻、电抗、变压器变比、电纳。（不考虑电导）
%B2——节点参数矩阵，其中第一列和第二列为节点编号和节点类型；第三列到第六列分别为：注入有功、注入无功、电压幅值、电压相位。
%节点类型分类如下：“0”为平衡节点“1”为PQ，“2”为PV节点；“3”为PQ（V）节点，“4”为PI节点。
n=33 ;       
n1=32;
isb=1;
H=32;                %%%%%%%%%%%%%18节点加DG    PQV处理
pr=0.001;   
v_amp=0;
Vimax=1.15;Vimin=0.85;
Z1=1;
 
B2=[1 0 0 0 1.05 0;
    2 1 -0.01 -0.006 1 0;
    3 1 -0.009 -0.004 1 0;
    4 1 -0.012 -0.008 1 0;
    5 1 -0.006 -0.003 1 0;
    6 1 -0.006 -0.002 1 0;
    7 1 -0.02 -0.01 1 0;
    8 1 -0.02 -0.01 1 0;
    9 1 -0.006 -0.002 1 0;
    10 1 -0.006 -0.0035 1 0;
    11 1 -0.0045 -0.003 1 0;
    12 1 -0.006 -0.0035 1 0;
    13 1 -0.006 -0.0035 1 0;
    14 1 -0.012 -0.008 1 0;
    15 1 -0.006 -0.001 1 0;
    16 1 -0.006 -0.002 1 0;
    17 1 -0.006 -0.002 1 0;
    18 1 -0.009 -0.004 1 0;
    19 1 -0.009 -0.004 1 0;
    20 1 -0.009 -0.004 1 0;
    21 1 -0.009 -0.004 1 0;
    22 1 -0.009 -0.004 1 0;
    23 1 -0.009 -0.005 1 0;
    24 1 -0.042 -0.02 1 0;
    25 1 -0.042 -0.02 1 0;
    26 1 -0.006 -0.0025 1 0;
    27 1 -0.006 -0.0025 1 0;
    28 1 -0.006 -0.002 1 0;
    29 1 -0.012 -0.007 1 0;
    30 1 -0.02 -0.06 1 0;
    31 1 -0.015 -0.007 1 0;
    32 1 -0.021 -0.01 1 0;
    33 1 -0.006 -0.004 1 0];

B2（54）=B2（24）+x（1）;
B2（154）=B2（154）+x（2）;
B2（314）=B2（314）+x（3）;



Y=zeros（n）;               %zeros就是生成一个全0的矩阵
Times=1;                  %置迭代次数为初始值
%创建节点导纳矩阵   电导统一为零


for i=1:n1
        p=B1（i1）;
        q=B1（i2）;
        Y（pq）=Y（pq）-1/（（B1（i3）+B1（i4））*B1（i5））;
        Y（qp）=Y（pq）;
        Y（pp）=Y（pp）+1/（B1（i3）+B1（i4））+0.5*B1（i6）;%低压侧不变
        Y（qq）=Y（qq）+1/（（B1（i3）+B1（i4））*B1（i5）^2）+0.5*B1（i6）;%高压侧阻抗乘以变比平方  输入时注意低压侧在前
end
%disp（‘节点导纳矩阵：‘） ;
Y;
G=real（Y）;
B=imag（Y）;
OrgS=zeros（2*n-21）;
DetaS=zeros（2*n-21）;   %将OrgS、DetaS初始化
%创建OrgS，用于存储初始功率参数
Q=0;
PQV=0;
x=1.655;         %%%x1=6.7x2=9.85x=x1+x2；其中x1为定子漏抗，x2为转子漏抗
xp=18.8;         %%%xc=xm=xp=xc*xm/（xc-xm）；其中xc为机端并联电容器电抗，xm为激磁电抗
h=0;

for i=1:n    %对PQ（V）节点的处理
    h=h+1;
    if i~=isb&&B2（i2）==3   
         Q（i）=-（B2（i5））^2/xp+（-（B2（i5））^2+sqrt（（B2（i5））^4-4*（B2（i3））^2*x^2））/2*x;
         B2（i4）=Q（i）;
         B2（i2）=1; 
         PQV=h;
    end  
end
       %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Ig=0.01;
Q=0;
PI=0;
h=0;
for i=1:n    %对PI节点的处理
    h=h+1;
    if i~=isb&&B2（i2）==4   
         Q（i）=sqrt（Ig^2*（（B2（i5））^2）-B（i3）^2）;        %e=B2（i5），f=0，e^2+f^2=B2（i5））^2，其中e和f为光伏发电系统接入节点电压的实部和虚部
         B2（i4）=Q（i）;
         B2（i2）=1; 
         PI=h;
    end  
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
h=0;
j=0;
for i=1:n            %对PQ节点的处理
    if i~=isb&&B2（i2）==1   
        h=h+1;
        for j=1:n
            OrgS（2*h-11）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      14512  2019-05-14 19:04  基于粒子群算法的配电网无功优化\fitness11.m

     文件      15782  2019-05-23 15:49  基于粒子群算法的配电网无功优化\powerieee.m

     文件       4403  2019-05-23 15:36  基于粒子群算法的配电网无功优化\QO.m

     文件        161  2019-05-14 12:52  基于粒子群算法的配电网无功优化\必读.txt

     目录          0  2019-12-27 16:28  基于粒子群算法的配电网无功优化

----------- ---------  ---------- -----  ----

                34858                    5