正交实验设计与数据处理.pptx

《正交实验设计与数据处理.pptx》由会员分享，可在线阅读，更多相关《正交实验设计与数据处理.pptx（70页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、试验优化设计Experimental design以概率论和数理统计为理论基础经济地、科学地安排试验正确地分析试验数据尽快地获得优化方案一项实用技术第1页/共70页试验的必要性生产实践中，试制新产品、改革新工艺、寻求好的生产条件等，这些都需要先做试验试验需要花费时间、消耗人力、物力人们希望做试验的次数尽量少，而得到的结果尽可能好 第2页/共70页完整的试验应包括三个方面：第一，试验的设计第二，试验的实施第三，试验结果的分析第3页/共70页试验设计的概念即试验的最优化设计，是研究如何合理而有效地获得数据资料的方法。其主要内容就是讨论如何合理地安排试验、取得数据，然后进行综合的科学分析，以达到尽快

2、获得最优方案的目的。是数理统计学的应用方法之一。第4页/共70页实际问题很复杂，对试验有影响的因素往往是多方面的。在多因素、多水平试验中，若对每个因素的每个水平都互相搭配进行全面试验，需要做的试验就会很多。第5页/共70页例如两个7水平的因素，如果两因素的各个水平都互相搭配进行全面试验，要做 7249次试验，3个7水平的因素进行全面试验，就要做 73343次试验6个7水平的因素进行全面试验，要做 76117649次试验第6页/共70页试验次数多，要花费大量的人力、物力，还要用相当长的时间由于时间长，条件改变，还会使试验失效长期的实践表明，要得理想的结果，并不需要做全面试验。尤其对那些试验费用很

3、高，或具有破坏性的试验，更不要做全面试验第7页/共70页应当在不影响试验效果的前提下，尽可能减少试验次数第8页/共70页试验设计的意义根据专业知识和统计学知识，制定出一个包含多因素多水平有机组合的试验计划；清楚地表达出每次试验涉及到各因素的哪个水平，各种水平组合下是否需要做重复试验，全部试验共需进行多少次；是试验的行动指南；使试验过程有据可查、有章可循，并能有效地、科学地解决试验研究中的复杂问题。第9页/共70页试验设计的基本步骤提出拟考察的试验因素及水平考虑现有的人力、物力和试验条件，有能力承担多大规模的试验研究说明需要观测哪些指标，是定量的还是定性的或两种兼而有之选择合适的设计类型并给出全

4、部因素水平组合的方案第10页/共70页正交试验设计是解决这类问题的有效方法之一；正交试验设计的主要工具是正交表，用正交表安排试验是一种较好的方法，在实践中已得到广泛的应用正交试验设计是用一套规格化的表格来安排试验，这种表格叫做正交表第11页/共70页正交表简介是一种特制的表格，每个表都有一个记号，如L9(34)，L8(27)，就是两个最常用的正交表；符号说明：L正交表 L下角的9、8正交表的行数 括号里的3、2因素所取的水平数，指数4、7正交表的列数表内的数字1、2、3因素的水平第12页/共70页二水平的正交表还有L16(215)、L12(211)，三水平的正交表还有L18(37)，L27(3

5、13)，四水平的正交表还有L16(45)等等。第13页/共70页正交表L9(34)试验号 列号1234111112122231333421235223162312731328321393321第14页/共70页正交表记法一般正交表记为Ln(mk)，n是表的行数，是要安排的试验次数；k表中列数，表示因素的个数；m是各因素的水平数。第15页/共70页Ln(Km)是一个n行m列矩阵（1）每个因素 ；（2）任一列中每个不同元素出现的次数都等于 ；（3）H中 任 意 二 列 同 行 上 的 有 序 数 偶：（1,1）（1,k），（2,1）（k,1）（k,k），出现的次数都等于n/k2。第16页/共70页

6、正交表的特点（或性质）-正交性每列中不同数字出现的次数是相等的，如L9(34)中不同的数字是1,2,3，它们各出现3次。(n/k)在任意两列中，将同一行的两个数字看成有序数对出现的次数也是相等的，如L9(34)，有序数对共有9个：(1,1),(1,2),(1,3),(2,1),(2,2),(2,3),(3,1),(3,2),(3,3)，它们各出现一次。(n/k2)第17页/共70页正交性的优点保证每个因素的各个水平参加试验的次数一样多，每两个因素之间实现了完全水平组合试验因素的各种水平的搭配是均衡的正交表所需要试验的次数在多因素多水平的完全组合试验中具有较好的代表性正交试验中的最优水平组合在完

