matlab程序实现模拟退火算法程序函数求极值引用后修改 模拟退火算法 matlab

%matlab 程序实现 模拟退火算法程序 函数求极值(引用后修改,感谢 ARMYLAU)

%使用模拟退火法求函数f(x,y) = 3*cos(xy) + x + y2的最小值
%解:根据题意,我们设计冷却表进度表为:
%即初始温度为30
%衰减参数为0.95
%马可夫链长度为10000
%Metropolis的步长为0.02
%结束条件为根据上一个最优解与最新的一个最优解的之差小于某个容差。
%使用METROPOLIS接受准则进行模拟, 程序如下

%* 日期:2012-11-29
%* 作者:steven
%* EMAIL:hxs2004@126.com
%* 结束条件为两次最优解之差小于某小量

function [BestX,BestY]=SimulateAnnealing1
clear;
clc;
%// 要求最优值的目标函数,搜索的最大区间
XMAX= 4;
YMAX = 4;
%冷却表参数
MarkovLength = 10000; %// 马可夫链长度
DecayScale = 0.95; %// 衰减参数
StepFactor = 0.02; %// 步长因子
matlab程序实现模拟退火算法程序函数求极值(引用后修改) 模拟退火算法 matlab
Temperature=30; %// 初始温度
Tolerance = 1e-8; %// 容差
AcceptPoints = 0.0; %// Metropolis过程中总接受点
rnd =rand;

% 随机选点 初值设定
PreX = -XMAX * rand ;
PreY = -YMAX * rand;
PreBestX = PreX;
PreBestY = PreY;

PreX = -XMAX * rand ;
PreY = -YMAX * rand;
BestX = PreX;
BestY = PreY;
% 每迭代一次退火一次(降温), 直到满足迭代条件为止
mm=abs( ObjectFunction( BestX,BestY)-ObjectFunction (PreBestX,PreBestY));
while mm > Tolerance

Temperature=DecayScale*Temperature;
AcceptPoints = 0.0;
% 在当前温度T下迭代loop(即MARKOV链长度)次
for i=0:MarkovLength:1
% 1) 在此点附近随机选下一点
p=0;
while p==0
NextX = PreX + StepFactor*XMAX*(rand-0.5);
NextY = PreY+ StepFactor*YMAX*(rand-0.5);
if p==(~(NextX >= -XMAX && NextX <= XMAX && NextY>= -YMAX && NextY <= YMAX))
p=1;
end
end

% 2) 是否全局最优解
if (ObjectFunction(BestX,BestY) >ObjectFunction(NextX,NextY))
% 保留上一个最优解
PreBestX =BestX;
Pr eBestY = BestY;
% 此为新的最优解
BestX=NextX;
BestY=NextY;
end
% 3) Metropolis过程

if( ObjectFunction(PreX,PreY) - ObjectFunction(NextX,NextY) >0 )
%// 接受, 此处lastPoint即下一个迭代的点以新接受的点开始
PreX=NextX;
PreY=NextY;
AcceptPoints=AcceptPoints+1;
else
changer = -1 * ( ObjectFunction(NextX,NextY) -ObjectFunction(PreX,PreY) ) / Temperature ;
rnd=rand;
p1=exp(changer);
double (p1);

if p1 >rand%// 不接受, 保存原解
PreX=NextX;
PreY=NextY;
AcceptPoints=AcceptPoints+1;
end
end
end
mm=abs( ObjectFunction( BestX,BestY)-ObjectFunction (PreBestX,PreBestY));
end

disp('最小值在点:');
BestX
BestY
disp( '最小值为:{0}');
ObjectFunction(BestX, BestY)
end

****************************************************

子函数,目标函数值计算

function value=ObjectFunction(x,y)
value=3*cos(x*y)+x+y*y;
end

  

爱华网本文地址 » http://www.aihuau.com/a/25101015/256637.html

更多阅读

怎样使用iTunes来安装应用程序 itunes看不到应用程序

iTunes是苹果所有移动设备的PC套件,不管是iPod、iTouch、iPhone还是今天要说的iPad,都要使用iTunes来安装应用程序。 下载完毕后,安装好下载的iTunes,把iPad用数据线连上电脑,iTunes就会识别了。怎样使用iTunes来安装应用程序——步骤/

怎样修改IOS程序名称 修改程序名称

1)iphone上显示的名字修改info.plist生成的app名字修改Target->Get info->build里面对应环境->packing里面的Productname2)这个功能太常用了:程序员经常以一个名字开始编写一个应用;但在发布之前,常常需要为应用更名.或者想从一个

IME消息函数常量 函数返回常量引用

IME消息  以下列出IME中用到的消息。WM_IME_CHAR(IME得到了转换结果中的一个字符)WM_IME_COMPOSITION(IME根据用户击键的情况更改了按键组合状态)WM_IME_COMPOSITIONFULL(IME检测到按键组合窗口的区域无法继续扩展)WM_IME_CONTROL(由应用

在STEP7V5.5中如何实现数值微分?算法 数值微分画线算法

在 STEP 7 V5.5 中如何实现数值微分?【算法】描述在数学中,微分是一个函数变化率的处理。不过在实践中,通常没有数学值,而 (举例来说)只有一个随时间变化的值。这个基于 STEP 7 Basic V5.5的程序根据当前值和前一函数值(一阶导数)计算

声明:《matlab程序实现模拟退火算法程序函数求极值引用后修改 模拟退火算法 matlab》为网友爺全球瞩目分享!如侵犯到您的合法权益请联系我们删除