clear?all
close?all
p=textread（‘data.txt‘）;
p1=p‘;
t=p1;
[pnminpmaxptnmintmaxt]?=?premnmx（p1t）;
p_train=[pn（1:）;pn（2:）;pn（3:）;pn（4:）;pn（5:）;pn（6:）;pn（7:）;pn（8:）;pn（9:）;pn（10:）;pn（11:）;pn（12:）;pn（13:）;...
????pn（14:）;pn（15:）;pn（16:）;pn（17:）;pn（18:）;pn（19:）;pn（20:）;pn（21:）;pn（22:）;pn（23:）;pn（24:）;pn（25:）;pn（26:）;pn（27:）;...
????pn（28:）;pn（29:）;pn（30:）;pn（31:）;pn（32:）;pn（33:）;pn（34:）;pn（35:）;pn（36:）;pn（37:）;pn（38:）;pn（39:）]‘
t_train=[pn（2:）;pn（3:）;pn（4:）;pn（5:）;pn（6:）;pn（7:）;pn（8:）;pn（9:）;pn（10:）;pn（11:）;pn（12:）;pn（13:）;pn（14:）;pn（15:）;pn（16:）;...
????pn（17:）;pn（18:）;pn（19:）;pn（20:）;pn（21:）;pn（22:）;pn（23:）;pn（24:）;pn（25:）;pn（26:）;pn（27:）;pn（28:）;pn（29:）;pn（30:）;...
????pn（31:）;pn（32:）;pn（33:）;pn（34:）;pn（35:）;pn（36:）;pn（37:）;pn（38:）;pn（39:）;;pn（40:）]‘;
save?p_train
save?t_train
threshold=minmax（p_train）;?
%創建網絡
net=newff（threshold[122512]{‘tansig‘‘tansig‘‘purelin‘}‘trainlm‘）;
%下面使用遺傳算法對網絡進行優化
P=p_train;
T=t_train;
R=size（P1）;
S2=size（T1）;
S1=25;%隱含層節點數
S=R*S1+S1*S2+S1+S2;%遺傳算法編碼長度
aa=ones（S1）*[-11];
popu=50;%種群規模
initPpp=initializega（pop