第八章混料试验的设计 《试验的设计与建模》课件.pptVIP

例 8.4 (续) - 第八章 - 混料试验设计 8.1 引言 合金钢材：以铁为基础再加适量的镊、鉻、錳、碳等成分组成 氖灯灯泡：有氦、氖、氩、氙等多种惰性气体 饮料：含有多种果汁及糖份和水； 中药、酒类产品的制造则更复杂，多种的物质混制而成 例 8.1 咖啡面包 面粉 水 白糖 植物 起酥油 椰子汁 盐 发酵粉 乳化剂 丙酸钙 咖啡粉 香精 食用色素 如何决定各成分的比率? 试验 混料试验设计 设有 s 个因素: X1, ?, Xs 并满足 Xi ? 0, i = 1, ?, s 且 X1 + ? + Xs = 1. 试验区域为单纯形： Ts = {(x1, ?, xs): xi ? 0, i = 1, ?, s , x1 + ? + xs = 1. } 单纯形格子点设计 (Scheffé, 1958) 单纯形重心设计 (Scheffé, 1963) Cox 设计 (1971) 及轴设计 (Cornell, 1975) 常见的混料设计方法: 极大极小或极小极大距离设计 (Johnson, 1990) 模型 一次型 二次型 三次型 由于存在各成分之和为1的限制条件，混料试验的回归方程中中没有常数项、二次项、三次项等，而只有一次项和交互项，采用Scheffé 典型多项式 ? = b1x1 + b2x2 + ? + bs xs 二维及三维标准单纯形格子点集 B. 单纯形重心设计 s 个设计点为单纯形各顶点，即 (1, 0, · · · , 0), · · · , (0, · · · , 0, 1); 个设计点采自单纯形各边的中点，坐标为 (1/2, 1/2, 0, · · · , 0) 的置换; · · · 个设计点采自单纯形各r ? 1 维面的重心，坐标为 (1/r, 1/r, · · · , 0) 的置换; · · · 最后一个设计点为标准单纯形的重心，坐标为 (1/s, 1/s, · · · , 1/s)。 三因素单纯形重心设计试验点 各坐标为 (1, 0, 0), (0, 1, 0), (0, 0, 1), (1/2, 1/2, 0), (1/2, 0, 1/2), (0, 1/2, 1/2), (1/3, 1/3, 1/3)， C. 最优回归设计 一阶多项式模型： 连续设计：设计点为全部 s 个顶点，且权重为1/s， 确定性设计：下列集合中任一设计 二阶模型：当 s=3， 连续设计：三顶点加三条边中点，权重都为1/6 确定性设计：与上面一样 (同单纯形格子点设计) 三阶模型： 十个设计点，不同于单纯形格子点设计 D. Scheffé 型设计 单纯形格子点设计和单纯形重心设计在边界有许多试验点，以致不能做具体的试验。 Scheffé 型设计：将这些设计往单纯形的重心压缩 例 s=3 例 s=3 8.3 混料均匀设计 将 n 个试验点，即 n 种不同的试验配方均匀地散布在Ts 内，而不存在边界上的设计点 怎样设计这些试验点? 逆变换方法 条件法 假设随机向量c = (c1, c2, · · · , cs?1) 服从 s?1 维立方体 Cs?1 上的均匀分布，则从 Cs?1 到标准单纯形 Ts 上的变换 所得的随机向量 x = (x1, x2, · · · , xs) 服从 s ? 1 维标准单纯形 Ts 上的均匀分布。 逆变换方法的理论基础 逆变换方法 给定在超立方体 C s-1 上的均匀设计： {ck = (ck1, ?,ck,s-1), k = 1, ?,n} 做如下的变换: 则 {xk = (xk1, ?,xks), k = 1, ?,n} 为 Ts 上面的均匀设计. 例 8.2 构造一个 n = 12, s = 3 的混料均匀设计。 此时变换 (8.8) 变为: (8.9) 几何意义： T3 x1 x2 x3 给定 (x1, x2, x3)，原坐标 (c1, c2) 如下 c1 c2 不同设计变换后的比较 原来的均匀性好 原来的均匀性差 试验结果 筛选后的模型： ? = 0.317 + 1.823x1 ? 2.568x12 + 1.257x1x2 ?0.947x22, 其 R2 = 0.995，σ2 = 0.0002436 最优解：x1 = 0.4238，x2 =0.2813 时，y 达到最大值 0.7033，此时 x3 = 0.2949 。 在Ts 上的均匀性测度 F-偏差：单纯形上偏差与变换前的偏差一样 均方距离 均方根偏差 最大距离偏差 DM2偏差：直接定义在 Ts 上 8.4 有限制的混料均匀设计