7、全水平组合试验中仍是名列前茅第18页/共70页混合水平正交表水平数不等的正交表，称为混合水平正交表如L8(4124)是由一个4水平的列，4个2水平的列组成，表示用该表设计试验时最多可安排一个4水平的因素，4个2水平的因素，需要试验的总次数为8次其它如L18(2137)，L32(814626)等等，都有类似的含义。第19页/共70页交互作用表需要考虑因素的交互作用时，许多正交表都配有一张交互作用表从交互作用表上可以查出相应正交表中任何两列之间的交互作用所占的列比如对L8(27)，由相应的交互作用表可以查到L8(27)中第1,2两列的交互作用应占第3列，而第1,4两列的交互作用应占第5列，第2,4

8、两列的交互作用应占第6列等等。而三水平正交表中任意两列的交互作用都占有两列。第20页/共70页正交试验设计应用范围常常用来解决二水平或三水平或二、三水混合水平的多因素设计问题；适用于需要考察的交互作用不多、也不太复杂的多因素试验研究的场合；通过方差分析鉴别各因素对试验指标的影响。第21页/共70页正交试验设计步骤首先要根据试验目的，确定要观察的因素确定每个因素的水平然后选用适当的正交表安排试验。第22页/共70页例一为提高某化工酸洗过程中产品的收率，选择了四个有关因素：反应温度（A）、反应时间（B），用酸量（C）和酸浓度(D)，由生产实践经验及专业知识选取的水平如下：表1 因素水平表水平 反应

9、温度（A）反应时间（B）加酸量（C）酸的浓度（D）1 85 90 1:1 10 2 90 120 2:1 15 3 95 150 3:1 20 第23页/共70页用正交表安排试验的步骤：首先选择适合试验的正交表，这里是四因素三水平试验，用L9(34)比较适合，这时只要做9次试验，每一行代表一次试验的条件。把A、B、C、D四个因素放到L9(34)的四列上，每一行所代表的试验条件，即把A、B、C、D对应的“1”“2”“3”翻译成具体的水平，如表2所示，并将试验结果填写入相应的试验序号中，然后进行计算和分析。第24页/共70页表2 试验方案与结果列号试验号1A2B3C4D列号试验号1A2B3C4D1

10、111111(85)1(90)1(1:1)1(10)21222212(120)2(2:1)2(15)31333313(150)3(3:1)3(20)4212342(90)123522315223162312623127313273(95)13283213832139332193321第25页/共70页正交实验结果分析对试验结果的分析：直观分析 方差分析第26页/共70页表2 正交试验直观分析计算结果11111A1B1C1D12221222A1B2C2D25231333A1B3C3D34342123A2B1C2D35852231A2B2C3D16162312A2B3C1D26173132A3B1

11、C3D25583213A3B2C1D37093321A3B3C2D164K1117135153147K2180183174168K3189168159171k139455149k260605856k363565357R241578第27页/共70页说明与结果分析：K为各因素相同水平的指标之和；极差R为相同因素不同水平的指标平均值中最大值与最小值之差；并作出因子水平与R值间关系图：由图可以看出因子与指标的变化规律，明确因子影响指标的主次：ABCD1 2 3 1 2 3 1 2 3 1 2 3 k A C B 因素 D 选定最优方案为A3B2C2D3，并验证该方案。第28页/共70页第29页/共7

12、0页第30页/共70页第31页/共70页第32页/共70页第33页/共70页第34页/共70页第35页/共70页第36页/共70页第37页/共70页第38页/共70页第39页/共70页第40页/共70页第41页/共70页第42页/共70页第43页/共70页第44页/共70页第45页/共70页第46页/共70页第47页/共70页第48页/共70页第49页/共70页第50页/共70页第51页/共70页第52页/共70页正交试验直观分析：看出因子与指标的变化规律；明确因子影响指标的主次：ABCD选定最优方案,取务因子中有最高指标的水平组合1 2 3 1 2 3 1 2 3 1 2 3 k A C B

