PGSA主代码

clear all; 
clc; 
tic 
global mvab kvb busdata linedata nf ng Ntiesw Lsw Lswcnt Fbest; 
data16; 
Xb=[5 11 16]; 
[TPLVolt]=loadflow（Xb）; 
TPLold=TPL; 
Fbest=TPL; 
for i=1:ng %maxno.of iterations     
    [SolutSet]=searchdomain（Xb）; 
     u=length（SolutSet（:1））; 
     for j=1:u %length（SolutSet（:1）） 
         clear TPL Volt;  
         [TPLVolt]=loadflow（SolutSet（j:））; 
         g1（j）=TPL; 
         g2（j:）=Volt; 
     end 
     g1; 
     fbest=min（g1）; 
     k=min（find（fbest==g1））; 
     Xbest=SolutSet（k:）; 
     if（fbest         Fbest=fbest; 
         Xb=Xbest;