粒子群算法无功优化程序

基于粒子群算法的配电网无功优化程序，利用了IEEE30节点进行MATLAB编程

function [ V ] = act（）
digits（6）;
%初始化，系统数据录入
clear ObjV;
clear objv;
[YnmPpvPpqQpqVpvVpqcita]=initial（）;
%定义遗传算法参数
NIND=40;           %个体数目
global MAXGEN;        %最大遗传代数
NVAR=6;            %变量维数
PRECI=30;          %变量的二进制位数
GGAP=0.5;          %代沟
global trace;
trace=zeros（MAXGEN2）;
%建立区域描述器
FieldD=[rep（[PRECI][1NVAR]）;rep（[0.9700;1.0300][1NVAR]）;rep（[1;0;1;1][1NVAR]）];
Chrom=crtbp（NINDNVAR*PRECI）;             %创建初始种群
gen=0;                                    %代计数器
%计算初始种群的目标函数值
%-------------------------------------------------------
%潮流计算
%n--节点总数，n-m--PV节点数量，m-1--PQ节点数量
%[PcitaV]=pf（）;
global V;
global cita;
global P;

detaP=zeros（n-11）;
detaQ=zeros（m-11）;
%P=[Ppq;Ppv];
%V=[Vpq;Vpv;1];
Q=Qpq;
G=real（Y）;
B=imag（Y）;
cita1=[cita];
n=30;ObjV=0;
Vpv0=bs2rv（ChromFieldD）;             %控制变量转为十进制
for k=1:NIND;
    V（n-6）=Vpv0（k1）;
    V（n-5）=Vpv0（k2）;
    V（n-4）=Vpv0（k3）;
    V（n-3）=Vpv0（k4）;
    V（n-2）=Vpv0（k5）;
    V（n-1）=Vpv0（k6）;
%求detaP
for i=1:n-1;
    s=0;
    for j=1:n;
        s=s+V（j）*（G（ij）*cos（cita1（i）-cita1（j））+B（ij）*sin（cita1（i）-cita1（j）））;
    end
    detaP（i）=P（i）-V（i）*s; 
end
ObjV（k）=0;
ObjV（k）=ObjV（k）+detaP（i）;
ObjV（k）=abs（ObjV（k））;                   %得到初始种群目标函数值
end

global ObjV;
ObjV=reshape（ObjVNIND1）;
global objv;
objv = ObjV;
while gen    FitnV=ranking（ObjV）;                  %分配适应度值
    SelCh=select（‘sus‘ChromFitnVGGAP）; %选择
    SelCh=recombin（‘xovsp‘SelCh0.7）;    %重组
    SelCh=mut（SelCh）;                     %变异
    Vpv0=bs2rv（SelChFieldD）;             %控制变量转为十进制
   % disp（Vpv0）;
%计算子代目标函数值
    for m=1:NIND*GGAP;
    V（n-6）=Vpv0（m1）;
    V（n-5）=Vpv0（m2）;
    V（n-4）=Vpv0（m3）;
    V（n-3）=Vpv0（m4）;
    V（n-2）=Vpv0（m5）;
    V（n-1）=Vpv0（m6）;
%求detaP
for i=1:n-1;
    s=0;
    for j=1:n;
        s=s+V（j）*（G（ij）*cos（cita1（i）-cita1（j））+B（ij）*sin（cita1（i）-cita1（j）））;
    end
    detaP（i）=P（i）-V（i）*s;
end
    ObjVSel（m）=0;
    ObjVSel（m）=ObjVSel（m）+detaP（i）;
    ObjVSel（m）=abs（ObjVSel（m））;
end
    ObjVSel=reshape（ObjVSelNIND*GGAP1）;
    [Chrom ObjV]=reins（ChromSelCh11ObjVObjVSel）;
    %disp（ObjVSel）;
    gen=gen+1;
    %输出最优解及其对应的30个自变量的十进制值，Y为最优解，I为种群的序号
    trace（gen1）=min（ObjV）;               %遗传算法性能跟踪
    trace（gen2）=sum（ObjV）/length（ObjV）;
end
    % plot（ObjV）;hold on;
    % plot（ObjV‘b*‘）;grid;
   % plot（trace（:1））; hold on;
   
   % plot（trace（:2）‘-.‘）;grid;
   % 
%legend（‘种群均值的变化‘‘解的变化‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-05-17 20:04  程序\
     文件        2751  2014-05-17 10:12  程序\act.asv
     文件        2737  2014-05-17 10:22  程序\act.m
     文件        1231  2014-05-08 21:01  程序\B1.txt
     文件        1906  2014-05-09 11:07  程序\B2.txt
     文件         160  2013-12-11 14:36  程序\detacitaV.m
     文件         626  2013-12-16 09:10  程序\detaPQ.m
     文件        1156  2014-05-15 21:32  程序\initial.asv
     文件        1153  2014-05-10 22:23  程序\initial.m
     文件        2183  2014-05-10 10:15  程序\jacobi.asv
     文件        2183  2014-05-10 10:18  程序\jacobi.m
     文件       12071  2014-05-17 10:17  程序\mai.asv
     文件        9984  2014-05-17 21:02  程序\mai.fig
     文件       12071  2014-05-17 10:21  程序\mai.m
     文件        1091  2014-05-10 22:47  程序\pf.asv
     文件        1106  2014-05-10 23:36  程序\pf.m
     文件         317  2014-05-15 21:32  程序\read_B1.asv
     文件         314  2014-05-09 16:50  程序\read_B1.m
     文件         314  2014-05-09 16:58  程序\read_B2.m
     文件          88  2013-12-10 09:45  程序\update.m