工程试验设回归正交试验设计.doc

第一节 一次回归正交设计 一 正交设计和回归设计的特点 1 正交设计的特点 正交设计是一种很实用的试验设计方法，它利用较少的试验次数获得较好的试验结果；但是通过正交设计得到的优方案只是局限在确定的水平组合中，而不是一定试验范围内的最优方案。 2 回归设计 回归分析是一种有效的数据处理方法，通过所确定的回归方程，可对试验结果进行预测和控制；但是，它只能对试验数据进行被动的分析和处理，不涉及对试验设计的要求。 如果把两者的优势统一起来，不仅有合理的试验设计和较少的试验次数，还能建立有效的数学模型，这就是回归正交设计方法。 二 一次回归正交设计基本方法 一次回归正交设计就是利用回归正交设计原理，建立试验指标y与m个因素x1、x2、…、xm之间的一次回归方程： （k=1，2，…，m） 如果不考虑交互作用，则一次回归方程为 一次回归正交设计的基本步骤如下： 1 确定因素的变化范围 根据指标y，确定需要考察的m个因素xj（j=1，2，…，m），并确定每个因素的取值范围。 设：xj的变化范围为[xj1，xj2]，分别称xj1和xj2为因素xj的下水平和上水平，并将其算术平均值称为零水平，即 上水平与零水平之差或零水平与下水平之差称为xj的变化间距，即 例如，某试验中温度的变化范围为30-90℃，则其上水平为xj2=90℃，xj1=30℃，零水平xj0=60℃，变化间距△j=30℃。 2 因素水平的编码 编码（coding）就是将xj的各水平进行线性变换，即 式中，zj——xj的编码。 显然，xj1、xj0、xj2的编码分别为-1、0、+1，即zj1=-1，zj0=0，zj=+1。一般，xj称为自变量，zj称为规范变量。 因素水平的编码见表8-1规范变量zj 自然变量xj x1 x2 … xm 下水平（-1） 零水平（0） 上水平（+1） 变化区间△j x11 x10 x12 △1 x21 x20 x22 △2 xm1 xm0 xm2 △m 对因素水平进行编码的目的，是为了使每个因素的每个水平在编码空间是“平等”的，即规范变量zj的取值范围在[-1，+1]内变化，不会受到自然变量xj的单位和取值大小的影响。所以编码能将试验结果y与因素xj各水平之间的回归问题，转换为试验结果y与因素zj之间的回归问题，从而简化了回归计算量。 3 一次回归正交设计编码表 将二水平正交表中的2”用-1”代换，就可得到一次回归正交设计编码表。 例如，L8(27)经过变换后得到的回归正交设计表如表8-2所示。 代换后，正交表中的编码不仅表示因素的不同水平，也表示因素水平数值上的大小。从表8-2可看出回归正交设计表具有如下特点： （1）任意一列编码的和为零，即 （2）任意两列编码的乘积之和为零，即 这些特点说明了转换之后的正交表同样具有正交性。 表8-2 一次回归正交设计编码表 试验号 列号 1 2 3 4 5 6 7 1 2 3 4 5 6 7 8 1 1 1 1 -1 -1 -1 -1 1 1 -1 -1 1 1 -1 -1 1 1 -1 -1 -1 -1 1 1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 -1 1 -1 1 1 -1 -1 1 1 -1 -1 1 1 -1 -1 1 -1 1 1 -1 4 试验方案的确定 在确定试验方案时，也要将规范变量zj安排在一次回归正交编码表相应的列中，即进行表头设计。 例如，考察因素x1、x2和x3，可选用L8(27)，根据L8(27)的表头设计表，应将x1、x2和x3分别安排在第1、2和4列，即将z1、z2和z3安排在表8-2的第1、2和4列上。如果要考虑交互作用x1x2和x1x3，也可参考L8(27)的交互作用表，将z1z2和z1z3分别安排在表8-2的第3、5列上。表头设计结果见表8-3。 表8-3 三因素一次回归正交表 试验号 1 2 3 4 5 z1 z1 z1z2 z1 z1z3 1 2 3 4 5 6 7 8 9 10 1 1 1 1 -1 -1 -1 -1 0 0 1 1 -1 -1 1 1 -1 -1 0 0 1 1 -1 -1 -1 -1 1 1 0 0 1 -1 1 -1 1 -1 1 -1 0 0 1 -1 1 -1 -1 1 -1 1 0 0 从表8-3中看出，第3列的编码等于第1列和第2列编码的乘积，即交互作用列的编码等于表中对应两因素列编码的乘积。 表8-3中的第9、10号试验成为零水平试验或中心试验。安排零水平试验的目的是为了进行更精确的统计分析，得到精度较高的回归方程。 三 一次回归方程的建立 建立回归方程，关键是确定回归系数。 设总试验次数为n，其包括mc次二水平试验和m0次零水平试验，即 如果试验结果为yi（i=1，2，…，n），根据最小二乘