试验设计讲义

1、Design Of ExperimentsDOE培训纲要 TOC o 1-5 h z HYPERLINK l bookmark34 o Current Document 1简介 4 HYPERLINK l bookmark37 o Current Document 什么是DOE? 4 HYPERLINK l bookmark40 o Current Document 试验设计的作用：4 HYPERLINK l bookmark47 o Current Document 试验设计的分类：4 HYPERLINK l bookmark53 o Current Document 计划和实施试验的系统方

2、法：5 HYPERLINK l bookmark56 o Current Document 试验的基本原则：6 HYPERLINK l bookmark63 o Current Document 一个例子：8 HYPERLINK l bookmark66 o Current Document 2常用的试验设计在实际工作中个的应用范围及步骤： 9 HYPERLINK l bookmark70 o Current Document 析因阶段9 HYPERLINK l bookmark73 o Current Document 分析探索阶段9 HYPERLINK l bookmark76 o Cur

3、rent Document 验证阶段及改进阶段9 HYPERLINK l bookmark79 o Current Document 3常用的统计分析方法： 9 HYPERLINK l bookmark82 o Current Document 显著性检验9 HYPERLINK l bookmark86 o Current Document 回归及相关分析10 HYPERLINK l bookmark96 o Current Document FULL FACTORIAL DESIGN 完全析因试验 11 HYPERLINK l bookmark99 o Current Document 2的

4、k次完全析因设计 11 HYPERLINK l bookmark112 o Current Document 3的k次完全析因设计 13 HYPERLINK l bookmark117 o Current Document FACTIONAL FACTORIALDESIGN 部分析因试验 16 HYPERLINK l bookmark120 o Current Document 试验准则16 HYPERLINK l bookmark125 o Current Document 一个实例17 HYPERLINK l bookmark129 o Current Document 6响应曲面法: 2

5、0 HYPERLINK l bookmark135 o Current Document 基础知识：21 HYPERLINK l bookmark140 o Current Document 中心复合设计(central COMPOSITE DESIGN)： 216.3.8统一设计和正交设计.23 HYPERLINK l bookmark148 o Current Document Box-Behnken 设计： 23 HYPERLINK l bookmark151 o Current Document 田口设计(TAGUCHIDESIGN) 23 HYPERLINK l bookmark15

6、4 o Current Document 试验目的：23试验输入的分类： 24信噪比(SIGNAL-TO-NOISE RATIO 简称 SN): 24田口设计的两步程序：24田口设计的适用范围：24一个例子25 HYPERLINK l bookmark164 o Current Document 混合设计(MIXTURE DESIGN) 25 HYPERLINK l bookmark167 o Current Document 适用范围：25 HYPERLINK l bookmark170 o Current Document 混合设计的参数：25 HYPERLINK l bookmark17

7、6 o Current Document 后记 251简介1.1什么是DOE ?DOE是翻译为试验设计。1.2试验设计的作用：处理比较。主要目的是比较几种不同的处理并 选择最好的。变量筛选。一个系统中有许多变量，但是通常 只有一小部分是重要的，筛选试验可以用来识别 这些重要变量。响应曲面探察。当少数的重要变量被识别出来 后，它们在响应上的效应就需要深入探察。系统优化。在许多研究中，兴趣在于优化系统。 若一个响应曲面被识别了，它就可以用来进行优 化。1.2.5系统文件型。在质量改进中提高系统的抗干扰水平组合来达到目的）。1.3试验设计的分类:1.3.1若目的为处理比较使用常用的统计分析方法。1.

8、3.2若目的为变量筛选。用完全析因试验或部分析 因试验(full factorial design 和 factional factorial design)。1.3.3若目的为响应曲面探察。用响应曲面设计 (Response Surface Design)。1.3.4若目的为了在质量改进中提高系统的抗干扰 能力，则使用田口设计(Taguchi)设计。1.3.5若为了配比混合物则可以使用混合设计 (Mixture Design )。1.4计划和实施试验的系统方法：1.4.1目标阐述。1.4.2选择响应。所谓响应就是实验的观察或结果。 所选的响应可以根据阐述的目标分类。主要考虑 望目(nomin

9、al-the-best)、望大(larger the-best)和望小(smaller-the-better)o1.4.3选择因子和水平。因子(factor )试验中所研 究的变量。为了研究因子对响应的影响，需要用 到因子的两个或多个不同的取值，这些取值称为 因子的水平(level )或设置(Setting )。处理 (treatment )是指因子的水平组合，当只有一个 因子时候，它的水平就是处理。因子可以是定量 的也可以是定性的。如温度、时间、压力等可以 在一个连续区域内取值的因子叫定量因子。定性 因子是在离散点上取值的因子。一般来说，定量 因子的水平值布置应该足够分散，这样效应才可 能检

