1、第四章第四章 正交试验设计与数据处理正交试验设计与数据处理 在生产实践中，试制新产品、改革工艺、寻求好的生产条件等，在生产实践中，试制新产品、改革工艺、寻求好的生产条件等，这些都需要做试验，而试验总是要花费时间，消耗人力、物力。因这些都需要做试验，而试验总是要花费时间，消耗人力、物力。因此，试验的次数应尽可能少。此，试验的次数应尽可能少。 全面试验：全面试验： 如如 4 个个 3 水平的因素，要做水平的因素，要做 3481 次试验；次试验； 6 个个 5 水平的因素，要做水平的因素，要做 5615625次试验。非常困难。次试验。非常困难。 能否减少试验次数，而又不影响试验效果呢？能否减少试验次
2、数，而又不影响试验效果呢？有有4.1 正交表及其用法正交表及其用法 正交表的记号：正交表的记号：L9（34）表示表示 4 个因素，每个因素取个因素，每个因素取 3 个个水平的正交表。格式如水平的正交表。格式如表表41所示。所示。2水平正交设计的方差分析水平正交设计的方差分析4.1 正交表及其用法正交表及其用法 正交表记为正交表记为 Ln（mk），），m 是各因素的水平，是各因素的水平，k （列数）是因列数）是因素的个数，素的个数，n 是安排试验的次数（行数）。是安排试验的次数（行数）。 L9（34）4因素因素3水平正交试验，共做水平正交试验，共做9次试验，面全面试验要次试验，面全面试验要做做
3、34=81 次，减少了次，减少了72次。次。 L25（56） 6因素因素5水平正交试验，共做水平正交试验，共做25次试验，面全面试验次试验，面全面试验要做要做 56=15625 次，减少了次，减少了15600次。次。 正交表的两条重要性质：正交表的两条重要性质： （1）每列中不同数字出现的次数是相等的，如）每列中不同数字出现的次数是相等的，如 L9（34），），每每列中不同的数字是列中不同的数字是1，2，3。它们各出现三次。它们各出现三次。 （2）在任意两列中，将同一行的两个数字看成有序数对时，）在任意两列中，将同一行的两个数字看成有序数对时，每种数对出现的次数是相等的，如如每种数对出现的次数
4、是相等的，如如 L9（34），），有序数对共有有序数对共有9个：个：（1，1），（），（1，2），（），（1，3），（），（2，1），（），（2，2），（），（2，3），），（3，1），（），（3，2），（），（3，3），），它们各出现一次它们各出现一次。4.1 正交表及其用法正交表及其用法 由于正交表的性质，用它来安排试验时，各因素的各种水平由于正交表的性质，用它来安排试验时，各因素的各种水平是搭配均衡的。是搭配均衡的。 下面通过具体例子来说明如何用正交表进行试验设计。下面通过具体例子来说明如何用正交表进行试验设计。 例例1 某炼铁厂为了提高铁水温度，需要通过试验选择最好的某炼铁厂为了提高铁
5、水温度，需要通过试验选择最好的生产方案，经初步分析，主要有生产方案，经初步分析，主要有 3 个因素影响铁水温度，它们是个因素影响铁水温度，它们是焦比、风压和底焦高度，每个因素都考虑焦比、风压和底焦高度，每个因素都考虑 3 个水平，具体情况如个水平，具体情况如表表42，问对这，问对这 3 个因素的个因素的 3 个水平如何安排，才能获得最高的个水平如何安排，才能获得最高的铁水温度？铁水温度？ 解：解： 在这个问题中，人们关心的是铁水温度，我们称它为试在这个问题中，人们关心的是铁水温度，我们称它为试验指标，如何安排试验才能获得最高的铁水温度，这只有通过试验验指标，如何安排试验才能获得最高的铁水温度，
6、这只有通过试验才能解决，这里有才能解决，这里有3个因素，每个因素有个因素，每个因素有3个水平，是一个个水平，是一个3因素、因素、3水平的问题，如果每个因素的每个水平都互相搭配着进行全面试验，水平的问题，如果每个因素的每个水平都互相搭配着进行全面试验，必须做试验必须做试验33=27次，将这种搭配试验编号列在次，将这种搭配试验编号列在表表43中。中。例例1 进行进行27次试验次试验 要花很多时间，耗费不少人力、物力，为了减少要花很多时间，耗费不少人力、物力，为了减少试验次数，但又不能影响试验的效果，因此，不能随便地减少试验，试验次数，但又不能影响试验的效果，因此，不能随便地减少试验，应当把有代表性
7、的搭配保留下来，为此，应当把有代表性的搭配保留下来，为此，按按 L9（34）表中前表中前3列列的的情况从情况从27个试验中选取个试验中选取9个，个，它们的序号分别为它们的序号分别为1，5，9，11，15，16，21，22，26，将这，将这9个试验按新的编号个试验按新的编号19写出来，正好是正写出来，正好是正交表交表 L9（34）的前的前3列，如列，如表表44所示。所示。 为了便于分析计算，把考查指标（铁水温度）列于表为了便于分析计算，把考查指标（铁水温度）列于表44的右的右边，做成一个新的边，做成一个新的表表45，利用张表进行分析计算。，利用张表进行分析计算。 从表从表45中的数据处理与分析，
8、可以得出结论：各因素对考查中的数据处理与分析，可以得出结论：各因素对考查指标（铁水温度）的影响按大小次序来说应当是指标（铁水温度）的影响按大小次序来说应当是A（焦比）、焦比）、C（底焦高度）、底焦高度）、B（风压），最好的方案应当是风压），最好的方案应当是A3C2B2，即：即：例例1 C2：底焦高度，第底焦高度，第2水平，水平，1.5； A3：焦比，焦比， 第第3水平，水平，1：14； B2：风压，风压， 第第2水平，水平，230。 得出的最好方案在已经做过的得出的最好方案在已经做过的9次试验中没有出现，与它比较次试验中没有出现，与它比较接近的是第接近的是第9号试验，在第号试验，在第9号试验中
9、只有风压号试验中只有风压B不是处于最好水平，不是处于最好水平，而且风压对铁水温度的影响是而且风压对铁水温度的影响是3个因素中最小的。从实际做出的结个因素中最小的。从实际做出的结果看出果看出第第9号试验中的铁水温度是号试验中的铁水温度是14100C，是是9次试验中最高的次试验中最高的，这也说明我们找出的最好方案是符合实际的。这也说明我们找出的最好方案是符合实际的。 为了最终确定试验方案为了最终确定试验方案C2A3B2是不是最好方案，可以按这个是不是最好方案，可以按这个方案再试验一次，若比方案再试验一次，若比9号好，作为最好结果，若比号好，作为最好结果，若比9号差，则以号差，则以9号为最佳条件。如
10、出现后一结果，说明我们的理论分析与实践有号为最佳条件。如出现后一结果，说明我们的理论分析与实践有一定的差距，最终还是要接受实践的检验。一定的差距，最终还是要接受实践的检验。正交试验步骤归纳如下：正交试验步骤归纳如下：1、确定要考核的试验指标；、确定要考核的试验指标；2、确定要考察的因素和各因素的水平；、确定要考察的因素和各因素的水平； 以上两条要实践经验来决定。以上两条要实践经验来决定。3、选用合适的正交表，一般只要正交表中的因素个数比试验要考、选用合适的正交表，一般只要正交表中的因素个数比试验要考察的因素的个数稍大或相等就行了。这样既保证了试验目的，而察的因素的个数稍大或相等就行了。这样既保
11、证了试验目的，而试验次数又不致太多，省工省时；试验次数又不致太多，省工省时；4、试验，测定试验指标；、试验，测定试验指标；5、试验结果分析计算，得出合理的结论。、试验结果分析计算，得出合理的结论。 以上的方法以上的方法直观分析法。简单、计算量小、很实用。直观分析法。简单、计算量小、很实用。 正交试验的主要分析工具是正交表，而在因素及其水平都确正交试验的主要分析工具是正交表，而在因素及其水平都确定的情况下，正交表并不是唯一的，常见的正交表见本书末定的情况下，正交表并不是唯一的，常见的正交表见本书末P214的附表的附表6。4.2 多指标的分析方法多指标的分析方法习题习题 4.1 在例在例1中，试验
12、指标只有一个，考察起来比较方便，但实际问中，试验指标只有一个，考察起来比较方便，但实际问题中，需要考察的指标往往不止一个，有时有两个、三个或更多。题中，需要考察的指标往往不止一个，有时有两个、三个或更多。如何评价考察指标呢？如何评价考察指标呢？两种方法两种方法。一、综合平衡法一、综合平衡法 通过具体的例子来加以说明。通过具体的例子来加以说明。 例例2 为了提高产品质量，要对生产的原料进行配方试验。要为了提高产品质量，要对生产的原料进行配方试验。要检验检验3项指标：抗压强度、落下强度和裂纹度，前两个指标越大越项指标：抗压强度、落下强度和裂纹度，前两个指标越大越好，第好，第3个指标越小越好。根据以
13、往的经验，配方有个指标越小越好。根据以往的经验，配方有3个重要因素：个重要因素：水分、粒度和碱度。它们各有水分、粒度和碱度。它们各有3个水平，具体数据如个水平，具体数据如表表46所示。所示。试进行试验分析，找出最好的配方方案。试进行试验分析，找出最好的配方方案。第四章第四章 习题习题4.1 为提高烧结矿的质量，做下面的配料试验。各因素水平如下为提高烧结矿的质量，做下面的配料试验。各因素水平如下 表（单位：表（单位：t）：）：因素因素水平水平A精矿精矿B生矿生矿C焦粉焦粉D石灰石灰E白云石白云石F铁屑铁屑18.05.00.82.01.00.529.54.00.93.00.51.0 反映质量好坏的
14、试验指标为含铁量，越高越好。从附表反映质量好坏的试验指标为含铁量，越高越好。从附表6中选中选择择L8（27）的正交表安排试验。各因素依次放在正交表的的正交表安排试验。各因素依次放在正交表的16列上，列上，8次试验所得含铁量（次试验所得含铁量（%）依次为：）依次为：50.9, 47.1, 51.4, 51.8, 54.3, 49.8, 51.5, 51.3。一、综合平衡法（例一、综合平衡法（例2的解）的解） 解解 3因素因素3水平，应选水平，应选L9（34）正交表来安排试验，将正交表来安排试验，将3个因个因素依次放在前素依次放在前3列（第列（第4列不要），得出一张具体的试验方案表，列不要），得出
15、一张具体的试验方案表，测出需要检验的指标结果，列于测出需要检验的指标结果，列于表表47中，然后用直观分析法对中，然后用直观分析法对每个指标分别进行计算分析。每个指标分别进行计算分析。 将将3 个指标分别进行计算分析后，得出个指标分别进行计算分析后，得出3个好的方案：对抗压个好的方案：对抗压强度是强度是A2B3C1；对落下强度是对落下强度是A3B3C2；对裂纹度是对裂纹度是A2B3C1，这这3个个方案不完全相同，对一个指标是好方案，而对另一个指标却不一方案不完全相同，对一个指标是好方案，而对另一个指标却不一定是好方案，如何找出对各个指标都较好的一个共同方案呢？定是好方案，如何找出对各个指标都较好
16、的一个共同方案呢？ 综合分析，将指标随因素水平变化的情况用图形表示出来，综合分析，将指标随因素水平变化的情况用图形表示出来，如如图图41所示（为了看得清楚，将各点用直线连接起来，实际上所示（为了看得清楚，将各点用直线连接起来，实际上并不一定是直线。并不一定是直线。 把图把图41和表和表47结合起来分析，看每一个因素对各指标的结合起来分析，看每一个因素对各指标的影响。影响。图图4110.19.011.29.26.814.312.78.39.38.63.15.91.13.313.36.98.32.42.03.71.73.02.71.71.72.73.00.02.04.06.08.010.012.0
17、14.016.07894681.11.31.5水分水分 粒度粒度 碱度碱度-抗压强度抗压强度落下强度落下强度裂纹度裂纹度一、综合平衡法（例一、综合平衡法（例2的解的综合分析）的解的综合分析）（1）粒度）粒度B对抗压强度和落下强度来讲，对抗压强度和落下强度来讲，极差最大极差最大，是最大的影，是最大的影响因素。从响因素。从图图41中看出三个指标中看出三个指标B均取均取8为最好为最好即取即取B3。（2）碱度）碱度C，极差不大极差不大，次要因素。由，次要因素。由图图41分析，取分析，取1.1时两个时两个指标好，指标好，1个指标稍差，对三个指标综合考虑，个指标稍差，对三个指标综合考虑，C取取1.1即取即
18、取C1。（3）水分）水分A，对对裂纹度影响极差最大裂纹度影响极差最大，A取取9最好，由最好，由图图41综合综合考虑考虑A取取9即取即取A2。 通过各因素对各指标影响的综合分析，得出较好的试验方案通过各因素对各指标影响的综合分析，得出较好的试验方案是：是： B3：粒度取第粒度取第3水平，水平，8； C1：碱度取第碱度取第1水平，水平，1.1； A2：水分取第水分取第2水平，水平，9。二、综合评分法二、综合评分法 对多指标的问题，真正做到好的综合平衡，有时很困难，这对多指标的问题，真正做到好的综合平衡，有时很困难，这是综合平衡法的缺点。综合评分法可以克服这个缺点。是综合平衡法的缺点。综合评分法可以
19、克服这个缺点。 例例3 某厂生产一种化工产品，需要检验两个指标：核酸纯度某厂生产一种化工产品，需要检验两个指标：核酸纯度和回收率，这两个指标都是越大越好。有影响的因素有和回收率，这两个指标都是越大越好。有影响的因素有4个，各有个，各有3个水平，具体情况如个水平，具体情况如表表48所示。试通过试验分析出较好方案，所示。试通过试验分析出较好方案，使产品的核酸含量和回收率都有提高。使产品的核酸含量和回收率都有提高。 解解 这是这是4因素因素3水平的试验，可以选用正交表水平的试验，可以选用正交表L9（34）安排出安排出试验方案（这里有试验方案（这里有4个因素，正好将表排满），进行试验，将得出个因素，正
20、好将表排满），进行试验，将得出的结果列入的结果列入表表49中。中。 综合评分法是根据各个指标的重要性的不同，按照得出的试验综合评分法是根据各个指标的重要性的不同，按照得出的试验结果综合分析，给每一个试验评出一个分数，作为这个试验的总指结果综合分析，给每一个试验评出一个分数，作为这个试验的总指标。根据这个总指标作进一步的分析。标。根据这个总指标作进一步的分析。二、综合评分法（例二、综合评分法（例3的解）的解） 这个方法的关键是如何评分。这个方法的关键是如何评分。 在这个试验中，两个指标的重要性是不同的，根据实践经验知在这个试验中，两个指标的重要性是不同的，根据实践经验知道，纯度的重要性大于回收率
21、，从实际分析，可以认为纯度是回收道，纯度的重要性大于回收率，从实际分析，可以认为纯度是回收率的率的4倍。也就是倍。也就是纯度占权数为纯度占权数为4，回收率占权数为，回收率占权数为1，按这个权数，按这个权数给出这个试验的总分为：给出这个试验的总分为：总分总分4纯度纯度1回收率回收率 由上式计算出这个试验的总分数，列于由上式计算出这个试验的总分数，列于表表49的最右边，再的最右边，再根据这个分数，用直观分析法进行分析。根据这个分数，用直观分析法进行分析。 从表从表49看出，看出，A、D两个因素的极差都很大，是对试验影响两个因素的极差都很大，是对试验影响很大的两个因素，很大的两个因素，A1、D1为好
22、；为好；B因素的极差比因素的极差比A、D的极差小，的极差小，对试验的影响比对试验的影响比A、D都小；都小；B因素取因素取B3为好；为好；C因素的极差最小，因素的极差最小，影响最小，影响最小，C取取C2为好。综合考虑，最好的试验方案应当是为好。综合考虑，最好的试验方案应当是A1B3C2D1，按影响大小次序排列为：按影响大小次序排列为：二、综合评分法（例二、综合评分法（例3的解）的解）A1：时间，时间， 25小时；小时；D1：加水量，加水量， 1：6；B3：料中核酸含量，料中核酸含量，6.0；C2：pH值，值， 6.0。 可以看出，这里分析出来的最好方案，在已经做过的可以看出，这里分析出来的最好方
23、案，在已经做过的9个试验个试验中是没有的，可以按这个方案再试验一次，中是没有的，可以按这个方案再试验一次，看能不能得出比第看能不能得出比第1号号试验更好的结果试验更好的结果，从而确定出真正最好的试验方案。，从而确定出真正最好的试验方案。 综合评分法是将多指标的问题，通过加权计算总分的方法化综合评分法是将多指标的问题，通过加权计算总分的方法化成一个指标的问题，这样对结果的分析计算都比较方便、简单。成一个指标的问题，这样对结果的分析计算都比较方便、简单。但如何合理地评分，是最关键的问题。这一点只能依据实际经验但如何合理地评分，是最关键的问题。这一点只能依据实际经验来解决，单纯从数学上是无法解决的。
24、来解决，单纯从数学上是无法解决的。习题习题42习题习题42 某厂生产液体葡萄糖，要对生产工艺进行优选试验。因素水某厂生产液体葡萄糖，要对生产工艺进行优选试验。因素水平如下表所示：平如下表所示：因因 素素 水水 平平A粉浆浓度（粉浆浓度（%）B粉浆酸度粉浆酸度C稳压时间稳压时间D工作压力（工作压力（105Pa）1161.502.22182.052.73202.5103.2 试验指标有两个试验指标有两个：(1)产量，越高越好；（产量，越高越好；（2）总还原糖，在）总还原糖，在3240%之间。用正交表之间。用正交表L9（34）安排试验，安排试验，9次试验所得结果如下：次试验所得结果如下：产量产量(k
25、g): 498，568，568，577，512，540，501，550，510；还原糖还原糖(%)：41.6, 39.4, 31.0, 42.4, 37.2, 30.2, 42.2, 40.4, 30.0。试用综合平衡法对结果进行分析。找出最好的生产方案。试用综合平衡法对结果进行分析。找出最好的生产方案。4.3 混合水平的正交试验设计混合水平的正交试验设计 在实际情况中，有时做试验时，每个因素的水平数是不同的在实际情况中，有时做试验时，每个因素的水平数是不同的混合水平混合水平。两种解决方案。两种解决方案。一、混合水平正交表及其用法一、混合水平正交表及其用法 混合水平正交表就是各因素的水平数不完
26、全相等的正交表。混合水平正交表就是各因素的水平数不完全相等的正交表。这种正交表有好多种。比如这种正交表有好多种。比如 L8（4124）就是一个混合水平的正交就是一个混合水平的正交表，如表，如表表410所示。所示。 其它混合水平的正交表还有很多，见附表所示，它们都有上其它混合水平的正交表还有很多，见附表所示，它们都有上面所说的两点。面所说的两点。 例例4 某农科站进行品种试验，具体试验因素及水平如某农科站进行品种试验，具体试验因素及水平如表表411所示。试验指标是产量，数值越大越好。试用混合正交表所示。试验指标是产量，数值越大越好。试用混合正交表安排试验，找出最好的试验方案。安排试验，找出最好的
27、试验方案。例例 4 的的 解解解解 这个问题中有这个问题中有4个因素，个因素，1个是个是4水平的，水平的，3个是个是2水平的，正好水平的，正好可以选用混合正交表可以选用混合正交表 L8（4123），），因素因素A为为4水平，放在第水平，放在第1列，列，其余其余3个个2水平的因素水平的因素B、C、D顺序放在顺序放在2、3、4列上，第列上，第5列不用。列不用。按这个方案进行试验，将得出的试验结果放在正交表的右边，然按这个方案进行试验，将得出的试验结果放在正交表的右边，然后进行分析，见后进行分析，见表表412。 经分析得最佳方案为：经分析得最佳方案为： A2B2C2D2。因为，从极差分析可知，因为，
28、从极差分析可知，因素因素D影响很小，这个方案影响很小，这个方案与第与第4号试验结果号试验结果A2B2C2D1很接近很接近，从，从试验试验 结果看出，第结果看出，第4号试验是号试验是8个试验中产量最高的，因此完全有个试验中产量最高的，因此完全有理由取第理由取第4号试验作为最好的试验方案加以推广。号试验作为最好的试验方案加以推广。二、拟水平法二、拟水平法例例5 今有某一试验，试验指标只有一个，它的数值越小越好，这今有某一试验，试验指标只有一个，它的数值越小越好，这个试验有个试验有4个因素个因素A、B、C、D，其中因素其中因素C是是2水平的，其余水平的，其余3个个因素都是因素都是3水平的，具体数值如
29、水平的，具体数值如表表413所示。试安排试验，并对所示。试安排试验，并对试验结果进行分析，找出最好的试验方案。试验结果进行分析，找出最好的试验方案。 解解 ： 4因素试验，因素试验，C为为2个水平，个水平，A、B、D为为3个水平。没有个水平。没有合适的正交表。设想：假若合适的正交表。设想：假若C有有3个水平，就变成个水平，就变成4因素因素3水平的问水平的问题了。如何将题了。如何将C变成变成3水平的因素呢？从水平的因素呢？从C中的中的1和和2水平中选一个水平中选一个水平让它重复一次作为第水平让它重复一次作为第3水平，这就叫虚拟水平。取哪一个水水平，这就叫虚拟水平。取哪一个水平作为第平作为第3水平
30、呢？一般说，都是要根据实际经验，选取一个较水平呢？一般说，都是要根据实际经验，选取一个较好的水平。比如，如果认为第好的水平。比如，如果认为第2水平比第水平比第1水平好，就选第水平好，就选第2水平水平作为第作为第3水平。这样因素水平水平。这样因素水平表表413就变为表就变为表414的样子，它的样子，它比比表表413多了一个虚拟的第多了一个虚拟的第3水平。水平。例例 5 的的 解解 这样就变成了一个这样就变成了一个4因素因素3水平的试验，可以按水平的试验，可以按 L9（34）表安表安排试验，并对正交表进行重构，测出结果，并进行分析，见排试验，并对正交表进行重构，测出结果，并进行分析，见表表415所
31、示。所示。 从表从表415的极差可以看出，因素的极差可以看出，因素D对试验的影响最大，取第对试验的影响最大，取第3水平最好；其次是因素水平最好；其次是因素A，取第取第3水平为好；再者是因素水平为好；再者是因素B，取第取第1水平为好；因素水平为好；因素C影响最小，取第影响最小，取第1水平为好。最优方案为：水平为好。最优方案为：A3B1D3C1。这个方案在这个方案在9个试验中没有。从试验结果看个试验中没有。从试验结果看8号试验为号试验为最好。这个试验只有最好。这个试验只有B不是处在最好情况，而因素不是处在最好情况，而因素B的影响是最小的影响是最小的。可以按这个方案再试验一次，看是否会得出比第的。可
32、以按这个方案再试验一次，看是否会得出比第8号试验更好号试验更好的结果，从而确定出真正的最优方案。的结果，从而确定出真正的最优方案。4.4 有交互作用的正交试验设计有交互作用的正交试验设计.:.:;BABA记记为为因因素素的的交交互互作作用用各各因因素素联联合合起起作作用用因因素素如如各各因因素素各各自自独独立立起起作作用用多多因因素素试试验验例例6 有有4块试验田，土质情况基本一样，种植同样的作物。现将块试验田，土质情况基本一样，种植同样的作物。现将氮肥、磷肥采用不同的方式分别加在氮肥、磷肥采用不同的方式分别加在4块地里，收获后算出平均亩块地里，收获后算出平均亩产，记在产，记在表表416中。中
33、。氮肥、磷肥交互作用的效果氮肥、磷肥的总效果（只加氮肥、磷肥交互作用的效果氮肥、磷肥的总效果（只加氮肥的效果只加磷肥的效果）氮肥的效果只加磷肥的效果） 在多因素的试验中，交互作用影响的大小参照实际经验。如在多因素的试验中，交互作用影响的大小参照实际经验。如果确有把握认定交互作用的影响很小，就可以忽略不计；如果不果确有把握认定交互作用的影响很小，就可以忽略不计；如果不能确认交互作用的影响很小，就应该通过试验分析交互作用的大能确认交互作用的影响很小，就应该通过试验分析交互作用的大小。小。一、交互作用表一、交互作用表 下面介绍交互作用表和它的用法，下面介绍交互作用表和它的用法，表表417就是正交表就
34、是正交表 L8（27）所对应的交互作用表。所对应的交互作用表。 P214列出了几个交互作用的正交表。列出了几个交互作用的正交表。正交表自由度的确定：正交表自由度的确定：（1）每列的自由度）每列的自由度 f列列水平数水平数1（2）两因素交互作用的自由度）两因素交互作用的自由度 fABfAfB （两因素自由度的乘积）两因素自由度的乘积）对对2因素因素2水平的正交表，因为：水平的正交表，因为：fAfB 21，每列只有一每列只有一个自由度；而个自由度；而fABfAfB 111，所以也占一列。，所以也占一列。一、交互作用表一、交互作用表 对于对于2 因素因素3水平，水平， fAfB 312，每列有每列有
35、2个自由度；个自由度；而而fABfAfB 224，由于交互作用列有，由于交互作用列有4个自由度，而每个自由度，而每列是列是2个自由度，因此个自由度，因此2个个3水平因素的交互作用列占水平因素的交互作用列占2列。列。 对于对于2因素因素n水平，水平， fAfB n1，每列有每列有n个自由度；个自由度；而两因素交互作用的自由度为：而两因素交互作用的自由度为：fABfAfB （n1）（）（n1），），所以交互作用列要占（所以交互作用列要占（n1）列。列。二、水平数相同的有交互作用的正交设计二、水平数相同的有交互作用的正交设计 例例7 某产品的产量取决于某产品的产量取决于3个因素个因素A、B、C，每个
36、因素都有每个因素都有两个水平，具体数值如两个水平，具体数值如表表418所示。每两个因素都有交互作用，所示。每两个因素都有交互作用，必须考虑。试验指标为产量，越高越好。试安排试验，并分析试必须考虑。试验指标为产量，越高越好。试安排试验，并分析试验结果，找出最好的方案。验结果，找出最好的方案。例例7的解的解 解解 这是这是3因素因素2水平的试验。水平的试验。3个因素个因素A、B、C要占要占3列，它列，它们之间的交互作用们之间的交互作用AB、B C、A C又各占又各占3列，共占列，共占6列，可列，可以用正交表以用正交表 L8（27）来安排试验。若将来安排试验。若将A、B放在第放在第1、2列，从列，从
37、表表417查出查出AB应在第应在第3列，因此列，因此C就不能放在第就不能放在第3列，否则就要和列，否则就要和AB混杂。现将混杂。现将C放在第放在第4列，由列，由表表417查出查出A C应在第应在第5列，列， B C应在第应在第6列。按这种安排进行试验。测出结果，用直观分析列。按这种安排进行试验。测出结果，用直观分析法进行分析，把交互作用当成新的因素看待。整个分析过程记录法进行分析，把交互作用当成新的因素看待。整个分析过程记录在在表表419中。中。 最后要说明一点，在这里只考虑两列间的交互作用最后要说明一点，在这里只考虑两列间的交互作用AB、B C、A C，3个因素的交互作用个因素的交互作用AB
38、C，一般影响很小，这一般影响很小，这里不去考虑它。里不去考虑它。习题习题 43 在梳棉机上纺粘棉混纱，为了提高质量，选了在梳棉机上纺粘棉混纱，为了提高质量，选了3个因素，每个个因素，每个因素有两个水平，因素有两个水平，3因素之间有一级交互作用。因素水平如下：因素之间有一级交互作用。因素水平如下：因素因素水平水平A金属针布金属针布B产量水平（产量水平（kg）C速度（速度（r/min）1甲地产品甲地产品62382乙地产品乙地产品10320 试验指标为棉结粒数，越小越好。用试验指标为棉结粒数，越小越好。用L8（27）安排试验，安排试验，8次次试验所得试验指标的结果依次为：试验所得试验指标的结果依次为
39、：0.3，0.55, 0.40, 0.30, 0.15, 0.40, 0.50, 0.35。试对结果进行分析，选出最佳工艺条件。试对结果进行分析，选出最佳工艺条件。4.5 正交表的构造法正交表的构造法 从前面的内容可以看出，正交表的用处和好处。那么正交表从前面的内容可以看出，正交表的用处和好处。那么正交表是如何得来的呢？下面就介绍两种正交表的构造方法。是如何得来的呢？下面就介绍两种正交表的构造方法。一、阿达玛矩阵法一、阿达玛矩阵法1、阿达玛矩阵、阿达玛矩阵 阿阵阿阵 定义：以定义：以1，1为元素，并且任意两列都是正交的矩为元素，并且任意两列都是正交的矩阵。阵。性质：性质： （1）每列元素个数都
40、）每列元素个数都 是偶数；是偶数； （2）任意两列（两行）交换后，仍为阿阵；）任意两列（两行）交换后，仍为阿阵；（3）任意一列（或行）乘）任意一列（或行）乘1以后，仍为阿阵。以后，仍为阿阵。 标准阿阵：第一列全为标准阿阵：第一列全为1列（用对行乘列（用对行乘1可得）。可得）。阿方阵：行、列相等阿方阵：行、列相等阿阵，偶阶方阵。阿阵，偶阶方阵。1、阿达玛矩阵、阿达玛矩阵 n 阶阿阵记为阶阿阵记为Hn。感兴趣：感兴趣：第一列，第一行全为第一列，第一行全为 1 的阿阵。例如：的阿阵。例如：1111111111111111,111142HH直积构造高阶阿阵的方法：直积构造高阶阿阵的方法：定义：设两个定
41、义：设两个2阶方阵阶方阵A、B2221121122211211,bbbbBaaaaA它们直积记为它们直积记为A B，定义如下：定义如下：BaBaBaBabbbbaaaaBA2221121122211211222112111、阿达玛矩阵、阿达玛矩阵2222212222212121122211221221112122122112221121111212111212111111babababababababababababababababa这是一个这是一个4阶方阵。阶方阵。有下面两个定理：有下面两个定理： 定理一定理一 设设2阶方阵阶方阵A、B如果它们中的两列是正交的，则它如果它们中的两列是正交的，
42、则它们的直积们的直积A B的任意两列也是正交的。的任意两列也是正交的。 定理二定理二 两个阿阵的直积是一个高阶阿阵两个阿阵的直积是一个高阶阿阵。 据此，可以用简单的低阶阿阵，用求直积的方法得出高阶阿据此，可以用简单的低阶阿阵，用求直积的方法得出高阶阿阵，例如有：阵，例如有：1、阿达玛矩阵、阿达玛矩阵422111111111111111111111111HHH依此类推有：依此类推有：,16441682, 842HHHHHHHHH或或一个固定阶的阿阵并不是唯一的。比如：一个固定阶的阿阵并不是唯一的。比如：1111,1111,1111,1111都是都是2 阶阿阵阶阿阵H2，但我们最感兴趣的是第一个但
43、我们最感兴趣的是第一个。2、2个水平正交表的阿达玛矩阵法个水平正交表的阿达玛矩阵法 有了第有了第1列第列第1行全为行全为1的标准阿阵，构造的标准阿阵，构造2水平的正交表就非水平的正交表就非常方便了。常方便了。（1）L4（23）正交表的构造正交表的构造 取标准阿阵取标准阿阵H4 如下：如下： 11111111111111114H 将全将全1列去掉，得出：列去掉，得出： 1111111111112、2个水平正交表的阿达玛矩阵法个水平正交表的阿达玛矩阵法 将将1 改写为改写为2，按顺序配上列号、行号，就得到，按顺序配上列号、行号，就得到2水平正交表水平正交表 L4（23），），见见表表428所示。所
44、示。（2） L8（27）正交表的构造法正交表的构造法 取标准阿阵取标准阿阵 H8 如下：如下： 1111111111111111111111111111111111111111111111111111111111111118H2、2个水平正交表的阿达玛矩阵法个水平正交表的阿达玛矩阵法 去掉全去掉全 1 列；列； 将将 1 改写为改写为2，并按顺序配上列号、行号，就得到正交表，并按顺序配上列号、行号，就得到正交表 L8（27），），见见表表429。总结：先取一个标准阿阵总结：先取一个标准阿阵 Hn ，去掉全去掉全1列，将列，将1 列改写为列改写为2， 配上列号、行号，就得正交表配上列号、行号，就
45、得正交表 Ln（2n1）。）。 上法只能构造上法只能构造 2 水平正交表，更多水平的正交表，用正交拉水平正交表，更多水平的正交表，用正交拉丁方的方法来解决。丁方的方法来解决。二、正交拉丁方的方法二、正交拉丁方的方法1、拉丁方、拉丁方定义：用定义：用 n 个不同的拉丁字母排成一个个不同的拉丁字母排成一个 n 阶方阵（阶方阵（n26），），如如 果每个字母在任一行、任一列中只出现一次，则称这种方果每个字母在任一行、任一列中只出现一次，则称这种方 阵为阵为 nn 拉丁方，简称为拉丁方，简称为 n 阶拉丁方。阶拉丁方。例如，用例如，用3个字母个字母 A、B、C 可排成：可排成：BACACBCBA33
46、拉丁方拉丁方用用4个字母个字母 A、B、C、D 可排成：可排成：ADCBBADCCBADDCBA44 拉丁方拉丁方这这两两个个拉拉丁丁方方不不是是唯唯一一的。的。1、拉丁方、拉丁方正交拉丁方。正交拉丁方。定义：设有两个同阶的拉丁方，如果对第一个拉丁方排列着相同定义：设有两个同阶的拉丁方，如果对第一个拉丁方排列着相同字母的各个位置上，第二拉丁方在同样位置上排列着不同字母，字母的各个位置上，第二拉丁方在同样位置上排列着不同字母，则称这两个拉丁方为互相正交的拉丁方。则称这两个拉丁方为互相正交的拉丁方。BACACBCBAACBBACCBA3阶阶拉拉丁丁方方与与是正交拉丁方。正交拉丁方只有两个。是正交拉
47、丁方。正交拉丁方只有两个。ABCDBADCCDABDCBABADCCDABABCDDCBA与与 4阶拉丁方中，正交拉丁方只有阶拉丁方中，正交拉丁方只有3个；个； 5阶拉丁方中，正交拉丁方只有阶拉丁方中，正交拉丁方只有4个；个； 6阶拉丁方中，正交拉丁方只有阶拉丁方中，正交拉丁方只有5个；个； 数学上已经证明：数学上已经证明：n 阶拉丁方的正交拉丁方阶拉丁方的正交拉丁方个数为：（个数为：（n1）个。个。1、拉丁方、拉丁方1、拉丁方、拉丁方 将字母拉丁方改写为数字拉丁方性质没有影响。比如将字母拉丁方改写为数字拉丁方性质没有影响。比如 3 阶拉阶拉丁方可写为：丁方可写为：21313232113221
48、3321与与为正交拉丁方。为正交拉丁方。其它的正交拉丁方见其它的正交拉丁方见P236附表附表8。2、3水平正交表的构造水平正交表的构造 首先考虑两个首先考虑两个 3 水平因素水平因素A、B，把它们所有水平搭配都写出把它们所有水平搭配都写出来来 32 9个，个，按下面的方式排成两列按下面的方式排成两列：3、4水平正交表水平正交表 因素因素A、B两个两个4水平的全排列水平的全排列 4216个，构成基本列；个，构成基本列； 三个正交拉丁方，按三个正交拉丁方，按1，2，3，4列分别按顺序排成列分别按顺序排成1列，共列，共3列，放在基本列右则，得列，放在基本列右则，得5列列16行矩阵。行矩阵。 得表得表
49、431，为，为L16（45）正交表正交表。四、混合型正交表的构造法四、混合型正交表的构造法 混合型正交表可以由一般水平数相等的正交表通过混合型正交表可以由一般水平数相等的正交表通过“并列法并列法”改造而成。举典型的例子加以说明。改造而成。举典型的例子加以说明。例例8 混合型正交表混合型正交表 L8（424）的构造法。的构造法。解：解：（1）先列出正交表）先列出正交表 L8（424），），见见表表432。（2）取出）取出表表432中的中的1，2列，将数对（列，将数对（1，1）、）、 （1，2）、）、 （2，1）、（）、（2，）与单数字，）与单数字1，2，3，4依次对应，作为依次对应，作为 新表第
50、新表第1列。列。（3）去掉去掉12的第的第3列列。交互作用。交互作用。（4）4，5，6，7列左移，依次变为新表的列左移，依次变为新表的2，3，4，5列列。列号列号行号行号11111111211122223122112241222211521212126212212172211221822121124567123表表432 L8（27）正交表正交表列号列号行号行号11111121222232112242221152121262212171122181211234567列号列号行号行号111111121122223221122422221153212126322121741122184121121
