等高图法在单纯形格子设计最优条件选择中的应用

山西医科大学本科生毕业论文基于单形格子混料设计的等高线寻优学生姓名王晓荣学号20060494指导教师仇丽霞专业预防医学班级306032011年6月目录中文摘要……………………………………………………………………………1英文摘要……………………………………………………………………………2前言…………………………………………………………………………………31原理与方法………………………………………………………………………41.1单纯形格子……………………………………………………………………51.1.1单纯形格子的原理.............................................61.1.2单纯形格子的模型建立.........................................71.2等高图法…………………………………………………………………………61.2.1等高图法的原理………………………………………………………………71.2.2等高图法寻优的步骤…………………………………………………………72实例分析……………………………………………………………………………72.1实例基本介绍……………………………………………………………………2.2实例结果…………………………………………………………………………2.2.1模型建立方法…………………………………………………………………2..2.2等高图法搜索最优配方…………………………………………………………………………3讨论………………………………………………………………………………104结论………………………………………………………………………………11参考文献……………………………………………………………………………12致谢…………………………………………………………………………………13等高图法在单纯形格子设计最优条件选择中的应用【摘要】目的了解单纯形格子实验设计的原理、方法、模型建立；应用等高图法寻找单纯形格子试验设计的最优条件；探讨等高图法的不足之处。方法认真阅读《实验优化设计与分析》，了解和掌握了单形格子实验设计的原理、方法以及模型建立，并通过查阅CNKI搜索有关单形格子实验设计的文章，例如《冲积土区域绿竹氮磷钾施肥反应及其最优施肥参数的研究》，并利用此文献中单形格子实验设计的数据，建立回归模型做为目标函数，用等高图法搜索最优试验条件。结果山地绿竹配方施肥三分量二阶多项式回归模型有统计学意义，等高图法确定的最优试验条件：结论以三分量二阶多项式回归模型为目标函数，利用等高图法确定的试验最优条件主观性强、准确度低，应用范围有限，只有当某一变量被选为固定值,才能做出其余两个变量与响应值的等高图,从而确定最优解,由于重要变量被锁定，只选择了一些次要变量进行优化，认识只能在低水平上循环，难以有大的突破。当单纯形格子实验设计的影响因素在3个以上时．应用等高图法将其统计模型优化分析却比较复杂，有待于进一步研究。【关键词】等高图法；单纯形格子试验设计（Simplex-latticedesign）；最优试验条件Useofthecontourplotinselectingtheoptimalconditionforsimplex-latticedesign【Abstract】objectiveToexploretheuseofcontourplotinselectingtheoptimalconditionforsimplex-latticedesign.MethodsUsingdatafromthefertilizerexperimentsofN,PandKconductedforsinocalamueoldhamiplantedinalluvionsoilarea,optimalexperimentalconditionwassearchedbythecontourplotbasedonthequadraticregressionmodelasatargetfunction.ResultsTheoptimalconditionsofextractionprocesswereasfollow:.ConclusionsTheoptimalexperimentalconditionsearchedbythecontourplotwasnotmorepreciseandobjective.【Keywords】Thecontourplot;Simplex-latticedesign;Thebestexperimentalcondition前言混料问题是近三十年来新发展起来的十分活跃的统计学分支，自1955年Claringbold到目前为止，许多统计学家提出适合于不同数据结构或各种统计目的所需要的混料模型【2】。Scheffe1958年创造了单纯形格子混料设计，1963年提出单纯形重心设计等多项式混料模型。这两种试验设计有三个突出的优点：(1)试验处理少，属饱和最优设计；(2)统计方法特别简单．由人工计算即可完成；(3)拟合精度高，理论值与实测值吻合得好．由于这三个突出的优点，使得这两种试验设计方法在配方研究上非常实用．Kenworthy和Snee等人使用了混料分量之比的一般多项式混料模型。为了研究当一个或多个分量的值趋于零(即趋于利益区域的边界)时，响应值产生的急剧变化，Draper和John在1977年首先在混料多项式(格子或中心多项式)中加入某些或全部分量的倒数项1/Xi(i＝1，2，，q)，提出具有倒数项的多项式混料模型。1986年朱伟勇等人提出带有对数项的多项式混料模型。混料试验设计是各组分之和为1的特殊试验设计。它主要用于配方或用料配比的研究．即研究各成分因子最佳含量的百分比．混料设计的理论和它的应用都有发展，人们针对各种数学模型、试验区域与各种意义下的撟钣判詳提出了各种设计方法与分析计算法.本文主要介绍单纯形格子设计.．单纯形格子设计是混料回归设计方案中最先出现的，也是最基本的设计方案，将试验点取在相应阶数的正规单纯形格子点上，这样的试验设计称为单纯形格子设计。Scheffe于1958创造了单纯形格子点设计，其方法如下：先确定一个正整数m，然后让每个原料取值例如当s=3，m=1时，只有3个试验点：（1，0，0），（0，1，0），（0，0，1），当s=3，m=2时，有6个试验点：（1，0，0），（0，1，0），（0，0，1），（1/2，1/2，0），（1/2，0，1/2），（0，1/2，1/2），当s=3，m=3时，有10个试验点：（1，0，0，），（0，1，0），（0，0，1），（1/3，2/3，0），（1/3，0，2/3），（0，1/3，2/3），（2/3，1/3，0），（2/3，0，1/3），（0，2/3，1/3），（1/3，1/3，1/3），一般记为{s，m}设计，一个{s，m}设计有个试验点。单纯形混料设计方法规范，试验点均匀、对称，适用于因素、水平不太多的混料场合，当因素或水平较多时，试验次数增加很快，此时应选择其它几种混料试验设计方法。对称单纯形设计的探索法用于解决3个自变量，出现摰憔奂瘮的试验设计场合效果很好，试验点比极端顶点设计还要少，而且均匀性很好。由于该法是使用探索法，在自变量多时，设计出适用的对称单纯形设计也不是一件易事。所以该法有很大的局限性。等高线图解分析法起源于1953年Loewe提出的等效线图，但是一直没有引起关注，直到1970年，Gessner等在研究水合氯醛和乙醇的合并应用时，发现不同水合氯醛和乙醇的比例产生不同的效应，有的具有相加作用，有的具有拮抗作用，有的具有协同作用，才引起人们的兴趣。但是这些数据虽然在药效图上具有明显的不同作用的分布，却没有显著的统计学差异，这些剂量,效应数据内在的误差意味着仅仅是“高于”或“低于”相加线是不够的，还缺少精确性。当等效线图加上统计分析时，被称为评价药物相互作用的“黄金标准”，1原理及方法1.1单纯形格子1.1.1单纯形格子的原理所谓单形格子设计，就是在给定的混料回归模型下，将试验点取在相应阶数的正单形格子点上的混料设计。单形是指顶点数与坐标空间维数相等的凸图形，如正三角形、正四面体等。p维单形即（1p）维单纯形。在单形混料设计中，各分量jx的变化范围可由高为1的p维正单形表示，顶点代表单一成分组成的混料，棱上的点代表2种成分的混料，面上的点代表多于2种而少于p种成分组成的混料，而单形内的点则代表全部p种成分组成的混料。对于p分量混料实验，p维正单形的p个顶点分别为：(1,0,0A1…,0),(0,1,0,A2…,0),…,pA(0,0,…,0,1)。设R,,(21xx…，)px为单形的点，定义,1x…，px分别表示R点到2A…PA面的距离，…，R点到1A…1PA面的距离。所以，满足混料条件的有序数组,,(21xx…，)px，总可以在单形上找到一点与之对应。反之，在单形上的任意一点R,,(21xx…，)px都满足混料条件。单形格子混料设计需满足：)3,2,1(,1,0321jxxxxj。在单形格子设计中，p维d阶格子点集﹛p，d﹜中有ddp1C个点，正好与所采用的d阶完全多项式回归模型中待估计参数的个数相等，故单形格子设计是饱和设计，是在“实验次数最少”意义下的最优设计。1.1.2单纯形格子的模型建立一般情况下，混料回归设计的p分量d次多项式回归方程常采用Scheffe多项式形式，当试验设计为{P，1}单形格子设计时．Scheffe多项式形式为线性规划问题；当试验设计为P因子2次以上单形格子时，模型(2)为非线性规划问题4。如：二次式（d=2）yˆ=jhjhhjjpjjxxbxb1(式2)不完全三次式（d=3）yˆ=jhjhhjjpjjxxbxb1+kjhkjhhjkxxxb（式3）完全三次式（d=3）yˆ=jhjhhjjpjjxxbxb1+)(jhjhjhhjxxxxr+kjhkjhhjkxxxb（式4）对于Scheffe多项式回归模型，当采用单形格子设计时，每个回归系数的值只取决于按一定规律对应的一些格子点上的实验值，而与其它设计点上的实验值无关，各回归系数都可表示成相应设计点上的实验值的简单线性组合。1.2等高图法1.2.1等高图法的原理等高线是地面高程相等的相邻各点连成的闭合曲线。如图A-2所示，设想有一座高出水面的小岛，与某一静止水面相交形成的水涯线为闭合曲线，曲线形状由小岛与水面相交的位置确定，曲线上各点的高程相等。例如，当水面高为70m时，曲线上任一点的高程均为70m；若水位继续升高至80m、90m，则水涯线的高程分别为80m、90m。将不同高程的水涯线垂直投影到水平面H上，按一定比例尺缩绘在图纸上，就可将小岛用等高线表示在地形图上。这些等高线的形状和高程，客观地反映了小岛的空间形态。这便是等高线的绘制原理。图A-2等高线的绘制原理1.2.2等高图法寻优的步骤为了使系统达到最优的目标所提出的各种求解方法称为最优化方法。最优化方法解决问题一般可以分为以下几个步骤：（1）提出需要进行最优化的问题，开始收集有关资料和数据；（2）建立求解最优化问题的有关数学模型，确定变量，列出目标函数和有关约束条件；（3）分析模型，选择合适的最优化方法；（4）求解方程。一般通过编制程序在电子计算机上求得最通过上述五个相互独立和互相渗透的步骤优解；（5）最优解的验证和实施。，最终求得系统的最优解。2实例分析2.1实例基本介绍引用文献《冲积土区域绿竹氮磷钾施肥反应及其最优施肥参数的研究》中实验设计数据进行分析。在冲积土区域绿竹氮磷钾施肥反应及其最优施肥参数的研究，绿竹的经济收益越高越好，每配方施肥经济收益按笋重（kg）*3.0元/kg+新竹重（kg）*0.12元/kg计。影响绿竹经济收益的因素主要有：专用复合肥配方中氮、磷、钾的含量，选择绿竹施肥{3，2}设计表进行单形格子混料设计，在每个实验条件下重复2次，施