版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖

我無(wú)恥又無(wú)賴(lài)

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 1184.9
散金: 500
帖子: 1402
在線: 224小時(shí)
蟲(chóng)號(hào): 5857479
注冊(cè): 2017-03-06
專(zhuān)業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)

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

小弟在用matlab優(yōu)化參數(shù)的時(shí)候遇到了一個(gè)問(wèn)題，我用的是fmincon函數(shù)，但是因?yàn)檫x不到合適的初值所以得到的結(jié)果不是很滿(mǎn)意，所以想著用全局優(yōu)化的算法來(lái)試一下，奈何不太會(huì)，希望論壇里的大佬們看一下能不能解決這個(gè)問(wèn)題。我直接把程序復(fù)制到下面了

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;%熔點(diǎn)
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

回復(fù)此樓

» 猜你喜歡

【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有93人回復(fù)
幫我的英語(yǔ)口語(yǔ)老師找學(xué)生已經(jīng)有4人回復(fù)
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時(shí)參數(shù)如何選擇？臨近原子數(shù)怎么取值已經(jīng)有7人回復(fù)
2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有0人回復(fù)

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

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

dingd

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

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

引用回帖:

9樓: Originally posted by 我無(wú)恥又無(wú)賴(lài) at 2021-09-17 13:11:22
老師好，您的結(jié)果也是通過(guò)設(shè)置全局變量得到的嗎，我剛才又試了下，第一個(gè)目標(biāo)函數(shù)的值到5點(diǎn)多就降不下來(lái)了...

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

贊一下

回復(fù)此樓

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

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

查看全部 13 個(gè)回答

曾經(jīng)落伍了

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

應(yīng)助: 0 (幼兒園)
金幣: 16291.8
散金: 2541
紅花: 6
帖子: 3321
在線: 575.9小時(shí)
蟲(chóng)號(hào): 4078713
注冊(cè): 2015-09-17
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)

★ ★
獨(dú)孤神宇: 金幣+2, 鼓勵(lì)交流 2021-09-15 08:37:52

用遺傳算法試試

發(fā)自小木蟲(chóng)Android客戶(hù)端

回復(fù)此樓

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

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

hzlhm

至尊木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 387 (碩士)
金幣: 17808.5
紅花: 53
帖子: 2879
在線: 606.3小時(shí)
蟲(chóng)號(hào): 1517335
注冊(cè): 2011-11-30
性別: GG
專(zhuān)業(yè): 常微分方程與動(dòng)力系統(tǒng)

【答案】應(yīng)助回帖

★ ★
感謝參與，應(yīng)助指數(shù) +1
獨(dú)孤神宇: 金幣+2, 鼓勵(lì)交流 2021-09-16 20:51:23

控制誤差的表達(dá)式，
f=100*sum(abs((yexp-ycal)./(yexp)))./n;
可能有問(wèn)題
因?yàn)榭刂普`差的表達(dá)式不同，得到的結(jié)果也不同。

贊一下

回復(fù)此樓

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

我無(wú)恥又無(wú)賴(lài)

QQ：2120156492

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

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

我無(wú)恥又無(wú)賴(lài)

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 1184.9
散金: 500
帖子: 1402
在線: 224小時(shí)
蟲(chóng)號(hào): 5857479
注冊(cè): 2017-03-06
專(zhuān)業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)

送紅花一朵

引用回帖:

3樓: Originally posted by hzlhm at 2021-09-15 22:37:40
控制誤差的表達(dá)式，
f=100*sum(abs((yexp-ycal)./(yexp)))./n;
可能有問(wèn)題
因?yàn)榭刂普`差的表達(dá)式不同，得到的結(jié)果也不同。

主要的目的是讓yexp和ycal更接近，所以用了個(gè)相對(duì)誤差的衡量方法，之前擬合過(guò)挺多數(shù)據(jù)，總的相對(duì)誤差達(dá)到5以下差不多就能用了，現(xiàn)在的主要問(wèn)題就是用fmincon函數(shù)，初值選取的不合適的話(huà)，相對(duì)誤差10以下都很難達(dá)到。

