24小时热门版块排行榜    

查看: 1102  |  回复: 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的回帖

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的回帖
查看全部 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

新虫 (小有名气)

不好意思,上面第二列是y,第三列是z @月只蓝
4楼2016-03-01 10:51:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

847396025

金虫 (小有名气)

你好 我用的是在网上下载的1stopt软件,拟合的时候有参数限制只能拟合4个参数,但是我至少需要拟合6个参数,请问你的软件可以拟合4个以上的参数吗,如果可以是否可以发我一份,谢谢,不胜感激。
简简单单比什么都好
5楼2016-03-01 14:38:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:45 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 01:14 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:14 by ZPa0EcMwuECS
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:12 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 01:04 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 01:03 by ZPa0EcMwuECS
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 01:03 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:52 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 00:31 by ZPa0EcMwuECS
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 jCd0dEvKHShX 2026-08-29 3/150 2026-08-29 23:14 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G6APbkg8SA6w 2026-08-29 4/200 2026-08-29 22:27 by ZPa0EcMwuECS
[基金申请] 面上意见出来了 +8 黄鸟于飞Chao 2026-08-29 16/800 2026-08-29 20:45 by 黄鸟于飞Chao
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见