| 查看: 2059 | 回复: 3 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？
|
|||
|
催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？ 我按照个基元反应，依据控制速率步骤，得到相应速率表达式，从而求的各个基元反应的速率常数，但是在设定初值后，结果再次迭代最后算的结果，但是把最后的结果当初值，同时改掉其中一个变量，有些变量对最后结果影响很大，有些却几乎不影响其他结果，求高手帮忙解决！  同时求其他更好的解动力学方程组的方法！ |
回复此楼» 猜你喜欢
0817一志愿苏州大学280复试调剂
已经有9人回复
-
0817找调剂
已经有6人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有90人回复
-
PDF在线压缩软件
已经有5人回复
-
材料/化学相关专业2026级学术/专业型硕士研究生
已经有19人回复
-
求催化作用基础第三版 （甄开吉) 课件
已经有3人回复
-
26年材料/化学相关专业博士研究生招生 温州大学化学与材料工程学院(26年秋季入学)
已经有0人回复
-
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有0人回复
-
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有0人回复
-
304求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Matlab同时拟合活化能，反应级数，速率常数
已经有26人回复
- 求助用matlab三次样条法算一个反应的反应速率
已经有13人回复
- 反应速率的计算方法
已经有3人回复
- 速率方程的问题求教！
已经有4人回复
- 有关反应速率常数的计算
已经有9人回复
- 酶反应动力学常数
已经有5人回复
- 阿伦尼乌斯公式求速率常数K
已经有5人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- 请教matlab反应动力学参数估计遇到的问题，谢谢
已经有15人回复
- Matlab如何实现这样的计算：赋值——计算——判断结果——重复
已经有10人回复
- Gaussian 计算二级反应速率常数
已经有4人回复
- 进行催化过程的动力学研究时怎样判断反应是否处于内外扩散区，该如何排除他们的影响
已经有3人回复
- 求反应动力学的反应熵变数据
已经有6人回复
- 求gaussian计算反应速率常数，指前因子方法。
已经有25人回复
- 求两个动力学常数的计算
已经有6人回复
- 请教：管式积分反应器做动力学时，反应速率如何计算？
已经有10人回复
- 询问固定床上反应速率常数以表观活化能的具体操作？
已经有11人回复
- 用matlab统计图片中某种颜色成分占总图片中所有颜色成分的比率
已经有4人回复
- 求DNA切割动力学常数kobs意义
已经有3人回复
- 如何用MS计算反应速率常数k
已经有3人回复
- 【求助】等温条件下速率常数与频率因子的关系是相等的吗？
已经有7人回复
- 【求助】反应速率常数怎么算
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
坐标上海，诚征女友，非常 着急，私信必回
+1/490
-
南京医科大学生殖医学与子代健康国重实验室-董飞宏课题组-学术硕士招生
+1/101
-
物理学 调剂
+1/84
-
何时使用 CODA™ 科里奥利质量流量仪表- 艾里卡特（Alicat）
+2/52
-
山东师范大学有机化学专业胡忠燕老师课题组招收2026届硕士研究生以及调剂生
+1/44
-
华南师范大学（211）博士招生- 电子、自动化、机械、生物学、物理相关专业
+2/42
-
广东工业大学管理科学与工程专业博士招生
+1/34
-
新疆大学招收学硕调剂
+1/24
-
厦门大学航空航天学院智能制造课题组招2026年申请审核制博士生1名
+1/22
-
南京林业大学-国家级青年人才团队 招2026级申请考核制博士（合成化学、植物化学方向）
+1/17
-
博士后招聘（高薪40万+）
+1/7
-
上海科技大学物质学院招收有机化学方向2026级工程博士
+1/7
-
211/双一流石河子大学化学化工学院-电池/电催化研究方向
+1/7
-
26年申博自荐-计算机视觉
+1/6
-
重庆大学药学院闫海龙课题组拟招收2026年申请考核制博士研究生
+1/5
-
深圳大学蔡兴科课题组博士后及研究助理招聘
+1/4
-
【2026年春季补录】厦门大学田中群院士/易骏教授纳米光子学团队招收2026年入学博士生
+1/4
-
天津大学药学院招聘AI生物医疗相关的研究助理、博士、博士后
+1/3
-
284 求调剂
+1/1
-
招收理论凝聚态物理/纳米光学/量子计算方向博士、硕士研究生/博士后
+1/1
3楼2011-03-05 22:40:16
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 ] |
2楼2011-03-03 20:00:42
4楼2011-03-06 08:37:35
