正交设计应用实例(毕业论文)

2正交实验设计2.1正交实验设计概述正交实验设计(Orthogonalexperimentaldesign)11是研究多因素多水平的一种设计方法，它是根据从全面实验中挑选出部分有代表性的点进行实验，正交实验设计又称正交设计或多因素优选设计，是一种合理安排、科学分析各实验因素的一种有效的数理统计方法。它是在实践经验和理论认识的基础上，借助一种规格化的“正交表”，从众多的实验条件中确定出若干个代表性较强的实验条件，科学地安排实验，然后对实验结果进行综合比较，统计分析，探求各因素水平的最佳组合，从而得到最优或较优实验方案的一种实验设计方法。正交实验设计的特点是用不太多的实验次数，找出实验因素的最佳水平组合，了解实验因素的重要性程度及交互作用情况，减少实验盲目性，避免资金浪费等。它能以较少的实验次数找到较好的实验（生产）方案，由正交实验寻找出的优化参数（条件）与全面实验所找出的最优条件有一致的趋势。正交实验设计具有正交性，使实验具备均衡分散和综合可比性。此法应用方便，准确性高，在多因素条件下应用有很大的优越性,是一种高效率、快速、经济的实验设计方法。日本著名的统计学家田口玄一将正交实验选择的水平组合列成表格，称为正交表。例如作一个三因素三水平的实验，按全面实验要求，须进行33=27种组合的实验，且尚未考虑每一组合的重复数。若按L9(3)3正交表按排实验，只需作9次，显然大大减少了工作量。因而正交实验设计在很多领域的研究中已经得到广泛应用。2．2正交实验设计基本程序正交设计中常用的术语有：指标、因子和水平。正交设计把实验设计要考虑的结果和评价准则称为指标，一般以yi表示第i次实验的指标值；把对实验结果和对评价指标可能产生影响且在实验中明确了条件加以对比的因素称为因子，一般以大写字母表示；把每个因子在实验中的具体条件称为因子的水平，简称水平，一般以表示因子的大写字母加上脚标来表示。对于多因素实验，正交设计是简单常用的一种设计方法，其设计程序12如图4所示。图4正交实验设计程序结论优组合因数主次顺序优水平计算极差R值计算K值绘制因数指标趋势图实验结果极差分析进行试验，记录实验结果列实验方案表头设计选择合适的正交表实验指标选因数、定水平因数筛选因数、水平确定因果分析实验目的与要求图4正交实验设计程序结论优组合因数主次顺序优水平计算极差R值计算K值绘制因数指标趋势图实验结果极差分析进行试验，记录实验结果列实验方案表头设计选择合适的正交表实验指标选因数、定水平因数筛选因数、水平确定因果分析实验目的与要求Figure4Theprogrammingofpotassiumdiformate2.3正交实验方案设计2.3.1明确实验目的，确定实验指标实验设计前必须明确实验目的，即本次实验要解决什么问题。实验目的确定后，对实验结果如何衡量，即需要确定出实验指标。实验指标可为定量指标，如强度、硬度、产量、出品率、成本等；也可为定性指标如颜色、口感、光泽等。一般为了便于实验结果的分析，定性指标可按相关的标准打分或模糊数学处理进行数量化，将定性指标定量化。2.3.2选择实验因素，确定实验水平，列出因素水平表根据专业知识、以往的研究结论和经验，从影响实验指标的诸多因素中，通过因果分析筛选出需要考察的实验因素。一般确定实验因素时，应以对实验指标影响大的因素、尚未考察过的因素、尚未完全掌握其规律的因素为先。实验因素选定后，根据所掌握的信息资料和相关知识，确定每个因素的水平，一般以2-4个水平为宜。对主要考察的实验因素，可以多取水平，但不宜过多（≤6），否则实验次数骤增。因素的水平间距，应根据专业知识和已有的资料，尽可能把水平值取在理想区域。2.3.3选择合适的正交表,进行表头设计正交表的选择是实验设计的首要问题。正交表选得太小，实验因素可能安排不下；正交表选得过大，实验次数增多，不经济。正交表的选择原则是在能够安排下实验因素和交互作用的前提下，尽可能选用较小的正交表，以减少实验次数。表头设计就是指将实验因素和交互作用合理地安排到所选正交表的各列中去的过程。若实验因素间无交互作用，各因素可以任意安排；若要考察因素间的交互作用，各因素应按相对应的正交表的交互作用列表来进行安排，以防止设计“混杂”。正交表是一整套规则的设计表格，表示方法记为Lm(rn)，其中L为正交表的代号，m为实验的次数，r为水平数，n为列数，也就是可能安排最多的因素个数。例如L9(33)，(表1)，它表示需作9次实验，最多可观察3个因素，每个因素均为3水平。一个正交表中也可以各列的水平数不相等，我们称它为混合型正交表。根据正交表的数据结构看出，正交表是一个m行n列的表，表11中，大写字母表示A、B、C三个因素，下标表示每个因素有1、2、3三个水平。表1三因数三水平正交实验表Table1Theorthogonaltesttableofthreefactorsandthreelevels因因素实验号ABC效益平均值Y1A1B1C12A1B2C23A1B3C34A2B1C25A2B2C36A2B3C17A3B1C38A3B2C19A3B3C2正交表具有以下两项性质：(1)每一列中，不同的数字出现的次数相等。例如在三水平正交表中，任何一列都有下标“1”、“2”、“3”，且在任一列的出现次数均相等。

