| 查看: 290 | 回复: 0 | |||
wodaifei银虫 (小有名气)
|
[交流]
【求助】各位师哥师姐帮小弟检查一下复杂参数回归的程序吧！！小弟万分感谢！！！
|
|
各位师哥师姐，小弟求助： dy1/dt=k1y1y2-k2y3 dy2/dt=k1y1y2-k2y3 dy3/dt=k2y3 y1(0)=10,y2(0)=8,y3(0)=0 小弟想用四阶-龙阁库塔积分和非线性最小二乘估计以上的动力学参数k1,k2,k3,数据如下：t y1 y2 y3 0 10 8 0 20 8 7 5 40 6 6 7 60 5 5 8 120 3 3 10 180 2 2 12 300 1 1 14 以下是小弟的程序： k0 = [1,2,3]; %设参数初值 y0 = [10,8,0]; tspan = [0,20,40,60,120,180,300]; %反应时间 y1 = [10,8,6,5,3,2,1]';%测得的第一组数据 y2 = [8,7,6,5,3,2,1]';%测得的第二组数据 y3 = [0,5,7,8,10,12,14]';%测得的第三组数据 yexp = [y1;y2;y3]; function dy = KineticEqs(t,y,k)%动力学参数方程 dy = zeros(3,1); dy(1)=k1*y(1)*y(2)-k2*y(3); dy(2)=k1*y(1)*y(2)-k2*y(3); dy(3)=k2*y(3); function f =OptObjFunc(k,y0,yexp) %目标函数 [t Y] = ode45(@KineticEqs,tspan,y,[],k); f1 = y(1) - y1; f2 = y(2) - y2; f3 = y(3) - y3; f = [f1; f2; f3]; [k,resnorm,resid,exitflag,output,lambda,jacobian] = ...%最小二乘法 lsqnonlin(@OptObjFunc,k0,[0 0],[],[],tspan,y0) ci = nlparci(k,resid,jacobian)%置信区间 fprintf('Estimated Parameters:\n') fprintf('\tk1 = %.4f ± %.4f\n',k1) fprintf('\tk2 = %.4f ± %.4f\n',k2) 可是运行不出来，小弟是新手，不知道错了哪里，希望各位高手帮我一下，小弟真的是急死了！！！！最好能给出正确的程序，小弟刚入论坛不久，金币不多，但小弟愿拿出一半给大家一表自己的心意。谢谢了！！ [ Last edited by wodaifei on 2011-3-26 at 10:07 ] |
回复此楼» 猜你喜欢
299求调剂
已经有8人回复
-
一志愿北京理工大学本科211材料工程294求调剂
已经有6人回复
-
300求调剂，材料科学英一数二
已经有8人回复
-
招收生物学/细胞生物学调剂
已经有5人回复
-
070305高分子化学与物理 304分求调剂
已经有7人回复
-
289求调剂
已经有13人回复
-
一志愿哈尔滨工业大学材料与化工方向336分
已经有9人回复
-
081200-11408-276学硕求调剂
已经有6人回复
-
调剂求院校招收
已经有5人回复
-
调剂310
已经有8人回复
找到一些相关的精华帖子，希望有用哦~
SFDA认证中心2010版药品GMP认证检查结果判定程序
已经有14人回复
- 求助格氏反应，小弟感谢！！！
已经有9人回复
- 小弟初学生物信息学，完全摸不着头脑，求启蒙
已经有7人回复
- 关于XPS的几个问题，小弟先再次谢过了
已经有8人回复
- 求助申请文部省奖学金的师兄师姐们帮小弟解解惑呀
已经有11人回复
- 请质谱高人帮忙解析质谱，万分感谢！
已经有14人回复
- 请教各位师兄师姐们了~~不胜感激
已经有5人回复
- 求助，各位出国读博的师兄师姐，请给我点建议，小弟在此拜谢！
已经有8人回复
- 【求助】各位师哥师姐帮小弟检查一下复杂参数回归的程序吧！！小弟万分感谢！！！
已经有3人回复
- 【求助】师兄师姐们啊 wingslib快搞死我了 帮帮小弟吧
已经有4人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
20