| 查看: 2603 | 回复: 6 | |||
[交流]
【求助】回归拟合多元二次方程 已有4人参与
|
|
我的实验是球面对称设计后，对数据进行回归拟合，数据如下： x1 x2 x3 y1 y2 1 18.45 0.17 1.85 109.6 18 2 41 0.17 1.85 128.4 38.04 3 18.45 0.48 1.85 98.7 43.05 4 41.55 0.48 1.85 124.6 55.71 5 18.45 0.17 4.15 129.6 22.61 6 41.55 0.17 4.15 128.8 49 7 18.45 0.48 4.15 103.2 60.19 8 41.55 0.48 4.15 141.7 78.71 9 10 0.33 3 105.8 57.66 10 50 0.33 3 196.1 64.53 11 30 0.05 3 102.3 30.5 12 30 0.6 3 119.5 50.23 13 30 0.33 1 118 65.02 14 30 0.33 5 121.3 82.77 15 30 0.33 3 130 87 目标是拟合成三元二次方程，具体形式如下： y1=b0+b1x1+b2x2+b3x3+b11x1^2+b12x1x2+b13x1x3+b22x2^2+b23x2x3+b33x3^2 y2=............. 通过拟合得到b0,b1,b2,b3,b11等常数的值 要求y1>100 y2>80%,解不等式求出x1,x2,x3, 并找出最适合的解。 手里面有SPSS，lingo和statistica，哪种软件更合适些？具体怎么操作？ 最好能给个模板，谢谢各位了， 我的邮箱liuyu0710@163.com |
回复此楼» 猜你喜欢
时间戳变了，能看出什么问题？
已经有11人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有109人回复
-
重要消息，中午系统在维护
已经有10人回复
-
两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练
已经有7人回复
-
明天放榜？
已经有4人回复
-
投稿咨询
已经有6人回复
-
时间戳又变了8-15
已经有26人回复
-
今天维护系统维护 祝所有人 高中
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有16人回复
-
2027广东省杰青
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Matlab多元回归！
已经有10人回复
- matlab 多元函数拟合
已经有3人回复
- 求助——求解三元二次多项式方程组
已经有9人回复
- 【求助】非线性方程组的求解问题
已经有6人回复
- 求助！！求解三元二次方程
已经有5人回复
- 麻烦高手帮忙解一下三元二次方程组！谢谢谢谢！
已经有4人回复
- 响应面中中心点的个数如何确定
已经有11人回复
- 回归和拟合有啥区别？
已经有13人回复
- MATLAB用于三元回归参数
已经有3人回复
- (1+x)^1/2次方能否二次项展开，怎么做呢？
已经有4人回复
- 【紧急求助】1的无穷大次方、无穷大的0次方什么意义？为什么是未定式？
已经有7人回复
- 求助SAS 多元非线性回归拟合编程
已经有3人回复
- 响应曲面法和多元回归法拟合的模型比较
已经有3人回复
- 【求助】三元二次方程求最大值
已经有3人回复
- 【求助】三元二次方程的回归拟合
已经有5人回复
- 【求助】请高人帮忙三元二次方程
已经有11人回复
dzl723
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1722.2
- 散金: 100
- 红花: 2
- 帖子: 473
- 在线: 111.2小时
- 虫号: 1062905
- 注册: 2010-07-23
- 专业: 中药制剂
2楼2010-07-31 01:03:56
mrzouhao
木虫之王 (文坛精英)
- PhEPI: 4
- 应助: 179 (高中生)
- 贵宾: 0.108
- 金币: 102385
- 散金: 11621
- 红花: 134
- 沙发: 23
- 帖子: 47341
- 在线: 4437.5小时
- 虫号: 201350
- 注册: 2006-02-28
- 专业: 药剂学
3楼2010-07-31 07:46:57
qiruiyun
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4706.3
- 散金: 66
- 红花: 1
- 帖子: 581
- 在线: 281.4小时
- 虫号: 744562
- 注册: 2009-04-09
- 性别: MM
- 专业: 光学和光电子材料
liuyu0710(金币+5): 2010-08-05 19:39:08
|
我用mathematica写了个代码 cc = Import["d:\\xmc.txt", "Table"]; y1 = Table[Delete[cc[], 5], {i, Length[cc]}]; y2 = Table[Delete[cc[], 4], {i, Length[cc]}]; fitY1 = Fit[ y1, {1, x1, x2, x3, x1^2, x1 x2, x1 x3, x2^2, x2 x3, x3^2}, {x1, x2, x3}] fitY2 = Fit[ y2, {1, x1, x2, x3, x1^2, x1 x2, x1 x3, x2^2, x2 x3, x3^2}, {x1, x2, x3}] 结果为： y1=79.7916 - 1.98694 x1 + 0.0439335 x1^2 + 90.7723 x2 + 3.26941 x1 x2 - 295.394 x2^2 + 25.4082 x3 - 0.0779321 x1 x3 + 1.30797 x2 x3 - 3.43553 x3^2 y2=-143.898 + 5.43815 x1 - 0.0813438 x1^2 + 502.35 x2 - 1.10043 x1 x2 - 696.748 x2^2 + 26.1101 x3 + 0.118026 x1 x3 + 16.5905 x2 x3 - 4.94622 x3^2 你参考一下吧。 |
4楼2010-07-31 08:24:23
5楼2010-07-31 08:36:54
6楼2014-07-02 09:43:24
7楼2018-11-27 16:38:37
