求一元线形回归方程与预报.docVIP

求一元线形回归方程及预报 一 功能 x为自变量，y为随机变量，给出一组n次观测值（Xi, Yi）, I=1, 2, … n, 求线形回归方程 Y=A+BX 来描述Y与X的变化规律。并用T检验法检验线形回归是否显著。若显著，可用求得的线形回归方程作预报，并给出预报值的置信区间。 二 算法简介 [16] , [15] 设X为可控制变量，Y是依赖于X的随机变量，取X不全相同的n个值X1, X2, …Xn 作n次独立实验，得到容量为n的样本值（X1, Y1）, (X1, Y2) , … , (Xn, Yn) 定Y的数学期望值存在（是X的函数） ， 记为 （x）, 称 (x) 为Y对X的。回归。现假定样本有线形模型 （即（x）为线形函数） Y1=A+BX+1, I=1,2,…,n (1) 其中A, B为待定函数；X为实验误差，且12….n. 独立同步分布：c1~N(0, 62); 对固有的X ,有Y~N (A+BX,)。 这相当于假设 (2) 其中A ，B ， 都不依赖于X 。式（2）称为一元线形回归模型。 在一元线形回归问题中，我们主要解决三个问题 如何用n对数据（Xi,,Yi）对A ,B , 作出点估计； 回归系数B进行假设检验； 给出X0 ，对Y作预报，并给出预报值Y0的置信区间。 先对A、B 作点估计。设估计值为， ，从而得经验回归直线方程 （3） 作离差平方和 （4） 用最小二乘法求A , B ,令 （5） 对该方程组变形，得 （6） 该方程组称为正规（则）方程组，解此方程组，得 . （7） (8) 由此即得经验回归直线方程（3） 实际计算时，可令 （9） 则得A , B的计算公式为 （10） （11） 还可得Q的计算公式[15] (12) 可得的估计。由于 于是由此得的无偏估计 （13） （14） 对回归系数B作T检验。可知 （15） 对回归系数作 假设 给定显著性水平，一次抽样后由（15）式计算得T , 若 则拒绝原假设，即认为线性回归显著；否则认为线形回归不显著。当线性回归显著时可作预报。 对Y作预报。设给出, 则Y的预报值为 = 记 （）=ta/2(n-2). （2） 则Y0的置信区间（置信度为1-a）为 （17） 三`程序使用说明 1`输入参数 N 实数`独立实验次数 R 实数，给定的显著性水平，取R=0.05等值 即ta/2(n-2),由查t分布表得到 X(N) 一维实数组 Y(N) 一维实数组，按(X1,Y1),(X2,Y2),…(Xn,Yn)存放X(N)与Y(N),也按这个顺序输入（出） 2`输出参数 上述1中的输入参数全部输出 经验直线回归方程 G 实数 G= T 实数，由（15）计算的T值 对回归系数作T检验的结果 对X0作预报的结果。 WW=0 表示预报的结束 10 ‘************************************************************************** 20‘ 906 求一元线性回归方程及