minitab正交分析、响应分析

1、Minitab实验之试验设计实验目的： 本实验主要引导学生利用Minitab统计软件进行试验设计分析，包括全因子设计、部分因子设计、响应曲面设计、混料设计、田口设计以及响应优化，并能够对结果做出解释。实验仪器：Miinitabb软件、计算算机实验原理：“全因子试验设设计”的定义是：所有因子的的所有水平的的所有组合都都至少要进行行一次试验的的设计。由于于包含了所有有的组合，全全因子试验所所需试验的总总次数会比较较多，但它的的优点是可以以估计出所有有的主效应和和所有的各阶阶交互效应。所所以在因子个个数不太多，而而且确实需要要考察较多的的交互作用时时，常常选用用全因子设计计。一般情况况下，当因子子水

2、平超过22时，由于试试验次数随着着因子个数的的增长而呈现现指数速度增增长，因而通通常只作2水水平的全因子子试验。进行2水平全因因子设计时，全全因子试验的的总试验次数数将随着因子子个数的增加加而急剧增加加，例如，66个因子就需需要64次试试验。但是仔仔细分析所获获得的结果可可以看出，建建立的6因子子回归方程包包括下列一些些项：常数项项、主效应项项有6项、二二阶交互作用用项15项、三三阶交互项220项，6阶交互互项1项，除除了常数项、主主效应项和二二阶交互项以以外，共有442项是3阶阶以及3阶以以上的交互作作用项，而这这些项实际上上已无具体的的意义了。部部分因子试验验就是在这种种思想下诞生生的，它

3、可以以使用在因子子个数较多，但但只需要分析析各因子和22阶交互效应应是否显著，并并不需要考虑虑高阶的交互互效应，这使使得试验次数数大大减少。在实际工作中，常常常要研究响响应变量Y是是如何依赖于于自变量，进进而能找到自自变量的设置置使得响应变变量得到最佳佳值（望大、望望小或望目）。如如果自变量的的个数较少（通通常不超过33个），则响响应曲面方法法（respponse surfaace meethodoology，RRSM）是最最好的方法之之一，本方法法特别适合于于响应变量望望大或望小的的情形。通常常的做法是：先用2水平平因子试验的的数据，拟合合一个线性回回归方程（可可以包含交叉叉乘积项），如如果

4、发现有弯弯曲的趋势，则则希望拟合一一个含二次项项的回归方程程。其一般模模型是（以两两个自变量为为例）：这些项比因子设设计的模型增增加了各自的的变量的平方方项。由于要要估计这些项项的回归系数数，原来因子子设计所安排排的一些设计计点就不够用用了，需要再再增补一些试试验点。这种种先后分两阶阶段完成全部部试验的策略略就是“序贯试验”的策略。适适用于这种策策略的方法有有很多种，其其中最常用的的就是中心复复合设计（ccentraal commpositte dessign，CCCD）。稳健参数设计（rrobustt paraameterr desiign）（也也称健壮设计计、鲁棒设计计，简称参数数设计）是

5、工工程实际问题题中很有价值值的统计方法法。它通过选选择可控因子子的水平组合合来减少一个个系统对噪声声变化的敏感感性，从而达达到减小此系系统性能波动动的目的。过过程的输入变变量有两类：可控因子和和参数因子。可可控因子是指指一旦选定就就保持不变的的变量，它包包括产品或生生产过程设计计中的设计参参数，而噪声声因子是在正正常条件下难难以控制的变变量。在做参参数设计时，就就是把可控因因子的设计当当做研究的主主要对象，与与此同时让噪噪声因子按照照设定的计划划从而系统改改变其水平的的方法来表示示正常条件下下的变化，最最终按照我们们预定的望大大、望小或望望目地目标选选出最佳设置置。田口玄一一博士在参数数设计方

6、法方方面贡献非常常突出，他在在设计中引进进信噪比的概概念，并以此此作为评价参参数组合优劣劣的一种测度度，因此很多多文献和软件件都把稳健参参数设计方法法称为田口方方法（Tagguchi desiggn）。在实际工作中，常常常需要研究究一些配方配配比试验问题题。这种问题题常出现在橡橡胶、化工、制制药、冶金等等课题中。例例如不锈钢由由铁、镍、铜铜和铬4种元元素组成；闪闪光剂由镁、硝硝酸钠、硝酸酸锶及固定剂剂组成；复合合燃料、复合合塑料、混纺纺纤维、混泥泥土、粘结剂剂、药品、饲饲料等都是由由多种成分按按相应比例而而不是其绝对对数值；而且且显然所有分分量之和总是是为1的。对对于这种分量量之和总是为为1的

