pso-bp代码

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