试验设计与数据处理作业----333333

1、学生姓名陈柏娥试验设计与数据处理正交实验方差分析法确定优方案化工081班20084540104学院名称指导教师化学化工学院范明舫2011年04月20日实验设计与数据处理课程的收获与体会实验设计与数据处理课程具有公式多、计算多、图表多等特点，涉及较多概率论基 础知识， 课程本身的繁杂性决定了理解和掌握起来难度较大。 一开始的时候， 我还有点担心 这一门课会学不好， 因为我的概率论和数理统计的知识基础薄弱， 可能会对里面的内容产生 难以理解的心理， 有点感觉他是郁闷枯燥乏味的课程。 不过， 在老师的指导下我否认了之前 的观点。这门课的安排很合理，从简单到复杂，由浅入深的思维发展规律，现将单因素试验

2、、 双因素试验、 正交试验、 均匀实验设计等常用实验设计方法及常规数据处理方法、 理论、方差分析、回归分析等数据处理的理论知识、最后讲得出的方差分析、回归分析 结论和处理方法直接应用到实验设计方法。老师也让我们先熟悉实验设计方法， 数据处理方法，使我较早的感受到应用试验设计方法指导实践的“收获” 学习兴趣。通过学习，我初步认识了这一门课。这门课是研究如何合理而有效地获得数据资料的 方法。讨论如何合理安排实验、 取得数据、然后进行综合的科学分析， 从而达到尽快获得最 优方案的目的， 即实验的最优设计。 实验设计方法是数据统计学的应用方法之一。 一般的数 据统计方法主要是对已获得的数据资料尽可能精

3、确的判断。如果试验安排得好且分析得当， 就能以较少的试验次数、较短的试验时间、较低的费用，得到较满意的实验结果；反之，如 果试验安排的不得当，分析不得当，则试验次数增加，试验时间延长，浪费人力、物力、财 力，难以达到预期的结果，甚至导致实验失败。通过这门课程的学习，是我对误差理论、方 差分析、正交试验设计与应用、回归分析都有了一个很好的理解，并且将它们做了笔记。比如方差分析的理解：方差分析市实验设计中的重要分析方法，应用非常广泛，它是 将不同因素， 不同水平组合下的实验数据作为不同总体的样本数据， 进行统计分析， 找出对 实验结果影响大的因素及其影响程度。 对于单因素试验的数据进行统计分析，

4、找出对试验指 标影响大的因素及其影响程度。对于单因素试验的方差分析，主要步骤如下： 1，建立线性 统计模型，提出需要检验的假设。 2，总离差平方和的分析与计算。 3，统计分析，列出方差 分析表。对于双因素试验的方差分析，分为两种， 一种无交互作用的方差分析，另一种有交 互作用的方差分析，对于这两种类型分别有各自的设计方法， 但是总体步骤都和单因素试验 的方差分析一样。我们又通过正交试验设计合理安排实验，他是尽快有效的获得最优方案的一种设计方 法。了解了他是避免做全面试验， 再多因素多水平实验中选择最有代表性的搭配。 否则花费 时间过长，人力，物力，财力消耗太多。尤其是一些长周期、高费用或破坏性

5、试验，更不要 做全面性试验。我觉得学习了这门课我得到了很大的收获，特别是在分析方面。我学习到的科学方法 相对于以前高中学到的观点有了很大的进步，扩展了我认识的视野。而且从这一门课程中，再讲误差 等 并掌握常规 ，从而激发并维持最大的乐趣就是他很多都能应用到实践中。在实验的课程中能和同学们一起发现问题和解决问题，更加深了我对这门课的认识。1.为了通过正交试验寻找从某矿物中提取稀土元素的最优工艺条件，使稀土 元素提取率最高，选取的因素水平如下（见表1-1）表1-1因素表需要考虑的交互作用有A X B,A X C,BX C,如果将A,B,C分别安排在正交表L8（27）的1,2,4列上，试验结果（提取

6、量/ml）依次 为:1_1.01,1.13,1.13,1.06,1.03,0.80,0.76,0.56试用方差分析法（a =0.05）分析试验 结果，确定较优工艺条件解：列出正交表L8（27）和试验结果，见表1-2。表1-2例1试验结果及分析试验1234567提取量号ABAX BCAX CBX C空列yi/(ml)111111111.01211122221.3331221122106521212121.03621221210.8722112210.76822121120.56K4.5304.1703.6603.9303.5003.6603.630T=7.68K23.1

