2021-03-13

我用遗传算法优化BP神经网络运行很正常，然后用粒子群算法优化BP神经网络就出现索引超出数组元素的数目(21)。
出错 fun (line 13)
w1=x(1:inputnum*hiddennum);
数据都一样，请问这是为啥
%% 该代码为基于PSO和BP网络的预测
%
%
%

该案例作者申明：

1：本人长期驻扎在此板块里，对该案例提问，做到有问必答。本套书籍官方网站为：video.ourmatlab.com

2：点此从当当预定本书：《Matlab神经网络30个案例分析》。

3：此案例有配套的教学视频，视频下载方式video.ourmatlab.com/vbuy.html。

4：此案例为原创案例，转载请注明出处（《Matlab神经网络30个案例分析》）。

5：若此案例碰巧与您的研究有关联，我们欢迎您提意见，要求等，我们考虑后可以加在案例里。

%
%% 清空环境
clc
clear

%读取数据
load data input output

%节点个数
inputnum=6;
hiddennum=8;
outputnum=1;

%训练数据和预测数据
input_train=input(1:108,:)’;
input_test=input(109:120,:)’;
output_train=output(1:108)’;
output_test=output(109:120)’;

%选连样本输入输出数据归一化
[inputn,inputps]=mapminmax(input_train);
[outputn,outputps]=mapminmax(output_train);

%构建网络
net=newff(inputn,outputn,hiddennum);

% 参数初始化
%粒子群算法中的两个参数
c1 = 1.49445;
c2 = 1.49445;

maxgen=50; % 进化次数
sizepop=30; %种群规模

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

for i=1:sizepop
pop(i,:)=5*rands(1,21);
V(i,:)=rands(1,21);
fitness(i)=fun(pop(i,:),inputnum,hiddennum,outputnum,net,inputn,outputn);
end

% 个体极值和群体极值
[bestfitness bestindex]=min(fitness);
zbest=pop(bestindex,:); %全局最佳
gbest=pop; %个体最佳
fitnessgbest=fitness; %个体最佳适应度值
fitnesszbest=bestfitness; %全局最佳适应度值

%% 迭代寻优
for i=1:maxgen
i;

for j=1:sizepop
    
    %速度更新
    V(j,:) = V(j,:) + c1*rand*(gbest(j,:) - pop(j,:)) + c2*rand*(zbest - pop(j,:));
    V(j,find(V(j,:)>Vmax))=Vmax;
    V(j,find(V(j,:)<Vmin))=Vmin;
    
    %种群更新
    pop(j,:)=pop(j,:)+0.2*V(j,:);
    pop(j,find(pop(j,:)>popmax))=popmax;
    pop(j,find(pop(j,:)<popmin))=popmin;
    
    %自适应变异
    pos=unidrnd(21);
    if rand>0.95
        pop(j,pos)=5*rands(1,1);
    end
  
    %适应度值
    fitness(j)=fun(pop(j,:),inputnum,hiddennum,outputnum,net,inputn,outputn);
end

for j=1:sizepop
%个体最优更新
if fitness(j) < fitnessgbest(j)
    gbest(j,:) = pop(j,:);
    fitnessgbest(j) = fitness(j);
end

%群体最优更新 
if fitness(j) < fitnesszbest
    zbest = pop(j,:);
    fitnesszbest = fitness(j);
end

end

yy(i)=fitnesszbest;  %每代最优值记录到yy数组中  

end

%% 结果分析
plot(yy)
title(['适应度曲线 ’ ‘终止代数＝’ num2str(maxgen)]);
xlabel(‘进化代数’);ylabel(‘适应度’);

x=zbest;
%% 把最优初始阀值权值赋予网络预测
% %用遗传算法优化的BP网络进行值预测
w1=x(1:inputnumhiddennum);
B1=x(inputnumhiddennum+1:inputnumhiddennum+hiddennum);
w2=x(inputnumhiddennum+hiddennum+1:inputnumhiddennum+hiddennum+hiddennumoutputnum);
B2=x(inputnumhiddennum+hiddennum+hiddennumoutputnum+1:inputnumhiddennum+hiddennum+hiddennumoutputnum+outputnum);

net.iw{1,1}=reshape(w1,hiddennum,inputnum);
net.lw{2,1}=reshape(w2,outputnum,hiddennum);
net.b{1}=reshape(B1,hiddennum,1);
net.b{2}=B2;

%% BP网络训练
%网络进化参数
net.trainParam.epochs=100;
net.trainParam.lr=0.1;
%net.trainParam.goal=0.00001;

%网络训练
[net,per2]=train(net,inputn,outputn);

%% BP网络预测
%数据归一化
inputn_test=mapminmax(‘apply’,input_test,inputps);
an=sim(net,inputn_test);
test_simu=mapminmax(‘reverse’,an,outputps);
error=test_simu-output_test;