7、试验设设计，称为混混料设计（mmixturre dessign）。实验内容和步骤骤：实验之一：全因因子试验设计计:例：改进热处处理工艺提高高钢板断裂强强度问题。合合金钢板经热热处理后将提提高其断裂其其抗断裂性能能，但工艺参参数的选择是是个复杂的问问题。我们希希望考虑可能能影响断裂强强度的4个因因子，确认哪哪些因子影响响确实是显著著的，进而确确定出最佳工工艺条件。这这几个因子及及其试验水平平如下：A：加热温度，低低水平：8220，高水平平：860（摄摄氏度）B：加热时间，低低水平：2，高高水平：3（分分钟）C：转换时间，低低水平：1.4，高水平平：1.6（分分钟）D：保温时间，低低水平：500，

8、高水平：60（分钟钟）由于要细致考虑虑各因子及其其交互作用，决决定采用全因因子试验，并并在中心点处处进行3次试试验，一共119次试验。步骤1：全因子子设计的计划划（创建）选择统计=DOE=因子=创建因子设设计，单击打开开创建因子设设计对话框。, 选择两水平因子子（默认生成成元），在因因子数中选择择4，单击“设计”选项，弹出出“设计”选项对话框框。选择“全因子”试验次数为为16的那行行，并在“每个区组的的中心点数”中选择3，其其他项保持默默认（本例中中没有分区组组，各试验点点皆不需要完完全复制）。单单击确定。单击“因子”选选项打开，分分别填写四个个因子的名称称及相应的低低水平和高水水平的设置。单

9、单击确定。“选项”选项可可以使用折叠叠设计（这是是一种减少混混杂的方法）、指指定部分（用用于设计生成成）、使设计计随机化以及及在工作表中中存储设计等等；“结果”选项用于控控制会话窗口口中显示的输输出。本例中中这两项保持持默认。单击击确定，计算算机会自动对对于试验顺序序进行随机化化，然后形成成下列表格。在在表的最后一一列，写上响响应变量名（强强度），这就就完成了全部部试验的计划划阶段的工作作。步骤2：拟合选选定模型按照上图的试验验计划进行试试验，将结果果填入上表的的最后一列，则则可以得到试试验的结果数数据（数据文文件：DOEE_热处理（全全因），如如下：拟合选定模型的的主要任务是是根据整个试试验

10、的目的，选选定一个数学学模型。通常常首先可以选选定“全模型”，就是在模模型中包含全全部因子的主主效应及全部部因子的二阶阶交互效应。在在经过细致的的分析之后，如如果发现某些些主效应和二二阶交互效应应不显著，则则在下次选定定模型的时候候，应该将不不显著的主效效应和二阶交交互效应删除除。选择统计=DOE=因子=分析因子设设计，打开分析析因子设计对对话框。点击“项”选项项后，在“模型中包含含项的阶数”中选择2（表表示模型中只只包含2阶交交互作用和主主效应项，三三阶以上交互互作用不考虑虑），对默认认的“在模型中包包括中心点”保持不选。单单击确定。在“图形”选项项中，“效应图”中选择“正态”和“Paret

11、to”，“图中的标准准差”中选择“正规”，“残差图”中选择“四合一”，在“残差与变量量”图中将“加热温度”、“加热时间”、“转换时间”和“保温时间”选入，单击击确定。在“存储”选项项中，在“拟合值与残残差”中选定“拟合值”和“残差”，在“模型信息”中选定“设计矩阵”。单击确定定。结果如下：拟合因子: 强强度 与 加热温度, 加热时间, 转换时间, 保温时间 强度 的估计效效应和系数（已已编码单位） 系系数标项 效应应 系数 准误 T P常量 5441.6322 1.377 393.339 0.000加热温度 200.038 10.019 1.5000 6.68 0.0000加热时间 166.8

12、87 8.444 1.5000 5.63 0.0000转换时间 33.813 1.906 1.5000 1.27 0.2440保温时间 111.113 5.556 1.5000 3.70 0.0006加热温度*加热热时间 0.7377 00.369 1.5500 0.255 0.8812加热温度*转换换时间 -0.4877 -00.244 1.5500 -0.166 0.8875加热温度*保温温时间 3.0622 11.531 1.5500 1.022 0.3337加热时间*转换换时间 1.2633 00.631 1.5500 0.422 0.6685加热时间*保温温时间 7.1133 33.

13、556 1.5500 2.377 0.0045转换时间*保温温时间 0.8377 00.419 1.5500 0.288 0.7787S = 6.000146 PPRESS = 17778.45R-Sq = 92.499% RR-Sq（预预测） = 53.688% RR-Sq（调调整） = 83.111%强度 的方差分分析（已编码码单位）来源 自由度度 Seeq SS Adjj SS Adj MS F P主效应 4 33298.885 32298.855 8224.71 22.990 0.0002因子交互作用用 6 2522.17 252.17 42.003 11.17 0.4088残差误差

14、8 2288.144 2888.14 366.02 弯曲 1 9.92 9.992 9.922 0.25 00.633 失拟 5 169.72 169.772 33.944 0.63 00.709 纯误差 2 108.550 1108.500 554.25合计 18 3839.16强度 的估计系系数（使用未未编码单位的的数据）项 系数常量 9332.26加热温度 -0.250063加热时间 -111.2262转换时间 43.8812保温时间 -16.56637加热温度*加热热时间 0.0366875加热温度*转换换时间 -0.1211875加热温度*保温温时间 00.01533125加热时间*

