| 查看: 1763 | 回复: 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] 最好能附上程序。十分感谢~ |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 数据拟合的问题
已经有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人回复
匿名
用户注销 (著名写手)
- 应助: 84 (初中生)
- 金币: 5015.3
- 散金: 1030
- 红花: 35
- 帖子: 2148
- 在线: 235.4小时
- 虫号: 0
- 注册: 2011-10-29
- 性别: GG
- 专业: 伊斯兰教
6楼2013-11-29 14:29:58
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3881
- 散金: 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 (院士)
- 金币: 146664
- 红花: 1374
- 帖子: 93091
- 在线: 7695小时
- 虫号: 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
- 金币: 3881
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
4楼2013-11-28 11:03:13
25