| 查看: 2326 | 回复: 10 | ||
| 【悬赏金币】回答本帖问题，作者Evoly_z将赠送您 88 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求助代谢动力学系数模拟代码 1stOpt或者MATLAB 已有3人参与
|
||
 求代谢动力学系数模拟一级代谢物的动力学方程：dC2/dt=k1*C1-k2*C2 初始条件：t=0，C2=0 且C1=exp（-A*t）。A=0.2779 我尝试了用1stOpt(破解版)和MATLAB ODE方法，都没成功，想请教一下大神。 另外t不是严格的等差数列，取值如：t=0，1，2，4，6，10，15，24 1stOpt代码： Title Kinetic_ave Parameters k1[0,100], k2[0,100]; Variable t, C; StartProgram var i:integer; begin for i:=0 to DataLength -1 do begin if i ==0 C=0; else C:=C[i-1]+k1*(t-t[i-1])*exp(-0.2779 *t) - k2*C*(t-t[i-1]);end; EndProgram; Data; //t C 0 xxx 1 xxx 2 xxx 4 xxx 6 xxx 10 xxx 15 xxx 24 xxx Matlab代码： function ODE_ave clear all;clc format long aveall; t=T_h(  ;yexp=OLEave( ;k0=[1 1]; y0=0; lb=[0 0]; ub=[+inf +inf]; yy=[y0 yexp']; tspan=0:1:24; [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,y0,yexp); fprintf('\n\n使用函数lsqnonlin()估计得到的参数值为:\n') fprintf('\t待拟合参数 k1 = %.6f\n',k(1)) fprintf('\t待拟合参数 k2 = %.6f\n',k(2)) fprintf(' \t残差平方和= %.6f\n\n',resnorm) ts=0:1:24; [ts ys]=ode45(@KineticsEqs,ts,y0,[],k); [ttt XXsim] = ode45(@KineticsEqs,tspan,y0,[],k); y=XXsim(2:end); xexp=yexp; R2=1-sum((xexp-y).^2)./sum((xexp-mean(y)).^2); fprintf('\n\t决定系数R-Square = %.6f',R2); figure(1) plot(ts,ys,'b',tspan,yy,'or'),legend('计算值','实验值','Location','best'); yr=y-yexp; figure(2) plot(tspan(2:end),yr,'r*',[-1 15],[0 0]),axis([-1 15 -0.5 0.5]); figure(3) plot(yexp,y,'ro',[21 29],[21 29],'b-'); (作图这块儿是copy的，没有做修改) %--------------------------------------------------------- function f = ObjFunc(k,tspan,y0,yexp) [t Xsim] = ode45(@KineticsEqs,tspan,y0,[],k) ; ysim = Xsim(2:end); size(ysim); size(yexp); f=ysim(1,1)+ysim(2,1)+ysim(4,1)+ysim(6,1)+ysim(10,1)+ysim(15,1)+ysim(24,1) - sum(yexp(:,1)); %---------------------------------------------------------- function dydt = KineticsEqs(t,y,k) beta(1)=k(1); beta(2)=k(2); dydt = beta(1)*exp(-0.2779*t)-beta(2)*y; 求求啦，被这个问题卡了两个多月了，不知道怎么解出k1 k2 |
回复此楼» 猜你喜欢
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
物理学I论文润色/翻译怎么收费?
已经有170人回复
-
桂林理工大学物理学专业招收调剂
已经有18人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有11人回复
5楼2021-03-24 20:15:22
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31016.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 856.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
2楼2021-03-24 07:20:06
3楼2021-03-24 11:13:13
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2021-03-24 16:12:00
20