| 查看: 2262 | 回复: 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 |
回复此楼» 猜你喜欢
焦虑
已经有8人回复
-
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
- 金币: 100702
- 散金: 11601
- 红花: 134
- 沙发: 22
- 帖子: 46727
- 在线: 4404.9小时
- 虫号: 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
