24小时热门版块排行榜    

查看: 2736  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

一直在纠结

铁杆木虫 (正式写手)

[求助] 用MATLAB回归方程的参数,是个隐函数方程 已有2人参与

请大神帮忙回归一下方程的参数:
x=
53.428        54.076        55.046        57.571        78.753        79.154        79.485        80.022        80.488        80.938        81.162        81.387        81.613        82.068        82.296        82.524        83.209        83.436        84.116        84.344        84.809        85.048        85.544        85.989        86.577        87.222        87.987        88.951        89.272        89.832        91.809        101.409        138.539
y=
1.29E+01        1.30E+01        1.32E+01        1.37E+01        1.93E+01        1.93E+01        1.93E+01        1.93E+01        1.93E+01        1.93E+01        1.93E+01        1.92E+01        1.92E+01        1.92E+01        1.92E+01        1.92E+01        1.92E+01        1.91E+01        1.91E+01        1.91E+01        1.91E+01        1.91E+01        1.89E+01        1.85E+01        1.84E+01        1.83E+01        1.83E+01        1.83E+01        1.83E+01        1.84E+01        1.88E+01        2.14E+01        3.36E+01
方程如下:
logy=a1+a2/(a3+x)+a4*x+a5*x^2+a6*logx+a7*y/x^2
其中a1~a7是方程要回归的7个参数
能帮忙回归出来就可以,不限于用MATLAB
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

一直在纠结

铁杆木虫 (正式写手)

引用回帖:
9楼: Originally posted by hc315 at 2014-08-01 14:04:32
将原来的式子两边同时乘以a3+x展开,整理同类项,将合并系数视为整体,进行回归,从而形成新的回归方程。只不过该方程可能无法还原为原方程,但如果你只是要拟合数据的话,那就无所谓了。...

http://muchong.com/bbs/viewthread.php?tid=7590721这是我另外开的一个帖子,里面有个人帮我说了一个方法,但是我照着写了一下有错误,能不能麻烦你给看一下,这个上面忘记说要把程序写上了。。。
10楼2014-08-01 15:27:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
一直在纠结: 金币+10, ★★★很有帮助 2014-06-22 10:10:58
均方差(RMSE): 0.0267854899556039
残差平方和(SSE): 0.0236762615813381
相关系数(R): 0.99996676427131
相关系数之平方(R^2): 0.999933529647235
决定系数(DC): 0.999933528953622
卡方系数(Chi-Square): 0.000664060614024672
F统计(F-Statistic): 466342.333528455

参数        最佳估算
----------        -------------
a1                 -21.123329221464
a2                 -0.0843462681141015
a3                 -75.7548418376161
a4                 -0.0861965264233731
a5                 0.000196616383852714
a6                 6.50999607267736
a7                 402.908092685275
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2014-06-20 22:19:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
一直在纠结: 金币+10, ★★★很有帮助 2014-06-22 10:11:05
隐函数拟合难度一般比较大,下面是用1stOpt求得的一组结果,参考下,有时间多计算几次结果或许会更好。

y=exp(a1+a2/(a3+x)+a4*x+a5*x^2+a6*ln(x)+a7*y/x^2)

目标函数值(最小): 41.5483587726347
均方差(RMSE): 1.12206781545482
残差平方和(SSE): 41.5481940218251
相关系数(R): 0.939896703711614
确定系数(DC): 0.883353547084267

参数        最佳估算
----------        -------------
a1        -4.37936849843139
a2        119.930462705615
a3        2.63450951276076
a4        0.0319253374374684
a5        -9.63612717156336E-5
a6        0.911783693805508
a7        -41.3173340800678
用MATLAB回归方程的参数,是个隐函数方程
c105.jpg

3楼2014-06-20 22:40:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

一直在纠结

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by 月只蓝 at 2014-06-20 22:19:31
均方差(RMSE): 0.0267854899556039
残差平方和(SSE): 0.0236762615813381
相关系数(R): 0.99996676427131
相关系数之平方(R^2): 0.999933529647235
决定系数(DC): 0.999933528953622
卡方系数(Chi-Square): 0. ...

能麻烦您把写的程序发上来吗?  忘记写上了。。。
4楼2014-06-21 08:11:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[文学芳草园] 梦想 +5 myrtle 2026-08-26 7/350 2026-09-01 15:18 by myrtle
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 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
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见