| 查看: 1177 | 回复: 3 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
forestsong12银虫 (小有名气)
|
[求助]
1stopt软件拟合非线性方程，计算结果全为0，这是为什么 已有2人参与
|
|||
|
Parameters T; Variable x,y; Function y= 3.741832*10^-16/(exp(0.01438786/(T*x*10^-3))-1)/(x*10^-3)^5; Data; 405.425 19197 424.055 14239 441.669 11853 463.348 6832.2 541.255 4912.8 566.321 3806.9 609.001 2420.1 711.297 1356.1 744.153 810.13 764.477 613.14 804.786 412.09 852.208 315.52 900.307 463.06 迭代数: 33 计算用时(时:分:秒:微秒): 00:00:00:319 优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法 计算结束原因: 达到收敛判断标准 均方差(RMSE): 7875.40731737143 残差平方和(SSE): 806286525.388597 相关系数(R): 0.985793202644828 相关系数之平方(R^2): 0.971788238380747 决定系数(DC): -0.75817360362487 卡方系数(Chi-Square): 67231.0399999996 F统计(F-Statistic): 378.908299540322 参数 最佳估算 ---------- ------------- T 27.2788094020072 ====== 结果输出 ===== No 实测值y 计算值y 1 19197 0.0000000 2 14239 0.0000000 3 11853 0.0000000 4 6832.2 0.0000000 5 4912.8 0.0000000 6 3806.9 0.0000000 7 2420.1 0.0000000 8 1356.1 0.0000000 9 810.13 0.0000000 10 613.14 0.0000000 11 412.09 0.0000000 12 315.52 0.0000000 13 463.06 0.0000000 |
回复此楼» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有116人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有29人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有53人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 1stopt拟合的动力学参数相关系数不高，有哪些原因
已经有3人回复
- 1stopt拟合出参数为0，为什么？
已经有6人回复
- 求高手给做一道题，用MATLAB或1stOpt，请把结果直接给出来，我的软件就不行，谢谢
已经有5人回复
- 求大神帮忙拟合一个非线性方程，求出模型参数
已经有15人回复
- 为什么同样的设置，计算出的结果有时会相差很大呢？
已经有8人回复
- 用1stOpt进行复数曲线拟合时，结果出错。急求教！！！
已经有5人回复
- 用1stOpt进行多模型多变量非线性拟合得不到结果
已经有4人回复
- 1stopt拟合问题求助
已经有7人回复
- 采用1stopt拟合
已经有8人回复
- 求大神帮忙用1stOpt 拟合复数数据
已经有13人回复
- 紧急求助：求前辈帮忙用1stOpt拟合一下数据
已经有3人回复
- 1stopt中使用sharedmodel拟合曲线出现问题，求指示！！
已经有6人回复
- 用1stopt拟合非线性方程结果与其他软件拟合结果差异大
已经有5人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- 用1stopt或者origin非线性拟合磁导率
已经有5人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
forestsong12
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 120.3
- 散金: 36
- 红花: 1
- 帖子: 168
- 在线: 53.2小时
- 虫号: 2772604
- 注册: 2013-11-02
- 专业: 光学
3楼2014-07-07 09:48:43
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2014-07-07 09:39:07
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
感谢参与，应助指数 +1
|
公式本身不合适。 先来看你的数据 x y 405.425 19197 424.055 14239 441.669 11853 463.348 6832.2 541.255 4912.8 566.321 3806.9 609.001 2420.1 711.297 1356.1 744.153 810.13 764.477 613.14 804.786 412.09 852.208 315.52 900.307 463.06 你给出方程简化得： y=3.741832*0.1*(x^(-5)) / (exp(14.38786/T/x)-1) 注意x是几百数量级的，那么式中3.741832*0.1*(x^(-5))将是非常小的数，至少小于10^-11，4^5与9^5相差50倍；而y是数百到几万的数值，那么显然 exp(14.38786/T/x)-1 这一项应当是非常接近于0的数，也就是说 exp(14.38786/T/x)这一项应当非常接近于1，这要求T非常大；如前述，3.741832*0.1*(x^(-5))在x整个变化范围中数值变化达50倍，而至于指数函数而言，其变化范围是很窄的。 也就是，不可能存在合适的T值能满足高拟合的质量。 至于你说y的计算值为0的问题，在1stopt算法设置中，以简面体爬山法计算，可得： 相关系数(R): 0.985799655457293 相关系数之平方(R^2): 0.971800960699717 决定系数(DC): -0.758172807118177 卡方系数(Chi-Square): 67230.9893136268 F统计(F-Statistic): 379.084211127347 参数 最佳估算 ---------- ------------- T 117315023909293 ====== 结果输出 ===== No 实测值y 计算值y 1 19197 0.0038974 2 14239 0.0032564 3 11853 0.0027671 4 6832.2 0.0022845 5 4912.8 0.0012269 6 3806.9 0.0010237 7 2420.1 0.0007655 8 1356.1 0.0004113 9 810.13 0.0003434 10 613.14 0.0003083 11 412.09 0.0002510 12 315.52 0.0001996 13 463.06 0.0001603 从结果中可见，T的估算值非常大。 简而概之，就是公式不合适。此外，在拟合的时候，公式输入能化简就化简，有利于避免诸如过拟合之类的问题，从而得到高质量的拟合结果。 |
4楼2014-07-07 18:47:01
