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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

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

我無恥又無賴

木蟲 (著名寫手)

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

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

小弟在用matlab優(yōu)化參數(shù)的時候遇到了一個問題，我用的是fmincon函數(shù)，但是因為選不到合適的初值所以得到的結(jié)果不是很滿意，所以想著用全局優(yōu)化的算法來試一下，奈何不太會，希望論壇里的大佬們看一下能不能解決這個問題。我直接把程序復(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;%熔點
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ù)此樓

» 猜你喜歡

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

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

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

hzlhm

至尊木蟲 (著名寫手)

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

引用回帖:

8樓: Originally posted by 我無恥又無賴 at 2021-09-17 12:47:54
非常感謝，不過這個全局變量和約束函數(shù)怎么編寫我不是很清楚，老師能否給出代碼，我想學(xué)習(xí)學(xué)習(xí)

...

%目標(biāo)函數(shù)：
function f = myfun(k)
global xexp yexp T w n R
.....
end

%線性約束條件：
function [c,ceq] = mycon(k)
global xexp yexp T w n R
....
c=0.11-sum(abs((yexp-ycal).^2./(yexp)));
end

贊一下

回復(fù)此樓

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

我無恥又無賴

QQ：2120156492

11樓2021-09-17 18:43:53

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

查看全部 13 個回答

曾經(jīng)落伍了

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

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

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

用遺傳算法試試

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

回復(fù)此樓

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

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

hzlhm

至尊木蟲 (著名寫手)

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

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

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

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

贊一下

回復(fù)此樓

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

我無恥又無賴

QQ：2120156492

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

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

我無恥又無賴

木蟲 (著名寫手)

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

送紅花一朵

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

查看全部 13 個回答

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

[考研] 求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以 +3	小丁想進步 2026-03-11	3/150	2026-03-17 11:13 by 出門打醬油是蜜?/a>
[考研] 341求調(diào)劑 +5	搗蛋豬豬 2026-03-11	5/250	2026-03-17 06:56 by Muziwl
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運氣yunqi
[考研] 344求調(diào)劑 +3	knight344 2026-03-16	3/150	2026-03-16 09:42 by 無際的草原
[考研] 326求調(diào)劑 +4	上岸的小葡 2026-03-15	5/250	2026-03-16 08:39 by Linda Hu
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 0703化學(xué)調(diào)劑 +4	快樂的香蕉 2026-03-11	4/200	2026-03-13 22:41 by JourneyLucky
[考研] 0703，333分求調(diào)劑一志愿鄭州大學(xué)-物理化學(xué) +3	李魔女斗篷 2026-03-11	3/150	2026-03-13 22:24 by JourneyLucky
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4	寶貴熱 2026-03-13	4/200	2026-03-13 22:01 by 星空星月
[考研] （081700）化學(xué)工程與技術(shù)-298分求調(diào)劑 +12	11啦啦啦 2026-03-11	35/1750	2026-03-13 21:25 by JourneyLucky
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8	努力奮斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky
[考研] 工科調(diào)劑 +4	Jiang191123！ 2026-03-11	4/200	2026-03-13 15:15 by Miko19
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3	土木小天才 2026-03-12	3/150	2026-03-13 15:01 by JourneyLucky
[考研] 085600材料與化工 309分請求調(diào)劑 +7	dtdxzxx 2026-03-12	8/400	2026-03-13 14:43 by jxchenghu

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

24小時熱門版塊排行榜

我無恥又無賴

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

» 猜你喜歡

hzlhm

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

曾經(jīng)落伍了

hzlhm

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

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

我無恥又無賴

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