Mathematica实现模拟退火法

王茂南

3427
文章

75
评论

2018年1月2日17:05:26

评论3 1931字阅读6分26秒

摘要这一篇文章，我们将会讲一下在mathematica中实现模拟退火法的方式。但是，这一篇不是为了说明模拟退火法是如何实现的，而是为了说明如何编写和使用mathematica的函数包。

这一篇文章，我们将会讲一下在 Mathematica 中实现模拟退火法的方式。其实，在 Mathmatica 中最优化算法里是内置了模拟退火法的，我们可以看下面的例子。

NMinimize[{100 (y - x^2)^2 + (1 - x)^2, -2.084 <= x <= 2.084 && -2.084 <= y <= 2.084}, {x, y}, Method -> "SimulatedAnnealing"]
>> {0,{x->1,y->1}}

在上面指定 Method -> "SimulatedAnnealing" 就是使用模拟退火法。那么，为什么我还要来讲使用 Mathematica 来进行模拟退火法呢，其实我是想利用这个例子来说明一下如何写 Mathematica 的函数包，及 Mathematica 的函数包如何使用，好了，下面开始吧。

哦，还有个事情，就是如果想看 Mathematica 数值非线性最优化全局最优解，可以在帮助文档中搜索：tutorial/ConstrainedOptimizationGlobalNumerical

文章目录(Table of Contents)

模拟退火法例子

模拟退火函数包

BeginPackage["tuihuoLXS`"]
tuihuoLXS::usage = "tuihuoLXS[一组位置坐标]将给出经退火算法优化后的结果,起始概率默认为0.8,可用setInitial设定；停止概率默认为0.2,可通过setEnd设定;每次退火的大小默认为.01，可用setStep设定";
setInitial::usage ="setInitial[初始概率],如setInitial[.8]";
setEnd::usage = "setEnd[停止概率]";
setStep::usage = "setEnd[退火步数]";
evaluationFunction::usage="评估函数,即测试总距离"
Begin["`Private`"]
Ncity;(*城市个数*)
step=.01;
Tstart=.8;(*初始温度*)
Tend=.2;(*停止温度*)
setInitial[n_/;0<n<1]:=Tstart=n;(*设置初始概率*)
setEnd[n_/;0<n<1]:=Tend=n;
setStep[n_/;0<n<1]:=step=n;(*设置步数*)
acceptQ:=RandomChoice[{initialT,1}->{True,False}];(*是否接受当前结果*)
evaluationFunction[list_]:=
Total[EuclideanDistance@@@Partition[list,2,1]];(*计算某条路线的总距离*)
newRoute[list_]:=
ReplacePart[list,{#1->list[[#2]],#2->list[[#1]]}&@@RandomSample[Range[Ncity],2]];
mainfunction[list_]:=
(
routeNew=newRoute[list];
{old,new}=evaluationFunction/@{list,routeNew};
If[new < old, routeNew,If[acceptQ,initialT=initialT-step;routeNew,list]]
);(*迭代函数*)
tuihuoLXS[initialList_]/;MatrixQ[initialList]&&And@@(NumberQ/@Flatten[initialList]):=
(
initialT = (Tstart&)[];
Ncity = Length@initialList;
Nlist = N[initialList];
NestWhile[mainfunction,Nlist,(#;initialT>Tend)&]
);(*主函数*)
End[]
EndPackage[]