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

liangsumei

我运行的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 randomize([2:N]')' 1];
for j=1:i-1
while gpool(i,:)==gpool(j,:)
       gpool(i,:)=[1 randomize([2:N]')' 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
[costmin,idx]=min(cost);
tourmin=gpool(idx,:);
disp([num2str(increase) 'minmum trip length = ' num2str(costmin)])

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
[costmin,idx]=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
[csort,ridx]=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=[gpool(2*i-1,:)==gpool(2*i,:)];
diffidx=find(sameidx==0);
           if length(diffidx)<=2
                gpool(2*i,:)=[1 randomize([2:12]')' 1];
                           end
                               end
%===========
%cross gens in couples
           for i=1:ngpool/2
                  [gpool(2*i-1,:),gpool(2*i,:)]=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
[costmin,idx]=min(cost);
tourmin=gpool(idx,:);
disp([num2str(increase) 'minmum trip length = ' num2str(costmin)])
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,
              [m,n]=size(A);
              p=rand(m,1);
              [p1,I]=sort(p);
              B=A(I,:);
elseif rowcol==1,
          Ap=A';
          [m,n]=size(Ap);
          p=rand(m,1);
          [p1,I]=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
[m,n]=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 = [x(2:n) x(1)];
       elseif col==1,
             y = [x(2:n); x(1)]; % rotate up
end
   elseif dir==0, % default rotate right or down
              if col==0,
                    y = [x(n) x(1:n-1)];
             elseif col==1 % column vector
                       y = [x(n); x(1:n-1)];
                   end
             end
%==================================================
function [L1,L2]=crossgens(X1,X2)
% Usage:[L1,L2]=crossgens(X1,X2)
s=randomize([2:12]')';
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

查一查你少了多少end
出现下面加载文件错误，程序应该通了

1. ??? Error using ==> load
   
2. Unable to read file distTSP.txt: No such file or directory.
   
3. 
   
4. Error in ==> gatsp1 at 3
   
5. load distTSP.txt;

复制代码

1. function gatsp1()
   
2. clear;
   
3. load distTSP.txt;
   
4. distance=distTSP;
   
5. N=5;
   
6. ngen=100;
   
7. ngpool=10;
   
8. %ngen=input('# of generations to evolve = ');
   
9. %ngpool=input('# of chromosoms in the gene pool = '); % size of genepool
   
10. gpool=zeros(ngpool,N+1); % gene pool
    
11. for i=1:ngpool, % intialize gene pool
    
12.     gpool(i,:)=[1 randomize([2:N]')' 1];
    
13.     for j=1:i-1
    
14.         while gpool(i,:)==gpool(j,:)
    
15.             gpool(i,:)=[1 randomize([2:N]')' 1];
    
16.         end
    
17.     end
    
18. end
    
19. 
    
20. costmin=100000;
    
21. tourmin=zeros(1,N);
    
22. cost=zeros(1,ngpool);
    
23. increase=1;resultincrease=1;
    
24. for i=1:ngpool,
    
25.     cost(i)=sum(diag(distance(gpool(i,:)',rshift(gpool(i,:))')));
    
26. end
    
27. % record current best solution
    
28. [costmin,idx]=min(cost);
    
29. tourmin=gpool(idx,:);
    
30. disp([num2str(increase) 'minmum trip length = ' num2str(costmin)])
    
31. 
    
32. costminold2=200000;costminold1=150000;resultcost=100000;
    
33. tourminold2=zeros(1,N);
    
34. tourminold1=zeros(1,N);
    
35. resulttour=zeros(1,N);
    
36. while (abs(costminold2-costminold1):100)&(abs(costminold1-costmin):100)&(increase:500)
    
37.     costminold2=costminold1; tourminold2=tourminold1;
    
38.     costminold1=costmin;tourminold1=tourmin;
    
39.     increase=increase+1;
    
40.     if resultcost>costmin
    
41.         resultcost=costmin;
    
42.         resulttour=tourmin;
    
43.         resultincrease=increase-1;
    
44.     end
    
45.     for i=1:ngpool,
    
46.         cost(i)=sum(diag(distance(gpool(i,:)',rshift(gpool(i,:))')));
    
47.     end
    
48.     % record current best solution
    
49.     [costmin,idx]=min(cost);
    
50.     tourmin=gpool(idx,:);
    
51.     %==============
    
52.     % copy gens in th gpool according to the probility ratio
    
53.     % >1.1 copy twice
    
54.     % >=0.9 copy once
    
55.     % ;0.9 remove
    
56.     [csort,ridx]=sort(cost);
    
57.     % sort from small to big.
    
58.     csum=sum(csort);
    
59.     caverage=csum/ngpool;
    
60.     cprobilities=caverage./csort;
    
61.     copynumbers=0;removenumbers=0;
    
62.     for i=1:ngpool,
    
63.         if cprobilities(i) >1.1
    
64.             copynumbers=copynumbers+1;
    
65.         end
    
66.         if cprobilities(i) <0.9
    
67.             removenumbers=removenumbers+1;
    
68.         end
    
69.     end
    
70.     copygpool=min(copynumbers,removenumbers);
    
71.     for i=1:copygpool
    
72.         for j=ngpool:-1:2*i+2 gpool(j,:)=gpool(j-1,:);
    
73.         end
    
74.         gpool(2*i+1,:)=gpool(i,:);
    
75.     end
    
76.     if copygpool==0
    
77.         gpool(ngpool,:)=gpool(1,:);
    
78.     end
    
79.     %=========
    
80.     %when genaration is more than 50,or the patterns in a couple are too close,do mutation
    
81.     for i=1:ngpool/2
    
82.         %
    
83.         sameidx=[gpool(2*i-1,:)==gpool(2*i,:)];
    
84.         diffidx=find(sameidx==0);
    
85.         if length(diffidx)<=2
    
86.             gpool(2*i,:)=[1 randomize([2:12]')' 1];
    
87.         end
    
88.     end
    
89.     %===========
    
90.     %cross gens in couples
    
91.     for i=1:ngpool/2
    
92.         [gpool(2*i-1,:),gpool(2*i,:)]=crossgens(gpool(2*i-1,:),gpool(2*i,:));
    
93.     end
    
94.    
    
95.     for i=1:ngpool,
    
96.         cost(i)=sum(diag(distance(gpool(i,:)',rshift(gpool(i,:))')));
    
97.     end
    
98.     % record current best solution
    
99.     [costmin,idx]=min(cost);
    
100.     tourmin=gpool(idx,:);
     
101.     disp([num2str(increase) 'minmum trip length = ' num2str(costmin)])
     
102. end
     
103. 
     
104. disp(['cost function evaluation: ' int2str(increase) ' times!'])
     
105. disp(['n:' int2str(resultincrease)])
     
106. disp(['minmum trip length = ' num2str(resultcost)])
     
107. disp('optimum tour = ')
     
108. disp(num2str(resulttour))
     
109. end
     
110. %====================================================
     
111. function B=randomize(A,rowcol)
     
112. % Usage: B=randomize(A,rowcol)
     
113. % randomize row orders or column orders of A matrix
     
114. % rowcol: if =0 or omitted, row order (default)
     
115. % if = 1, column order
     
116. 
     
117. rand('state',sum(100*clock))
     
118. if nargin == 1,
     
119.     rowcol=0;
     
120. end
     
121. if rowcol==0,
     
122.     [m,n]=size(A);
     
123.     p=rand(m,1);
     
124.     [p1,I]=sort(p);
     
125.     B=A(I,:);
     
126. elseif rowcol==1,
     
127.     Ap=A';
     
128.     [m,n]=size(Ap);
     
129.     p=rand(m,1);
     
130.     [p1,I]=sort(p);
     
131.     B=Ap(I,:)';
     
132. end
     
133. end
     
134. %=====================================================
     
135. function y=rshift(x,dir)
     
136. % Usage: y=rshift(x,dir)
     
137. % rotate x vector to right (down) by 1 if dir = 0 (default)
     
138. % or rotate x to left (up) by 1 if dir = 1
     
139. if nargin ;2, dir=0; end
     
140. [m,n]=size(x);
     
141. if m>1,
     
142.     if n == 1,
     
143.         col=1;
     
144.     elseif n>1,
     
145.         error('x must be a vector! break');
     
146.     end % x is a column vectorelseif m == 1,
     
147.     if n == 1, y=x;
     
148.         return
     
149.     elseif n>1,
     
150.         col=0; % x is a row vector endend
     
151.         if dir==1, % rotate left or up
     
152.             if col==0, % row vector, rotate left
     
153.                 y = [x(2:n) x(1)];
     
154.             elseif col==1,
     
155.                 y = [x(2:n); x(1)]; % rotate up
     
156.             end
     
157.         elseif dir==0, % default rotate right or down
     
158.             if col==0,
     
159.                 y = [x(n) x(1:n-1)];
     
160.             elseif col==1 % column vector
     
161.                 y = [x(n); x(1:n-1)];
     
162.             end
     
163.         end
     
164.     end
     
165. end
     
166. end
     
167. 
     
168. % ==================================================
     
169. function [L1,L2]=crossgens(X1,X2)
     
170. % Usage:[L1,L2]=crossgens(X1,X2)
     
171. s=randomize([2:12]')';
     
172. n1=min(s(1),s(11));n2=max(s(1),s(11));
     
173. X3=X1;X4=X2;
     
174. for i=n1:n2
     
175.     for j=1:13
     
176.         if X2(i)==X3(j)
     
177.             X3(j)=0;
     
178.         end
     
179.         if X1(i)==X4(j)
     
180.             X4(j)=0;
     
181.         end
     
182.     end
     
183. end
     
184. j=13;k=13;
     
185. for i=12:-1:2,
     
186.     if X3(i)~=0,
     
187.         j=j-1;
     
188.         t=X3(j);X3(j)=X3(i);X3(i)=t;
     
189.     end
     
190.     if X4(i)~=0,
     
191.         k=k-1;
     
192.         t=X4(k);X4(k)=X4(i);X4(i)=t;
     
193.     end
     
194. end
     
195. for i=n1:n2
     
196.     X3(2+i-n1)=X2(i);
     
197.     X4(2+i-n1)=X1(i);
     
198. end
     
199. L1=X3;L2=X4;
     
200. end
     

复制代码