工程数学第21讲.pptVIP

工程数学第21讲 第九章 方差分析及回归分析 §1 单因素试验的方差分析 (一)单因素试验 在科学试验和生产实践中, 影响一事物的因素往往是很多的.例如, 在化工生产中,有原料成份, 原料剂量, 催化剂, 反应温度, 压力, 溶液浓度, 反应时间, 机器设备及操作人员的水平等因素. 每一因素的改变都有可能影响产品的数量和质量. 有些因素影响较大, 有些较小. 为了使生产过程得以稳定, 保证优质, 高产, 就有必要找到对产品质量有显著影响的那些因素. 为此, 我们需进行试验, 方差分析就是根据试验结果进行分析, 鉴别各个有关因素的作用的有效方法. 在试验中, 将要考察的指标称为试验指标. 影响试验指标的条件称为因素. 因素可分为两类, 一类是人们可以控制的(可控因素); 一类是人们不能控制的. 例如, 反应温度, 原料剂量, 溶液浓度等是可以控制的, 而测量误差, 气象条件等一般是难以控制的. 以下所说的因素都是可控因素. 因素所处的状态, 称为该因素的水平. 如果一项试验中只有一个因素在改变称为单因素试验, 如果多于一个因素在改变称为多因素试验. 例1 设有三台机器, 生产规格相同的铝合金薄板. 取样,测量薄板的厚度得结果如下表所示: 表9.1 铝合金板的厚度 这里, 试验的指标是薄板的厚度. 机器为因素, 不同的三台机器就是这个因素的三个不同的水平. 假定除机器这一因素外, 材料的规格, 操作人员的水平等其它条件都相同. 这是单因素试验. 试验的目的是为了考察各台机器所生产的薄板的厚度有无显著的差异, 即考察机器这一因素对厚度有无显著的影响. 如果厚度有显著差异, 就表明机器这一因素对厚度的影响是显著的. 例2 下面给出了随机选取的, 用于计算器的四种类型的电路的响应时间(以毫秒计). 表9.2 电路的响应时间 例3 一火箭使用四种燃料, 三种推进器作射程试验. 每种燃料与每种推进器的组合各发射火箭两次, 得射程如下(以海里计): 表9.3 火箭的射程 这里的试验指标是射程. 推进器和燃料是因素, 它们分别有3个, 4个水平. 这是一个双因素试验. 试验的目的在于考察在各种因素的各个水平下射程有无显著的差异, 即考察推进器和燃料这两个因素对射程是否有显著的影响. 本节限于讨论单因素试验. 就例1讨论. 在例1中, 在因素的每一个水平下进行独立试验, 其结果是一个随机变量. 表中数据可看成来自三个不同总体(每个水平对应一个总体)的样本值. 将各个总体的均值依次记为m1,m2,m3. 按题意需检验假设 H0:m1=m2=m3, H1:m1,m2,m3不全相等.进而假设各总体均为正态变量, 且各总体的方差相等, 但参数均未知. 这就是一个检验同方差多个正态总体均值是否相等的问题. 设因素A有s个水平A1,A2,...,As, 在水平Aj(j=1,2,...,s)下, 进行nj次独立试验, 得到如下表的结果. 假定:各个水平Aj(j=1,2,...,s)下的样本X1j,X2j ... 2,...,s)的正态总体N(mi,s2), mi与s2未知. 且设不同水平Aj下的样本之间相互独立. 由于Xij~N(mj,s2), 即有Xij-mj~N(0,s2), 故Xij-mj可看成随机误差. 记Xij-mj=eij, 则 其中mj与s2均为未知参数, (1.1)式称为单因素试验方差分析的数学模型. 这是本节研究的对象.方差分析的任务是对于模型(1.1),1o 检验s个总体N(m1,s2),...,N(ms,s2)的均值是否相等, 即检验假设 H0:m1=m2=...=ms, H1:m1,m2,...,ms不全相等. (1.2)2o 作出未知参数m1,m2,...,ms, s2的估计. 为了将问题(1.2)写成便于讨论的形式, 令n=n1+n2+...+ns, 再令 m称为总平均. 再引入 dj=mj-m, j=1,2,...,s, (1.4) 此时有n1d1+n2d2+...+nsds=0, dj表示水平Aj下的总体平均值与总平均的差异, 习惯上将dj称为水平Aj的效应. 利用这些记号, 模型(1.1)可改写成 而假设(1.2)等价于假设 H0:d1=d2=...=ds=0, H1:d1,d2,...,ds不全为零. (1.2) (二)平方和的分解 引入总偏差平方和 是数据的总平均. ST能反映全部试验数据之间的差异, 因此ST又称为总变差. 又记水平Aj下的样本平均值为`X?j, 即 将ST写成 注意最后一项交叉项等于0. 即交叉项为 于是将ST分解成为 上述SE的各项(Xij-`X?j)2表示在水平Aj下, 样本观察值与样本均值的差异, 这是由随机误差所引起的. S