10、测出来，但是也不要布置太分散以至与各种 物理机械因素都会包括进来。对于定性因子的水 平选择或缺乏灵活。1.4.4选择试验计划。1.4.6分析数据。1.4.7推断和建议1.4.8有时需要做验证性的试验(Confirmation experiment)，以验证所推荐的水平组合是否真的 好。有时需要给是否要做跟随试验(follow-up experiment)的建议。例如，若根据试验结果显 示两组模型所得出的结果同样的好，单只能选择 一个作为最优模型，那么必须做一个跟随试验。1.5试验的基本原则：重复(replication)、随机化(randomization) 和分区组(blocking)。1.

11、5.1试验单元(experimental unit)是试验载体的通用术语，是指诸如材料、动物、人、机器或 时间段等这样的在试验中一个处理所施用于的基 本单元。所谓重复(replication)，是指每个处 理都使用于数个试验单元，这些单元是试验推断 都适用的单元总体的代表。它可以使得试验误差范围得以估计(即误差方差的估计)，并以此判断出处理间的差异。增加重复的次数可以减小处理匹出处理间的差异。增加重复的次数可以减小处理匹1=1效应估计的方差，并可更有效地检测出处理间的 差异。不同单元重复和同单元重复加以区别。来 自同一试验单元的三个读数则是同一个单元的重 复，而来自不同的试验单元的读数上不同单

12、元的 重复。因而称前者为同单元重复试验(为区别起 见，也称复试)，而称后者为不同单元重复试验。1.5.2随机化(randomization)它用于试验中的处 理在单元上的分配，也用于处理在执行试验中使 用的顺序和响应的测试顺序。它可以抵抗那些试 验者未知但可能会对响应产生影响的变量。它也 可以减少在处理分配上不希望的由主观判断而造 成的影响。此外，随机化还可保证试验误差估计 的合理性，从而为分析试验进行判断提供一个基 础。一组齐性(homogeneous)单元称为一个区组(block )。要使分区组有效，需把实验单元进行 合理的安排使得区组内单元的差异远比区组间单 元的差异小。通过在同一区组内

13、比较处理，区组 效应在处理效应比较中得以消除，从而使得试验 更有效。1.5.4若分区组有效，则其就可以由来消除组与组之 间的差异，然后把随机化用于在区组内处理到单 元的安排，以进一步减小未知变量的影响。1.6 一个例子：假设有两种键盘(A、B)，需要比较它们的打字效率, 现在有6份不同的稿件给同以打字员打印。原先的试验设计1：打印顺序如下：1。A，B。2。A，B。 3。 A， B。 4。 A， Bo 5。 A， B。 6。 A， B。尽管此试验重复了 6次。也用了分区组在同一份稿件 上来比较两种键盘，但这个设计存在一个严重的缺陷， 在A键盘上打完稿件以后，当打字员再在B键盘上打 同一份稿件时已

14、经熟悉该稿件的内容。这种“学习效 应”将会不公平地帮助B键盘的表现。优化的试验设计：1o A，Bo 2。B，Ao 3。A，B。4。B，Ao 5。 A， B。 6o B， Ao2常用的试验设计在实际工作中个的应用范围及步骤:2.1析因阶段2.1.1完全析因设计2.1.2不完全析因设计2.2分析探索阶段2.2.1响应曲面法2.2.2混合设计2.3验证阶段及改进阶段田口设计3常用的统计分析方法：3.1显著性检验配对 t 检验(Matched Pairs)o3.1.2两样本均数比较(非配对t检验)(MatchedPairs)o3.1.3方差分析3.1.4卡方检验3.1.5非参检验3.1.6检验适用范Y

15、的类型X的个数X的类型检验方式1检验方式2连续单个二分类T检验配对t检验配对T检单因素方差连续单个多分类验分析连续多个分类方差分析分类单个二分类卡方检验分类单个多分类卡方检验3.2回归及相关分析3.2.1线性回归Logistic 回归3.2.3相关分析3.2.4：适用范围应变量自变量个数变量类型检验方式连续单个连续线性回归连续多个连续相关分析分类连续连续Logistic 回归分类多个分类Logistic 回归分类多个连续Logistic 回归4 full factorial design 完全析因试验2的k次完全析因设计2的k次完全析因设计指它研究二水平因子的所有可能组合。其关键性质是均衡性（

16、balance ） 和正交性（orthogonality）。均衡性指的是每个 因子的水平出现的次数相同。正交性指的是所有因子水平组合出现的次数相同。一个设计是正交 的如果它的所有因子对都是正交的。4.1.2望目特征。有一个工程设计要求的目标值。4.1.3般步骤:4.1.3般步骤:4.1.3.1选择某些因子的水平变化最小。4.1.3.2选择一个不再（1）中的因子使水平接近 目标值。调节因子(adjustment factor)：步骤2中的因子称为调节因子，如果它对均值有显 著的影响而对方差没有影响。调节因子的存在使得 利用两步程序最小化二次损失成为可能。否则，在 步骤2中调节因子的水平时均值和方

