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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 工科数学 » 非线性拟合
41/1返回列表
查看: 872  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

RABBIT

铜虫 (小有名气)

[求助] 非线性拟合

小白有四列数据

x1        x2        x3        y
40        0.7        45.47         3.87
40        0.7        76.62         6.37
40        0.7        106.28         8.28
40        0.7        135.67         9.23
40        0.7        235.56         9.90
40        0.7        268.20         11.03
40        0.7        384.58         11.79
40        0.7        503.57         12.72
40        0.7        586.93         13.39
40        0.7        682.19         14.58
40        0.7        752.19         15.60
40        0.8        45.81         3.00
40        0.8        73.78         5.38
40        0.8        94.77         7.02
40        0.8        135.19         8.54
40        0.8        214.05         9.09
40        0.8        331.18         10.09
40        0.8        488.09         10.27
40        0.8        519.04         10.71
40        0.8        561.47         11.42
40        0.8        695.22         12.36
40        0.8        824.44         14.21
40        0.9        46.18         2.67
40        0.9        69.34         4.71
40        0.9        108.64         6.64
40        0.9        132.45         7.26
40        0.9        241.82         7.94
40        0.9        294.09         8.70
40        0.9        491.76         9.39
40        0.9        603.92         10.12
40        0.9        729.84         11.02
40        0.9        825.79         11.88
40        0.9        900.19         13.04
50        0.7        45.47         3.13
50        0.7        75.17         4.73
50        0.7        106.28         6.41
50        0.7        142.14         8.34
50        0.7        227.50         9.07
50        0.7        302.03         10.00
50        0.7        402.87         10.55
50        0.7        449.17         11.00
50        0.7        545.41         11.33
50        0.7        586.38         12.00
50        0.7        634.31         13.21
50        0.8        45.81         2.47
50        0.8        72.99         4.03
50        0.8        106.36         6.10
50        0.8        144.29         7.43
50        0.8        229.08         8.13
50        0.8        310.76         8.98
50        0.8        405.17         9.33
50        0.8        468.47         9.90
50        0.8        557.94         11.00
50        0.8        594.68         11.30
50        0.8        647.29         11.87
50        0.9        46.18         1.67
50        0.9        77.84         3.34
50        0.9        109.13         5.19
50        0.9        138.23         6.00
50        0.9        235.56         6.67
50        0.9        313.86         7.67
50        0.9        440.95         8.53
50        0.9        486.01         8.90
50        0.9        591.26         9.35
50        0.9        606.10         9.80
50        0.9        674.13         10.65
60        0.7        45.09         2.67
60        0.7        78.39         3.00
60        0.7        104.07         4.37
60        0.7        136.89         6.67
60        0.7        239.04         8.67
60        0.7        299.79         9.00
60        0.7        371.65         10.33
60        0.7        440.04         10.47
60        0.7        558.66         11.00
60        0.7        621.48         11.79
60        0.7        671.71         13.14
60        0.8        45.68         2.00
60        0.8        79.90         2.67
60        0.8        104.94         4.67
60        0.8        134.02         5.33
60        0.8        272.64         6.33
60        0.8        369.07         7.33
60        0.8        464.95         8.42
60        0.8        501.91         8.83
60        0.8        532.78         10.05
60        0.8        610.20         10.95
60        0.8        681.87         11.76
60        0.9        46.39         1.32
60        0.9        80.95         1.67
60        0.9        105.59         3.24
60        0.9        137.26         3.67
60        0.9        260.49         5.00
60        0.9        376.01         6.47
60        0.9        446.68         8.03
60        0.9        533.76         8.33
60        0.9        637.88         9.00
60        0.9        690.83         9.31
60        0.9        703.69         10.33
x1、x2、x3、y
单独分析满足如下关系
y＝a1*x1^b1,
y=a1*x2^b2
y=z3*lnx3+b3
经过数据整理，y与x1 x2 x3之间应当满足如下关系
lny=Aln（x1）+Bln（x2）+1/3*ln(C*ln（x3）+D)
excel  我只会二元拟合，
matlab 小白我只会regress，
origin 1stopt 暂时没有学会