7、503.5104.0203.7504.1804.0204.050Q=7.7816(A)酸 用量/ml(B)水 用量/ml（C）反应时间/h125201220402:求K:K仁SUMIF(B$3:B$10,1,$l$3:$1$10),选中该公式，然后水平拖 动填充柄，就可计算出后六列的Ki值；K2=SUMIF(B$3 : B$10,2,$I$3:$I$10),选中该公式，然后水平拖动填充柄，就可计算出后六列的K2值。2求k:首先选中单元各区域B11:H12，在该区域的左上角第一个单元格即B11中或在编辑栏中输入：=B17:H18/4,然后在同时按“Shift+Ctrl+Enter”，即可 在B1

8、3:H14中显示k值。3求极差R:在B15中选中该单元格输入：=MAX(B11:B12)-MIN(B11:B12)，按下Enter键，然后选中该单元格，向右拖动填充柄就可计算出后六列的极 差Ro (也可用下式计算R的值：R=MAX(B13:B14)-MIN(B13:B14)(2)T=SUM(l3:l10)=7.68;Q=SUM(yi)2)=7.7816(i=1,2,3,4,5,6,7,8);P=T2/n=7.682/8=7.3728;计算力差平方和：总离差平方和:SST=Q-P=7.7816-7.3728=0.4088;因素与交互作用的离差平方和：SSA=SS1=(1/8)RA2=(1/n)(

9、K1-K2)2=(1/8)(4.53-3.15)2=0.23805;SS3=SS2=(1/8)RB2=(1/n)(K1-K2)2=(1/8)(4.17-3.51)2=0.05445;2 2 2SS(AXB)=SS3=(1/8)R(AxB)2=(1/n)(K1-K2)2=(1/8)(3.66-4.02)2=0.0162;SSC=SS4=(1/8)RC2=(1/n)(K1-K2)2=(1/8) (3.93-3.75)2=0.00405;SS(AXC)=SS5=(1/8)R(AxC)2=(1/n)(K1-K2)2=(1/8)(3.5-4.18)2=0.0578;SS(BXC)=SS3=(1/8)R(

10、BxC)2=(1/n)(K1-K2)2=(1/8)(3.66-4.02)2=0.0162;误差的离差平方和：SSe=SS7=(1/8)R72=(1/n)(K1-K2)2=(1/8)(3.63-4.05)2=0.02205或:SS=SST-(SSA+SSB+SS(AXB)+SSC+SS：AXC)+SS(BXc)=0.4088-(0.23805+0.05445+0.0162)=0.02205计算自由度：总自由度:dfT= n-1=8-1=7;各因素自由：dfA=dfB=dfc=r-1=2-1=1; df(AxB)=dfAX dfB=1X1=1或df(AXB)=df3=r-1=2-1=1;同理:df

11、(AXc)=dfAX dfc=df5=1 X 1=1;df(Bxc)=dfBX dfc=df6=1 X 1=1;误差自由度：dfe=df7=r-1=2-1=1；或dfe=dfT-(dfA+dfB+df(AxB)+dfc+df(Axc)+df(Bxc)=7(1+1+1+1+1+1)=1(4)计算均方： 由于各因素、交互作用和误差的自由度都为1,所以它们的均方应该等于它们各自的离差平方和，即：MSA=SSA=0.23805;MSB=SSB=0.05445;MS(AXB)=SS(AXB)=0.0162;MSC=SSC=0.00405;MS(AXC)=SSAXC)=0.0578;MS(BXC)=SSB

12、Xc)=0.0162;MSe=SSe=0.02205;k11.1331.0430.9150.9830.8750.9150.908k20.7880.8781.0050.9381.0451.0051.013极差R1.3800.6600.3600.1800.6800.3600.420因素 主次A,B X C,B优方案A1B1OP=7.3728可知：MS(AxB)vMSe、MScvMSe、MS(BXc)Mse,这说明因素C和交互因素A X B、交互因素B X C对试验结果的影响较小，为次要因素。所以，可将它们归入误差，这样误差的离差平方和、自由度、和和均方都 会随之变化，即新误差平方和：SSe=SS3