15、转换换时间 12.66250加热时间*保温温时间 1.422250转换时间*保温温时间 0.833750结果分析：分析要点一：分分析评估回归归的显著性。包包含三点：（1）看方差分分析表中的总总效果。方差差分析表中，主主效应对应的的概率P值为为0.0000小于显著性性水平0.005，拒绝原原假设，认为为回归总效果果是显著的。（2）看方差分分析表中的失失拟现象。方方差分析表中中，失拟项的的P值为0.709，无无法拒绝原假假设，认为回回归方程并没没有因为漏掉掉高阶交互作作用项而产生生失拟现象。（3）看方差分分析表中的弯弯曲项。方差差分析表中，弯弯曲项对应的的概率P值00.633，表表明无法拒绝绝原假

16、设，说说明本模型中中没有弯曲现现象。分析要点二：分分析评估回归归的总效果（1）两个确定定系数R-SSq与R-Sq（调调整），计算算结果显示，这这两个值分别别为92.449%和833.11%，二二者的差距比比较大，说明明模型还有待待改进的余地地。（2）对于预测测结果的整体体估计。计算算结果显示RR-Sq和R-Sq（预测）分别为922.49%和和53.688%，二者差差距比较大；残差误差的的SSE为2288.144，PRESSS 为 17788.45，两两者差距也比比较大；说明明在本例中，如如果使用现在在的模型，则则有较多的点点与模型差距距较大，模型型应该进一步步改进。分析要点三：分分析评估各项项

17、效应的显著著性。计算结结果显示，44个主效应中中，加热温度度、加热时间间和保温时间间是显著的，只只有转换时间间不显著；66个2因子水水平交互效应应中，只有加加热时间*保保温时间是显显著的。说明明本例中还有有不显著的自自变量和2因因子交互作用用，改进模型型时应该将这这些主效应和和交互作用删删除。对于各项效应的的显著性，计计算机还输出出了一些辅助助图形来帮助助我们判断和和理解有关结结论。Pareto图图是将各效应应的t检验的的t值的绝对对值作为纵坐坐标，按照绝绝对值的大小小排列起来，根根据选定的显显著性水平，给给出t值的临临界值，绝对对值超过临界界值的效应将将被选中，说说明这些效应应是显著的。从从

18、图中可以看看到，加热温温度、加热时时间、保温时时间以及加热热时间*保温温时间是显著著的。正态效应图，凡凡是因子效应应离直线不远远者，就表明明这些效应是是不显著的；反之，则是是显著的。从从图中可以看看到，加热温温度、加热时时间、保温时时间以及加热热时间*保温温时间是显著著的。步骤3：残差诊诊断 残差差诊断的主要要目的是基于于残差的状况况来诊断模型型是否与数据据拟合得比较较好。如果数数据和模型拟拟合得比较好好，则残差应应该是正常的的。残差分析析包括四个步步骤： （1）在在“四合一”图的右下角角图中，观察察残差对于以以观测值顺序序为横轴的散散点图，重点点考察此散点点图中，各点点是否随机地地在水平轴上

19、上下无规则的的波动着。（2）在“四合合一”图的右上角角图中，观察察残差对于以以响应变量拟拟合预测值为为横轴的散点点图，重点考考察此散点图图中，残差是是否保持等方方差性，即是是否有“漏斗型”或“喇叭型”。（3）在“四合合一”图的左上角角正态概率图图（或右下角角的直方图）中中，观察残差差的正态检验验图，看残差差是否服从正正态分布。（4）观察残差差对于以各自自变量为横轴轴的散点图，重重点观察此散散点图中是否否有弯曲趋势势。 从上面这些图可可以看到，这这些图形都显显示残差是正正常的。步骤4：判断模模型是否需要要改进 这一步步需要综合前前面的分析：包括残差诊诊断和显著性性分析。从上上面的分析我我们得知，

20、在在模型中包含含不显著项，应应该予以删除除，所以需要要建立新的模模型。选择统计=DOE=因子=分析因子设设计，打开分析析因子设计对对话框。主要要是修改“项”选项中的设设置，在选取取的项中将加加热温度、加加热时间和保保温时间保留留，其他项皆皆删去，操作作中的其余各各项都保持不不变。单节确确定。结果如下：拟合因子: 强强度 与 加热温度, 加热时间, 保温时间 强度 的估计效效应和系数（已已编码单位） 系系数标项 效应应 系数 准误 T P常量 5441.3199 1.363 397.227 0.000加热温度 200.038 10.019 1.3663 7.35 0.0000加热时间 166.8

