| 查看: 2219 | 回复: 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 |
回复此楼» 猜你喜欢
寻合作：应力腐蚀多尺度模拟
已经有6人回复
-
考研交流
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有288人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
帮我的英语口语老师找学生
已经有3人回复
-
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有0人回复
-
重庆交通大学光子学微结构与器件课题组2026年招收调剂硕士研究生信息
已经有0人回复
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31016.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 856.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
|
这是免费的试用版，你这被人忽悠了啊，最新版已经到9.0了，而且只能去官网购买。 均方差(RMSE): 0.0171691634594409 残差平方和(SSR): 0.00324258191286701 相关系数(R): 0.802484090026993 相关系数之平方(R^2): 0.64398071474645 修正R平方(Adj. R^2): 0.554975893433063 确定系数(DC): 0.468028249120648 F统计(F-Statistic): 23.5213130305514 参数 最佳估算 -------------------- ------------- k1 0.13256115104059 k2 0.358504426638858 ====== 结果输出 ====== 文件: 数据文件-1 No x 目标 y 计算 y 1 0.33 0.061043 0.0298491702527491 2 1 0.03675 0.053602315170092 3 2 0.05932 0.0683875960100214 4 4 0.095993 0.079360928514414 5 6 0.072057 0.0789300799888936 6 10 0.05085 0.0667093831255373 7 15 0.04678 0.0492934543907035 8 24 0.047673 0.02721397228202 9 30 0.034973 0.0182160615899268 10 48 0.030375 0.00545427669955797 11 72 0 0.00109238646097193 |
10楼2021-03-30 08:47:04
独孤神宇
版主 (知名作家)
- 应助: 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
