matlab达人，这里可以打擂了！（多元函数二阶拟合）

a88wangwei

我有如下的一个数表，里面包含x1-x6的自变量的值，以及对应的y的值。如何将其进行二阶拟合呢？？

x1	x2	x3	x4	x5	x6	y
500	22.2	44.4	33.3	33.3	88.9	1.4
555.6	277.8	77.8	44.4	44.4	266.7	10
611.1	88.9	22.2	300	166.7	122.2	9.1
666.7	100	211.1	266.7	22.2	44.4	5.2
722.2	0	144.4	122.2	111.1	11.1	3.4
777.8	111.1	122.2	166.7	0	100	3.6
833.3	300	233.3	244.4	200	0	12.3
888.9	166.7	11.1	22.2	188.9	277.8	11.3
944.4	55.6	266.7	11.1	211.1	200	6.7
1000	288.9	277.8	133.3	122.2	144.4	11.1
1055.6	177.8	288.9	288.9	300	188.9	27
1111.1	222.2	244.4	211.1	66.7	244.4	12.8
1166.7	211.1	66.7	88.9	244.4	66.7	12.1
1222.2	133.3	0	77.8	222.2	111.1	9.8
1277.8	244.4	177.8	0	277.8	233.3	19.9
1333.3	77.8	188.9	233.3	288.9	155.6	12.7
1388.9	66.7	133.3	200	55.6	288.9	5.6
1444.4	200	100	155.6	233.3	177.8	16.3
1500	266.7	222.2	100	133.3	166.7	10.9
1555.6	155.6	155.6	55.6	11.1	33.3	2
1611.1	33.3	300	144.4	88.9	55.6	5.1
1666.7	11.1	55.6	177.8	177.8	133.3	3.9
1722.2	255.6	88.9	111.1	155.6	300	19.7
1777.8	188.9	33.3	66.7	77.8	77.8	4.7
1833.3	122.2	200	188.9	255.6	255.6	16.1
1888.9	144.4	255.6	277.8	100	22.2	7.5
1944.4	233.3	111.1	255.6	266.7	222.2	29
2000	44.4	166.7	222.2	144.4	211.1	6.7

如果各位觉得上面的数据看着不舒服，附件中有excel格式的。
另，请达人们将matlab的实现程序也贴出来，越精炼越好。
现在计时开始。

a88wangwei

半个小时过去:time:。。。。。。

stats01

这个问题有点复杂，可能的结果如下：
Yhat=1.8043+.008347X3-.02035X5-.0001149X2^2+.0002858X2X5+.0001778X2X6+.00012407X4X5
RSS=27.6957, MSe=1.3188
R^2=0.97904

TBE_Legend

题目没半点难度不说。

自己想让大家帮忙，就不要以什么擂台的名义了， 很不厚道， 会让人瞧不起的。

1. -4.63062 + 0.0548905 x[1] - 0.0000528502 x[1]^2 -
   
2. 1.87352*10^-9 x[1]^3 + 1.80421*10^-11 x[1]^4 -
   
3. 4.81012*10^-15 x[1]^5 - 0.169085 x[2] + 0.000498351 x[2]^2 +
   
4. 8.32498*10^-6 x[2]^3 - 3.77583*10^-8 x[2]^4 +
   
5. 3.50405*10^-11 x[2]^5 - 0.182677 x[3] + 0.00122014 x[3]^2 +
   
6. 6.82775*10^-6 x[3]^3 - 6.93016*10^-8 x[3]^4 +
   
7. 1.35503*10^-10 x[3]^5 + 0.0758468 x[4] - 0.00107314 x[4]^2 +
   
8. 6.84147*10^-6 x[4]^3 - 1.98795*10^-8 x[4]^4 +
   
9. 2.49796*10^-11 x[4]^5 - 0.0325794 x[5] + 0.0000833505 x[5]^2 -
   
10. 7.92855*10^-6 x[5]^3 + 7.0895*10^-8 x[5]^4 - 1.4634*10^-10 x[5]^5 -
    
11. 0.236513 x[6] + 0.00417562 x[6]^2 - 0.0000193053 x[6]^3 +
    
12. 1.22548*10^-8 x[6]^4 + 5.3544*10^-11 x[6]^5

复制代码

stats01

典型的过拟合（over-parameterization）。

liuyalong008

非常赞同四楼观点
别把大家当小丑
有问题不丢人，但是装逼就比较可恨了

ps: nlinfit可以解决你的问题
下次要摆擂，弄个有难度一点的，免得贻笑大方

qibbxxt

用

1. regstats

复制代码

函数应该就可以搞定了