版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

我無恥又無賴

木蟲 (著名寫手)

應助: 0 (幼兒園)
金幣: 1184.9
散金: 500
帖子: 1402
在線: 224小時
蟲號: 5857479
注冊: 2017-03-06
專業(yè): 化工熱力學和基礎數據

[求助] matlab全局優(yōu)化參數已有2人參與

小弟在用matlab優(yōu)化參數的時候遇到了一個問題，我用的是fmincon函數，但是因為選不到合適的初值所以得到的結果不是很滿意，所以想著用全局優(yōu)化的算法來試一下，奈何不太會，希望論壇里的大佬們看一下能不能解決這個問題。我直接把程序復制到下面了

clear all
clc
data=[288.05 0.000257 0.3003
293.05 0.000274 0.3003
298.2 0.000281 0.3003
303.25 0.000294 0.3003
308.15 0.000305 0.3003
313.25 0.000326 0.3003
318.2 0.000345 0.3003
323.15 0.00039 0.3003
288.05 0.000493 0.3995
293.05 0.00051 0.3995
298.2 0.000532 0.3995
303.25 0.000562 0.3995
308.15 0.00061 0.3995
313.25 0.000641 0.3995
318.2 0.000698 0.3995
323.15 0.000845 0.3995
288.05 0.000877 0.501
293.05 0.00092 0.501
298.2 0.00099 0.501
303.25 0.00115 0.501
308.15 0.00127 0.501
313.25 0.00139 0.501
318.2 0.00157 0.501
323.15 0.00189 0.501
288.05 0.0016 0.5998
293.05 0.00172 0.5998
298.2 0.00196 0.5998
303.25 0.00232 0.5998
308.15 0.00252 0.5998
313.25 0.00281 0.5998
318.2 0.0032 0.5998
323.15 0.00358 0.5998
288.05 0.00293 0.7009
293.05 0.00309 0.7009
298.2 0.00357 0.7009
303.25 0.00424 0.7009
308.15 0.00474 0.7009
313.25 0.00535 0.7009
318.2 0.00611 0.7009
323.15 0.00711 0.7009
288.05 0.0048 0.8003
293.05 0.00548 0.8003
298.2 0.00638 0.8003
303.25 0.0078 0.8003
308.15 0.00863 0.8003
313.25 0.00988 0.8003
318.2 0.0112 0.8003
323.15 0.0128 0.8003
288.05 0.0103 0.8999
293.05 0.0112 0.8999
298.2 0.0121 0.8999
303.25 0.014 0.8999
308.15 0.0152 0.8999
313.25 0.0173 0.8999
318.2 0.0193 0.8999
323.15 0.022 0.8999];

T=data(:,1);
xexp=data(:,2);
w=data(:,3);

deltaH=46266;%熔融焓
R=8.314;
Tm=532.34;%熔點
yexp=exp((deltaH./(R.*T)).*(T./Tm-1)-log(xexp));

n=size(T,1);
k0=[-3125.70317518255,-2378.98766584565,5092.47952427953,-718.563001024490,-1495.43799248924,-1965.61341412590,0.899991420305235,0.100010845691102,0.219058309329965];
lb=[-inf -inf -inf -inf -inf -inf 0.1 0.1 0.1];
ub=[inf inf inf inf inf inf 0.9 0.9 0.9];
[k,fval,flag] = fmincon(@myfun,k0,[],[],[],[],lb,ub,[],[],xexp,yexp,T,w,n,R);
display(fval);
display(k);

function f = myfun(k,xexp,yexp,T,w,n,R)
%a12=k(1)
%a21=k(2)
%b12=k(3)
%b21=k(4)
%c12=k(5)
%c21=k(6)
%o=k(7)
%p=k(8)
%q=k(9)
tao12=k(1)./(R*T);
tao21=k(2)./(R*T);
tao13=k(3)./(R*T);
tao31=k(4)./(R*T);
tao23=k(5)./(R*T);
tao32=k(6)./(R*T);
G12=exp(-k(7).*tao12);
G21=exp(-k(7).*tao21);
G23=exp(-k(8).*tao23);
G32=exp(-k(8).*tao32);
G13=exp(-k(9).*tao13);
G31=exp(-k(9).*tao31);
x2=(1-xexp).*w;
x3=(1-xexp).*(1-w);
temp1=x2.*G21+x3.*G32;
temp2=tao21.*G21.*x2+tao31.*G31.*x3;
temp3=(xexp+x2.*G21+x3.*G31).^2;
temp4=tao12.*G12.*(x2.^2);
temp5=G12.*G32.*x2.*x3.*(tao12-tao32);
temp6=(x2+xexp.*G12+x3.*G32).^2;
temp7=tao13.*G13.*(x3.^2);
temp8=G13.*G23.*x2.*x3.*(tao13-tao23);
temp9=(x3+x2.*G23+xexp.*G13).^2;
lngamma1=((temp1.*temp2)./temp3)+((temp4+temp5)./temp6)+((temp7+temp8)./temp9);
ycal=exp(lngamma1);
f=100*sum(abs((yexp-ycal)./(yexp)))./n;
end

回復此樓

» 猜你喜歡

