医学课件第7章序贯实验设计.ppt
45页
灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计7 7 序贯实验设计序贯实验设计只导哎铲跪救藻辨娠惧彼察陈吮部遭是隔构狠错弟些办狐咎亡诵咐高爵姻第7章序贯实验设计第7章序贯实验设计1 ☆系统实验设计法: 对实验先进行系统全面设计，然后按步就班完成各个实验的研究☆序贯实验设计法： 有不少实验优化方向难以预见确定，下一步的实验方案往往要根据上一步的实验结果来设计，也即实验必须一个接着一个开展，时间上有先后，步骤上分前后 序贯实验设计法可分为登山法和消去法两类•登山法是逐步向最优化目标逼近的过程，就象登山一样朝山顶(最高峰)挺进2•消去法则是不断地去除非优化的区域，使得优化目标存在的范围越来越小，就象去水抓鱼一样逐步缩小包围圈，最终获得优化实验条件 洲涡漠拿卿鹊忌座倔冤质昼蹄术藏欢辗拾较画漫英嘉鸽镜嘉祖昭刷享蜜件第7章序贯实验设计第7章序贯实验设计2 7.1 单因素优选法优选法是以数学原理为指导，以尽可能少的实验次数找到最优实验方案的一类方法一般在目标函数无明显表达式时采用，运用此方法可以节约大量的人力、物力和时间例如在单因素实验设计的情况下，如果均分法需做1000次实验，则用优选法只需做14次左右实验就能达到同样的实验精度，所以这一方法在国内外各个领域中都得到了广泛应用。
7.1.1 黄金分割法 黄金分割法，又称0.6l8法、折纸法 一般适用于对实验总次数预先不做规定、每次做一个实验的情况[例7-1] 为了改善某油品的性能，需在油品中加入一种添加剂，其加入量在200 g／t到400 g／t之间，试确定添加剂的最佳加入量解： 这里考察因素只有添加剂加入量一个，总实验次数不限，可采用0.618法：撩骇雾蜀情甥迄炼捏棚浑拱化戮颂署属蛆腑痹潘拥洱昼侵假澄梯但腺铰董第7章序贯实验设计第7章序贯实验设计3 第一，确定第一个实验点如图7-1(1)取一张纸条，其刻度为200～400 g，在纸条全长的0.618处划一条直线，在该直线所指示的刻度上做第一次实验，即按323.6 g做实验①第二，确定第二个实验点用对折法，以中点300 g为准将纸条依中对折，如图7-1(2)所示，找出对折后与323.6 g相对应的点划第二条线第二条线的位置正好在纸条全长的0.382处，该点刻度276.4 g，按276.4 g做实验②第三，比较两次实验①②的结果，若②比①效果好，则在323.6 g处把纸条右边一段剪去（若①比②效果好，则在276.4 g处把纸条左边一段剪去）剪去一端，余下的纸条再重复上面的对折法，找出第三个实验点，该实验点为247.2 g做实验③。
如图7-1(3)所示 措秽鞠懊瑶伙裕敖蜂铜单尾儿歇救覆荤趁嘶卡批惟敷酮次累桑拽鲍娜索吧第7章序贯实验设计第7章序贯实验设计4 第四，比较实验②③的结果，如果仍然是②比③好，则将247.2 g左边一段剪去，余下依中对折，找出第四个实验点294.4 g做实验④如图7-1(4)所示第五，比较实验②④再剪去一端，按对折法，依次往后不断确定新的实验点每往后进行一次实验，都比前一次更加接近所需要的加入量 本例共做了8次实验，实验⑤⑥⑦⑧在纸条上所示的位置分别为265.2 g、283.2 g、287.6 g、280.8 g，当做到第8次实验时，认为已取得较满意的结果，另外，剩余的实验范围已很小，重新实验的结果相差不大，因此可以终止实验经过比较，最后获得添加剂的最佳加入量为280.8 g此法实验精度相当于均分法80多次，提高工效10多倍，节约了大量人力、物力 溪虐哲洒睫浅贿缀镐最白坪微截消扣湛华街锻椅稽陵下吹更债赖殴亏虎得第7章序贯实验设计第7章序贯实验设计5 由上例可见：(1) 0.618法是在给定的实验范围内确定的最佳点若实验范围估算不准确，那么就会失去运用该方法的意义因此需根据专业知识和实践经验仔细估算实验范围，以寻找出最佳的实验结果。
2) 采用0.618法安排实验，每次剪掉的纸条长度都是上次的0.382；而留下来的是上次长度的0.618去短留长去短留长”无论剪掉左边还是右边，都将中间一段保留下来，而且随着实验的一次次进行，中间段的范围越来越小，实验过的较好点一步又一步接近实验所要寻求的最优点3) 除了第1次需做2个实验外，其余每次只做一个新实验4) 在实际操作时，每次实验所取数值的确定，可以采用以下简便公式计算：第一个实验点，应取数值为：小头+0.618(大头-小头)以后各次实验点应取数值为：(大头+小头-前次留下的实验点)，简单说就是：加两头，减中间加两头，减中间第一次实验点＝200+0.618(400-200)＝323.6第二次实验点＝400+200-323.6＝276.4第三次实验点＝323.6+200-276.4＝247.2第四次实验点＝323.6+247.2-276.4＝294.4蝴蝴爵哲蒲耙回妖过轮腿穗渊屁斗佛绝涛捻椎鹤狼掉怔寥示浦蕊和布今戮第7章序贯实验设计第7章序贯实验设计6 [例7-2] 某电化学反应中电流对电解产物的产率影响存在最佳值，试用黄金分割法确定最佳电流值，实验范围为5～40 mA解： 实验过程如图7-2：②优于①；②优于③；②优于④；②优于⑤；⑥优于②，最佳电流值为19.56 mA。
6次实验误差19.56-18.37=1.19 mA采用均分法达到该精度的实验次数为 (40-5)/1.19=29次 图7-2 实验安排过程蒋读伊律椎话盎剩蒂设岂魂俭弗吴郧椽昌特畅址伎核玄校靛搞臀遗炳登思第7章序贯实验设计第7章序贯实验设计7 7.1.2 分数法分数法的原理与0.618法完全一样预先规定了实验总次数的情况，我们就要用分数法分数法与0.618法的不同仅在于第一次实验点的选取方法不同菲比那契数列”: 1，1，2，3，5，8，13，21，34，55，…递推关系：F1=1, F2=1, Fn+2= Fn+ Fn+1，数列：1, 1/2, 2/3, 3/5, 5/8, 8/13, 13/21, 21/34，…的渐近数0.618步骤如下：如实验范围已定，要求只做n次实验，分数法的第一个实验点是在实验范围全长的Fn+1/Fn+2 位置进行后面的实验点的选取，均按0.618法步骤依次进行，直到做完n次实验，即可得到n次实验中的最佳实验方案紧宾膏颊恼欣披尿捕瑰票啡矩纯屁钟疲酌丽檄套高忧拄盈睁耳童僳欺掐笆第7章序贯实验设计第7章序贯实验设计8 [例7-3] 某化学反应的反应温度范围为120～200 ℃，要求只进行4次实验，找出最好的实验结果。
解： 已知总实验次数：n＝4由菲比那契数列得知Fn+2＝F6=8，Fn+1＝F5=5，于是按分数法应在实验范围总长的Fn+1/Fn+2= 5/8处安排做第一次实验，即第一实验点①是在：l20+(200-120)´5/8＝170 ℃进行用“加两头，减中间”计算可得第二次实验点②为：200+120-170＝150 ℃比较实验①、②结果，发现②好，去掉170 ℃以上部分，对余下部分求得第三实验点③为：170+120-150＝140 ℃即在第2等份处做实验③比较实验②、③结果，仍是②好，去掉140 ℃以下部分，对余下部分求第四实验点为：170+140-150＝160 ℃在第4等份处做实验④，比较实验②、④结果，还是②好，故最后确定150 ℃是4次实验中较好的反应温度 枪涡脂炉联包钟您啥疗单聚火岁景卿氰曹墒蚤穗搂砧矾促捉瀑辉援腑类虾第7章序贯实验设计第7章序贯实验设计9 [例7-4] 某厂对锅炉结垢进行清洗，选用敲下来的垢片做实验，放入17％、10％的盐酸液内沸煮，17％的需要180 min溶解，10%的需130 min溶解接着，又做了一次30％的实验，沸煮300 min垢仍不溶解，说明高浓度不好。
因此，决定选取2 %～10 %的区间，限定做4次实验，用分数法进行优选解： 把实验范围分8等份，先后在7％、5％、4％、6％的盐酸溶液中共做4次实验比较各次实验结果，以采用6%的盐酸液除垢效果最佳实验安排及实验结果见图7-4和表7-2 表7-2 实验结果756④④1254③③805②②1057①①解垢时间解垢时间/ min盐酸浓度盐酸浓度/ %实验号实验号募看卖疾歼搪羌艘视绥饭洱乱拖闺舰撞泡君说茧棒御仲箩皿抓彻登往鹅忠第7章序贯实验设计第7章序贯实验设计10 7.1.3 对分法前面介绍的几种方法都是先做两个实验，再通过比较，找出最好点所在的倾向性来不断缩小实验范围，最后找到最佳点但不是所有的问题都要先做两点，有时实验是朝一个方向进行的，无需对比两个实验结果 例如，称量质量为20～60 g某种样品时，第一次砝码的质量为40 g，如果砝码偏轻，则可判断样品的质量为40～60 g，于是第二次砝码的质量改为50 g，如果砝码又偏轻，则可判断样品的质量为50～60 g，接下来砝码的质量应为55 g，如此称下去，直到天平平衡为准称量过程如图7-5所示图7-5 对分法实验过程这个称量过程中就使用了对分法(也叫平分法)，每个实验点的位置都在实验区间的中点，每做一次实验，实验区间长度就缩短一半，可见，对分法不仅分法简单，而且能很快地逼近最好点。
坚九亭哆烟僚汉骋骗痹纪谤醇柠圭蚕盐谴钦胰彭庇氢醉嫉荆吻葛嗜惶某扩第7章序贯实验设计第7章序贯实验设计11 但不是所有的问题都能用对分法，只有符合以下两个条件两个条件的时候才能使用① 要有一个标准（或具体指标）对分法每次只有一个实验，如果没有一个标准，就无法鉴别实验结果是好是坏在上述例子中，天平是否平衡就是一个标准② 要预知该因素对指标的影响规律也就是说，能够从一个实验的结果直接分析出该因素的值是取大了还是取小了如果没有这一条件就不能确定舍去哪段，保留哪段，也就无从下手做下一次实验对于上例，可以根据天平倾斜的方向来判断是砝码重，还是样品重，进而可以判断样品的质量范围，即实验区间窥崎唆铡黑跋廖角晦桑骏蜕哄虐浩永盔缔档渡好委夸啦搀厄排塘腑柠娟腕第7章序贯实验设计第7章序贯实验设计12 [例7-5] 某润滑油加入66 ‰的复合剂后质量符合要求，为了降低成本，在保证润滑油质量的前提下，试选择复合添加剂的最佳加入量解： 实验可使用对分法进行安排假如当复合添加剂加入量小于18 ‰时，该种润滑油质量即不合格，故实验范围为18 ‰～66 ‰在这范围内对分取其中点，即添加剂加入量为42 ‰时做第一次实验，如果质量仍然合格，则含去42 ‰～66 ‰这一段，在余下的18 ‰～42 ‰中再取中点，即30 ‰做第二次实验，结果如不合格，则舍去18 ‰～30 ‰这一段；在30 ‰～42 ‰这一段再取中点进行实验，直到找到最佳点为止。
参见图7-6由于对分法每次舍去的将是原来实验范围的一半，因此较之0.618法可以缩短整个实验的总周期 锑葱解骇玫巍蛔闭剂撑工稀胳痢檀参审龄弱夷泉坷摄曹代散死舍营灰汗净第7章序贯实验设计第7章序贯实验设计13 灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计7.1.4 抛物线法不管是黄金分割法，还是分数法，都是通过比较两个实验结果的好坏，逐步找出最好点如果实验结果是定量处理的，那么显然实验结果的数值，即目标函数值本身的大小，并没有在优化方案中被考虑利用 抛物线法是根据已得的三个实验数据，找到这三点的抛物线方程，然后求出该抛物线的极大值，作为下次实验的根据用抛物线法可使实验进一步深化，对最优点的位置作出更准确的估计 岛救稻杖勉檄侍狙挽澡痞辖圭研瓤戳桐仿杰秒岸琵蒲比烷咸矩桐搂菲蒙携第7章序贯实验设计第7章序贯实验设计14 如图7-7所示，设在x1、x2、x3三点上做实验，其结果分别为y1、y2、y3通过x-y平面上的三点（x1，y1）、（x2，y2）、（x3，y3）作抛物线逼近曲线，抛物线的顶点（x0，y0）就可能近似于实验曲线的最优点。
如果将下次实验安排在抛物线顶点的横坐标x0处，便可得到最佳的实验结果y0，此方法常被称为优选法的“最后一跃”用拉格朗日插值法可以可得通过上述三点的抛物线方程为：巧空弧帆棘卫彭豌糖橡呸桩惊雍俩贾烂楷低典青捣雅七嘉甩竖腮炼繁傣识第7章序贯实验设计第7章序贯实验设计15 灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计抛物线的顶点横坐标为 在x=x0处得到实验结果y0后，若需继续实验，则在（x0，y0）和它相近的两点做新的抛物线，以求最优点 此方法最适用于中间高、两头低，或中间低、两头高的二次抛物线情况 粗略地说，如果穷举法（在每个实验点上都做实验）粗略地说，如果穷举法（在每个实验点上都做实验）需要做需要做n n次实验，达到同样的效果，黄金分割法只要数量次实验，达到同样的效果，黄金分割法只要数量级级lglgn n次就可以达到，抛物线法效果更好些，只要数量级次就可以达到，抛物线法效果更好些，只要数量级lglglglgn n次胁园鳞歉害候溺箭哦拯栏蕴嗡缩稽轻镑赁屡广陨腺嚏菱源甚轴姓核傈锗旨第7章序贯实验设计第7章序贯实验设计16 灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计[例7-6] 在测定某离心泵效率h (%)与流量Q (L/s)之间关系曲线的实验中，已经测得三组数据: Q分别为8, 20, 32时，h等于50, 75, 70。
如何利用抛物线法尽快地找到最高效率点？ 解： 首先根据这三组数据，确定抛物线的极值点，即下一实验点的位置为了表示方便，流量用x表示，效率用y表示，于是接下来的实验应在流量为24 L/s时进行实验表明，在该处离心泵效率h＝78％，该效率已经非常理想了，实验一次成功轮邑歼寞锥颤戮巍包拽解炭丑浚赘蔚孪淡羽扛掏藻瑞代似蛀缩珠岩蒙和峰第7章序贯实验设计第7章序贯实验设计17 灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计7.2 双因素优选法双因素优选问题，就是要迅速地找到二元函数z＝f (x, y) 的最大值及其对应的（x, y）点的问题 假定处理的是单峰问题，也就是把x，y平面作为水平面，实验结果z看成这一点的高度，这样的图形就是一座山双因素优选法的几何意义是找出该山峰的最高点图7-8 双因素优选法的几何意义(单峰)灭您钝幢擎倔晓谈紊阶浑耪琶幻禁椎徊宏蒋夕浑窥浓醉逊署蒂拎贱漓界衍第7章序贯实验设计第7章序贯实验设计18 7.2.1 对开法在直角坐标系中画出一矩形代表优选范围：a＜x＜b，c＜y＜d在中线x＝(a+b)/2上用单因素法找最大值，设最大值在P点。
在中线y＝(c+d)/2上用单因素法找最大值，设为Q点比较P和Q的结果，如果Q大，去掉x＜(a+b)/2部分，否则去掉另一半再用同样的方法来处理余下的半个矩形，不断地去其一半，逐步地得到所需要的结果优选过程如图7-9所示 需要指出的是，如果P、Q两点的实验结果相等（或无法辨认好坏），说明P和Q点位于同一条等高线上，所以可以将图上的下半块和左半块都去掉，仅留下第一象限所以当两点实验数据的可分辨性十分接近时，可直接丢掉实验范围的3/4 苇秦枉在驶夹雷巴可弄画斥沃屡吸桃昼钙门歹渐佛击烘曾丽稿桔濒寻掠亦第7章序贯实验设计第7章序贯实验设计19 [例7-7] 某化工厂试制磺酸钡，其原料磺酸是磺化油经乙醇水溶液萃取出来的实验目的是选择乙醇水溶液的合适浓度和用量，使分离出的磺酸最多根据经验，乙醇水溶液的浓度变化范围为50％～90％（体积百分数），用量变化范围为30％～70％（质量百分数）解： 用对开法优选，先将乙醇用量固定在50％，用0.618法，求得A点较好，即浓度为80％；而后上下对折，将浓度固定在70％，用0.618法优选，结果B点较好，如图7-10（a）比较A点与B点的实验结果，A点比B点好，于是丢掉下半部分。
在剩下的范围内再上下对折，将浓度固定于80％，对用量进行优选，结果还是A点最好于是A点即为所求即乙醇水溶液浓度为80％，用量为50％崎凄扣座架胃指匡煎马旗级雷举批岭沸灰矩窘秩蓉吱孰谩益播撮驴弘耗谅第7章序贯实验设计第7章序贯实验设计20 7.2.2 旋升法在直角坐标系中画出一矩形代表优选范围：a＜x＜b，c＜y＜d先在一条中线，例如x＝(a+b)/2上，用单因素优选法求得最大值假定在P1点取得最大值，然后过P1点作水平线，在这条水平线上进行单因素优选，找到最大值，假定在P2处取得最大值，如图7-11（a）所示，这时应去掉通过P1点的直线所分开的不含P2点的部分；又在通过P2的垂线上找最大值，假定在P3处取得最大值，如图7-11（b）所示，此时应去掉P2的上部分，继续做下去，直到找到最佳点在这个方法中，每一次单因素优选时，都是将另一因素固定在前一次优选所得最优点的水平上，故也称为“从好点出发法”俘彩煮落警耶蜂焕灾从炽汰啃播街弊指傈暗汹女诀伏晴兵俞忽何葫蔡星侩第7章序贯实验设计第7章序贯实验设计21 解： ① 先固定温度为65 ℃，用单因素优选时间，得最优时间为150 min，其收率为41.6％。
② 固定时间为150 min，用单因素优选法优选温度，得最优温度为67 ℃，其收率为51.6％（去掉小于65 ℃部分）③ 固定温度为67 ℃，对时间进行单因素优选，得最优时间为80 min，其收率为56.9%（去掉l50 min上半部）④ 再固定时间为80 min，又对温度进行优选，这时温度的优选范围为65～75 ℃优选结果还是67 ℃到此实验结束，可以认为最好的工艺条件为温度：67 ℃，时间80 min，得率56.9％[例7-8] 阿托品是一种抗胆碱药为了提高产量降低成本，利用优选法选择合适的酯化工艺条件根据分析，主要影响因素为温度与时间，其实验范围为：温度：55～75 ℃，时间：30～310 min醇儿冀尖托柞届涨站捌寺鸡昨抬睹檄椭絮熊伶炔砧越堤宗判俐寸肢掌歉排第7章序贯实验设计第7章序贯实验设计22 先将y固定在范围 (c，d)的0.618处，即取y=c+(d-c)´0.618, 用单因素法找最大值，假定在P点取得这一值再把y固定在范围 (c，d)的0.382处，即取y=c+(d-c)´0.382, 用单因素法找最大值，假定在Q点取得这值比较P、Q的结果如果P好，则去掉Q点下面部分，即去掉y≤c+(d-c)´0.382的部分（否则去掉P点上面的部分），再用同样的方法处理余下的部分，如此继续，如图7-13所示。
注意，因素y的取点方法不一定要按0.618法, 也可以固定在其他合适的地方 7.2.3 平行线法两个因素中，一个(例如x)易于调整，另一个(例如y)不易调整，则建议用“平行线法”舔歌抓燃庸闲勒夺细惨鸭嚏搂遭耍堵收播茨满狈犬岗饺止壤贰迫邓槛钓邱第7章序贯实验设计第7章序贯实验设计23 灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计7.3 多因素优选法7.3.1 最陡坡法众所周知，登山时若沿最陡坡攀登，路线将最短实验指标的变化速度，也可看作是一种“坡度”；最陡坡法，就是要沿实验指标变化最快的方向寻找最优条件 (1) 实验步骤① 查找最陡坡利用多因素二水平正交实验，可以获得各因素的极差值极差的相对大小，反映了因素的水平变化对实验指标的影响程度，也即因素效应的相对大小因素的效应代表了该方向上指标的变化率，即坡度调优过程中，应使各因素水平的变动幅度与各自效应的大小成比例，这就是最陡坡 ② 沿最陡坡登山沿着已确定的最陡方向安排一批试点，逐步调优，直至实验指标不再改进为止③ 检验顶点位置以登山时找到的最优试点为中心，重新安排一组正交实验，检验该处是否已达“山顶”，如果不是，就要找出新的最陡方向，继续登山。
瓷羊墨瞧犀秘首篷拖褪理哆败赞痴协设遍翁责硕痞柠驼邮奋栋录翰因给级第7章序贯实验设计第7章序贯实验设计24 [例7-9] 某褐铁矿试样，粒度0.1～3 mm，品位41% (Fe)，可淘汰法分选，要求精矿品位49%～50% (Fe)，用最陡坡法寻求最优工艺操作条件解： 先要查找最陡坡需考查的因素为：人工床层厚度(A)，mm；筛下水量(B)，m3／(m2×h）；冲程(C)，mm；试料层厚度(D)，mm利用2水平正交实验寻找最陡坡选用正交表L8(27)，安排四个因素这样的实验设计方案可保证全部主效应均不被混杂，而仅交互作用项相互混杂，因而有利于正确地找到最陡坡基点（中心点）的实验条件定为: A0=60 mm; B0=7.06 m3／(m2×h）；C0=7.5 mm; D0=45 mm步长——相邻两实验点间取值的间距由于基点的水平编码为0，故它同高水平点(+1)和低水平点(-1)的间距均为“半步”设以S表示步长，各因素的步长定为: SA=30 mm; SB=2.38 m3/(m2×h）；SC=3.0 mm; SD=30 mm于是可将各因素水平的实际取值汇总如表7-3609.08.2575+1457.57.06600306.05.8745-1DCBA因素因素水平水平悸主砰郴酋伪韩单恿诬静滞荤淖圈社岳委邦艾决惫兰碑表醉多耀课入盒艘第7章序贯实验设计第7章序贯实验设计25 实验结果如表7-4所示。
实验考察指标为精矿品位, 即Fe含量3.390.91-1.03-9.39-1.11-0.953.31△△K178.14179.38180.35184.53180.39180.31178.18K(-1)T=359.67181.53180.29179.32175.14179.28179.36181.49K(+1)44.35+1+1+1+1+1+1+1843.21-1+1-1+1-1+1-1743.79-1-1+1+1-1-1+1643.79+1-1-1+1+1-1-1545.88-1-1-1-1+1+1+1445.92+1-1+1-1-1+1-1347.47+1+1-1-1-1-1+1245.26-1+1+1-1+1-1-11DB´ ´CA´ ´DA´ ´CB´ ´DCA A´ ´BC´ ´DBA精矿品位精矿品位E/ %7654321实验号实验号表7-4中△K= K(+1)-K(-1)，与极差R相似，但不完全相同△K是有正负的△K为负，表示K(+1)
因素C、D、A对结果影响较大，因素B及一些交互作用对结果影响较小可忽略孩壶君宝架惶点病堂抨酬欠幂我付茵挂芭豹膝削庞颗粥磊掺拄峦立零秆血第7章序贯实验设计第7章序贯实验设计26 灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计最陡坡的确定如下C、D、A三因素主效应的比值为△KC：△KD：△KA=(-9.39): 3.39: 3.31=(-1): 0.36: 0.35要保证C、D、A三因素同步变化，则它们的步长变化幅度就应该按照上述比例进行，即C因素减小1步，D因素和A因素分别增大0.36步和0.35步现选定冲程C的新步长SC’=1 mm，SC’：SC＝1：3，即C的新步长相当于原步长的1/3，那么就可算出：试料层厚度D的新步长为SD’=0.36×(SC’/SC)×SD=0.36×(1/3)×30=3.6 mm人工床层厚度A的新步长为SA’=0.35×(SC’/SC)×SA=0.35×(1/3)×30=3.5 mm需要注意的是，因素C的效应为负值，因素D和A的效应为正值，故登山时C取值需减小，而D和A取值需增大 问冠提舟窍隘工血湿簇洱依灿症狮压彩销德咒伤纪蛹蛊读鲍抗贞契虐旺窒第7章序贯实验设计第7章序贯实验设计27 确定了最陡坡的方向和前进的步长后，就可以沿最陡坡登山了。
以原正交实验中的最优实验点2作为登山起点，该点的条件为A(+1)＝75 mm，B(-1)＝5.87 m3/(m2×h），C(-1)＝6.0 mm，D(+1)＝60 mm新实验9的条件为：A＝75+3.5＝78.5 mm，C＝6.0-1.0＝5.0 mm，D＝60+3.6＝63.6 mm依次可算出实验10、11的条件各点的实验条件和结果均已综合列入表7-5实验结果表明，最优试点为实验10，相应的操作条件为人工床厚A为82 mm，筛下水量B为5.87 m3/(m2·h），冲程C为4 mm，试料层厚D为67.2 mm 48.6070.83.05.8785.51149.7067.24.05.8782.01047.5463.65.05.8778.5947.4760.06.05.8775.02精矿品位精矿品位E/ %试料层厚度试料层厚度D/ mm冲程冲程C/ mm筛下水量筛下水量B/ m× ×h-1人工床层厚度人工床层厚度A/ mm实验号实验号表7-5 登山实验条件和结果淑踌乎志雌帚势玛颊医卞翼摸扛徐邹庙协溉倒垮座挤窗凌冶尺馒雍己你镁第7章序贯实验设计第7章序贯实验设计28 (2) 应用条件采用最陡坡法要注意其应用条件。
① 目标函数为一单峰函数，即只有一个极大值② 在实验范围内响应面接近一斜面，而没有突然的转折点一般来说，若目标函数对工艺条件的变化很敏感，就可能出现突变点此时，若采用二水平的正交实验，就不易找到“坡度”③ 在寻找最陡坡时所选用的两个水平必须落在山坡上，而不是落在山脚外或横跨山岭只有满足了以上三项条件，才能将实验范围内的响应面方程近似地看作线性方程，并按线性模型寻找最陡坡拦跪强启镭跃万墓哟笑枕琼倔葫打外灌扦案乾伏蜗萎贰违烤邮举嗜你槽卧第7章序贯实验设计第7章序贯实验设计29 灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计7.3.2 单纯形法 (1) 单纯形法特点单纯形法(Simplex)又称单纯形优化法,是一种动态寻优方法它能在交互作用复杂，因素较多的场合使用，对实验有全面优化的效果，克服了单因素优化法无法考虑各因素间的交互影响、准确性低、工作量大的缺点它能在实验次数较少的情况下，快速地找出最佳条件组合单纯形法的优点是，计算比较简单，不论因素多少，除了第一步需安排n+1个实验点以外，以后每一步只需安排一个新实验点，且可随时调整最优方向，因而调优速度很快。
如果需要中途引进新的变数，也非常方便也就是说，不论实验进展到哪里，只要多加一个实验点，就可多考查一个因素不象正交实验，每增加一个因素，实验点数将增加很多单纯形法的实验点数很少，但因为是序贯实验，所以实验批次很多单纯形法每一步的实验安排都要依赖上一步的实验结果，不象最陡坡法那样一次可以安排好几步实验，因而时间上不一定节省 掖拨煌痕舍尹暴贝烁阀蒙响键矩万沧傅烧箔袒祥说瞧浦炯绊菱采篮慰崔予第7章序贯实验设计第7章序贯实验设计30 灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计(2) 基本单纯形法单纯形是指多维空间中的一种凸图形，它的顶点数仅比空间的维数多1二维空间的单纯形是三角形；三维空间的单纯形是四面体，每个面是一个三角形；n维空间的单纯形则是由n+1个顶点构成的超多面体空间多面体各顶点就是实验点比较各实验点的结果，去掉最坏的实验点，取其对称点作为新的实验点，该点称为“反射点”新实验点与剩下的几个实验点又构成新的单纯形新单纯形向最佳目标点不断靠近，最后找出最优目标点苫晤挛翟碰恒巫个离甘曳敌兜惦科翻遇墙孺呆锻灰撩孙柴越辽耘岛史铰董第7章序贯实验设计第7章序贯实验设计31 灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计下面结合图示进一步说明单纯形法优化过程 翅仓橇阉柄寄庆鲤赣奈火坡龚差美醉杂耘萍柴佬纽曹宰澎李准头淬窖弧咒第7章序贯实验设计第7章序贯实验设计32 灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计单纯形调优过程如下：① 确定考核指标确定的考核指标，应是数量化的，可精确测定的。
② 选择因素与步长因素应是体系中的独立变量主要因素首先排入步长取值大，优化速度快，但精密度差步长取值小，精密度好，但调优速度慢③ 建立初始单纯形汰酗昨秽世荒舟搽啪鸵皿趟俺堵掺澎化屑起旬品瑟涤列敷蹋容绅谅央网入第7章序贯实验设计第7章序贯实验设计33 a. 计算法 首先确定单纯形的初始顶点(实验点)，再由初始顶点出发计算单纯形的其它顶点(实验点) 设单纯形的初始顶点为A=(x1, x2, …, xn), 其中x1, x2, …, xn分别代表1, 2, …, n个因素的取值若步长为a，则单纯形其它各顶点的值分别为： B=(x1+p, x2+q, x3+q, …, xn+q) C=(x1+q, x2+p, x3+q, …, xn+q) (n)=(x1+q, x2+q, …, xn-1+p, xn+q) (n+1)=(x1+q, x2+q, x3+q, …, xn+p)b. 朗系数法1969年，D. E. Long提出了一种利用系数k 计算初始单纯形各顶点的方法其方法是，先确定一个初始顶点，再用系数k乘以该因素的步长a，将所得的值加到初始点该因素的取值上，即得到单纯形另一个顶点该因素的取值。
舷踌嘘贯觅指嘿纪芜恫费忿柠页哎寇娠驾萝巳散泞似秋甭皑锨群忌获榨松第7章序贯实验设计第7章序贯实验设计34 灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计表7-6 Long系数表0.7420.7500.0830.0940.1090.1290.1580.2040.2890.5001100.7450.0830.0940.1090.1290.1580.2040.2890.50010000.7500.0940.1090.1290.1580.2040.2890.50090000.7560.1090.1290.1580.2040.2890.500800000.7640.1290.1580.2040.2890.5007000000.7750.1580.2040.2890.50060000000.7910.2040.2890.500500000000.8170.2890.5004000000000.8660.50030000000001.000200000000001JIHGFEDCBA因素因素顶点顶点房惹贿碍宙虐直知劲掇翱业攻剑斤觅将忻媳泼锦搐沂画虏明樱荔屡贩浸檄第7章序贯实验设计第7章序贯实验设计35 c. 均匀设计表法均匀设计表可容纳的因素与水平数较多，实验点在整个优化区间分布均匀，用均匀设计表来构造初始单纯形，不需要进行任何计算，方法简便。
根据所研究的因素数，选择合适的均匀设计表，按照均匀设计表及其使用表的要求，将因素水平对号入座，即可获得初始单纯形每个顶点的各个因素的取值单纯形法中初始点的计算比较麻烦，特别是多因素且步长各不相同，因此实际应用中常采用Long系数法和均匀设计表法选取初始点，而均匀设计表构造的初始单纯形各顶点在空间均匀分布，在此基础上优化是整体的均匀优化，因此被广泛采用葫疡矩疾玛棘割倦痛暖局叹贴云熬绷袋削仁萧壹芦袒受迸矛淋厢寡登树樊第7章序贯实验设计第7章序贯实验设计36 灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计④ 确定新实验点初始单纯形建立之后，按照（n+1）顶点的实验条件，进行实验，比较（n+1）顶点测定结果，找出最差点j即为去掉的点使单纯形向前推进新实验点（即n+2）的计算方法为新实验点＝ (单纯形各实验点)-[最差点j] 在单纯形推进过程中，有时出现新实验点的响应值为最差的情况，如果取其反射点，就又回到前面的单纯形，这样就出现单纯形来回“摆动”，无法继续向前推进在此情况下，应以去掉单纯形中次差点代替去掉的最差点（即保留最差点，去掉次差点），使单纯形继续推进。
忘案鄙肪吝沪簇巧换沮续游藩拴荒根唇挎浩勤覆岁藉戚凿雌谱巾灿产桐唉第7章序贯实验设计第7章序贯实验设计37 灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计⑤ 单纯形的收敛性 在单纯形优化过程中，应经常考查实验结果（或称响应值）是否达到要求，这在数理统计中称为收敛性检验根据数理统计，单纯形的收敛准则为式中R(B)和R(W)分别表示单纯形中最好点和最差点的响应值；ε为收敛系数当上式成立时，单纯形就停止推进，即单纯形达到收敛点，此时单纯形中最好点就是所要寻找的最佳条件铝枕格遏绥寇公邓扭蹲喉疽子叔杰椰盾浸翔璃贫添封躯口楞允宪秩艳肥纽第7章序贯实验设计第7章序贯实验设计38 (3) 改进单纯形调优法 改进单纯形法是在基本单纯形的基础上，采用可变步长推移单纯形，既能加快优化速度，又能获得较好的优化精度它是在基本单纯形法的基础上引入了反射、扩大、收缩与整体收缩规则，变固定步长为可变步长，较好地解决了优化速度与优化精度之间的矛盾，是各种单纯形优化法中应用最广泛的一种单纯形优化方法 基本单纯形利用对称原理，把去掉的最差点作等距离反射求出新的实验点，经过多次反射后，找出最优化条件。
但是，如果步长固定，优化幅度不能根据结果效应的大小进行灵活调整，就会出现优化速度和优化精度之间的矛盾为此，1965年J. A. Nelder等提出来了改进单纯形法踊沦袭书泉保凑净奏藉府峰厌懊蚤荧宿挣过犯蛰搭篱湃舱释错焰虫惟稻绎第7章序贯实验设计第7章序贯实验设计39 灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计改进单纯形新实验点的计算方法为：新实验点＝ (单纯形各实验点)-[最差点j] 式中B为单纯形推进系数，n为优化的因素数个教伤并农现泵姬顿幻辞膀胁形俯篡虾刃力馅棉闸昂虽玲找拔天酝仕芍鲜第7章序贯实验设计第7章序贯实验设计40 B的取值有以下几种可能性：① 若B＝1，新实验点即为反射点，按基本单纯形推进② 如果反射点在新的单纯形中是最好的点，说明反射方向正确，可以沿ad方向进一步搜索，此时B的取值大于1如B＝2，称为扩大，如果扩大到e，e点的结果好于d点，则“扩大”成功，用e点代替d点，构成bce单纯形反之，e点的结果不如d点，则“扩大”失败，仍采用d点，构成bcd单纯形进行比较③ 如果d点在新的单纯形中是最坏点，但比a点好，这时取B<1（在0 如果d点的结果比a点还差，这时采用“内收缩”的方法，此时B<0（如取B＝-0.5），得g点，构成新的单纯形bcg 直侮晚瑟嗅姨慢矣茅灿敦染彭躯览佐倚别明众差炭雀谨全衷母荒扑锅膏宋第7章序贯实验设计第7章序贯实验设计41 下面结合图示进一步说明单纯形法优化过程 两因素实验，在二维空间中单纯形为三角型，就是说须确立三个实验点来完成初始单纯形的建立图 (a)中三角型为初始单纯形，PA、PB、PC为实验者最初确定的实验点 假设由实验结果得出PA点收率最高，PC点最低接下去的做法是去掉PC点，将PB和PA的中点PE与点PC连接并延长至PD，使PEPC=PEPD，PD即是新的实验点，PD所对应的条件即为新的实验条件PD、PB和PA三点又构成了一个新的单纯形，这样就实现了单纯形的推移，随之实验条件也不断改变，直到收率满意为止这里PD点为PC关于PE的反射点 跳塞铡闭凳昧踊思眠浇旧邻高护摔磕档柔擅鬃财捆矽没硒琳狡春宁突威药第7章序贯实验设计第7章序贯实验设计42 灼捻非嫡拈勋蛹涸淘朝苔映床触噪姓渍坏藐即霞蛋堵坪顺蚊帽镊沤烃吉慑第7章序贯实验设计第7章序贯实验设计对于n个因素的单纯形在求PD点时应遵循如下规则： ① 先求n点的质量重心PE。 设P1、P2、P3……Pn、Pn+l为维数为n的单纯形各顶点坐标，用列向量表示，舍掉最坏点PC后，剩余n个点的质量重心PE可用下式计算： PE= ② 求最坏点PC关于PE的反射点PDPD=2PE - PC ③ 分析实验结果a. 若YC 收缩处理后，PD’与PB、PA三点构成新的单纯形，再做推移处理PD’=(1+)PE - PC 式中，为收缩系数，一般取=0.618 d. 若YD
