精品回归正交试验设计.docx

1、精品回归正交试验设计【关键字】精品8 回归正交试验设计 前面介绍的正交试验设计是种很实用的试验设计方法，它能利用较少的试验次数获得较好的试验结果，但是通过正交设计所得到的优方案只能限制在已定的水平上，而不是一定试验范围内的最优方案；返回分析是一种有效的数据处理方法，通过所确立的返回方程，可以对试验结果进行预测和优化，但返回分析往往只能对试验数据进行被动的处理和分析，不涉及对试验设计的要求。如果能将两者的优势统一起来，不仅有合理的试验设计和较少的试验次数，还能建立有效的数学模型，这正是我们所期望的。返回正交设计(orthogonal regression design)就是这样一种试验设计方法，

2、它可以在因素的试验范围内选择适当的试验点，用较少的试验建立一个精度高、统计性质好的返回方程，并能解决试验优化问题。8.1 一次返回正交试验设计及结果分析一次返回正交设计就是利用返回正交设计原理，建立试验指标(y)与m个试验因素x1，x2，xm，之间的一元返回方程： (81)或者 k=1，2，m1（jk） (82) 一次返回正交设计的基本方法 (1)确定因素的变化范围根据试验指标y，选择需要考察的m个因素xj (j1，2，m)，并确定每个因素的取值范围。设因素xj的变化范围为xj1，xj2，分别称xj1和xj2为因素xj的下水平和上水平，并将它们的算术平均值称作因素xj的零水平，用xj0。表示。

3、 (83) 上水平与零水平之差称为因素xj的变化间距，用j表示，即： (84) 或 (85) (2)因素水平的编码编码(coding)是将xj的各水平进行线性变换，即： (86) 式(86)中zj就是因素xj的编码，两者是一一对应的。显然，与xj1，xj0和xj2的编码分别为1，0和1，即zj11，zj20，zj21。一般称xj为自然变量，zj为规范变量。因素水平的编码结果可表示成表81。对因素xj的各水平进行编码的目的，是为了使每个因素的每个水平在编码空间是“平等”的，即规范变量zj的取值范围都在1，1内变化，不会受到自然变量xj的单位和取值大小的影响。所以编码能将试验结果y与因素zj (j

4、1，2，m)各水平之间的返回问题，转换成试验结果y与编码值zj之间的返回问题，从而大大简化了返回计算量。表81 因素水平编码表规范变量zj自然变量xjx1x2xm下水平(1)x11x21xm1上水平(1)x12x22xm2零水平(0)x10x20xm0变化间距j12m (3)一次返回正交设计表 将二水平的正交表中“用“代换，就可以得到一次返回正交设计表。例如正交表L8(27)经过变换后得到的返回正交设计表如表82。表82 一次返回正交设计表试验号列 号l23456711111ll121l111113111111l411111115111l1116111111171lll11181111ll1

5、代换后，正交表中的编码不仅表示因素的不同水平，也表示了因素水平数值上的大小。从表82可以看出返回正交设计表具有如下特点：1任一列编码的和为0，即： (87)所以有 ， (88)2任两列编码的乘积之和等于零，即：， (89) 这些特点说明了转换之后的正交表同样具有正交性，可使返回计算大大简化。(4)试验方案的确定 与正交试验设计类似，在确定试验方案之前，要将规范变量zj安排在一元返回正交表相应的列中，即表头设计。 例如，需考察三个因素x1、x2、x3，可选用L8(27)进行试验设计，根据正交表L8(27)的表头设计表，应将x1、x2、x3分别安排在第1、2和4列，也就是将zl，z2，z3安排在表

6、82的第1、2和4列上。如果还要考虑交互作用x1 x2、x1x3，也可参考正交表L8(27)的交互作用表，将zlz2和z2z3，分别安排在表82的第3、5列上，表头设计结果见表83。每号试验的方案由zl，z2，z3z对应的水平确定，这与正交试验是一致的。表8-3 三因素一次返回正交表试验号12345Z1Z2Z1 Z2Z3Z1 Z3123451111111111111111111111l116789101ll001110011100111001l10O从表83可以看出，第3列的编码等于第1，2列编码的乘积，同样第5列的编码等于第1，4列编码的乘积，即交互作用列的编码等于表中对应两因素列编码的乘积

