一元线性回归的matlab实现(含R方和F分布检验).doc

一元线性回归的 matlab 实现(含置信度检验) 说明:正文中命令部分可以直接在 M atlab 中运行,作者(Yangfd09_LZU) 在 MATLAB R2009a() 中运行通过% 求一元线性回归方程% 数据要求:两行。第一行存放 x 的观察值,第二行存放 y 的观察值% 数据文件名: ;变量名: test load N=length(test(1,:)); % 注:也可以用[M,N]=size(test) % 但不能用 N=size(test(1,:)) sx=0;sx2=0;sy=0;sy2=0;sxy=0; for i=1:N sx=sx+test(1,i); sx2=sx2+test(1,i)^2; sy=sy+test(2,i); sy2=sy2+test(2,i)^2; sxy=sxy+test(1,i)*test(2,i); end r=[N,sx;sx,sx2]\[sy;sxy]; a=r(1);b=r(2); %F 分布检验 U=b*sxy; Q=sy2-U; F=(N-2)*U/Q; % 命令窗口中限制回归方程 str=[blanks(5), 'y=' , '(' ,num2str(a), ')' , '+' , '(' ,num2str(b), ')' , '*x' ]; disp( '') disp( ' 回归方程为:') disp(str) disp( 'F- 分布检验: ') stri=[ 'F 计算值: ' ,num2str(F),blanks(4), ' 自由度: f1=1,f2=' ,num2str(N-2)]; dis