tic;
close?all;
clear;
clc;
format?compact;
%%
train=xlsread（‘C:\Users\XD\Desktop\date2.xlsx‘‘train‘‘A2:G1101‘）;
train_out=xlsread（‘C:\Users\XD\Desktop\date2.xlsx‘‘train‘‘I2:I1101‘）;
test=xlsread（‘C:\Users\XD\Desktop\date2.xlsx‘‘test‘‘A2:G181‘）;
test_out=xlsread（‘C:\Users\XD\Desktop\date2.xlsx‘‘test‘‘I2:I181‘）;
%%歸一化
[train_data?pstrain0]?=?mapminmax（train‘01）;????????
[test_data]?=?mapminmax（‘a(chǎn)pply‘test‘pstrain0）;???????
[train_resultpstrain1]?=?mapminmax（train_out‘01）;???
[test_result]?=?mapminmax（‘a(chǎn)pply‘test_out‘pstrain1）;??

train_data?=?train_data‘;?????????
train_result=train_result‘;
test_data?=?test_data‘;
test_result=test_result‘;
indices?=?crossvalind（‘Kfold‘length（X）3）;%交叉驗(yàn)證函數(shù)，這里的5表示5折，即將length（p_train）個(gè)樣本分成5組。
MAXGEN?=?200;
NIND?=?25;
NVAR?=?2;
PRECI?=?20;
GGAP?=?0.95;
trace?=?zeros（MAXGEN2）;
FieldID?=?...
[rep（[PRECI][1NVAR]）;[00;120200];...
??[11;00;11;11]];
Chrom?=?crtbp（NINDNVAR*PRECI）;
gen?=?1;
BestMSE?=?inf;
Bestgam?=?0;
Bestsig2?=?0;
%%
cg?=?bs2rv（ChromFieldID）;
for?nind?=?1:NIND
????gam=cg（nind1）;
????sig2=cg（nind2）;

????ObjV（nind1）=CVerror1（XYgamsig2indices）;%計(jì)算適應(yīng)值
%?????ObjV（nind1）=CVmse（XYgamsig2）;
end
[BestMSEI]?=?min（ObjV）;
Bestgam?=?cg（I1）;
Bestsig2?=?cg（I2）;

%%
while?gen<=MAXGEN??
????str?=?sprintf（?‘迭代數(shù)為?=?%g‘gen）;
????disp（str）;
????FitnV?=?ranking（ObjV）;%計(jì)算適應(yīng)度
???SelCh=select（‘sus‘ChromFitnVGGAP）;??????????????%選擇
???SelCh=recombin（‘xovsp‘SelCh0.7）;??????????????????%重組
???SelCh=mut（SelCh）;??
????cg?=?bs2rv（SelChFieldID）;
????for?nind?=?1:size（SelCh1）
????????gam=cg（nind1）;
????????sig2=cg（nind2）;
????????ObjVSel（nind1）=CVerror1（XYgamsig2indices）;%計(jì)算子代目標(biāo)函數(shù)值
%?????????ObjVSingle（nind1）=CVmse（XYgamsig2）;
????end
????[ChromObjV]?=?reins（ChromSelCh11ObjVObjVSel）;?%重插入
????[NewBestMSEI]?=?min（ObjV）;
????cg_temp?=?bs2rv（ChromFieldID）;
????temp_NewBestMSE?=?NewBestMSE;
????NewwrostMSE?=?max（ObjV）;
????avgfift=sum（ObjV）/length（ObjV）;
????if?NewBestMSE????????BestMSE?=?NewBestMSE;
???????Bestgam?=?cg