 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1891 | 回復: 12 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
我無恥又無賴木蟲 (著名寫手)
|
[求助]
matlab全局優(yōu)化參數(shù) 已有2人參與
|
|
小弟在用matlab優(yōu)化參數(shù)的時候遇到了一個問題，我用的是fmincon函數(shù)，但是因為選不到合適的初值所以得到的結果不是很滿意，所以想著用全局優(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 |
回復此樓» 猜你喜歡
德國亥姆霍茲Hereon中心招收兩位醫(yī)用鎂合金腐蝕與LPSO相變方向2026公派博士生
已經(jīng)有0人回復
-
推薦一款可以AI輔助寫作的Latex編輯器SmartLatexEditor，超級好用，推薦試試
已經(jīng)有11人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有79人回復
-
2026-CJ開始申報了
已經(jīng)有1人回復
-
西安電子科學大學杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲能器件方向博士
已經(jīng)有13人回復
-
重慶交通大學光子學微結構與器件課題組招收2026年碩士研究生信息
已經(jīng)有1人回復
-
一志愿鄭大材料學碩298分，求調(diào)劑
已經(jīng)有6人回復
-
尋合作：應力腐蝕多尺度模擬
已經(jīng)有2人回復
hzlhm
至尊木蟲 (著名寫手)
- 應助: 387 (碩士)
- 金幣: 17805.5
- 紅花: 53
- 帖子: 2879
- 在線: 606.3小時
- 蟲號: 1517335
- 注冊: 2011-11-30
- 性別: GG
- 專業(yè): 常微分方程與動力系統(tǒng)
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
獨孤神宇: 金幣+2, 鼓勵交流 2021-09-17 07:49:14
我無恥又無賴: 金幣+100, ★★★很有幫助 2021-09-17 12:49:15
獨孤神宇: 金幣+2, 鼓勵交流 2021-09-17 07:49:14
我無恥又無賴: 金幣+100, ★★★很有幫助 2021-09-17 12:49:15
|
將 f=100*sum(abs((yexp-ycal)./(yexp)))./n 改為 f=100*sum(abs((yexp-ycal).^2./(yexp)))./n; 使用全局變量 global xexp yexp T w n R 增加一個約束條件函數(shù) mycon(k) 123.jpg |
我無恥又無賴
6樓2021-09-16 21:07:22
曾經(jīng)落伍了
鐵桿木蟲 (職業(yè)作家)
- 應助: 0 (幼兒園)
- 金幣: 16291.8
- 散金: 2541
- 紅花: 6
- 帖子: 3321
- 在線: 575.9小時
- 蟲號: 4078713
- 注冊: 2015-09-17
- 性別: GG
- 專業(yè): 機械動力學
2樓2021-09-15 00:45:41
hzlhm
至尊木蟲 (著名寫手)
- 應助: 387 (碩士)
- 金幣: 17805.5
- 紅花: 53
- 帖子: 2879
- 在線: 606.3小時
- 蟲號: 1517335
- 注冊: 2011-11-30
- 性別: GG
- 專業(yè): 常微分方程與動力系統(tǒng)
3樓2021-09-15 22:37:40
我無恥又無賴
木蟲 (著名寫手)
- 應助: 0 (幼兒園)
- 金幣: 1184.9
- 散金: 500
- 帖子: 1402
- 在線: 224小時
- 蟲號: 5857479
- 注冊: 2017-03-06
- 專業(yè): 化工熱力學和基礎數(shù)據(jù)
4樓2021-09-15 23:19:45
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 學碩材料275調(diào)劑 +4 | 路三三 2026-03-03 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 290求調(diào)劑 +9 | ErMiao1020 2026-03-02 | 9/450 |
|
|
|
[碩博家園] 2025屆雙非化工碩士畢業(yè)，申博 +4 | 更多的是 2026-02-27 | 5/250 |
|
|
|
[考研] 283求調(diào)劑 +9 | 鹿沫笙 2026-03-02 | 9/450 |
|
|
|
[考研]
材料270求調(diào)劑
 6+6
|
Eiiiio 2026-03-01 | 11/550 |
|
|
|
[考研] 267求調(diào)劑 +6 | 釣魚佬as 2026-03-02 | 6/300 |
|
|
|
[考研] 中國科學技術大學材料與化工281求調(diào)劑，有科研和獲獎經(jīng)歷 +6 | wsxw 2026-03-02 | 7/350 |
|
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
|
[考研] 材料類考研調(diào)劑 +6 | gemmgemm 2026-03-01 | 7/350 |
|
|
|
[考研] 298求調(diào)劑 +10 | 人間唯你是清歡 2026-02-28 | 14/700 |
|
|
|
[考研] 材料化工調(diào)劑 +13 | 今夏不夏 2026-03-01 | 16/800 |
|
|
|
[考研] 275求調(diào)劑 +7 | 明遠求學 2026-03-01 | 7/350 |
|
|
|
[考博] 26超級電容器申博 +3 | dhdjdjend 2026-02-25 | 3/150 |
|
|
|
[考研] 265分求調(diào)劑不調(diào)專業(yè)和學校有行學上就 +6 | 禮堂丁真258 2026-02-28 | 9/450 |
|
|
|
[基金申請] 此成果不能導入原因：元數(shù)據(jù)必填信息不完整，可 進行補充。 +4 | Kittylucky 2026-03-02 | 5/250 |
|
|
|
[基金申請] 成果系統(tǒng)訪問量大，請一小時后再嘗試。---NSFC啥時候好哦，已經(jīng)兩天這樣了 +4 | NSFC2026我來了 2026-02-28 | 4/200 |
|
|
|
[考研] 291分工科求調(diào)劑 +9 | science餓餓 2026-03-01 | 10/500 |
|
|
|
[考研] 313求調(diào)劑 +3 | 水流年lc 2026-02-28 | 3/150 |
|
|
|
[考研] 調(diào)劑 +3 | 簡木ChuFront 2026-02-28 | 3/150 |
|
|
|
[論文投稿]
求助coordination chemistry reviews 的寫作模板
10+3
|
ljplijiapeng 2026-02-27 | 4/200 |
|