7、，所以用返回正交表安排交互作用时，可以不参考正交表的交互作用表，直接根据这一规律写出交互作用列的编码，这比原正交表的使用更方便。 表8-3中的第9，10号试验称为零水平试验或中心试验。安排零水平试验的目的是为了进行更精确的统计分析(如返回方程的失拟检验等)，得到精度较高的返回方程。当然，如果不考虑失拟检验，也可不安排零水平试验。 一次返回方程的建立建立返回方程，关键是确定返回系数。设总试验次数为n，其中包括mc次二水平试验(原正交表所规定的试验)和m0次零水平试验，即： (810) 如果试验结果为yi (i1，2，n)，根据最小二乘法原理和返回正交表的两个特点，可以得到一次返回方程系数的计算公

8、式如下(证明略)： (811)， (812)， ， (813) 上述式中，zji表示zi列各水平的编码，(zjzk)i表示zjzk列各水平的编码。 需要指出的是，如果一次返回方程中含有交互作用项zjzk (jk)，则返回方程不是线性的，但交互作用项的返回系数的计算和检验与线性项zj是相同的，这是因为交互作用对试验结果也有影响，可以被看作是影响因素。 通过上述方法确定偏返回系数之后，可以直接根据它们绝对值的大小来判断各因素和交互作用的相对重要性，而不用转换成标准返回系数。这是因为，在返回正交设计中，所有因素的水平都经过了无因次的编码变换，它们在所研究的范围内都是“平等的”，因而所求得的返回系数不

9、受因素的单位和取值的影响，直接反映了该因素作用的大小。另外，返回系数的符号反映了因素对试验指标影响的正负。 返回方程及偏返回系数的方差分析 (1)无零水平试验时首先计算各种平方和及自由度。总平方和为： (814)其自由度为。根据式(444)和式(812)，推导出一次项偏返回平方和的计算公式为： (444)， (815)同理可以得到交互项偏回归平方和的计算公式：， (816) 各种偏回归平方和的自由度都为1。一次项偏回归平方和与交互项偏回归平方和的总和就是回归平方和： (817)所以回归平方和的自由度也是各偏回归平方和的自由度之和： (818)于是残差平方和为： (819)其自由度为： (820

10、)如果考虑了所有的一次项和交互项，则可参照表84进行方差分析。表8-4 一次回归正交设计的方差分析差异源111111回归残差总和 在实际做试验时，往往只需要考虑几个交互作用，或者可以不考虑交互作用，所以在计算回归和残差自由度时应与实际情况相符。如果不考虑交互作用，。值得注意的是，无论是否考虑交互作用，都不影响偏回归系数的计算公式。 经偏回归系数显著性检验，证明对试验结果影响不显著的因素或交互作用，可将其从回归方程中剔除，而不会影响到其他回归系数的值，也不需要重新建立回归方程。但应对回归方程再次进行检验，将被剔除变量的偏回归平方和、自由度并入到残差平方和与自由度中，然后再进行相关的分析计算。 (

11、2)有零水平试验时 如果零水平试验的次数m02，则可以进行回归方程的失拟性(1ack of fit)检验。前面对回归方程进行的显著性检验，只能说明相对于残差平方和而言，各因素对试验结果的影响是否显著。即使所建立的回归方程是显著的，也只反映了回归方程在试验点上与试验结果拟合得较好，不能说明在整个研究范围内回归方程都能与实测值有好的拟合。为了检验一次回归方程在整个研究范围内的拟合情况，则应安排m0(m02)次零水平试验，进行回归方程的失拟性检验，或称拟合度检验(test of goodness of fit)。 设m0次零水平试验结果为y01，y02，y0 m0，根据这m0次重复试验，可以计算出重

12、复试验误差为： (821)试验误差对应的自由度为： (822)由前述知，只有回归系数与零水平试验次数m0有关，其他各偏回归系数都只与二水平试验次数m0有关，所以增加零水平试验后回归平方和SSR没有变化，于是定义失拟平方和为： (823)或 (824)可见，失拟平方和表示了回归方程未能拟合的部分，包括未考虑的其他因素及各xj的高次项等所引起的差异。它对应的自由度为： (825)所以有： (826) (827)这时 (828)服从自由度为(，)的F分布。对于给定的显著性水平(一般取0.1)，当 时，就认为回归方程失拟不显著，失拟平方和是由随机误差造成的，否则说明所建立的回归方程拟合得不好，需要进一