13、+SS(AXB)+SSC+SS(Bc)=0.02205+0.0162+0.00405+0.0162=0.0585;误差自由度：B)+dfc+df(BXC)=1+1+1+1=4;新误差均方：MSe二SSe/dfe=0.0585/4=0.014625(5)计算F值：FA=MSA/MSe=0.23805/0.014625=16.2764;岳=MSB/MSed=0.05445/0.014625=3.7230;FAXC)=MS(AXC)/MSe=0.0578/0.014625=3.9531;由于：A X B、C、B X C已并入误差，所以就不需要计算它们的F值。(6)F值检验：查得临界值F0.05(1,

14、4)=7.71,F0.01(1,4)=21.20,所以对于给定 显著水平a =0.05，因素A对试验结果有显著影响，而B,A X C对试验结果影响 不显著(对于主要因素，一定要按有利于指标的要求选取最好的水平；而对于不重要的因素，由于其水平改变对试验结果的影响较小，则可以根据有利于降低消耗、提高效率等目的来考虑别的水平。)最后将结果列于表1-3例1方差分析表 方差分析表中(见表1-3)。差异源SSdfMSF显著性A0.23810.23816.276*B0.05410.0543.726A X C0.05810.0583.953A X B0.0161C0.00410.015B X C0.0161误

15、差 e0.0221总和0.4097从表中F值的大小可看出因素主次顺序为：A、A X C、B,这与极差分析结果 是一致的。(7)优方案的确定：因为因素A对试验结果有显著影响，因素B、C和交互 作用AXC,A X B对试验结果无显著影响，所以从有利利于降低消耗，提高效 益等目的来考虑别的水平，所以，因素B取Bi(用水量20ml)、因素C取C1(反 应时间1h)。因素A取K值对应的最大水平，即Ai酸用量25ml。所以最优方案 为：AiBiCi。2、为了提高陶粒混凝土的抗压强度，考察了A,B,C,D,E,F六个因素，每个 因素有六个水平，因素水平表如下(见表2-1):表2-1因素水平表水平(A)水泥标

