| 查看: 1915 | 回复: 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); |
回复此楼» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有226人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
找到一些相关的精华帖子，希望有用哦~
用matlab，非线性拟合两个参数之间的数量级差距比较大，参数怎么设置权重？
已经有8人回复
- 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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
30