解非线性互补问题的约束积分水平集算法.pdfVIP

维普资讯 2004年 l2月 应用数学与计算数学学报 第 18卷第2期 Dec．，2004 C0MM．oNAPPL．MATH．ANDCOMPUT Vo1．18 NO．2 解非线性互补问题的约束积分水平集算法 俞武扬 邬冬华 上海大学理学院，上海 200436 摘要：本文考虑有约束的非线性互补问题的全局最优化问题，在文 [1】[5】的基础上， 利用数论中一致分布佳点集列，给出了以数论方法代替Monte~Caclo投点的实现算法，并 证明了所给实现算法的全局收敛性．最后所给出的两个数值例子表明本算法对求非线性互 补问题是有效的． 关键词：约束总极值；积分一水平集；非线性互补问题；一致分布佳点集 1．引 言 求解非线性互补问题即寻找一个X∈R 使得它满足： (NCP) F()≥0， 0，xTF(x)=0 这里F()是R 一R 上的连续函数． 解非线性互补问题的实质仍是求解全局最优化的全局最优点，这一问题显然可化 为如下有约束的全局优化问题： min{xTF()IF() 0， 0}． 定义 设F为Rn—Rn上的连续函数，S为R 中正的子空间，若满足 ，()=XTF(x)=0 其中 ，F(x)∈S，称，()为F()的互补函数． 我们考虑有约束的全局最优化问题的模型为： (P) c ，() 其中SCRn是一个紧的约束集，f：R 一R 上的连续函数，c为，()在S中的总 极值． 对这一问题最常用的处理方法是利用罚函数将有约束问题转化为无约束问题，再 结合求无约束总极值问题的方法解决．而我们在文 [1】中以另一种途径，通过构造一 个简单的新函数9()，结合积分水平集方法提出了一种新的概念性算法并证明了其收 敛性． 本文 2003年 4月4日收到． 维普资讯 应用数学与计算数学学报 18卷 设DDS是R“中的有界闭箱，算法如下： 算法A 步1，取co吨n，()，￡0，k：=0； 酬 = 一 HCk-~ 嘶 )c) 作函数g (X)满足： ㈤= ∈∈ 计算均值： c gc() 这里 (D)指集合D上的勒贝格测度． 步3，计算均方差： ／d(()一c) 若Vk+1 ￡，令k：=k+1，转步2；否则转步4； 步4，令C=Ok+1且H =Hc 则H 为f(x)在S中的近似总极值点集，C 为近似总极值；算法终止． 注意，其中co可以这样取：任取一xoES，然后再取cof(xo)即可． 类似于文 2【】，我们利用数论中的一致分布佳点集代替Monte-Caclo投点构造相应 的实现算法．首先介绍数论中的有关概念． 数论中的相关概念 令t为一正整数，及 (s)=(i(s)，·一，n￡(s))， 1 8 t 表示 中的一个点集，此处0 (s) 1，(1 i 几)， 为几维单位闭箱． 对于任意7=(71，72，…，)E ，设 (71，，…，)表示 (s)(1 8 t)中适 合不等式 0 xi(k)71，…，0 n￡() 的个数．若 l ～… 则称点集有偏差 ()．