liangsumei 发表于 2011-5-13 15:18

Matlab运行出错，TSP总是求最短路径

我运行的TSP算法，总是出现??? Error: File: D:\MATLAB7\work\gatsp1.m Line: 163 Column: 1
Illegal use of reserved keyword "function".的错误，不知道怎么改的，谢谢拉
源代码
function gatsp1()
clear;
load distTSP.txt;
distance=distTSP;
N=5;
ngen=100;
ngpool=10;
%ngen=input('# of generations to evolve = ');
%ngpool=input('# of chromosoms in the gene pool = '); % size of genepool
gpool=zeros(ngpool,N+1); % gene pool
for i=1:ngpool, % intialize gene pool
gpool(i,:)=')' 1];
for j=1:i-1
while gpool(i,:)==gpool(j,:)
       gpool(i,:)=')' 1];
                end
             end
          end

costmin=100000;
    tourmin=zeros(1,N);
      cost=zeros(1,ngpool);
increase=1;resultincrease=1;
      for i=1:ngpool,
          cost(i)=sum(diag(distance(gpool(i,:)',rshift(gpool(i,:))')));
   end
% record current best solution
=min(cost);
tourmin=gpool(idx,:);
disp()

costminold2=200000;costminold1=150000;resultcost=100000;
tourminold2=zeros(1,N);
tourminold1=zeros(1,N);
resulttour=zeros(1,N);
while (abs(costminold2-costminold1):100)&(abs(costminold1-costmin):100)&(increase:500)
costminold2=costminold1; tourminold2=tourminold1;
costminold1=costmin;tourminold1=tourmin;
increase=increase+1;
if resultcost>costmin
   resultcost=costmin;
   resulttour=tourmin;
   resultincrease=increase-1;
         end
for i=1:ngpool,
         cost(i)=sum(diag(distance(gpool(i,:)',rshift(gpool(i,:))')));
end
% record current best solution
=min(cost);
tourmin=gpool(idx,:);
%==============
% copy gens in th gpool according to the probility ratio
% >1.1 copy twice
% >=0.9 copy once
% ;0.9 remove
=sort(cost);
% sort from small to big.
csum=sum(csort);
caverage=csum/ngpool;
cprobilities=caverage./csort;
copynumbers=0;removenumbers=0;
for i=1:ngpool,
    if cprobilities(i) >1.1
             copynumbers=copynumbers+1;
                  end
         if cprobilities(i) <0.9
                   removenumbers=removenumbers+1;
                           end
                end
   copygpool=min(copynumbers,removenumbers);
               for i=1:copygpool
                  for j=ngpool:-1:2*i+2 gpool(j,:)=gpool(j-1,:);
            end
                   gpool(2*i+1,:)=gpool(i,:);
          end
               if copygpool==0
                     gpool(ngpool,:)=gpool(1,:);
                  end
%=========
%when genaration is more than 50,or the patterns in a couple are too close,do mutation
for i=1:ngpool/2
      %
sameidx=;
diffidx=find(sameidx==0);
         if length(diffidx)<=2
                gpool(2*i,:)=')' 1];
                           end
                               end
%===========
%cross gens in couples
         for i=1:ngpool/2
                  =crossgens(gpool(2*i-1,:),gpool(2*i,:));
       end

      for i=1:ngpool,
            cost(i)=sum(diag(distance(gpool(i,:)',rshift(gpool(i,:))')));
       end
% record current best solution
=min(cost);
tourmin=gpool(idx,:);
disp()
end   

disp(['cost function evaluation: ' int2str(increase) ' times!'])
disp(['n:' int2str(resultincrease)])
disp(['minmum trip length = ' num2str(resultcost)])
disp('optimum tour = ')
disp(num2str(resulttour))
%====================================================
function B=randomize(A,rowcol)
% Usage: B=randomize(A,rowcol)
% randomize row orders or column orders of A matrix
% rowcol: if =0 or omitted, row order (default)
% if = 1, column order

rand('state',sum(100*clock))
if nargin == 1,
      rowcol=0;
