现代试验方法与数据处理作业 3.docVIP

PAGE PAGE 6 第三次作业 4.用某种菌生产酯类风味物质，为了寻找最优的发酵工艺条件，重点考察了葡萄糖用量的x1（50~150g/L)和蛋白胨用量x2(2~10g/L)影响，试验指标为菌体生长量y(g/L),其他的发酵条件不变。试验方案结果如下： 试验号 Z1 Z2 y 1 1 1 9.61 2 1 -1 9.13 3 -1 1 9.37 4 -1 -1 8.57 5 1.078 0 9.34 6 -1.078 0 8.97 7 0 1.078 10.21 8 0 -1.078 9.48 9 0 0 10.24 10 0 0 10.33 试用二次回归正交设计在试验范围内建立二次回归方程； 对回归方程和回归系数进行显著性检验； 失拟性检验； 试验范围内最优试验方案的确定。 解：(1)确定因素的变化范围： 由于因素数为2，取零水平实验次数，根据星号臂长r的计算公式或查表得r=1.078 因为葡萄糖用量（X1r）上限为150 g/L ，下限（X-1r）为50 g/L 零水平X10=(X1r+X-1r)/2=(150+50)/2=100 g/L , 变化间距 =(X1r-X10)/r=（150-100）/1.078=46.382， 上水平X12=X10+ =100+46.382=146.382 g/L , 下水平X11= X10- =100-46.382=53.618 g/L。 X12对应的编码为Z12=(X12-X10)/ =(146.382-100)/46.382=1, 同理可以求出其它因素的水平及编码如下： 因素水平编码表 规范变量Zj 自然变量Xj X1/(g/L) X2(g/L) 上星号臂r 150 10 上水平1 146.382 9.711 零水平0 100 6 下水平-1 53.618 2.289 下星号臂-r 50 2 变化间距 46.382 3.711 正交组合设计 本实验因素数m=2，参照正交表可选取 进行变换，二水平实验次数，星号实验为2m=4，并取零水平实验数，实验方案见下表： 实验号 Z1 Z2 葡萄糖用量 X1（g/L） 蛋白胨用量 X2（g/L） 1 1 1 146.382 9.711 2 1 -1 146.382 2.289 3 -1 1 53.618 9.711 4 -1 -1 53.618 2.289 5 1.078 0 150 6 6 -1.078 0 50 6 7 0 1.078 100 10 8 0 -1.078 100 2 9 0 0 100 6 10 0 0 100 6 二元二次回归正交组合设计表及实验结果： 实验号 Z1 Z2 Z1Z2 Z1*Z1 Z2*Z2 Z1’ Z2’ Y 1 1 1 1 1 1 0.368 0.368 9.61 2 1 -1 -1 1 1 0.368 0.368 9.13 3 -1 1 -1 1 1 0.368 0.368 9.37 4 -1 -1 1 1 1 0.368 0.368 8.57 5 1.078 0 0 1.162 0 0.530 -0.632 9.34 6 -1.078 0 0 1.162 0 0.530 -0.632 8.97 7 0 1.078 0 0 1.162 -0.632 0.530 10.21 8 0 -1.078 0 0 1.162 -0.632 0.530 9.48 9 0 0 0 0 0 -0.632 -0.632 10.24 10 0 0 0 0 0 -0.632 -0.632 10.33 ∑ 95.25 回归方程的建立 所以规范变量与实验指标y之间的回归关系为： y=9.525+0.1896Z1+0.3268Z2-0.08Z1Z2-0.83Z1’-0.236Z2’ （2）回归方程及偏回归系数的显著性检验 方差分析结果如下表： 差异源 SS df MS F 显著性 Z1 0.277 1 0.227 6.91 Z2 0.675 1 0.657 20.55 * Z1Z2 0.0256 1 0.0256 0.78 Z1’ 1.861 1 1.861 56.65 ** Z2’ 0.15 1 0.15 4.57 回归 2.94 5 0.588 17.89 ** 残差e 0.1314 4 0.03285 总和 3.07 n-1=9 由表可知因素Z1、Z1Z2、Z1’ 对实验指标影响不显著，因此可以将其平方和及自由度归入残差，再进行一次方差分析，分析结果如下表： 二次方差分析表 差异源 SS df MS F 显著性 Z1 Z1’ 0.675 1.861 1 1 0.657 1.861 9.25 25.49 * ** 回归 残差e 2.536 0.534