Matlab一元一次函数拟合及显著性检验

X=xlsread('E:\DATA\data1.xlsx','sheet1','B2:B13');%读取X数据
Y=xlsread('E:\DATA\data1.xlsx','sheet1','C2:C13');%读取Y数据

%函数拟合如下
Xsum=0;
Ysum=0;
mol=0;
den=0;
for i=1:length(X)
    Xsum=Xsum+X(i);
    Ysum=Ysum+Y(i);
end
Xmean=Xsum/12;
Ymean=Ysum/12;
for i=1:length(X)
   mol=mol+(X(i)-Xmean)*(Y(i)-Ymean);
   den=den+(X(i)-Xmean)^2;
end
b=mol/den;
a=Ymean-b*Xmean;%求出a，b值，即可得到拟合函数方程y=a+bx;

%显著性F检验如下
Yfore=[];
for i=1:length(X)
   Yfore(i)=-0.65*X(i)+3.85;%-0.65为上面求出的b值，3.85为上面求出的a值
end
Q=0;
U=0;
for i=1:length(X)
    Q=Q+(Y(i)-Yfore(i))^2;
    U=U+(Yfore-Ymean)^2;
end
F=U/(Q/(12-2));%查F表进行显著性检验

 
 
想了解更多，请关注微信公众号“GAO戏精光芒万丈”：