 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 726 | 回復(fù): 1 | |||
| 【懸賞金幣】回答本帖問題，作者臣本布衣將贈(zèng)送您 150 個(gè)金幣 | |||
臣本布衣新蟲 (初入文壇)
|
[求助]
求助matlab進(jìn)行反應(yīng)動(dòng)力學(xué)擬合，大神幫忙看看代碼哪里有問題 已有1人參與
|
||
|
動(dòng)力學(xué)方程如下： (dC_A)/dt=-(k_1+k_2)*K_A*C_A/(1+K_A*C_A+K_B*C_B+K_C*C_C+K_D*C_D+K_E*C_E) (dC_B)/dt=k_1*K_A*C_A/(1+K_A*C_A+K_B*C_B+K_C*C_C+K_D*C_D+K_E*C_E) (dC_C)/dt=(k_2*K_A*C_A-(k_3+k_4)*K_C*C_C+k_5*K_E*C_E)/(1+K_A*C_A+K_B*C_B+K_C*C_C+K_D*C_D+K_E*C_E) (dC_D)/dt=k_3*K_C*C_C/(1+K_A*C_A+K_B*C_B+K_C*C_C+K_D*C_D+K_E*C_E) (dC_E)/dt=(k_4*K_C*C_C-k_5*K_E*C_E)/(1+K_A*C_A+K_B*C_B+K_C*C_C+K_D*C_D+K_E*C_E) 用matlab語言的ode45解算子求上述微分方程組，以各物質(zhì)濃度計(jì)算值和實(shí)驗(yàn)值的殘差平方和為目標(biāo)函數(shù)，用lsqnonline解算子求解上述最優(yōu)化問題。 代碼如下： function piadatfit clc; clear all; close all; %己烯異構(gòu)化反應(yīng)動(dòng)力學(xué)數(shù)據(jù)擬合％輸入實(shí)驗(yàn)數(shù)據(jù) tspan=[0 30 36 45 60 90 180]; cexp=[0.84691 0 0 0 0; 0.65164 0 0.10305 0.02491 0.00485; 0.65409 0 0.11407 0.02995 0.00681; 0.60309 0 0.14508 0.03510 0.00949; 0.08397 0 0.65655 0.04303 0.01346; 0.05952 0.00752 0.61224 0.05047 0.01814; 0 0.05346 0.53471 0.06414 0.02753]; c0=[0.84691 0 0 0 0]; %輸入?yún)?shù)初值及范圍 k0=[0.1 2 0.05 0.05 0.05 0.3 0.3 0.3 0.1 0.1]; lb=[0 0 0 0 0 0 0 0 0 0]; ub=[10 10 10 10 10 10 10 10 10 10]; %非線性最小二乘擬合 options.algorithm = 'levenberg-marquardt'; [k,resnorm,residual,exitflag,output,lambda,jacobian]=lsqnonlin(@objpia,k0,lb,ub,options,cexp,tspan,c0) %擬合結(jié)果標(biāo)繪 [tplot,cplot]=ode45(@piakin,tspan,c0,[],k); plot(tspan,cexp(:,1),'bx',tplot,cplot(:,1),'b-',tspan,cexp(:,2),'ko',tplot,cplot(:,2),'k-',tspan,cexp(:,3),'g*',tplot,cplot(:,3),'g-',tspan,cexp(:,4),'rs',tplot,cplot(:,4),'r-',tspan,cexp(:,5),'cd',tplot,cplot(:,5),'c-') xlabel('time(min)') ylabel('concentration(mol/l)') legend('caexp','cacal','cbexp','cbcal','ccexp','cccal','cdexp','cdcal','ceexp','cecal') function f=objpia(k0,cexp,tspan,c0) %最小二乘擬合目標(biāo)函數(shù) [t,c]=ode45(@piakin,tspan,c0,[],k0); f1=c(:,1)-cexp(:,1); f2=c(:,2)-cexp(:,2); f=[f1;f2]; function dcdt=piakin(t,c,k) %反應(yīng)動(dòng)力學(xué)方程 r1=k(1)*k(6)*c(1)/(1+k(6)*c(1)+k(7)*c(2)+k(8)*c(3)+k(9)*c(4)+k(10)*c(5)); r2=k(2)*k(6)*c(1)/(1+k(6)*c(1)+k(7)*c(2)+k(8)*c(3)+k(9)*c(4)+k(10)*c(5)); r3=k(3)*k(8)*c(3)/(1+k(6)*c(1)+k(7)*c(2)+k(8)*c(3)+k(9)*c(4)+k(10)*c(5)); r4=k(4)*k(8)*c(3)/(1+k(6)*c(1)+k(7)*c(2)+k(8)*c(3)+k(9)*c(4)+k(10)*c(5)); r5=k(5)*k(10)*c(5)/(1+k(6)*c(1)+k(7)*c(2)+k(8)*c(3)+k(9)*c(4)+k(10)*c(5)); dcdt(1)=-r1-r2; dcdt(2)=r1; dcdt(3)=r2-r3-r4+r5; dcdt(4)=r3; dcdt(5)=r4-r5; dcdt=dcdt'; matlab小白，程序是模仿教材寫的，應(yīng)該存在不少問題，請(qǐng)大神幫忙修改一下！ |
回復(fù)此樓» 猜你喜歡
德國(guó)亥姆霍茲Hereon中心招收兩位醫(yī)用鎂合金腐蝕與LPSO相變方向2026公派博士生
已經(jīng)有0人回復(fù)
-
推薦一款可以AI輔助寫作的Latex編輯器SmartLatexEditor，超級(jí)好用，推薦試試
已經(jīng)有11人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有278人回復(fù)
-
2026-CJ開始申報(bào)了
已經(jīng)有1人回復(fù)
-
西安電子科學(xué)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲(chǔ)能器件方向博士
已經(jīng)有13人回復(fù)
-
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有1人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有3人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
dingd
鐵桿木蟲 (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
微分方程擬合問題，1stOpt試試，更簡(jiǎn)單、方便： Sum Squared Error (SSE): 0.257211769123463 Root of Mean Square Error (RMSE): 0.0925944147205909 Correlation Coef. (R): 0.867645804757432 R-Square: 0.752809242513172 Parameter Best Estimate --------- ------------- k1 0.00075219124176959 k2 0.0093680714188618 k3 0.0883840107805217 k4 9.99999999961713 k5 9.99999999989503 k6 9.99999999999999 k7 4.61483024965692E-14 k8 0.0268675102514994 k9 5.40027468998149E-19 k10 0.307240694154773 |
2樓2022-08-13 12:58:37
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 臣本布衣 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 能動(dòng)297求調(diào)劑，本科川大 +3 | 邵11 2026-03-04 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|
|
|
[考研] 080500材料科學(xué)與工程 +5 | 202114020319 2026-03-03 | 5/250 |
|
|
|
[考研] 求調(diào)劑院校 +6 | 云朵452 2026-03-02 | 11/550 |
|
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
|
[考研] 江蘇省農(nóng)科院招調(diào)劑1名 +5 | Qwertyuop 2026-03-01 | 5/250 |
|
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|
|
|
[考研]
 計(jì)算機(jī)學(xué)碩 分?jǐn)?shù)285求調(diào)劑
