[交流] 【求助】催化反應(yīng)動力學(xué)matlab計算各基元反應(yīng)的速率常數(shù)時，該如何避免較小量被忽略？

催化反應(yīng)動力學(xué)matlab計算各基元反應(yīng)的速率常數(shù)時，該如何避免較小量被忽略？
我按照個基元反應(yīng)，依據(jù)控制速率步驟，得到相應(yīng)速率表達(dá)式，從而求的各個基元反應(yīng)的速率常數(shù)，但是在設(shè)定初值后，結(jié)果再次迭代最后算的結(jié)果，但是把最后的結(jié)果當(dāng)初值，同時改掉其中一個變量，有些變量對最后結(jié)果影響很大，有些卻幾乎不影響其他結(jié)果，求高手幫忙解決！

同時求其他更好的解動力學(xué)方程組的方法！

回復(fù)此樓

» 猜你喜歡

0817一志愿蘇州大學(xué)280復(fù)試調(diào)劑已經(jīng)有9人回復(fù)
0817找調(diào)劑已經(jīng)有6人回復(fù)
化學(xué)工程及工業(yè)化學(xué)論文潤色/翻譯怎么收費? 已經(jīng)有190人回復(fù)
PDF在線壓縮軟件已經(jīng)有5人回復(fù)
材料/化學(xué)相關(guān)專業(yè)2026級學(xué)術(shù)/專業(yè)型碩士研究生已經(jīng)有19人回復(fù)
求催化作用基礎(chǔ)第三版（甄開吉) 課件已經(jīng)有3人回復(fù)
26年材料/化學(xué)相關(guān)專業(yè)博士研究生招生溫州大學(xué)化學(xué)與材料工程學(xué)院(26年秋季入學(xué)) 已經(jīng)有0人回復(fù)
求助電催化劑制備，怎么在碳布表面生長PdCo合金已經(jīng)有0人回復(fù)
求助電催化劑制備，怎么在碳布表面生長PdCo合金已經(jīng)有0人回復(fù)
304求調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

Matlab同時擬合活化能，反應(yīng)級數(shù)，速率常數(shù) 已經(jīng)有26人回復(fù)
求助用matlab三次樣條法算一個反應(yīng)的反應(yīng)速率已經(jīng)有13人回復(fù)
反應(yīng)速率的計算方法已經(jīng)有3人回復(fù)
速率方程的問題求教！已經(jīng)有4人回復(fù)
有關(guān)反應(yīng)速率常數(shù)的計算已經(jīng)有9人回復(fù)
酶反應(yīng)動力學(xué)常數(shù) 已經(jīng)有5人回復(fù)
阿倫尼烏斯公式求速率常數(shù)K 已經(jīng)有5人回復(fù)
matlab 擬合反應(yīng)動力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a> 已經(jīng)有14人回復(fù)

請教matlab反應(yīng)動力學(xué)參數(shù)估計遇到的問題，謝謝已經(jīng)有15人回復(fù)
Matlab如何實現(xiàn)這樣的計算：賦值——計算——判斷結(jié)果——重復(fù) 已經(jīng)有10人回復(fù)
Gaussian 計算二級反應(yīng)速率常數(shù) 已經(jīng)有4人回復(fù)
進行催化過程的動力學(xué)研究時怎樣判斷反應(yīng)是否處于內(nèi)外擴散區(qū)，該如何排除他們的影響已經(jīng)有3人回復(fù)
求反應(yīng)動力學(xué)的反應(yīng)熵變數(shù)據(jù) 已經(jīng)有6人回復(fù)
求gaussian計算反應(yīng)速率常數(shù)，指前因子方法。已經(jīng)有25人回復(fù)
求兩個動力學(xué)常數(shù)的計算已經(jīng)有6人回復(fù)
請教：管式積分反應(yīng)器做動力學(xué)時，反應(yīng)速率如何計算？已經(jīng)有10人回復(fù)
詢問固定床上反應(yīng)速率常數(shù)以表觀活化能的具體操作？已經(jīng)有11人回復(fù)
用matlab統(tǒng)計圖片中某種顏色成分占總圖片中所有顏色成分的比率已經(jīng)有4人回復(fù)
求DNA切割動力學(xué)常數(shù)kobs意義已經(jīng)有3人回復(fù)
如何用MS計算反應(yīng)速率常數(shù)k 已經(jīng)有3人回復(fù)
【求助】等溫條件下速率常數(shù)與頻率因子的關(guān)系是相等的嗎？已經(jīng)有7人回復(fù)
【求助】反應(yīng)速率常數(shù)怎么算已經(jīng)有5人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-03-02 22:43:34

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

kylexy3792

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 413.2
帖子: 150
在線: 91.6小時
蟲號: 704366

jinkai838: 什么反應(yīng)？能具體么？可以給我發(fā)文獻(xiàn)，可以聯(lián)系我 2011-03-06 08:15:20

沒有人幫忙看看，解決一下么？

