第9章 响应面设计分析.ppt

《第9章 响应面设计分析.ppt》由会员分享，可在线阅读，更多相关《第9章 响应面设计分析.ppt（157页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第9章响应面设计分析 在响应分析中，观察值y可以表述为：其中 是自变量 的函数，是误差项。在响应面分析中，首先要得到回归方程，然后通过对自变量 的合理取值，求得使 最优的值，这就是响应面设计试验的目的。确信或怀疑因素对指标存在非线性影响；确信或怀疑因素对指标存在非线性影响；因素个数因素个数2-72-7个，一般不超过个，一般不超过4 4个；个；所有因素均为计量值数据；所有因素均为计量值数据；试验区域已接近最优区域；试验区域已接近最优区域；基于基于2 2水平的全因子正交试验。水平的全因子正交试验。适用范围适用范围2023/4/1982023/4/1992023/4/19102023/4/19112

2、023/4/19122023/4/1913王玉顺：试验设计与SAS分析2023/4/1914王玉顺：试验设计与SAS分析2023/4/19152023/4/19162023/4/19172023/4/19182023/4/19192023/4/19202023/4/1921王玉顺：试验设计与SAS分析2023/4/19222023/4/19232023/4/19242023/4/19252023/4/19262023/4/19272023/4/1928王玉顺：试验设计与SAS分析2023/4/19292023/4/19302023/4/1931响应面设计与分析响应面设计与分析2023/4/19

3、32(1)响应面模型9 响应面设计与分析响应面模型是多元二次多项式2023/4/1933(1)响应面模型 连续定量的变量间关系，可用下面的二次响应面模型描述：10 响应面设计与分析2023/4/1934王玉顺：试验设计与SAS分析(1)响应面模型 考察第i次试验，用响应面模型表述响应变量y与自变量(x1,x2,xp)样本观测间的关系：10 响应面设计与分析2023/4/1935 所有n个变量观测间的关系，用矩阵形式的二次响应面模型表述如下：(1)响应面模型设计矩阵10 响应面设计与分析2023/4/1936(1)响应面模型响应向量回归参数向量10 响应面设计与分析2023/4/1937(1)响

4、应面模型残差向量响应面模型10 响应面设计与分析2023/4/1938(2)响应面设计响应面设计(Response Surface Design)是一种探析变量间相关关系的试验设计方法，亦称回归设计。它从控制响应精度分布特性的角度出发，着眼于实现诸多优良统计特性，设计并优化变量的水平组合(处理)或试验点。与其它试验设计方法不同，它以“果”的要求设计“因”的水平组合(处理)，而不是以效应要求设计“因”的水平组合。10 响应面设计与分析2023/4/1939(2)响应面设计 响应面设计的特点：(1)正交化设计，试验布点均衡分散，代表性强；(2)统计控制响应的精度，优化试验的响应特性；(3)力图实现

5、尽可能少的试验次数；(4)实现能对因素效应进行独立解释的试验方案；(5)实现能求解最佳试验处理的试验方案。响应面设计最终可归结为对设计矩阵的优化设计，按试验点的分布及特性，可分为正交设计、中心组合设计、旋转设计和最优设计等。10 响应面设计与分析2023/4/1940(3)响应面分析 响应面分析(Response Surface Analysis)是与响应面设计配套的数据分析方法,主要任务如下：(1)回归方程估计和检验；(2)回归参数估计和检验；(3)模型拟合不足检验；(4)因素效应检验；(5)一次项、二次项及交叉项模型效应检验；(6)典型分析：求解响应面驻点及其响应；10 响应面设计与分析2

6、023/4/1941王玉顺：试验设计与SAS分析(7)岭脊分析：它以中心试验点为原点，在试验范围内以到原点的距离为半径搜索响应面上的最大或最小响应点，最后获得由一系列距离及其最佳响应点构成的集合，可据此选择需要的最佳水平组合；(8)等值线分析：垂直于不同坐标轴剖切响应面，获得响应面的等值线，可从不同剖面观察试验的响应特性。(3)响应面分析10 响应面设计与分析2023/4/1942王玉顺：试验设计与SAS分析本章内容本章内容10.1 一般性问题10.1.1 案例和解决方案10.1.2 设计变量和响应变量10.1.3 设计点类型和空间分布10.1.4 设计点类型和数目10.1.5 设计编码和归一

