24小时热门版块排行榜    

查看: 1910  |  回复: 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的回帖

xujince

木虫 (正式写手)

【答案】应助回帖

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

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-18 10:30 by yS4CkeWDhD6V
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 s3fFTmArrBt6 2026-09-14 6/300 2026-09-18 10:16 by yS4CkeWDhD6V
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 6/300 2026-09-18 09:18 by yS4CkeWDhD6V
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 C79jjtjAKjEn 2026-09-15 4/200 2026-09-18 06:44 by manzzl
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 C79jjtjAKjEn 2026-09-15 3/150 2026-09-18 03:42 by cITaGg3p5edr
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 5BDX0d0WFp7t 2026-09-15 3/150 2026-09-18 03:16 by cITaGg3p5edr
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-18 01:04 by cITaGg3p5edr
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-17 22:39 by cuPZTDXS3VOd
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-17 14:27 by ObVzqQQh1rrE
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 0pYnUiPDfdnk 2026-09-14 4/200 2026-09-17 14:15 by ObVzqQQh1rrE
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-17 13:51 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-17 13:15 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 mibUvS8DDCwf 2026-09-15 4/200 2026-09-17 06:26 by k06BKNblrGRH
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-17 03:14 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-17 02:38 by Tql5LQhh5rLK
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-17 01:16 by Tql5LQhh5rLK
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 16:00 by i3KBRD19ERkG
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 QUjhNVAcOSff 2026-09-13 7/350 2026-09-16 14:36 by i3KBRD19ERkG
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
信息提示
请填处理意见