 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1163 | 回復(fù): 5 | |||||
[求助]
matlab lsqnonlin函數(shù)進(jìn)行非線性動力學(xué)參數(shù)擬合，求高手指點. 已有1人參與
|
|
function kinetics1 clear all clc a0=[10^11 81 0.3]; lb=[0 0 0];ub=[+inf +inf +inf]; c0=[0.0525564 100]; c16exp=... [0.0525564 100 0.0403533 73.3083 0.0267938 50.3759 0.015037 42.1053 0.0082474 24.812 0.00415801 13.9098 0.00278195 11.6541 0.00185881 10.1504 0.00183435 9.77444 0.0015879 9.77444]; c14exp=... [0.0525564 100 0.0493758 89.8496 0.0441658 81.203 0.0378303 72.9323 0.0357841 67.6692 0.0337337 61.6541 0.0276196 54.8872 0.0239832 47.7444 0.018766 40.6015 0.0160359 35.7143 0.0153455 34.5865]; c12exp=... [0.0526 100 0.0516 98.1203 0.0507 94.7368 0.0484 93.2331 0.0464 89.8496 0.0439 80.4511 0.0429 80.8271 0.0407 77.8195 0.0393 76.3158 0.0381 75.5639 0.0365 72.1805]; [a,resnorm,residual,exitflag,output,lambda,jacobian]=lsqnonlin(@objlm,a0,lb,ub,[],c0,c16exp,c14exp,c12exp); ci=nlparci(a,residual,jacobian); fprintf('\n\n使用函數(shù)lsqnonlin()估計得到的參數(shù)值為：\n'),a function f=objlm(a,c0,c16exp,c14exp,c12exp) %160度 tspan1=[0 0.1005 0.2143 0.3360 0.4365 0.5661 0.6667 0.7593 0.8677 0.9603]; [t,c16]=ode45(@kinetic16,tspan1,c0,[],a); f1=c16(:,1)-c16exp(:,1); f2=c16(:,2)-c16exp(:,2); %140度 tspan2=[0 0.1005 0.1984 0.2857 0.3571 0.4471 0.5529 0.6746 0.8042 0.9074 0.9683]; [t,c14]=ode45(@kinetic14,tspan2,c0,[],a); f3=c14(:,1)-c14exp(:,1); f4=c14(:,2)-c14exp(:,2); %120度 tspan3=[0 0.1005 0.2011 0.2937 0.4021 0.5000 0.5979 0.7090 0.8069 0.8862 0.9894]; [t,c12]=ode45(@kinetic12,tspan3,c0,[],a); f5=c12(:,1)-c12exp(:,1); f6=c12(:,2)-c12exp(:,2); f=[f1;f2;f3;f4;f5;f6]; %構(gòu)造目標(biāo)函數(shù) function dc16dt=kinetic16(t,c16,a) %160度 T1=433.15;R=0.008314; dc16dt=... [-(a(1).*exp(-a(2)./(R.*T1)).*((c16(2)./100).^(a(3))).*c16(1).*0.943) -(7.*a(1).*exp(-a(2)./(R.*T1)).*((c16(2)./100).^(a(3))).*c16(1).*0.943)]; function dc14dt=kinetic14(t,c14,a) %140度 T2=413.15;R=0.008314; dc14dt=... [-(a(1).*exp(-a(2)./(R.*T2)).*((c14(2)./100).^(a(3))).*c14(1).*0.926) -(7.*a(1).*exp(-a(2)./(R.*T2)).*((c14(2)./100).^(a(3))).*c14(1).*0.926)]; function dc12dt=kinetic12(t,c12,a) %120度 T3=393.15;R=0.008314; dc12dt=... [-(a(1).*exp(-a(2)./(R.*T3)).*((c12(2)./100).^(a(3))).*c12(1).*0.907) -(7.*a(1).*exp(-a(2)./(R.*T3)).*((c12(2)./100).^(a(3))).*c12(1).*0.907)]; |
回復(fù)此樓» 猜你喜歡
食品專碩 一志愿雙一流 328
已經(jīng)有4人回復(fù)
-
材料調(diào)劑
已經(jīng)有6人回復(fù)
-
300分，材料，求調(diào)劑，英一數(shù)二
已經(jīng)有5人回復(fù)
-
求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有5人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有13人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0856材料專碩353求調(diào)劑
已經(jīng)有6人回復(fù)
-
07化學(xué)280分求調(diào)劑
已經(jīng)有6人回復(fù)
-
086003食品工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有16人回復(fù)
Mr__Right
專家顧問 (著名寫手)
-
專家經(jīng)驗: +31 - 應(yīng)助: 317 (大學(xué)生)
- 金幣: 14456.3
- 散金: 500
- 紅花: 54
- 帖子: 2716
- 在線: 950.6小時
- 蟲號: 1972612
- 注冊: 2012-09-04
- 性別: GG
- 專業(yè): 應(yīng)用數(shù)學(xué)方法
- 管轄: 程序語言
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
898766282: 金幣+5, ★有幫助 2016-01-12 14:02:40
898766282: 金幣+5, ★有幫助 2016-01-12 16:41:14
898766282: 金幣+10, ★有幫助 2016-01-13 13:59:33
感謝參與，應(yīng)助指數(shù) +1
898766282: 金幣+5, ★有幫助 2016-01-12 14:02:40
898766282: 金幣+5, ★有幫助 2016-01-12 16:41:14
898766282: 金幣+10, ★有幫助 2016-01-13 13:59:33
|
這類問題“嘗試”的工作量大。只給一點提示： 1。擬合問題最終數(shù)值上都是“最優(yōu)化”問題，這個是相對復(fù)雜的非線性非凸優(yōu)化問題 2。這類優(yōu)化問題的求解，常用的迭代方法或基于牛頓法延伸出來的方法，通常無法找到好的最優(yōu)解 理論上，這類問題，只有唯一的全局最優(yōu)解的情況下，沒有數(shù)值算法能夠確保必然收斂到這個點 現(xiàn)在通用的方法是： 先用全局優(yōu)化算法（differential-evolution差分演化，最常用，粒子群PSO）作初步篩查，然后再用nonlincon，fmincon之類函數(shù)自帶的局部優(yōu)化算法 通常要多算幾次做比較，可能還要調(diào)整優(yōu)化的參數(shù)，直到找到最佳結(jié)果 |
2樓2016-01-12 13:03:33
3樓2016-01-12 16:41:04
4樓2016-01-13 14:04:11
Mr__Right
專家顧問 (著名寫手)
-
專家經(jīng)驗: +31 - 應(yīng)助: 317 (大學(xué)生)
- 金幣: 14456.3
- 散金: 500
- 紅花: 54
- 帖子: 2716
- 在線: 950.6小時
- 蟲號: 1972612
- 注冊: 2012-09-04
- 性別: GG
- 專業(yè): 應(yīng)用數(shù)學(xué)方法
- 管轄: 程序語言
5樓2016-01-13 19:12:34
中計量張興賜
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 27.5
- 帖子: 10
- 在線: 22.6小時
- 蟲號: 20222575
- 注冊: 2019-12-21
- 性別: GG
- 專業(yè): 計算機(jī)硬件技術(shù)
|
您好，我現(xiàn)在研一，也在做matlab非線性擬合求解動力學(xué)參數(shù)，當(dāng)待求解的參數(shù)多了之后，求解的結(jié)果總是限入局部最優(yōu)，請問您怎么解決的��？謝謝了！ 發(fā)自小木蟲Android客戶端 |
6樓2019-12-22 18:21:34
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 898766282 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 材料與化工考研調(diào)劑
+7
|
孅華 2026-03-22 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 08工學(xué)調(diào)劑 +11 | 用戶573181 2026-03-20 | 16/800 |
|
|
|
[考研] 求調(diào)劑323材料與化工 +3 | 1124361 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
|
[考研] 收08調(diào)劑生 +6 | komorebi69 2026-03-18 | 6/300 |
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|
|
|
[考研] 321求調(diào)劑 +4 | Ymlll 2026-03-24 | 4/200 |
|
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|
|
[考研] 材料專業(yè)求調(diào)劑 +11 | hanamiko 2026-03-18 | 11/550 |
|
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
|
[考研] 303求調(diào)劑 +4 | 元夕元 2026-03-20 | 4/200 |
|
|
|
[考研] 工科0856求調(diào)劑 +5 | 沐析汀汀 2026-03-21 | 5/250 |
|
|
|
[考研] 333求調(diào)劑 +6 | 87639 2026-03-21 | 10/500 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
|
[考研] 考研調(diào)劑 +3 | 呼呼？~+123456 2026-03-21 | 3/150 |
|
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
 40+4
|
babangida 2026-03-19 | 9/450 |
|
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|