正交设计-方差分析.ppt

11.1正交试验方差分析的数学模型(一)数学模型根据一般线性模型的假定,若9次试验结果(如例10.2中的转化率)以x1、x2,x表示,我们首先假定:(1)三个因素间没有交互作用。(2)为9个数据可分解为:x1=+a1+b1+c1+1x2=+a1+b2+c2+2x3=+a1+b3+c3+3x4=+a2+b1+c2+4x5=+a2+b+c3+5,x6=+a2+b3+c1+6x7=+a3+b1+c3+7x8=+a3+b2+c1+8x9=+a+b+c+9其中:一般平均;估计=xi=x1+x2+x9叫全部数据的总体平均值。a1、a2、a3表示A在不同水平时的效应。b1、b2、b3表示B在不同水平时的效应。c1、c2、c3表示C在不同水平时的效应。(3)各因素的效应为零,或者,各因素的效应的加和为零ai=0bi=0ci=0,(4)i是试验误差,它们相互独立,且遵从标准正态分布N(0,1),所以多个试验误差的平均值近似等于零。(二)参数估计有了数学模型,还应通过子样的实测值,对以上的各个参数作出估计。由数理统计知识E（）E（）表示的数学期望。即，是的一个无偏估计量。可表示为：,11.2正交试验的方差分析法,一、方差分析的必要性极差分析不能估计试验中以及试验结果测定中必然存在的误差大小。为了弥补这个缺点，可采用方差分析的方法。方差分析法是将因素水平（或交互作用）的变化所引起的试验结果间的差异与误差波动所引起的试验结果间的差异区分开来的一种数学方法所谓方差分析，就是给出离散度的各种因素将总变差平方和进行分解，而你还进行统计检验的一直数学方法。,二、单因素方差分析法方差分析法的基本思路：（1）由数据中的总变差平方和中分出组内变差平方和、组间变差平方和，并赋予它们的数量表示；（2）用组间变差平方和与组内变差平方和在一定意义下进行比较，如两者相差不大，说明因素水平的变化对指标影响不大；如两者相差较大，组间变差平方和比组内变差平方和大得多，说明因素水平的变化影响很大，不可忽视；（3）选择较好的工艺条件或进一步的试验方向。,例11.1考察温度对一化工产品的得率的影响，选了五种不同的温度，同一温度做了三次试验，结果如下：,11.3.1：正交设计方差分析的步骤11.3.2：3水平正交设计的方差分析11.3.3：混合型正交设计的方差分析11.3.4：拟水平法的方差分析11.3.5：重复试验的方差分析,11.3：正交试验方差分析,计算离差的平方和:设用正交表安排m个因素的试验，试验总次数为n,试验的结果分别为x1,x2,xn.假定每个因素有na个水平，每个水平做a次试验，则n=ana.1)总离差的平方和ST记:记为其中ST反映了试验结果的总差异，它越大，说明各次试验的结果之间的差异越大。试验结果之所以有差异，一是由因素水平的变化所引起的，二是因为有试验误差。,11.3.1正交设计方差分析的步骤,2)各因素离差的平方和由单因素的方差分析记为其中Ki表示因素的第i个水平a次试验结果的和。SA反映了因素A的水平变化时所引起的试验结果的差异，即因素A对试验结果的影响。对于两因素的交互作用，我们把它当作一个新的因素。如果交互作用占两列，则交互作用的离差的平方和等于这两列的离差的平方和之和。比如SAxB=S(AxB)1+S(AxB)2,3)试验误差的离差的平方和SE设S因+交为所有因素以及要考虑的交互作用的离差的平方和，因为ST=S因+交+SE,所以SE=ST-S因+交计算自由度:试验的总自由度f总=试验总次数-1=n-1各因素的自由度f因=因素的水平数-1=na-1两因素交互作用的自由度等于两因素的自由度之积fAxB=fAXfB试验误差的自由度fE=f总-f因+交,计算平均离差平方和(均方):在计算各因素离差平方和时，我们知道，它们都是若干项平方的和,它们的大小与项数有关，因此不能确切反映各因素的情况。为了消除项数的影响，我们计算它们的平均离差的平方和。因素的平均离差平方和=(因素离差的平方和)/因素的自由度=S因/f因试验误差的平均离差平方和=(试验误差的离差的平方和)/试验误差的自由度=SE/fE求F比:将各因素的平均离差的平方和与误差的平均离差平方和相比，得出F值。这个比值的大小反映了各因素对试验结果影响程度的大小。,对因素进行显著性检验:给出检验水平，从F分布表中查出临界值F(f因，fE)。将在“求F比”中算出的F值与该临界值比较，若FF(f因，fE)，说明该因素对试验结果的影响显著，两数差别越大，说明该因素的显著性越大。,11.3.2：3水平正交设计的方差分析,例1(无交互作用):磁鼓电机是彩色录像机磁鼓组件的关键部件之一，按质量要求其输出力矩应大于210g.cm。某生产厂过去这项指标的合格率较低，从而希望通过试验找出好的条件，以提高磁鼓电机的输出力矩。根据工程技术人员的经验，取试验因素和相应水平如下表：,解：(选用正交表L9(34)表头设计：试验计划与试验结果：,详细计算如下：,列方差分析表如下：最佳条件的选择：对显著因子应取最好的水平，对不显著因子的水平可以任意选取；在实际中通常从降低成本操作方便等角度加以选择，上面的例子中对因子A与B应选择A2B2，因子C可以任选，譬如为节约材料可选择C1,验证试验：对A2B2C1进行三次试验，结果为：234，240，220，平均值为231.3.此结果是满意的,（二）有交互作用的正交试验的方差分析当任意两因素之间(如A与B）存在交互作用而且显著时，则不论因素A、B本身的影响是否显著，A和B的最佳因素都应从A与B的搭配中去选择例2某分析试验，起测定值受A、B、C三种因素的影响，每因素去两个水平，由于因素间存在交互作用，在设计试验方案时，可选用L8(27)表，试验安排结果如表（试验指标要求越小越好）,因素,试验号,12345678K1K2QiSi,A1,B1,AB3,C4,BC6,误差7,试验指标（经简化后）,AC5,11112222-506.253.125,11221122+10-1581.2578.125,112222110-56.253.125,12121212-3025706.25703.125,2121212125-30256.25253.15,12211221-506.253.125,122121125-1031.2528.125,05-100-510-1510,正交试验结果计算表,结果表明B、C、AC对试验结果影响最大，B可取B2，而A和C见存在显著的交互作用，可通过二元表和二元图来确定其最优水平,由图可知，A2C1最好,故最佳试验条件为A2B2C1，这正好是第7号试验。事实上，从试验结果看，它的效果也最好。,说明：对二水平因素，平方和的计算有一个简单的公式设计算方法对任何二水平的因素都是适用的，设共做了n次试验，某一列是二水平，相应的K值是K1和K2则该列的平方和S为:,例4（无交互作用）某钢厂生产一种合金，为便于校直冷拉，需要进行一次退火热处理,以降低合金的硬度。根据冷加工变形量，在该合金技术要求范围内,硬度越低越好。试验的目的是寻求降低硬度的退火工艺参数。考察的指标是洛氏硬度(HR)，经分析研究，要考虑的因素有3个：退火温度A，保温时间B，冷却介质C。,11.3.3：混合型正交设计的方差分析,解：,方差分析表：从F值和临界值的比较看出，各因素均无显著影响，相对来说，B的影响大些。为提高分析精度，我们只考虑因素B，把因素A，C都并入误差。这样一来，SE就变成SA+SC+S4+S5=0.445+0.18+1.125+0.500=2.250，再列方差分析表。,方差分析表(2)：临界值F0.05(1，6)=5.99,F0.01(1,6)=13.75从F值和临界值的比较来看，因素B就是显著性因素了。因素影响从大到小的顺序为BCA，选定的最优方案应为A2B2C1,例6:钢片在镀锌前要用酸洗的方法除锈。为了提高除锈效率，缩短酸洗时间，先安排酸洗试验。考察指标是酸洗时间。在除锈效果达到要求的情况下，酸洗时间越短越好。要考虑的因素及其水平如表：选取正交表L9(34)，将因素C虚拟1个水平。据经验知，海鸥牌比OP牌的效果好，故虚拟第2水平并安排在第1列。,11.3.4拟水平法的方差分析,解：虚拟水平的因素C的第1水平重复3次，第二水平重复6次。因此，离差平方和为：其余因素的离差平方和为误差的离差平方和为：,方差分析表：从F值和临界值比较看出，各因素均无显著影响，相对来说，因素D