| 查看: 2076 | 回复: 3 | |||
[交流]
【求助】催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？
|
|
催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？ 我按照个基元反应，依据控制速率步骤，得到相应速率表达式，从而求的各个基元反应的速率常数，但是在设定初值后，结果再次迭代最后算的结果，但是把最后的结果当初值，同时改掉其中一个变量，有些变量对最后结果影响很大，有些却几乎不影响其他结果，求高手帮忙解决！  同时求其他更好的解动力学方程组的方法！ |
回复此楼» 猜你喜欢
国家级人才团体课题组招收2026届化学类学硕，专硕。
已经有1人回复
-
河南理工大学绿色催化课题组招收材料，化工，化学等方向研究生
已经有0人回复
-
专业论文润色/翻译怎么收费?
已经有62人回复
-
国家级人才课题组招研究生 （催化方向）调剂
已经有13人回复
-
招收硕士研究生
已经有24人回复
-
[天津]中国民航大学 招收材料类研究生
已经有0人回复
-
材料与化工专业有大量调剂指标
已经有8人回复
-
国家大人才课题组招收化学、材料化工类2026年调剂生
已经有0人回复
-
085600调剂
已经有4人回复
-
浙师大化与材学院2026年有大量硕士研究生调剂名额 （线上 调剂）
已经有6人回复
-
吉林工程技术师范学院 材料与化工 （0856） 招收调剂
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
» 抢金币啦！回帖就可以得到:
-
广西师范大学-石伟课题组2026年招收药学学硕 、生物与医药专硕调剂生
+1/482
-
江西科技师范大学 材料与能源招收调剂研究生
+1/92
-
安徽省大别山中医药研究院 接收中药调剂生（中药、药学等10开头代码均可）
+1/77
-
重庆工商大学 教育部研究中心招收化工/环境/材料/石油/机械方向2026年调剂研究生
+1/53
-
2026年河北大学无机化学专业招收调剂生
+1/40
-
安徽工业大学有调剂名额
+1/38
-
国家级人才课题组招研究生 （催化方向）调剂
+1/37
-
南京信息工程大学 化材院接受调剂研究生
+1/34
-
生物与医药专业硕士调剂（河北大学 化学与材料科学学院）
+2/28
-
招收地理类、环境类、地化类、生物类、土壤类、农业资源与环境等硕士调剂生
+1/11
-
武科大冶能学院招收材料、冶金、能源环境、机械等相关专业调剂生
+1/9
-
广州医科大学生物医学工程生物技术与工程（086001）和制药工程（086002）接收调剂
+1/9
-
齐齐哈尔大学李莉课题组诚招2026级考研调剂生（学硕和专硕）
+1/8
-
徐工-环境工程学院-招收调剂硕士
+1/6
-
福建理工大学招收专业代码08开头学术型硕士研究生
+1/5
-
渤海大学材料与环境工程学院发光及透明陶瓷材料、地聚物吸波混凝土方向接收研究生调剂
+1/5
-
南昌大学工程建设学院2026年招聘2名科研助理
+1/4
-
调剂招生化学化工（0817、0856）
+1/3
-
齐齐哈尔大学化学与化学工程学院招收专硕调剂
+1/1
-
北京高校副校长团队招收机械类，环境类学硕和专硕
+1/1
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
3楼2011-03-05 22:40:16
4楼2011-03-06 08:37:35