7、化编码10.1.6 试验方案和设计矩阵10 响应面设计与分析2023/4/1943王玉顺：试验设计与SAS分析本章内容本章内容10.2 响应面设计10.2.1 因素水平编码10.2.2 中心组合正交设计10.2.3 中心组合可旋转设计10.2.4 中心组合正交旋转设计10.2.5 Box-Behnken Design10.2.6 Hybrid Design 10.2.7 设计方案的响应精度10.2.8 响应面设计的SAS实现10 响应面设计与分析2023/4/1944本章内容本章内容10 响应面设计与分析10.3 响应面分析10.3.1 数据整理10.3.2 SAS编程10.3.3 简单统计1

8、0.3.4 响应面回归分析10.3.5 响应面典型分析10.3.6 响应面岭脊分析10.3.7 响应面图形分析2023/4/194510.1 一般性问题Universality Problems10 响应面设计与分析2023/4/1946(1)案例和解决方案10.1 一般性问题 响应面设计与分析的过程陈述，将通过下面的案例展开。案例案例：某温室废弃物制沼气试验，试验因素为配料浓度(ferment)、PH值(PH)和配料比(scale)，检测指标为产气量(gasvalue)和甲烷含量(CH4)。试完成试验设计、试验实施和试验结果分析，找出最佳处理，以提高产气量为目标确定下一步试验的方向。2023

9、/4/194710.1 一般性问题 解决问题选用响应面设计和分析，主要依据是：(1)与线性模型相比，响应面模型是二次曲面，它能更逼近地拟合实际问题中的变量间关系；(2)与其它设计相比，响应面设计是较强效应检验和较小试验方案的完美结合；(3)典型分析确定响应面的驻点，并判断是鞍点还是最优点，为确定最佳处理提供信息；(1)案例和解决方案2023/4/1948(4)岭脊分析给出从中心点出发到试验边界终止的最优点轨迹（岭脊），为确定最佳处理提供多种选择；(5)网格图和等值线图展示响应面形状、变化规律和变化趋势，为搜寻最佳处理的下一步试验提供方向性信息。(1)案例和解决方案10.1 一般性问题2023/

10、4/1949(2)设计变量和响应变量10.1 一般性问题 明确试验目的及确定设计变量和响应变量是首先要解决的问题。首先要理清欲考察谁对谁的影响或谁对谁的关系，从机理层面上分析因果关系，确定哪些是自变量，哪些是因变量。根据试验目的，发酵浓度、pH值及原料配比应是设计变量(试验因素)，沼气产量和甲烷含量应是响应变量(试验指标)，欲研究因素对指标的影响程度、变量间关系及关系密切程度。2023/4/1950(3)设计点类型和分布10.1 一般性问题 因素的水平组合称作处理，亦称试验点或设计点，用几何形式表述主要有析因设计点、主轴设计点和中心设计点三种类型，还有棱边设计点和其它设计点等。2023/4/1

11、951(4)设计点类型和数目10.1 一般性问题自变量个数p，水平数a，回归参数个数q，试验容量n,析因点数np,主轴点数n,中心点数n0ap-k部分析因设计点、主轴设计点和中心设计点2023/4/1952 试验设计若满足基本要求，则全部回归估计和假设检验可进行；否则，只有部分回归估计和假设检验可进行。响应面设计一般可能存在下述三种情形：(1)nq，称作过饱和设计。此时若三种设计点齐全，则可进行全部回归估计和假设检验。若中心设计点只重复1次，则不能进行拟合不足检验(lack fit test)；若没有主轴设计点，则多数二次项回归系数不能估计和检验。析因设计点主轴设计点中心设计点2023/4/1

12、954(5)设计编码10.1 一般性问题2023/4/1955(6)归一化编码10.1 一般性问题2023/4/1956(7)试验方案和设计矩阵10.1 一般性问题 响应面设计的试验方案和设计矩阵可用编码矩阵表达。为便于理解设计原理，后面的讨论将以三因素二次响应面设计为例。由编码书写的三因素响应面中心组合设计的试验方案如右面所示：2023/4/195710.1 一般性问题(7)试验方案和设计矩阵设计矩阵X:2023/4/195810 响应面设计与分析10.2 响应面设计Response Designs2023/4/195910.2.1 因素水平编码10.2.2 中心组合正交设计10.2.3 中

