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

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

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

回復(fù)此樓

» 猜你喜歡

0817一志愿蘇州大學(xué)280復(fù)試調(diào)劑已經(jīng)有9人回復(fù)
0817找調(diào)劑已經(jīng)有6人回復(fù)
物理化學(xué)論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有66人回復(fù)
PDF在線壓縮軟件已經(jīng)有5人回復(fù)
材料/化學(xué)相關(guān)專業(yè)2026級(jí)學(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)價(jià)值貼推薦，對(duì)您同樣有幫助:

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

請(qǐng)教matlab反應(yīng)動(dòng)力學(xué)參數(shù)估計(jì)遇到的問題，謝謝已經(jīng)有15人回復(fù)
Matlab如何實(shí)現(xiàn)這樣的計(jì)算：賦值——計(jì)算——判斷結(jié)果——重復(fù) 已經(jīng)有10人回復(fù)
Gaussian 計(jì)算二級(jí)反應(yīng)速率常數(shù) 已經(jīng)有4人回復(fù)
進(jìn)行催化過程的動(dòng)力學(xué)研究時(shí)怎樣判斷反應(yīng)是否處于內(nèi)外擴(kuò)散區(qū)，該如何排除他們的影響已經(jīng)有3人回復(fù)
求反應(yīng)動(dòng)力學(xué)的反應(yīng)熵變數(shù)據(jù) 已經(jīng)有6人回復(fù)
求gaussian計(jì)算反應(yīng)速率常數(shù)，指前因子方法。已經(jīng)有25人回復(fù)
求兩個(gè)動(dòng)力學(xué)常數(shù)的計(jì)算已經(jīng)有6人回復(fù)
請(qǐng)教：管式積分反應(yīng)器做動(dòng)力學(xué)時(shí)，反應(yīng)速率如何計(jì)算？已經(jīng)有10人回復(fù)
詢問固定床上反應(yīng)速率常數(shù)以表觀活化能的具體操作？已經(jīng)有11人回復(fù)
用matlab統(tǒng)計(jì)圖片中某種顏色成分占總圖片中所有顏色成分的比率已經(jīng)有4人回復(fù)
求DNA切割動(dòng)力學(xué)常數(shù)kobs意義已經(jīng)有3人回復(fù)
如何用MS計(jì)算反應(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小時(shí)
蟲號(hào): 704366

沒人幫忙解決！

回復(fù)此樓

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

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

查看全部 4 個(gè)回答

kylexy3792

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 413.2
帖子: 150
在線: 91.6小時(shí)
蟲號(hào): 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ù)誤差
%%不知道你說的是哪個(gè)誤差。如果是‘結(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的回帖

zt_chem

木蟲 (正式寫手)

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

★ ★
kylexy3792(金幣+5): 謝謝，我去試試看，如果可以計(jì)算出結(jié)果，再給你追加金幣！ 2011-03-06 14:02:49
jinkai838(金幣+2): 不錯(cuò)！可以具體發(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ā)表

[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 304求調(diào)劑 +8	小熊joy 2026-03-14	8/400	2026-03-17 17:29 by ruiyingmiao
[考研] 302求調(diào)劑 +9	負(fù)心者當(dāng)誅 2026-03-11	9/450	2026-03-17 17:13 by ruiyingmiao
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[基金申請(qǐng)] 國自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 283求調(diào)劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	5/250	2026-03-14 16:30 by JourneyLucky
[考研] 材料371求調(diào)劑 +9	鱷魚? 2026-03-11	11/550	2026-03-13 22:53 by JourneyLucky
[考研] 一志愿中科院，化學(xué)方向，295求調(diào)劑 +4	一氧二氮 2026-03-11	4/200	2026-03-13 22:35 by JourneyLucky
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 【考研調(diào)劑求收留】 +3	Ceciilia 2026-03-11	3/150	2026-03-13 20:18 by JourneyLucky
[考研] 一志愿山大07化學(xué) 332分四六級(jí)已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky
[考研] 277求調(diào)劑 +4	anchor17 2026-03-12	4/200	2026-03-13 11:15 by 白夜悠長