应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 用MATLAB回归方程的参数，是个隐函数方程
51/1返回列表
查看: 2711  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

一直在纠结

铁杆木虫 (正式写手)

[求助] 用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朵）

笑靥如花?

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-06-20 21:25:28
已阅   回复此楼   关注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的回帖
查看全部 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
赞一下(1人) 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2014-06-20 22:19:31
已阅   回复此楼   关注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的回帖

一直在纠结

铁杆木虫 (正式写手)
引用回帖:
3楼: Originally posted by dingd at 2014-06-20 22:40:40
隐函数拟合难度一般比较大，下面是用1stOpt求得的一组结果，参考下，有时间多计算几次结果或许会更好。

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

目标函数值(最小): 41.5483587726347
均方差(RMS ...

呃。。跟上面的回归结果差距好大啊，这个还真是不好弄啊，辛苦您了，能麻烦您把写的程序发一下吗？
赞一下 回复此楼
5楼2014-06-21 08:12:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +9 晴天加油 2026-08-26 10/500 2026-08-27 16:20 by 443933501
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +6 苏东坡二世 2026-08-23 6/300 2026-08-27 13:57 by zhangjianl
[基金申请] 青基如何破局 +3 lyz123lyz 2026-08-27 4/200 2026-08-27 13:08 by kulium
[基金申请] 基金不中，共勉 +8 eulota 2026-08-26 8/400 2026-08-27 12:34 by lice_1987
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 17/850 2026-08-27 08:20 by 大萍1987
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 国际合作可查了，中了面上 (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
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