13、心组合可旋转设计10.2.4 中心组合正交旋转设计10.2.5 Box-Behnken Designs10.2.6 Hybrid Designs 10.2.7 设计方案的响应精度10.2.8 响应面设计的SAS实现10.2 响应面设计本节内容本节内容2023/4/1960fermentPHscale设计码归一码x-j450.2-1xminj5.625.81.2-1-x0j872.600 xmaxj10.388.24+1+x j1295+1Sj2.381.21.4=1.681793表10.1 温室废弃物制沼气试验的因素水平编码表(1)因素水平编码10.2 响应面设计2023/4/1961(2)中

14、心组合正交设计10.2 响应面设计 中心组合正交设计(Orthogonal Central Composite Design)，“正交”指设计矩阵可正交化，“中心”指中心设计点，“组合”指试验方案由三种设计点构成。响应面设计的正交性具体指设计矩阵的可正交化，即设计矩阵具有能够变换为正交矩阵的性质，它是一个优良性质，在回归估计和分离效应时很有用，它还可提高估计精度和数值计算稳定性。用设计矩阵可正交化确定中心点的个数。析因设计点主轴设计点中心设计点2023/4/1962(2)中心组合正交设计10.2 响应面设计正交条件2023/4/1963(2)中心组合正交设计10.2 响应面设计Orthogon

15、alCentral Composite Design正交条件选定23析因设计点+2p主轴设计点+4中心设计点2023/4/1964(3)中心组合可旋转设计10.2 响应面设计 中心组合可旋转设计(Rotatable Central Composite Design)。所谓“可旋转”指位于同一球面上的所有试验点的响应方差相同。调整值，就可实现二次响应面中心组合设计的可旋转性，值仅与析因设计点的个数有关，值称作轴的缩放比例(axial scaling)。旋转条件2023/4/196510.2 响应面设计(3)中心组合可旋转设计精度均衡中心组合设计Uniform Precision Central

16、Composite Design 23析因设计点+2p主轴设计点+3中心设计点2023/4/196610.2 响应面设计(3)中心组合可旋转设计精度均衡中心组合设计Uniform Precision Central Composite Design 2023/4/196710.2 响应面设计(3)中心组合可旋转设计小试验量中心组合设计Small Central Composite Design23-1部分析因设计点+2p主轴设计点+5中心设计点2023/4/196810.2 响应面设计(3)中心组合可旋转设计小试验量中心组合设计Small Central Composite Design 20

17、23/4/196910.2 响应面设计(4)中心组合正交旋转设计 同时满足正交性和可旋转性的中心组合设计，称作中心组合正交旋转设计(Orthogonal and Rotatable Central Composite Design)。设计的基本步骤是：(1)调节 值实现旋转设计；(2)通过调整中心设计点个数n0决定试验容量n实现正交设计。注意旋转或正交设计单独均易实现，同时满足两者设计的试验容量可能不是整数，需要圆整。2023/4/197010.2 响应面设计(4)中心组合正交旋转设计2023/4/197110.2 响应面设计(4)中心组合正交旋转设计中心组合正交旋转设计Orthogonal

18、and Rotatable Central Composite Design2023/4/197210.2 响应面设计(5)Box-Behnken Design Box-Behnken Design，以立方体的12条棱边的中心作析因设计点，取消主轴点，选较少中心点，实际上是3水平正交设计。该设计的特点是，试验次数较少，均衡分散，代表性较好，响应回归精度均衡，只能进行回归估计和部分检验。2023/4/197310.2 响应面设计(5)Box-Behnken Design2023/4/197410.2 响应面设计(6)Hybrid Design Hybrid Design是一种具有最少试验次数的响

