跳到主要内容

用全局收敛算法求解具有无界非凸约束集的非线性规划问题

摘要

本文提出了一组无界条件,在此条件下,我们可以通过组合同伦内点算法求解一类无界非凸集合中的非线性规划问题。我们还获得了组合同伦内点算法的全局收敛性,并详细分析了该算法的实现。

1引言

考虑约束最小化问题:

最小值 x个 n个 (f) ( x个 ) 科学技术 u个 ( x个 ) 0 , = 1 , , , 不锈钢 v(v) j个 ( x个 ) = 0 , j个 = 1 , , ,
(1)

哪里(f): n个 1 , u个 : n个 1 ,=1,, v(v) j个 : n个 1 ,j个=1,,是连续可微的三倍。

我们称之为要点 x个 是(1)的Karush-Kuhn-Tucker(K-K-T)点,并且 , z(z) 是对应的拉格朗日乘子向量,如果( x个 , , z(z) )满足

(f) ( x个 ) + u个 ( x个 ) + v(v) ( x个 ) z(z) = 0 , v(v) ( x个 ) = 0 , Y(Y) u个 ( x个 ) = 0 , u个 ( x个 ) 0 , 0 ,
(2)

哪里 ,z(z) ,(f)(x个)= ( (f) ( x个 ) / x个 ) T型 n个 ,u个(x个)=( u个 1 (x个),, u个 (x个)) n个 × ,v(v)(x个)=( v(v) 1 (x个),, v(v) (x个)) n个 × ,u个(x个)= ( u个 1 ( x个 ) , , u个 ( x个 ) ) T型 ,v(v)(x个)= ( v(v) 1 ( x个 ) , , v(v) ( x个 ) ) T型 Y(Y)=诊断() × .

自从Kellogg等。(请参见[1])和Smale(参见[2])提出了著名的同伦方法,该方法已成为处理各种非线性问题的有力工具,例如映射的零点或不动点(参见[5],等。以及其中的参考)。然而,直到1988年Megiddo(参见[6])和小岛等。(请参见[7])发现线性规划的Karmarkar内点法是一种路径允许方法。从那时起,数学规划的中心路径允许方法成为一个活跃的研究课题。此外,它最近被扩展到凸非线性规划问题(参见[811],等。). 但它们的所有收敛结果都是在对数障碍函数严格凸且解集非空且有界的假设下得到的。

最近,在中提出了一种求解非线性规划问题的组合单纯形内点方法(为了方便起见,称为CHIP方法)[12,13](详细摘要于[14]). 在[13]与中心路径允许方法相比,作者消除了对数障碍函数的凸性以及解集的有界性和非空性。这表明,CHIP方法可以解决内部路径允许方法无法解决的问题。在[15]采用分段技术,在常用条件下,给出了CHIP方法的多项式,这表明CHIP方法也具有很好的效率。在[16],我们介绍了 C类 2 功能 α (x个) n个 ,=1,, β j个 (x个) n个 ,j个=1,,将结果扩展到[12,17]到更一般的非凸集。然而,文献中没有关于[16]扩展到无界非凸集。在本文中,我们试图完成这项工作。为此,通过使用一些不等式技术和在[18,19],我们发展了一组新的无界条件。在这些条件下,我们得到了CHIP方法的全局收敛结果,从而扩展了[16]到无界非凸集。

2主要成果

在本节中 表示为 + + + 分别是。我们还表示为B类(x个)={{1,,}: u个 (x个)=0}活动集位于x个此外,让X(X)={x个 n个 : u个 (x个)0,=1,,, v(v) j个 (x个)=0,j个=1,,}作为可行集,让 X(X) 0 ={x个 n个 : u个 (x个)<0,=1,,, v(v) j个 (x个)=0,j个=1,,}成为严格可行的集合X(X)=X(X) X(X) 0 是的边界集X(X).

在[16],为了求解(2),构造了以下同伦:

H(H) ( w个 , w个 ( 0 ) , λ ) = ( ( 1 λ ) [ (f) ( x个 ) + u个 ( x个 ) + λ ξ ( x个 ) 2 ] + [ v(v) ( x个 ) + λ ( η ( x个 ) v(v) ( x个 ) ) ] z(z) + λ ( x个 x个 ( 0 ) ) v(v) ( x个 ) Y(Y) u个 ( x个 ) λ Y(Y) ( 0 ) u个 ( x个 ( 0 ) ) ) = 0 ,
(3)

哪里w个=(x个,,z(z)) n个 + + , w个 ( 0 ) =( x个 ( 0 ) , ( 0 ) ,0) X(X) 0 × + + × 、和 Y(Y) ( 0 ) =诊断( ( 0 ) ) × .

在本节中,我们使用无限解的概念[18,19]从而给出了一组新的无界条件。如果存在序列,则称非线性规划问题在无穷远处有解{ x个 ( k个 ) }令人满意的{ x个 ( k个 ) }X(X), x个 ( k个 ) 作为k个,对于任何给定的ζX(X),存在 ( k个 ) + z(z) ( k个 ) 这样的话

