24小时热门版块排行榜    

查看: 434  |  回复: 1
当前主题已经存档。

kky258

木虫 (正式写手)

公鸡爱吃虫

[交流] 【求助】如何在 matlab 中使用模拟退火的方法求全局最优解?

有个优化过程,我原本使用的是单纯形法,但最优解和初值有很大关系,而且很有可能不是全局最优解。
想请高手给我弄个模拟退火求最优解的程序。不胜感激
具体实例如下:

方程为:
x^4+a*x^3+b*x^2+c*x+d=0;                    ------------------ 方程1

其中
a = Aa2 + (Aa1-Aa2)/(1 + exp((T-xa0)/dax))

b = yb0 + (2*Ab/3.14)*(wb/(4*(T-xcb)^2 + wb^2))

c= Ac1*exp(T/tc1) + yc0

d=exp(Ad*T)
-----------------
要解决的问题是

在给定 T 的情况下 ;如何用模拟退火的方法优化a b c d 中的参数值,使得 方程1 的最大解和最小解得平方和最小

a b c d 中的参数是[ A1a A2a xa0 dax yb0 Ab wb xcb Ac1 tc1 yc0 Ad]
(注意其中如果 abcd 的值给的不好的话,方程有可能无解,该如何避免?)


希望有高手能给个程序 ,重谢
谢谢啊
回复此楼
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)

★ ★ ★ ★
kuhailangyu(金币+2,VIP+0):辛苦了!呵呵 7-17 16:26
kky258(金币+2,VIP+0):谢谢啊,课关键是我不知道模拟退火怎么弄啊,程序不会啊 7-17 19:08
方程1的四个解应该可以求得,表示成a,b,c,d的函数,进一步就是[ A1a A2a xa0 dax yb0 Ab wb xcb Ac1 tc1 yc0 Ad]的函数,为了保证方程有解,进一步可以得出a,b,c,d的不等式关系,也就得到了参数[ A1a A2a xa0 dax yb0 Ab wb xcb Ac1 tc1 yc0 Ad]的若干个不等式,然后代入模拟退火的循环里面。
2楼2009-07-17 15:52:08
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kky258 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 22:24 by xiangy672
[基金申请] 国社科又开始会评了,不知道这次命运如何 +5 雨打竹帘 2026-08-30 9/450 2026-08-30 22:23 by 余韵清
[基金申请] 能否申诉? +4 echo8914667 2026-08-30 5/250 2026-08-30 21:11 by echo8914667
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见