%?NP---?????????
%?G----?????
%ITERmax--??????
%?Zmax---???????
%?Zmin---???????
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%???????????%%%%%%%%%%%%%%%%%%%%%%%
function?syst
clear?all;
clc;
%format?long;
NP=40;
G=10;
InitialTON=[8800000000];
%InitialTON=[8855555555];
%InitialTOFF=[0055633111];
InitialTOFF=[0000000000];
MDT=[8855633111];%min?down?time
MUT=[8855633111];%min?up?time
x0=[1100000000];%初始機組方案
%%pwind=[42?37?60?19?20?26?40?37?21?16?23?24?41?46?80?99?72?50?38?59?66?86?71?61];
pwind=[32.29?30.87?37.94?34.37?36.09?41.64?50.15?44.28?38.46?33.38?25.87?25.01?22.67?54.13?54.00?59.29?53.84?45.48?47.61?47.71?59.43?41.16?46.22?59.12];
T=24;
PMAX=[4554551301301628085555555];
PMIN=[1501502020252025101010];
PD=[700750850950100011001150120013001400145015001400130012001050100011001200140013001100900800];
SR=0.1*PD;


a=[1000970700680450370480660665670];
b=[16.1917.2616.616.519.722.2627.7425.9227.2727.79];
c=[0.000480.000310.0020.002110.003980.007120.000790.004130.002220.00173];
SUH=[45005000550560900170260303030];%hot?start?cost
SUC=[900010000110011201800340520606060];%?cold?start?cost
Tcold=[5544422000];%cold?start?hours

ITERmax=10;
Zmax=0.05*pi;
Zmin=0.01*pi;
[Yh]=sort（PMAX‘descend‘）;%降序排序，Y是排序后的向量，h是排序的順序

swarm=zeros（NPG*T）;
ab=1/sqrt（2）*ones（2G*T）;
%***************???NP??????**********************
for?n=1:NP
????%--??
????swarml=initial（GT）;
????%??????????????
????[swarm（n:）swarmTOFF（::n）]=Verify（swarmlGTx0MUTMDTInitialTONInitialTOFFPMAXPDSRh）;
????%*****************?????????***********************
????[fitness（n）]=objectFitness（swarm（n:）GTx0MDTPMAXPMINPDabcSUHSUCTcoldswarmTOFF（::n））;
end
%*************?????????????????*********************
Gmbest=swarm（1:）;%將swarm矩陣中的第一行附給Gmbest
Gfitness=fitness（1）;%Gfitness=fitness矩陣中的第一個元素
for?n=1:NP
????Pmbest（n:）=swarm（n:）;
????Pfitness（n）=fitness（n）;
????if?Gfitness>Pfitness（n）
????????Gmbest=Pmbest（n:）;
????????Gfitness=Pfitness（n）;
????end
end

%**************???????Updata********************
for?n=1:NP
????ab1（::n）=ab;
end
for?k=1:ITERmax
????for?n=1:NP
????????%??????
????????[swarm（n:）ab1（::n）]=Update（swarm（n:）GTPmbest（n:）Gmbestfitness（n）Pfitness（n）GfitnessZmaxZminkITERmaxab1（::n））;
???????%?????
????????[swarm（n:）swarmTOFF（::n）]=Verify（swarm（n:）GTx0MUTMDTInitialTONInitialTOFFPMAXPDSRh）;
????????%??????
????????[fitness（n）]=objectFitness（swarm（n:）GTx0MDTPMAXPMINPDabcSUHSUCTcoldswarmTOFF（::n））;
???%?end
???%?for?n=1:NP
????????if?Pfitness（n）>=fitness（n）
????????????Pmbest（n:）=swarm（n:）;
???????????%?Ptmbest（::n）=Mp（::n）;
????????????Pfitness（n）=fitness（n）;
????????end
????????if?Gfitness>Pfitness（n）
????????Gmbest=Pmbest（n:）;
???????%?Gptbest=Ptmbest（::n）;
????????Gfitne