21、87 8.444 1.3663 6.20 0.0000保温时间 111.112 5.556 1.3663 4.08 0.0001加热时间*保温温时间 7.1133 33.556 1.3363 2.611 0.0022Ct Pt 1.9981 3.4299 00.58 0.5733S = 5.445038 PPRESS = 7244.350R-Sq = 89.944% RR-Sq（预预测） = 81.133% RR-Sq（调调整） = 86.077%强度 的方差分分析（已编码码单位）来源 自由度度 Seeq SS Adjj SS Adj MS F P主效应 3 33240.771 32240.7

22、11 10880.24 36.336 0.0002因子交互作用用 1 2022.35 202.35 202.335 66.81 0.0222 弯曲 1 9.92 9.992 9.922 0.33 00.573残差误差 113 3386.199 3886.19 299.71 失拟 3 151.52 151.552 50.511 2.15 00.157 纯误差 10 234.667 2234.677 223.47合计 18 3839.16强度 的估计系系数（使用未未编码单位的的数据）项 系数常量 212.788加热温度 0.5009338加热时间 -661.35000保温时间 -22.445000

23、加热时间*保温温时间 1.422250Ct Pt 1.981255结果分析：从方差分析表中中可以看到，主主效应和2阶阶交互作用对对应的概率都都小于显著性性水平0.005，应该拒拒绝原假设，认认为本，本模模型总的来说说是有效的；失拟值和弯弯曲对应的概概率分别为00.157和和0.5733，都大于显显著性水平，不不应拒绝原假假设，说明本本模型删除了了很多项之后后，并没有造造成失拟的现现象。再看删减后的模模型是否比原原来的有所改改进。从上述述表中，可以以看到，由于于模型的项数数减少了6项项，R-Sqq通常都会有有微小的降低低（本例由00.92499降到0.88968），但但关键还是要要看调整的RR-

24、Sq（调调整）是否有有所提高，本本例中，该值值从0.83311提高到到0.86773，可见删删除不显著的的效应之后，回回归的效果明明显好了；而而s的值有66.001446降为5.319133，PRESSS由17778.45降降到704.408，再再次证明删除除不明显的主主效应和交互互效应后，回回归的结果更更好了。步骤5：对选定定的模型进行行分析解释经过前三步的多多次反复以后后，我们可以以获得一个满满意的回归方方程：对选定的模型进进行分析，主主要是在拟合合选定模型后后输出更多的的图形和信息息，并做出有有意义的解释释。主要包括括下面四个方方面：（1）再次进行行残差诊断。具体做法是：选选择统计=DO

25、E=因子=分析因子设设计，打开分析析因子设计对对话框。点击击“图形”窗口后，在在“图中的残差差”中选择“标准化”，在“残差图”中，在单独独视图下选择择“直方图”，单击确定定。点击“存储”窗口后，在在拟合值与残残差中，选择择“标准化残差差”和“删后”。单击确定定。结果如下：从得出的直方图图可知，残差差及所有残差差数据都是正正常的。（2）确认主效效应及交互作作用的显著性性，并考虑最最优设置通过输出各因子子的主效应图图和交互效应应图来判定。具具体做法是：选择统计计=DOE=因子=因子图，打开因子子图对话框。选选定“主效应图”和“交互作用图图”，在图中使使用的均值类类型中选择“数据均值”。在主效应应图

26、的设置中中，将“强度”选入到响应应中，将可用用中的所有项项选入所选中中；在交互作作用图的设置置中，重复前前面主效应图图设置的步骤骤。单击确定定。结果如下：从主效应图中可可以看到，加加热温度、加加热时间和保保温时间三者者的回归线比比较陡，顾主主效应影响确确实显著，而而转换时间的的回归线较平平，故主效应应影响不显著著；为了使断断裂强度达到到最大，三因因子都是取值值越大越好，即即加热温度应应取上限8660摄氏度，加加热时间应取取上限3分钟钟，保温时间间应取上限660分钟。从从交互作用图图可以看出，只只有加热时间间和保温时间间二者效应线线明显不平行行，说明二者者交互作用显显著。（3）输出等值值线图、响

27、应应曲面图等以以确认最佳设设置本例中，只有加加热时间和保保温时间的交交互作用显著著，因此绘制制这组等值线线图和响应曲曲面图，而设设定另一个影影响显著的变变量（加热温温度）为最佳佳设置。具体体操作为：选选择统计=DOE=因子=等值线/曲曲面图，打开等值线线/曲面图对对话框。选定定“等值线图”和“曲面图”。在等值线线图设置中，在在因子中，XX轴选为加热热时间，Y轴轴选为保温时时间，在设置置中，选择保保留附加因子子在高设置，并并在加热时间间中设置8660，单击确确定；在曲面面图设置中，XX轴中选择加加热时间，YY轴中选择保保温时间，单单击确定。结果如下： 从等值线图和和曲面图可以以看出，断裂裂强度的

