%生成初始解,求目标函数f(x)=x1^2+x2^2+8在x1^2-x2>0;-x1-x2^2+2=0约束下的最小值问题
sol_new2=1;%(1)解空间(初始解)
sol_new1=2-sol_new2^2;
sol_current1 = sol_new1;
sol_best1 = sol_new1;
sol_current2 = sol_new2;
sol_best2 = sol_new2;
E_current = inf;
E_best = inf;
rand('state',sum(clock)); %初始化随机数发生器
t=90; %初始温度
tf=89.9; %结束温度
a = 0.99; %温度下降比例
while t>=tf%(7)结束条件
for r=1:1000 %退火次数
%产生随机扰动(3)新解的产生
sol_new2=sol_new2+rand*0.2;
sol_new1=2-sol_new2^2;
%检查是否满足约束
if sol_new1^2-sol_new2>=0 && -sol_new1-sol_new2^2+2==0 && sol_new1>=0 &&sol_new2>=0
else
sol_new2=rand*2;
sol_new1=2-sol_new2^2;
continue;
end
%退火过程
E_new=sol_new1^2+sol_new2^2+8;%(2)目标函数
if E_new<E_current%(5)接受准则
E_current=E_new;
sol_current1=sol_new1;
sol_current2=sol_new2;
if E_new<E_best
%把冷却过程中最好的解保存下来
E_best=E_new;
sol_best1=sol_new1;
sol_best2=sol_new2;
end
else
if rand<exp(-(E_new-E_current)/t)%(4)代价函数差
E_current=E_new;
sol_current1=sol_new1;
sol_current2=sol_new2;
else
sol_new1=sol_current1;
sol_new2=sol_current2;
end
end
plot(r,E_best,'*')
hold on
end
t=t*a;%(6)降温
end
disp('最优解为:')
disp(sol_best1)
disp(sol_best2)
disp('目标表达式的最小值等于:')
disp(E_best)
智能算法之模拟退火算法.rar
版权申诉
74 浏览量
2023-07-31
20:24:42
上传
评论
收藏 841B RAR 举报
小正太浩二
- 粉丝: 195
- 资源: 5919
最新资源
- C程序语言设计 课程设计-万年历程序
- 基于Vue和CSS的泛康云诊所前端设计源码
- stc15xx-STC8-sch-pcb-library-ver1.0 STC单片机的原理图和PCB库文件资料
- 基于Scala的微服务架构后台服务设计源码 - biBACK
- 基于Java的日志推送appender设计源码 - logback-push-appender
- MedievalV1PathTracingBack00000-1920x1080-7e4d63a5fecdc81444c4cacd46d528da.avif
- 基于Vue的医疗消毒检测管理系统设计源码
- 老公-2405282109.awb
- 老公-2405282044.awb
- 8bccc0032ca49b59886db7a4af668bf8.zip
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