| 查看: 1859 | 回复: 0 | ||
[求助]
用Matlab求反应动力学参数
|
|
中间先拟合为一段一级和一段二级，并求出各点的导数，最后拟合出需要的参数。拟合过程中出现导数为负数，求各位大神指导！ 代码如下： clc; clear all; format short g %%%分段求解 一级 expE1=[0 2054.794521 300 903.6835616 600 449.3047945]; t1=expE1(:,1); C0=expE1(1,2); C1=expE1(:,2); LCO=log(C0); LC=log(C1); d=LCO-LC; K0=[0.02]; K=lsqcurvefit(@FunS,K0,t1,d); mdcdt1=K.*C1; %%%%分段求解 二级 expE2=[900 200.3150685 1200 107.0773973 1500 78.67123288 2100 57.00136986 2700 47.67328767 3300 38.26164384 3900 28.4 4500 21.78219178]; C02=expE2(1,2); C2=expE2(:,2); t2=expE2(:,1); RC2=1./C2; d2=RC2-1/C02; K20=[0.02]; K2=lsqcurvefit(@FunS2,K20,t2,d2); mdcdt2=K2.*(C2.^2); %%%%合并求参数 mdcdt=[mdcdt1;mdcdt2]; C=[C1;C2]; %%%% T=250C下，mnts初值P0 P0=[89 0.0328 110 2284.62585]; P=lsqcurvefit(@FunP,P0,C,mdcdt) mdcdt3=(P(1)+P(2).*C)./(P(3)+1/P(4)./C./C) figure plot(C,mdcdt,'r-'); hold on; plot(C,mdcdt3,'b*') function f=FunS(K,t1) f=K(1).*t1; function f=FunS2(K2,t2) f=K2(1).*t2; function f=FunP(P,C) f=(P(1)+P(2).*C)./(P(3)+1/P(4)./C./C); |
回复此楼» 猜你喜欢
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
物理学I论文润色/翻译怎么收费?
已经有126人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有0人回复
找到一些相关的精华帖子，希望有用哦~
使用MATLAB拟合动力学参数问题
已经有18人回复
- 用matlab，非线性拟合两个参数之间的数量级差距比较大，参数怎么设置权重？
已经有8人回复
- Matlab同时拟合活化能，反应级数，速率常数
已经有26人回复
- matlab求模型参数有错误，求指导
已经有5人回复
- 使用matlab最优化方法拟合获得多个动力学参数中的问题
已经有4人回复
- matlab动力学模型编程求助
已经有6人回复
- matlab拟合参数显示初始点是一个局部最小值
已经有18人回复
- 拜求！matlab拟合动力学参数
已经有14人回复
- 求助用matlab拟合动力学方程
已经有9人回复
- 求助用matlab拟合动力学方程
已经有13人回复
- 酶动力学参数计算，Kcat
已经有5人回复
- 关于matlab的参数估计
已经有15人回复
- 谁有成熟的求动力学参数的MATLAB的程序。
已经有11人回复
- 动力学参数拟合
已经有26人回复
- "ploylog()"MATLAB求教！！！
已经有6人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- 请教matlab反应动力学参数估计遇到的问题，谢谢
已经有15人回复
- matlab拟合拟均相动力学参数
已经有6人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 【求助】使用Matlab拟合反应动力学方程问题
已经有7人回复
- 【求助】使用Matlab预估动力学方程问题
已经有13人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
精华III:
20