| 查看: 1869 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
自定义函数拟合 请假大侠帮忙！非常非常感激不尽！！！
|
|||
|
因为不知道初值，origin自定义函数怎么也拟合不出来，4stopt拟合还在摸索中，请假哪位大侠帮我拟合一下吧！小女子在此十万分的感谢了！ x,y的数据如下： x=0.16382 0.2513 0.33647 0.36263 0.45197 0.4947 0.53649 0.56201 0.59333 y=420 492 550 570 608 612 620 629 656 自定义函数：y=a0+a1*(x+a2)^a3 |
回复此楼» 收录本帖的淘帖专辑推荐
 SCI论文 |
» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有295人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有1stOpt4.0版本的朋友帮忙拟合一下复数吧 感激不尽呐
已经有4人回复
- Matlab 的自定义函数拟合
已经有3人回复
zzcentaurus
铁虫 (初入文坛)
- 应助: 9 (幼儿园)
- 金币: 48.2
- 帖子: 34
- 在线: 30.5小时
- 虫号: 1450805
- 注册: 2011-10-19
- 性别: GG
- 专业: 粒子物理与核物理实验方法
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
快乐的乐乐: 金币+10, 应该是自定义函数不合适，a3的初值应在0.1-0.3之间，其它初值还确定不了，其实最终目的是根据这些点回归出来一个接近且准确的幂函数方程，现在这个自定义的函数是他们之前回归出来的，看来是完全错误的。请问怎么根据这些数据点回归出来一种幂函数方程？这个你能帮我一下吗？或者能教我一下吗？目前我还只会操作origin,实在是很感谢你！！ 2013-03-28 10:07:45
csgt0: 金币+1, 谢谢 2013-03-28 17:52:51
快乐的乐乐: 金币+10, 应该是自定义函数不合适，a3的初值应在0.1-0.3之间，其它初值还确定不了，其实最终目的是根据这些点回归出来一个接近且准确的幂函数方程，现在这个自定义的函数是他们之前回归出来的，看来是完全错误的。请问怎么根据这些数据点回归出来一种幂函数方程？这个你能帮我一下吗？或者能教我一下吗？目前我还只会操作origin,实在是很感谢你！！ 2013-03-28 10:07:45
csgt0: 金币+1, 谢谢 2013-03-28 17:52:51
|
嗯，而且试了一下，给定不一样的初始值，最后拟合出来的结果也不同，参数的结果差异非常大。。。如果没有一个比较靠谱的初始值的话，最后拟合出来估计也不可信。。。 |
5楼2013-03-27 23:38:59
xujince
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2516.7
- 红花: 1
- 帖子: 548
- 在线: 376.9小时
- 虫号: 1926989
- 注册: 2012-08-06
- 性别: GG
- 专业: 固体力学
3楼2013-03-27 22:22:28
zzcentaurus
铁虫 (初入文坛)
- 应助: 9 (幼儿园)
- 金币: 48.2
- 帖子: 34
- 在线: 30.5小时
- 虫号: 1450805
- 注册: 2011-10-19
- 性别: GG
- 专业: 粒子物理与核物理实验方法
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
快乐的乐乐: 金币+10, ★★★很有帮助, 非常感谢，可能是自定义的函数有问题 2013-03-28 10:02:33
csgt0: 金币+1, 谢谢 2013-03-28 17:52:42
感谢参与，应助指数 +1
快乐的乐乐: 金币+10, ★★★很有帮助, 非常感谢，可能是自定义的函数有问题 2013-03-28 10:02:33
csgt0: 金币+1, 谢谢 2013-03-28 17:52:42
|
刚才用matlab的cftool帮你拟合了一下，感觉按照你的这个函数拟合出来结果不怎么样啊。。。我不会传图片，下面是拟合结果： General model: f(x) = a0+a1*(x+a2)^a3 Coefficients (with 95% confidence bounds): a0 = -65.19 (-2.799e+005, 2.798e+005) a1 = 0.455 (-7803, 7804) a2 = 5.475 (-1.08e+004, 1.081e+004) a3 = 4.084 (-5690, 5698) Goodness of fit: SSE: 3554 R-square: 0.9219 Adjusted R-square: 0.875 RMSE: 26.66 你对这些参数有没有些概念，比如大概是多少，什么范围，一般来说，给出这些之后拟合也才更加准确，反正我觉得我刚才的拟合结果我自己是不能接受的。。。 |
4楼2013-03-27 23:23:19
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
★
csgt0: 金币+1, 谢谢 2013-03-28 17:52:56
csgt0: 金币+1, 谢谢 2013-03-28 17:52:56
|
用1stOpt吧，不需初值，结果稳定唯一： 均方差(RMSE): 6.62844319454137 残差平方和(SSE): 395.426332649356 相关系数(R): 0.995644022086177 相关系数之平方(R^2): 0.991307018715939 决定系数(DC): 0.991307018715939 卡方系数(Chi-Square): 0.320745467120516 F统计(F-Statistic): 191.725555618373 参数 最佳估算 ---------- ------------- a0 772.095936316939 a1 -177.457418189015 a2 0.566934026145924 a3 -2.19432258957609  t3.jpg |
6楼2013-03-28 08:52:19
10