应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】各位师哥师姐帮小弟检查一下复杂参数回归的程序吧！！小弟万分感谢！！！
41/1返回列表
查看: 783  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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)
可是运行不出来，小弟是新手，不知道错了哪里，希望各位高手帮我一下，小弟真的是急死了！！！！最好能给出正确的程序，小弟刚入论坛不久，金币不多，但小弟愿拿出一半给大家一表自己的心意。谢谢了！！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
要想拥有一切，就要让自己变得足够优秀！！！
1楼 2011-03-26 10:18:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
xiegangmai(金币+1): 谢谢应助！ 2011-03-27 09:05:07
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);
%%%
%这个函数不对，输入参数是k，你用的是k1k2k3
赞一下(2人) 回复此楼
2楼2011-03-26 23:01:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wodaifei

银虫 (小有名气)
引用回帖:
Originally posted by 信彼南山 at 2011-03-26 23:01:03:
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);
%%%
%这个函数不对，输入参数是k，你用的是k1k2k3

那应该用什么?？？在函数里面输入k1k2k3？?
赞一下 回复此楼
要想拥有一切，就要让自己变得足够优秀！！！
3楼2011-03-28 08:19:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
xiegangmai(金币+1): 谢谢应助！ 2011-03-29 22:56:43
function dy = KineticEqs(t,y,k1,k2,k3)%动力学参数方程
赞一下(2人) 回复此楼
4楼2011-03-29 11:16:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wodaifei 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 oepdf5z972 2026-06-20 3/150 2026-06-21 21:01 by zo4669e5n4
[基金申请] E0414, 我的本子有没有希望？ +7 布布和一二 2026-06-17 9/450 2026-06-21 19:29 by 布布和一二
[考研] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 19:18 by q1p9ixelzp
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 18:58 by q1p9ixelzp
[考研] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 4/200 2026-06-21 17:45 by q1p9ixelzp
[论文投稿] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 17:40 by q1p9ixelzp
[论文投稿] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 17:37 by q1p9ixelzp
[公派出国] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 4/200 2026-06-21 17:25 by ky8v3skvft
[考博] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 5/250 2026-06-21 17:20 by ky8v3skvft
[博后之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 6/300 2026-06-21 17:05 by ky8v3skvft
[公派出国] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 5/250 2026-06-21 17:00 by ky8v3skvft
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +5 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:57 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 vlfdw50o0u 2026-06-20 5/250 2026-06-21 16:45 by ky8v3skvft
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:40 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 vlfdw50o0u 2026-06-20 9/450 2026-06-21 16:37 by ky8v3skvft
[基金申请] 青A35岁以下通知答辩了吗 +3 暨阳一只柴 2026-06-17 3/150 2026-06-21 09:39 by kudofaye
[考博] 【全奖博士/科研助理/博后招生】新加坡南洋理工大学机械与航空航天学院 +3 robohaha 2026-06-15 4/200 2026-06-21 00:36 by lanbing1
[基金申请] 希望面上有个好结果 +7 碧水00 2026-06-16 7/350 2026-06-18 12:18 by wuke100666
[论文投稿] 三区计算机方向期刊推荐 +5 1457340941 2026-06-15 5/250 2026-06-17 13:16 by 会议编辑
[论文投稿] sci论文二审求助 +5 潘倍倍 2026-06-15 5/250 2026-06-16 10:15 by xs74101122
信息提示
关闭
请填处理意见
关闭