广义Newton法求解非线性互补问题

１引言考虑非线性互补问题ＮＣＰ（ｆ）：的求解，即我们要寻求某ｘ∈Ｒｎ，使其满足（１．１）．其中映射ｆ：Ｒｎ→Ｒｎ为具有连续Ｆ－导数的非线性映射．众所周知，问题（１．ｌ）可以等价地转化为Ｂ－可微方程组：求解，其中：容易证明，由（１．３）定义的映射Ｇ处处Ｂ－可微，且其在点ｘ∈Ｒｎ处的Ｂ－导数ＢＧ（ｘ）为而对于问题（１．２）（１．３），我们希望直接用经典的广义Ｎｅｗｔｏｎ法进行求解．但是，由于由（１．３）定义映射Ｇ在（１．１）的解ｘ∈Ｒｎ处，没有可逆的强Ｆ－导数存在，因此，关于算法（１．５）（１．６）…     

求解不等式约束最优化问题的MBF—ODE方法

１引言考虑用基于修正内罚函数的常微分方程（ＭＢＦ－ＯＤＥ）方法求解下列不等式约束极小化问题：其中ｆi∈ｃ２：Ｒ,ｉ=0,１,…,ｍ．求解无约束极小化问题的ＯＤＥ的一般形式是其中,φ（ｘ）∈Ｃ１：ΩＲｎ→Ｒ;ｓ（ｘ）∈Ｃ１：ΩＲｎ→Ｒｎ且满足φ（ｘ）＞０,ｓＴ（ｘ）ｆ（ｘ）＜０,ｆ（ｘ）∈Ｃ１：Ｒｎ→Ｒ为目标函数．为便于用ＯＤＥ方法求解（１．ｌ）,可藉助于罚函数将（１．ｌ）变换为无约束极小化问题（见［７］．但由于经典罚函数（ＣＢＦ）在计算上有较大的困难,我们采用修正内罚函数（ＭＢＦ）．其基本思想是用…     

定步长的连续极小化方法...

求解非线性方程组Ｆ（ｘ）＝０可转化为求非线性最小二乘问题ｍｉｎ１／２Ｆ（ｘ）＾ＴＦ（ｘ）的极小点。文章提出了一种求解上述非线性最小二乘问题的连续极小化方法，方法给出确定的步长，并证明具有整体和线性的收敛性。两个数值例子说明了方法的优越性。     

周期问题分类浅析

近年在各级各类的试卷中,常见与周期有关的试题．由于周期函数在中学课本中仅介绍定义．练习中也只有求三角函数的周期,除此之外,在课本中再也找不到与周期有关的问题．因此,涉及周期问题的求解,学生颇感困难,本文以近年各类试题为例,对周期问题从周期给出的方式入手进行分类浅析．１　直接给出周期周期为Ｔ的函数,直接给出周期的方式有：ｆ（ｘ＋Ｔ）＝ｆ（ｘ）、ｆ（ｘ＋Ｔ２）＝－ｆ（ｘ）、ｆ（ｘ＋Ｔ２）＝１ｆ（ｘ）、ｆ（ｘ＋Ｔ２）＝－１ｆ（ｘ）等．例１　若ｆ（ｘ）（ｘ∈Ｒ）是以２为周期的偶函数,当ｘ∈［０,１］时,…     

非光滑约束问题的既约次梯度法

１引言 对带约束的不可微的非线性规划问题,由于不能使用梯度,求极小点就比较困难．本文给出解决此问题的一种有效的算法． ２　非光滑约束问题的既约次梯度法 １）非线性规划问题的Ｌａｅｒａｎｅ对偶理论 考虑下面非线性规划问题其中ｇ（ｘ）＝（ｇ１（ｘ）,…,ｇｒ（ｘ））Ｔ,ｈ（ｘ））＝（ｈ１（ｘ）,…,ｈｍ（ｘ））Ｔ,ｆ（ｘ）＝　　　　　　Ｒｎ中是Ｌｉｓｐｓｃｈｉｔｚ连续的ｉ＝１,２,…,ｒ,ｊ＝１,２,…,ｍ相应的Ｌａｇｒａｎｇｅ对偶问题为其中　　（ｕ,　）＝ｉｎｆＬ（ｘ;ｕ,ｖ）＝ｉｎｆ（ｆ（ｘ）＋ｕＴ…     

再论求导数零点的二次收敛迭代法

一维搜索是最优化理论数值计算的一个基本问题,它可归结为求定义在开凸区域Ｄ上的可微函数 ｆ的导数零点．若用 Ｎｅｗｔｏｎ法求导数零点,则涉及到二阶导数的计算．若用带导数的三次插值法则需要开平方的计算［１］．为了克服上述问题,本文作者之一在 １９７９年［２］首次提出了下述具有二阶收敛速度的迭代法：通常,我们称迭代法（０．１）为基于信息集（ｆ（ｘｎ）,ｆ’（ｘｎ）,ｆ（ｘｎ－１）,ｆ’（ｘｎ－１）｝的迭代法,而δ（ｆｘｙ）是基于信息集｛ｆ（ｘ）,ｆ＇（ｘ）,ｆ（ｙ）,Ｆ＇（ｙ））｝的三次插值多项式在ｘ处…     

