求助 哪有问题啊！！！！！

zyli

cities= [18 87 74 71 25 58 4 13 18 24 71 64 68 83 58 54 51 37 41 2 7 22 25 62 87 91 83 41 45 44;54 76 78 71 38 35 50 40 40 40 42 44 60 58 69 69 62 67 84 94 99 64 60 62 32 38 46 26 21 35] ;
t=30+1;   %城市的数目是30个
s=10;   %样本的数目是1400个
  ]1 i! i0 d: S/ U+ V4 Q
x=cities(1,;
y=cities(2,;
3 h1 f! c$ \) w% ]9 e
; x! l. Q3 P9 C% cpop=zeros(s,t);  %得初始的pop矩阵，矩阵的最后一列表示所在行的样本的路径距离
for i=1:s
( }: h3 `* S- ?  F4 P$ `( I   pop(i,1:t-1)=randperm(t-1);   %随机产生1—（t-1）的t-1个打乱的随机数
end
M=100;
for l=1:1:2for i=1:2:s-1
$ T- E2 V' r5 X; F2 I%随机选择一个交叉点
9 J' `6 b* O& {0 O6 p3 b   m=randperm(t-3)+1;
4 p( ?, M8 Q- M9 M   crosspoint=min(m(1),m(2));
%任意两行交叉   
   x1=pop(i,;
   x2=pop(i+1,;
% n* a2 j0 j+ J6 W* J5 W+ h   k=0;
  for K=1:10
+ U# g! P3 }8 w6 }  e" k1 I     while k==0
+ k+ c$ M+ q  h% e       if sqrt((x(x1(crosspoint))-x(x2(crosspoint)))^2+(y(x1(crosspoint))-y(x2(crosspoint)))^2)<M|sqrt((x(x2(crosspoint))-x(x1(crosspoint)))^2+(y(x2(crosspoint))-y(x1(crosspoint)))^2)<M
            pop(i,=[pop(i,1:crosspoint),pop(i+1,crosspoint+1:t)];
            pop(i+1,=[pop(i+1,1:crosspoint),pop(i,crosspoint+1:t)];
6 t2 ~0 q. [9 ^' Q            %检查x1左边的重复性并得到x1的左边
- S0 a/ D9 `3 H! a- H          for j=1:crosspoint
0 o: t( W' B7 u0 c. N* F           while find(x1(crosspoint+1:t)==x1(j))
              zhi=find(x1(crosspoint+1:t)==x1(j)); %确定重复位置
              temp=x2(crosspoint+zhi);
              x1(j)=temp;
           end
          end
         for j=1:crosspoint
: y+ o  U- W1 x+ w+ t. b! U& N8 H3 P           while find(x2(crosspoint+1:t)==x2(j))
% ?1 N3 ^7 u0 K  L7 H- O              zhi=find(x2(crosspoint+1:t)==x2(j)); %确定重复位置
              temp=x1(crosspoint+zhi);
' G$ l, o4 p0 `+ L! Q              x2(j)=temp;
! b% v2 r( E- S+ e' N' P           end
          end
, G- ?. {) Y6 `: n& O( e5 A4 F            k=1;
' m" ]  M  W/ P! @6 Y       else
3 V. Q6 k8 M7 w+ A8 P1 t            m=randperm(t-3)+1;
) {: o( Q+ n, s) ^$ Z            crosspoint=min(m(1),m(2));
3 x& q# N$ _% k0 |% i# [            k=0;
+ Y) v5 A; H) {' A% E8 I       end
- d9 u6 x0 _* z$ H8 V     end
   end
end
endbestL=min(pop(:,t))
5 p; j/ E1 m8 Z% lJ=pop(:,t);
$ u2 H! E# K1 Afi=1./J;

4 v" _* M* x' q* ~/ e4 n! x/ O[Oderfi,Indexfi]=sort(fi);   %对于fi进行排序
BestS=pop(Indexfi(s),;     %得到最短路
  r; }' C5 v+ P! D8 c
$ }, H4 I+ {" k  U9 J/ g) w  B: CI=BestS;
1 j3 X' Y* [1 K, l
for i=1:1:t-1
    x1(i)=x(I(i));
/ \: L6 L' H8 m0 p6 x8 r0 \    y1(i)=y(I(i));
/ F+ s) l% l' K- vend
4 l/ z- ?2 f4 F& ?x1(t)=x(I(1));
$ O- p; A; j* ~$ w* }4 wy1(t)=y(I(1));

' x- C& r+ J7 x8 y; }! t+ X+ }" Dcities_new=[x1;y1];
disp('Best Route is:');disp(cities_new);
pos=[cities_new cities_new(:,1)];

' e6 y1 Z8 l: ?( {2 A7 z" y1 [lentemp=0;
9 L- ?; ^. R* |" w' Dfor i=1:1:t-1
    temp=sqrt((pos(1,i)-pos(1,i+1))^2+(pos(2,i)-pos(2,i+1))^2);
8 ~0 w0 T/ @: e8 k    lentemp=lentemp+temp;
6 u& B& r  X0 _6 H, K4 nend
disp('Shortest Length is:');disp(lentemp);
: d3 K" ^+ ]$ U# ?( A3 t
. e: m& `; h! L% Tfigure(1);
subplot(1,2,1);       %窗口分割的左边部分
x(t)=x(1);y(t)=y(1);
6 I* A) I# c: @7 Nplot(x,y,'-or');
xlabel('X axis'), ylabel('Y axis'), title('原始路径');
* V& i" \8 ^4 O5 Raxis([0,1,0,1]);
* Z3 ^+ E- A& k/ raxis([0,100,0,100]);
axis on
$ P3 i  R7 n4 X  U  V& g2 Mhold on;
subplot(1,2,2);        %窗口分割的右边部分
  e; J0 |* q7 _, S+ Cplot(x1,y1,'-or');
5 _9 U+ y' u+ q1 S/ _( }xlabel('X axis'), ylabel('Y axis'), title('最新的路径');
/ g- T# C: W% y& V9 ~  f  Caxis([0,1,0,1]);
axis([0,100,0,100]);
axis on

- B: Z! w& e4 `! j' w这只是遗传算法的交叉部分，现在有点问题，就是为什么去掉红色字体就可以运行出来，而加上后程序就没反应了呢？
那位高人能找出原因，麻烦看看啊 ，不胜感激啊！