PSO_dg_prog

PARTICLE SWARM OPTIMIZTION FOR DG PLACEMENT

% The implementation of PSO program to find the location and size of DG. 
% The user can make m-files for bus data line data and additional file myloss.m 
% before using this program as per the user requirement with support files.

% The basic reffence papers for this program  “James Kennedy and Russell Eberhart“ on
% PSO
% Author: AASHISH KUMAR BOHRE MANIT BHOPAL MP INDIA

clear all
clc

[basemva busdata linedata]= Mydgalldata;
up=28; low=0.001; accuracy = 0.001;  maxiter = 10;
qwe=busdata;
for bss=4:1:16
    aasuu=aasuu1;
    aasuu（bss1）=2;
             n=length（qwe（:1））;
         Dn= n-6;           
         itera_max  = 10; c1 =2; c2 = 2; R1 = 0.5; R2 = 0.5;
           Pdg=rand（Dn1）.*（up-low）+low;
            current_pop = [Pdg];
           current_fitnsssss=ones（Dn1）.*inf;                     
           local_best_fitnsssss  = current_fitnsssss;           
           [global_best_fitnsssssg] = min（local_best_fitnsssss）;
           local_best_pop  = current_pop;
           globl_best_pop = local_best_pop（g:）;
            velocity = rand（Dn1）.*0.5;
        itera = 0 ;
      while  （ itera < itera_max ）
        itera = itera + 1;
        w=（itera_max-itera）/itera_max;
         velocity = w *velocity + c1*rand.*（local_best_pop-current_pop）+...
             c2*rand.*（（ones（Dn1）*globl_best_pop）-current_pop）;
         current_pop = current_pop + velocity;        
              for i=1:Dn
                  if current_pop（i:）                      current_pop（i:）=low;
                  elseif current_pop（i:）>up
                      current_pop（i:）=up;
                  end
              end
              p_g=gar;
            for dgb=1:1:Dn
                p_g（bss1）=current_pop（dgb1）;
                aai=p_g（bss1）;
                
                myloss 
                ajj = real（myloss）;
             
              fitP_LOSS（dgb）=ajj;
            end
            current_fitnsssss=fitP_LOSS‘;
        for i = 1 : Dn
                if current_fitnsssss（i:） < local_best_fitnsssss（i:）
                   local_best_fitnsssss（i:）  = current_fitnsssss（i:）;  
                   local_best_pop（i:） = current_pop（i:）;
                end   
        end
         [current_global_best_fitnsssssg] = min（local_best_fitnsssss）;
        if current_global_best_fitnsssss < global_best_fitnsssss
           global_best_fitnsssss = current_global_best_fitnsssss;
        end   
         globl_best_pop = local_best_pop（g:）;
        minfit（itera）=global_best_fitnsssss;
        itera
      end   
        global_best_fitnsssss   
        globl_best_pop         
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1509  2020-09-19 08:46  license.txt
     文件        2620  2020-09-19 08:46  PSO_dg_prog.m