| 查看: 2306 | 回复: 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 |
回复此楼» 猜你喜欢
070200求调剂，一志愿某211，288分
已经有18人回复
-
法国博士后职位
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有275人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有7人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有2人回复
|
我想买一个5.0以上版本的，结果入手了一个5.0的破解版，这个程序好像也运行不了。 哥你能帮我算一次吗？ Variable x, y; ODEFunction y'= k1*(0.701*exp(-3.211*x)+0.299*exp(-0.067*x))-k2*y; Data; //x, y 0 0 0.33 0.061043 1 0.03675 2 0.05932 4 0.095993 6 0.072057 10 0.05085 15 0.04678 24 0.047673 30 0.034973 48 0.030375 72 0 我按教程写成这样了，应该没问题吧？就是算不出结果。 |
9楼2021-03-30 01:33:28
独孤神宇
版主 (知名作家)
- 应助: 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
