 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (4876)
- >考研 (898)
- >導(dǎo)師招生 (747)
- >文獻(xiàn)求助 (278)
- >蟲(chóng)友互識(shí) (178)
- >考博 (170)
- >博后之家 (127)
- >碩博家園 (106)
- >休閑灌水 (90)
- >基金申請(qǐng) (77)
- >論文投稿 (75)
- >公派出國(guó) (43)
- >招聘信息布告欄 (42)
- >教師之家 (30)
- >找工作 (30)
- >SciFinder/Reaxys (23)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 645 | 回復(fù): 0 | ||
叮咣三腳貓新蟲(chóng) (初入文壇)
|
[求助]
求1stOpt程序代跑，求解微分方程組系數(shù)，可有償，捉急捉急
|
|
需要擬合k1-k8的參數(shù)值 以下是我編寫(xiě)的matlabm文件: function QQQ_4_13_NJ_cs_2 format long clear all clc tspan = [0 0.25 1 3 9 22 35 91 135 ]; x0 = [19750 0 0 0 5250]; k0 = [0.147;0.19;0.03 ;0.02;0.01;0.000001184;3.7;13]; lb = [0;0;0;0;0;0;0;0]; ub = []; data=... [ 0.25 21300 100.225875 46.63575 16.168 3536.970375 1 22050 40.657785 14.123025 8.6688 2874.46475 3 21750 101.43456 48.06252 19.9383 3052.32539 9 19500 279.2715 190.106 23.3415 4908.7372 22 20200 390.3344 403.5 38.5308 3759.5435 35 18550 245.784 262.899 11.2602 5756.7424 91 17950 296.3025 363.12 20.31 6015.16875 135 17500 279.4176 467.418 11.8272 6891.531067 ]; yexp = data(:,2:6); [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,x0,yexp); ci = nlparci(k,residual,jacobian); fprintf('\n\n使用函數(shù)lsqnonlin()估計(jì)得到的參數(shù)值為:\n') fprintf('\tk(1) = %.9f ± %.9f\n',k(1),ci(1,2)-k(1)) fprintf('\tk(2) = %.9f ± %.9f\n',k(2),ci(2,2)-k(2)) fprintf('\tk(3) = %.9f ± %.9f\n',k(3),ci(3,2)-k(3)) fprintf('\tk(4) = %.9f ± %.9f\n',k(4),ci(4,2)-k(4)) fprintf('\tk(5) = %.9f ± %.9f\n',k(5),ci(5,2)-k(5)) fprintf('\tk(6) = %.9f ± %.9f\n',k(6),ci(6,2)-k(6)) fprintf('\tk(7) = %.9f ± %.9f\n',k(7),ci(7,2)-k(7)) fprintf('\tk(8) = %.9f ± %.9f\n',k(8),ci(8,2)-k(8)) fprintf('The sum of the squares is:%.9e\n\n',resnorm) function f = ObjFunc(k,tspan,x0,yexp) % 目標(biāo)函數(shù) [t Xsim] = ode45(@KineticsEqs,tspan,x0,[],k); Xsim1=Xsim(:,1); Xsim2=Xsim(:,2); Xsim3=Xsim(:,3); Xsim4=Xsim(:,4); Xsim4=Xsim(:,5); ysim(:,1) = Xsim1(2:end); ysim(:,2) = Xsim2(2:end); ysim(:,3) = Xsim3(2:end); ysim(:,4) = Xsim4(2:end); ysim(:,5) = Xsim4(2:end); size(ysim(:,1)); size(ysim(:,2)); size(ysim(:,3)); size(ysim(:,4)); size(ysim(:,5)); size(yexp(:,1)); size(yexp(:,2)); size(yexp(:,3)); size(yexp(:,4)); size(yexp(:,5)); f1=ysim(:,1)-yexp(:,1); f2=ysim(:,2)-yexp(:,2); f3=ysim(:,3)-yexp(:,3); f4=ysim(:,4)-yexp(:,4); f5=ysim(:,5)-yexp(:,5); f = [f1;f2;f3;f4;f5]; function dCdt = KineticsEqs(t,C,k) % ODE模型方程 dCAdt =-(k(3)+k(5)+k(7))*C(1)+k(1)*C(2)+k(2)*C(3)+k(2)*C(4)+k(8)*C(5); dCBdt = k(3)*C(1)-(k(1)+k(6))*C(2)+k(4)*C(3)+k(4)*C(4); dCCdt = k(5)*C(1)+k(6)*C(2)-(k(2)+k(4))*C(3); dCDdt = k(5)*C(1)+k(6)*C(2)-(k(2)+k(4))*C(4); dCEdt = k(7)*C(1)-k(8)*C(5); dCdt = [dCAdt;dCBdt;dCCdt;dCDdt;dCEdt]; |
回復(fù)此樓» 猜你喜歡
幫我的英語(yǔ)口語(yǔ)老師找學(xué)生
已經(jīng)有4人回復(fù)
-
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時(shí)參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有7人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有55人回復(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)有2人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有2人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？
已經(jīng)有0人回復(fù)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 叮咣三腳貓 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 紡織、生物、化學(xué)、材料相關(guān)專業(yè)招生了 +3 | 耶耶業(yè) 2026-03-09 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 304求調(diào)劑（085602一志愿985） +9 | 化工人999 2026-03-09 | 9/450 |
|
|
|
[考研] 材料工程085601調(diào)劑求老師收留 +9 | 強(qiáng)木木木 2026-03-07 | 11/550 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-10 | 9/450 |
|
|
|
[考研] 哈工大材料324求調(diào)劑 +6 | 閆旭東 2026-03-10 | 8/400 |
|
|
|
[考研] 0817學(xué)碩華東區(qū)求調(diào)劑 +3 | 30660438 2026-03-08 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 三dd. 2026-03-10 | 3/150 |
|
|
|
[考研] 278求調(diào)劑 +8 | Gale1314 2026-03-06 | 8/400 |
|
|
|
[考博] 求材料讀博院校 +9 | yanglei131 2026-03-08 | 9/450 |
|
|
|
[考研] 一志愿南大化學(xué)339分求調(diào)劑，四六級(jí)已過(guò)，有比賽，有文章 +7 | Gallantzhou 2026-03-07 | 7/350 |
|
|
|
[考研] 337求調(diào)劑 +3 | 睡醒，。 2026-03-09 | 3/150 |
|
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +5 | 加油呀fxy 2026-03-07 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +3 | 拾柒12。 2026-03-08 | 3/150 |
|
|
|
[考研] 347求調(diào)劑 +4 | 浮云滿足 2026-03-07 | 4/200 |
|
|
|
[考研] 材料化工求調(diào)劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
|
[考研] 材料科學(xué)（0805）338 求調(diào)劑 +7 | xiaokang3286 2026-03-07 | 7/350 |
|
|
|
[考研] 一志愿鄭大071000分?jǐn)?shù)282求調(diào)劑 +3 | 研研顏 2026-03-05 | 7/350 |
|
|
|
[考研] 第一志愿上海大學(xué)，專業(yè)化學(xué)工程與技術(shù)，總分288，求調(diào)劑 +3 | 1829197082 2026-03-07 | 3/150 |
|
|
|
[考研] 275求調(diào)劑 +4 | 大爆炸難民 2026-03-06 | 5/250 |
|
|
|
[考博] 2026申博自薦 六級(jí)440電催化方向 +4 | 櫻落成影花成雙 2026-03-05 | 4/200 |
|