19、应面设计方法，设计点的水平编码既有大于或小于“1”的情况又有大于或小于“-1”的情况，仅选用两个对称或不对称的主轴点，有1个或没有中心点。设计点主要考虑均衡分散和代表性，设计点的响应回归精度均衡，分布半径范围或收缩或伸张以满足不同试验目的。该设计能进行全部回归估计，能部分或不能进行假设检验，一般用于筛选处理的探索性试验。2023/4/197510.2 响应面设计(6)Hybrid Design 2023/4/197610.2 响应面设计(6)Hybrid Design 2023/4/197710.2 响应面设计(6)Hybrid Design 2023/4/1978OBSR310R310AR3

20、10BSCCDBBDUPCCD ORCCD11.00000.75000.80000.90740.75000.66980.669721.00000.75000.80000.90740.75000.66980.669731.00000.93750.92500.90740.75000.66980.669741.00000.93750.92500.90740.75000.66980.669751.00000.93750.92500.90740.75000.66980.669761.00000.93750.92500.90740.75000.66980.669771.00000.93750.92500.

21、90740.75000.66980.669781.00000.93750.92500.90740.75000.66980.669791.00000.93750.92500.90740.75000.60730.6072101.00000.93750.92500.90740.75000.60730.6072111.00001.00000.18520.75000.60730.6072120.18520.75000.60730.6072130.18520.33330.60730.6072140.18520.33330.60730.6072150.18520.33330.16630.1110160.16

22、630.1110170.16630.1110180.16630.1110190.16630.1110200.16630.1110210.1110220.1110230.1110(7)各个设计方案的响应精度比较2023/4/197910.2 响应面设计(8)响应面设计的SAS实现Step1点击SolutionsAnalysisDesign of Experiments出现ADX试验设计与分析试验设计与分析界面Analysis and Design of Experiments2023/4/198010.2 响应面设计(8)响应面设计的SAS实现响应面设计按钮Step2点击Cerate Respo

23、nse-Surface Design按纽出现New Design窗口2023/4/1981王玉顺：试验设计与SAS分析10.2 响应面设计(8)响应面设计的SAS实现定义变量按钮Step3点击New Design窗口的右侧按钮Define Varibles2023/4/1982王玉顺：试验设计与SAS分析10.2 响应面设计(8)响应面设计的SAS实现Step4点击Factor页的右侧按钮Add，选，选32023/4/1983王玉顺：试验设计与SAS分析10.2 响应面设计(8)响应面设计的SAS实现Step5点击Factor页下方中部的OK按钮，选yes2023/4/1984王玉顺：试验设计

24、与SAS分析10.2 响应面设计(8)响应面设计的SAS实现Step6在New Design页点击Select Design按钮在Design List Options页选定响应面设计方案2023/4/1985王玉顺：试验设计与SAS分析10.2 响应面设计(8)响应面设计的SAS实现Step7点击Design Details按钮浏览Design Information页和Design List页2023/4/1986王玉顺：试验设计与SAS分析10.2 响应面设计(8)响应面设计的SAS实现Step8关闭Design Details页，点击yes保存设计New Design窗口显示完成的试验

25、设计方案2023/4/1987王玉顺：试验设计与SAS分析10.2 响应面设计(8)响应面设计的SAS实现Step9点击New Design窗口的Report进入ADX:Report窗口在Design Items列表框中选择Design Points项点击Generate Report在ADX Repot窗口输出设计方案 2023/4/1988王玉顺：试验设计与SAS分析10.2 响应面设计(9)响应面设计注意事项 就响应面设计而言，要注意平衡试验者各种矛盾的需求，在试验目的、试验次数、试验效率、试验成本、响应精度及响应特性之间取得平衡。筛选因素或寻求最佳处理的探索性试验，可选择具有较少试验次

26、数的Hybrid设计、BOX设计或小型中心组合设计。探明变量关系及其规律的分析性试验，可选择均衡精度或正交可旋转的中心组合设计。选择试验次数较少的设计必会以削弱精度和损失部分检验项目为代价。响应面设计，实质上是对设计矩阵的设计，可采用编码方法使设计变得简便易行。2023/4/198910 响应面设计与分析10.3 响应面分析Response Analysis2023/4/199010.3 响应面分析响应面分析的实现 根据响应面设计方案实施试验和获得试验结果，创建SAS数据表，用试验数据对响应面模型进行回归估计、效应检验、典型分析、岭脊分析等，称作响应面分析。响应面分析，利用 RSREG 响应面