k个 ( ζ x个 ( k个 ) ) T型 ( (f) ( x个 ( k个 ) ) + u个 ( x个 ( k个 ) ) ( k个 ) + v(v) ( x个 ( k个 ) ) z(z) ( k个 ) ) 0.

然后我们假设存在光滑映射 ξ (x个) 1 ,=1,, η j个 (x个) 1 ,j个=1,,这样:

(A)1) X(X) 0 非空;非线性规划问题没有无限解。

(A)2)x个X(X),如果

B类 ( x个 ) ( ( x个 ) + α ξ ( x个 ) ) +v(v)(x个)z(z)+η(x个)β=0, 0, α 0,z(z) ,β ,

然后 =0, α =0z(z)=β=0,B类(x个).

(A))x个X(X),我们有

{ x个 + B类 ( x个 ) ξ ( x个 ) + z(z) η ( x个 ) : 0 , B类 ( x个 ) , z(z) } X(X)={x个},

哪里η(x个)=( η 1 (x个),, η (x个)).

(A)4)x个X(X),v(v)(x个)列秩满v(v) ( x个 ) T型 η(x个)是非奇异的。

然而,仅将无限解技术应用于项目(f)(x个),u个(x个)、和v(v)(x个),我们无法将结果扩展到[16]到无界情况,因为项目的存在ξ(x个)η(x个)对于这些项目ξ(x个)η(x个),我们需要使用其他技术。在本文中,我们仍然需要以下假设:

(A)5)ζX(X),

u个 ( ζ ) T型 u个 ( x个 ) T型 ( ζ x个 ) T型 ξ (x个) ,=1,,,v(v) ( ζ ) T型 v(v) ( x个 ) T型 ( ζ x个 ) T型 η(x个).

对于给定的 w个 ( 0 ) ,重写H(H)(w个, w个 ( 0 ) ,λ)作为 H(H) w个 ( 0 ) (w个,λ)。的零点集 H(H) w个 ( 0 )

H(H) w个 ( 0 ) 1 (0)= { ( w个 , λ ) X(X) × + × × ( 0 , 1 ] : H(H) w个 ( 0 ) ( w个 , λ ) = 0 } .

逆像定理告诉我们,如果0是映射的正则值 H(H) w个 ( 0 ) ,那么 H(H) w个 ( 0 ) 1 (0)由一些平滑曲线组成。的规律性 H(H) w个 ( 0 ) 可以通过以下引理得到。

引理2.1(参数化萨尔德定理)

V(V) n个 ,U型 成为开放集, Φ:V(V)×U型 ( k个 ) C类 第页 地图,哪里 第页>最大值{0,k个}.如果 0 ( k个 ) 是的常规值Φ,那么几乎所有人 V(V), 0是的常规值 Φ Φ(,).

引理2.2 H(H) 定义如下(3).此外,let假设(A)1)-(A)5)持有, u个 (x个),=1,, v(v) j个 (x个),j个=1,, C类 功能,然后让 ξ (x个),=1,, η j个 (x个),j个=1,, C类 2 功能.然后,几乎所有人 w个 ( 0 ) X(X) 0 × + + × , 0是的常规值 H(H) w个 ( 0 ) , H(H) w个 ( 0 ) 1 (0) 由一些平滑曲线组成,其中有一条平滑的曲线 Γ w个 ( 0 ) 从开始 ( w个 ( 0 ) ,1).

引理2.3 H(H) 定义如下(3).此外,let假设(A)1)-(A)5)持有, u个 (x个),=1,, v(v) j个 (x个),j个=1,, C类 功能,然后让 ξ (x个),=1,, η j个 (x个),j个=1,, C类 2 功能.然后,几乎所有 w个 ( 0 ) X(X) 0 × + + × ,光滑曲线的投影 Γ w个 ( 0 ) 上的 x个-平面有界.

证明如果不是,则存在一系列点 { ( x个 ( k个 ) , ( k个 ) , z(z) ( k个 ) , λ k个 ) } k个 = 1 这样的话 x个 ( k个 ) 作为k个.

很容易证明以下不等式成立:

x个 ζ 2 x个 ( 0 ) ζ 2 2 ( x个 ζ ) T型 ( x个 x个 ( 0 ) ) .
(4)

通过同伦方程(),我们有

H(H) ( w个 ( k个 ) , w个 ( 0 ) , λ k个 ) = ( ( 1 λ k个 ) [ (f) ( x个 ( k个 ) ) + u个 ( x个 ( k个 ) ) ( k个 ) + λ k个 ξ ( x个 ( k个 ) ) ( ( k个 ) ) 2 ] + λ k个 ( x个 ( k个 ) x个 ( 0 ) ) + [ v(v) ( x个 ( k个 ) ) + λ k个 ( η ( x个 ( k个 ) ) v(v) ( x个 ( k个 ) ) ) ] z(z) ( k个 ) v(v) ( x个 ( k个 ) ) Y(Y) ( k个 ) u个 ( x个 ( k个 ) ) λ k个 Y(Y) ( 0 ) u个 ( x个 ( 0 ) ) ) = 0 .
(5)

