MATLAB一元线性回归方程的计算和检验.doc

input 语句键盘输入一组数据( x i,y i), i=1,2, …n。 2. 计算一元线性回归方程 y=ax+b 的系数 a和b ,用两种方法计算: 一是公式: xaybxx yyxxa i ii????????,)( ) )(( 2 ; 二是用最小二乘法的公式求出最小值点( a,b ),使???? 2)( min },(baxybaQ ii 3. 检验回归方程是否有效(用 F 分布检验)。 4. 把散列点( x i,y i )和回归曲线 y=ax+b 画在一个图上。 5. 每种计算法都要有计算框图,且每种计算法都要编成一个自定义函数。 function yiyuan clc; disp(' 从键盘输入一组数据:'); x=input('please Input data x: '); y=input('please Input data y: '); disp(' 一元线性回归的计算和检验:'); disp('1. 公式法'); disp('2. 最小二乘'); disp('3. 检验'); disp('0. 退出'); global a0 b0; while 3 num=input(' 选择求解的方法:'); switch num case 1 [a0,b0]=huigui(x,y) case 2 [a0,b0]=zxec(x,y) case 3 break; case 0 return; otherwise disp(' 输入错误,请重先输入! '); end end X=x';Y=y'; X=[ones(size(X)),X];alpha=; [b,bint,e,rint,stats]=regress(Y,X) if stats(3)<alpha disp(' 有效的 x') end n=[min(x)::max(x)]; f=a0*n+b0; xlabel('x','b');ylabel('y','r');legend(' 散点','k'); end %................................. function [a0,b0]=huigui(x,y) n=length(x); x1=0;y1=0; for i=1:n x1=x1+x(i); y1=y1+y(i); end x0=x1/n; y