 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 867 | 回復(fù): 3 | |||
| 【懸賞金幣】回答本帖問(wèn)題，作者自制奶酪將贈(zèng)送您 30 個(gè)金幣 | |||
自制奶酪銅蟲(chóng) (小有名氣)
|
[求助]
用matlab求動(dòng)力學(xué)參數(shù)k1-k7,金幣沒(méi)有，太寒磣請(qǐng)見(jiàn)諒 已有1人參與
|
||
|
希望木蟲(chóng)大神幫忙用Matlab解解方程，可以酬金支付， t/h C1 C2 C3 C4 C5 0 100 0 0 0 0 0.5 18.7 48.9 19.3 9.7 0.2 1 6.5 31.3 34.7 19.6 2.6 1.5 2.7 15.3 41.2 28.7 5.7 2 2.1 6.9 40.9 34.5 9.3 3 1.7 2.2 37.2 38.4 17.1 4 1.5 1.8 31.7 38.2 24.8 5 1.2 1.7 26.1 36.3 30.5 6 0.2 0.5 19.4 33.4 36.7 微分方程組 dC1dt = -k1*C1-k7*C1; dC2dt = k1*C1-k2*C2-k5*C2; dC3dt = k2*C2-k3*C3-k6*C6; dC4dt = k3*C3-k4*C4; dC5dt = k4*C4; dC6dt = k5*C2+k6*C3+k7*C1; 我網(wǎng)上依葫蘆畫(huà)瓢填的報(bào)錯(cuò)結(jié)果出不來(lái) 代碼： function odes_fit format long clear all clc k0 = [0 0 0 0 0 0]; lb = -[1 1 1 1 1 1]*1e9; ub = [1 1 1 1 1 1]*1e9; data=... [0 100 0 0 0 0 0 0.5 18.7 48.9 19.3 9.7 0.2 0.0; 1 6.5 31.3 34.7 19.6 2.6 0.0; 1.5 2.7 15.3 41.2 28.7 5.7 0.05; 2 2.1 6.9 40.9 34.5 9.3 0.1; 3 1.7 2.2 37.2 38.4 17.1 0.12; 4 1.5 1.8 31.7 38.2 24.8 0.16; 5 1.2 1.7 26.1 36.3 30.5 0.22; 6 0.2 0.5 19.4 33.4 36.7 0.3; ]; x0=data(1,2:end); tspan = [data(:,1)']; yexp = [data(2:end,2) data(2:end,3) data(2:end,4) data(2:end,5) data(2:end,6) data(2:end,7)]; [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('\tk1 = %.9f \n',k(1)) fprintf('\tk2 = %.9f \n',k(2)) fprintf('\tk3 = %.9f \n',k(3)) fprintf('\tk4 = %.9f \n',k(4)) fprintf('\tk5 = %.9f \n',k(5)) fprintf('\tk5 = %.9f \n',k(6)) fprintf('\tk5 = %.9f \n',k(7)) figure(1) ts=0  (max(tspan)-min(tspan))/100):max(tspan);[ts ys] = ode45(@KineticsEqs,ts,x0,[],k); yy = [data(:,2) data(:,3) data(:,4) data(:,5) data(:,6) data(:,7)]; figure(1) plot(ts,ys(:,1),'b',tspan,yy(:,1),'bo'); figure(2) plot(ts,ys(:,2),'r',tspan,yy(:,2),'ro'); figure(3) plot(ts,ys(:,3),'k',tspan,yy(:,3),'ko'); figure(4) plot(ts,ys(:,4),'g',tspan,yy(:,4),'go'); figure(5) plot(ts,ys(:,5),'m',tspan,yy(:,5),'mo'); figure(6) plot(ts,ys(:,6),'h',tspan,yy(:,6),'ho'); figure(7) plot(ts,ys(:,1),'b',tspan,yy(:,1),'bo',ts,ys(:,2),'r',tspan,yy(:,2),'ro',ts,ys(:,3),'k',tspan,yy(:,3),'ko',ts,ys(:,4),'g',tspan,yy(:,4),'go',ts,ys(:,5),'m',tspan,yy(:,5),'mo',ts,ys(:,5),'h',tspan,yy(:,5),'ho'), legend('C1的計(jì)算值','C1的實(shí)驗(yàn)值','C2的計(jì)算值','C2的實(shí)驗(yàn)值','C3的計(jì)算值','C3的實(shí)驗(yàn)值','C4的計(jì)算值','C4的實(shí)驗(yàn)值','C5的計(jì)算值','C5的實(shí)驗(yàn)值','C6的計(jì)算值','C6的實(shí)驗(yàn)值','Location','best'); 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); Xsim5=Xsim(:,5); Xsim6=Xsim(:,6); ysim(:,1) = Xsim1(2:end); ysim(:,2) = Xsim2(2:end); ysim(:,3) = Xsim3(2:end); ysim(:,4) = Xsim4(2:end); ysim(:,5) = Xsim5(2:end); ysim(:,6) = Xsim6(2:end); f = [(ysim(:,1)-yexp(:,1)) (ysim(:,2)-yexp(:,2)) (ysim(:,3)-yexp(:,3)) (ysim(:,5)-yexp(:,5)) (ysim(:,6)-yexp(:,6)) ... (ysim(:,6)-yexp(:,6))]; function dCdt = KineticsEqs(t,C,k) % ODE模型方程 C1=C(1);C2=C(2);C3=C(3);C4=C(4); k1=k(1);k2=k(2);k3=k(3);k4=k(4);k5=k(5);k6=k(6);k7=k(7); dC1dt = -k1*C1-k7*C1; dC2dt = k1*C1-k2*C2-k5*C2; dC3dt = k2*C2-k3*C3-k6*C6; dC4dt = k3*C3-k4*C4; dC5dt = k4*C4; dC6dt = k5*C2+k6*C3+k7*C1; dCdt = [dC1dt;dC2dt;dC3dt;dC4dt;dC5dt;dC6dt;]; |
回復(fù)此樓» 猜你喜歡
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有2人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有258人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有6人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
幫我的英語(yǔ)口語(yǔ)老師找學(xué)生
已經(jīng)有3人回復(fù)
hzlhm
至尊木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 387 (碩士)
- 金幣: 17807
- 紅花: 53
- 帖子: 2879
- 在線: 606.3小時(shí)
- 蟲(chóng)號(hào): 1517335
- 注冊(cè): 2011-11-30
- 性別: GG
- 專(zhuān)業(yè): 常微分方程與動(dòng)力系統(tǒng)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
獨(dú)孤神宇: 金幣+10, 鼓勵(lì)交流 2021-01-11 20:05:59
獨(dú)孤神宇: 金幣+10, 鼓勵(lì)交流 2021-01-11 20:05:59
|
如方程沒(méi)有錯(cuò)誤的話，其k值為 k1=2.9355； k2=1.318； k3=0.44525； k4=0.2476； k5=2.2315e-14； k6=2.2214e-14； k7=0.47562 C1的擬合精度R^2=0.99696 C2的擬合精度R^2=0.98746 C3的擬合精度R^2=0.74181 C4的擬合精度R^2=0.71926 C5的擬合精度R^2=0.94598 |
2樓2021-01-11 19:19:58
梅子炒米粉
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 127
- 帖子: 52
- 在線: 12.2小時(shí)
- 蟲(chóng)號(hào): 22571169
- 注冊(cè): 2020-06-03
3樓2021-04-24 21:45:22
梅子炒米粉
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 127
- 帖子: 52
- 在線: 12.2小時(shí)
- 蟲(chóng)號(hào): 22571169
- 注冊(cè): 2020-06-03
4樓2021-04-24 21:48:28
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 自制奶酪 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 歡迎211本科同學(xué)，過(guò)A區(qū)國(guó)家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組 +12 | lisimayy 2026-03-04 | 22/1100 |
|
|---|---|---|---|---|
|
|
[考研] A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向 +17 | lisimayy 2026-03-04 | 25/1250 |
|
|
|
[考研] 085602 293分求調(diào)劑 +3 | SivanNano. 2026-03-05 | 3/150 |
|
|
|
[考研] 一志愿武理085601專(zhuān)碩347分求調(diào)劑 +4 | 啊歐歐歐 2026-03-04 | 5/250 |
|
|
|
[考研] 0856材料與化工求調(diào)劑！ +5 | 化工考生111 2026-03-04 | 9/450 |
|
|
|
[考研] 313求調(diào)劑 +3 | 水流年lc 2026-02-28 | 4/200 |
|
|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +7 | xiiiia 2026-03-04 | 7/350 |
|
|
|
[考研] 085600材料調(diào)劑 總分330 +5 | 池池丶 2026-03-03 | 5/250 |
|
|
|
[考研] 0856材料與化工，270求調(diào)劑 +17 | YXCT 2026-03-01 | 20/1000 |
|
|
|
[考研] 材料與化工，291分，求調(diào)劑 +9 | 咕嚕咕嚕123123 2026-03-04 | 11/550 |
|
|
[論文投稿]
 主編終審卡者不動(dòng)
 100+4
|
Stray2021 2026-03-03 | 4/200 |
|
|
|
[考研] 材料專(zhuān)碩346求調(diào)劑 +3 | 旺一下 2026-03-04 | 3/150 |
|
|
|
[考研] 292求調(diào)劑 +9 | yhk_819 2026-02-28 | 9/450 |
|
|
|
[考研] 291求調(diào)劑 +3 | 咕嚕咕嚕123123 2026-03-02 | 4/200 |
|
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
|
[考研] 環(huán)境調(diào)劑 +8 | chenhanheng 2026-03-02 | 8/400 |
|
|
|
[考研] 321求調(diào)劑一志愿東北林業(yè)大學(xué)材料與化工英二數(shù)二 +5 | 蟲(chóng)蟲(chóng)蟲(chóng)蟲(chóng)蟲(chóng)7 2026-03-01 | 9/450 |
|
|
|
[考研] 275求調(diào)劑 +7 | 明遠(yuǎn)求學(xué) 2026-03-01 | 7/350 |
|
|
|
[考研] 哈工大計(jì)算機(jī)劉劼團(tuán)隊(duì)招生 +4 | hit_aiot 2026-03-01 | 6/300 |
|
|
|
[考研] 0856材料求調(diào)劑 +4 | 麻辣魷魚(yú) 2026-02-28 | 4/200 |
|