clc;clear;
format long
C=[18000 572.7587379 17427.24126 503.100394 63.00150801 5.429614942 1.227220877 5.34957E-05 31.80663631 28.64788272 2.468940458 0.558038702
18000 442.6621111 17557.33789 391.4982651 46.31421905 3.97596704 0.873659836 2.99887E-05 25.4952724 22.99518737 1.974082883 0.433775459
18000 395.8533613 17604.14664 353.0097887 38.82819644 3.315565494 0.699810656 2.01511E-05 23.85648271 21.55796881 1.840849281 0.388544864
18000 345.603434 17654.39657 311.2417414 31.19218479 2.643727575 0.525780251 1.48443E-05 21.89345442 19.83639301 1.681255082 0.334365283
18000 326.7498543 17673.25015 295.4076551 28.47206709 2.405087713 0.465044362 2.99887E-05 21.11229134 19.15131299 1.617746524 0.312805182
18000 268.2211638 17731.77884 245.416994 20.77111267 1.733166874 0.299890213 2.01511E-05 18.68215469 17.02008314 1.420176413 0.24573341
18000 1386.82292 16613.17708 1331.139206 101.6813129 8.431488388 1.370722461 1.48443E-05 88.51218746 80.81389331 6.701146786 1.089417668
18000 1109.037017 16890.96298 1063.169464 83.66013195 6.960106488 1.169991466 2.99887E-05 72.90022276 66.48318761 5.531070232 0.929770971
18000 557.2888809 17442.71112 532.9063777 44.3791472 3.714372299 0.661563426 2.01511E-05 38.75233847 35.2670053 2.951719351 0.525728013
18000 276.5175885 17723.48241 264.074297 22.62549864 1.899122798 0.347306091 1.48443E-05 19.78050922 17.98334006 1.509472637 0.276047995
];
P_P=C(:,4);
P_O=C(:,5);
P_D=C(:,8);
P_H=C(:,3);
P_A=C(:,7);

W_P=C(:,9);
W_O=C(:,10);
W_D=C(:,11);
W_A=C(:,12);

% W_P=(K_1.*P_P-K_a.*P_O.*P_H)./A;
% W_O=(K_1.*P_P-K_a.*P_O.*P_H-K_2.*P_O+K_b.*P_D.*P_H)./A;
% W_D=(K_2.*P_O-K_b.*P_D.*P_H)./A;

xx=[P_P P_O P_D P_H W_P W_O W_D W_A P_A];

yy=[0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0];