将(5)中的第一个方程乘以 ( x个 ( k个 ) ζ ) T型 ,我们得到

( 1 λ k个 ) ( x个 ( k个 ) ζ ) T型 (f) ( x个 ( k个 ) ) + ( 1 λ k个 ) ( x个 ( k个 ) ζ ) T型 u个 ( x个 ( k个 ) ) ( k个 ) + ( 1 λ k个 ) λ k个 ( x个 ( k个 ) ζ ) T型 ξ ( x个 ( k个 ) ) ( ( k个 ) ) 2 + ( x个 ( k个 ) ζ ) T型 v(v) ( x个 ( k个 ) ) z(z) ( k个 ) + ( x个 ( k个 ) ζ ) T型 λ k个 ( η ( x个 ( k个 ) ) v(v) ( x个 ( k个 ) ) ) z(z) ( k个 ) + λ k个 ( x个 ( k个 ) ζ ) T型 ( x个 ( k个 ) x个 ( 0 ) ) = 0 ,
(6)

,

λ k个 ( x个 ( k个 ) ζ ) T型 ( x个 ( k个 ) x个 ( 0 ) ) = ( 1 λ k个 ) ( x个 ( k个 ) ζ ) T型 (f) ( x个 ( k个 ) ) ( 1 λ k个 ) ( x个 ( k个 ) ζ ) T型 u个 ( x个 ( k个 ) ) ( k个 ) ( 1 λ k个 ) λ k个 ( x个 ( k个 ) ζ ) T型 ξ ( x个 ( k个 ) ) ( ( k个 ) ) 2 ( x个 ( k个 ) ζ ) T型 [ v(v) ( x个 ( k个 ) ) + λ k个 ( η ( x个 ( k个 ) ) v(v) ( x个 ( k个 ) ) ) ] z(z) ( k个 ) = ( 1 λ k个 ) ( x个 ( k个 ) ζ ) T型 (f) ( x个 ( k个 ) ) ( 1 λ k个 ) ( x个 ( k个 ) ζ ) T型 u个 ( x个 ( k个 ) ) ( k个 ) ( 1 λ k个 ) λ k个 ( x个 ( k个 ) ζ ) T型 ξ ( x个 ( k个 ) ) ( ( k个 ) ) 2 ( 1 λ k个 ) ( x个 ( k个 ) ζ ) T型 v(v) ( x个 ( k个 ) ) z(z) ( k个 ) λ k个 ( x个 ( k个 ) ζ ) T型 η ( x个 ( k个 ) ) z(z) ( k个 ) .
(7)

所以

λ k个 ( x个 ( k个 ) ζ 2 x个 ( 0 ) ζ 2 ) 2 λ k个 ( x个 ( k个 ) ζ ) T型 ( x个 ( k个 ) x个 ( 0 ) ) = 2 ( 1 λ k个 ) ( x个 ( k个 ) ζ ) T型 (f) ( x个 ( k个 ) ) 2 ( 1 λ k个 ) ( x个 ( k个 ) ζ ) T型 u个 ( x个 ( k个 ) ) ( k个 ) 2 ( 1 λ k个 ) λ k个 ( x个 ( k个 ) ζ ) T型 ξ ( x个 ( k个 ) ) ( ( k个 ) ) 2 2 ( 1 λ k个 ) ( x个 ( k个 ) ζ ) T型 v(v) ( x个 ( k个 ) ) z(z) ( k个 ) 2 λ k个 ( x个 ( k个 ) ζ ) T型 η ( x个 ( k个 ) ) z(z) ( k个 ) = 2 ( 1 λ k个 ) ( x个 ( k个 ) ζ ) T型 ( (f) ( x个 ( k个 ) ) + u个 ( x个 ( k个 ) ) ( k个 ) + v(v) ( x个 ( k个 ) ) z(z) ( k个 ) ) + 2 ( 1 λ k个 ) λ k个 ( ζ x个 ( k个 ) ) T型 ξ ( x个 ( k个 ) ) ( ( k个 ) ) 2 + 2 λ k个 ( ζ x个 ( k个 ) ) T型 η ( x个 ( k个 ) ) z(z) ( k个 ) 2 ( 1 λ k个 ) ( x个 ( k个 ) ζ ) T型 ( (f) ( x个 ( k个 ) ) + u个 ( x个 ( k个 ) ) ( k个 ) + v(v) ( x个 ( k个 ) ) z(z) ( k个 ) ) + 2 ( 1 λ k个 ) λ k个 ( u个 ( ζ ) T型 ( k个 ) u个 ( x个 ( k个 ) ) T型 ( k个 ) ) + 2 λ k个 ( v(v) ( ζ ) T型 z(z) ( k个 ) v(v) ( x个 ( k个 ) ) T型 z(z) ( k个 ) ) .
(8)

