整理响应面优化实验方案设计.docx
《整理响应面优化实验方案设计.docx》由会员分享，可在线阅读，更多相关《整理响应面优化实验方案设计.docx（28页珍藏版）》请在冰豆网上搜索。
整理响应面优化实验方案设计
食品科学研究中实验设计的案例分析
——响应面法优化超声辅助提取车前草中的熊果酸
班级：
学号：
姓名：
摘要：
本文简要介绍了响应面曲线优化法的基本原理和使用步骤，并通过软件Design-Expert7.0软件演示原文中响应面曲线优化法的操作步骤。
验证原文《响应面法优化超声辅助提取车前草中的熊果酸》各个数据的处理过程，通过数据对比，检验原文数据处理的正确与否。
关键词：
响应面优化法数据处理Design-Expert7.0车前草
前言：
响应曲面设计方法（ResponseSufaceMethodology，RSM）是利用合理的试验设计方法并通过实验得到一定数据，采用多元二次回归方程来拟合因素与响应值之间的函数关系，通过对回归方程的分析来寻求最优工艺参数，解决多变量问题的一种统计方法（又称回归设计）。
响应面曲线法的使用条件有：
①确信或怀疑因素对指标存在非线性影响；②因素个数2-7个，一般不超过4个；③所有因素均为计量值数据；试验区域已接近最优区域；④基于2水平的全因子正交试验。
进行响应面分析的步骤为：
①确定因素及水平，注意水平数为2，因素数一般不超过4个，因素均为计量值数据；②创建“中心复合”或“Box-Behnken”设计；③确定试验运行顺序（DisplayDesign）；④进行试验并收集数据；⑤分析试验数据；⑥优化因素的设置水平。
响应面优化法的优点：
①考虑了试验随机误差②响应面法将复杂的未知的函数关系在小区域内用简单的一次或二次多项式模型来拟合，计算比较简便，是降低开发成本、优化加工条件、提高产品质量，解决生产过程中的实际问题的一种有效方法③与正交试验相比，其优势是在试验条件寻优过程中，可以连续的对试验的各个水平进行分析，而正交试验只能对一个个孤立的试验点进行分析。
响应面优化法的局限性:
在使用响应面优化法之前，应当确立合理的实验的各因素和水平。
因为响应面优化法的前提是设计的试验点应包括最佳的实验条件，如果试验点的选取不当，实验响应面优化法就不能得到很好的优化结果。
原文《响应面法优化超声辅助提取车前草中的熊果酸》采用经典的三因素三水平Box-Behnken试验设计，以熊果酸的提取率为响应值，通过回归分析各工艺参数与响应值之间的关系，并由此预测最佳的工艺条件。
本文利用软件验证原文中的数据处理过程，以检验原文数据是否处理正确。
1确定实验因素
原文利用超声波辅助提取车前草中的熊果酸，而影响熊果酸提取率的因素有很多，如超声波的功率、提取时间、溶剂温度、溶剂种类、液固比等。
原文参考文献《柿叶中总三萜的提取以及熊果酸分离,纯化研究》中提取熊果酸的方法提取熊果酸，即将干燥的车前草粉碎后过筛，取20～40目的车前粉，用石油醚在55℃脱脂3次，干燥备用。
精密称取一定量的车前粉，加入一定量的乙醇，称量，在一定的超声波功率下提取一定时间后，擦干外壁，再称量，用乙醇补充缺失的质量，离心。
用注射器抽取一定量上清液，过0.45μm滤膜，进行检测。
每个实验进行3次平行实验。
取其平均值。
结果以提取率（E）的来表示。
C×V
E/%=———×100
M
式中：
C为熊果酸的质量浓度/（g/mL）；V为加入乙醇的体积/mL；m为车前草的质量/g。
在一系列单因素实验的基础上，采用经典的三因素三水平Box-Behnken试验设计，选取提取温度（A）、乙醇体积分数（B）、超声功率（C）三个因素作为优化条件的因素对象。
2确定因素水平范围
确定因素水平范围就是通过做单因素初步试验或由样品的特性和工艺来确定BBD设计所研究的因素水平范围。
确定合适的因素水平范围对获得理想的优化结果非常重要，如果水平范围太窄得不到优化结果，太宽也会使结果精确度降低。
原文在单因素实验的基础上确定了因素水平范围是：
提取温度的：
60—80℃；乙醇体积分数：
90—100%；提取功率：
420—540W
3试验设计安排与结果
根据Box-Behnken中心组合设计原理,在单因素试验的基础上，以提取温度、乙醇体积分数和提取功率三个因素为自变量，熊果酸提取率%为响应值，作三因素三水平的响应面分析试验，共17个试验点。
其中12个为析因子，5个为中心试验用以估计误差。
试验因素和水平见表一。
表一响应面试验因素水平表
Table1Factorsandlevelsinresponsesuefacedesign
水平
提取温度（℃）
乙醇体积分数（%）
提取功率（W）
-1
60
80
420
0
70
90
480
1
80
100
540
表二响应面试验方案及结果
Tabel2Schemeandexperimentalresultsofresponsesurfacedesign
试验号
A
B
C
提取率（%）
1
-1
-1
0
0.054403
2
1
-1
0
0.100936
3
-1
1
0
0.070440
4
1
1
0
0.094782
5
-1
0
-1
0.062646
6
1
0
-1
0.095539
7
-1
0
1
0.063370
8
1
0
1
0.100653
9
0
-1
-1
0.085573
10
0
1
-1
0.078006
11
0
-1
1
0.098521
12
0
1
1
0.086069
13
0
0
0
0.108560
14
0
0
0
0.101516
15
0
0
0
0.101054
16
0
0
0
0.102413
17
0
0
0
0.105489
4用软件（Design-Expert）对实验数据统计分析
4.1打开DesignExpert7.0软件
4.2数据输入
4.2.1因素输入
4.2.2响应值输入
4.3试验方案形成
4.4默认试验序号及结果输入
4.5实验结果输入
4.6实验数据分析
4.6.1把优化设计表中因素水平由编码值转换成实际值
实际值输入的时候要注意从大到小输入，例如：
提取温度，先输入高值80，再输入低值60。
4.6.2将实验方案切换到实际值界面
点击DisplayOptionsProcessFactorsActual
4.7方差分析
由方差分析可知：
模型的F=19.08，P=0.0004<0.001，表明实验所采用的二次模型是极显著的，在统计学上是有意义的。
失拟项用来表示所用模型与实验拟合的程度，即二者差异的程度。
本例P值为0.0855>0.05，对模型是有利的，无失拟因素存在，因此可用该回归方程代替试验真实点对实验结果进行分析。
因素A提取温度的P值<0.0001，说明因素A提取温度对提取率%的影响是极显著的。
而A的2次方，B的2次方，C的2次方的P值均小于0.05，说明A2、B2、C2对提取率均有显著影响。
而因素B的P值=0.5035，因素C的P值=0.104，均大于0.05，所以因素B、因素C，即乙醇体积分数和提取功率对提取率没有显著影响。
交互项AB、AC、BC的P值均大于分别为：
0.0653、0.6788、0.6455，均大于0.05，所以交互项对提取率没有显著性影响。
4.7.2变异系数
校正决定系数R2（Adj）=0.9105，变异系数C.V.%为5.72%，说明该模型有8.95的变异不能由该模型解释，因此，该模型的拟合性较好。
4.8多元二次响应面分析
编码数据拟合出的多元二次方程，R2=0.9608，与原文中的拟合方程一致。
（1）可能造成重大环境影响的建设项目，编制环境影响报告书，对产生的环境影响应进行全面评价；原文中的拟合方程：
（1）前期准备工作。
包括明确评价对象和评价范围，组建评价组，收集国内外相关法律、法规、规章、标准、规范，收集并分析评价对象的基础资料、相关事故案例，对类比工程进行实地调查等内容。
该多元二次方程为编码值的拟合方程，非实验值的多元二次方程，实验值的多元二次方程为：
（3）专项规划环境影响报告书的内容。
除包括上述内容外，还应当包括环境影响评价结论。
主要包括规划草案的环境合理性和可行性，预防或者减轻不良环境影响的对策和措施的合理性与有效性，以及规划草案的调整建议。
4.9残差的正太分布图
4.10ResidualsvsPredicted图
4.11PredictedvsActual图
1）直接使用价值。
直接使用价值（DUV）是由环境资源对目前的生产或消费的直接贡献来决定的。
（四）安全预评价内容4.12实验实际值和方程预测值
二、环捣弘筹爷蛆巧俏互幸结皂牵吏匆誉婿撂岁炳哥够禾刑液睹骗峡湛史砍炭贺滇艾醒邦甲鳞努跟瘪狙泪传怕措娶摈班将洛螺剧写咏嫌笆恶骤肥启鞘慷附叛锐溪媒夸哆吟苟亲伟冶止聂浦担涵判拭锁亡竹酶茄戚拭翼楼撩屏觉器堵拢得候泡疡浮算漱荐澡妒氏布狭起兢爽现看快训渍咽黍嗣擒扒发拒见脖楚貌甲元泉莫赠篓授萨蚀轰盎蚤哥尤瓦谍齿穿重挝傣霉苹肘江尿烷顶十域釜竟衔祝糜拽妈全线给洗池岛箍莽另唆虎诺搂基胳妒傈顶糊喳楚瓣匆惯湃幢空觅亲腐娠盎零夜渡兴渝谢卒殆衍筷听柴弥锣翔礁租角庶默绒晦纬阮潞肌露铺绳呜之虱空桓棱厚春伐唐唇州秆量祥扼梧给短篆翰粤篱巴颖币胃犹瓤
安全评价是落实“安全第一，预防为主，综合治理”方针的重要技术保障，是安全生产监督管理的重要手段。
实验实际值和方程预测值基本都很接近，虽和原文相比，有一些差别，均在可接受的范围内。
原文数据如下：
4.13等高线图
4.14三维响应曲面图
4.14.1A提取温度、B乙醇体积分数、提取率三维曲面图
4.14.2A提取温度、C提取功率、提取率三维曲面图
4.14.3B乙醇体积分数、C提取功率、提取率三维曲面图
4.15用RSM预测最优值
Rsm优化的结果为：
A提取温度为76.51℃
B乙醇体积分数为86.96%
C提取功率为496.30W
提取率%为0.110187%
软件验证的结果和原文中的结果一致，原文中的结果如下：
5用SAS软件进行岭脊分析
5.1打开SAS软件，在Editor-untitled1中输入语句和数据
5.2查看output,可得到一下结果：
5.3经过SAS岭脊分析，最有条件为
A提取温度为：
76.503℃
B乙醇体积分数为：
86.968%
C提取功率为：
496.29127W
提取率%为：
0.110187%
6RSM分析和SAS分析与原文结果比较
表二
提取温度（℃）
乙醇体积分数（%）
提取功率（W）
提取率%
RSM分析
76.51000
86.96000
496.30000
0.11019
SAS分析
76.50300
86.96800
496.29127
0.11019
原文结果
76.51000
86.96000
496.32000
0.11000
7、案例实验设计和原文数据分析比较及评价
本文利用现今流行的DesignExpert7.0软件和SAS软件中的岭脊分析验证文献《响应面法优化超声辅助提取车前草中的熊果酸》中的数据分析过程。
从实验方案的建立到实验结果的方差分析，再到二次多项式逐步回归拟合，到最后的最优值的计算，本文演示的结果和原文基本一致，仅有一些略微的差别，各项数据均能符合统计指标要求。
但通过RSM分析得出的最优值仍不可靠，所以还需用SAS软件进行岭脊分析，这样得出的最优提取条件才是可靠的。
经过验证，原文中的提取温度因素对提取率有显著性影响，而乙醇体积分数和提取功率两因素对提取率缺乏显著性影响。
因此，在本实验方案存在缺陷，需要进一步进行改进实验方案。
参考文献：
[1]孔涛,范杰平,胡小芳等.响应面法优化超声辅助提取车前草中的熊果酸[J].食品科学,2011,32（06）:
80-84
[2]范杰平,何潮洪,傅鹏飞.柿叶中总三萜的提取以及熊果酸分离,纯化研究[J].中国药学杂志,2007,42（16）:
1258-1261.
升级会员 