小弟写的源代码如下:
clc;
close;
clear all;
tic
warning off
NIND=35;
MAXGEN=500;
NVAR=14;
PRECI=15;
GGAP=0.7;
trace=zeros(MAXGEN,2);
FieldD=[rep([PRECI],[1,NVAR]);rep([0;1],[1,NVAR]);rep([1;0;1;1],[1,NVAR])];
Chrom=crtbp(NIND,NVAR*PRECI);
gen=0;
x_value=bs2rv(Chrom,FieldD);
a=x_value(1)+x_value(2)+x_value(3)+x_value(4)+x_value(5)+x_value(6)+x_value(7)+x_value(8)+x_value(9)+x_value(10)+x_value(11)+x_value(12)+x_value(13)+x_value(14);
if a==1
objV=entropy(x_value(:,1),x_value(:,2),x_value(:,3),x_value(:,4),x_value(:,5),x_value(:,6),x_value(:,7),x_value(:,8),x_value(:,9),x_value(:,10),x_value(:,11),x_value(:,12),x_value(:,13),x_value(:,14));
else
x_value=bs2rv(Chrom,Fiel