| 查看: 3552 | 回复: 10 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
matlab 求指点 动力学方程拟合过程中导数的获取
|
|||
|
本人matlab不通，勉强做了一个程序能够达到预期目的，可是老师想要中间运算的导数值，也就是dxdt的值，不知道如何填补表达式实现，望各位高手出手相助，不胜感激。 主程序： format short global Umax Ks Kp Y1 Y2 Uo Ko D Ux Kx Ku Ki; Ki=0.0001;Ku=0.0006; k0=[0.2 0.5 0.1 10 10 0.01 0.01 0.1 0.01 0.01]; %参数的初始值 x0=[0.0514 0.250 0]; %菌体浓度、底物浓度和产物浓度的初始值 t1=[0 6 12 18 24 30 36 42 48 54 60 66]';%发酵时间的实验数据 tspan=[0 6 12 18 24 30 36 42 48 54 60 66 ]';%时间间隔 %yexp实验数据[菌体浓度x1、底物浓度x2和产物浓度x3] yexp=[[0.0514 0.0711 0.0997 0.1887 0.3478 0.6634 0.8711 1.1023 1.2336 1.4239 1.5849 1.6739];[0.25 0.2447 0.2366 0.2258 0.2196 0.203 0.1872 0.1698 0.1453 0.1124 0.0855 0.0608];[0 0.0012 0.0021 0.0039 0.0088 0.0137 0.0183 0.0268 0.0361 0.0438 0.0557 0.0643]]'; lb=[0 0 0 0 0 0.0216 0.001 0 0.001 0.001];ub=[1 10 10 50 50 2 1 1 0.1 0.1];%参数的下、上限 [k,resnorm,residual,exitflag,output,lambda,jacobian]=lsqnonlin(@ObjFunc4LNL103,k0,lb,ub,[],x0,yexp);%非线性拟合 y1=[yexp(:,1)];y2=[yexp(:,2)];y3=[yexp(:,3)]'; [t4plot,x4plot]=ode45(@kineticseqs103,[0 100],x0,[],k); plot(t1,y1,'b*',t4plot,x4plot(:,1),'k-'),xlabel('T(h)'),ylabel('Cx/（g/L）'); figure plot(t1,y2,'g*',t4plot,x4plot(:,2),'k-'),xlabel('T(h)'),ylabel('Cs/(mmol/L)'); figure plot(t1,y3,'r*',t4plot,x4plot(:,3),'k-'),xlabel('T(h)'),ylabel('Cp/(mmol/L)'); disp(k) 子程序： 1） function f=ObjFunc4LNL103(k,x0,yexp) tspan=[0 6 12 18 24 30 36 42 48 54 60 66]; [t1,x]=ode45(@kineticseqs103,tspan,x0,[],k); y(:,1)=x(:,1);y(:,2)=x(:,2);y(:,3)=x(:,3); f1=y(: ,1)-yexp(: ,1);f2=y(: ,2)-yexp(: ,2);f3=y(: ,3)-yexp(: ,3); f=[f1*1 f2*5 f3*10]; 2） function dxdt=kineticseqs103(t,x,k) %模型方程 global Umax Ks Kp Y1 Y2 Uo Ko D Ux Kx Ku Ki Umax=k(1); Ks=k(2); Kp=k(3); Y1=k(4); Y2=k(5); Ux=k(6); Kx=k(7); D=k(8); Uo=k(9); Ko=k(10); dxdt1=Umax*x(2)*x(1)*(1-Kp*x(3))/(Ks+x(2)); if dxdt1<0 dxdt1=0; end dxdt2=-1/Y1*dxdt1-Ux*x(2)*x(1)/(Kx+x(2))*(Ki/(Ki+x(3))); dxdt3=D*(-dxdt2-1/Y2*dxdt1-Uo*x(2)*x(1)/(Ko+x(2))*(Ku/(Ku+x(3)))); if dxdt3<0 dxdt3=0; end dxdt=[dxdt1;dxdt2;dxdt3]; 只要能显示对应时间t1,这些点的dxdt值就可以了。谢谢 |
回复此楼» 收录本帖的淘帖专辑推荐
 matlab |
» 本帖已获得的红花（最新10朵）
shidoudou» 猜你喜欢
帮我的英语口语老师找学生
已经有4人回复
-
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
物理学I论文润色/翻译怎么收费?
已经有209人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有4人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab拟合拟均相动力学参数
已经有6人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 吸附的动力学方程拟合
已经有18人回复
- matlab可否求离散三维函数（即一组（x，y，z）数据）的偏导数（如：Dz/Dy）？
已经有5人回复
- 【求助】拟合动力学方程求助
已经有13人回复
- origin，sigmaplot等软件自定义曲线拟合（主要是一些动力学方程的模拟）
已经有2人回复
- 【求助】光催化动力学方程拟合
已经有7人回复
- 【求助】使用Matlab拟合反应动力学方程问题
已经有7人回复
- 【求助】使用Matlab预估动力学方程问题
已经有13人回复
715211229
木虫 (正式写手)
- 应助: 52 (初中生)
- 金币: 3652.8
- 散金: 20
- 红花: 6
- 帖子: 339
- 在线: 249.2小时
- 虫号: 1357002
- 注册: 2011-07-27
- 性别: GG
- 专业: 传动机械学
8楼2012-02-11 00:49:28
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2012-02-09 17:16:45
3楼2012-02-09 20:30:07
hhucy
实习版主 (著名写手)
DOE锅炉工
- 应助: 53 (初中生)
- 贵宾: 0.264
- 金币: 7775.6
- 散金: 1861
- 红花: 13
- 沙发: 7
- 帖子: 2588
- 在线: 1275.5小时
- 虫号: 1079650
- 注册: 2010-08-22
- 性别: GG
- 管辖: 能源
4楼2012-02-09 20:40:21
5