求解双障碍问题的一个内点迭代算法.pdfVIP

第 32卷 第 6期 武汉理工大学学报(鸯鍪差) Vo1．32 No．6 2O08年 12月 JournalofW uhanUniVersityofTechnology Dec．20O8 (Transp0rtationScienceEngineering) 求解双障碍问题的一个内点迭代算法* 吴绪权 (武汉理工大学理学院 武汉 430063) 摘要：介绍一种解决双障碍问题的迭代算法，该算法称之为内点迭代算法．首先将原问题转化为等 价非光滑方程组 ，然后构造新的光滑函数来逼近非光滑方程组．文中证明了惩罚参数可能需要修 正有限次情形下算法的全局收敛性和超线性收敛性．经数值实验表明，该算法是有效的． 关键词：双障碍问题；内点迭代算法；收敛性定理 中图法分类号 ：O224 考虑尺“中的双障碍问题 ：求z∈ ，使 强F导数 ，并且它是处处方向可微的，在向量 处 (／(‘)， — z)≥ 0，V ∈X (1) 沿d的导数由下式给出． 式 中：厂： 一R，I连续可微 ； 一 {z∈尺 l口≤ ≤ g(，d)一 ∑ (z一口)+ ∑ (一口)× f∈， () ∈ ( ，一1，…，}；f， 为常数． 。 式(1)是变分不等式的重要的类型之一，主要 min(， ()d)+ ∑ ) )d+ i∈I z’ 来源于经济领域 中的实际问题 (如平衡问题等)和 ∑ (一z)min(一 ，一V ()d)+ 物理、力学及工程领域 中的有关问题 (一般为无穷 ∈ () 维问题的离散化[1]．由于其应用的广泛性 ，对双障 ∑ (z一 ) 。 f∈ ( 碍问题 的数值解法的研究近年来很受重视．构造 式 中： ()一 {： ()zf一 日f}；J ()= {： 数值解法的基本途径之一是将方程组的解法直接 ()一zf一口f)；，，()= {：z一6 ()zf一 推广到对双障碍 问题的求解 [2]，如线性化牛顿法 口)； ()一 {： ()一z一6f}；，6()一 {：，l() 和拟牛顿法等．本文试图在此基础上提出与式(1) z一 }．令 等价的非光滑优化问题的内点迭代法．设文中使 ()一 L()U L () 用的内积和范数都是欧氏空间的． ()= L()U () 1 算 法 则若 ∈ ，有 g(，)≤ ∑ (五一日)+ ∑ ()× f∈J口() ∈，，(’ 若记 口一(口1，…，以)，西一(6，…，6)贝Ⅱ式