牙膏的销量matlab操作命令

?%价格差 x1=[-0.05 0.25 0.60 0 0.25 0.20 0.15 0.05 -0.15 0.15 0.20 0.10 0.40 0.45 0.35 0.30 0.50 0.50 0.40 -0.05 -0.05 -0.10 0.20 0.10 0.50 0.60 -0.05 0 0.05 0.55];
%广告费 x2=[5.50 6.75 7.25 5.50 7.00 6.50 6.75 5.25 5.25 6.00 6.50 6.25 7.00 6.90 6.80 6.80 7.10 7.00 6.80 6.50 6.25 6.00 6.50 7.00 6.80 6.80 6.50 5.75 5.80 6.80];
%销售量 y=[7.38 8.51 9.52 7.50 9.33 8.28 8.75 7.87 7.10 8.00 7.89 8.15 9.10 8.86 8.90 8.87 9.26 9.00 8.75 7.95 7.65 7.27 8.00 8.50 8.75 9.21 8.27 7.67 7.93 9.26];



%图一和图二,需要分别执行
plot(x1,y,'r*');
plot(x2,y,'r*');

x3=x2.*x2; %这里有.,说明是数组之间的乘法
e=ones(30,1); %这是一个30行1列的全为1的向量
x=[e x1' x2' x3'];
alpha=0.05;
[b,bint,r,rint,stats]=regress(y',x,alpha)



%改进,分开执行
x4=x1.*x2;
x5=[e x1' x2' x3' x4'];
[b,bint,r,rint,stats]=regress(y',x5,alpha)


%图四
y2=29.1133+11.1342*x1-7.6080*6.5+0.6712*6.5^2-1.4777*6.5*x1;
plot(x1,y2)


%最后一个图

x6=[x1' x2'];
rstool(x6,y,'quadratic')