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

  • 147
  • 0
  • 约1.32万字
  • 约 25页
  • 2017-09-07 发布于陕西
  • 举报
第一节 一次回归正交设计 一 正交设计和回归设计的特点 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),根据最小二乘

文档评论(0)

1亿VIP精品文档

相关文档