投影算子的一种简单算法.docxVIP

投影算子的一种简单算法何松年，赵子祎（中国民航大学理学院，天津 300300）摘 要 ： 提 出 了 投 影 算 子 的 一 种 Halpern 型 的 松 弛 算 法 ， 由 于 这 种 算 法 把 计 算 关 于 一 个 凸 函 数 水 平 集 的 投 影 转 化为 计 算 关 于 一 列 包 含 水 平 集 的 半空 间 的 投 影 ，因 而 算 法 容易 实 现 ，并 且 证 明 了 算 法 的 强 收 敛 性 。关 键 词 ： 投 影 ；半空 间 ； 强 收 敛 ；Hilbert 空 间中 图 分 类 号 ： O177文 献 标 志 码 ： A文 章 编 号 ： 1674－5590（2014）04－0052－03Simple algorithm for projection operatorHE Song-nian， ZHAO Zi-yi（College of Science， CAUC， Tianjin 300300， China）Abstract：A relaxed Halpern’s projection algorithm is proposed. Since this algorithm computes the projection onto levelset of a convex function by computing the projection onto a series of half-spaces containing a level set， it is easy to be implemented. Strong convergence of this algorithm is proved.Key words：projection； half-space； strong convergence； Hilbert space设 H 是一实 Hilbert 空间，其内积与范数分别表示为·，·和‖·‖，C奂H 非空闭凸。从 H 到其非空闭 凸子集 C 上的投影 PC 定义为：对于任意 x∈H，必有 C 中唯一一点，记作 PC x，满足分别表示为·，·和|·| ||，C奂H 非空闭凸，T：C→C 非扩张并且满足 F（T）≠覫，又设1）λn→0（n→∞）；∞2）Σλn = ∞；∞PC x：= arg min ‖x - y‖2y∈C且 PC 有如下特征：对 x∈H，有x∈Hn=1λn3）Σ│λn+1n- λ │ +∞，lim= 1n→∞ λx- PC x，y - PC x≤0（1）坌y∈Cn=1则序列（xn）强收敛到 PFix( T )u。n+1投影算子在很多算法中都会涉及到：如不动点的求解，凸优化问题的求解，变分不等式[1]以及分裂可行 性问题[2-5]的求解等。由于对一般闭凸子集 C，投影算子 PC 没有显式表达式（除非 C 是闭球或者半空间等简单 情形），所以投影算子的计算一般难以实现。假设 T：C→C 为非扩张映像，其不动点集 Fix（T） 非空，众所周知 Fix（T）闭凸的。计算 T 的不动点的Halpern 迭代格式为xn+1 = λn u +（1 - λn）Txn其中：u∈C 取定，x0∈C 任意取定，{λn}奂（0，1）。关于Halpern 迭代的收敛性，有如下基本结果：定理 1设 H 是一实 Hilbert 空间，其内积与范数假设 c：H→R 为凸函数，本文讨论关于凸函数 c的水平集投影的计算方法，其中水平集为C = {x∈H│c（x）≤0}本文提出投影算子的一种 Halpern 型松弛算法， 具体如下xn+1 = λn u +（1 - λn）PC xnn其 中：取 u ∈H，任意取定初值 x0 ∈H，Cn 劢C（n = 1，2，…）是一列半空间（Cn 的确切构造见第 2 节），由于 这种算法是把计算关于一个凸函数水平集的投影转 化为计算关于一列包含水平集的半空间的投影，所以 此算法容易实现。将在第 2 节证明此算法产生的序列基 金 项 目 ： 中央高校基本科研业务费专项（3122013SY30）收 稿 日 期 ： 2013-07-12； 修 回 日 期 ： 2013-09-02作 者 简 介 ： 何松年（1963—），男，山西太原人，教授，博士，研究方向为非线性分析理论、算法及其应用.- 53 -第 32 卷 第 4 期何松年，赵子祎：投影算子的一种简单算法在一定条件下强收敛于 PC u。R 为一凸函数，总是假设 c 在 H 上是次可微的，并且鄣c 是有界算子（也就是在有界集上是有界的）。值得注 意的是定义在有限维 Hilbert 空间中的每一个凸函数都是次可微的，并且其次可微算子是有界的[8]。假设算法第 n 步迭代 xn 已经得到