13、步改进回归模型，如引入别的因素或建立更高次的回归方程。只有当回归方程显著、失拟检验不显著时，才能说明所建立的回归方程是拟合得很好的。 最后需要指出的是，回归正交试验得到的回归方程是规范变量与试验指标之间的关系式，还应对回归方程的编码值进行回代，得到自然变量与试验指标的回归关系式。 【例81】 用石墨炉原子吸收分光光度计法测定食品中的铅，为提高测定灵敏度，希望吸光度(y)大。为提高吸光度，讨论了xl(灰化温度)、x2(原子化温度)和x3(灯电流mA)三个因素对吸光度的影响，并考虑交互作用x12，x13。 已知x1300700，x218002400，x3810mA。试通过回归正交试验确定吸光度与三

14、个因素之间的函数关系式。 解：(1)因素水平编码 因x1300700，所以其上水平x12700，下水平x11300，零水平，变化间距，以x11300为例，对应的编码。同理可对其他因素水平进行编码，编码结果见表85。表85 例81因素水平编码表编码zj灰化温度xl原子化温度x2灯电流x3mA上水平(1)下水平(1)零水平(0)变化间距j70030050020024001800210030010891 (2)正交表的选择和试验方案的确定依题意，可以选用正交表L8(27)，经编码转换后，得到表82所示的回归正交表。如表86所示，将zl、z2、z3分别安排在第1，2和4列，则第3列和第5列分别为交互作

15、用z1z2、z2z3；列。不进行零水平试验，故总试验次数n8，试验结果也列在表86中(注：本例的试验方案和试验结果与例65是完全一样的)。表8-6 例81三元一次回归正交设计试验方案及试验结果试验号z1z2z1 z2z3z1 z3灰化温度xl原子化温度x2灯电流x3mA吸光度yi1111l17002400100.55221l1l1700240080.55431ll117001800100.4804111ll700180080.4725111113002400100.516611l1l300240080.5327l-11113001800100.44881-1111300180080.484 (

16、3)回归方程的建立 依题意，m00，nmc8。根据回归系数的计算公式，将有关计算列在表87中。表87 例81三元一次回归正交设计计算表试验号z1z2z1z2z3z1 z3yy2z1yz2yz3y(z1z2)y(z1 z3)yl23456781l1l1111111l1111111111111111111111111l110.5520.5540.4800.4720.5160.5320.4480.4840.3047040.3069160.2304000.222784O.266256O.2830240.2007040.2342560.5520.5540.4800.4720.5160.5320.4480

17、.4840.5520.554O.4800.4720.5160.5320.4480.4840.5520.5540.4800.4720.5160.5320.4480.4840.5520.5540.4800.4720.5160.5320.4480.4840.5520.5540.4800.4720.5160.5320.4480.4844.0382.0490440.0780.2700.0460.0380.058由表87得：所以回归方程为 由该回归方程中偏回归系数绝对值的大小，可以得到各因素和交互作用的主次顺序为：x2xlx1x3x3x1x2，这与例65中正交试验的分析结果是一样的。(4)方差分析方差分析

18、结果见表88。由表88，对于显著性水平0.05，只有因素z2对试验指标y有非常显著的影响，其他因素和交互作用对试验指标都无显著影响，所以应将z1，z3，zlz3，zlz2的平方和及自由度并入残差项，然后再进行方差分析。这时的方差分析为一元方差分析，分析结果见表89。表88 例81方差分析表差异源SSdfMSF显著性z10.00076l10.00076112.27z20.009113l0.009113146.98z30.00026510.0002654.27z1 z30.000l8l10.00018l2.92zl z20.00042110.0004216.79回归0.01074150.00214