对于上述的四元简单非线性拟合，小白无助了。

请 各位大神多多帮助！！
  真诚希望 各位帮忙！！

ps：希望各位能够告知拟合后结果
pps：最好能够告诉拟合方法，因为我还有好多数据和公式模型需要拟合。

再次  真诚感谢，
  期待帮助
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 四列数据.xlsx
    • 2013-04-15 09:18:38, 13.05 K

» 收录本帖的淘帖专辑推荐

数据处理资源

» 猜你喜欢
1楼 2013-04-15 09:23:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
1stOpt代码：
CODE:
Function y=exp(A*ln(x1)+B*ln(x2)+1/3*ln(C*ln(x3)+D));
Data;
x1        x2        x3        y
40        0.7        45.47        3.87
40        0.7        76.62        6.37
40        0.7        106.28        8.28
40        0.7        135.67        9.23
40        0.7        235.56        9.90
40        0.7        268.20        11.03
40        0.7        384.58        11.79
40        0.7        503.57        12.72
40        0.7        586.93        13.39
40        0.7        682.19        14.58
40        0.7        752.19        15.60
40        0.8        45.81        3.00
40        0.8        73.78        5.38
40        0.8        94.77        7.02
40        0.8        135.19        8.54
40        0.8        214.05        9.09
40        0.8        331.18        10.09
40        0.8        488.09        10.27
40        0.8        519.04        10.71
40        0.8        561.47        11.42
40        0.8        695.22        12.36
40        0.8        824.44        14.21
40        0.9        46.18        2.67
40        0.9        69.34        4.71
40        0.9        108.64        6.64
40        0.9        132.45        7.26
40        0.9        241.82        7.94
40        0.9        294.09        8.70
40        0.9        491.76        9.39
40        0.9        603.92        10.12
40        0.9        729.84        11.02
40        0.9        825.79        11.88
40        0.9        900.19        13.04
50        0.7        45.47        3.13
50        0.7        75.17        4.73
50        0.7        106.28        6.41
50        0.7        142.14        8.34
50        0.7        227.50        9.07
50        0.7        302.03        10.00
50        0.7        402.87        10.55
50        0.7        449.17        11.00
50        0.7        545.41        11.33
50        0.7        586.38        12.00
50        0.7        634.31        13.21
50        0.8        45.81        2.47
50        0.8        72.99        4.03
50        0.8        106.36        6.10
50        0.8        144.29        7.43
50        0.8        229.08        8.13
50        0.8        310.76        8.98
50        0.8        405.17        9.33
50        0.8        468.47        9.90
50        0.8        557.94        11.00
50        0.8        594.68        11.30
50        0.8        647.29        11.87
50        0.9        46.18        1.67
50        0.9        77.84        3.34
50        0.9        109.13        5.19
50        0.9        138.23        6.00
50        0.9        235.56        6.67
50        0.9        313.86        7.67
50        0.9        440.95        8.53
50        0.9        486.01        8.90
50        0.9        591.26        9.35
50        0.9        606.10        9.80
50        0.9        674.13        10.65
60        0.7        45.09        2.67
60        0.7        78.39        3.00
60        0.7        104.07        4.37
60        0.7        136.89        6.67
60        0.7        239.04        8.67
60        0.7        299.79        9.00
60        0.7        371.65        10.33
60        0.7        440.04        10.47
60        0.7        558.66        11.00
60        0.7        621.48        11.79
60        0.7        671.71        13.14
60        0.8        45.68        2.00
60        0.8        79.90        2.67
60        0.8        104.94        4.67
60        0.8        134.02        5.33
60        0.8        272.64        6.33
60        0.8        369.07        7.33
60        0.8        464.95        8.42
60        0.8        501.91        8.83
60        0.8        532.78        10.05
60        0.8        610.20        10.95
60        0.8        681.87        11.76
60        0.9        46.39        1.32
60        0.9        80.95        1.67
60        0.9        105.59        3.24
60        0.9        137.26        3.67
60        0.9        260.49        5.00
60        0.9        376.01        6.47
60        0.9        446.68        8.03
60        0.9        533.76        8.33
60        0.9        637.88        9.00
60        0.9        690.83        9.31
60        0.9        703.69        10.33

