24小时热门版块排行榜    

CyRhmU.jpeg
查看: 340  |  回复: 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 的主题更新
信息提示
请填处理意见