24小时热门版块排行榜    

查看: 1874  |  回复: 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论文

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国社科又开始会评了,不知道这次命运如何 +3 雨打竹帘 2026-08-30 4/200 2026-08-30 12:35 by 白大多亮
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:14 by l0VvVHGBGRLv
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 12:13 by l0VvVHGBGRLv
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:41 by l0VvVHGBGRLv
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 11:32 by l0VvVHGBGRLv
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:30 by l0VvVHGBGRLv
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 17/850 2026-08-30 10:28 by 孤独的英雄6
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 06:55 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 06:22 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 06:11 by ZPa0EcMwuECS
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 6/300 2026-08-29 18:10 by manplx
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见