| 查看: 998 | 回复: 9 | ||||
水中的梦88木虫 (正式写手)
|
[求助]
1stOpt拟合没有反应，迭代数为0，求大神帮忙看看已有1人参与
|
|
如题： 代码如下： Parameters De0<1e-8,kr>0; Variables t,y; OdeFunction y'=kr*3*(1-y)^(2/3)/(2-2*De0*ln(1/y)*(1-y)^(1/3)); Data; 0 0.448761513 50 0.468487883 100 0.481780381 150 0.493508003 187.5 0.499387463 求大神帮忙看看。以前用1.5版本的算过没有问题，现在下载的好像是5.0版本，算了之后就没反应。 |
回复此楼» 猜你喜欢
法国博士招聘
已经有0人回复
-
法国博士招聘
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有284人回复
-
法国博士招聘
已经有0人回复
-
法国博士招聘
已经有0人回复
-
超导三元凸包图
已经有0人回复
-
最近遇到这样一个问题
已经有10人回复
-
论文接收，散金币了！
已经有175人回复
-
vasp算居里温度
已经有0人回复
-
郑州大学光电器件方向招收一名博士
已经有6人回复
2楼2018-01-13 20:21:27
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15356.8
- 散金: 101
- 红花: 234
- 帖子: 3409
- 在线: 1222.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
水中的梦88: 金币+10 2018-01-14 14:17:00
感谢参与，应助指数 +1
水中的梦88: 金币+10 2018-01-14 14:17:00
|
正版7.0计算没任何问题： 均方差(RMSE): 0.00357407148406669 残差平方和(SSE): 5.10959478928747E-5 相关系数(R): 0.995876310044685 相关系数之平方(R^2): 0.991769624908217 修正R平方(Adj. R^2): 0.975308874724651 确定系数(DC): 0.90874075449774 F统计(F-Statistic): 35.1948517749889 参数 最佳估算 -------------------- ------------- de0 9.99995331641504E-9 kr 0.000298523773919182 ====== 结果输出 ====== 文件: 数据文件-1 No t 目标 y 计算 y 1 50 0.468487883 0.46367711931354 2 100 0.481780381 0.47832120807955 3 150 0.493508003 0.492696273365113 4 187.5 0.499387463 0.503302526007435 |
» 本帖已获得的红花（最新10朵）
水中的梦883楼2018-01-13 20:22:12
水中的梦88
木虫 (正式写手)
- 应助: 52 (初中生)
- 金币: 2378.9
- 散金: 26
- 红花: 5
- 帖子: 323
- 在线: 195.7小时
- 虫号: 1308847
- 注册: 2011-05-28
- 性别: MM
- 专业: 结构工程
4楼2018-01-13 20:23:50
水中的梦88
木虫 (正式写手)
- 应助: 52 (初中生)
- 金币: 2378.9
- 散金: 26
- 红花: 5
- 帖子: 323
- 在线: 195.7小时
- 虫号: 1308847
- 注册: 2011-05-28
- 性别: MM
- 专业: 结构工程
5楼2018-01-14 14:16:53
lipenggg
专家顾问 (正式写手)
-
专家经验: +65 - 应助: 122 (高中生)
- 金币: 4841.3
- 红花: 43
- 帖子: 722
- 在线: 435.7小时
- 虫号: 3929201
- 注册: 2015-06-17
- 性别: GG
- 专业: 流体力学
- 管辖: 计算模拟
6楼2018-01-15 09:22:18
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 67794.9
- 散金: 1938
- 红花: 442
- 沙发: 4
- 帖子: 4373
- 在线: 3284.9小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
复制你的代码，一点不改，用5.5试用版计算，尽管提示有语法错误，但是仍能计算，并得到跟3楼7.0正式版非常接近的结果，如下图： 用5.5试用版的时候，parameters的约束，用以下的写法更为兼容： Parameters De0[,1e-8],kr[0,]; Variables t,y; OdeFunction y'=kr*3*(1-y)^(2/3)/(2-2*De0*ln(1/y)*(1-y)^(1/3)); Data; 0 0.448761513 50 0.468487883 100 0.481780381 150 0.493508003 187.5 0.499387463 此外你说“以前用1.5版本的算过没有问题”，这是不可能的，因为1.5版本压根就没有常微分方程拟合的功能，不支持关键词OdeFunction。 |
» 本帖已获得的红花（最新10朵）
7楼2018-01-15 10:01:12
水中的梦88
木虫 (正式写手)
- 应助: 52 (初中生)
- 金币: 2378.9
- 散金: 26
- 红花: 5
- 帖子: 323
- 在线: 195.7小时
- 虫号: 1308847
- 注册: 2011-05-28
- 性别: MM
- 专业: 结构工程
8楼2018-01-15 10:53:54
水中的梦88
木虫 (正式写手)
- 应助: 52 (初中生)
- 金币: 2378.9
- 散金: 26
- 红花: 5
- 帖子: 323
- 在线: 195.7小时
- 虫号: 1308847
- 注册: 2011-05-28
- 性别: MM
- 专业: 结构工程
9楼2018-01-15 10:54:08
追星的吐司
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1563.4
- 散金: 15
- 红花: 2
- 帖子: 214
- 在线: 145.8小时
- 虫号: 3777120
- 注册: 2015-03-31
- 专业: 法学其他学科
10楼2018-01-16 10:49:05
2