oor为向负无穷取整函数
bfloorNra
d11
marri1a2
参与交叉的个体标记初值
marri1b3
ifmarri1a1marri1b1ab
marri1a1
marri1b1
参与交叉的个体标记为1
pair
opair
o1
locatio
floor
ra
d11用随机数确定个体中单交叉点位置
l10
l20
forilocatio
step

以下按步长step进行交叉
forj1

用for确定交叉位置
ifgroup1iagroup1jb
l1j
e
d
e
d
forj1

ifgroup1ibgroup1ja
l2j
e
d
e
d
tempgroup1ia
group1iagroup1l2a
group1l2atemp
tempgroup1ib
group1ibgroup1l1b
group1l1btemp
e
d
fe
d
e
d

变异操作
pb01
个体变异率
b
umpbN
变异个体数
fori1b
um
逐个取个体，随机选择位置进行变异
a1floor
ra
d11
a2floor
ra
d11
bfloorNra
d11
tempgroup1a1b
group1a1bgroup1a2b
group1a2btemp
e
d

fori1N
forj1

groupjigroup1ji构造经过复制，交叉，变异后的矩阵，group准备下次循环
e
d
e
d
e
d
disp最优路径为：fori1
fpri
tfdgroupi1ifremi100
fpri
tf
e
de
dfigure1lu
11maxlu
mi
dismi
disiii1lu
plotlu
mi
distitle每一代种群最短距离的收敛过程xlabel遗传代数ylabel每一代种群最短距离
xzeros1
1yzeros1
1fori1
kgroupi1x1icityk1y1icityk2e
dx1
1x11y1
1y11figureplotxygtitle闭合曲线即为最优路径xlabelxylabely
fu
ctio
Ccity_xy
if
10
C00123252584533172571515152525154112elseif
30
C41943784546725627642996858714454628369646018542260834691382538244258697171747887761840134082762325835452141264435450
elseif
50C31323239403037692768375238463162304821472555165717634241173325325648521242738525101745354257322227235637524149495848575839104610591551214828523358276133626320265613132110301536166242636952644367
e
d
fr