24小时热门版块排行榜    

查看: 2332  |  回复: 14
【悬赏金币】回答本帖问题,作者figo2015将赠送您 500 个金币

figo2015

木虫 (正式写手)

[求助] 求助:多变量曲线拟合已有3人参与

各路大神:小弟遇到一个多元曲线拟合的问题,请大神们多多指教!5个自变量,1个因变量,希望通过多变量拟合出经验公式和图像。求助有拟合过程的教程。
编号        x1        x2        x3        x4        x5        y
1        3        1        3        0.15        0.10        100
2        5        1        3        0.15        0.10        97
3        7        1        3        0.15        0.10        75
4        9        1        3        0.15        0.10        50
5        7        1.5        3        0.15        0.10        55
6        7        2        3        0.15        0.10        36
7        7        2.6        3        0.15        0.10        52
8        7        1        2        0.15        0.10        38
9        7        1        3        0.15        0.10        65
10        7        1        4        0.15        0.10        62
11        7        1        6        0.15        0.10        41
12        4        1        2        0.10        0.10        124
13        4        1        2        0.15        0.10        102
14        4        1        2        0.20        0.10        55
15        4        1        2        0.25        0.10        37
16        4        1        2        0.15        0.05        120
17        4        1        2        0.15        0.15        79
18        4        1        2        0.15        0.20        60
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

shikang999

新虫 (小有名气)

【答案】应助回帖

引用回帖:
8楼: Originally posted by figo2015 at 2020-02-18 14:34:18
最好都是幂次方的函数,因为自变量都不是角度,三角函数的话不太可信...

排除三角函数的话,下面是在限定条件下搜索到的一个结果。
CODE:
y = (a1 + a2 * x1 + a3 * x1^2 + a4 * x4^2) / (1 + a5 * sqrt(x1) + a6 * sqrt(x2) + a7 * x3^2 + a8 * x1 * x2 * x3 * x4 * x5) + (a9 + a10 * x3 + a11 * x3^2 + a12 * sqrt(x3)) / (1 + a13 * sqrt(x1) + a14 * x2 + a15 * sqrt(x4) + a16 * x1 * x2 * x3 * x4 * x5)

下面是变量的值
a1 = 41.4358532428512
a2 = -9.3450084621907
a3 = 0.503316542998572
a4 = -118.134503443035
a5 = -0.283266564550067
a6 = -0.35919173090838
a7 = -0.00577064680795153
a8 = 0.417014822906771
a9 = -116.495224207427
a10 = -74.2926531230855
a11 = 3.58481686172072
a12 = 177.341976667852
a13 = -0.216714046182733
a14 = -0.299062603053879
a15 = -0.93176728225343
a16 = 0.823131693849135
求助:多变量曲线拟合
result.PNG

我还是那个我,过去如此,未来亦如此!
9楼2020-02-18 19:25:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

代码小站

铜虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
两个问题想问下,1 有经验公式的表达形式吗?2 用于拟合的数据集只有列出来的这么多吗?

发自小木虫Android客户端
qq:3200130194,有问题留言
2楼2020-02-12 18:45:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

figo2015

木虫 (正式写手)

引用回帖:
2楼: Originally posted by 代码小站 at 2020-02-12 18:45:57
两个问题想问下,1 有经验公式的表达形式吗?2 用于拟合的数据集只有列出来的这么多吗?

没有经验公式,只有这么多数据

发自小木虫IOS客户端
3楼2020-02-13 00:28:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangxd0228

木虫 (著名写手)

★ ★
独孤神宇: 金币+2, 鼓励交流 2020-02-15 22:11:03
直接拟合线性模型吧,大体看了一下,有两个因子是不显著的。剩下的三个因子线性显著

发自小木虫Android客户端
呵呵
4楼2020-02-14 09:31:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shikang999

新虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
独孤神宇: 金币+5, 鼓励交流 2020-02-15 22:10:48
1. 一个结果
找到的一个公式,不知是否满足你的需求:

y = a6 + a7 * Sin((a5 + a2 / sqrt(x1) + a3 * sqrt(x4) + a4 * log(x5)) / (a1 - x1 * sqrt(x1))^2) + (a8 * x2 + a9 /x2^2) / (1  + a10 * Acosh(x2))


误差平方和 = 1085.52481200043
R2相关系数 = 0.95967005360806
a1 = 19.8131655361676
a2 = 625.512731990641
a3 = -2478.95724020620
a4 = -182.208407316690
a5 = 292.495799014147
a6 = 79.8975187939561
a7 = 43.1450726807051
a8 = 0.07116260928769
a9 = 0.96038348046488
a10= -0.77644590040080

2. 至于教程,等我自己写的免费查找公式的软件稳定点过后,会专门写一个搜索公式的教程

PS:小木虫发图片按钮找半天,居然没有找到!

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

我还是那个我,过去如此,未来亦如此!
5楼2020-02-15 21:38:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

【答案】应助回帖