28、最大大值确实在加加热时间为33分钟，保温温时间为600分钟，加热热温度固定在在860摄氏氏度时达到最最大。（4）实现最优优化Minitabb软件中有专专门的响应变变量优化器窗窗口。具体做做法：统计计=DOE=因子=响应优化器器，打开响应优优化器对话框框。将“可用项”中的强度选选入到“所选项”中；点击“设置”窗口，根据据本例的要求求，在“目标”中选择“望大”，在“下限”中填入5660（这个值值是在做过的的试验中已经经实现了的），在在“望目”中填入6000（这个值值是在做过的的试验中未能能达到的，是是较高理想），上上限留为空白白。结果如下：这个图中共有33列，分别为为选中的自变变量。最上端端列出各

29、变量量的名称、取取值范围以及及最优设置，上上半图是合意意值d的取值值情况，下半半图是最优化化结果：最大大值在加热温温度取8600摄氏度、加加热时间取33分钟、保温温时间取600分钟达到，断断裂强度最终终可以达到5569.20066。合意意度d为0.230166。步骤6：进行验验证试验通常的做法是在在先算出在最最佳点的观测测值的预测值值及其变动范范围，然后再再最佳点做若若干次验证试试验，如果验验证试验结果果的平均值落落在事先计算算好的范围内内，则说明一一切正常，模模型是正确的的，预测结果果可信；否则则就要进一步步分析发生错错误的原因，改改进模型，再再重新验证，以以求得符合实实际数据的统统计模型。

30、具具体做法是：选择统计计=DOE=因子=分析因子设设计，打开分析析因子设计对对话框。在前前面建立的模模型的基础上上，即在“项”中已经将最最终选定的模模型中包括了了加热温度、加加热时间、保保温时间以及及加热时间和和保温时间的的交互作用项项。再打开“预测”窗口，在“因子”中按顺序设设定各个主效效应的最优值值，分别为8860 3 60。单击击确定。结果如下：根据该模型在新新设计点处对对 强度 的预测响应应 拟合值点 拟合值值 标准误 95% 置信区间间 955% 预测区区间 1 5669.2077 2.926 (562.931, 575.4483) (556.186, 582.2227)结果解释：最

31、左左侧给出的拟拟合预测值是是569.2207，就是是将自变量值值代入回归方方程所得的结结果，这与最最优值的预测测是一致的。拟拟合值标准误误为2.9226，是拟合合值的标准差差，此值在作作进一步计算算时还有用。预预测值平均值值置信区间的的结果是（5562.9331,5755.438），具具体的理解可可以是：当加加热温度取8860摄氏度度，加热时间间取3分钟，保保温时间取660分钟时，我我们有95%的把握断言言，断裂强度度平均值将落落入（5622.931,575.4438）之内内。95%的的预测区间是是将来一次验验证试验时将将要落入的范范围，可供做做验证试验时时使用，具体体的理解是：当加热温度度取

32、860摄摄氏度，加热热时间取3分分钟，保温时时间取60分分钟时，我们们有95%的的把握断言，任任何一块钢板板的断裂强度度将落入（5556.1886,5822.227）之之内。试验之二：部分分因子试验设设计部分因子试验设设计与全因子子试验设计的的不同之处在在于大大减少少了试验的次次数，具体表表现在试验设设计创建阶段段的不一致，下面主要就部分因子试验设计的创建进行讲述。步骤1：部分因因子试验的计计划（创建）默认生成元的计划例：用自动刨床床刨制工作台台平面的工艺艺条件试验。在在用刨床刨制制工作台平面面试验中，考考察影响其工工作台平面光光洁度的因子子，并求出使使光洁度达到到最高的工艺艺条件。共考察6个

33、因子子：A因子：进刀速速度，低水平平1.2，高高水平1.44（单位：mmm/刀）B因子：切屑角角度，低水平平10，高水水平12（单单位：度）C因子：吃刀深深度，低水平平0.6，高高水平0.88（单位：mmm）D因子：刀后背背角，低水平平70，高水水平76（单单位：度）E因子：刀前槽槽深度，低水水平1.4，高高水平1.66（单位：mmm）F因子：润滑油油进给量，低低水平6，高高水平8（单单位：毫升/分钟）要求：连中心点点在内，不超超过20次试试验，考察各各因子主效应应和2阶交互互效应AB、AAC、CF、DDE是否显著著。由于试验验次数的限制制，我们在因因子点上只能能做试验166次，另4次次取中心

34、点，这这就是的试验，通过过查部分因子子试验分辨度度表可知，可可达分辨度为为的设计。具具体操作为：选择 统统计=DOE=因子=创建因子设设计，单击打开开创建因子设设计对话框。在在“设计类型”中选择默认认2水平因子子（默认生成成元），在“因子数”中选定6。单击“显示可用用设计”就可以看到到下图的界面面，可以确认认：用16次次试验能够达达到分辨度为为的设计。单击“设计”选选项，选定11/4部分实实施，在每个个区组的中心心点数中设定定为4，其他他的不进行设设定，单击确确定。单击“因子”选选项，设定各各个因子的名名称，并设定定高、低水平平值。点击确确定。再点击确定后，就就可以得到试试验计划表，如如下：与