求解简单界约束优化问题的一种逐次逼近法

１引言考虑变量带简单界约束的非线性规划问题：其中二阶连续可微,ａ＝（ａ１,ａ２,…,ａｎ）,ｂ＝（ｂ１,ｂ２,…,ｂｎ）,＋ｉ＝１,２,…,ｎ．问题（１）不仅是实际应用中出现的简单界约束最优化问题,而且相当一部分最优化问题可以把变量限制在有意义的区间内（参见［１］）．因此无论在理论方面还是在实际应用方面,都有研究此类问题并给出简便而有效算法的必要．假设ｆ是凸函数,记ｇ（ｘ）＝ｆ（ｘ）,则由Ｋ－Ｔ条件,问题（１）可化为求解下面的非光滑方程组：显然,（２）等价于易证,（３）等价于求解下面的非光滑方程…     

新题征展(3)

题组新编１．（１）设Ｍ＝｛ｘ｜ｆ（ｘ）＝０｝、Ｎ＝｛ｘ｜ｇ（ｘ）＝０｝,则｛ｘ｜ｆ（ｘ）·ｇ（ｘ）＝０｝为（　　）;（Ａ）Ｍ　（Ｂ）Ｎ　（Ｃ）Ｍ∪Ｎ　（Ｄ）以上都不对（２）设ｆ（ｘ）＝ｘ－１ｘ＋３,ｇ（ｘ）＝ｘ＋３ｘ－１,则集合｛ｘ｜ｆ（ｘ）·ｇ（ｘ）＝０｝＝　　;（３）设函数ｆ（ｘ）、ｇ（ｘ）的定义域依次是Ｆ、Ｇ,且Ｍ＝｛ｘ｜ｆ（ｘ）＝０｝、Ｎ＝｛ｘ｜ｇ（ｘ）＝０｝,则｛ｘ｜ｆ（ｘ）·ｇ（ｘ）＝０｝＝　　．２．（１）设ｍ、ｋ∈Ｎ,则Ｃｎｎ＋Ｃｎｎ＋１＋Ｃｎｎ＋２＋…＋Ｃｎｎ＋ｋ＝　　;（２）求…     

关于微积分基本公式的一点注记

大多数分析教材将微积分基本公式叙述为：定理１　（ｉ）ｆ（ｘ）在［ａ,ｂ］上连续;（ｉｉ）Ｆ（ｘ）是ｆ（ｘ）的任意一个原函数,则　　∫ｂａｆ（ｘ）ｄｘ＝Ｆ（ｂ）－Ｆ（ａ）．某些教科书将定理１的条件减弱,改述定理１为：定理２　（ｉ）ｆ（ｘ）在［ａ,ｂ］上可积;（ｉｉ）存在Ｆ（ｘ）在［ａ,ｂ］上连续,在［ａ,ｂ］－Ａ（Ａ为［ａ,ｂ］的一有限子集）上Ｆ′（ｘ）＝ｆ（ｘ）,则∫ｂａｆ（ｘ）ｄｘ＝Ｆ（ｂ）－Ｆ（ａ）．我们知道,黎曼函数Ｒ（ｘ）＝１ｑ,ｘ＝ｐｑ,ｑ＞０,ｐ,ｑ互质,０,ｘ为无理数．在［ａ,ｂ］…     

复Schr(o)dinger场和实Klein-Gordon场相互作用下一类方程组守恒差分格式的收敛性和稳定性

１引言关于Ｚａｋｈａｒｏｖ方程各种定解问题的讨论,近年来已发表了许多的文献［１－５］,很明显,这是复非线性Schrod-inger（NLS）方程和实非线性波动方程耦合的一类方程组．文献［６］考察了复ＮＬＳ方程和实非线性Ｋｌｅｉｎ－Ｇｏｒｄｏｎ（Ｇ－Ｋ）方程耦合的一类方程组的各种定解问题,其中φ（ｘ,ｔ）是复值函数,ｕ（ｘ,ｔ）是实值函数,ｆ（ｓ）∈（－∞,∞）,μ,ｇ,ｍ均为正常数,ｉ２＝－１．文中证明了整体解的存在性、唯一性和存在如下两个守恒律其中Ｆ（ｓ）＝ｆ（ｚ）ｄｚ,Ｅ００,Ｅ０１均为仅与初始条件有关…     

半无限规划问题的一个有效解法