u个 ( ζ ) T型 0 ( k个 ) 0,那么u个 ( ζ ) T型 ( k个 ) 0此外,根据(5)中的第三个方程,我们得出u个 ( x个 ( k个 ) ) T型 ( k个 ) = λ k个 u个 ( x个 ( 0 ) ) T型 ( 0 ) 因此,以下不等式成立:

λ k个 ( x个 ( k个 ) ζ 2 x个 ( 0 ) ζ 2 ) 2 ( 1 λ k个 ) ( x个 ( k个 ) ζ ) T型 ( (f) ( x个 ( k个 ) ) + u个 ( x个 ( k个 ) ) ( k个 ) + v(v) ( x个 ( k个 ) ) z(z) ( k个 ) ) 2 ( 1 λ k个 ) λ k个 u个 ( x个 ( k个 ) ) T型 ( k个 ) = 2 ( 1 λ k个 ) ( x个 ( k个 ) ζ ) T型 ( (f) ( x个 ( k个 ) ) + u个 ( x个 ( k个 ) ) ( k个 ) + v(v) ( x个 ( k个 ) ) z(z) ( k个 ) ) 2 ( 1 λ k个 ) λ k个 2 u个 ( x个 ( 0 ) ) T型 ( 0 ) .
(9)

那么到了(9),我们已经

( ζ x个 ( k个 ) ) T型 ( (f) ( x个 ( k个 ) ) + u个 ( x个 ( k个 ) ) ( k个 ) + v(v) ( x个 ( k个 ) ) z(z) ( k个 ) ) λ k个 2 ( 1 λ k个 ) ( x个 ( k个 ) ζ 2 x个 ( 0 ) ζ 2 ) + λ k个 2 u个 ( x个 ( 0 ) ) T型 ( 0 ) .
(10)

什么时候? x个 ( k个 ) ,到(10),我们已经

k个 ( ζ x个 ( k个 ) ) T型 ( (f) ( x个 ( k个 ) ) + u个 ( x个 ( k个 ) ) ( k个 ) + v(v) ( x个 ( k个 ) ) z(z) ( k个 ) ) k个 λ k个 2 ( 1 λ k个 ) ( x个 ( k个 ) ζ 2 x个 ( 0 ) ζ 2 ) + λ k个 2 u个 ( x个 ( 0 ) ) T型 ( 0 ) 0 ,
(11)

这与假设(A)相矛盾1). □

定理2.1 H(H) 定义如下(3), (f)(x个), u个 (x个),=1,, v(v) j个 (x个),j个=1,, 是三次连续可微函数,let假设(A)1)-(A)5)持有,然后让 ξ (x个),=1,, η j个 (x个),j个=1,, 是两倍连续可微函数.然后几乎所有人 w个 ( 0 ) X(X) 0 × + + × ,存在一个 C类 1 曲线 (w个(),λ()) 尺寸的1这样的话

H(H) ( w个 ( ) , w个 ( 0 ) , λ ( ) ) =0, ( w个 ( 0 ) , λ ( 0 ) ) = ( w个 ( 0 ) , 1 ) .
(12)

什么时候? λ()0,w个() 趋于某一点 w个 =( x个 , , z(z) ).特别地,组件 x个 属于 w个 是K-K(K)-问题的T点(1).

证明根据引理2.2,必须有一个 C类 1 曲线(w个(),λ())尺寸1(表示为 Γ w个 ( 0 ) )这样的话

H(H) ( w个 ( ) , w个 ( 0 ) , λ ( ) ) =0, ( w个 ( 0 ) , λ ( 0 ) ) = ( w个 ( 0 ) , 1 ) .

根据一维光滑流形的分类定理, Γ w个 ( 0 ) 既不同于单位圆,也不同于单位区间。对于任何 w个 ( 0 ) X(X) 0 × + + × , H(H) w个 ( 0 ) ( w个 ( 0 ) ,1)/w个是非奇异的,所以 Γ w个 ( 0 ) 不能与单位圆相异。那就是, Γ w个 ( 0 ) 与单位区间不同。

( w个 , λ )是的一个极限点 Γ w个 ( 0 ) .如果( w个 , λ ) X(X) 0 × + + × ×(0,1),因为0是的常规值 H(H) w个 ( 0 ) ,( w个 , λ ) H(H) w个 ( 0 ) 1 (0),以及的雅可比矩阵H(H)( w个 , λ )是满行秩的,那么根据隐函数定理, Γ w个 ( 0 ) 可以在以下位置扩展( w个 , λ )这一结果与以下事实相矛盾:( w个 , λ )是的极限点 Γ w个 ( 0 ) .

( w个 , λ )=( x个 , , z(z) , λ )因此,( w个 , λ )(X(X)× + × ×(0,1])可能出现以下三种情况:

  1. (a)

    ( w个 , λ )=( x个 , , z(z) , λ )X(X)× + × ×{0}.

  2. (b)

    ( w个 , λ )=( x个 , , z(z) , λ ) X(X) 0 × + + × ×{1}.

  3. (c)

    ( w个 , λ )=( x个 , , z(z) , λ )(X(X)× + × )×(0,1].