35、全因子设计不不同的是，我我们不能肯定定这个试验计计划表一定能能满足要求，因因为部分因子子试验中一定定会出现混杂杂，这些混杂杂如果破坏了了试验要求，则则必须重新进进行设计，从从运行窗中可可以看到下列列结果：设计生成元: E = AABC, FF = BCCD别名结构I + ABCCE + AADEF + BCDFFA + BCEE + DEEF + AABCDFB + ACEE + CDDF + AABDEFC + ABEE + BDDF + AACDEFD + AEFF + BCCF + AABCDEE + ABCC + ADDF + BBCDEFF + ADEE + BCCD + AABCE

36、FAB + CEE + ACCDF + BDEFAC + BEE + ABBDF + CDEFAD + EFF + ABBCF + BCDEAE + BCC + DFF + ABBCDEFAF + DEE + ABBCD + BCEFBD + CFF + ABBEF + ACDEBF + CDD + ABBDE + ACEFABD + AACF + BEF + CDEABF + AACD + BDE + CEF从此表得知，计计算机自己选选择的生成元元是：E=AABC,F=BCD。后后面的别名结结构中列出了了交互作用项项的混杂情况况，即每列中中互为别名的的因子有哪些些；从上表可可以看出，主主效应

37、与三阶阶及四阶交互互作用混杂，二二阶交互作用用与四阶交互互作用混杂，三三阶交互作用用与四阶交互互作用混杂；关键是要检检查一下题目目所要求的22阶交互作用用情况，将33阶以上的交交互作用忽略略不计，混杂杂的情况有： AB=CCE,AC=BE,ADD=EF, AF=DEE,AE=BBC=DF,BD=CFF,BF=CCD。本例中中所要求的44个2阶交互互作用是ABB,AC,CCF,DE，显显然可以看到到，这四个22阶交互作用用均没有混杂杂。因此可以以看到此试验验计划是可行行的。步骤2：指定生生成元的部分分因子试验计计划例：和前面的例例子是一样的的，考察的是是各因子主效效应和2阶交交互效应ABB,AC

38、,CCE和DE是是否显著。从上例的别名结结构表中可以以看出，ABB与CE是相相互混杂，因因此用默认的的生成元构造造的试验计划划是不能满足足要求的。指定生成元的步步骤：由要求求条件可知，AAB,AC,CE和DEE不能混杂，这这相当于ABBCE,ABBDE,ACDE,运用用移项法则，变变形后可知，即即EABC,EEABD,EEACD.对对于分辨度为为的设计生成成元中，只能能含3个字母母。而试验次次数为16的的的各列中，字字母个数为33的项只有44个：ABCC,ABD,ACD以及及BCD。既既然给定条件件中有3个选选择不可接受受，因此，生生成元只能选选择E=BCCD，试验计计划对于F没没有要求，因因

39、此F可以任任选，取F=ABC。具体操作为：选选择 统计计=DOE=因子=创建因子设设计，单击打开开创建因子设设计对话框。在在“设计类型”中选择2水水平因子（指指定生成元），在在“因子数”中选定4（这这是基本设计计的因子数，其其他两个因子子是通过指定定生成元加入入的）。打开“因子”对对话框，选定定全因子，并并在“每个区组的的中心点数”中选择4。打打开“生成元”选项，在“通过列出生生成元将因子子添加到基本本设计中”中填写生成成元：E=BBCD F=ABC,单单击确定。单击确定后，得得到的结果如如下：设计生成元: E = BBCD, FF = ABBC别名结构（直到到 4阶项）I + ABCCF +

40、 AADEF + BCDEEA + BCFF + DEEF + AABCDEB + ACFF + CDDE + AABDEFC + ABFF + BDDE + AACDEFD + AEFF + BCCE + AABCDFE + ADFF + BCCD + AABCEFF + ABCC + ADDE + BBCDEFAB + CFF + ACCDE + BDEFAC + BFF + ABBDE + CDEFAD + EFF + ABBCE + BCDFAE + DFF + ABBCD + BCEFAF + BCC + DEE + ABBCDEFBD + CEE + ABBEF + ACDFBE

41、+ CDD + ABBDF + ACEFABD + AACE + BEF + CDFABE + AACD + BDF + CEF从上面的结果可可以看出，AAB,AC,CE和DEE均没有相互互混杂，此设设计满足原定定的要求。部分因子试验的的分析步骤总总体来说与全全因子试验设设计是一致的的。但是有一一个要注意的的地方：在第第一步选定模模型中显著的的主效应和22阶交互作用用时，当某些些2阶交互作作用效用显著著时，不能仅仅从表面上的的结果来定取取舍，要仔细细分析混杂结结构，查看在在结构表中，此此显著项是与与哪个（或哪哪些）2阶交交互作用效应应相混杂的，再再根据背景材材料予以判断断，最终决定定入选。比如