数据基本符合下列数学模型
y=6316.3138 + 6597.0761 x1 + 28281.2976 x2 + 9200.4959 x3 + -121.8516 x4 + 660.554 x5 + -7286.06 x1^0.975 -30123.2864 x2^0.975 -10078.5308 x3^0.975 -9667.6374 x4^0.15 -10611.4405 x5^0.15 + 15865.22 (x1*x2*x3*x4*x5)^0.075
决定系数: R^2= 0.983691506031832530
求助:多变量曲线拟合-1
123.jpg

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

QQ:2120156492
6楼2020-02-16 17:27:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

figo2015

木虫 (正式写手)

送红花一朵
引用回帖:
6楼: Originally posted by hzlhm at 2020-02-16 17:27:01
数据基本符合下列数学模型
y=6316.3138 + 6597.0761 x1 + 28281.2976 x2 + 9200.4959 x3 + -121.8516 x4 + 660.554 x5 + -7286.06 x1^0.975 -30123.2864 x2^0.975 -10078.5308 x3^0.975 -9667.6374 x4^0.15 -10611 ...

这个是怎么算出来的?
7楼2020-02-18 14:32:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

figo2015

木虫 (正式写手)

送红花一朵
引用回帖:
5楼: Originally posted by shikang999 at 2020-02-15 21:38:58
1. 一个结果
找到的一个公式,不知是否满足你的需求:

y = a6 + a7 * Sin((a5 + a2 / sqrt(x1) + a3 * sqrt(x4) + a4 * log(x5)) / (a1 - x1 * sqrt(x1))^2) + (a8 * x2 + a9 /x2^2) / (1  + a10 * Acosh(x2))
...

最好都是幂次方的函数,因为自变量都不是角度,三角函数的话不太可信
8楼2020-02-18 14:34:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

引用回帖:
7楼: Originally posted by figo2015 at 2020-02-18 14:32:06
这个是怎么算出来的?...

用matlab软件计算的

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

QQ:2120156492
10楼2020-02-18 20:52:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 figo2015 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 换导师 +5 是柠檬呀! 2024-05-18 10/500 2024-05-18 11:16 by 是柠檬呀!
[教师之家] 来用亲身体会,一起说说年轻老师的辛苦 +17 zylfront 2024-05-17 21/1050 2024-05-18 10:25 by songwz
[有机交流] 求助在硝基苯上进行F-C反应甲基化的方案 105+3 epj006144 2024-05-14 8/400 2024-05-18 10:17 by zyp0009928
[论文投稿] 推荐转投( transfer pending)是否有用? 50+3 lily5289 2024-05-17 3/150 2024-05-18 10:14 by bobvan
[教师之家] 白天不懂夜的黑,90后青椒与60后老板 +3 zylfront 2024-05-18 4/200 2024-05-18 10:12 by zylfront
[教师之家] 博士去高校就是为了有寒暑假吗? +16 wenwen0825 2024-05-16 22/1100 2024-05-18 09:57 by 鱼翔浅底1
[硕博家园] 又见奇葩学生,可咋办 +14 appleapple2 2024-05-17 18/900 2024-05-18 09:14 by 斯坦福01
[硕博家园] 耐高温垫片求购 +6 Sexyflea 2024-05-16 9/450 2024-05-18 08:31 by nanbeiY
[找工作] 浙江的高校现在门槛都这么高吗 +7 dadqweq_qw 2024-05-16 7/350 2024-05-18 07:41 by kangshisan
[硕博家园] 民办高校入职后稳定吗? +8 905452934 2024-05-16 13/650 2024-05-18 06:54 by dong5391
[基金申请] 有知道工材e01情况的吗? +9 xiaopang8958 2024-05-17 12/600 2024-05-17 19:27 by 嗨皮一下
[基金申请] 南林海归博士首聘期 没拿到国自然,被降级后自杀 +32 babu2015 2024-05-13 34/1700 2024-05-17 18:52 by 清晨和风一起
[教师之家] 南京林业大学老师首聘期考核不过,服药自杀 +30 babu2015 2024-05-12 35/1750 2024-05-17 10:55 by hunt_a_job
[教师之家] 青椒的未来在哪里? +14 sxhuabin 2024-05-15 14/700 2024-05-17 08:58 by Quakerbird
[考博] 寻找高分子领域博导 +8 小政雅图 2024-05-14 11/550 2024-05-16 22:10 by 单调的胖子
[考博] 24应届生寻博导 +3 8139173 2024-05-15 4/200 2024-05-16 21:09 by 不容易456
[考博] 深度学习申博 +4 酸菜鱼. 2024-05-15 4/200 2024-05-16 11:12 by 安塔瓦拉多
[教师之家] 学生家长私下联系老师修改成绩不成,唆使19名学生联名要求复核成绩 +23 sjtu2012 2024-05-11 26/1300 2024-05-16 07:57 by zhangysbad
[基金申请] 青年拔尖出炉 +5 babu2015 2024-05-14 6/300 2024-05-15 21:27 by msjy
[考博] 24/25申博,科研助理找导师 +9 cc毕业了吗 2024-05-12 9/450 2024-05-13 14:31 by zxl_1105
信息提示
请填处理意见