%% 遗传算法优化
popu=20;  
bounds=ones(S,1)*[0,1];
% 产生初始种群
% initPop=crtbp(popu,S);
initPop=randint(popu,S,[0 1]);
% 计算初始种群适应度
initFit=zeros(popu,1);
for i=1:size(initPop,1)
    initFit(i)=de_code(initPop(i,:));
end
initPop=[initPop initFit];
gen=100;
% 优化计算
[X,EndPop,BPop,Trace]=ga(bounds,'fitness',[],initPop,[1e-6 1 0],'maxGenTerm',...
    gen,'normGeomSelect',0.09,'simpleXover',2,'boundaryMutation',[2 gen 3]);
[m,n]=find(X==1);
disp(['优化筛选后的输入自变量编号为:' num2str(n)]);
% 绘制适应度函数进化曲线
figure
plot(Trace(:,1),Trace(:,3),'r:')
hold on
plot(Trace(:,1),Trace(:,2),'b')
xlabel('进化代数')
ylabel('适应度函数')
title('适应度函数进化曲线')
legend('平均适应度函数','最佳适应度函数')
xlim([1 gen])
%% 新训练集/测试集数据提取
p_train=zeros(size(n,2),size(T_train,2));
p_test=zeros(size(n,2),size(T_test,2));
for i=1:length(n)
    p_train(i,:)=P_train(n(i),:);
    p_test(i,:)=P_test(n(i),:);
end
t_train=T_train;
%% 创建优化BP网络
t=cputime;
net_ga=newff(minmax(p_train),[s1,1],{'tansig','purelin'},'trainlm');
% 训练参数设置
net_ga.trainParam.epochs=1000;
net_ga.trainParam.show=10;
net_ga.trainParam.goal=0.1;
net_ga.trainParam.lr=0.1;
net_ga.trainParam.showwindow=0;
%% 训练优化BP网络
net_ga=train(net_ga,p_train,t_train);
%% 仿真测试优化BP网络
tn_ga_sim=sim(net_ga,p_test);
% 反归一化
T_ga_sim=postmnmx(tn_ga_sim,mint,maxt);
e=cputime-t;
T_ga_sim(T_ga_sim>1.5)=2;
T_ga_sim(T_ga_sim<1.5)=1;
result_ga=[T_ga_sim' T_test'];
%% 结果显示（优化BP网络）
number_b_sim=length(find(T_ga_sim==1 & T_test==1));
number_m_sim=length(find(T_ga_sim==2 &T_test==2));
disp('(2)优化BP网络的测试结果为：');
disp(['良性乳腺肿瘤确诊：' num2str(number_b_sim)...
      '  误诊：' num2str(number_B-number_b_sim)...
      '  确诊率p1=' num2str(number_b_sim/number_B*100) '%']);
disp(['恶性乳腺肿瘤确诊：' num2str(number_m_sim)...
      '  误诊：' num2str(number_M-number_m_sim)...
      '  确诊率p2=' num2str(number_m_sim/number_M*100) '%']);
disp(['建模时间为：' num2str(e) 's'] );
以上程序段，在加粗的两个地方提示错误，请教高手如何解决？急急急！