24小时热门版块排行榜    

查看: 1961  |  回复: 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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

★
csgt0: 金币+1, 谢谢 2013-03-28 17:52:56
用1stOpt吧,不需初值,结果稳定唯一:
CODE:
Function y=a0+a1*(x+a2)^a3;
Data;
0.16382        420
0.2513        492
0.33647        550
0.36263        570
0.45197        608
0.4947        612
0.53649        620
0.56201        629
0.59333        656

均方差(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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

xujince

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
试一下对自定义函数两边取lg。
3楼2013-03-27 22:22:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzcentaurus

铁虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzcentaurus

铁虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
快乐的乐乐: 金币+10, 应该是自定义函数不合适,a3的初值应在0.1-0.3之间,其它初值还确定不了,其实最终目的是根据这些点回归出来一个接近且准确的幂函数方程,现在这个自定义的函数是他们之前回归出来的,看来是完全错误的。请问怎么根据这些数据点回归出来一种幂函数方程?这个你能帮我一下吗?或者能教我一下吗?目前我还只会操作origin,实在是很感谢你!! 2013-03-28 10:07:45
csgt0: 金币+1, 谢谢 2013-03-28 17:52:51
引用回帖:
4楼: Originally posted by zzcentaurus at 2013-03-27 23:23:19
刚才用matlab的cftool帮你拟合了一下,感觉按照你的这个函数拟合出来结果不怎么样啊。。。我不会传图片,下面是拟合结果:
General model:
     f(x) = a0+a1*(x+a2)^a3
Coefficients (with 95% confidence boun ...

嗯,而且试了一下,给定不一样的初始值,最后拟合出来的结果也不同,参数的结果差异非常大。。。如果没有一个比较靠谱的初始值的话,最后拟合出来估计也不可信。。。
5楼2013-03-27 23:38:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 12/600 2026-10-13 20:40 by Merryya
[硕博家园] 极限拉扯,如何抉择 +9 Isak0320 2026-10-12 13/650 2026-10-13 11:57 by Isak0320
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见