16、号(B)水泥用量/kg(C)陶粒用量/kg(D)含砂率/%(E)养护方式(F)搅拌时间/h130018015038空气1240019018040水1.5Xdfe =dfe+df(AX213500200200421蒸气2根据经验还要考察交互作用A X B,A X C,BX C。如果将A,B,C,D,E,F依次安排在 正交表L27（313）的1,2,5,9,12,13列上，试验结果（抗压强度/kg.cm-2）依次为100,98,97,95,96,99,94,99,101,85,82,98,85,90,85,91,89,80,73,90,77,84,80,76,89,78,85试用方差分析法（a =

17、0.05） 分析试验结果，确定较有水平组合。解：（1）实验设计：本实验要考虑六个因素和三种交互作用，且每种交互作用占两列，这样因素和交互作用在正交表中总共占有12列，所以应该选择正交表L27（313）。根据L27（313）的交互作用表和题目要求进行表头设计，然后进 行试验，得到试验结果yi（1,2,，27）。试验设计及结果列于表2-2。表2-2试验设计及结果12345678910111213抗压强度(A(A(A(A(B空 列(Byi/(kg.cm-2)ABXXCXXXDXEFB)1B)2C)1C)2C)1C)2111111111111110011112222222229811113333333

18、339712221112223339512222223331119612223331112229913331113332229413332221113339913333332221111012123123123123852123231231231822123312312312982231123231312852231231312123902231312123231852312123312231912312231123312892312312231123803132132132132733132213213213903132321321321771234567891011121314151617

19、181920试验号(1)求Ki:在公式，然后水平拖动填充柄，就可计算出后12列的Ki值。同理，求K2：在B38中输入：=SUMIF(B$10:B$36,2,$0$10:$0$36),然后水平拖动填充柄，就可计算出 后12列的K2值。求K3:在B39中输入：=SUMIF(B$10:B$36,3,$0$10:$0$36)， 然后水平拖动填充柄，就可计算出后12列的K2值。2求k1：首先选中单元分区域B40：N42，在该区域的左上角第一个单元格输 入：=B37：N39/3，再同时按住“Shift+Ctrl+Enter,即可在B40：N42显示出结果。3求极值R：在B43单元格中输入：=MAX(B37

20、:B39)-MIN(B37:B39)，回车， 然后选中该单元格，选中该填充柄，然后水平拖动该填充柄，就可计算出后12列的极差R值。(2)计算离差平方和：T=SUM(010:036)=(y1+y2+y27)=2396,Q=SUM(y12+y22+y27)=214418, P=T2/n=23962/27=212622.82所以离差平方和：SST=Q-P=214418-212622.82=1796.739SSB=1/9(K12+K22+K32)=1/9(8002+7902+8062)-212622.82=14.518，SS(AXB)1=1/9(K12+K22+K32)=1/9(7952+8072+7

21、842)-212622.82=1 1 .629，SS(AXB)2=1/9(K12+K22+K32)=1/9(8072+7902+7992)-212622.82=16,074，SSC=1/9(K12+K22+K32)=1/9(7962+8022+7982)-212622.82=2.074，SS(AXC)1=1/9(K12+K22+K32)=1/9(8002+7922+8042)-212622.82=8.296，SS(AXC)2=1/9(K12+K22+K223213132213321842332132133211328024321332113221376253321132321213892633

22、212131323217827332132121313285K1879800795807796800788784836791790824794K2785790807790802792802800799800804785801K3732806794799798804806812761805802787801k1293.0266.7265.0269.0265.3266.7262.7261.3278.7263.7263.3274.7264.7k2261.7263.3269.0263.3267.3264.0267.3266.7266.3266.7268.0261.7267.0k3244.0268.72

23、64.7266.3266.0268.0268.7270.7253.7268.3267.3262.3267.0T=2396，极值1471613176121828751414397P=212622.82 ，Q=214418因素主次优方案A,D,E,B X C,B,A X C,AX B,F,CA1B1C1D1E1F1中输入：=SUMIF(B$10:B$36,1,$0$10:$0$36),选中该B3732)=1/9(7882+8022+8062)-21 2622.82=1 9.852,SS(BXC)1=1/9(K12+K22+K32)=1/9(7842+8002+8122)-212622.82=43.

24、852，SSD=1/9(K12+K22+K32)=1/9(8362+7992+7612)-212622.82=312.518,SSe=1/9(K12+K22+K32)=1/9(7912+8002+8122)-212622.82=11.185,SS(BxC)2=1/9(K12+K22+K32)=1/9(7902+8042+8022)-212622.82=12.741,SSE=1/9(K12+K22+K32)=1/9(8242+7852+7872)-212622.82=107.182,SSF=1/9(K12+K22+K32)=1/9(7942+8012+8012)-212622.82=3.629所

25、以：SS(AXB)=SS(AxB)1+SS(AxB)2=11.629+16.074=27.703, SSAxC)=SS(AxC)1+SS：AxC)2=8.296+19.852=28.148, SSBxC)=SSBxC)1+SSBXC)2=43.852+12.741=56.593(3)计算自由度：总自由度：dfT= n-1=27-1=26,各因素自由度：dfA=dfB=dfc=dfD=dfE=dfF=r-1=3-1=2,交互作用自由度：df(AxB)=dfAx dfB=2x 2=4或df(AxB)=df(AxB)1+df(AXB)2=2+2=4,同理：df(Axc)=dfAx dfc=2 x 2

26、=4,df(Bxc)=dfBx dfc=2 x 2=4,误差自由度：dfe=r-1=3-1=2或dfe=dfT-dfA-dfB-dfc-dfD-dfE-dfF-df(AxB)-df(AxC)-df(BxC)=26-2-2-2-2-2-2-4-4-4=2;(4) 计算均方：各因素和交互作用的均方为：MSA=SSA/dfA=1231.629/2=615.815,MS3=SSB/dfB=14.518/2=7.259,MS(AXB)=SS：AxB)/df(AxB)=27.703/4=6.92,MS(AXC)=SSA AXc)/df(AxC)=28.148/4=7.037,MSc=SSc/dfc=2.0

27、74/2=1.037, MSD=SSD/dfD=312.518/2=156.259,MSE=SSE/dfE=107.182/2=53.591,MSF=SSF/dfF=3.629/2=1.81,MS(Bxc)=SSBxC)/df(BxC)=56.593/4=14.148,但误差的均方为：MSe=SSe/dfe=11.185/2=5.593。计算到这里,我们发现MScMSe,MSFMSe,这说明因素C、F对试验结果的影响较小,可以将它们归 入误差，这样新误差离差平方和为：SSe=SSc+SSF+SSe=3.629+11.185+2.074=16.888,新误差自由度为：dfe=dfF+dfc+dfe=2+2+2=6,新误差均方为：MSe=SSef=16.888/6=2.815(5) 计算F值：FA=MSA/MSe=615.815/2.815=218.762,FB=MSB/MSe=7.259/2.815=2.579,xB)/MSe=6.926/2.815=2.460,xC)/M