回归分析(7.19).ppt

2004年全国数模竞赛的B题 “电力市场的 输电阻塞管理” 第一个问题： 某电网有8台发电机组，6条主要线路，表 1和表2中的方案0给出了各机组的当前出力和 各线路上对应的有功潮流值，方案1~32给出了 围绕方案0的一些实验数据，试用这些数据确 定各线路上有功潮流关于各发电机组出力的近 似表达式。 由电网的拓扑结构，线路上的有功潮流由机 组出力决定。又根据功率的叠加原理，各线路 上有功潮流应为各发电机组出力的线性组合， 考虑对所有实验数据采用最小二乘法进行线性 拟合，从而得到各线路有功潮流关于各发电机 组出力的近似表达式。 根据表1和表2围绕方案0的1--32组实验数 据，可以列出关于未知数的32个方程的方程 组，利用SAS或Matlab编程求解方程组，得 还需要根据样本值运用假设检验来判断， 以确定求得的回归方程是否有价值。 （Ⅰ）F 检验法 （Ⅱ）r 检验法 (残差平方和) 3、多元线性回归中的检验与预测 2）预测 （1）点预测 （2）区间预测 （4）“有进有出”的逐步回归分析. （1）从所有可能的因子（变量）组合的回归方程中选择最优者； （2）从包含全部变量的回归方程中逐次剔除不显著因子； （3）从一个变量开始，把变量逐个引入方程； 选择“最优”的回归方程有以下几种方法： “最优”的回归方程就是包含所有对Y有影响的变量, 而不包含对Y影响不显著的变量回归方程. 以第四种方法，即逐步回归分析法在筛选变量方面较为理想. 4、逐步回归分析 这个过程反复进行，直至既无不显著的变量从回归方程中剔除，又无显著变量可引入回归方程时为止. 逐步回归分析法的思想： 从一个自变量开始，视自变量Y对作用的显著程度，从大到小地依次逐个引入回归方程. 当引入的自变量由于后面变量的引入而变得不显著时，要将其剔除掉. 引入一个自变量或从回归方程中剔除一个自变量，为逐步回归的一步. 对于每一步都要进行Y值检验，以确保每次引入新的显著性变量前回归方程中只包含对Y作用显著的变量. 统计工具箱中的回归分析命令 1．多元线性回归 2．多项式回归 3．非线性回归 4．逐步回归 多元线性回归 b=regress( Y, X ) 1．确定回归系数的点估计值： 3．画出残差及其置信区间： rcoplot（r，rint） 2．求回归系数的点估计和区间估计、并检验回归模型： [b, bint,r,rint,stats]=regress(Y,X,alpha) 回归系数的区间估计 残差 用于检验回归模型的统计量， 有三个数值：相关系数r 2、 F值、与F 对应的概率p 置信区间 显著性水平 （缺省时为0.05） 例1 解： 1．输入数据： x=[143 145 146 147 149 150 153 154 155 156 157 158 159 160 162 164]; X=[ones(16,1) x]; Y=[88 85 88 91 92 93 93 95 96 98 97 96 98 99 100 102]; 2．回归分析及检验： [b,bint,r,rint,stats]=regress(Y,X) b,bint,stats To MATLAB(liti11) 题目 3．残差分析，作残差图： rcoplot(r,rint) 从残差图可以看出，除第二个数据外，其余数据的残差离零点均较近，且残差的置信区间均包含零点，这说明回归模型 y=-16.073+0.7194x能较好的符合原始数据，而第二个数据可视为异常点. 4．预测及作图： z=b(1)+b(2)* plot(x,Y,k+,x,z,r) To MATLAB(liti12) 多 项 式 回 归 （一）一元多项式回归 （1）确定多项式系数的命令：[p，S]=polyfit（x，y，m） （2）一元多项式回归命令：polytool（x，y，m） 1．回归： y=a1xm+a2xm-1+…+amx+am+1 2．预测和预测误差估计： （1）Y=polyval（p，x）求polyfit所得的回归多项式 在x处 的预测值Y； （2）[Y，DELTA]=polyconf（p，x，S，alpha） 求polyfit所得的回归多项式在x处的预测值Y及预测值的显 著性为1-alpha的置信区间Y DELTA；alpha缺省时为0.5. ± 法一 直接作二次多项式回归： t=1/30:1/30:14/30; s=[11.86 15.67 20.60