由于方程式H(H)(w个, w个 ( 0 ) ,1)=0有独特的解决方案( w个 ( 0 ) ,1)在里面 X(X) 0 × + + × ×{1},情况(b)是不可能的。

通过同伦方程(),我们有

( 1 λ k个 ) [ (f) ( x个 ( k个 ) ) + u个 ( x个 ( k个 ) ) ( k个 ) + λ k个 ξ ( x个 ( k个 ) ) ( ( k个 ) ) 2 ] + λ k个 ( x个 ( k个 ) x个 ( 0 ) ) + [ v(v) ( x个 ( k个 ) ) + λ k个 ( η ( x个 ( k个 ) ) v(v) ( x个 ( k个 ) ) ) ] z(z) ( k个 ) = 0 ,
(13)
v(v) ( x个 ( k个 ) ) =0,
(14)
Y(Y) ( k个 ) u个 ( x个 ( k个 ) ) λ k个 Y(Y) ( 0 ) u个 ( x个 ( 0 ) ) =0.
(15)

(x个)= { = 1 , , : k个 ( k个 ) = + } ,J型(x个)= { j个 = 1 , , : k个 z(z) j个 ( k个 ) = + } .

如果J型(x个),那么

( 1 λ k个 ) [ (f) ( x个 ( k个 ) ) + ( x个 ) ( u个 ( x个 ( k个 ) ) ( k个 ) + λ k个 ξ ( x个 ( k个 ) ) ( ( k个 ) ) 2 ) ] + λ k个 ( x个 ( k个 ) x个 ( 0 ) ) + [ ( v(v) ( x个 ( k个 ) ) + λ k个 ( η ( x个 ( k个 ) ) v(v) ( x个 ( k个 ) ) ) ) z(z) ( k个 ) + ( 1 λ k个 ) ( x个 ) ( u个 ( x个 ( k个 ) ) ( k个 ) + λ k个 ξ ( x个 ( k个 ) ) ( ( k个 ) ) 2 ) ] = 0 .
(16)

因为X(X)(0,1]是有界的,通过假设(A2),(16)左侧的第三部分趋于无穷大,如下所示k个,但其他两部分是有界的,这是不可能的。因此,我们得出以下结论:平滑曲线的投影 Γ w个 ( 0 ) 上的z(z)-平面也是有界的。

在情况(c)中,首先,我们证明 + .如果 + ,那么就存在 0 {1,,}和一系列的点{( w个 ( k个 ) , λ k个 )} Γ w个 ( 0 ) 这样的话 0 ( k个 ) 0 =0作为k个+从(15)开始,我们有

0 ( k个 ) u个 0 ( x个 ( k个 ) ) = λ k个 0 ( 0 ) u个 0 ( x个 ( 0 ) ) .
(17)

因为X(X)(0,1]有界,当k个+,(17)的左侧趋于0。同时,(17)的右侧倾向于 λ 0 ( 0 ) u个 0 ( x个 ( 0 ) ),严格小于0。这一事实导致了矛盾。

那么我们只需要证明案例(c)的其余部分是不可能的。如果不是,则存在一系列点{( w个 ( k个 ) , λ k个 )} Γ w个 ( 0 ) 这样的话 u个 ( x个 ( k个 ) )0对一些人来说{1,,}作为k个+从(15)中,我们得到 ( k个 ) +.因为X(X)(0,1]有界,则存在一个点的子序列(也表示为{( w个 ( k个 ) , λ k个 )})这样的话 x个 ( k个 ) x个 , ( k个 ) +, z(z) ( k个 ) z(z) 、和 λ k个 λ 作为k个+.

什么时候? λ >0,从(15)开始,活动索引集为B类( x个 )=( x个 ).何时 λ =0,索引集( x个 )B类( x个 ).

(1) 何时 λ =1根据(16) ( k个 ) 每个都有界( x个 ),假设(A1)-(A)2),我们得出结论 k个 + (1 λ k个 ) ( k个 ) =0 k个 + (1 λ k个 ) ( ( k个 ) ) 2 = 因此,当k个+,(16)成为

η ( x个 ) z(z) + ( x个 ) ξ ( x个 ( k个 ) ) + x个 = x个 ( 0 ) ,
(18)

这与假设(A)相矛盾).

(2) 何时0< λ <1,将(16)重写为

( x个 ) ( u个 ( x个 ( k个 ) ) ( k个 ) + λ k个 ξ ( x个 ( k个 ) ) ( ( k个 ) ) 2 ) + λ k个 1 λ k个 ( x个 ( k个 ) x个 ( 0 ) ) = (f) ( x个 ( k个 ) ) ( x个 ) ( u个 ( x个 ( k个 ) ) ( k个 ) + λ k个 ξ ( x个 ( k个 ) ) ( ( k个 ) ) 2 ) 1 1 λ k个 ( v(v) ( x个 ( k个 ) ) + λ k个 ( η ( x个 ( k个 ) ) v(v) ( x个 ( k个 ) ) ) ) z(z) ( k个 ) .
(19)

