谁能给我举一个模拟退火算法MATLAB源代码的简单例子 15

就如题,记住,要简单一点的,就像求某函数的最小值之类的例子... 就如题,记住,要简单一点的,就像求某函数的最小值之类的例子 展开
 我来答
司其玄5u
2013-07-20 · 超过17用户采纳过TA的回答
知道答主
回答量:45
采纳率:0%
帮助的人:42.4万
展开全部
clear
clc
a = 0.95
k = [5;10;13;4;3;11;13;10;8;16;7;4];
k = -k; % 模拟退火算法是求解最小值,故取负数
d = [2;5;18;3;2;5;10;4;11;7;14;6];
restriction = 46;
num = 12;
sol_new = ones(1,num); % 生成初始解
E_current = inf;E_best = inf;
% E_current是当前解对应的目标函数值(即背包中物品总价值);
% E_new是新解的目标函数值;
% E_best是最优解的
sol_current = sol_new; sol_best = sol_new;
t0=97; tf=3; t=t0;
p=1;

while t>=tf
for r=1:100
%产生随机扰动
tmp=ceil(rand.*num);
sol_new(1,tmp)=~sol_new(1,tmp);

%检查是否满足约束
while 1
q=(sol_new*d <= restriction);
if ~q
p=~p; %实现交错着逆转头尾的第一个1
tmp=find(sol_new==1);
if p
sol_new(1,tmp)=0;
else
sol_new(1,tmp(end))=0;
end
else
break
end
end

% 计算背包中的物品价值
E_new=sol_new*k;
if E_new<E_current
E_current=E_new;
sol_current=sol_new;
if E_new<E_best
% 把冷却过程中最好的解保存下来
E_best=E_new;
sol_best=sol_new;
end
else
if rand<exp(-(E_new-E_current)./t)
E_current=E_new;
sol_current=sol_new;
else
sol_new=sol_current;
end
end
end
t=t.*a;
end

disp('最优解为:')
sol_best
disp('物品总价值等于:')
val=-E_best;
disp(val)
disp('背包中物品重量是:')
disp(sol_best * d)
本回答被网友采纳
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
可爱的任刚1314
2017-07-15 · TA获得超过3006个赞
知道小有建树答主
回答量:3482
采纳率:3%
帮助的人:734万
展开全部
th a concert at t
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
学员医学学生
2010-11-12 · TA获得超过7635个赞
知道大有可为答主
回答量:1371
采纳率:0%
帮助的人:2172万
展开全部
先学好C语言 里有很多教程和文章的 不懂问我
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
收起 1条折叠回答
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式