 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1167 | 回復(fù): 5 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
[求助]
matlab lsqnonlin函數(shù)進(jìn)行非線性動力學(xué)參數(shù)擬合，求高手指點(diǎn). 已有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()估計(jì)得到的參數(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ù)此樓» 猜你喜歡
283分材料與化工求調(diào)劑
已經(jīng)有12人回復(fù)
-
349求調(diào)劑
已經(jīng)有7人回復(fù)
-
296材料專碩求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿上海海洋大學(xué)083200食品學(xué)碩，求調(diào)劑，接受其他專業(yè)
已經(jīng)有7人回復(fù)
-
材料340分調(diào)劑
已經(jīng)有8人回復(fù)
-
316求調(diào)劑
已經(jīng)有13人回復(fù)
-
321求調(diào)劑
已經(jīng)有11人回復(fù)
-
315求調(diào)劑
已經(jīng)有6人回復(fù)
-
330分求調(diào)劑
已經(jīng)有13人回復(fù)
-
化工調(diào)劑303分，過四級
已經(jīng)有9人回復(fù)
3樓2016-01-12 16:41:04
Mr__Right
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +31 - 應(yīng)助: 317 (大學(xué)生)
- 金幣: 14456.3
- 散金: 500
- 紅花: 54
- 帖子: 2716
- 在線: 950.6小時(shí)
- 蟲號: 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
|
這類問題“嘗試”的工作量大。只給一點(diǎn)提示： 1。擬合問題最終數(shù)值上都是“最優(yōu)化”問題，這個(gè)是相對復(fù)雜的非線性非凸優(yōu)化問題 2。這類優(yōu)化問題的求解，常用的迭代方法或基于牛頓法延伸出來的方法，通常無法找到好的最優(yōu)解 理論上，這類問題，只有唯一的全局最優(yōu)解的情況下，沒有數(shù)值算法能夠確保必然收斂到這個(gè)點(diǎn) 現(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
4樓2016-01-13 14:04:11
Mr__Right
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +31 - 應(yīng)助: 317 (大學(xué)生)
- 金幣: 14456.3
- 散金: 500
- 紅花: 54
- 帖子: 2716
- 在線: 950.6小時(shí)
- 蟲號: 1972612
- 注冊: 2012-09-04
- 性別: GG
- 專業(yè): 應(yīng)用數(shù)學(xué)方法
- 管轄: 程序語言
5樓2016-01-13 19:12:34
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿山東大學(xué)，085600，344 +3 | 魏子per 2026-04-02 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 材料調(diào)劑 +6 | 一樣YWY 2026-04-02 | 6/300 |
|
|
[考研] 366求調(diào)劑一志愿東北大學(xué) +7 | 運(yùn)氣來得若有似?/a> 2026-04-02 | 7/350 |
|
|
|
[考研] 265求調(diào)劑11408 +4 | 劉小鹿lu 2026-03-27 | 4/200 |
|
|
|
[考研] 385分 生物學(xué)（071000）求調(diào)劑 +6 | qf626 2026-04-01 | 6/300 |
|
|
|
[考研] 266求調(diào)劑 +10 | 陽陽哇塞 2026-04-01 | 10/500 |
|
|
|
[考研] 求生物學(xué)調(diào)劑 +8 | 15172915737 2026-04-01 | 8/400 |
|
|
|
[考研] 295材料工程專碩求調(diào)劑 +19 | 1428151015 2026-03-27 | 19/950 |
|
|
|
[考研] 材料調(diào)劑 +11 | 一樣YWY 2026-03-31 | 11/550 |
|
|
|
[考研] 085410人工智能 初試316分 求調(diào)劑 +3 | 殘星拂曙 2026-03-31 | 3/150 |
|
|
|
[考研] 080200學(xué)碩，機(jī)械工程專業(yè)277分，求帶走！ +4 | 瓶子PZ 2026-03-31 | 4/200 |
|
|
|
[考研] 375求調(diào)劑 +7 | 雨夏整夜 2026-03-29 | 7/350 |
|
|
|
[考研] 254材料與化工求調(diào)劑 +3 | 翰冬林楠 2026-03-30 | 4/200 |
|
|
|
[考研] 材料與化工272求調(diào)劑 +25 | 阿斯蒂芬2004 2026-03-28 | 25/1250 |
|
|
|
[考研] 334求調(diào)劑 +7 | Trying] 2026-03-31 | 7/350 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +16 | kikiki7 2026-03-30 | 16/800 |
|
|
|
[考研] 085602化工求調(diào)劑（331分） +8 | 111@127 2026-03-30 | 8/400 |
|
|
|
[考研] 0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|
|
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|