| 查看: 1589 | 回复: 6 | |||||
[求助]
跪求数据拟合，形式是y=a*(x1)^b*(x2)^c*(x3)^d
|
|
数据拟合，最终拟合的公式是y=a*(x1)^b*(x2)^c*(x3)^d，已知y, x1,x2,x3,求a，b，c，d 其中 y=[0.145 0.086 0.074 0.181 0.097 0.084 0.219 0.126 0.103 0.24 0.152 0.124 0.205 0.117 0.1] x1=[0.0059 0.0059 0.0059 0.0107 0.0107 0.0107 0.0131 0.0131 0.0131 0.0161 0.0161 0.0161 0.0115 0.0115 0.0115] x2=[304 304 304 304 304 304 304 304 304 258 258 258 340 340 340] x3=[0.850 4.850 8.850 0.850 4.850 8.850 0.850 4.850 8.850 0.850 4.850 8.850 0.850 4.850 8.850] 最好能附上程序。十分感谢~ |
回复此楼» 猜你喜欢
材料学硕297已过四六级求调剂推荐
已经有8人回复
-
本人考085602 化学工程 专硕
已经有21人回复
-
265求调剂
已经有12人回复
-
求调剂
已经有5人回复
-
一志愿中南化学337求调剂
已经有6人回复
-
350求调剂
已经有3人回复
-
295复试调剂
已经有8人回复
-
0817 化学工程 299分求调剂 有科研经历 有二区文章
已经有22人回复
-
287求调剂
已经有7人回复
-
329求调剂
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 数据拟合的问题
已经有3人回复
- 怎样拟合数据求出x=f(A,B,C),y=g(A,B,C,D)中的A,B,C,D参数
已经有18人回复
- 关于origin数据拟合方程后，求一个y值对应的x值
已经有7人回复
- Matlab非线性拟合（y=sqrt(1-(x-a).^2)）
已经有5人回复
- 通过实验数据拟合，求解公式中的参数
已经有14人回复
- matlab数据拟合问题求助
已经有7人回复
- 数据点的二次曲线方程拟合
已经有6人回复
- 非线性最小二乘拟合（初值问题）lsqcurvefit 函数的设定
已经有14人回复
- 一阶微分方程求解，并优化参数值
已经有14人回复
- origin进行非线性拟合如何选择函数或者自定义合适的函数
已经有19人回复
- 求教origin拟合简单的渗流方程（又称逾渗方程），做复合材料的应该会遇到吧
已经有20人回复
- 菜鸟求助：关于Origin多元线性回归拟合y=x1*（a*x2+b*x3+c)
已经有10人回复
- 【求助】多元方程组求解及拟合
已经有9人回复
- 求助： 1stopt 进行非线性拟合
已经有5人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
l^o^ise: 金币+12, ★★★★★最佳答案 2013-12-02 09:36:04
感谢参与，应助指数 +1
l^o^ise: 金币+12, ★★★★★最佳答案 2013-12-02 09:36:04
|
Parameters a, b, c, d; Variable x1, x2, x3, y; Function y=a*(x1)^b*(x2)^c*(x3)^d; data; 0.0059 304 0.85 0.145 0.0059 304 4.85 0.086 0.0059 304 8.85 0.074 0.0107 304 0.85 0.181 0.0107 304 4.85 0.097 0.0107 304 8.85 0.084 0.0131 304 0.85 0.219 0.0131 304 4.85 0.126 0.0131 304 8.85 0.103 0.0161 258 0.85 0.24 0.0161 258 4.85 0.152 0.0161 258 8.85 0.124 0.0115 340 0.85 0.205 0.0115 340 4.85 0.117 0.0115 340 8.85 0.1 迭代数: 38 计算用时(时:分:秒:微秒): 00:00:01:641 优化算法: 最大继承法 计算结束原因: 达到收敛判断标准 均方差(RMSE): 0.00714506098071919 残差平方和(SSE): 0.000765778446272938 相关系数(R): 0.989972462281014 相关系数之平方(R^2): 0.980045476074734 决定系数(DC): 0.980045242663603 卡方系数(Chi-Square): 0.00318900654592896 F统计(F-Statistic): 638.481340706879 参数 最佳估算 ---------- ------------- a 3.49152900492693 b 0.531011058765137 c -0.0940361498458002 d -0.305225716804401 ====== 结果输出 ===== No 实测值y 计算值y 1 0.145 0.1404031 2 0.086 0.0825140 3 0.074 0.0686755 4 0.181 0.1926016 5 0.097 0.1131907 6 0.084 0.0942075 7 0.219 0.2144515 8 0.126 0.1260318 9 0.103 0.1048950 10 0.24 0.2429876 11 0.152 0.1428022 12 0.124 0.1188528 13 0.205 0.1980238 14 0.117 0.1163773 15 0.1 0.0968596 |
2楼2013-11-27 17:50:43
peterflyer 
木虫之王 (文学泰斗)
peterflyer
- 数学EPI: 10
- 应助: 20282 (院士)
- 金币: 146160
- 红花: 1374
- 帖子: 93091
- 在线: 7694.3小时
- 虫号: 1482829
- 注册: 2011-11-08
- 性别: GG
- 专业: 功能陶瓷
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
l^o^ise(feixiaolin代发): 金币+5, 又一种方法。 2013-11-27 23:16:12
l^o^ise: 金币+2, ★★★很有帮助 2013-12-02 09:36:15
感谢参与，应助指数 +1
l^o^ise(feixiaolin代发): 金币+5, 又一种方法。 2013-11-27 23:16:12
l^o^ise: 金币+2, ★★★很有帮助 2013-12-02 09:36:15
|
求解过程： 先对两边求对数，并令 fi(a,b,c,d)=Lnyi-Lna-b*Lnx1i-c*Lnx2i-d*Lnx3i (i=1~N,N为数据个数，此题中N=15） 再令F(a,b,c,d)=Σfi(a,b,c,d)^2 将F分别对Lna、b、c、d求偏导并分别令其为零，得到四元一次方程组，解之，得到Lna、b、c、d的最小二乘法拟合值。最后利用 a=exp[Lna]获得a的拟合值。 完毕。 |
3楼2013-11-27 22:30:09
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
4楼2013-11-28 11:03:13
peterflyer
木虫之王 (文学泰斗)
peterflyer
- 数学EPI: 10
- 应助: 20282 (院士)
- 金币: 146160
- 红花: 1374
- 帖子: 93091
- 在线: 7694.3小时
- 虫号: 1482829
- 注册: 2011-11-08
- 性别: GG
- 专业: 功能陶瓷
5楼2013-11-28 12:02:02
匿名
用户注销 (著名写手)
- 应助: 84 (初中生)
- 金币: 5015.3
- 散金: 1030
- 红花: 35
- 帖子: 2148
- 在线: 235.4小时
- 虫号: 0
- 注册: 2011-10-29
- 性别: GG
- 专业: 伊斯兰教
6楼2013-11-29 14:29:58
baconhu
金虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 1909.3
- 红花: 15
- 帖子: 632
- 在线: 31.5小时
- 虫号: 2800309
- 注册: 2013-11-14
- 专业: 高分子材料的加工与成型
【答案】应助回帖
★
感谢参与，应助指数 +1
l^o^ise: 金币+1, ★★★很有帮助 2013-12-02 09:36:40
感谢参与，应助指数 +1
l^o^ise: 金币+1, ★★★很有帮助 2013-12-02 09:36:40
| 看看这个帖子，里面全是绘图函数和程序，还有实例。http://muchong.com/bbs/viewthread.php?tid=6682326 |
7楼2013-11-29 15:11:58
