应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 化学化工区 » 催化 » 催化反应 » 【求助】催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？
41/1返回列表
查看: 2056  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

kylexy3792

银虫 (小有名气)

[交流] 【求助】催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？

催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？
我按照个基元反应，依据控制速率步骤，得到相应速率表达式，从而求的各个基元反应的速率常数，但是在设定初值后，结果再次迭代最后算的结果，但是把最后的结果当初值，同时改掉其中一个变量，有些变量对最后结果影响很大，有些却几乎不影响其他结果，求高手帮忙解决！

同时求其他更好的解动力学方程组的方法！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2011-03-02 22:43:34
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

kylexy3792

银虫 (小有名气)
jinkai838: 什么反应？能具体么？可以给我发文献，可以联系我 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是函数误差
%%不知道你说的是哪个误差。如果是‘结果误差’就把上句的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人) 回复此楼
2楼2011-03-03 20:00:42
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

kylexy3792

银虫 (小有名气)
没人帮忙解决！
回复此楼
3楼2011-03-05 22:40:16
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zt_chem

木虫 (正式写手)
★ ★
kylexy3792(金币+5): 谢谢，我去试试看，如果可以计算出结果，再给你追加金币！ 2011-03-06 14:02:49
jinkai838(金币+2): 不错！可以具体发帖 2011-03-06 14:42:48
试着用lsqnonlin优化效果是否能好些？
赞一下(1人) 回复此楼
4楼2011-03-06 08:37:35
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kylexy3792 的主题更新
41/1返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 化学调剂0703 +7 啊我我的 2026-03-11 7/350 2026-03-15 23:03 by 凌千颂111
[文学芳草园] 伙伴们，祝我生日快乐吧 +15 myrtle 2026-03-10 24/1200 2026-03-15 21:16 by 苏州_逗号
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 311求调剂 +5 冬十三 2026-03-15 5/250 2026-03-15 18:38 by 无际的草原
[考研] 311求调剂 +3 26研0 2026-03-15 3/150 2026-03-15 09:12 by JourneyLucky
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 一志愿华中农业大学071010，总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-10 3/150 2026-03-14 00:12 by JourneyLucky
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 工科调剂 +4 Jiang191123！ 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
信息提示
关闭
请填处理意见
关闭