K—S函数拟合最大差值极小响应面的研究.pdfVIP

K-S函数拟合最大差值极小响应面的研究 字慧平隋允康丁力 (北京工业大学，100124) 摘要：现有的响应曲面是按最小二乘法实现的。本文则提出了一种新的响应面拟合方法： 使响应面函数同样本值之间的最大距离极小．本文提出的途径不仅丰富了响应面的构造方 法，而且可以满足实际工程问题中希望响应函数与样本值最大距离极小化的需求． 关键词：响应面方法；最大值极小；Kreisselmerier-Steinhauser函数；最小二乘法 一、概论 现有的响应面方法通常采用最小二乘法，即拟合值与样本值在各点的残差的平方和 最小来确定响应系数。本文在以往研究的基础上，提出了响应面函数距样本值最大离差 值极小化的方法。本文提出的问题属于最大极小化(Minimax)问题，可以采用K-S函 数求解。如艾剑良等…将K-S函数法运用到已知多个目标函数求解最大值最小化，而本 文则首次运用该方法确定响应面函数的系数。 二、基本原理 2．1Kreisselmerier-Steinhauser函数 对于标量函数石(x)，五(x)，…，五(x)，构造下列标量函 数： 蹦，^。…，厶)-／。In【窆。戚(町卜}min (1) 式中，p0为熵参数，K为标量函数个数，《奠．／『； x=(薯，…，邑)7，函数sty,，五，…，五)称为关于 石，五，…，厶的K-S函数。可以证明K-S函数具有以下性质： (2) 哭二 m毗)s蹦，驴。fDsmax(A)+{lnK 一● 其物理意义可用图1描述。 图1卜S函数示意图 2．2数学模型 假定响应面函数为一次函数，定义为： f～(x,f1)=flo+∑历_ (3) J=l 若响应面函数为含交叉项的二次型 萝=属+∑岛■+∑∑岛‘■ (4’ jffil s=l』4 式中反——常数项待定系数；尼——一次项待定系数；岛——二次项待定系数， 381 共Q+1)Q+2)／2+。目的就是求解出待定系数。通过代换得到： 7：Fk-I口r (5’(5) 7=∑％而 根据响应面函数与样本值的差的最大值最小，有目标函数： max(f(x‘n，口)一f‘‘’)2专min(七=l，…，K) (6) p(氕。“一卜，‘”)2，一m；n (7) 转化为：s=吉·nc砉P 令响应面函数与样本值的差值为g(x‘n，口)，于是： s=扣喜砂m卅，一m；n (8) j=o…1．．，刀 (9) 令x：=l，有：窆g(x㈣，口弦p92(xO)a)一Jk=o 本文以最dx-乘法拟合得到的ao为初始展开点，采用一阶泰勒展开，于是式(9) 在口。点展开，有： 整理后得： z(x㈣，a。)+1)a? 窆g(x㈨，a。)P腰2(t(k，ao)．，k∑n∑gJ—k