| 查看: 2216 | 回复: 8 | |||
[交流]
个人发现软件1stopt的隐函数拟合功能存在严重问题 已有4人参与
|
|
关于隐函数公式的拟合，大家都说用1stopt来处理。本人发现1stopt的所谓隐函数拟合功能好像存在严重问题。 我们看一段简单的代码： Title "Type your title here"; Parameters a; Variable x,y; Function y=a*x-y; Data; 1 1.1 2 2.4 3 3.1 4 3.8 5 5.3 这个函数我们认为是隐函数，这个写法也就是隐函数拟合的写法。 结果如下： 参数 最佳估算 ---------- ------------- a 2.06909090943816 ====== 结果输出 ===== No 实测值y 计算值y 1 1.1 0.9690909 2 2.4 1.7381818 3 3.1 3.1072727 4 3.8 4.4763636 5 5.3 5.0454545 OK，可能很多人到此就认为参数a的估计值已经得到了，没问题了。但是我对软件的解隐式方程存疑。我们来验算一下。 对于实测值y=1.1这个点，拟合公式计算值是0.9690909，但是请用手算一下呢？y=a*x-y，令x=1，a=2.06909090943816，得到的应该是y=1.0345 那么这个“计算值”y是怎么算出来的？其实是直接将(x,y)=(1,1.1)给代入拟合公式的右边，y=a*x-y=2.06909090943816*1-1.1=0.9690909 所以，1stopt根本就没有进行隐式方程的求解！所谓的隐函数公式的拟合，里面所包含的解隐式方程得到计算值这一步都是有问题的，那么其拟合结果也是不正确的。 希望用过该功能的虫友能交流一下，1stopt到底能否做隐函数的拟合？ |
回复此楼» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有178人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 使用AdNDP方法以及ELF/LOL、多中心键级研究多中心键
已经有14人回复
- 回答一些关于Multiwfn的疑问以及未来Multiwfn的发展打算
已经有37人回复
- 这些网站你要是没进过？那你就Out了
已经有2人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2015-09-28 17:24:09
3楼2015-09-28 17:26:09
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2015-09-28 17:28:37
5楼2015-09-28 17:31:33
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
一楼代码6.0计算： 目标函数值(最小): 0.244363635752727 均方差(RMSE): 0.221071768953803 残差平方和(SSE): 0.244363635141818 相关系数(R): 0.989343958232479 确定系数(DC): 0.975095430580736 参数 最佳估算 ---------- ------------- a 2.06909090871176 ====== 结果输出 ====== 目标 y 计算 y 1.1 1.03454545451952 2.4 2.06909090953904 3.1 3.10363636305855 3.8 4.13818181657807 5.3 5.17272727209759 1.5版的，估计还是破解的，有些Out了吧，尤其是隐函数拟合，那个1.5版本本身就支持的不好。 |
6楼2015-09-28 20:01:19
7楼2015-09-28 20:19:14
nankailiang
铁虫 (正式写手)
- 应助: 24 (小学生)
- 金币: 185.1
- 散金: 91
- 红花: 5
- 帖子: 812
- 在线: 402.6小时
- 虫号: 1294784
- 注册: 2011-05-13
- 专业: 无机材料化学
8楼2015-09-30 16:53:54
皓小天
木虫之王 (职业作家)
- 应助: 19 (小学生)
- 贵宾: 0.1
- 金币: 76016.6
- 散金: 502
- 红花: 137
- 帖子: 3801
- 在线: 869小时
- 虫号: 1864745
- 注册: 2012-06-19
- 性别: GG
- 专业: 力学中的基本问题和方法
9楼2015-09-30 21:46:39