42、如：数据显示示B,C,DD以及AD是是显著的，但但是背景材料料又说明A和和D没有交互互作用，而AAD与BC是是相混杂的，这这个时候，应应该是B,CC,D以及BBC是显著的的。实验之三：响应应曲面设计 响应曲面设计包包括两种方法法：中心复合合设计和Boox-Behhnken设设计。在中心心复合设计中中，整个试验验由下面三部部分试验点组组成：（1）立方体点点或称角点点，各坐标标皆为1或-1，这是因因子试验的组组成部分；（2）中心点，各各点坐标皆为为0；（3）星号点或轴点，除除了一个自变变量的坐标为为外，其余余自变量皆为为0，在k个个因子的情况况下，共有22k个星号点点。中心复合设计，包包括三种设计

43、计：中心复合序贯设设计，当“水平定义”栏中选定“立方点”时，表示这时时设定的水平平作为立方点点，星号点将将超出立方体体。中心复合有界设设计，当“水平定义”栏中选定“轴点”时，表示这这时设定的水水平作为轴上上的星号点，立立方点将向内内收缩。中心复合表面设设计，意味着着将星号点的的位置向中心心收缩而设定定在立方体的的表面上。Box-Behhnken设设计，这种设设计是将各试试验点取在立立方体的棱的的中点上，除除非极端重视视试验次数，否否则通常不采采用这种设计计。步骤1：响应曲曲面的计划例：提高密封胶胶条黏合力试试验。影响黏黏合力的3个个因子是：AA:烘烤温度度（220-240摄氏氏度）、B:烘烤时

44、间（66-10秒）、C:黏合压力（100-140帕）。在因子设计中，分别取下列条件，安排了全因子试验：A:烘烤温度,低水平2220，高水平平240（摄摄氏度）B:烘烤时间，低低水平7，高高水平9（秒秒）C:黏合压力，低低水平1100，高水平1130（帕）试验后发现，数数据明显呈现现弯曲状况，希希望进一步安安排些实验以以拟合响应曲曲面方程。由由于要进行序序贯试验，最最好选中心复复合设计。具体做法是：选选择统计DOEE响应应曲面创建响应曲曲面设计，打打开创建响应应曲面设计对对话框。在“设计类型”中选择“中心复合”，在“因子数”中设定为33。打开“显示可用用设计”对话框，可可以看到未划划分区组时试试

45、验次数为220。打开“设计”后，本例中中需要的试验验次数为200次，这是可可行的，因此此不必修改，中中心点数也不不用另设；但但是选取哪种种中心复合设设计，需要考考虑更多条件件，由于在烘烘烤温度上，原原来的试验温温度条件已经经取在边界上上了，不允许许再超界因而而不能使用中中心复合序贯贯设计，但是是又考虑到要要保持序贯性性，只能放弃弃中心复合边边界设计（没没有序贯性），因而选用用中心复合表表面设计，即即在Alphha值中选择择表面中心；在“因子”中，选择“立方点”，并填写各各因子的名称称及水平；在在“选项”中，为了看看清楚结构，暂暂时先删除随随机化。单节节确定。结果如下：在表中的20次次试验中，第

46、第1至第8号号因子点以及及第15至117号中心点点，已经在因因子设计阶段段获得了数据据，只要将这这些结果填在在后面第8列列上，然后再再补充其他99个试验，及及可以完成全全部响应曲面面的试验任务务。步骤2：响应曲曲面设计的分分析例：提高烧碱纯纯度问题。在在烧碱生产过过程中，经过过因子的筛选选，最后得知知反应炉内压压力及温度是是两个关键因因子。在改进进阶段进行全全因子试验，因因子A压力的的低水平和高高水平分别取取为50帕和和60帕，因因子B反应温温度的低水平平和高水平分分别取为2660及3200摄氏度，在在中心点处也也作了3次试试验，试验结结果在数据文文件：DOEE_烧碱纯度度（响应1）。对于这批

47、数据按按全因子试验验进行分析，具具体操作为：选择统计计=DOE=因子=分析因子设设计，打开分析析因子设计对对话框。首先先将全部备选选项列入模型型，删除在模模型中包括中中心点，在“图形”中的残差与与变量下将压压力和温度选选入进去。得得到的结果如如下：纯度 的效应和和系数的估计计（已编码单单位）项 效应 系数 系数标准准误 T P常量 96.9661 0.44150 233.663 0.000压力 -2.665 -1.3332 0.55490 -2.443 0.094温度 -0.765 -0.3882 0.55490 -0.770 0.536压力*温度 0.0335 00.018 0.54990

48、0.03 0.9777S = 1.009803 PPRESS = 1344.203R-Sq = 68.011% RR-Sq（预预测） = 0.00% R-Sq（调整整） = 336.01%对于 纯度 方方差分析（已已编码单位）来源 自由度 Seq SSS AAdj SSS Addj MS F P主效应 2 77.68744 7.668745 3.844372 3.199 0.11812因子交互作用用 1 0.00012 00.001223 0.001233 0.00 00.977残差误差 33 3.6170 3.611701 1.205567 弯曲 1 3.51778 3.517811 3.5