(2)任意两列中数字的排列方式齐全而且均衡。例如在三水平情况下，任何两列(同一横行内)有序对共有9种，1.1、1.2、1.3、2.1、2.2、2.3、3.1、3.2、3.3，且每对出现数也均相等。

以上两点充分的体现了正交表的两大优越性，即“均匀分散性，整齐可比”。通俗的说，每个因素的每个水平与另一个因素各水平各碰一次，这就是正交性。2.3.4编制实验方案，按方案进行实验，记录实验结果根据实验方案进行实验，记录实验结果。2.4正交实验结果分析极差分析法（R法）又称直观分析法，此法计算简便而直观，简单、易懂，是正交实验结果分析最常用的方法。2.4.1确定实验因素的优水平和最优水平组合最佳水平是指每个因子的各水平中使指标达最佳的水平。为确定因子的最佳水平，必须确定该因子各水平对指标的影响。分析A因素各水平对实验指标的影响。A因素的i水平所对应的实验指标平均值用kAi表示（Kij—第j个因素第i个水平的所有实验结果指标值的均值）。根据正交设计的特性，对A1、A2、A3来说，三组实验的实验条件是完全一样的（综合可比性），可进行直接比较。如果因素A对实验指标无影响时，那么kA1、kA2、kA3应该相等，但kA1、kA2、kA3实际上不相等。说明，A因素的水平变动对实验结果有影响。因此，根据kA1、kA2、kA3的大小可以判断A1、A2、A3对实验指标的影响大小，k值大的对实验指标的影响大，为A因素的优水平。同理，可以计算并确定B、C因素的优水平。三个因素的优水平组合为实验的最优水平组合。2.4.2确定因素的主次顺序各因子对指标的影响是不同的，其重要性也各不相同。为了评价各因子的重要性，需拟定一评价指标。通常采用均方和S或极差R作为评价指标。Sj—均方和Sj=[(K1j-K2j)2+(K1j-K3j)2+(K2j-K3j)2]MAX=MAX(K1j,，K2j，K3j)MIN=MIN(K1j,，K2j，K3j)极差=MAX－MIN根据极差Rj的大小，可以判断各因素对实验指标的影响主次。极差Rj各数据间的差距越大，说明该因子各水平相差悬殊，对指标的影响大，反之则小。因此以均方和或极差可粗略的揭示出各因子的重要性。2.4.3因子显著性的检验因子的重要性只说明该因子相对其他因子的重要程度，而未说明该因子对指标影响的显著程度。如果某因子对指标的作用不显著，则可排除该因子而使决策简化。经显著性检验之后，可确定对指标有显著影响的因子、排除对指标影响不显著的因子。在此基础之上可选择与确定最佳方案。当计算出的F值大于临界值,k因子在α水平下作用显著，否则作用不显著。最佳方案的确定方法是选择对指标有显著影响的因子中的最佳水平，对于对指标无显著影响的因子可不考虑，或根据实际情况决定。2.4.4绘制因素与指标趋势图以各因素水平为横坐标，实验指标的平均值（kjm）为纵坐标，绘制因素与指标趋势图。由因素与指标趋势图可以更直观地看出实验指标随着因素水平的变化而变化的趋势，可为进一步实验指明方向。2．5正交实验总结1在实际生产中，最优条件的确定是灵活的。对于主要的影响因素，一定要选最优水平，而对于次要因素，则应权衡利弊，综合考虑来选取优水平。2极差分析得到的最优工艺条件并不一定在所实施的正交实验方案中。为了考察最优条件的再现性，应追加验证性实验，从而进一步判断研究所找出的生产工艺条件是否最优。3由极差分析得出的最优工艺条件，只有在实验所考察的范围内有意义。4正交实验的极差分析简明、直观，是有效的数理统计方法，应提倡大力推广。3.1材料与方法3.1.1原料与仪器甲酸（化学纯，85%）、无水碳酸钾（化学纯，99%）、蒸馏水、邻苯二甲酸氢钾(基准试剂，99.8%)、氢氧化钠（化学纯，96%）、酚酞；三颈烧瓶，冷凝管，恒温水槽，磁力搅拌器，抽滤漏斗，真空干燥箱，电子天平。3.1.2实验原理甲酸与碳酸钾按4：1的物质的量比在一定温度、时间条件下反应后，通过降温结晶，分离，干燥，得到的白色晶体二甲酸钾。其化学反应4HCOOH+K2CO3=2KH(HCOO)2+CO2+H3O二甲酸钾的生产工艺流程图如下无水碳酸钾合成降温结晶无水碳酸钾合成降温结晶分离干燥甲酸产品二甲酸钾纯度的测定有反相高效液相色谱法和酸碱滴定法。反相高效液相色谱法是以二甲酸钾标准品的HPLC图谱和二甲酸钾产品HPLC图谱比较，运用峰高外标定量法即可计算出产品纯度；酸碱滴定法是用已知浓度的氢氧化钠溶液滴定二甲酸钾，其反应如下:HCOOH-HCOOK+NaOH=NaCOOH-HCOOK+H2O二甲酸钾与氢氧化钠是一比一反应，称取一定量二甲酸钾加水溶解，用已知浓度的氢氧化钠溶液滴定，滴入两滴酚酞做指示剂，到滴定终点记下所用氢氧化钠的体积，算出所用氢氧化钠的量就能算出二甲酸钾的量，从而可算出二甲酸钾的纯度。综合考虑到实验条件和实验准确性，本实验用滴定法测定二甲酸钾纯度。3.2实验步骤3.2.1配制标准氢氧化钠溶液用电子天平称取二甲酸钾产品，从准确性考虑：量不能太少。若定为取产品1.5g，二甲酸钾分子量为130.12，则滴定所用氢氧化钠的量为1.5g÷130.12g/mol=0.011528mol,用50ml碱式滴定管滴定则需用25.00ml，所配氢氧化钠溶液的浓度应该为0.01537mol÷0.025l=0.4611mol/l,九次实验，每次产品滴定3次至少要滴定27次，一次滴定要25.00ml，需用1000ml的容量瓶配制。氢氧化钠纯度为96%，总共需称取粗氢氧化钠：0.4611mol/L×1L×40g/mol÷0.96=19.21g.3.2.1.1配制浓度为0.5mol/L的氢氧化钠溶液把烧杯放在天平上，回零后称取氢氧化钠约20g,用蒸馏水溶解后，冷却到室温，用玻璃棒引流转移到1000ml容量瓶，玻璃棒和烧杯用蒸馏水洗三遍移入容量瓶，然后加蒸馏水至刻度线，塞紧瓶塞后摇匀即可。3.2.1.2标定氢氧化钠溶液取在105℃干燥至恒重的基准邻苯二甲酸氢钾约2g,精密称定质量为m=2.013克,加蒸馏水振摇,使其尽量溶解,加酚酞指示液2滴,用本液滴定,在接近终点时,溶液由无色显粉红色。记录氢氧化钠溶液的体积v=22.30ml，计算氢氧化钠的准确浓度为0.44113mol/l,公式为CNaOH=m×99.8%÷M×1000÷做三次平行实验，把实验数据及计算结果填入表2。表2氢氧化钠标定Table2Thecalibrationofsodiumhydroxide邻苯二甲酸氢钾质量m(g)氢氧化钠溶液体积V(ml)氢氧化钠溶液浓度C(mol/l)2.01322.300.441132.03722.750.437562.02622.700.43616计算氢氧化钠溶液平均浓度：CNaOH=(0.44113mol/l+0.43756mol/l+0.43616mol/l)÷3=0.4383mol/l3.2.2合成二甲酸钾本实验考察甲酸与碳酸钾合成二甲酸钾，考察指标为二甲酸钾的收率。查二甲酸钾合成专利及各类文献，可知甲酸与碳酸钾合成二甲酸钾影响二甲酸钾收率的主要因数为反应摩尔比、反应温度、反应时间，其最优条件范围为反应摩尔比（甲酸：碳酸钾=4.0～5.0），反应温度（50℃～60℃），反应时间（1.5h～3h）表3实验因数与水平Table3experimentalfactorsandlevels因因子水平摩尔比（甲酸：碳酸钾）反应温度(℃)反应时间(h)14.0:1.0501.524.5:1.0552.035.0:1.0602.5表4.摩尔比例与质量对应表摩尔比甲酸用量(85%)(g)碳酸钾用量(99%)（g）4.0:1.064.941.8824.5:1.073.141.8825.0:1.081.241.882具体实验正交表安排如下表：表5实验正交表排列Table5Thetableoforthogonaltest列列号水平实验号A摩尔比B反应时间(h)C反应温度(℃)收率Y14.0:1.01.55024.0:1.02.05534.0:1.02.56044.5:1.01.55554.5:1.02.06064.5:1.02.55075.0:1.01.56085.0:1.02.05095.0:1.02.555为了更好更快完成实验，排除外界环境因数对实验产生误差，先把以上九次实验按反应温度分成三组：表6.实验分组表因子分组A摩尔比B反应时间（h）C反应温度(℃)第一组4.0:1.01.5504.5:1.02.55.0:1.02.0第二组4.0:1.02.0554.5:1.01.55.0:1.02.5第三组4.0:1.02.5604.5:1.02.05.0:1.01.53.2.2.1实验步骤（一）接通电源，打开恒温水浴，把温度设定到50℃（二）分别取三个大烧杯三个小烧杯编号并称重，用小烧杯称取41.882克99%的碳酸钾试剂三份，用大烧杯称取64.9克、73.1克、81.2克85%的甲酸试剂各一份。（用天平称取干燥时的烧杯重后，加入所需称量的试剂至所需量。）（三）向甲酸里缓慢加入碳酸钾，为减少甲酸挥发用玻片遮着点;反应过程中不断搅拌，一是可以使反应充分，二是可以使温度不太高。（四）反应后，把三个烧杯一起放入恒温水浴，按表6所列反应时间依次取出各烧杯，在同一条件下冷却到室温结晶。（五）安装好吸滤瓶与布氏漏斗，剪两张比漏斗内颈稍小的圆形滤纸。室温为16℃3.2.2.2用氢氧化钠滴定产品干燥大约24小时，称量产品重量直至重量不再减少，可认为产品已经干燥，记录产品重量；由反应式计算出二甲酸钾理论值，然后计算出二甲酸钾的收率,见表7。收率=实际产品值÷理论产品值×100%用电子天平准确称取制得的二甲酸钾晶体1.2-1.3g，溶解在80ml蒸馏水中，然后用已标定浓度的氢氧化钠溶液滴定，酚酞做指示剂，计算出二甲酸钾的纯度,见表8。表7二甲酸钾收率计算结果Table7Thereceivingrateofpotassiumdiformate’spurities实验号原烧杯重(g)烘干后总重(g)实际产品值(g)理论产品值(g)收率1165.201226.01860.81778.7300.772289.333144.90155.56878.7220.7063151.875203.67551.80078.7300.6584190.037237.97147.93478.7300.609573.881120.28946.40878.7040.5906136.303182.38946.08678.7800.5857179.215213.61334.39878.7070.437881.389114.93133.54278.7370.4269194.110226.48032.37078.7600.411表8二甲酸钾纯度的计算结果Table8Thecalculationofpotassiumdiformate’spurities实验号二甲酸钾质量（g）氢氧化钠体积（ml)二甲酸钾纯度收率平均纯度11.32222.460.9710.7730.9681.33025.060.9651.32627.600.96821.31825.180.9730.7060.9721.32322.760.9701.31925.320.97331.27721.900.9780.6580.9761.26521.660.9771.29822.140.97341.21020.600.9710.6090.9711.20620.480.9681.21320.700.97351.19220.480.9800.590.9811.22621.100.9821.27822.000.98261.17920.280.9810.5850.9831.16920.160.9841.23221.280.98571.21820.940.9800.4370.9821.23621.240.9801.22221.100.98581.32427.880.9690.4260.9711.31825.460.9731.32130.040.97191.28622.080.9790.4110.9791.29522.220.9791.31322.560.9804实验结果分析4.1对正交实验表分析把实验结果填入正交实验表，用Excel计算各数据值。表9正交实验分析结果Table9Theanalysisoforthogonaltest因因数水平实验号AHCOOH:K2CO3B反应时间(h)C反应温度(℃)D收率(%)14.0:1.01.55077.324.0:1.02.05570.634.0:1.02.56065.844.5:1.01.55560.954.5:1.02.06059.064.5:1.02.55058.575.0:1.01.56043.785.0:1.02.05042.695.0:1.02.55541.1K171.2360.6359.47K259.4757.4057.53K342.4755.1356.17S1254.9845.8416.50MAX71.23860.6359.47MIN42.4755.1356.17R28.775.503.30优水平A1B1C1表中K1,K2.K3，表示每个因素所进行的3次实验(3水平)所得二甲酸钾收率的平均值，即在代表因素A(摩尔比)的第一列中，将与水平1对应的l,2,3共3个二甲酸钾收率相加除3所得，其他类推。Ri为平均数的极差(K中最大值与最小值之差)。Ri愈大表明变化幅度大，说明该因素的水平变化对二甲酸钾收率的影响愈大，该因素即是最大因素。注：针对每个因素，其哪个水平的K值大，该水平即为该因素的最优水平。其中，针对每个因素Kij—第j个因素第i个水平的所有试验结果指标值的均值；Sj—均方和Sj=[(K1j-K2j)2+(K1j-K3j)2+(K2j-K3j)2]