( k个 ) 是一个有分量的向量 ( ( k个 ) ) 2 ,( x个 ).然后设置

ρ ( k个 ) = ( ( k个 ) ) 2 ( k个 ) , ( x个 ) .
(20)

请注意0 ρ ( k个 ) 1; 然后存在一个子序列{ ρ ( k个 ) },仍表示为{ ρ ( k个 ) },因此 ρ ( k个 ) ρ 对于每个( x个 )作为k个+此外,包含分量的向量 ρ ,( x个 )表示为 ρ ; 因此, ρ =1将(19)的两边除以 ( k个 ) 并让k个+,我们有

( x个 ) λ ρ ξ ( x个 ) =0,

这与假设(A)相矛盾2).

(3) 何时 λ =0,因为非空索引集( x个 )B类( x个 ),(3)的证明与(2)的证明相似。

从上述讨论中,我们得出结论,(a)是唯一可能的情况。因此x个-的组件 w个 是(1)的K-K-T点。□

3算法分析

几乎所有人 w个 ( 0 ) =( x个 ( 0 ) , ( 0 ) ,0) X(X) 0 × + + × ,根据定理2.1,同伦生成 C类 1 曲线 Γ w个 ( 0 ) 通过微分(12)中的第一个方程,我们得到了以下定理。

定理3.1 同伦路径 Γ w个 ( 0 ) 由以下常微分方程的初值问题决定:

D类 H(H) w个 ( 0 ) ( w个 ( ) , λ ( ) ) ( w个 ˙ ( ) λ ˙ ( ) ) =0, ( w个 ( 0 ) , λ ( 0 ) ) = ( w个 ( 0 ) , 1 ) ,
(21)

哪里 是曲线的弧长 Γ w个 ( 0 ) .

关于如何数值跟踪同伦路径,已有许多预测-校正算法,请参阅[20],等。以供参考。因此,我们在本文中省略了它们。在算法的实现中,通常我们需要致力于在 Γ w个 ( 0 ) 保持行列式的符号 | D类 H(H) w个 ( 0 ) ( w个 , λ ) 第页 T型 | 不变量。在第一次迭代中,符号由以下引理确定。

引理3.1 如果 Γ w个 ( 0 ) 是平滑的,然后是正向 第页 ( 0 ) 在初始点 ( w个 ( 0 ) ,1) 满足 签名 | D类 H(H) w个 ( 0 ) ( w个 ( 0 ) , 1 ) 第页 ( 0 ) T型 | = ( 1 ) + + 1 .

证明

=(f) ( x个 ( 0 ) ) u个 ( x个 ( 0 ) ) +(12λ)ξ ( x个 ( 0 ) ) ( ( 0 ) ) 2 ,b条= Y(Y) ( 0 ) ( x个 ( 0 ) ) ,

然后

H(H) w个 ( 0 ) ( w个 ( 0 ) , 1 ) ( w个 , λ ) = ( 0 η ( x个 ( 0 ) ) v(v) ( x个 ( 0 ) ) T型 0 0 0 Y(Y) ( 0 ) u个 ( x个 ( 0 ) ) T型 诊断 ( u个 ( x个 ( 0 ) ) ) 0 b条 ) =( M(M) 1 , M(M) 2 ),

哪里 M(M) 1 ( n个 + + ) × ( n个 + + ) , M(M) 2 ( n个 + + ) × 1 .切线向量 第页 ( 0 ) 属于 Γ w个 ( 0 ) ( w个 ( 0 ) ,1)满足

( M(M) 1 , M(M) 2 ) 第页 ( 0 ) =( M(M) 1 , M(M) 2 ) ( 第页 1 ( 0 ) 第页 2 ( 0 ) ) =0,
(22)

哪里 第页 1 ( 0 ) n个 + + , 第页 2 ( 0 ) 1 .

它很容易得到 第页 1 ( 0 ) = M(M) 1 1 M(M) 2 第页 2 ( 0 ) ,那么

| D类 H(H) w个 ( 0 ) ( w个 ( 0 ) , 1 ) 第页 ( 0 ) T型 | = | M(M) 1 M(M) 2 第页 1 ( 0 ) T型 第页 2 ( 0 ) T型 | = | M(M) 1 M(M) 2 M(M) 2 T型 M(M) 1 T型 1 | 第页 2 ( 0 ) = | M(M) 1 | 第页 2 ( 0 ) ( 1 + M(M) 2 T型 M(M) 1 T型 M(M) 1 1 M(M) 2 ) = | 0 η ( x个 ( 0 ) ) v(v) ( x个 ( 0 ) ) T型 0 0 Y(Y) ( 0 ) u个 ( x个 ( 0 ) ) T型 诊断 ( u个 ( x个 ( 0 ) ) ) 0 | 第页 2 ( 0 ) ( 1 + M(M) 2 T型 M(M) 1 T型 M(M) 1 1 M(M) 2 ) = | 诊断 ( u个 ( x个 ( 0 ) ) ) | | η ( x个 ( 0 ) ) v(v) ( x个 ( 0 ) ) T型 0 | 第页 2 ( 0 ) ( 1 + M(M) 2 T型 M(M) 1 T型 M(M) 1 1 M(M) 2 ) = ( 1 ) | 诊断 ( u个 ( x个 ( 0 ) ) ) | | v(v) ( x个 ( 0 ) ) T型 η ( x个 ( 0 ) ) | 第页 2 ( 0 ) ( 1 + M(M) 2 T型 M(M) 1 T型 M(M) 1 1 M(M) 2 ) .

