对互补问题不可行内点算法探讨.docVIP

对互补问题不可行内点算法探讨 　　摘要：对一类非单调线性(P-矩阵)互补问题窄邻域和宽邻域的几种不可行内点算法做数值计算，验证其收敛性及实际表现，本文就对互补问题的不可行内点算法起进行简单的探讨。 　　关键词：互补；非线性；算法，不可行内点 　　【中图分类号】G642 　　一、互补问题简介 　　非线性互补问题，记为NCP(F),是指：求（x,s)∈Rn*Rn,使得s=F(X),(X,S)0≥0,xTs=0,其中F（x）=（f1(x),f2(x)，……fn(x)）T:∈Rn―Rn是连续映射。 　　记Ω=｛（x,s）|s=F(x),(x,s)≥0},Ω+={(x,s)|(x,s)∈Ω,(x,s)0},本文主要是针对互补问题的不可行内点算法做数值实验，而不可行内点算法是相对可行内点算法而言的.简言之，如果在迭代过程中能保证sk=F(xk),k=0,1……，则称为可行内点法；否则称为不可行内点算法.在初始点的选取方面不可行内点算法有明显的优势，因为内点算法对于初始点如何选取至今没有一个完整的方法，尤其是大规模问题，更是难上加难，就实际效果而言，不可行内点算法要比可行内点算法表现好现.。 　　（i）当F（x）为一致P函数时，分析结果表明，该方法是一多项式时间算法,（ii）当F(x)=Mx+q(M为n*n矩阵，q为n维列向量)时，NCP(F)退化为线性互补问题，记为LCP(q,M)。众所周知，当Ω+非空时，单调线性互补问题可在多项式时间内求解，而且人们已经设计出了多种求解单调互补问题的有效内点算法,对于非单调互补问题内点算法可以说是刚刚开始.本文假设当M是P矩阵时，不难验证它是一致P-函数的特殊情况，一种不可行内点算法是特殊情况，其中这种算法应用在窄邻域中。当为一致P函数时，互补问题宽邻域不可行内点算法。本文提供的算法进行编程研究，用数值实验的方法验证不可行内点算法的收敛性。 　　二、互补问题符号约定及基本概念 　　本文约定当英文小写字母例如：x,dx表示向量时，对应的大写字母表示对角矩阵X=diag(x), Dx=diag(dx). e=(1,...,1)T, E为n 阶的单位矩阵， 分别表示2-范数和∞-范数。 x=(x1,x2,...xn)T Rn, 记向量x+=( )T, x-=（ )T其中记 用 表示矩阵M+MT的最小特征值并假定 （当 时，LCP(q,M)问题成为单调线性互补问题）。 　　定义1.1称n阶方阵 为 　　（i） 矩阵，若M的每一个主子式均非负； 　　（ii） 矩阵，若M的每一个主子式均为正。 　　当M为 矩阵时，还可以使用另一个定义方式，即所谓 矩阵是满足下面的条件的矩阵： 记 则有 ， ， 称为 矩阵M的正惯指数，显然 矩阵的概念是非对称矩阵正定性推广。 　　定义1.2设 F： ,称F为 　　（i） 函数，若对任意 , ,存在 使得 且 ； 　　（ii）一致 函数，若存在 ,使 ， 　　三、互补问题的不可行内点算法描述 　　本文首先采用2-范数对一类非单调线性互补问题的不可行内点算法进行描述，对于问题LCP(q,M),假定M是P-矩阵，记s=Mx+q, 　　r=r(x,s)=s-Mx-q,(x,s)∈Rn*Rn,C={(x,s)∈Ω|Xs=μe,0μ≤μ0, 　　 　　 　　 　　其中 而β∈(0,1)及0,ρ，μ0是适当选取的常数；称r为问题LCP(q,M)在（x,s）处的的不可行残向量；称C为问题LCP(q,M)的中心路径；称N2（β）， （β）为中心路径的各个领域，用N（β）泛指N2（β）， 都适用。 　　取x0=e,s0=e,r0=s0-Mx0-q,μ0=[(x0)Ts0]/n=1,p=‖r0‖/ ，则容易验证（x0,s0）∈N(β)。 　　记 　　求解线性互补问题的不可行内点算法的基本思想是通过产生N（β）中的迭代点列（xk,sk）逼近问题的解，即使其满足 。 　　具体实施时，对问题LCP(q,M)的近似解（x,s)∈N2β，在当前迭代首先求解方程组 : 　　ds-Mdx=-(1-μ)r,Xds+Sdx=rμe-Xs，其中η，r∈(0,1)是适当选取的常数，以下取: 　　 ,然后记 x(θ)=x+θdx,s(θ)=s+θdx,r(θ)=s(θ)-Mx(θ)-q,并求θ1使其满足θ1=max{θ-∈(0,1)}|(x(θ),s(θ))∈N(β),θ∈(0,θ-)}。最后求Φ*使其满足Φ（x(θ),s(θ*)）= 　　并取x+=x(θ*),s+=S(θ*),r+=r(θ+)作为问题LCP(q,M)的新近似解及不可行残量。 　　综合上述可得如下的求解问题LCP(q,M)的不可行内点算法的基本步骤： 　　（i）给定数,ε0，令k=0； 　　（