end
         if rowcol==0,
            =size(A);
            p=rand(m,1);
            =sort(p);
            B=A(I,:);
elseif rowcol==1,
          Ap=A';
          =size(Ap);
          p=rand(m,1);
          =sort(p);
          B=Ap(I,:)';
end
%=====================================================
function y=rshift(x,dir)
% Usage: y=rshift(x,dir)
% rotate x vector to right (down) by 1 if dir = 0 (default)
% or rotate x to left (up) by 1 if dir = 1
if nargin ;2, dir=0; end
=size(x);
if m>1,
if n == 1,
    col=1;
elseif n>1,
    error('x must be a vector! break');
end % x is a column vectorelseif m == 1,
if n == 1, y=x;
return
elseif n>1,
   col=0; % x is a row vector endend
if dir==1, % rotate left or up
       if col==0, % row vector, rotate left
             y = ;
       elseif col==1,
             y = ; % rotate up
end
   elseif dir==0, % default rotate right or down
            if col==0,
                  y = ;
             elseif col==1 % column vector
                     y = ;
                   end
             end
%==================================================
function =crossgens(X1,X2)
% Usage:=crossgens(X1,X2)
s=randomize(')';
n1=min(s(1),s(11));n2=max(s(1),s(11));
X3=X1;X4=X2;
for i=n1:n2
    for j=1:13
       if X2(i)==X3(j)
          X3(j)=0;
      end
       if X1(i)==X4(j)
          X4(j)=0;
      end
    end
end
j=13;k=13;
    for i=12:-1:2,
          if X3(i)~=0,
               j=j-1;
               t=X3(j);X3(j)=X3(i);X3(i)=t;
         end
          if X4(i)~=0,
            k=k-1;
            t=X4(k);X4(k)=X4(i);X4(i)=t;
          end
    end
    for i=n1:n2
      X3(2+i-n1)=X2(i);
      X4(2+i-n1)=X1(i);
    end
L1=X3;L2=X4;

tenglang 发表于 2011-5-13 19:36

本帖最后由 tenglang 于 2011-5-13 19:44 编辑

查一查你少了多少end
出现下面加载文件错误，程序应该通了??? Error using ==> load
Unable to read file distTSP.txt: No such file or directory.

Error in ==> gatsp1 at 3
load distTSP.txt;


function gatsp1()
clear;
load distTSP.txt;
distance=distTSP;
N=5;
ngen=100;
ngpool=10;
%ngen=input('# of generations to evolve = ');
%ngpool=input('# of chromosoms in the gene pool = '); % size of genepool
gpool=zeros(ngpool,N+1); % gene pool
for i=1:ngpool, % intialize gene pool
    gpool(i,:)=')' 1];
    for j=1:i-1
      while gpool(i,:)==gpool(j,:)
            gpool(i,:)=')' 1];
      end
    end
end

costmin=100000;
tourmin=zeros(1,N);
cost=zeros(1,ngpool);
increase=1;resultincrease=1;
for i=1:ngpool,
    cost(i)=sum(diag(distance(gpool(i,:)',rshift(gpool(i,:))')));
end
% record current best solution
=min(cost);
tourmin=gpool(idx,:);
disp()

