| 查看: 707 | 回复: 3 | |||
wodaifei银虫 (小有名气)
|
[交流]
【求助】各位师哥师姐帮小弟检查一下复杂参数回归的程序吧！！小弟万分感谢！！！ 已有2人参与
|
|
各位师哥师姐，小弟求助： 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) 可是运行不出来，小弟是新手，不知道错了哪里，希望各位高手帮我一下，小弟真的是急死了！！！！最好能给出正确的程序，小弟刚入论坛不久，金币不多，但小弟愿拿出一半给大家一表自己的心意。谢谢了！！ |
回复此楼» 猜你喜欢
0855求调剂材料
已经有12人回复
-
297分083200求助
已经有7人回复
-
070300化学学硕311分求调剂
已经有10人回复
-
专硕0854初试考材科基，求调剂
已经有8人回复
-
0857大类环境工程B区求调剂
已经有3人回复
-
0817化学工程与技术求调剂，一志愿中海洋319
已经有14人回复
-
0703化学调剂325分
已经有11人回复
-
286求调剂
已经有8人回复
-
求调剂
已经有5人回复
-
生物学308分求调剂（一志愿华东师大）
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- SFDA认证中心2010版药品GMP认证检查结果判定程序
已经有14人回复
- 求助格氏反应，小弟感谢！！！
已经有9人回复
- 小弟初学生物信息学，完全摸不着头脑，求启蒙
已经有7人回复
- 关于XPS的几个问题，小弟先再次谢过了
已经有8人回复
- 求助申请文部省奖学金的师兄师姐们帮小弟解解惑呀
已经有11人回复
- 请质谱高人帮忙解析质谱，万分感谢！
已经有14人回复
- 请教各位师兄师姐们了~~不胜感激
已经有5人回复
- 求助，各位出国读博的师兄师姐，请给我点建议，小弟在此拜谢！
已经有8人回复
- 【求助】师兄师姐们啊 wingslib快搞死我了 帮帮小弟吧
已经有4人回复
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
2楼2011-03-26 23:01:03
wodaifei
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3339.1
- 散金: 572
- 红花: 8
- 帖子: 248
- 在线: 351.1小时
- 虫号: 683543
- 注册: 2008-12-28
- 性别: GG
- 专业: 化工系统工程
3楼2011-03-28 08:19:53
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
4楼2011-03-29 11:16:40