19、834.65残差e0.00012320.000062总和0.010864n17注：F0.05(1，2)18.51，F0.01 (1，2)=98.49， F0.05 (5，2)19.30，F0.01 (5，2)=99.30。表89 例81第二次方差分析表差异源SSdfMSF显著性回归(z2)0.009113l0.00911331.21残差e0.00175160.000292总和0.010864n17注：F0.05(1，6)5.99，F0.01(1，6) 13.74。 由表89可知，因素z2对试验指标y有非常显著的影响，因此原回归方程可以简化为：可见，只有原子化温度z2对吸光度有显著影响，两者之间

20、存在显著的线性关系，而且原子化温度取上水平时试验结果最好。根据编码公式，将上述线性回归方程进行回代：整理后得到： 【例82】 从某种植物中提取黄酮类物质，为了对提取工艺进行优化，选取三个相对重要的因素：乙醇浓度(x1)、液固比(x2)和回流次数(x3)进行了回归正交试验，不考虑交互作用。已知xl6080，x2812，x313次。试通过回归正交试验确定黄酮提取率与三个因素之间的函数关系式。 解：(1)因素水平编码及试验方案的确定表8-10 例8-2因素水平编码表编码zj乙醇浓度液固比回流次数16081070102180123j1021由于不考虑交互作用，所以本例要求建立一个三元线性方程。因素水平

21、编码如表810所示。选正交表L8(27)安排试验，将三个因素分别安排在回归正交表的第1、2、4列，试验方案及试验结果见表811，表中的第9、10、11号试验为零水平试验。表8-11 例8-2试验方案及试验结果试验号z1z2z3乙醇浓度液固比回流次数提取率y1111801238.02111801217. 3311180836. 94l1l808l6.45111601236.96111601216.5711l60836.0811160815.19000701026.610000701026.511000701026.6回归方程的建立 将有关计算过程列在表812中。表812 例82试验结果及计算表试

22、验号z1z2z3提取率yy2z1yz2yz3y11118.064.008.08.08.021117.353.297.37.37.33l116.947.616.96.96.941116.440.966.46.46.451116.947.616.96.96.961116.542.256.56.56.571116.036.006.06.06.081115.126.015.15.15.190006.643.560.00.00.0100006.542.250.00.00.0110006.643.560.00.00.072.8487.14.14.32.5由表812计算回归方程为 由该回归方程偏回归系数绝对

23、值的大小，可以得到各因素的主次顺序为：x2x1x3，即液固比乙醇浓度回流次数。又由于各偏回归系数都为正，所以这些影响因素取上水平时，试验指标最好。 (3)回归方程显著性检验有关平方和的计算如下：方差分析结果见表813。表8-13 例8-2方差分析表差异源SSdfMSF显著性z12.10112.101142.9z22.31112.311157.2z30.78110.78153.1回归5.19331.731117.8残差0.10370.0147总和5.296n110 注：F0.01(1，7)12.25，F0.01(3，7) 8.45 可见，三个因素对试验指标都有非常显著的影响，所建立的回归方程也非

24、常显著。 (4)失拟性检验本例中，零水平试验次数m03，可以进行失拟性检验，有关计算如下。检验结果表明，失拟不显著，回归模型与实际情况拟合得很好。 (5)回归方程的回代根据编码公式：，代人上述回归方程得：整理后得到：8.2 二次回归正交组合设计 在实际生产和科学试验中，试验指标与试验因素之间的关系往往不宜用一次回归方程来描述，所以当所建立的一元回归方程经检验不显著时，就需用二次或更高次方程来拟合。8.2.1 二次回归正交组合设计表 组合设计试验方案的确定假设有m个试验因素(自变量)zj (j1，2，m)，试验指标为因变量y，则二次回归方程的一般形式为：，其中，为回归系数，可以看出该方程共有项，要使回归系数的估算成为可能，必要条件为试验次数；同时，为了计算出二次回归方程的系数，每个因素至少要取3个水平，所以用一元回归正交设计的方法来安排试验，往往不能满足这一条件。例如，当因素数m3时，二次回归方程的项数为10，要求试验次数n10，如果用正交表L9(34)安排试验，则试验次数不符合要求，如果进行全面试验，则试验次数为3327次，试验次数又偏多。为解决这一矛盾，可以在一次回归正交试验设计的基础上再增加一些特定