clc
clear
load?data1
load?data2
data（1:14:）=data1（1:14:）;
data（15:69:）=data2（1:55:）;
input=data（:1:2）;
output?=data（:3）;
input_train=input（（1:60）:）‘;
output_train=output（（1:60）:）‘;
input_test=input（（61:69）:）‘;
output_test=output（（61:69）:）‘;
[inputninputps]=mapminmax（input_train）;
[outputnoutputps]=mapminmax（output_train）;
net=newff（inputnoutputn8）;
c1?=?2;
c2?=?2;

maxgen=100;???
sizepop=100;??

Vmax=1;
Vmin=-1;
popmax=5;
popmin=-5;


for?i=1:sizepop
???
????pop（i:）=5*rands（12）;??
????V（i:）=rands（12）;?
??
????fitness（i）=fun（pop（i:））;?
end


[bestfitness?bestindex]=min（fitness）;
zbest=pop（bestindex:）;??
gbest=pop;??
fitnessgbest=fitness;??
fitnesszbest=bestfitness;??


for?i=1:maxgen
????
????for?j=1:sizepop
????????
?????
????????V（j:）?=?V（j:）?+?c1*rand*（gbest（j:）?-?pop（j:））?+?c2*rand*（zbest?-?pop（j:））;
????????V（jfind（V（j:）>Vmax））=Vmax;
????????V（jf