27、分析过程编程实现，也可利用ADX界面的菜单操作实现。2023/4/199110.3.1 数据整理10.3.2 SAS编程10.3.3 简单统计10.3.4 响应面回归分析10.3.5 响应面典型分析10.3.6 响应面岭脊分析10.3.7 响应面图形分析本节内容本节内容10.3 响应面分析2023/4/199210.3 响应面分析10.3.1 数据整理Data Management2023/4/1993试验号编码设计方案因素设计方案 响应值fermentPHscalefermentPHscaleGasvalueCH41-1-1-15.625.81.219.443.32-1-1+15.625.8

28、418.150.93-1+1-15.628.21.225.347.54-1+1+15.628.242748.95+1-1-110.385.81.236.744.86+1-1+110.385.8425.647.37+1+1-110.388.21.233538+1+1+110.388.2430.947.79-1.681800472.618.748.610+1.6818001272.635.747.5110-1.68180852.622.749.0120+1.68180892.628.450.01300-1.6818870.232.951.11400+1.681887528.651.51500087

29、2.624.747.916000872.626.548.217000872.629.851.618000872.628.249.219000872.627.848.520000872.626.749.321000872.625.250.322000872.626.350.423000872.624.949.3表10.2 温室废弃物制沼气试验的响应面设计方案和试验结果2023/4/1994fermentPHscalegasvalueCH45.625.81.219.443.35.625.8418.150.95.628.21.225.347.55.628.242748.910.385.81.236.

30、744.810.385.8425.647.310.388.21.2335310.388.2430.947.7472.618.748.61272.635.747.5852.622.749892.628.450870.232.951.187528.651.5872.624.747.9872.626.548.2872.629.851.6872.628.249.2872.627.848.5872.626.749.3872.625.250.3872.626.350.4872.624.949.3表10.3 温室废弃物制沼气试验的SAS数据表(sasuser.gasch4)2023/4/199510.3 响

31、应面分析10.3.2 SAS编程Programing by SAS2023/4/1996(1)SAS编程实现响应面分析10.3.2 SAS编程title 中心组合正交旋转设计响应面分析;proc rsreg data=sasuser.gasch4;model CH4 gasvalue=ferment PH scale/lackfit press;ridge max min center=8 7 2.6 outr=aa;run;quit;问题分析：响应面模型是多元二次多项式，SAS提供了一个专门按此模型进行响应面分析的过程rsregrsreg，主要功能是参数估计和检验、模型（回归方程）检验、回归

32、分项检验、失拟检验、因素效应检验、典型分析、岭脊分析、决定系数及其它统计量的计算。建模语句model CH4 gasvalue=ferment PH scale为两个响应变量CH4和gasvalue分别指定一个响应面模型，等号右端只需列出自变量，变量相互之间用空格隔开；建模语句后面的两个选项lackfit和press分别指定失拟检验和输出残差平方和。岭脊分析语句ridge max min center=8 7 2.6 outr=aa中，关键词max指定计算最大响应岭脊点；关键词min指定计算最小响应岭脊点；关键词center=8 7 2.6指定岭脊搜索中心，一般用试验设计的中心点，若数据表中的

33、因素水平值是编码数据，则关键词为center=0 0 0；关键词outr=aa指定将岭脊分析结果存盘为临时库的数据文件aa。2023/4/1997(2)SAS编程生成响应面绘图数据/*计算并创建绘响应面图的数据表GasValue*/data GasValue;x2=9;/*x2分别取5,7,9三个水平*/do x1=4 to 12 by 0.1;do x3=0.2 to 5 by 0.1;z=-40.1388+7.3477*x1+9.6052*x2-6.4517*x3-0.5777*x1*x2-0.5102*x1*x3+0.8929*x2*x3+0.0016*x1*2-0.4080*x2*2+