1．引言在计算机辅助设计和工程设计中,经常遇到下面的两类优化问题1,2］．1．无约束半无限极大极小问题．其中外x）二——x。。Im。x。。。Yi夕（x;N）这里J二（】,2,·,｝对任何7E八岁：R－xR”。＋R是连续可微的函数,X是R”。中的一个紧子集,且VYj）一O,这里问h）表示X体积．2．约束半无限代化问题．其中I一（1,2,·．小记L二《0｝UI对任何jCL冲’（x）一max。。。Yi夕（x,yi）·这里拉：PX＊n＋R是连续可微函数,X是”。中的一个紧子集,且NU）一0·注．设Y（Z,一二切EyW一叫卜4．今后对本文用到的紧子集地做…     

一种改进的混沌序列产生方法

根据灰度图像的二维直方图的特点,在已有的二维Arnold混沌系统的基础上,结合Bernstein形式的Bézier曲线的生成算法,给出了一种基于生成Bézier曲线的de Casteljau算法构造伪随机序列的方法,实验结果表明生成的二维序列不仅具有伪随机性,而且还具有在近似圆盘中随机分布的性质,这使得该伪随机序列更适合对灰度图像的二维灰度直方图进行基于混沌优化的图像分割.在此基础上,给出了一种基于混沌优化的二维最大熵的灰度图像分割算法,该算法对于含噪图像取得了良好的分割效果.     

多目标minimax问题的极大熵逼近收敛性

本文利用极大熵逼近函数，展开了多目标ｍｉｎｉｍａｘ问题的逼近方法的研究，并讨论了该逼近方法的收敛性，所得结果是目前已有的结果进一步拓广．     

基于离差最大化的挤奶机方案综合评价

针对以往挤奶设备方案评价过程中权重确定存在的主观随意性,提出了一种新的基于离差最大化方法的挤奶机方案综合评价方法.首先从经济性、技术性和质量三个角度建立了多层次的挤奶机方案综合评价指标体系,然后介绍了属性权重完全未知时利用离差最大化方法确定属性权重的基本原理.在此基础上,建立了基于离差最大化方法的挤奶机方案综合评价模型.最后通过仿真实例对三型挤奶机方案进行了评价和分析,验证了方法的可行性和有效性.     

A NEW NUMERICAL METHOD FOR TWO-PHASE IMMISCIBLE INCOMPRESSIBLE PROBLEM

Two-phase, immiscible, incompressible flow in porous media is governed by a system of nonlinear partial differential equations. In most practical applications convection physically dominates diffusion, and the object of this paper is to develop a finite difference method combined with the method of characteristics and the lumped mass method to treat the parabolic equation of the differential system. This method is shown satisfy the maximum principle and its error analysis is presented.     

工程设计中约束规划熵方法的收敛性分析

极大熵方法在工程设计优化中得到成功的应用,但它的收敛性分析一直没有得到很好的解决．本文讨论了这个有意义的问题,在一般连续条件下解决了工程设计中的外点极大熵方法和内点极大熵方法的收敛性．     

Binary clustering with missing data

A clustering method is presented for analysing multivariate binary data with missing values. When not all values are observed, Govaert³ has studied the relations between clustering methods and statistical models. The author has shown how the identification of a mixture of Bernoulli distributions with the same parameter for all clusters and for all variables corresponds to a clustering criterion which uses L₁ distance characterizing the MNDBIN method (Marchetti⁸). He first generalized this model by selecting parameters which can depend on variables and finally by selecting parameters which can depend both on variables and on clusters. We use the previous models to derive a clustering method adapted to missing data. This method optimizes a criterion by a standard iterative partitioning algorithm which removes the necessity either to ignore objects or to substitute the missing data. We study several versions of this algorithm and, finally, a brief account is given of the application of this method to some simulated data.     

交替迭代最大似然估计量的偏倚校正

In this paper, we give a definition of the alternating iterative maximum likelihood estimator (AIMLE) which is a biased estimator. Furthermore we adjust the AIMLE to result in asymptotically unbiased and consistent estimators by using a bootstrap iterative bias correction method as in Kuk (1995). Two examples and simulation results reported illustrate the performance of the bias correction for AIMLE.     

On a class of logarithmic integrals

In this note we give closed forms for a class of logarithmic integrals in terms of Bernoulli polynomials. This provides a method for unifying a large class of definite integrals.     

A polynomial-time simplex method for the maximumk-flow problem

A generalization of the maximum-flow problem is considered in which every unit of flow sent from the source to the sink yields a payoff of $k. In addition, the capacity of any arce can be increased at a per-unit cost of $c _e . The problem is to determine how much arc capacity to purchase for each arc and how much flow to send so as to maximize the net profit. This problem can be modeled as a circulation problem. The main result of this paper is that this circulation problem can be solved by the network simplex method in at mostkmn pivots. Whenc _e = 1 for each arce, this yields a strongly polynomial-time simplex method. This result uses and extends a result of Goldfarb and Hao which states that the standard maximum-flow problem can be solved by the network simplex method in at mostmn pivots.Research partially supported by Office of Naval Research Grant N00014-86-K-0689 at Purdue University.