| 查看: 805 | 回复: 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的：后加了空格是因为它会自动变成微笑表情。。 |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
物理学I论文润色/翻译怎么收费?
已经有125人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有2人回复
找到一些相关的精华帖子，希望有用哦~
使用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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