34、0.6241*x3*2;output;end;end;run;10.3.2 SAS编程2023/4/1998/*从绘图数据表中选择部分数据创建网格图数据表gridfile*/proc g3grid data=GasValue out=gridfile;grid x3*x1=z/naxis1=20 naxis2=20;run;(3)SAS编程生成网格图数据表10.3.2 SAS编程2023/4/1999/*绘图设置*/goptions reset=all ftext=swiss htext=1.5 colors=(black r b black r b black r b black r b);

35、axis1 label=(f=宋体 c=black h=15pt 配料浓度()c=blue width=1.5;axis2 label=(A=90 f=宋体 c=black h=15pt 配料比)c=blue width=1.5;symbol height=1.5 width=2.0;(4)SAS编程实现绘图设置10.3.2 SAS编程2023/4/19100/*根据数据表gridfile绘网格图*/proc g3d data=gridfile;plot x3*x1=z/caxis=blue ctext=black cbottom=red ctop=black tilt=70 rotate=7

36、0 xticknum=7 yticknum=7 zticknum=7 grid;run;/*根据数据表GasValue绘等值线图*/proc gcontour data=GasValue;plot x3*x1=z/levels=23 to 37 by 2 autolabel nolegend haxis=axis1 vaxis=axis2 caxis=b grid;run;quit;(5)SAS编程绘制网格图和等值线图10.3.2 SAS编程2023/4/1910110.3 响应面分析10.3.3 简单统计Simple Statistic2023/4/19102王玉顺：试验设计与SAS分析(1

37、)因素水平编码系数表10.3.3 简单统计 编码系数表，列出了因素水平归一化编码运算的减数和除数，因素水平减去减数再除以除数就等于编码值，最大值1，最小值-1，该编码用于参数估计、典型分析和岭脊分析。-The RSREG ProcedureCoding Coefficients for the Independent Variables Factor Subtracted off Divided by ferment 8.000000 4.000000 PH 7.000000 2.000000 scale 2.600000 2.400000-2023/4/19103(2)简单统计量的计算结果R

38、esponse Surface for Variable CH4:CH4 Response Mean 48.947826 Root MSE 1.438121 R-Square 0.7398 Coefficient of Variation 2.9381 Sum of Squared Residuals 26.886487802 Predicted Residual SS(PRESS)133.11245810.3.3 简单统计2023/4/19104王玉顺：试验设计与SAS分析(2)简单统计量的计算结果Response Surface for Variable gasvalue:gasvalue

39、 Response Mean 27.091304 Root MSE 1.548607 R-Square 0.9392 Coefficient of Variation 5.7163 Sum of Squared Residuals 31.176382471 Predicted Residual SS(PRESS)99.50462317810.3.3 简单统计2023/4/1910510.3 响应面分析10.3.4 响应面回归分析Regression Analysis of Response Surface2023/4/19106(1)回归参数估计和检验10.3.4 响应面回归分析表10.4 C

40、H4响应面参数估计及其t检验 ParameterDFEstimatet ValuePr|t|Coded EstimateIntercept111.0629420.620.543049.443717Ferment11.5512980.930.36900.043252PH14.4414171.130.27841.532787Scale110.3417393.930.00170.850588Ferment*Ferment1-0.146894-2.300.0384-2.350309PH*Ferment10.2801121.570.13962.240896PH*PH1-0.229113-0.900.38

41、37-0.916453Scale*Ferment1-0.442677-2.900.0124-4.249700Scale*PH1-1.041667-3.440.0044-5.000000Scale*Scale10.1626450.910.37880.9368372023/4/19107(1)回归参数估计和检验ParameterDFEstimate t Value Pr|t|Coded EstimateIntercept1-40.1388-2.100.055326.6941Ferment17.34474.090.00138.0001PH19.60522.270.04073.1844Scale1-6

42、.4517-2.280.0403-2.4918Ferment*Ferment10.00160.020.98150.0259PH*Ferment1-0.5777-3.010.0100-4.6218PH*PH1-0.4080-1.490.1598-1.6322Scale*Ferment1-0.5102-3.100.0084-4.8980Scale*PH10.89292.740.01694.2857Scale*Scale10.62413.250.00643.5946表10.5 GasValue响应面参数估计及其t检验10.3.4 响应面回归分析2023/4/19108王玉顺：试验设计与SAS分析(2

43、)响应面模型及其分项检验表10.6 CH4响应面模型及其分项检验RegressionDFTypeSSR SquareF ValuePrFLinear314.8908150.14412.400.1148Quadratic314.5150890.1405 2.340.1210Crossproduct347.0250000.4552 7.580.0035Total Model976.4309040.73984.11 0.010910.3.4 响应面回归分析2023/4/19109王玉顺：试验设计与SAS分析RegressionDFTypeSSR SquareF ValuePrFLinear3388.

44、0729480.756553.94.0001Quadratic330.8689300.06024.290.0260Cross product362.9000000.12268.740.0020Total Model9481.8418780.939222.32 FLack of Fit515.6775993.1355202.240.1489Pure Error811.2088891.401111Total Error 1326.8864882.068191表10.8 CH4响应面模型失拟检验10.3.4 响应面回归分析2023/4/19111王玉顺：试验设计与SAS分析ResidualDFSum

45、 of SquaresMean SquareF Value Pr FLack of Fit58.4208271.6841650.590.7078Pure Error822.7555562.844444Total Error 1331.1763822.398183表10.9 GasValue响应面模型失拟检验(3)响应面模型拟合不足检验10.3.4 响应面回归分析2023/4/19112王玉顺：试验设计与SAS分析(4)响应面因素效应检验FactorDFSum of SquaresMean SquareF Value Pr FFerment433.5102168.3775544.050.0239

46、PH442.76596610.6914925.170.0103Scale447.03835011.7595875.690.0072表10.10 CH4响应面因素效应检验10.3.4 响应面回归分析2023/4/19113王玉顺：试验设计与SAS分析FactorDFSum of SquaresMean SquareF ValuePr FFerment4354.16853888.54213436.92.0001PH494.59684523.6492119.860.0007Scale495.71488923.9287229.980.0006表10.11 GasValue响应面因素效应检验(4)响应面

47、因素效应检验10.3.4 响应面回归分析2023/4/19114王玉顺：试验设计与SAS分析式中：z1甲烷含量；z2产气量；x1配料浓度；x2PH值；x3配料比(5)响应面方程10.3.4 响应面回归分析2023/4/19115(6)回归分析结论 回归分析结果表明，产气量和甲烷含量的响应面模型决定系数分别为0.9392和0.7398，显著性水平达0.0001和0.0109；产气量响应面模型的线性项、平方项和交叉项的效应均显著，显著性水平分别达0.0001、0.0260和0.0020；甲烷含量响应面模型只有交叉项效应显著，显著性水平达0.0035；拟合不足检验均不显著，说明模型拟合适当。因素效应

48、检验均显著，尤以对产气量的效应极强，显著性水平最大值达0.0007。由此可见，产气量与配料浓度、PH值和配料比存在较密切的关系，而甲烷含量与三个因素的关系相对较松散，响应面回归有效，模型可用。10.3.4 响应面回归分析2023/4/1911610.3 响应面分析10.3.5 响应面典型分析Canonical Analysis of Response Surface2023/4/19117稳定点(驻点)向量：响应面方程稳定点(驻点)响应：(1)确定响应面稳定点的方法10.3.5 响应面典型分析2023/4/19118王玉顺：试验设计与SAS分析矩阵A特征值矩阵A性质稳定点性质均取负值负定极大值

49、点均取正值正定极小值点有正有负不定鞍点含有0值在一平面上稳定点性质与矩阵A特征值和矩阵A性质的关系(2)判断稳定点性质的方法10.3.5 响应面典型分析2023/4/19119王玉顺：试验设计与SAS分析(3)CH4典型分析结果The RSREG ProcedureCanonical Analysis of Response Surface Based on Coded Data for Variable CH4Critical ValueFactor Coded Uncoded Labelferment -0.034634 7.861464 ferment PH 0.271341 7.542

50、682 PH scale 0.191567 3.059760 scalePredicted value at stationary point:49.73239510.3.5 响应面典型分析2023/4/19120王玉顺：试验设计与SAS分析EigenvectorsEigenvalues ferment PH scale 3.629175 0.371042 0.516053 -0.772021-2.557687 -0.313999 0.852120 0.418684-3.401413 0.873918 0.087065 0.478213Stationary point is a saddle