u个( x个 ( 0 ) )<0,(1+ M(M) 2 T型 M(M) 1 T型 M(M) 1 1 M(M) 2 )>0、和 第页 2 ( 0 ) <0,所以

| D类 H(H) w个 ( 0 ) ( w个 ( 0 ) , 1 ) 第页 ( 0 ) T型 |

( 1 ) + + 1 . □

以下伪代码描述了通用预测-校正方法的基本步骤。

算法3.1(欧拉-纽顿法)

步骤0。提供初步猜测( w个 ( 0 ) ,1),初始步长 小时 0 >0,和三个小正数 ϵ 1 >0, ϵ 2 >0、和 ϵ >0.设置k个=0.

步骤1。计算方向 θ ( k个 ) 预测步骤。

  1. (a)

    计算单位切线向量 第页 ( k个 ) .

  2. (b)

    确定方向 θ ( k个 ) 预测步长如下:

如果行列式的符号 | D类 H(H) w个 ( 0 ) ( w个 ( k个 ) , λ k个 ) 第页 ( k个 ) T型 | ( 1 ) + + 1 ,那么 θ ( k个 ) = 第页 ( k个 ) .

如果行列式的符号 | D类 H(H) w个 ( 0 ) ( w个 ( k个 ) , λ k个 ) 第页 ( k个 ) T型 | ( 1 ) + ,那么 θ ( k个 ) = 第页 ( k个 ) .

第2步。计算校正点( w个 ( k个 + 1 ) , λ k个 + 1 ).

( w个 ¯ ( k个 ) , λ ¯ k个 ) = ( w个 ( k个 ) , λ k个 ) + 小时 k个 θ ( k个 ) , ( w个 ( k个 + 1 ) , λ k个 + 1 ) = ( w个 ¯ ( k个 ) , λ ¯ k个 ) D类 H(H) w个 ( 0 ) ( w个 ¯ ( k个 ) , λ ¯ k个 ) + H(H) w个 ( 0 ) ( w个 ¯ ( k个 ) , λ ¯ k个 ) .

如果 H(H) w个 ( 0 ) ( w个 ( k个 + 1 ) , λ k个 + 1 ) ϵ 1 ,然后让 小时 k个 + 1 =最小值{ 小时 0 ,2 小时 k个 },然后转至步骤3。

如果 H(H) w个 ( 0 ) ( w个 ( k个 + 1 ) , λ k个 + 1 )( ϵ 1 , ϵ 2 ),然后让 小时 k个 + 1 = 小时 k个 ,然后转至步骤3。

如果 H(H) w个 ( 0 ) ( w个 ( k个 + 1 ) , λ k个 + 1 ) ϵ 2 ,然后让 小时 k个 + 1 =最大值{ 2 25 小时 0 ,( 小时 k个 /2)},k个=k个+1,然后转至步骤2。

步骤3。如果 λ k个 + 1 ϵ ,然后停止。否则,k个=k个+1,然后转至步骤1。

4结论

本文给出了一组无界条件,在此条件下,我们可以通过组合同伦内点算法求解一类无界非凸集上的非线性规划问题。本文提出的算法的主要优点是它是一种全局收敛的算法,与局部收敛的算法相比,它的初始点更容易选择。由于非线性规划问题在工程、管理、经济等领域有着广泛的应用,我们的结果可能有助于为处理这些非线性问题提供强大的解决工具。未来,我们将致力于提出新的技术来解决更广泛的无界非凸集类中的非线性规划问题。

