第五篇-单因素实验设计及实验因素水平确定方法ppt课件

实验因素与水平 实验因素与水平 因素 在实验中 影响试验考核指标的量称为因素 水平 水平是试验中各因素的不同取值 实际中一般试验设计的因素水平均取2或3水平 一般用 号或1 2 3 来表示因素的不同水平 当因素只有高低两个水平时 用 号代表高 号代表水平 当因素有3个以上水平时 用1 2 3来依次表示从低到高的水平 在同一试验表中 只能出现同类符号 而不可混用 因素的选取 1 选择依据 专业知识 以往的研究结论 经验教训 最重要的是在阅读文献基础上结合自身实际情况选择 2 一般原则 尽可能全面地考虑到影响试验指标的各个因素 根据实验要求和尽量少选因素 3 首先选对试验指标影响大的因素 尚未完全掌握其规律的因素和未曾被考察研究过的因素 那些对试验指标影响较小的因素 对试验指标影响规律已完全掌握的因素 应当少选或不选 但要作为可控的条件因素参加试验 4 在初步筛选因素时 可以考虑多安排一些因素 因素 在试验设计时 试验因素 输入变量 有两种 一种是在试验时我们可以人为进行控制的可控因素 一种是人为无法控制的噪声 随机 因素 可控因素是在试验过程中我们可以设置和保持其在一个希望的水平上的因子 它应具有以下特征 1 根据经验和以往数据可以确信其对指标Y有重要影响 2 在试验过程中可以比较容易地进行人为改变 可控因子对Y的影响愈大 则潜在的改善机会愈大 噪声因素是试验过程中可使试验结果发生偏差 且无法对其进行控制的因子 它具有以下特征 1 使试验结果偏离目标 2 无法或很难人为控制 当试验中存在噪声因素时 有两种方法可以进行改善 1 首先确认此因素对指标Y的影响程度 如影响大 则须对其进行中和 即直接控制或降低其对Y的影响 2 通过重复精确试验来确定可控因素的最佳水平 当可控因素的水平足够好时 即可得到可靠的设计 对噪声因素不敏感 因素 1 水平有两种 量的变化 数量因素 和质的变化 质量因素 2 数量因素水平水平范围要足够宽 否则就可出现缩小甚至抵削变量影响 同时也看不出因素间交互作用对输出的影响 3 水平设置也不可过宽 否则同样可能缩小此因素的影响 或将其它因素的影响掩盖掉 过宽还可能超出允许操作范围 造成意外损失 一般要求3个以上 4 依据 专业知识 以往的研究结论 经验教训 最重要的是在阅读文献基础上结合自身实际情况选择 水平的选取 水平 确定实验因素 在对实验背景 实验条件 实验预期结果充分了解的基础上结合自身研究内容 现实条件 预期效果确定实验因素 方法 大量阅读文献及总结 确定实验因素水平 文献结合实际 参考单因素优选法 优选法就是根据生产和科研中的不同问题 利用数学原理 合理地安排试验点 减少试验次数 迅速地找到最佳点的一类科学方法 优选法可以解决那些试验指标与因素间不能用数学形式表达 或虽有表达式但很复杂的那些问题 假定f x 是定义在区间 a b 上的函数 但f x 的表达式是并不知道的 只有从试验中才能得出在某一点x0的数值f x0 应用单因素优选法 就是用尽量少的试验次数来确定f x 的最佳点 单因素优选法 单因素优选法 1均分法 在试验范围 a b 内 根据精度要求和实际情况 均匀地排开试验点 在每一个试验点上进行试验 并相互比较 以求得最优点 作法 如试验范围L b a 试验点间隔为N 则试验点n为 包含两个端点 例2 1对采用新钢种的某零件进行磨削加工 砂轮转速范围为420转 分 720转 分 拟经过试验找出能使光洁度最佳的砂轮转速值 N 30转 分 试验转速 420 450 480 510 540 570 600 630 660 690 720 1均分法 使用范围 这种方法的特点是对所试验的范围进行 普查 常常应用于对目标函数的性质没有掌握或很少掌握的情况 即假设目标函数是任意的情况 其试验精度取决于试验点数目的多少 1均分法 适用于试验范围 a b 内 目标函数为单调 连续或间断 的情况下 求最优点的方法 前提是有一个具体指标作为标准 2平分法 每次选取因素所在试验范围 a b 的中点处C做试验 每试验一次 试验范围缩小一半 重复做下去 直到找出满意的试验点为止 根据试验结果 如下次试验在高处 取值大些 就把此试验点 中点 以下的一半范围划去 如下次试验在低处 取值小些 就把此试验点 中点 以上的一半范围划去 2平分法 本方法是在试验范围 a b 内 首先安排两个试验点 再根据两点试验结果 留下好点 去掉不好点所在的一段范围 再在余下的范围内寻找好点 去掉不好的点 如此继续地作下去 直到找到最优点为止 3黄金分割法 0 618法 黄金分割 0 618法要求试验结果目标函数f x 是单峰函数 即在试验范围内只有一个最优点d 其效果f d 最好 比d大或小的点都差 且距最优点d越远的试验效果越差 3黄金分割法 0 618法 设x1和x2是因素范围 a b 内的任意两个试点 C点为问题的最优点 并把两个试点中效果较好的点称为好点 把效果较差的点称为差点 则 最优点与好点必在差点同侧 因而我们把因素范围被差点所分成的两部分中好点所在的那部分称为存优范围 即可以去掉不包含好点的一段 只留下存优范围 如何安排两个试验点 合理地缩小存优范围 3黄金分割法 0 618法 由来 安排两个试点时应该使两个试点关于因素范围的中点对称 对称原则 则无论哪点差 划去的长度都一样 最好每次舍去的区间都能占舍去前全区间同样的比例数 成比例地舍去 原则 由来 3黄金分割法 0 618法 两个试验点位置的确定 设第一次舍去的长度为x 则 3黄金分割法 0 618法 由来 它的左边是第一次舍去的比例数 右边是第二次舍去的比例数 对这个等式进行变形可得 整理可得x 0 382 b a 0 618或 1 0 618 0 382这正是黄金分割常数 3黄金分割法 0 618法 由来 第一个试验点x1设在范围 a b 的0 618位置上 第二个试验点x2取成x1的对称点 则 x1 大 小 0 618 小 b a 0 618 ax2 大 小 第一点 即前一点 b a x1第三个试验点的安排有三种情形 3黄金分割法 0 618法 做法 1 x1是好点 则划去 a x2 保留 x2 b x1的对称点x3 在x3安排第三次试验 x3 大 小 前一点 b x2 x1 由来 3黄金分割法 0 618法 2 x2是好点 则划去 x1 b 保留 a x1 第三个试验点x3应是好点x2的对称点 x3 大 小 前一点 x1 a x2 3黄金分割法 0 618法 3 如果f x1 和f x2 一样 则应该具体分析 看最优点可能在哪边 再决定取舍 一般情况下 可以同时划掉 a x2 和 x1 b 仅留中点的 x2 x1 把x2看成新a x1看成新b 然后在范围 x2 x1 内0 382 0 618处重新安排两次试验试验 3黄金分割法 0 618法 无论何种情况 在新的范围内 又有两次试验可以比较 根据试验结果 再去掉一段或两段试验范围 在留下的范围中再找好点的对称点 安排新的试验 这个过程重复进行下去 直到找出满意的点 得出比较好的试验结果 或者留下的试验范围已很小 再做下去 试验差别不大时也可终止试验 3黄金分割法 0 618法 例 炼某种合金钢 需添加某种化学元素以增加强度 加入范围是1000 2000克 求最佳加入量 第一步先在试验范围长度的0 618处做第 1 个试验x1 大 小 0 618 小 a b a 0 618 1000 2000 1000 0 618 1618克第二步第 2 个试验点由公式x2 大 小 第一点 2000 1000 1618 1382克 3黄金分割法 0 618法 x1 1618克x2 1382克第三步比较 1 与 2 两点上所做试验的效果 现在假设第 1 点比较好 就去掉第 2 点 即去掉 1000 1382 那一段范围 留下 1382 2000 则 x3 大 小 第一点 1382 2000 1618 1764克第四步比较在上次留下的好点 即第 1 处和第 3 处的试验结果 看那个点好 然后就去掉效果差的那个试验点以外的那部分范围 留下包含好点在内的那部分范围作为新的试验范围 如此反复 直到得到较好的试验结果为止 3黄金分割法 0 618法 可以看出每次留下的试验范围是上一次长度的0 618倍 随着试验范围越来越小 试验越趋于最优点 直到达到所需精度即可 3黄金分割法 0 618法 适用范围 试验要求预先给出试验总数 或者知道试验范围和精确度 这时试验总数就可以算出来 在这种情况下 用分数法比0 618法方便 且同样适合单峰函数的方法 裴波那契数列 1 1 2 3 5 8 13 21 34 55 89 144 用F0 F1 F2 依次表示上述数串 它们满足递推关系 Fn Fn 1 Fn 2 n 2 4分数法 所有可能的试验总数正好是某个Fn 1 第一步 先在第Fn 1 Fn 2点上做试验 比较这两个试验的结果 如果第Fn 1点好 划去第Fn 2点以下的试验范围 如果第Fn 2点好 划去Fn 1点以上的试验范围 4分数法 所有可能的试验总数正好是某个Fn 1 第二步 在留下的试验范围中 还剩下Fn 1 1个试验点 重新编号 其中第Fn 2和Fn 3个分点 有一个是刚好留下的好点 另一个是下一步要做的新试验点 两点比较后同前面的做法一样 从坏点把试验范围切开 短的一段不要 留下包含好点的长的一段 这时新的试验范围就只有Fn 2 1试验点 以后 以后的试验 照上面的步骤重复进行 直到试验范围内没有应该做的好点为止 4分数法 例 卡那霉素生物测定培养温度试验 卡那霉素发酵液测定 国内外都规定培养温度为37 1 培养时间在16h以上 某制药厂为缩短时间 决定进行试验 试验范围为29 50 精确度要求 1 中间试验点共有20个 用分数法安排试验 解 由题意可知 试验总次数为20次 正好等于F7 1 4分数法 卡那霉素发酵液测定 国内外都规定培养温度为37 1 培养时间在16h以上 某制药厂为缩短时间 决定进行试验 试验范围为29 50 精确度要求 1 中间试验点共有20个 用分数法安排试验 1 第 个试验点选在13个分点42 第 个试验点在第8个分点37 发现 点好 划去8分点以下的 再重新编号 2 和 比较 好 划去8分点以上的 再重新编号 3 和 比较 好 划去3分点以下的 再重新编号 4 和 比较 好 划去2分点以下的 再重新编号 5 和 比较 好 试验结束 定下43 1 只需8 10h 说明 F7 21 因而最多只需做7 1 6次试验 4分数法 所有可能的试验总数大于某个Fn 1而小于Fn 1 1 只在要试验范围之外虚设几个试验点 凑成Fn 1 1个试验 对于虚设点 并不真正做试验 直接判断其结果比其他点都坏 试验往下进行 很明显 这种虚设点 并不增加实际试验次数 4分数法 例假设某混凝沉淀试验 所用的混凝剂为某阳离子型聚合物与硫酸铝 硫酸铝的投入量恒定为10mg L 而某阳离子聚合物的可能投加量分别为0 10 0 15 0 20 0 25 0 30mg L 试利用分数法来安排试验 确定最佳阳离子型聚合物的投加量 4分数法 解 根据题意可知 可能的试验总次数为5次 由裴波那契数列可知 F5 1 8 1 7F4 1 5 1 4故F4 1 4 5 F5 1 7 1 首先需要增加两个虚设点 使其可能的试验总次数为7次 虚设点可以安排在试验范围的一端或两端 假设安排在两端 即一端一个虚设点 如下图所示 4分数法 4分数法 2 第 个试验点选在第5个分点0 25mg L 第 个试验点在第3个分点0 15mg L 假设 点好 划去3分点以下的 再重新编号 3 和 比较 假设 好 划去2分点以下的 再重新编号 4 此时第 个试验点为虚设点 直接认定它的效果比 差 即 好 试验结束 定下该阳离子型聚合物的最佳投加量为0 30mg L 分数法与0 618法的区别只是用Fn 1和Fn 2代替0 618和0 382来确定试验点 以后的步骤相同 一旦用Fn 1确定了第一个试验点 则以后根据公式确定其余的试验点 也会得出完全一样的试验序列来 4分数法 不管是0 618法 还是分数法 都只是比较两个试验结果的好坏 而不考虑试验的实际值 即目标函数值 抛物线法是根据已得的三个试验数据 找到这三点的抛物线方程 然后求出该抛物线的极大值 作为下次试验的根据 5抛物线法 1 在三个试验点 x1 x2 x3 且x1 x2 x3 分别得试验值y1 y2 y3 根据拉格朗日插值法可以得到一个二次函数 2 设上述二次函数在x4取得最大值 这时 5抛物线法 3 在x x4处做试验 得试验结果y4 如果假定y1 y2 y3 y4中的最大值是由x4给出的 除x4之外 在x1 x2 x3和x4中取较靠近x4的左右两点 将这三点记为x1 x2 x3 此处x1 x2 x3 若在x1 x2 x3 处的函数值分别为y1 y2 y3 则根据这三点又可得到一条抛物线方程 如此继续下去 直到找到函数的极大点 或它的充分邻近的一个点 被找到为止 5抛物线法 某离心泵效率与流量之间关系曲线的试验中 已经测得三组数据如下表所示 求流量为多大时 效率最高 看抛物线法做的过程 5抛物线法 三个试验点得到一条抛物线 该抛物线有一极值点 可根据抛物线方程求出该极值点X座标 X4 该流量处对应效率可能是最大的 在X4对应流量处安排试验 测定其效率 假定测得效率值为78 即X4流量对应效率最大 则下一步试验 X4对应流量最大 因而取较靠近X4的左右两点 再做抛物线 即以三点 20 75 24 78 32 70 做抛物线 求得一极值点 何时试验结束 粗略地说 如果穷举法 在每个试验点上都做试验 需要做n次试验 对于同样的效果 黄金分割法只要数量级lgn次就可以达到 抛物线法 效果更好些 只要数量级lglgn次 原因就在于黄金分割法没有较多地利用函数的性质 做了两次试验 比一比大小 就把它舍掉了 抛物线法则对试验结果进行了数量方面的分析 5抛物线法 抛物线法常常用在0 618法或分数法取得一些数据的情况 这时能收到更好的效果 此外 建议做完了0 618法或分数法的试验后 用最后三个数据按抛物线法求出x4 并计算这个抛物线在点x x4处的数值 预先估计一下在点x4处的试验结果 然后将这个数值与已经试得的最佳值作比较 以此作为是否在点x4处再做一次试验的依据 5抛物线法 在生产和科学实验中 为加速试验的进行 常常采用一批同时做几个试验的方法 即分批试验法 若可同时做多个试验 则一次做一个好 还是一次同时做几个好 一个试验后试验范围缩小多少 时间 成本 以下介绍的几种方法同样适用于单峰函数 6分批实验法 1 预给要求法 预给要求法是分批试验的一种方法 如能预先确定总的可能的试验个数 换句话说 知道了试验范围和要求的精密度 或事先限定试验的批数和每批的个数 就可以采用这种方法 1 预给要求法 每批做偶数个试验 若只做一批试验 每批两个试验 把试验范围平分3等份 在每个分点上做试验 如下所示 若做两批试验 每批两个试验 把试验范围分为7等份 在第3 4两点做第一批试验 如第4点好 再做5 6两点 如第3点好 则做1 2两点 若做三批试验 每批两个试验 把试验范围分为15等份 在第7 8两点做第一批试验 如第7点好 则把第8点以上的范围划去 如第8点好 则把7点以下的划去 在余下的部分做第二批试验 如下所示 依此可以推出做更多批数试验的情形来 1 预给要求法 每批做偶数个试验 每批做4个 6个或更多个试验的情形原理相同 容易推出 若每批做2k个试验 共作n批 则应将试验范围等分份 第一批试验点是 试验结果的精确度是 L b