2026年循環(huán)經濟功能材料國際會議（ICFMCE 2026）已經有0人回復
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經有0人回復
物理學I論文潤色/翻譯怎么收費? 已經有108人回復
哈爾濱理工大學物理系招收物理學考研調劑已經有5人回復
0702一志愿吉大B區(qū)求調劑已經有5人回復
求調劑已經有0人回復
0702一志愿吉大B區(qū)求調劑有論文已經有0人回復
請問還有沒有用Latex寫文章的小伙伴們？已經有0人回復
光學工程學碩調劑信息已經有24人回復
歡迎加入課題組已經有0人回復

1樓 2021-09-14 21:34:59

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計算強帖: 4
應助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

引用回帖:

9樓: Originally posted by 我無恥又無賴 at 2021-09-17 13:11:22
老師好，您的結果也是通過設置全局變量得到的嗎，我剛才又試了下，第一個目標函數的值到5點多就降不下來了...

用1stOpt計算試試，Matlab的全局優(yōu)化工具箱優(yōu)化效果不行。

贊一下

回復此樓

10樓2021-09-17 17:37:45

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 13 個回答

曾經落伍了

鐵桿木蟲 (職業(yè)作家)

應助: 0 (幼兒園)
金幣: 16291.8
散金: 2541
紅花: 6
帖子: 3321
在線: 575.9小時
蟲號: 4078713
注冊: 2015-09-17
性別: GG
專業(yè): 機械動力學

★ ★
獨孤神宇: 金幣+2, 鼓勵交流 2021-09-15 08:37:52

用遺傳算法試試

發(fā)自小木蟲Android客戶端

回復此樓

2樓2021-09-15 00:45:41

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

hzlhm

至尊木蟲 (著名寫手)

應助: 387 (碩士)
金幣: 17808.5
紅花: 53
帖子: 2879
在線: 606.3小時
蟲號: 1517335
注冊: 2011-11-30
性別: GG
專業(yè): 常微分方程與動力系統(tǒng)

【答案】應助回帖

★ ★
感謝參與，應助指數 +1
獨孤神宇: 金幣+2, 鼓勵交流 2021-09-16 20:51:23

控制誤差的表達式，
f=100*sum(abs((yexp-ycal)./(yexp)))./n;
可能有問題
因為控制誤差的表達式不同，得到的結果也不同。

贊一下

回復此樓

» 本帖已獲得的紅花（最新10朵）

我無恥又無賴

QQ：2120156492

3樓2021-09-15 22:37:40

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

我無恥又無賴

木蟲 (著名寫手)

應助: 0 (幼兒園)
金幣: 1184.9
散金: 500
帖子: 1402
在線: 224小時
蟲號: 5857479
注冊: 2017-03-06
專業(yè): 化工熱力學和基礎數據

送紅花一朵

引用回帖:

3樓: Originally posted by hzlhm at 2021-09-15 22:37:40
控制誤差的表達式，
f=100*sum(abs((yexp-ycal)./(yexp)))./n;
可能有問題
因為控制誤差的表達式不同，得到的結果也不同。

主要的目的是讓yexp和ycal更接近，所以用了個相對誤差的衡量方法，之前擬合過挺多數據，總的相對誤差達到5以下差不多就能用了，現在的主要問題就是用fmincon函數，初值選取的不合適的話，相對誤差10以下都很難達到。

贊一下

回復此樓

4樓2021-09-15 23:19:45

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 13 個回答

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 材料工程專碩274一志愿211求調劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學員h26Tkc
[考研] 302求調劑 +4	小賈同學123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學123
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 考研調劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 0703化學336分求調劑 +3	zbzihdhd 2026-03-15	4/200	2026-03-16 23:52 by zbzihdhd
[考研] 藥學383 求調劑 +3	藥學chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 304求調劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 085600材料與化工求調劑 +13	enenenhui 2026-03-13	14/700	2026-03-16 15:19 by 了了了了。。
[考研] 311求調劑 +6	冬十三 2026-03-15	6/300	2026-03-16 08:00 by wang_dand
[考研] 機械專碩調劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 085601材料工程315分求調劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 復試調劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 求材料調劑 085600英一數二總分302 前三科235 精通機器學習一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 281求調劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 310求調劑 +3	【上上簽】 2026-03-11	3/150	2026-03-13 16:16 by JourneyLucky
[考研] 【0856】化學工程（085602）313 分，本科學科評估A類院�；瘜W工程與工藝，誠求調劑 +7	小劉快快上岸 2026-03-11	7/350	2026-03-13 16:06 by ruiyingmiao
[考研] 307求調劑 +5	超級伊昂大王 2026-03-12	5/250	2026-03-13 15:56 by 棒棒球手
[考研] 295求調劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考研] 274求調劑 +3	S.H1 2026-03-12	3/150	2026-03-13 15:15 by JourneyLucky
[考研] 0817化學工程與技術考研312分調劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

我無恥又無賴

[求助] matlab全局優(yōu)化參數 已有2人參與

» 猜你喜歡

dingd

曾經落伍了

hzlhm

【答案】應助回帖

» 本帖已獲得的紅花（最新10朵）

我無恥又無賴

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] matlab全局優(yōu)化參數已有2人參與