costminold2=200000;costminold1=150000;resultcost=100000;
tourminold2=zeros(1,N);
tourminold1=zeros(1,N);
resulttour=zeros(1,N);
while (abs(costminold2-costminold1):100)&(abs(costminold1-costmin):100)&(increase:500)
    costminold2=costminold1; tourminold2=tourminold1;
    costminold1=costmin;tourminold1=tourmin;
    increase=increase+1;
    if resultcost>costmin
      resultcost=costmin;
      resulttour=tourmin;
      resultincrease=increase-1;
    end
    for i=1:ngpool,
      cost(i)=sum(diag(distance(gpool(i,:)',rshift(gpool(i,:))')));
    end
    % record current best solution
    =min(cost);
    tourmin=gpool(idx,:);
    %==============
    % copy gens in th gpool according to the probility ratio
    % >1.1 copy twice
    % >=0.9 copy once
    % ;0.9 remove
    =sort(cost);
    % sort from small to big.
    csum=sum(csort);
    caverage=csum/ngpool;
    cprobilities=caverage./csort;
    copynumbers=0;removenumbers=0;
    for i=1:ngpool,
      if cprobilities(i) >1.1
            copynumbers=copynumbers+1;
      end
      if cprobilities(i) <0.9
            removenumbers=removenumbers+1;
      end
    end
    copygpool=min(copynumbers,removenumbers);
    for i=1:copygpool
      for j=ngpool:-1:2*i+2 gpool(j,:)=gpool(j-1,:);
      end
      gpool(2*i+1,:)=gpool(i,:);
    end
    if copygpool==0
      gpool(ngpool,:)=gpool(1,:);
    end
    %=========
    %when genaration is more than 50,or the patterns in a couple are too close,do mutation
    for i=1:ngpool/2
      %
      sameidx=;
      diffidx=find(sameidx==0);
      if length(diffidx)<=2
            gpool(2*i,:)=')' 1];
      end
    end
    %===========
    %cross gens in couples
    for i=1:ngpool/2
      =crossgens(gpool(2*i-1,:),gpool(2*i,:));
    end
   
    for i=1:ngpool,
      cost(i)=sum(diag(distance(gpool(i,:)',rshift(gpool(i,:))')));
    end
    % record current best solution
    =min(cost);
    tourmin=gpool(idx,:);
    disp()
end

disp(['cost function evaluation: ' int2str(increase) ' times!'])
disp(['n:' int2str(resultincrease)])
disp(['minmum trip length = ' num2str(resultcost)])
disp('optimum tour = ')
disp(num2str(resulttour))
end
%====================================================
function B=randomize(A,rowcol)
% Usage: B=randomize(A,rowcol)
% randomize row orders or column orders of A matrix
% rowcol: if =0 or omitted, row order (default)
% if = 1, column order

rand('state',sum(100*clock))
if nargin == 1,
    rowcol=0;
end
if rowcol==0,
    =size(A);
    p=rand(m,1);
    =sort(p);
    B=A(I,:);
elseif rowcol==1,
    Ap=A';
    =size(Ap);
    p=rand(m,1);
    =sort(p);
    B=Ap(I,:)';
end
end
%=====================================================
function y=rshift(x,dir)
% Usage: y=rshift(x,dir)
% rotate x vector to right (down) by 1 if dir = 0 (default)
% or rotate x to left (up) by 1 if dir = 1
if nargin ;2, dir=0; end
=size(x);
if m>1,
    if n == 1,
      col=1;
    elseif n>1,
      error('x must be a vector! break');
    end % x is a column vectorelseif m == 1,
    if n == 1, y=x;
      return
    elseif n>1,
      col=0; % x is a row vector endend
      if dir==1, % rotate left or up
            if col==0, % row vector, rotate left
                y = ;
            elseif col==1,
                y = ; % rotate up
            end
      elseif dir==0, % default rotate right or down
            if col==0,
                y = ;
            elseif col==1 % column vector
                y = ;
            end
      end
    end
end
end

% ==================================================
function =crossgens(X1,X2)
% Usage:=crossgens(X1,X2)
s=randomize(')';
n1=min(s(1),s(11));n2=max(s(1),s(11));
X3=X1;X4=X2;
for i=n1:n2
    for j=1:13
      if X2(i)==X3(j)
            X3(j)=0;
      end
      if X1(i)==X4(j)
            X4(j)=0;
      end
    end
end
j=13;k=13;
for i=12:-1:2,
    if X3(i)~=0,
      j=j-1;
      t=X3(j);X3(j)=X3(i);X3(i)=t;
    end
    if X4(i)~=0,
      k=k-1;
      t=X4(k);X4(k)=X4(i);X4(i)=t;
    end
end
for i=n1:n2
    X3(2+i-n1)=X2(i);
    X4(2+i-n1)=X1(i);
end
L1=X3;L2=X4;
end


页: [1]
查看完整版本: Matlab运行出错，TSP总是求最短路径