24小时热门版块排行榜    

查看: 1105  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gxy2011

新虫 (小有名气)

[求助] 1stOpt非线性拟合问题 已有1人参与

各位大神,在拟合微分方程时出错,方程为dy/dx= 3.01*z/(8+z)-a*y/(b+y),运行错误,是哪出问题了呢?

Parameter a,b;
variable x,z,y
ODEFunction y'=3.01*z/(8+z)-a*y/(b+y);
Data;
//x      z       y
0        6.6136        92.9004
5        19.5602        77.9298
10        29.3171        67.1342
15        37.5293        55.1337
20        49.2231        35.2848
25        52.6523        27.7856
30        59.0752        17.7677
35        62.693        7.8028
40        62.4031        1.3144
45        57.7333        0.9989
50        50.1079        0.7628
55        45.3587        0.8256
60        39.7552        0
70        26.9176        0
80        18.5399        0
90        8.9469        0
100        1.2214        0
110        1.2586        0
120        1.1816        0
140        1.2251        0
160        1.2177        0
180        1.2535        0
210        1.2146        0
240        1.2054        0
270        1.1833        0

1stOpt非线性拟合问题
微分方程拟合.png
回复此楼

» 本帖@通知

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

847396025

金虫 (小有名气)

你好 我用的是在网上下载的1stopt软件,拟合的时候有参数限制只能拟合4个参数,但是我至少需要拟合6个参数,请问你的软件可以拟合4个以上的参数吗,如果可以是否可以发我一份,谢谢,不胜感激。
简简单单比什么都好
5楼2016-03-01 14:38:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
z和x有函数关系吗?
如果有,dy/dx就不对了,而是偏y偏x了,原方程就不是常微分方程了。
据我了解,在常微分方程拟合的算法中,是会包含龙格-库塔法解常微分方程步骤的,如果问题本身就不是常微分方程,显然拟合无法开展。

如果z和x没有函数关系,仅仅是变的常数项,那试用版的1stOpt做不了的,求助dingd用正版的1stOpt做一下吧。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2016-03-01 10:38:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gxy2011

新虫 (小有名气)

引用回帖:
2楼: Originally posted by 月只蓝 at 2016-03-01 10:38:35
z和x有函数关系吗?
如果有,dy/dx就不对了,而是偏y偏x了,原方程就不是常微分方程了。
据我了解,在常微分方程拟合的算法中,是会包含龙格-库塔法解常微分方程步骤的,如果问题本身就不是常微分方程,显然拟合无 ...

确实z与x,y没有函数关系,仅仅是变的常数项。谢谢您的解答!
3楼2016-03-01 10:46:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gxy2011

新虫 (小有名气)

不好意思,上面第二列是y,第三列是z @月只蓝
4楼2016-03-01 10:51:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:29 by ZPa0EcMwuECS
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:22 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:18 by ZPa0EcMwuECS
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:25 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:14 by ZPa0EcMwuECS
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 00:53 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 00:52 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:41 by ZPa0EcMwuECS
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见