工具书类

  1. Kellogg RB、Li TY、Yorke JA:Brouwer不动点定理的构造性证明和计算结果。 SIAM J.数字。分析。1976,13:473–483. 10.1137/0713041

    第条 数学科学网 数学 谷歌学者 

  2. 斯梅尔S:价格调整的收敛过程和全局牛顿法。 数学杂志。经济。1976,三:1–14. 10.1016/0304-4068(76)90002-1

    第条 数学科学网 数学 谷歌学者 

  3. Chow SN、Mallet-Paret J、Yorke JA:寻找映射零点:概率为1的同伦方法。 数学。计算。1978,32:887–899. 10.1090/S0025-5718-1978-0492046-9

    第条 数学科学网 数学 谷歌学者 

  4. 威斯康星州赞维尔加西亚CB:同伦与度理论的探讨。 数学。操作。物件。1979,4:390–405. 10.1287/门.4.3.90

    第条 数学科学网 数学 谷歌学者 

  5. 李毅、林志浩:Poincaré-Birkhoff定理的构造性证明。 事务处理。美国数学。Soc公司。1995,347:2111–2126.

    第条 数学 谷歌学者 

  6. 梅吉多N:线性规划中最优集的路径。数学规划、内点及相关方法的进展.编辑:Megiddo N.Springer,纽约;1988:131–158.

    谷歌学者 

  7. 小岛M、瑞野S、吉时A:线性规划的原对偶内点算法。内点及相关方法.编辑:Megiddo N.Springer,纽约;1988:29–47.

    谷歌学者 

  8. Kortanek KO、Potra F、Ye Y:非线性凸规划的一些有效内点算法。 线性代数应用。1991,152:169–189.

    第条 数学科学网 数学 谷歌学者 

  9. 麦考密克集团总经理:凸规划的投影SUMT方法。 数学。操作。物件。1989,14:203–223. 10.1287/门14.2.203

    第条 数学科学网 数学 谷歌学者 

  10. 蒙泰罗RDC,阿德勒一世:将Karmarkar型算法推广到一类具有全局线性收敛速度的凸可分规划问题。 数学。操作。物件。1990,15:408–422. 10.1287/门15.3.408

    第条 数学科学网 数学 谷歌学者 

  11. 朱JA:一类凸规划问题的路径跟踪算法。 ZOR-Methods模型操作。物件。1992,36:359–377. 2007年10月10日/BF01416235

    第条 数学科学网 数学 谷歌学者 

  12. 冯总、林志浩、于斌:非凸规划问题的Karush-Kuhn-Tucker点内部路径的存在性。 非线性分析。1998,32(6):761–768. 10.1016/S0362-546X(97)00516-6

    第条 数学科学网 数学 谷歌学者 

  13. Lin ZH、Yu B、Feng GC:求解凸非线性规划的一种组合的内点法。 申请。数学。计算。1997,84:193–211. 10.1016/S0096-3003(96)00086-0

    第条 数学科学网 数学 谷歌学者 

  14. 冯GC,余B:非线性规划问题的组合同伦内点法。数值与应用分析课堂讲稿14。数值数学的进展:第二届中日数值数学研讨会论文集.编辑:Fujita H,Yamaguti M.Kinokuniya,东京;1995:9–16.

    谷歌学者 

  15. 于斌、徐强、冯总经理:关于凸规划的组合同伦内点方法的复杂性。 J.计算。申请。数学。2007,200(1):32–46. 2016年10月10日/j.cam.2005.12.018

    第条 数学科学网 数学 谷歌学者 

  16. Su ML、Lv XR:利用同伦延拓方法求解一类非线性规划问题。 东北部。数学。J。2008,24(3):265–274.

    数学科学网 谷歌学者 

  17. Liu QH、Yu B、Feng GC:拟正规锥条件下非凸非线性规划问题的组合同伦内点方法。 数学学报。申请。罪。2003,26:372–377.

    谷歌学者 

  18. Xu Q、Yu B、Feng GC、Dang CY:无界集上变分不等式问题的同伦方法的全局收敛条件。 最佳方案。方法软件。2007,22:587–599. 10.1080/10556780600887883

    第条 数学科学网 数学 谷歌学者 

  19. 徐Q、林志浩:无界集上的组合同伦收敛。 数学学报。申请。罪。2004,27:624–631.

    数学科学网 谷歌学者 

  20. Allgower EL,Georg K:数值延拓方法简介SIAM,费城;2003

     数学 谷歌学者 

下载参考资料

致谢

本研究得到了NSFC-河南联合科学基金(No.U1304103)和河南省自然科学基金(No 122300410261)的资助。

作者信息

作者和附属机构

作者

通讯作者

与的通信苏梦龙.

其他信息

竞争性利益

作者声明,他们没有相互竞争的利益。

作者的贡献

MS完成了主要工作并起草了手稿。JL参与完成了引理2.3的证明。ZX参与完成了定理2.1的证明。所有作者阅读并批准了最终手稿。

权利和权限

开放式访问本文根据Creative Commons Attribution 4.0 International License的条款分发(https://creativecommons.org/licenses/by/4.0)允许以任何媒体或格式使用、复制、改编、分发和复制,只要您对原始作者和来源给予适当的信任,提供知识共享许可的链接,并指出是否进行了更改。

转载和许可

关于这篇文章

检查更新。通过CrossMark验证货币和真实性

引用这篇文章

Su,M.,Li,J.&Xu,Z。通过全局收敛算法求解具有无界非凸约束集的非线性规划问题。J不平等申请 2014, 281 (2014). https://doi.org/10.1186/1029-242X-2014-281

下载引文

  • 收到:

  • 认可的:

  • 出版:

  • 内政部:https://doi.org/10.1186/1029-242X-2014-281

关键词