| 查看: 829 | 回复: 0 | ||
泡花alkali新虫 (小有名气)
|
[求助]
关于自由基聚合参数拟合的一个matlab代码，求大神看看错哪了？
|
|
我在跑一篇文献上给出的动力学模型，但是最后结果不太对，怀疑是其中的参数有问题，于是我打算用文献中的实验数据去拟合参数，模型中的物质有6种，文献中数据有浓度值的有一种（即单体转化率），参数一共有9个，所以我打算固定其中8个参数不变，去拟合1个参数kp，代码如下，但是跑不出来，求大神改改。 function KineticsMX180 clear all;clc global t kp kd ki kt ka km a b f kp0 tspan t=[0 7 10 20 30]; CAm=[4.45 3.231 2.759 1.936 1.317]; CA0=[4.45 0 0 0 0 0]'; kd=2.37e-4; % 初级自由基生成速率常数,1/min ki=41.22; % 链引发速率常数,L/(mol·min) kt=2.83e11; % 链终止速率常数,L/(mol·min) ka=18.03; % 初级自由基迁入聚合物相速率,1/min km=42.83; % 链转移速率常数,L/(mol·min) a=0.4; b=0.4; f=0.8; kp0=[4.25e8]; %beta0=[kd ki kt ka km a b f kp]; tspan=t; [kp,resnorm,resid,exitflag,output,lambda,jacobian]=lsqnonlin(@OptObjFunc,kp0,[],[],[],tspan,CA0,CAm); ci=nlparci(kp,resid,jacobian); [t4plot CA4plot]=ode23(@KineticsEqs,[tspan(1) tspan(end)],CA0,[],kp); [t CAc]=ode23(@KineticsEqs,tspan,CA0,[],kp); function f=OptObjFunc(kp,tspan,CA0,CAm) global t kp kd ki kt ka km a b f kp0 tspan [t CAc]=ode23(@KineticsEqs,tspan,CA0,[],kp); CAc f=(CAm'-CAc); function dCAdt=KineticsEqs(t,CA,kp) global t kp kd ki kt ka km a b f kp0 tspan kd=2.37e-4; % 初级自由基生成速率常数,1/min ki=41.22; % 链引发速率常数,L/(mol·min) kt=2.83e11; % 链终止速率常数,L/(mol·min) ka=18.03; % 初级自由基迁入聚合物相速率,1/min km=42.83; % 链转移速率常数,L/(mol·min) a=0.4; b=0.4; f=0.8; CA c1=CA(1,: ); c2=CA(2,: ); c3=CA(3,: ); % C=k2*c1+k3*c2+k4*c3+k5*c4 Co=c6+c7+c8+c9 c4=CA(4,: ); c5=CA(5,: ); c6=CA(6,: ); dc1dt=-kd*c1^a*c6^b; dc2dt=f*kd*c1^a*c6^b-ka*c2; dc3dt=ka*c2-ki*c3*c5-kt*c3*c4; dc4dt=ki*c3*c5-kt*c3*c4; dc5dt=kp*c5*c4; dc6dt=-3*kd*c1^a*c6^b; dCAdt=[dc1dt;dc2dt;dc3dt;dc4dt;dc5dt;dc6dt]; 其中c1~c6的：后加了空格是因为它会自动变成微笑表情。。 |
回复此楼» 猜你喜欢
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
物理学I论文润色/翻译怎么收费?
已经有103人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有0人回复
找到一些相关的精华帖子，希望有用哦~
使用MATLAB拟合动力学参数问题
已经有18人回复
- 如何用Matlab曲线拟合一个多输入参数，而且其中一个输入参数与另一个参数有关呢？
已经有6人回复
- 拜求！matlab拟合动力学参数
已经有14人回复
- ]紧急求助复介电常数的模型拟合-Matlab or 1stOpt
已经有9人回复
- 关于matlab两个自变量、一个因变量曲线拟合问题
已经有3人回复
- 求助用matlab 多元非线性拟合
已经有4人回复
- matlab 拟合求助
已经有8人回复
- 求助matlab拟合三维曲面
已经有10人回复
- 求助，matlab自定义微分方程拟合实验数据来求方程中的参数
已经有4人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- Matlab 含有两个参数的公式拟合
已经有4人回复
- matlab数据拟合问题求助
已经有7人回复
- 求助 用matlab拟合三元线性回归方程及分析
已经有9人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- matlab拟合求值
已经有17人回复
- matlab 求指点 动力学方程拟合过程中导数的获取
已经有10人回复
- matlab非线性参数拟合问题
已经有7人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 【求助】使用MATLAB拟合出按照公式的系数
已经有19人回复
- 【求助】matlab曲面拟合表达式
已经有10人回复
- 【求助】如何使得matlab拟合的曲线强制经过一个点
已经有8人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
精华III:
20