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

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +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的回帖

847396025

金虫 (小有名气)

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

gxy2011

新虫 (小有名气)

引用回帖:
5楼: Originally posted by 847396025 at 2016-03-01 14:38:18
你好 我用的是在网上下载的1stopt软件,拟合的时候有参数限制只能拟合4个参数,但是我至少需要拟合6个参数,请问你的软件可以拟合4个以上的参数吗,如果可以是否可以发我一份,谢谢,不胜感激。

你好,我这个最多也只能拟合4个参数
6楼2016-03-01 14:53:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gxy2011 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 帮忙看看fileCode +4 wwncly 2026-08-10 8/400 2026-08-10 22:07 by demonstreets
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +3 Tide man 2026-08-10 3/150 2026-08-10 21:48 by GZUZM
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 15/750 2026-08-10 19:29 by 要乐观耀哥
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +4 布布和一二 2026-08-10 4/200 2026-08-10 14:50 by 医学老男孩
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
信息提示
请填处理意见