混料回归设计.doc

几种常用的混料回归设计返回主页 针对各种回归模型和试验区域以及各种“最优性”要求，人们提出了许多种混料回归设计方案，下面介绍几种常用的方法和应用实例。一、单纯形格子设计 将试验点取在相应阶数的正规单纯形格子点上，这样的试验设计称为单纯形格子设计。 单纯形格子设计是混料回归设计方案中最先出现的，也是最基本的设计方案，很多其它设计方案的构成要用到单纯形格子设计。对于由约束条件(1)构成的正规单纯形因子空间，当采用(4)、(5)、(6)模型形式的完全形规范多项回归模型时，试验点可以取在正规单纯形的格子点上，构成单纯形格子设计。它可以保证试验点分布均匀，而且计算简单、准确，回归系数只是相应格子点的响应值的简单函数。 下面介绍格子点概念与计算公式。 将下图中高为1的等边三角形(a)三条边各二等分，则此三角形(b)的三个顶点与三个边中点的总体称为二阶格子点集，记为3，2，3表示正规单纯形的顶点个数，2表示每边的等分数。将等边三角形(c)各边三等分，对应分点连成与一边平行的直线，在等边三角形上形成许多格子，则这些小等边三角形的顶点，即这些格子的顶点的总体称为三阶格子点集，记为3，3。前面的3指明正规单纯形顶点个数，后面的3指明了每边的等分数。 用类似的方法,可做出其它各种格子点集。三顶点正规单纯形的四阶格子点集记为3，4，总共有15个点。 四顶点正规单纯形(d)的二阶和三阶格子点集分别用4，2和4，3表示，如图(e)和(f)所示。 变量试验X1X2X311002010300141/21/2051/201/2601/21/2图b 3，2 单纯形格子设计 变量试验X1X2X311002010300142/31/3051/32/3062/301/371/302/3802/31/3901/32/3101/31/31/3 图c 3，4 单纯形格子设计 下面一般地介绍(n1)维正规单纯形(有n个顶点)d阶格子点集n，d中各格子点的正规单纯形坐标(重心坐标)的计算法。取n个互相正交的单位向量a1(1，0，.，0)，a2(0，1，.，0)，.，an(0，0，.，1)，则这n个单位向量的顶点便围成一个(n1)维正规单纯形，此正规单纯形上任一点都可以表示为i1a1i2a2.inan (i1，i2，.，in)， (10)其中i1，i2，.，in 0， i1i2.ip1 ， (11)当i1，i2，.，in都取分母是d的分数时，即 则式(11)所确定的总体就是(n1)维正规单纯形的d阶格子点集n，d，也就是说可以算出n，d中各点的单纯形坐标系的坐标。下面我们以n4为例，算出4，2，4，3各点的坐标。 4，2 单纯形格子设计变量试验X1 X2X3X41100020100300104000151/21/20061/201/2071/2001/2801/21/20901/201/210001/21/2 4，3 单纯形格子设计 变量试验X1X2X3X41100020100300104000152/31/30061/32/30072/301/3081/302/3092/3001/3101/3002/31102/31/301201/32/301302/301/31401/302/315002/31/316001/32/3171/31/31/30181/31/301/3191/301/31/32001/31/31/3(1) p4，d2.： 此时a1、a2、a3、a4只有两种取法： 某个a为2，其余为零，有4 个点； 某两个a为1，其余者为零，有6 个点。4，2的10个点如表所示。(2) p4，d3。 此时a1、a2、a3、a4有三种取法： 某个a为3，其余者为零，有4个点； 某个a为2，另外一个为1，其余两个为零，有12个点； 某三个a为1，剩下那个a为零，有4个点。 4，3的20个点如表所示。 单纯形格子设计法 设有n分量系统，各分量Xi(i1，2，.，n)的变化范围满足约束条件(1)，当采用d阶完全型规范多项式回归模型时，试验点选为n，d的C个格子点。 在单纯形格子设计中，n分量d阶格子点集n，d中有C个点，正好与所采用的d 阶完全型规范多项式回归方程中待估计的回归系数的个数相等，故单纯形格子设计是饱和设计，是在“试验次数最少”意义下的最优设计。 常用的单纯形格子设计的试验点数及相应的完全型规范多项式回归方程阶数d 之间的关系如表所示。 单纯形格子设计回归系数的计算 单纯形格子设计中，每个回归系数的值，取决于对应的一些格子点上的响应值，而与其它设计点上的响应值无关，各回归系数都可以表达成相应设计点上响应值的简单线性组合，用计算机建模更方便。下面是三分量混料设计的二阶多项式回归方程。三分量二阶多项式回归方程的规范形式为： YX1X2X3X1X2X1X3X2X3 非单纯形格子设计回归系数的计算 由于实际的试验区域有时是没有规则的几何形状，对于三分量的混料试验，当有上、下界约束条件时，有可能出现如下图所示的几种情况，可用