MATLAB回归分析求系数问题！！！！

求教大神：y = a0+{a1+a2*(1-m)/m+a3*n/m^2}*(x1+x2)/2+{a4+a5*(1-m)/m+a6*n/m^2}*（x1-x2）/2
怎么用MATLAB求回归系数？？？？
x1 = [281.1 283.0 284.8 286,7 288.7 290.7 292.7]
x2 = [282.5 285.4 288.3 291.3 294.3 297.4 300.6]
y = [280 285 290 295 300 305]
m = 0.972 n = -0.013
求a0-a6的值！

f=@(a,x)a(1)+(a(2)+a(3)*(1-m)/m+a(4)*n/m^2).*(x(1,:)+x(2,:))/2+(a(5)+a(6)*(1-m)/m+a(7)*n/m^2).*(x(1,:)-x(2,:))/2;
x = [281.1 283.0 284.8 286.7288.7 290.7 292.7; 282.5 285.4 288.3 291.3 294.3 297.4 300.6];
y = [280 285 290 295 300 305310];
m = 0.972;n = -0.013;
[a,e]=lsqcurvefit(f,ones(1,7),x,y);
disp(' a0 a1 a2 a3 a4 a5 a6 ');disp(a)
disp('cha_wu');disp(e)

首先，有个问题：系数a和b有什么区别吗？既然都是系数项。 如果不考虑系数项b，则可以这么解： 整理12组x1, x2, x3为如下格式的数据，其中括号里为数据组数： | 1 x1(1) x2(1) x3(1) | x = | 1 x1(2) x2(2) x3(2) | | ... ... ... ... | | 1 x1(12) x2(12) x3(12) | 整理12组y为如下格式： y = [ y(1); y(2); y(3); ... ; y(12) ]; 则可用如下代码求得回归系数项： aa = inv(x' * x) * x' * y; 其中所得系数项为如下格式： aa = [a;a1;a2;a3];