xdata = xx;
ydata = yy;
x0=[9.99E+09 1.410771372584031e+006 1.250036334940072e+010 1.550570624167701e+007 2.97E+10 8.250392906170209e+002 3.468987506591585e+002 2.790032030520837e+002 4.789732864204271e+002  9.528956235795578e+003
]
%options = optimset('TolFun',1e-3);%TolFun是函數(shù)誤差
%%不知道你說的是哪個誤差。如果是‘結(jié)果誤差’就把上句的TolFun給改成TolX x = lsqcurvefit(fun,x0,xdata,ydata,options);
%$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
x=lsqcurvefit(@myfun310,x0',xdata,ydata)
ww10=(x(1).*x(6).*(xdata(:,1)-xdata(:,2).*xdata(:,4)./x(1).*x(2)))./(1+x(6).*xdata(:,1)+x(7).*xdata(:,2)+x(9).*xdata(:,4)+x(8).*xdata(:,3)).^2-xdata(:,5);
ww20=(x(1).*x(6).*(xdata(:,1)-xdata(:,2).*xdata(:,4)./x(1).*x(2)))./(1+x(6).*xdata(:,1)+x(7).*xdata(:,2)+x(9).*xdata(:,4)+x(8).*xdata(:,3)).^2-((x(3).*x(7).*(xdata(:,2)-xdata(:,3).*xdata(:,4)./x(3).*x(4)))./(1+x(6).*xdata(:,1)+x(7).*xdata(:,2)+x(9).*xdata(:,4)+x(8).*xdata(:,3)).^2)-xdata(:,6);
ww30=1*((x(3).*x(7).*(xdata(:,2)-xdata(:,3).*xdata(:,4)./x(3).*x(4)))./(1+x(6).*xdata(:,1)+x(7).*xdata(:,2)+x(9).*xdata(:,4)+x(8).*xdata(:,3)).^2-xdata(:,8)-xdata(:,7));
x=lsqcurvefit(@myfun310,x0',xdata,ydata)%,options)
ww1=(x(1).*x(6).*(xdata(:,1)-xdata(:,2).*xdata(:,4)./x(1).*x(2)))./(1+x(6).*xdata(:,1)+x(7).*xdata(:,2)+x(9).*xdata(:,4)+x(8).*xdata(:,3)).^2-xdata(:,5);
ww2=(x(1).*x(6).*(xdata(:,1)-xdata(:,2).*xdata(:,4)./x(1).*x(2)))./(1+x(6).*xdata(:,1)+x(7).*xdata(:,2)+x(9).*xdata(:,4)+x(8).*xdata(:,3)).^2-((x(3).*x(7).*(xdata(:,2)-xdata(:,3).*xdata(:,4)./x(3).*x(4)))./(1+x(6).*xdata(:,1)+x(7).*xdata(:,2)+x(9).*xdata(:,4)+x(8).*xdata(:,3)).^2)-xdata(:,6);
ww3=1*((x(3).*x(7).*(xdata(:,2)-xdata(:,3).*xdata(:,4)./x(3).*x(4)))./(1+x(6).*xdata(:,1)+x(7).*xdata(:,2)+x(9).*xdata(:,4)+x(8).*xdata(:,3)).^2-xdata(:,8)-xdata(:,7));
if abs(ww1-ww10)>10^(-6) |abs(ww2-ww20)>10^(-6)| abs(ww3-ww30)>10^(-6)
      ww10=ww1;
      ww20=ww2;
      ww30=ww3;
      x=lsqcurvefit(@myfun310,x0',xdata,ydata)%,options)
end
K_1=x(1),K_a=x(2),K_2=x(3),K_b=x(4),K=x(5),K_P=x(6),K_O=x(7),K_D=x(8),K_H=x(9),K_A=x(10)
KK_1=x(1)/x(2)
KK_2=x(3)/x(4)
x0=[x(1) x(2) x(3) x(4) x(5) x(6) x(7) x(8) x(9) x(10)]
%$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$%$$$$$$$$$$
方程為：
function F=myfun310(x,xdata)
format long
F=[(x(1).*x(6).*(xdata(:,1)-xdata(:,2).*xdata(:,4)./x(1).*x(2)))./(1+x(6).*xdata(:,1)+x(7).*xdata(:,2)+x(9).*xdata(:,4)+x(8).*xdata(:,3)).^2-xdata(:,5),(x(1).*x(6).*(xdata(:,1)-xdata(:,2).*xdata(:,4)./x(1).*x(2)))./(1+x(6).*xdata(:,1)+x(7).*xdata(:,2)+x(9).*xdata(:,4)+x(8).*xdata(:,3)).^2-((x(3).*x(7).*(xdata(:,2)-xdata(:,3).*xdata(:,4)./x(3).*x(4)))./(1+x(6).*xdata(:,1)+x(7).*xdata(:,2)+x(9).*xdata(:,4)+x(8).*xdata(:,3)).^2)-xdata(:,6),1*((x(3).*x(7).*(xdata(:,2)-xdata(:,3).*xdata(:,4)./x(3).*x(4)))./(1+x(6).*xdata(:,1)+x(7).*xdata(:,2)+x(9).*xdata(:,4)+x(8).*xdata(:,3)).^2-xdata(:,8)-xdata(:,7))]
有沒有好心人幫忙看看啊？

[ Last edited by kylexy3792 on 2011-3-3 at 20:02 ]

贊一下(1人)

回復(fù)此樓

2樓2011-03-03 20:00:42

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

kylexy3792

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 413.2
帖子: 150
在線: 91.6小時
蟲號: 704366

沒人幫忙解決！

回復(fù)此樓

3樓2011-03-05 22:40:16

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zt_chem

木蟲 (正式寫手)

應(yīng)助: 12 (小學(xué)生)
金幣: 2147.8
帖子: 872
在線: 981.6小時
蟲號: 286920

★ ★
kylexy3792(金幣+5): 謝謝，我去試試看，如果可以計算出結(jié)果，再給你追加金幣！ 2011-03-06 14:02:49
jinkai838(金幣+2): 不錯！可以具體發(fā)帖 2011-03-06 14:42:48

試著用lsqnonlin優(yōu)化效果是否能好些？

贊一下(1人)

回復(fù)此樓

4樓2011-03-06 08:37:35

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 301求調(diào)劑 +6	yy要上岸呀 2026-03-17	6/300	2026-03-17 23:58 by 星空星月
[考研] 268求調(diào)劑 +7	好運連綿不絕 2026-03-12	8/400	2026-03-17 20:28 by xilongliang
[考研] 0703化學(xué)336分求調(diào)劑 +4	zbzihdhd 2026-03-15	5/250	2026-03-17 17:33 by ruiyingmiao
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	5/250	2026-03-17 17:27 by Little-xue
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 工科材料085601 279求調(diào)劑 +3	困于星晨 2026-03-17	3/150	2026-03-17 14:08 by ms629
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +5	薛云鵬 2026-03-13	5/250	2026-03-17 10:15 by Sammy2
[考研] 材料與化工304求B區(qū)調(diào)劑 +7	邱gl 2026-03-11	8/400	2026-03-17 09:36 by 努力學(xué)習(xí)賺彩禮
[考研] 274求調(diào)劑 +5	時間點 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 機械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5	d如愿上岸 2026-03-12	8/400	2026-03-16 15:19 by peike
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 四川大學(xué)085601材料工程專碩初試294求調(diào)劑 +4	祝我們好在冬天 2026-03-11	4/200	2026-03-13 21:39 by peike
[考研] 333求調(diào)劑 +3	球球古力 2026-03-11	3/150	2026-03-13 21:27 by JourneyLucky
[考研] 314求調(diào)劑 +7	無懈可擊的巨人 2026-03-12	7/350	2026-03-13 15:40 by JourneyLucky
[考研] 277求調(diào)劑 +4	anchor17 2026-03-12	4/200	2026-03-13 11:15 by 白夜悠長