![稳健参数设计](https://wfqqreader-1252317822.image.myqcloud.com/cover/826/43806826/b_43806826.jpg)
2.2.2 案例2-1
有四种不同品牌的汽车轮胎:A1—韩国制造;A2—本厂制造;A3—国内甲厂制造;A4—国内乙厂制造。
为了比较四种轮胎的耐磨损性,每种轮胎各取6个进行强化磨损实验,测量磨损量数据,见表2-2:
表2-2 磨损量数据表 (单位:mm)
![](https://epubservercos.yuewen.com/9A77D1/23083821501912306/epubprivate/OEBPS/Images/35_01.jpg?sign=1739401307-8yIEmmNv6e7FlkjqJEQAubnZVIuVfny0-0-3f17eddf58da39ab6a8b3b6d5261a06f)
试分析四种轮胎的耐磨损性有无显著性差异。
1. 方差分析表
应用Minitab中的方差分析模块,可以输出方差分析表,见表2-3。
表2-3 磨损量的方差分析表
![](https://epubservercos.yuewen.com/9A77D1/23083821501912306/epubprivate/OEBPS/Images/35_02.jpg?sign=1739401307-JaSuNAjp4NCTzTVjYSDDqPeNAtoQoLQX-0-dea41c3f46086aaa1a66204b9d5649ff)
观察方差分析表,判定因子A对输出特性Y的影响是否显著有两种方法:
(1)F检验法。在一定的显著水平下(通常取5%或1%),从F分布表查取临界值Fα,如果方差分析表中的F值(本例F值等于24.27)大于临界值Fα,则判定因子A对输出特性Y有显著性影响,否则判定因子A对输出特性Y无显著性影响。
(2)P值法。如果P<1%,则认为因子A对输出特性Y有高度显著性影响;如果1%<P<5%,则认为因子A对输出特性Y有显著性影响;如果P>5%,则认为因子A对输出特性Y无显著性影响。
上述两种方法是等价的,两者的判定结果是一致的。
2. 置信区间图
应用Minitab中的方差分析模块,可以输出磨损量的置信区间图,见图2-1。
![](https://epubservercos.yuewen.com/9A77D1/23083821501912306/epubprivate/OEBPS/Images/36_01.jpg?sign=1739401307-poSlvxCuvuqR8kmdjm8vXJeiSrTVkbII-0-79badc95bbacdb69163a94a030e9cc2a)
图2-1 磨损量的置信区间图
从上面的置信区间图可以看出,如果置信区间不相交,则说明有显著性差异,否则说明无显著性差异。本例的分析结论是:本厂制造的轮胎与韩国制造的轮胎的耐磨性无显著性差异,但是,明显优于国内甲、乙两个工厂制造的轮胎。
3. 因子主效应图
应用Minitab中的方差分析模块,可以输出磨损量的主效应图,见图2-2。
![](https://epubservercos.yuewen.com/9A77D1/23083821501912306/epubprivate/OEBPS/Images/36_02.jpg?sign=1739401307-MBh7ZGYXHAwAt1fQKfW26wSLvY4uLMv0-0-203fefa153b7c7054d18d93cbb7416a7)
图2-2 磨损量的主效应图
在图2-2中,横坐标表示因子A的水平,纵坐标表示与因子A的水平对应的输出特性Y的平均值。
从图2-2中可以看出,A2(本厂制造)的轮胎平均磨损量与A1(韩国制造)的轮胎平均磨损量较接近,但是,明显优于国内甲、乙两厂的轮胎平均磨损量。
4. 残差分析图
应用Minitab中的方差分析模块,可以输出磨损量的残差分析图,见图2-3。
![](https://epubservercos.yuewen.com/9A77D1/23083821501912306/epubprivate/OEBPS/Images/37_01.jpg?sign=1739401307-bh0Kk94CtJNVBkMwbJcb98UDuvQxERsF-0-4ee3515217bb7eca3860feab93060ae3)
图2-3 磨损量的残差分析图
所谓残差是指Y的观察值与行平均值之差,即
![](https://epubservercos.yuewen.com/9A77D1/23083821501912306/epubprivate/OEBPS/Images/37_02.jpg?sign=1739401307-NEEvBj8CzNClLJ7VC2ADGekD1pRdRyM5-0-bae65fe52ce9c854a7d23e77d1738efd)
残差分析是为了检验残差eij是否独立同分布于正态分布N(0,σ2)。实际上,在进行方差分析时,残差应当满足以下假定条件:
1)残差服从正态分布,这一点可以从图2-3中左侧的正态概率图和直方图得到验证;
2)残差平均值为0,这一点可以从图2-3中右上角的图得到验证(横坐标轴上下方的点数差不多);
3)残差数据相互独立,这一点可以从图2-3中右下角的图得到验证(数据点随机排列,没有缺陷);
4)残差方差相等,这一点可以通过对残差进行等方差检验来验证。
应用Minitab中的方差分析模块,可以进行等方差检验,并可输出图2-4。
![](https://epubservercos.yuewen.com/9A77D1/23083821501912306/epubprivate/OEBPS/Images/37_03.jpg?sign=1739401307-xb8uFKIBbFI0GKfnzdR2DT4w16Xx0y25-0-f32f66b7fa25f5ea80a788dc8a57d3a0)
图2-4 磨损量的等方差检验
从图2-4可以看出,本例满足等方差的假定。
5. 贡献率分析
上述传统的方差分析方法,要求在因子的不同水平下残差服从随机、相互独立、平均值为0、等方差的正态分布,在满足这些假定条件下才能进行F检验和计算P值。这种方法虽然理论上十分严谨,但是在实际应用时,这些假定条件未必能全部满足。田口博士提出,只要不进行F检验和计算P值,上述假定条件就没有必要。其实,只要利用波动平方和的分解公式,计算纯波动与贡献率,就完全可以判定因子A对输出特性的影响大小。
从表2-4可以看出,轮胎品牌对磨损量的贡献率为75.22%,除轮胎品牌以外,误差对磨损量的贡献率为24.78%。
表2-4 磨损量的贡献率分析表
![](https://epubservercos.yuewen.com/9A77D1/23083821501912306/epubprivate/OEBPS/Images/38_01.jpg?sign=1739401307-sZQZzhyezKDor1HNyKL5d2gyozpRaqCm-0-4036c5a305a35478266aa1c620261cc2)
注:1)因子A的纯波动445.69,等于因子A的波动464.83减去3(A的自由度)倍误差方差6.38;
2)因子A的贡献率75.22%,等于因子A的纯波动445.69除以总波动592.50;
3)误差的贡献率24.78%,等于100%减去A的贡献率75.22%。