小木虫

当前位置：首页 > 计算模拟 >求助，1stopt拟合微分方程图像结果不对

求助，1stopt拟合微分方程图像结果不对

作者 考考考考: 来源: 小木虫 500 10 举报帖子

用的1stopt 5.5试用版想要出的图像是以t，s为横纵轴，但实际出来的横轴是1,2,3,4,5,6
这是我写的代码
parameters a,b,c;
variable s,x,t;
odefunction s'=-a*s*x/(b+s+(651.42-s)^2/c);
data;
// s x t
658 0 0
651.42 0 48
625.1 15 72
559.3 45 84
460.6 90 96
282.94 171 108
72.38 267 120
0 300 132
请师兄师姐们帮忙看一下问题出在哪里了。。。谢谢！

结果.png

方程.png 返回小木虫查看更多

今日热帖

精华评论

独孤神宇

横坐标表示数据的个数，这个效果很好了。计算完了后将输出的数据用其他软件绘图。
dingd

均方差(RMSE): 3.73781051354546
残差平方和(SSR): 97.7985920461968
相关系数(R): 0.999889749386679
相关系数之平方(R^2): 0.999779510928556
修正R平方(Adj. R^2): 0.999669266392834
确定系数(DC): 0.999767604712921
F统计(F-Statistic): 8638.03334493322

参数最佳估算
-------------------- -------------
a 0.232564800688783
b 745.118198663314
c -706.354495111857
考考考考

引用回帖:
2楼: Originally posted by 独孤神宇 at 2019-02-25 11:52:32
横坐标表示数据的个数，这个效果很好了。计算完了后将输出的数据用其他软件绘图。

感谢老师回复，但输出数据也是错的，应该所有参数都是正值，但拟合的c值是负的。。。。
还有老师，我想把横坐标变成t应该怎么操作额
考考考考

引用回帖:
3楼: Originally posted by dingd at 2019-02-25 13:58:52
均方差(RMSE): 3.73781051354546
残差平方和(SSR): 97.7985920461968
相关系数(R): 0.999889749386679
相关系数之平方(R^2): 0.999779510928556
修正R平方(Adj. R^2): 0.999669266392834
确定系数(DC): 0.9997 ...

谢谢老师回复，我的拟合结果和您的结果差不多，但c的值应该是正值，所以请问是不是我的代码有哪里出错了。。。
考考考考

引用回帖:
6楼: Originally posted by 独孤神宇 at 2019-02-25 16:50:41
那就约束参数的范围
加一句 ParameterDomain=；
...

老师您好，我按您说的加了那一句，发现结果还是不对，值太大了，而且我发现软件对s和t分别做了两个曲线，不知道是怎么回事。。。我想要的是s和t的关系图。。。

1.png

2.png
dingd

版本问题吧。下面参考：
均方差(RMSE): 4.82770003217269
残差平方和(SSR): 163.146813204482
相关系数(R): 0.999846823935424
相关系数之平方(R^2): 0.999693671333755
修正R平方(Adj. R^2): 0.999540507000632
确定系数(DC): 0.999612320078465
F统计(F-Statistic): 5190.56058397433

参数最佳估算
-------------------- -------------
a 0.123246443246998
b 138.260772724859
c 2.39655917250376E20

====== 结果输出 ======

文件: 数据文件-1
No t 目标 s 计算 s
1 48 651.42 658
2 72 625.1 621.544225367598
3 84 559.3 567.579784257519
4 96 460.6 462.738957222542
5 108 282.94 279.544586815861
6 120 72.38 72.0786281228572
7 132 0 4.7389090429929，