13、 因素 D 选定最优方案为A3B2C2D3，并验证该方案。第53页/共70页正交试验的方差分析直观分析不能给出误差的大小，也不能知道结果的精度；方差分析能够反应数据的波动性，即数据的分散性，方差大小表明数据变化的显著程度，也表明因素对指标影响大小。第54页/共70页例二为提高某化工产品的转化率，选择 了三个有关因素：反应温度（A）、反应时间（B），用碱量（C），选取的水平如下：因素水平 1 2 3 温度（A）80 85 90 时间（B）90分 120分 150分 用碱量（C）5%6%7%第55页/共70页用正交表安排试验的步骤：首先选择适合试验的正交表，这里是三因素三水平试验，用L9(34)比

14、较适合，这时只要做9次试验把A、B、C三个因素放到L9(34)的任意三列上，例如放在前三列，把A、B、C对应的“1”“2”“3”翻译成具体的水平，如表2所示：第56页/共70页表2 试验方案与结果列号试验号1A2B3C4列号试验号1A2B3C指标（转化率）1111111(80)1(90)1(5)31(%)2122221(80)2(120)2(6)543133331(80)3(150)3(7)384212342(85)1(90)2(6)535223152(85)2(120)3(7)496231262(85)3(150)1(5)427313273(90)1(90)3(7)578321383(90)

15、2(120)1(5)629332193(90)3(150)2(6)64第57页/共70页对试验结果进行方差分析设9次试验结果以y1,y2,y9表示，我们假定yr r+r (r=1,2,9)r的结构式为 y1+a1+b1+c1+1 y2+a2+b2+c2+2 y9+a9+b9+c9+9 r N(0,2)r=1,2,9（1）第58页/共70页其中r(r=1,2,9)相互独立，即各次试验结果理论值r 的总平均，ai,bj,ck(i,j,k=1,2,3)分别表示因素A,B,C对指标的效应，即 ai=（水平Ai下各次试验结果值r的平均）i=1,2,3 bj=（水平Bj下各次试验结果值r的平均）j=1,2

16、,3 ck=（水平Ck下各次试验结果值r的平均）k=1,2,3第59页/共70页且满足要检验的假设为：H01：a1=a2=a3=0 H02：b1=b2=b3=0 H03：c1=c2=c3 =0 第60页/共70页令其中于是9次试验结果的总偏差平方和为：经过初等的代数运算，可以得到如下的分解式：第61页/共70页其中（2）第62页/共70页SA反映了因素A各水平效应引起的差异，它正好等于正交表L9(34)中第一列各水平的偏差平方和S1；SB反映了因素B各水平效应引起的差异，它正好等于正交表L9(34)中第二列各水平的偏差平方和S2；SC反映了因素C各水平效应引起的差异，它正好等于正交表L9(34

17、)中第三列各水平的偏差平方和S3；第63页/共70页 由此可见，正交表L9(34)的总偏差平方和恰好等于各列偏差平方和之和：ST=S1+S2+S3+S4 相应自由度也有fT=f1+f2+f3+f4第64页/共70页对（3）第65页/共70页按上页的公式，对表2中数据计算可得：ST=982 SA=618 SB=114 SC=234SE=18fA=fB=fC=fE=3-1=2ST 自由度fT=（总试验次数1）8在原假设：“第i列中每个水平的效应为0”成立时，有统计量 i=1,2,3第66页/共70页例题的方差分析表如下：方差来源平方和自由度均方Fi临界值A618230934.33*F0.05(2,

18、2)B1142576.33=19C234211713.00误差1829总和9848第67页/共70页作显著性检验时，常常取 0.01,0.05,0.10,0.20 定出四个临界值，则该因素对指标影响特别显著 ，则该因素对指标影响显著 ，则该因素对指标影响较为显著 ，则该因素对指标有一定影响 ，则该因素对指标影响不显著在本例中，F0.01=99,F0.1=9,F0.2=4，则各因素影响指标的主次顺序为：ACB第68页/共70页u在试验范围内，A,B,C的优水平组合为A3B2C2，即反应温度为90，反应时间为120分，用碱量6%u但是这个水平组合不在9次试验之中，它是否为优，还要通过试验检验，也可以根据数据结构式（1)对这个水平组合下，指标yr的理论值r作出估计。第69页/共70页感谢您的观看！第70页/共70页