均方差(RMSE): 1.26174829708875
残差平方和(SSE): 157.60886775543
相关系数(R): 0.934897845124764
相关系数之平方(R^2): 0.874033980818927
决定系数(DC): 0.856234052147805
卡方系数(Chi-Square): 12.4338955887395
F统计(F-Statistic): 220.265418478811

参数        最佳估算
----------        -------------
a        -0.627484130653009
b        -1.02179742605419
c        305413.535181483
d        -1158318.09682216

t1.jpg
赞一下(1人) 回复此楼
2楼2013-04-15 09:59:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

RABBIT

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2013-04-15 09:59:17
1stOpt代码：

Function y=exp(A*ln(x1)+B*ln(x2)+1/3*ln(C*ln(x3)+D));
Data;
x1        x2        x3        y
40        0.7        45.47        3.87
40        0.7        76.62        6.37
40        0.7        106.28        8.28
40        0.7        135.67        9.23
40        0.7        235.56        9.90
40        0.7        268. ...

dingd
非常感谢！
我使用你帮忙拟合的 相关式计算后发现，平均相对偏差71.5%呢！
你能帮忙看一下吗？
理论上讲，你帮忙拟合的这一相关式更合理（因为y与x3之间明显为对数关系）
但是，这一公式虽然合理，却相对偏差还不如
y=a1* x1^b1*x2^b2*x3^b3

这其中原因可能是什么？
你能帮忙分析一下吗？

可能表达得不好，我再描述一下这个问题：
这个问题是这样的：
对于y,x1,x2,x3
有两种拟合方式：
1、y=a1*x1^b1*x2^b2*x3^b3
2、y=a1*x1^b1*x2^b2*(lnx3+b3)
使用方式1，拟合较好，但不合理；
使用方式2，拟合不好，但更合理。
这是为什么？

谢谢
赞一下 回复此楼
3楼2013-04-15 11:18:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
RABBIT: 金币+15, ★★★★★最佳答案, 非常感谢！！ 2013-04-15 15:12:33
方法1：

均方差(RMSE): 0.858822599150106
残差平方和(SSE): 73.0200494242834
相关系数(R): 0.966657040009446
相关系数之平方(R^2): 0.934425832999825
决定系数(DC): 0.9333933631578
卡方系数(Chi-Square): 7.35569053843467
F统计(F-Statistic): 475.428097979022

参数        最佳估算
----------        -------------
a1        4.31783270994662
b1        -0.52270862849314
b2        -1.14363688062844
b3        0.425888718982471

方法2：
均方差(RMSE): 0.754727096141286
残差平方和(SSE): 56.391685975336
相关系数(R): 0.974048757610354
相关系数之平方(R^2): 0.948770982202275
决定系数(DC): 0.948561243407903
卡方系数(Chi-Square): 4.76101947178774
F统计(F-Statistic): 615.618820606017

参数        最佳估算
----------        -------------
a1        22.3458878154313
b1        -0.560984296857699
b2        -1.08987846862157
b3        -2.97267635021781

很明显，方法2效果更好啊！
赞一下 回复此楼
4楼2013-04-15 11:54:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 RABBIT 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 基金不中，共勉 +5 eulota 2026-08-26 5/250 2026-08-26 23:33 by young90
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国际合作可查了，中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