xmccy
这个根据因素数，因素多的话，水平不宜多，最好不超过5个。还要考虑交互项的影响。一般就3到5个吧。 
正交试验法的评价正交试验法的理论基础是正交拉丁方理论与群论。在工作中可用的多因素寻优工作方法，一类是从优选区某一点开始试验，一步一步到达较优点，这类实验方法叫序贯试验法，如因素轮换法、爬山法等；另一类是，在优选区内一次布置一批试验点，通过对这批试验结果的分析，逐步缩小优选范围从而达到较优点，如正交试验法等。科研中普遍采用正交试验法，因其具有如下优点：①实用上按表格安排试验，使用方便； ②布点均衡、试验次数较少； ③在正交试验法中的最好点，虽然不一定是全面试验的最好点，但也往往是相当好的点。特别在只有一两个因素起主要作用时，正交试验法能保证主要因素的各种可能都不会漏掉。这点在探索性工作中很重要，其他试验方法难于作到； ④正交试验法提供一种分析结果（包括交互作用）的方法，结果直观易分析。且每个试验水平都重复相同次数，可以消除部分试验误差的干扰； ⑤因其具有正交性，易于分析出各因素的主效应。 但其也有一些缺点：它提供的数据分析方法所获得的优选值，只能是试验所用水平的某种组合，优选结果不会超越所取水平的范围；另外，也不能给进一步的试验提供明确的指向性，使试验仍然带很强的摸索性色彩，不很精确。这样，正交试验法用在初步筛选时显得收敛速度缓慢、难于确定数据变化规律，增加试验次数。尤其在试验工作烦琐、费用昂贵的情况更显突出。正交试验法的代数学基础 对试验的寻优工作，用数学语言可描述为求多维连续空间上的最大或最小值（极值）。但现实的试验工作又往往没有可用的数学模型，不能确切知道数据变化的数学规律。故处理上可以先求出其数学模型，再计算极值；或直接从试验点的组合中推算出一个较好值作为较优解。 实际上，在高等数学上的泰勒级数：f(x)=∑f(n)(x0)·(x-x0)n/n!，n=0~+∞，用于求复杂函数的近似解时，就利用了其收敛性原理。在试验寻优时也可同理可将描述试验对象的“数学函数”运用泰勒级数的数项代数式近似地拟合试验规律，代入各次试验的结果获得一组线性方程，用解析法求出方程的优值，即曲线的极点。 显然，如果数据量大，则可以一次性用较高阶幂级数解出很精确的拟合曲线或函数。只消再用极少的进一步试验印证或寻优即可完成试验工作。缺点：代数学方法不能提供一种较好的试验安排方案，不如正交试验法规范直观，数据不易处理，故不常为人们使用。 以上分析可以看到二者均是寻优，各有所长。如何取长补短？二者联合运用原理 以L9(34)三水平标准正交表为例加以说明。L9(34)有九次试验，如对每个因素均使用二次方程拟合，因素间无交互作用即理解为各因素独立对结果作贡献。用代数方程表达有：f(x1，x2，x3，x4)=a0+a1x1+a2x12+a3x2+a4x22+a5x3+a6x32+a7x4+a8x42 代入各xi的值，获得含九个未知数的九个线性方程，可求解出各ai 。再运用多元函数求极值的方法，可以获得较优值。这样，一组试验便获得较优结果，而且不受水平取值范围的限制，并对进一步试验有较强的指向性。 显然，从正交表得到的线性方程组的系数矩阵均是满秩的，对应的线性方程组有唯一解，即试验的优值是唯一的： AX=B X=A-1B 其次，对因素间的交互作用，在代数处理上，可用因素间相乘的项来表达：xim·xjn。（显然，代数法还可用来分析多元交互作用等问题。） 以上分析说明正交试验法有其代数学基础。 很易看到，对标准正交表，各因素的基本拟合级数的最高幂次为其水平数减一。[备注]关于均匀设计 在均匀设计中，如3因素7水平的试验方案，只消7次试验，即U7（73）。这样的均匀设计的7次试验也可代数法处理求解，各因素分配二次幂，另设一常数项，但数据的正交性不易保证。而正交试验法的处理很明确易理解。觉得好的话去这个网
清访问：_asp?title=%D5%FD%BD%BB%CA%D4%D1%E9第5章 正交试验设计方法5．1 试验设计方法概述 试验设计是数理统计学的一个重要的分支。多数数理统计方法主要用于分析已经得到的数据，而试验设计却是用于决定数据收集的方法。试验设计方法主要讨论如何合理地安排试验以及试验所得的数据如何分析等。 例5-1 某化工厂想提高某化工产品的质量和产量，对工艺中三个主要因素各按三个水平进行试验（见表5-1）。试验的目的是为提高合格产品的产量，寻求最适宜的操作条件。对此实例该如何进行试验方案的设计呢？很容易想到的是全面搭配法方案（如图5-1所示）： 此方案数据点分布的均匀性极好，因素和水平的搭配十分全面，唯一的缺点是实验次数多达33＝27次（指数3代表3个因素，底数3代表每因素有3个水平）。因素、水平数愈多，则实验次数就愈多，例如，做一个6因素3水平的试验，就需36＝729次实验，显然难以做到。因此需要寻找一种合适的试验设计方法。试验设计方法常用的术语定义如下。试验指标：指作为试验研究过程的因变量，常为试验结果特征的量（如得率、纯度等）。例1的试验指标为合格产品的产量。因素：指作试验研究过程的自变量，常常是造成试验指标按某种规律发生变化的那些原因。如例1的温度、压力、碱的用量。水平：指试验中因素所处的具体状态或情况，又称为等级。如例1的温度有3个水平。温度用T表示，下标1、2、3表示因素的不同水平，分别记为T1、T2、T3。常用的试验设计方法有：正交试验设计法、均匀试验设计法、单纯形优化法、双水平单纯形优化法、回归正交设计法、序贯试验设计法等。可供选择的试验方法很多，各种试验设计方法都有其一定的特点。所面对的任务与要解决的问题不同，选择的试验设计方法也应有所不同。由于篇幅的限制，我们只讨论正交试验设计方法。5．2 正交试验设计方法的优点和特点用正交表安排多因素试验的方法，称为正交试验设计法。其特点为：①完成试验要求所需的实验次数少。②数据点的分布很均匀。③可用相应的极差分析方法、方差分析方法、回归分析方法等对试验结果进行分析，引出许多有价值的结论。从例1可看出，采用全面搭配法方案，需做27次实验。那么采用简单比较法方案又如何呢？先固定T1和p1，只改变m，观察因素m不同水平的影响，做了如图2-2（1）所示的三次实验，发现 m＝m2时的实验效果最好（好的用¬¬ □ 表示），合格产品的产量最高，因此认为在后面的实验中因素m应取m2水平。固定T1和m2，改变p的三次实验如图5-2（2）所示，发现p＝p3时的实验效果最好，因此认为因素p应取p3水平。 固定p3和m2，改变T 的三次实验如图5-2（3）所示，发现因素T 宜取T2水平。因此可以引出结论：为提高合格产品的产量，最适宜的操作条件为T2p3m2。与全面搭配法方案相比，简单比较法方案的优点是实验的次数少，只需做9次实验。但必须指出，简单比较法方案的试验结果是不可靠的。因为，①在改变m值（或p值，或T值）的三次实验中，说m2（或p3或T2 ）水平最好是有条件的。在T ≠T1，p ≠p1时，m2 水平不是最好的可能性是有的。②在改变m的三次实验中，固定T ＝T2，p ＝p3 应该说也是可以的，是随意的，故在此方案中数据点的分布的均匀性是毫无保障的。③用这种方法比较条件好坏时，只是对单个的试验数据进行数值上的简单比较，不能排除必然存在的试验数据误差的干扰。运用正交试验设计方法，不仅兼有上述两个方案的优点，而且实验次数少，数据点分布均匀，结论的可靠性较好。正交试验设计方法是用正交表来安排试验的。对于例1适用的正交表是L9（34），其试验安排见表5-2。所有的正交表与L9（34）正交表一样，都具有以下两个特点： （1） 在每一列中，各个不同的数字出现的次数相同。在表L9（34）中，每一列有三个水平，水平1、2、3都是各出现3次。 （2） 表中任意两列并列在一起形成若干个数字对，不同数字对出现的次数也都相同。在表L9（34）中，任意两列并列在一起形成的数字对共有9个：（1，1），（1，2），（1，3），（2，1），（2，2），（2，3），（3，1），（3，2），（3，3），每一个数字对各出现一次。表5－2 试验安排表试验号 列号 1 2 3 4 因素 温度℃ 压力Pa 加碱量kg 符号 T p m 123456789 1（T1）1（T1）1（T1）2（T2）2（T2）2（T2）3（T3）3（T3）3（T3） 1（p1）2（p2）3（p3）1（p1）2（p2）3（p3）1（p1）2（p2）3（p3） 1（m1）2（m2）3（m3）2（m2）3（m3）1（m1）3（m3）1（m1）2（m2） 123312231 这两个特点称为正交性。正是由于正交表具有上述特点，就保证了用正交表安排的试验方案中因素水平是均衡搭配的，数据点的分布是均匀的。因素、水平数愈多，运用正交试验设计方法，愈发能显示出它的优越性，如上述提到的6因素3水平试验，用全面搭配方案需729次，若用正交表L27（313）来安排，则只需做27次试验。在化工生产中， 因素之间常有交互作用。 如果上述的因素T的数值和水平发生变化时，试验指标随因素p变化的规律也发生变化，或反过来，因素p的数值和水平发生变化时，试验指标随因素T变化的规律也发生变化。这种情况称为因素T、p间有交互作用，记为T×p 。5． 3 正交表使用正交设计方法进行试验方案的设计，就必须用到正交表。正交表请查阅有关参考书。1 各列水平数均相同的正交表各列水平数均相同的正交表，也称单一水平正交表。这类正交表名称的写法举例如下： 各列水平均为2的常用正交表有：L4（23），L8（27），L12（211），L16（215），L20（219），L32（231）。各列水平数均为3的常用正交表有：L9（34），L27（313）。各列水平数均为4的常用正交表有：L16（45）各列水平数均为3的常用正交表有：L25（56）2 混合水平正交表各列水平数不相同的正交表，叫混合水平正交表，下面就是一个混合水平正交表名称的写法： L 8（41×24）常简写为L 8（4×24）。此混合水平正交表含有1 个4水平列，4个2水平列，共有1＋4＝5列。3 选择正交表的基本原则一般都是先确定试验的因素、水平和交互作用，后选择适用的L表。在确定因素的水平数时，主要因素宜多安排几个水平，次要因素可少安排几个水平。 （1）先看水平数。若各因素全是2水平，就选用L(2＊)表；若各因素全是3水平，就选L(3＊)表。若各因素的水平数不相同，就选择适用的混合水平表。 （2）每一个交互作用在正交表中应占一列或二列。要看所选的正交表是否足够大，能否容纳得下所考虑的因素和交互作用。为了对试验结果进行方差分析或回归分析，还必须至少留一个空白列，作为“误差”列，在极差分析中要作为“其他因素”列处理。 （3）要看试验精度的要求。若要求高，则宜取实验次数多的L表。 （4）若试验费用很昂贵，或试验的经费很有限，或人力和时间都比较紧张，则不宜选实验次数太多的L表。 （5）按原来考虑的因素、水平和交互作用去选择正交表，若无正好适用的正交表可选，简便且可行的办法是适当修改原定的水平数。（6）对某因素或某交互作用的影响是否确实存在没有把握的情况下，选择L表时常为该选大表还是选小表而犹豫。若条件许可，应尽量选用大表，让影响存在的可能性较大的因素和交互作用各占适当的列。某因素或某交互作用的影响是否真的存在，留到方差分析进行显著性检验时再做结论。这样既可以减少试验的工作量，又不致于漏掉重要的信息。4 正交表的表头设计 所谓表头设计，就是确定试验所考虑的因素和交互作用，在正交表中该放在哪一列的问题。 （1）有交互作用时，表头设计则必须严格地按规定办事。因篇幅限制，此处不讨论，请查阅有关书籍。（2）若试验不考虑交互作用，则表头设计可以是任意的。如在例5-1中，对L 9（3 4）表头设计，表5-3所列的各种方案都是可用的。但是正交表的构造是组合数学问题，必须满足2中所述的特点。对试验之初不考虑交互作用而选用较大的正交表，空列较多时，最好仍与有交互作用时一样，按规定进行表头设计。只不过将有交互作用的列先视为空列，待试验结束后再加以判定。5．4 正交试验的操作方法（1）分区组。对于一批试验，如果要使用几台不同的机器，或要使用几种原料来进行，为了防止机器或原料的不同而带来误差，从而干扰试验的分析，可在开始做实验之前，用L表中未排因素和交互作用的一个空白列来安排机器或原料。与此类似，若试验指标的检验需要几个人（或几台机器）来做，为了消除不同人（或仪器）检验的水平不同给试验分析带来干扰，也可采用在L表中用一空白列来安排的办法。这样一种作法叫做分区组法。 （2）因素水平表排列顺序的随机化。如在例5-1中，每个因素的水平序号从小到大时，因素的数值总是按由小到大或由大到小的顺序排列。按正交表做试验时，所有的1水平要碰在一起，而这种极端的情况有时是不希望出现的，有时也没有实际意义。因此在排列因素水平表时，最好不要简单地按因素数值由小到大或由大到小的顺序排列。从理论上讲，最好能使用一种叫做随机化的方法。所谓随机化就是采用抽签或查随机数值表的办法，来决定排列的别有顺序。 （3）试验进行的次序没必要完全按照正交表上试验号码的顺序。为减少试验中由于先后实验操作熟练的程度不匀带来的误差干扰，理论上推荐用抽签的办法来决定试验的次序。 （4）在确定每一个实验的实验条件时，只需考虑所确定的几个因素和分区组该如何取值，而不要（其实也无法）考虑交互作用列和误差列怎么办的问题。交互作用列和误差列的取值问题由实验本身的客观规律来确定，它们对指标影响的大小在方差分析时给出。 （5）做实验时，要力求严格控制实验条件。这个问题在因素各水平下的数值差别不大时更为重要。例如，例5-1中的因素（加碱量）m的三个水平：m1＝0，m2=5，m3=0，在以m＝m2=5为条件的某一个实验中，就必须严格认真地让m2=5。若因为粗心和不负责任，造成m2=2或造成m2=0，那就将使整个试验失去正交试验设计方法的特点，使极差和方差分析方法的应用丧失了必要的前提条件，因而得不到正确的试验结果。5．5 正交试验结果分析方法 正交试验方法之所以能得到科技工作者的重视并在实践中得到广泛的应用，其原因不仅在于能使试验的次数减少，而且能够用相应的方法对试验结果进行分析并引出许多有价值的结论。因此，有正交试验法进行实验，如果不对试验结果进行认真的分析，并引出应该引出的结论，那就失去用正交试验法的意义和价值。1 极差分析方法下面以表5-4为例讨论L4（23）正交试验结果的极差分析方法。极差指的是各列中各水平对应的试验指标平均值的最大值与最小值之差。从表5-4的计算结果可知，用极差法分析正交试验结果可引出以下几个结论：（1）在试验范围内，各列对试验指标的影响从大到小的排队。某列的极差最大，表示该列的数值在试验范围内变化时，使试验指标数值的变化最大。所以各列对试验指标的影响从大到小的排队，就是各列极差D的数值从大到小的排队。 （2）试验指标随各因素的变化趋势。为了能更直观地看到变化趋势，常将计算结果绘制成图。（3）使试验指标最好的适宜的操作条件（适宜的因素水平搭配）。（4）可对所得结论和进一步的研究方向进行讨论。2 方差分析方法1 计算公式和项目试验指标的加和值＝ ，试验指标的平均值 ，以第j列为例：⑴ Ⅰj ¬＿＿ “1”水平所对应的试验指标的数值之和⑵ Ⅱj —— “2”水平所对应的试验指标的数值之和⑶ ……⑷ kj —— 同一水平出现的次数。等于试验的次数除以第j列的水平数⑸ Ⅰj/ kj —— “1”水平所对应的试验指标的平均值⑹ Ⅱj/ kj —— “1”水平所对应的试验指标的平均值⑺ ……以上7项的计算方法同极差法（见表5-4）。⑻ 偏差平方和 ⑼ fj ——自由度。fj ＝第j列的水平数－1。⑽ Vj ——方差。Vj ＝Sj ／fj 。⑾ Ve ——误差列的方差。Ve ＝Se ／fe 。式中，e为正交表的误差列。⑿ Fj ——方差之比 Fj ＝Vj ／Ve 。⒀ 查F分布数值表（F分布数值表请查阅有关参考书）做显著性检验。⒁ 总的偏差平方和 ⒂ 总的偏差平方和等于各列的偏差平方和之和。即 式中，m为正交表的列数。 若误差列由5个单列组成，则误差列的偏差平方和Se等于5个单列的偏差平方和之和，即：Se ＝Se1 ＋Se2 ＋Se3 ＋Se4 ＋Se5 ；也可用Se ＝S总 ＋S，，来计算，其中S，，为安排有因素或交互作用的各列的偏差平方和之和。2 可引出的结论 与极差法相比，方差分析方法可以多引出一个结论：各列对试验指标的影响是否显著，在什么水平上显著。在数理统计上，这是一个很重要的问题。显著性检验强调试验在分析每列对指标影响中所起的作用。如果某列对指标影响不显著，那么，讨论试验指标随它的变化趋势是毫无意义的。因为在某列对指标的影响不显著时，即使从表中的数据可以看出该列水平变化时，对应的试验指标的数值与在以某种“规律”发生变化，但那很可能是由于实验误差所致，将它作为客观规律是不可靠的。有了各列的显著性检验之后，最后应将影响不显著的交互作用列与原来的“误差列”合并起来。组成新的“误差列”，重新检验各列的显著性。5．6 正交试验方法在化工原理实验中的应用举例例5-2 为提高真空吸滤装置的生产能力，请用正交试验方法确定恒压过滤的最佳操作条件。其恒压过滤实验的方法、原始数据采集和过滤常数计算等见《过滤实验》部分。影响实验的主要因素和水平见表5-5（a）。表中Δp为过滤压强差；T为浆液温度；w为浆液质量分数；M为过滤介质（材质属多孔陶瓷）。 解：（1）试验指标的确定：恒压过滤常数K（m2/s） （2）选正交表：根据表5-5（a）的因素和水平，可选用L 8（4×24）表。 （3）制定实验方案：按选定的正交表，应完成8次实验。实验方案见表5-5（b）。（4）实验结果：将所计算出的恒压过滤常数K（m2/s）列于表5-5（b）。 表5-5（a） 过滤实验因素和水平因素 压强差／kPa 温度／℃ 质量分数 过滤介质符号 Δp T w M水平 1234 88 （室温）18（室温＋15）33 稀（约5%）浓（约10%） G2* G3* * G2 、G3为过滤漏斗的型号。过滤介质孔径：G2 为30~50μm、G3为16~30μm。表2-5（b）正交试验的试验方案和实验结果列号 j＝1 2 3 4 5 6因素 Δp T w M e K（m2/s）试验号 水 平12345678 11223344 12121212 12122121 12211221 12212112 01×10－93×10－21×10－55×10－83×10－02×10－11×10－10×10－3（5）指标K的极差分析和方差分析： 分析结果见表5-5（c）。以第2列为例说明计算过程：Ⅰ2 ＝01×10－4＋21×10－4＋83×10－4＋11×10－4＝92×10－3Ⅱ2 ＝93×10－4＋55×10－4＋02×10－3＋10×10－3＝97×10－3 k2＝4 Ⅰ2/ k2＝92×10-3／4＝79×10－4 Ⅱ2/ k2＝97×10－3／4＝42×10－4 D2＝42×10－4 - 79×10－4＝63×10－4 ∑K＝88×10－3 11×10－4 S2＝k2（Ⅰ2/ k2－ ）2＋k2（Ⅱ2/ k2－ ）2 ＝4（79×10－4－11×10－4 ）2 ＋4（42×10－4－11×10－4 ）2 ＝38×10－7 f2＝第二列的水平数－1＝2－1＝1V2＝S2／f2＝38×10－7／1＝38×10－7Se＝S5＝k5（Ⅰ5/ k5－ ）2＋k5（Ⅱ5/ k5－ ）2 ＝4（22×10－4－11×10－4 ）2 ＋4（99×10－4－11×10－4 ）2 ＝06×10－9 fe＝f5＝1Ve＝Se／fe＝06×10－9／1＝06×10－9F2 =V2／Ve＝38×10－7／06×10－9＝2 查《F 分布数值表》可知：F（а＝01，f1＝1，f2＝1）＝4052 > F2F（а＝05，f1＝1，f2＝1）＝4 >F2F（а＝10，f1＝1，f2＝1）＝9 < F2F（а＝25，f1＝1，f2＝1）＝83 < F2 （其中：f1 为分子的自由度，f2 分母的自由度） 所以第二列对试验指标的影响在 ＝10水平上显著。其他列的计算结果见表2-5（c）。表5－5（c） K的极差分析和方差分析列号 j＝1 2 3 4 5 6因素 Δp T w M e K（m2/s）项目 ⅠjⅡjⅢjⅣjkjⅠj/ kjⅡj/ kjⅢj/ kjⅣj/ kjDjSjfjVjFjF01F05F10F25显著性 94×10－08×10－50×10－61×10－47×10－38×10－52×10－06×10－59×10－65×10－84×10－202*(10) 92×10－97×10－79×10－42×10－63×10－38×10－38×10－832*(10) 04×10－84×10－61×10－61×10－00×10－80×10－80×10－833*(05) 54×10－35×10－35×10－86×10－85×10－70×10－70×10－830*(25) 49×10－40×10－22×10－99×10－30×10－06×10－06×10－00 ∑K＝88×10－3（m2/s）11×10－4（m2/s）（6）由极差分析结果引出的结论：请同学们自己分析。 （7）由方差分析结果引出的结论。 ① 第1、2列上的因素 Δp、T 在 ＝10水平上显著；第3列上的因素w在 ＝05水平上显著；第4列上的因素M 在 ＝25水平上仍不显著。 ② 各因素、水平对K的影响变化趋势见图5-3。图5-3是用表5-5（a）的水平、因素和表5-5（c）的Ⅰj/ kj 、Ⅱj/ kj 、Ⅲj/ kj 、Ⅳj/ k值来标绘的。从图中可看出： A．过滤压强差增大，K值增大； B．过滤温度增大，K值增大； C．过滤浓度增大，K值减小； D．过滤介质由1水平变为2水平，多孔陶瓷微孔直径减小， K值减小。因为第4列对K值的影响在 ＝25水平上不显著，所以此变化趋势是不可信的。 ③ 适宜操作条件的确定。由恒压过滤速率议程式可知，试验指标K值愈大愈好。为此，本例的适宜操作条件是各水平下K的平均值最大时的条件： 过滤压强差为4水平，88kPa 过滤温度为2水平，33℃ 过滤浆液浓度为1水平，稀滤液 过滤介质为1水平或2水平（这是因为第4列对K值的影响在 ＝25水平上不显著。为此可优先选择价格便宜或容易得到者）。上述条件恰好是正交表中第8个试验号。
