| 查看: 2061 | 回复: 3 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？
|
|||
|
催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？ 我按照个基元反应，依据控制速率步骤，得到相应速率表达式，从而求的各个基元反应的速率常数，但是在设定初值后，结果再次迭代最后算的结果，但是把最后的结果当初值，同时改掉其中一个变量，有些变量对最后结果影响很大，有些却几乎不影响其他结果，求高手帮忙解决！  同时求其他更好的解动力学方程组的方法！ |
回复此楼» 猜你喜欢
0817一志愿苏州大学280复试调剂
已经有9人回复
-
0817找调剂
已经有6人回复
-
专业论文润色/翻译怎么收费?
已经有112人回复
-
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/168
-
深圳技术大学集成电路与光电芯片学院郝俊杰课题组诚聘博士后、研究助理、访问学生
+2/132
-
招收苏州国家实验室和苏州大学联合培养博士（2026年9月入学）
+1/82
-
福建师范大学招收2026年化学、材料硕士3-4名
+1/82
-
【武汉高校】【新能源材料与器件课题组】 招收优秀硕士调剂生！！！
+5/35
-
教研论文SCI期刊投稿选刊
+1/34
-
【沈阳农业大学机械专硕调剂】
+2/30
-
2026年工科硕士调剂-上海大学全国重点实验室团队-材料数据挖掘方向-研究生3-5人
+1/30
-
中山大学院士团队微纳器件、脑机接口方向博士招生、博后招聘
+1/22
-
课题组招收环境及相关专业调剂硕士研究生（欢迎优秀学生加入）
+1/7
-
教育部长江学者和创新团队发展计划”入选团队招收 材料，化学与化工博士研究生
+1/7
-
福州大学国家级人才汤育欣教授课题组2026年招收”申请考核制“储能电池方向博士研究生
+1/6
-
海南大学徐月山老师招生第二批博士名额2~3个，2026年9月份入学（高端设备开发方向）
+1/6
-
发paper
+1/5
-
药学专硕一志愿苏大本科苏大寻求B区高校调剂名额
+1/2
-
西安交通大学生命学院智能化诊疗装备研究中心2026年博士生招生
+1/2
-
新加坡国立大学刘小钢课题组招聘博士后
+1/2
-
层流压差式流量测控产品在环保与分析仪器领域的应用与优势
+1/2
-
五邑大学环境与化学化工学院欢迎化学、化工、材料背景考生报考
+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
