妖丶刺尾蝶
2012-11-23, 20:22
>> clear
x1=[1 2 3 4 5 6 7 8 9 10 ];
x2=[0.75 1.25 1.75 2.25 2.75 0.5 1 1.5 2 2.5];
x3=[4.5 5.5 4 5 6 4 5 6 4.5 5.5 ];
x4=[28 30 27 30 27 29 26 29 26 28 ];
x5=[140 120 100 80 60 140 120 100 80 60];
X=[ones(10,1)x1' x2' x3' x4' x5' (x1.^2)' (x2.^2)' (x3.^2)' (x4.^2)' (x5.^2)'];
Y=[3.167367 4.693917 5.0277 11.1097 11.29748 3.038083 3.275633 5.6382 7.9124 10.63053]
[b,bint,r,rint,stats]=regress(Y,X)
b,bint,stats
其实是想解如下方程的参数y=a+bx1+cx2+dx3+ex4+fx5+g(x1^2)+h(x1^2)+i(x1^2)+j(x1^2)+k(x1^2)+l(x1*x2+)+m(x2*x5)+n(x1*x5)
急求~
x1=[1 2 3 4 5 6 7 8 9 10 ];
x2=[0.75 1.25 1.75 2.25 2.75 0.5 1 1.5 2 2.5];
x3=[4.5 5.5 4 5 6 4 5 6 4.5 5.5 ];
x4=[28 30 27 30 27 29 26 29 26 28 ];
x5=[140 120 100 80 60 140 120 100 80 60];
X=[ones(10,1)x1' x2' x3' x4' x5' (x1.^2)' (x2.^2)' (x3.^2)' (x4.^2)' (x5.^2)'];
Y=[3.167367 4.693917 5.0277 11.1097 11.29748 3.038083 3.275633 5.6382 7.9124 10.63053]
[b,bint,r,rint,stats]=regress(Y,X)
b,bint,stats
其实是想解如下方程的参数y=a+bx1+cx2+dx3+ex4+fx5+g(x1^2)+h(x1^2)+i(x1^2)+j(x1^2)+k(x1^2)+l(x1*x2+)+m(x2*x5)+n(x1*x5)
急求~