| 查看: 4151 | 回复: 26 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[交流]
动力学参数拟合
|
||||
|
CO变换反应，CO+H2O=CO2+H2，用固定床管式积分反应器，要拟合动力学方程中的参数！动力学方程如下图。转化后的东西学方程为dy/dw=-22.4./v*k(1)*exp(-k(2)/R.*T).*yco^k(3).*(yh2oin-(ycoin-yco))^k(4).*(yco2in+(ycoin-yco))^k(5).*(yh2in+(ycoin-yco))^k(6).*(1-β)，其中β是逆反应程度，要拟合指前因子、反应活化能和四个反应级数（k1~k6）。整理后的数据是已知反应器进口一氧化碳、二氧化碳、氢气、水蒸气的湿基组成和反应器出口一氧化碳湿基组成，主要是一个微分方程dy/dw=f(y),其中f(y)中含有要拟合的6个参数，涉及到微分方程的初值问题和最小二乘法问题。自己编了一个m文件，可是运行不了，纠结好长时间也没有发现问题，想大家帮我看一下！代码如下 function kinetics clear all clc global t T R kp ycoin yh2oin yco2in yh2in v t=[448.69 449.03 450.18 450.71 450.47 430.01 430.31 430.03 430.17 427.29 410.13 410.08 410.37 410.00 409.96 390.00 390.00 389.99 390.00 390.00];%反应温度 T=t+273.15;%绝对温度K R=8.314;%气体常数 kp=exp(4575./T-4.33);%平衡常数 ycoin=[0.09865 0.13296 0.18299 0.23883 0.13809 0.13918 0.15477 0.14327 0.16374 0.13077 0.10624 0.11253 0.19904 0.20900 0.19895 0.12235 0.17431 0.14712 0.18064 0.16864];%进口CO含量 yh2oin=[0.52349 0.45073 0.39759 0.33124 0.45811 0.30844 0.39573 0.51067 0.54215 0.32203 0.48650 0.54179 0.32176 0.42416 0.32176 0.41059 0.32107 0.50444 0.48270 0.31474];%进口H2O含量 yco2in=[0.09440 0.08103 0.07302 0.05747 0.09028 0.11332 0.11931 0.04328 0.05115 0.10654 0.06181 0.03400 0.13821 0.09236 0.13688 0.04807 0.07740 0.07265 0.09696 0.08260];%进口CO2含量 yh2in=[0.09363 0.14281 0.17864 0.22005 0.13013 0.20040 0.20556 0.10088 0.10997 0.20593 0.17479 0.13200 0.16644 0.09602 0.17151 0.14214 0.13702 0.18113 0.16127 0.13905];%进口H2含量 v=[21.3952 18.2224 15.0234 13.5257 19.5599 14.5256 18.8683 19.0081 20.6600 16.2300 18.4185 20.6734 13.9253 15.8444 13.9233 16.3679 13.9538 20.7244 17.0155 14.2337];%总原料流量 ycoout=[0.07304 0.08526 0.12500 0.16207 0.08445 0.10695 0.13233 0.10264 0.11417 0.09878 0.09192 0.08657 0.17122 0.17377 0.17204 0.10501 0.14913 0.13370 0.15282 0.14068];%出口CO含量 k0=[3.5e5 8.4e4 0.9 0.6 -0.5 -0.04]; % k(指前因子、活化能、反应级数)初始向量 lb=[1e5 3e4 0 0 -1 -0.1]; ub=[1e6 1.5e5 1 1 0.5 0.5];% k上下限 % 使用函数lsqnonlin()进行参数估计 [k,resnorm,residual,exitflag]=lsqnonlin(@ObjFunc,k0,lb,ub,[],ycoin,ycoout); % ------------------------------------------------------------------ function f=ObjFunc(k,ycoin,ycoout) % 目标函数 wspan=[0 0.6091]; % 即催化剂装填量/g ycow=zeros(1,length(ycoout)); [w yco] =ode45(@Euqations,wspan,ycoin,[],k);%4,5阶龙格库塔法求解 ycow=yco(length(yco));%出口CO含量 f=ycoout-ycow; % --------此方程定义微分方程等号右端项 function dydw=Euqations(w,yco,k) % w催化剂质量 global T R kp ycoin yh2oin yco2in yh2in v dydw=-22.4./v*k(1)*exp(-k(2)/R.*T).*yco^k(3).*(yh2oin-(ycoin-yco))^k(4).*(yco2in+(ycoin-yco))^k(5).*(yh2in+(ycoin-yco))^k(6).*(1-(yco2in-(ycoin-yco).*yh2in-(ycoin-yco))./(kp.*yco.*yh2oin-(ycoin-yco)));% 反应速率  [ 来自科研家族 化工工程工艺 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 matlab |
» 本帖已获得的红花（最新10朵）
dingd» 猜你喜欢
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有147人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
-
计算与实验科学工程国际会议(8月7-11日,香港) 二维材料分会S5与声学换能器分会S8
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教matlab反应动力学参数估计遇到的问题，谢谢
已经有15人回复
- matlab 求指点 动力学方程拟合过程中导数的获取
已经有10人回复
- matlab拟合拟均相动力学参数
已经有6人回复
- 最小二乘拟合 初值 运算时间
已经有12人回复
- 吸附的动力学方程拟合
已经有18人回复
- 固定床中气-固相脱水反应本征动力学
已经有5人回复
- 【讨论】关于腐蚀速度与动力学参数
已经有15人回复
- origin，sigmaplot等软件自定义曲线拟合（主要是一些动力学方程的模拟）
已经有2人回复
- 【求助】如何从极化曲线求动力学参数
已经有7人回复
» 抢金币啦！回帖就可以得到:
-
【招聘】温州医科大学/附属医院/瓯江实验室李校堃院士团队诚聘博士后和科研人员
+1/187
-
耐热氧化耐热黄变效果好的聚醚多元醇制备的聚氨酯都有哪些潜在应用领域？
+2/142
-
Postdoc position in Hangzhou Normal University （陈健翔组博后招聘）
+1/83
-
资源与环境专硕
+1/72
-
北京-89175-事业单位-诚征女友
+1/63
-
武汉工程大学化环学院程老师课题组招收化学类学硕、化学工程与环境工程专硕调剂
+2/50
-
青岛大学固态电池研究团队招收材料、化学、物理等方向博士研究生
+2/44
-
天津城建大学招收物理、化学或者生物调剂生一名，前往天津大学联合培养
+1/38
-
招聘︱国家纳米科学中心丁宝全杰青团队诚聘博士后（特别研究助理）
+2/34
-
国家双一流高校-国家级青年人才课题组博士招生
+2/30
-
招聘科研助理——储能方向（储能、流体传热、CFD模拟）
+1/18
-
双一流大学湘潭大学“化工过程模拟与强化”国家地方联合工程研究中心招收博士生
+1/17
-
博士招生 | 青岛科技大学（高分子、化学、材料、力学、航天、机电工程、计算机仿真）
+1/13
-
青岛科技大学0860 调剂招生
+1/11
-
西华大学材料学院表面科学与工程技术科研团队2026年招收研究生
+1/11
-
同济大学环境学院 肖倩特聘研究员课题组 招聘硕士/博士（长期有效）
+1/6
-
NTU招聘生物微机电系统/微流控/3D打印/可穿戴器件方向全奖博士/博士后
+1/6
-
山东济南省属重点高校生态环境类硕士调剂信息
+1/2
-
中国科学技术大学纳米科学技术学院张振教授/江雷院士诚招2026级博士
+1/1
-
南开大学环境科学与工程学院王鑫教授课题组招收申请考核制博士生
+1/1
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
csgt0: 金币+2, 谢谢 2012-11-17 09:08:43
小木虫: 金币+0.5, 给个红包，谢谢回帖
csgt0: 金币+2, 谢谢 2012-11-17 09:08:43
|
Optimization terminated: relative function value changing by less than OPTIONS.TolFun. k = 1.0e+005 * Columns 1 through 4 1.000003573678879 0.762988537687594 0.000005843208004 0.000005981921222 Columns 5 through 6 -0.000000744621941 -0.000000996174141 resnorm = 5.657234679833130e-004 residual = -0.006805767672396 0.002960729069534 -0.006244673627901 -0.001713302178392 0.010027828622409 0.004669076015188 -0.006971457219631 -0.002080652235825 0.005813395104906 0.007827253485404 -0.008382737359081 0.000669267176915 -0.001261508186825 -0.002006057821641 -0.001995402177405 -0.001123333955809 0.003628239382343 -0.004182989575642 0.004948094412108 0.007850845335425 exitflag = 3 |
25楼2012-11-15 16:59:53
8楼2012-06-01 10:18:18
9楼2012-06-01 12:01:47
11楼2012-06-01 14:47:56
