%?改進(jìn)的快速粒子群優(yōu)化算法?（APSO）
function?apso

clear?all;
Lb=[0.9?0.9?0?0.75?0?0];?%下邊界
Ub=[1?1?1?1?12?2500];?%上邊界
%?默認(rèn)參數(shù)
para=[25?250?0.95];?%[粒子數(shù)，迭代次數(shù)，gama參數(shù)]

%?APSO?優(yōu)化求解函數(shù)
[gbestfmin]=pso_mincon（@cost@constraintLbUbpara）;

%?輸出結(jié)果
Bestsolution=gbest?%?全局最優(yōu)個(gè)體
fmin

%%?目標(biāo)函數(shù)
function?f=cost（x）
%s=[0;0;0;0;0;0;1;1;1;1;1;1;1;1;1;1;1;1;1;1;0;0;0;0];%內(nèi)燃機(jī)工作狀態(tài)
%a=[1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1];%制冷劑工作12小時(shí)
%c=[0.524;0.524;0.524;0.524;0.524;0.524;0.524;0.524;0.875;0.875;0.875;1.217;1.217;1.217;1.217;1.217;0.875;0.875;0.875;1.217;1.217;1.217;0.875;0.875];%電價(jià)
f=x（1）*0.059*700+x（2）*0.039*4000*5+x（3）*0.0013*13800+x（4）*0.013*1200+x（5）*1000*0.43+x（6）*0.5;
%?非線性約束
function?[ggeq]=constraint（x）
%?不等式限制條件
%Pd=xlsread（‘春夏負(fù)荷.xlsx‘‘sheet1‘‘C1:C24‘）;
%Qh=xlsread（‘春夏負(fù)荷.xlsx‘‘sheet1‘‘B1:B24‘）;
%Cload=xlsread（‘春夏負(fù)荷.xlsx‘‘sheet1‘‘A1:A24‘）;
%?如果沒有等式約束，則置geq=[]；
Pd=9.95;
Qh=17.06;
Cload=1.8;
%Pd=xlsread（‘春夏負(fù)荷.xlsx‘‘sheet1‘‘C1:C24‘）;
%Qh=xlsread（‘春夏負(fù)荷.xlsx‘‘sheet1‘‘B1:B24‘）;
%Cload=xlsread（‘春夏負(fù)荷.xlsx‘‘sheet1‘‘A1:A24‘）;
%?如果沒有等式約束，則置geq=[]；
g（1）=x（5）+x（1）*0.7+x（2）*4-Pd*1.0001-3.94*x（3）;
g（2）=Pd*0.9999-（x（5）+x（1）*0.7+x（2）*4-3.94*x（3））;
g（3）=x（1）+x（2）*4+0.8*x（6）/120-x（4）-Qh*1.0001;
g（4）=0.9999*Qh-（x（1）+x（2）*4+0.8*x（6）/120-x（4））;
g（5）=1.2*x（4）+13.8*x（3）-Cload*1.0001;
g（6）=Cload*0.9999-（1.2*x（4）+13.8*x（3））;
%if?x（1）<0.5?x（1）=0;elsex（1）=1;end
geq=[];


%%??APSO?Solver
function?[gbestfbest]=pso_mincon（fhandlefnonlinLbUbpara）
if?nargin<=6
????para=[25?250?0.95];
end
n=para（1）;%?粒子種群大小
time=para（2）;?%時(shí)間步長，迭代次數(shù)
gamma=para（3）;?%gama參數(shù)
?scale=abs（Ub-Lb）;?%取值區(qū)間
%?驗(yàn)證約束條件是否合乎條件
if?abs（length（Lb）-length（Ub））>0
????disp（‘Constraints?must?have?equal?size‘）;
????return
end

??alpha=0.2;?%?alpha=[01]粒子隨機(jī)衰減因子
??beta=0.5;??%?收斂速度（0->1）=（slow->fast）;

%?初始化粒子群
best=init_pso（nLbUb）;