17、差都会发生变 化，这时为了减小方差还需要将步骤1设定的因子 水平进行调整；因此步骤1和步骤2之间可能需要 进行几次重复。由此可见，如果没有调节因子的话，两步程序上难以奏效的。一个例子晶体外延层生长试验。外延层的厚度目标是14.5微米。其误差为0.5微 米。即实际的厚度尽可能接近14.5微米并在14-15 微米范围内。试验因子：感应旋转方法，喷嘴位置，沉积温度，沉积时间。因子水平-+A感应旋转方法连续震荡B喷嘴位置26C沉积温度12101220D沉积时间低高4.1.6 主效应(main effect)：4.1.7交互效应为度量两个因子，比如A和B的联合效应，A*B交 互效应(也称为交互作用)(i

18、nteraction effect 或 interaction)o4.1.8因子效应的基本原则：效应排序原则(Effect Hierarchy Principle)低阶效应比高阶效应更重要。 同阶效应的重要性是相同的。效应稀疏原则(Effect Sparsity Principle )在析因试验中相对重要的因子效应数目是小的效应遗 传原则(Effect Heredity Principle)要使一个交互作用是显著的，至少它的严格亲本 因子应该是显著的。4.1.9利用回归模型分析试验结果。3的k次完全析因设计3的k次完全析因设计。每个因子都有三个水 平因子，因而有两个自由度而且对交互效应进行 参

19、数化的常用系统。即正交成分系统(orthogonal components system )o4.2.2望大和望小特征问题及其解决步骤。 望大的两步程序：1寻找最大化均值的因子设置。2寻找其它的最小化标准差的因子设置。望小的两步程序：1寻找最小化均值的因子设置。2寻找其它的最小化标准差的因子设置。要使任意两步程序有效，步骤2中所使用的因子应 对步骤2中考虑的特征有一个显著的影响但是对1 中的特征影响不显著。4.2.3考虑使用三水平的主要原因。4.2.3.1定量因子如温度可能非单调的形式影响响应，仅有两个温度设置不能分析温度与响应 间的曲线关系，要探测一个曲线效应需要三个 或更多的设置。4.2.

20、3.2定性因子可能有几个水平，如生产某种零件的机器有三种类型等。如果当前的生产过 程中三种机器类型都要用到，那么在研究中它 们都应该进行比较。人们不能仅仅基于两种机器类型的研究来比较三种机器类型。而且，如 果首先在二水平试验中研究两种机器类型然后在跟随试验中在增加第三种类型，这样会更匹后在跟随试验中在增加第三种类型，这样会更匹1=1浪费更多的时间和财力。4.2.3.3如果一个产品或工艺的目标设置或当前 设置，比如一个零件的尺寸或者注射压力，通 常会研究其附近的两个设置。如果目标时通过 改变X来改善产品或工艺的性能，并且低一些 或高一些的设置可能实现这个目标，此时这种浪费更多的时间和财力。4.2

21、.3.3如果一个产品或工艺的目标设置或当前 设置，比如一个零件的尺寸或者注射压力，通 常会研究其附近的两个设置。如果目标时通过 改变X来改善产品或工艺的性能，并且低一些 或高一些的设置可能实现这个目标，此时这种方法特别适用。4.2.4 一个例子。安全带试验 两个质量特征：卷曲抗张强度（磅），它的最小值4000 磅;:瞬，时制动（毫米）。它是卷曲过程中的过剩金属，不能超过14.00毫米。强度必须最大化 保证安全带在事故中不会断裂，瞬时制动效应英 最小化因为它会引起后一个工序阶段装配的塑 料护套破裂。起初建议的设置为压力1700磅， 卷曲长度23毫米。试验的目标是在上述两个约 束条件下寻找最大化抗

22、张强度和最小化瞬时制 动的因子设置。总的期望目标是：强度均值为 7000磅，标准差为250磅;瞬时制动均值为13.50 毫米，标准差为0.2毫米。其设计是望大（Larger-the-better）和望小 （Small-the-better）特征问题。因子水平012A压力（磅/平方英寸）110014001700B模板（毫米）10.010.210.4C卷曲长度（毫米）182327D制动位置P74P75P764.2.5多个响应的分析策略：限外概率。限外概率(out-of-spec probability )5 factional factorial design 部分析因试验5.1试验准则5.1.1

