模拟退火matlab程序错误请求支援！

概率论 · 发表于 2005-6-27 18:29:05

<

>各位大侠：
<

>我这里有个模拟退火的简单易懂的matlab程序，运行时总是说函数Annealing的第12行有问题，我觉得程序没有任何问题，不知道问题出在哪，请各位高手指点。
<

>主程序如下：
%%%%%%%模拟退火tsp主程序%%%%%%%%%%%%%%%% % 选取初始点并初始化变量 Location=[0.37 0.75 0.45 0.76 0.71 0.07 0.42 0.59 0.32 0.60 0.30 0.67 0.62 0.67 0.20 0.35 0.27 0.94 0.82 0.37 0.61 0.42 0.60 0.39 0.53 0.40 0.63 0.50 0.98 0.68; 0.91 0.87 0.85 0.75 0.72 0.74 0.71 0.69 0.64 0.64 0.59 0.59 0.55 0.55 0.50 0.45 0.43 0.42 0.38 0.27 0.26 0.25 0.23 0.19 0.19 0.13 0.08 0.04 0.02 0.85]'; NumCity=size(Location,1); for I=1:NumCity for J=1:NumCity Distance(I,J)=norm(Location(I,

-Location(J,

); end end Path=randperm(NumCity); Energy=sum(Distance((Path-1)*NumCity+[Path(2:NumCity) Path(1)]));
%global MaxTrialN,NumCity MaxTrialN=30; NumCity=length(Location); StopTolerance=0.01; TempRatio=0.5; Temp=70; MinE=inf; MaxE=-1; while (MaxE-MinE)/MaxE>StopTolerance, [Energy,MaxE,MinE,Path]=Annealing(Path,Energy,Temp,MaxE,MinE); Temp=Temp*TempRatio; end

模拟退火函数如下：
%固定温度的模拟退火子函数 function [Energy,MaxE,MinE,Path]=Annealing(Path,Energy,Temp,MaxE,MinE) TrialN=0; MaxTrialN=30; NumCity=30; while TrialN<MaxTrialN NewPath=Path; Index=ceil(rand(2,1)*NumCity); Temp1=NewPath(Index'); NewPath(Index(1,1))=Temp1(1,2); NewPath(Index(2,1))=Temp1(1,1); NewEnergy=sum(Distance((NewPath-1)*NumCity+[NewPath(2:NumCity) NewPath(1)])); if rand<1/(1+exp((NewEnergy-Energy)/Temp)) Energy=NewEnergy; Path=NewPath; MinE=min(MinE,Energy); MaxE=max(MaxE,Energy); end TrialN=TrialN+1; end

		自动登录	找回密码
密码			注-册-帐-号