贊一下

回復(fù)此樓

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

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

查看全部 13 個(gè)回答

普通表情龍兔虎貓

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

[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-08	6/300	2026-03-09 11:49 by 勇敢太監(jiān)王公公
[考研] 中科大材料299求調(diào)劑 +10	DAIjiayo 2026-03-05	16/800	2026-03-09 10:45 by 斬魂滴兔子！
[考研] 一志愿東北大學(xué)080500，英一數(shù)二343分求調(diào)劑 +7	Seele_v5 2026-03-04	7/350	2026-03-09 09:55 by EBSD
[考研] 0817化學(xué)工程319求調(diào)劑 +3	lv945 2026-03-08	3/150	2026-03-09 05:23 by houyaoxu
[考研] 286求調(diào)劑 +12	Faune 2026-03-06	14/700	2026-03-08 14:54 by hmn_wj
[考研] 22408 275分求調(diào)劑 +3	宇智波比 2026-03-03	4/200	2026-03-08 12:47 by AaAa7420
[基金申請(qǐng)] 進(jìn)入個(gè)人成果庫(kù)好難，一下午都沒(méi)進(jìn)去 +6	mi_dilee 2026-03-05	6/300	2026-03-08 11:43 by lan111
[考研] 材料工程專(zhuān)碩調(diào)劑自薦信，初試323 +4	李白26 2026-03-07	4/200	2026-03-07 21:59 by 星空星月
[考研] 第一志愿上海大學(xué)，專(zhuān)業(yè)化學(xué)工程與技術(shù)，總分288，求調(diào)劑 +3	1829197082 2026-03-07	3/150	2026-03-07 19:14 by houyaoxu
[考研] 求調(diào)劑 +4	呼呼？~+123456 2026-03-05	5/250	2026-03-06 23:15 by L135790
[考研] 求調(diào)劑一志愿蘇州大學(xué)，0856化工323分 \| 本科應(yīng)化 \| 有專(zhuān)利/競(jìng)賽/科研助手經(jīng)歷 \| +5	橙子cyx 2026-03-06	6/300	2026-03-06 23:13 by L135790
[考研] 085600材料調(diào)劑總分330 +6	池池丶 2026-03-03	6/300	2026-03-06 15:50 by @颯颯颯颯
[考研] 材料與化工304求調(diào)劑 +7	邱gl 2026-03-05	10/500	2026-03-06 01:11 by 卷柏卷柏
[考研] 334求調(diào)劑 +3	Trying] 2026-03-03	4/200	2026-03-05 16:39 by zhukairuo
[考研] 316求調(diào)劑 +3	林小星發(fā)大財(cái) 2026-03-04	3/150	2026-03-05 07:49 by bxbo
[考研] 調(diào)劑材料學(xué)碩 +6	詞凝Y 2026-03-02	6/300	2026-03-05 07:10 by 昆工985
[考研] 293求調(diào)劑 +3	是樂(lè)渝哇 2026-03-04	3/150	2026-03-04 23:12 by wutongshun
[考研] 085601 材料305分求助 +4	泡泡郵件 2026-03-03	6/300	2026-03-04 13:00 by 無(wú)際的草原
[考研] 299求調(diào)劑 +5	kkcoco25 2026-03-02	9/450	2026-03-03 12:55 by 公瑾逍遙
[考研] 清華大學(xué) 材料與化工 353分求調(diào)劑 +5	awaystay 2026-03-02	6/300	2026-03-03 09:03 by houyaoxu

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

24小時(shí)熱門(mén)版塊排行榜

我無(wú)恥又無(wú)賴(lài)

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

» 猜你喜歡

dingd

曾經(jīng)落伍了

hzlhm

【答案】應(yīng)助回帖

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

我無(wú)恥又無(wú)賴(lài)

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