49、51781 70.992 0.014 纯误差 2 00.09922 0.009920 0.044960合计 6 111.30557从上述表中可以以看到，主效效应和2因子子交互作用对对应的概率PP值均大于00.1，说明明模型的总效效应不显著，而而且弯曲对应应的概率P值值为0.0114，拒绝原原假设，认为为存在明显的的弯曲趋势；R-Sq和R-Sq（预预测）的值都都比较小，说说明了模型的的总效果不显显著。从残差与各变量量的图也验证证了存在严重重的弯曲现象象。这些都表表明，对响应应变量单纯地地拟合一阶线线性方程已经经不够了，需需要再补充些些“星号点”，构成一个个完整的响应应曲面设计，拟拟合一个含二二阶

50、项的方程程就可能问题题了。补充的的4个星号点点的实验结果果见数据表：DOE_烧烧碱纯度（响响应2）。下面对全部111个点构成的的中心复合序序贯设计进行行分析，拟合合一个完整的的响应曲面模模型。分析如如下：第一步：拟合选选定模型。选择统计DOE响应曲曲面分分析响应曲面面设计，打打开分析响应应曲面设计对对话框。点击击窗口“项”以后，可以以看到模型中中将全部备选选项都列入了了模型，包括括A(压力)、B（温度度）以及它们们的平方项AAA、BB和和交互作用项项AB；打开开“图形”窗口，选定定“正规”、“四合一”以及残差与与变量，并将将压力和温度度都选入残差差与变量中；打开“储存”窗口，选定定“拟合值”、

51、“残差”以及“设计矩阵”。单击确定定。得到的结果如下下：纯度 的估计回回归系数项 系数 系数标准准误 TT P常量 97.7804 00.105002 9331.0666 0.0000压力 -1.8911 00.091114 -220.7500 0.0000温度 -0.6053 00.090992 -6.6577 0.0001压力*压力 -2.58822 0.115339 -16.835 0.0000温度*温度 -0.46615 0.115314 -3.014 0.0300压力*温度 0.03351 0.118253 0.192 0.8555S = 0.11819000 PRRESS = 0.

52、6993667R-Sq = 99.355% R-Sq（预测测） = 997.27% R-SSq（调整） = 98.770%对于 纯度 的的方差分析来源 自由由度 SSeq SSS Addj SS Adjj MS F P回归 5 25.22310 25.23310 55.046220 1552.51 0.0000 线性 2 15.77127 15.71127 77.856335 2337.44 0.0000 平方 2 9.55171 9.51171 44.758553 1443.82 0.0000 交互作用 1 0.00112 00.00122 0.000123 0.04 00.855残差误差

53、5 0.16554 00.16544 0.003309 失拟 3 0.00662 0.06662 00.022008 0.45 0.7447 纯误差 2 0.09992 0.09992 0.049600合计 10 25.33964结果解释：（1）看方差分分析表中的总总效果。在本本例中，回归归项的P值为为0.0000，表明应该该拒绝原假设设，认为本模模型总的来说说是有效的。看方差分析表中中的失拟现象象，本例中，失失拟项对应的的P值为0.747，明明显大于显著著性水平0.05，接受受原假设，认认为本模型中中不存在失拟拟现象。（2）看拟合的的总效果。本本例中，R-Sq与R-Sq（调调整）比较接接近，

54、认为模模型的拟合效效果比较好；R-Sq（预预测）比较接接近于R-SSq值且这个个值比较大，说说明将来用这这个模型进行行预测的效果果比较可信。（3）各效应的的显著性。从从表中可以看看到，压力、温温度以及它们们的平方项对对应的概率值值都小于显著著性水平，说说明这些效应应都是显著的的；而压力和和温度的交互互效应项对应应的概率值为为0.8555，显然大于于显著性水平平，认为该效效应项是不显显著的。第二步：进行残残差诊断利用自动输出的的残差图来进进行残差诊断断。从上述残差图中中可以看出，残残差的状况是是正常的。第三步：判断模模型是否需要要改进。根据第一步的分分析，我们得得知压力和温温度的交互作作用项是不

55、显显著的，应该该予以剔除，因因此需要重新新拟合新的模模型，使得新新的模型中不不包含交互作作用项。具体实现步骤是是：在项中将将交互作用项项剔除，在结结果中输出标标准化残差和和删后残差。得到的结果为：纯度 的估计回回归系数项 系数 系数标准准误 T P常量 97.7804 00.096222 10016.1777 0.000压力 -1.8911 00.083550 -22.6447 0.000温度 -0.6053 00.083331 -7.2665 0.000压力*压力 -2.58822 0.114054 -188.373 0.0000温度*温度 -0.46615 0.114031 -33.289 0.0117S = 0.11666655 PRRESS = 0.5446550R-Sq