23、定义关系(Defining Relation)即因子选取之 间的定义关系式。5.1.2字长：关系式的长度。5.1.3分辨度(resolution)：即定义字中最短的字 长5.1.4纯净的(Clear)：个主效应或二因子交互效 应如果它不与其它主效应和二因子交互效应别 名。5.1.5强纯净的(Strongly Clear)一个主效应或者 二因子交互效应如果它不与其它主效应和二因子 交互效应别名也不与三因子交互效应别名。最大分辨度(maximum resolution)准则。若 无先验知识，则选取分辨度较大的一种试验设计。最小低阶混杂准则(minimum aberration cirterion)

24、对任意两个2的k-p次设计D1和D2。 如果D1的低阶混杂比D2小，则成为D1有最小的低阶混杂。若无先验知识，则选取分低阶混杂较 小的一种试验设计。一个实例叶形弹簧试验。因子水平-+B高温加热温度（F）18401880C加热时间（秒）2325D传送时间（秒）1012E成形时间（秒）23成形时间（秒）130-150150-17016（2的4次）个水平组合中研究5个因子。称为 2的5-1次设计。其中5是因子个数，因为因子E 被安排到B，C，D的乘机上。即E是B，C，D列的 乘积。别名关系式子是：E=BCD或I二BCDE （表示截 距与BCDE交互作用的别名）I=I=BCDE称为定义关 系（Defi

25、ning Relation）其字长为4。所以具有分 辨度4。因子BCDEQ+其别名关系式:B=CDE, C=BDE, D=BCE, E=BCD,BC=DE, BD=CE, BE=CD,Q=BCDEQ, BQ=CDEQ, CQ=BDEQ, DQ=BCEQ,EQ=BCDQ, BCQ=DEQ, BDQ=CEQ, BEQ=CDQ由上式可看出。BC=DE则只有当先验知识认为DE可 以忽略的时候BC才是可估的。主效应B, C, D, E是纯净的。主效应Q和二因子 交互效应BQ, CQ, DQ, EQ是强纯净的。另一个部分析因设计：Q列为B, C, D, E列的乘积， 即Q=BCDE或I二BCDEQ。其定义

26、关系式长为5。所以 为分辨度为5的设计。因子BCDEQ+-+-+-+-+-+-+-+-+-+-+-+-+_+其别名关系式是：B=CDEQ, C=BDEQ, D=BCEQ, E=BCDQ, Q=BCDEBC=DEQ, BD=CEQ, BE=CDQ, BQ=CDE,CD=BEQ, CE=BDQ, CQ=BDE,DE=BCQ, DQ=BCE, EQ=BCD可见主效应均只与一个四因子交互作用别名。所有 的主效应都是强纯净的。类似的10个二因子交互 作用是纯净的。两个设计的优缺点：任意一个含Q的二因子交互作 用是强纯净的。弱点是其它二因子的交互作用是不 纯净的。第二个设计则平等地处理全部因子。因此 所有

27、的二因子交互作用是纯净的，但不是强纯净 的。因此若先验知识很强地表明了含Q的二因子交 互作用可能更重要，则应选择第一个设计，否则选 第二个。6响应曲面法：6.1用途：若一个过程或系统的研究常常集中在响应和输入因子 之间的关系上，其目标时优化响应或理解深层的机制。输入因子时定量的且只有少数几个因子，则响应曲面 法是研究这种关系的一个有效工具。6.2基础知识：立方体点(cube point)或角点(corner point)：位于立方体的角上的点。轴点(axial point)或星点(star point) 在坐标轴上的点。6.2.3中心点(center point)在设计区域的中心(即 原点)中

28、心复合设计(central composite design)：nf个立方体点(或角点)，它们组成设计的析 因部分。6.3.2nc个中心点6.3.32k个星点或轴点。6.3.4中心复合设计的三个问题。6.3.4.1选择设计的析因部分6.3.4.2确定星点或轴点的值6.3.4.3确定中心点的个数。6.3.5选择设计的析因部分6.3.6 星点或轴点的值一般在1和K (输入因子)开根号之间。6.3.6.2如果星点或轴点取值为1。则星点位于立 方体各个面的中心。因此成为此设计为中心立方体(face center cube)设计。这样的选择 有两个好处。1。它们是唯一的要求三水平的 中心复合设计。若有一个或多个因子是定性的 且只能取三水平，则5水平的中心符合设计就 不适用。2。若设计区域是立方体，则这种设 计是有效的。6.3.6.3可旋转的(rotatable)取值法。星点或 轴点取值为立方体点数开4次方。要求如果为不完全析因设计，要求分辨度要较高。6.3.6.4星点或轴点取值依赖于设计区域的几何6.3.6.4星点或轴点取值依赖于设计区域的几何性质和实际约束。若想收集靠近立方体面的信 息则应选择较小的值。若主要考虑的是估计效 率，则应选择较大的值。6.3.7中心点的试验次数