y=[-6.33E-06
9.95E-07
1.85E-05
1.09E-05];
x1=[-2.28E-05
-6.33E-06
9.95E-07
1.85E-05];
x2=[1.33E-05
6.84E-06
7.86E-06
2.65E-05];
% y=a+d1+d2+α1*x1+α2*x2+e
%根据你的公式 a+d1+d2+e 只能被看成是一个数
%公式变为 y=c1*x1+c2*x2+c3
c=regress(y,[x1 x2 ones(size(x1))]);
disp(' c1 c2 c3');
disp(c');
结果
c1 c2 c3
0.6910 -0.6789 0.0000
也就是回归结果是
y=0.6910*x1-0.6789*x2