+4
|
glwshine 2026-03-02 | 5/250 |
|
|
|
[考研] 267求調(diào)劑 +6 | 釣魚佬as 2026-03-02 | 6/300 |
|
|
|
[考研] 268求調(diào)劑 +6 | 好運(yùn)連綿不絕 2026-03-02 | 6/300 |
|
|
|
[考研] 清華大學(xué) 材料與化工 353分求調(diào)劑 +5 | awaystay 2026-03-02 | 6/300 |
|
|
|
[考研] 298求調(diào)劑 +7 | axyz3 2026-02-28 | 8/400 |
|
|
|
[考研] 【2026 碩士調(diào)劑】課題組 招收調(diào)劑生 +3 | 考研版棒棒 2026-03-02 | 5/250 |
|
|
|
[考研] 285求調(diào)劑 +9 | 滿頭大汗的學(xué)生 2026-02-28 | 9/450 |
|
|
|
[考研] 0856求調(diào)劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
|
[考研] 一志愿山東大學(xué)材料與化工325求調(diào)劑 +5 | 半截的詩(shī)0927 2026-03-02 | 5/250 |
|
|
|
[考研] 一志愿華南理工大學(xué)材料與化工326分，求調(diào)劑 +3 | wujinrui1 2026-02-28 | 3/150 |
|
|
|
[考研] 292求調(diào)劑 +7 | yhk_819 2026-02-28 | 7/350 |
|
|
|
[考研] 284求調(diào)劑 +10 | 天下熯 2026-02-28 | 11/550 |
|