求解非线性方程组的一种新的全局收敛的Levenberg-Marquardt算法

本文提出了求解非线性方程组的一种新的全局收敛的Levenberg-Marquardt算法,即μk=ακ(θ||F_k|| (1-θ)||J_k~TF_k||),θ∈[0,1],其中ακ利用信赖域技巧来修正.在不必假设雅可比矩阵非奇异的局部误差界条件下,证明了该算法是全局收敛和局部二次收敛的.数值试验表明该算法能有效地求解奇异非线性方程组问题.     

一种求解非线性互补问题的多步自适应Levenberg-Marquardt算法

本文采用Modulus-based变换将非线性互补问题转化为非光滑方程组,并将一种多步自适应Levenberg-Marquardt方法推广应用于求解所得的非光滑方程组,从而得到原问题的解.在适当条件下,本文证明了算法的全局收敛性.与一种已有的参数自适应Levenberg-Marquardt方法(PSA-LMM)相比较,...     

一类推广的拟牛顿多进程异步并行方法

本文通过近似雅可比矩阵Bk代替雅可比矩阵F′(xk),运用多进程异步并行方法求解非线性方程组。该方法在保持解的精度的情况下,缩短了运行时间和迭代步数。文中给出了算法收敛性的证明及八个非线性方程组的数值测试结果,表明该算法是可行的和快速的。     

对称不定问题的不精确Newton法

1.引 言 非线性方程组F（x）=0的数值求解,经典的算法是Newton迭代;xk 1=xk sk,k=0,1,2,…,（1.1）其中的sk满足F’（xk）sk=-F（xk）;k=0,1,2,….（1.2）这里x0为迭代的初始点,{xk}称为Newton迭代序列.当变量个数比较多时,每一步Newton迭代中计算Jacobi矩阵F’（xk）和求解线性方程组（1.2）的代价非常高;特别当xk远离方程组的解x＊时,高精度地求解线性方程组（1.2）     

ABS投影算法的两个改进

本文将求解线性方程的ABS投影算法进行两方面的改进和推广，一是使算法在第K次迭代产生的点xk＋1不仅满足前k个方程，还尽可能地使得在点xk处成立的方程j（j＞k）在xk＋1处仍成立，称之为强ABS投影算法，另外初始选代矩阵由非奇异的减弱为任意的．二是建立了系数矩阵有零子块的方程组的ABS投影算法，其存贮量和计算量比原ABS投影算法小．ABS算法可以作为这两种改进算法的特别情形．     

求解非线性方程组的L-M方法

本文研究了求解奇异非线性方程组的Levenberg-Marquardt方法的收敛性.利用选取新的迭代参数求解非线性方程组的L-M方法,获得点列的超线性收敛性和二阶收敛性,并把试验结果与文献[19,20]的结果进行了比较.     

一种解非线性方程组的不精确牛顿——兰索斯方法

本文提出一种不完全线搜索技术的不精确牛顿—克雷洛夫(Newton-Krylov)子空间方法解对称非线性方程组,其中克雷洛夫子空间方法采用的是兰索斯(Lanczos)类分解技术.迭代方向是通过使用兰索斯方法近似求解非线性方程组的牛顿方程获得的.在合理的假设条件下,分析了算法的全局收敛性和局部超线性收敛速率.最后,数值结果显示了该算法的有效性.     

解非线性最小二乘问题的ABS方法

在本文中,基于解非线性方程组的ABS方法的思想,我们对非线性最小二乘问题建立了一类新的算法。在类似于Gauss-Newton法的收敛条件下,我们证明了算法的局部收敛性。此外,在对算法结构进行深入分析的基础上,我们将新算法转化为一种近似Gauss-Newton法。并建立了它的Kantorovich型收敛定理。数值结果表明ABS算法是有效的,且在一定程度上优越于Gauss-Newton法。     

求解广义非线性互补问题的光滑化拟牛顿法

利用光滑对称扰动Fischer-Burmeister函数将广义非线性互补问题转化为非线性方程组,提出新的光滑化拟牛顿法求解该方程组.然后证明该算法是全局收敛的,且在一定条件下证明该算法具有局部超线性(二次)收敛性.最后用数值实验验证了该算法的有效性.     

一类非线性方程组的Newton-PSS迭代法

正定反Hermite分裂(PSS)方法是求解大型稀疏非Hermite正定线性代数方程组的一类无条件收敛的迭代算法.将其作为不精确Newton方法的内迭代求解器,我们构造了一类用于求解大型稀疏且具有非Hermite正定Jacobi矩阵的非线性方程组的不精确Newton-PSS方法,并对方法的局部收敛性和半局部收敛性进行了详细的分析.数值结果验证了该方法的可行性与有效性.     

A REPRESENTATION FORMULA RELATED TO SCHR(O)DINGER OPERATORS

Schr(o)dinger operator is a central subject in the mathematical study of quantum mechanics.Consider the Schrodinger operator H = -△ V on R, where △ = d2/dx2 and the potential function V is real valued. In Fourier analysis, it is well-known that a square integrable function admits an expansion with exponentials as eigenfunctions of -△. A natural conjecture is that an L2 function admits a similar expansion in terms of "eigenfunctions" of H, a perturbation of the Laplacian (see [7], Ch. Ⅺ and the notes), under certain condition on V.     

On a class of self-similar processes with stationary increments in higher order Wiener chaoses

We study a class of self-similar processes with stationary increments belonging to higher order Wiener chaoses which are similar to Hermite processes. We obtain an almost sure wavelet-like expansion of these processes. This allows us to compute the pointwise and local Hölder regularity of sample paths and to analyse their behaviour at infinity. We also provide some results on the Hausdorff dimension of the range and graphs of multidimensional anisotropic self-similar processes with stationary increments defined by multiple Wiener–Itô integrals.     

Qualitative Analysis of Bifurcating Solutions in the Lengyel-Epstein Model

One of the most fundamental problems in theoretical biology is to explain the mechanisms by which patterns and forms are created in the＇living world. In his seminal paper ＂The Chemical Basis of Morphogenesis＂, Turing showed that a system of coupled reaction-diffusion equations can be used to describe patterns and forms in biological systems. However, the first experimental evidence to the Turing patterns was observed by De Kepper and her associates（1990） on the CIMA reaction in an open unstirred reactor, almost 40 years after Turing＇s prediction. Lengyel and Epstein characterized this famous experiment using a system of reaction-diffusion equations. The Lengyel-Epstein model is in the form as follows     

(0,1 ;0)-INTERPOLATION ON INFINITE INTERVAL (-∞, +∞)

In this paper, we study the explicit representation and convergence of (0, 1; 0)-interpolation on infinite interval, which means to determine a polynomial of degree ≤ 3n - 2 when the function values are prescribed at two set of points namely the zeros of Hn(x) and H′n(x) and the first derivatives at the zeros of H′n(x).     

On a class of Riemann surfaces

It is considered the class of Riemann surfaces with dimT1 = 0, where T1 is a subclass of exact harmonic forms which is one of the factors in the orthogonal decomposition of the spaceΩH of harmonic forms of the surface, namely The surfaces in the class OHD and the class of planar surfaces satisfy dimT1 = 0. A.Pfluger posed the question whether there might exist other surfaces outside those two classes. Here it is shown that in the case of finite genus g, we should look for a surface S with dimT1 = 0 among the surfaces of the form Sg\K , where Sg is a closed surface of genus g and K a compact set of positive harmonic measure with perfect components and very irregular boundary.     

CONTENTS OF VOLUME 30

正~~     

Instructions for contributions

正Applied Mathematics-A Journal of Chinese Universities,Series B(Appl.Math.J.Chinese Univ.,Ser.B)is a comprehensive applied mathematics journal jointly sponsored by Zhejiang University,China Society for Industrial and Applied Mathematics,and Springer-Verlag.It is a quarterly journal with     

Journal of Mathematical Research with Applications Guide for Authors

正Journal overview:Journal of Mathematical Research with Applications(JMRA),formerly Journal of Mathematical Research and Exposition(JMRE)created in 1981,one of the transactions of China Society for Industrial and Applied Mathematics,is a home for original research papers of the highest quality in all areas of mathematics with applications.The target audience comprises:pure and applied mathematicians,graduate students in broad fields of sciences and technology,scientists and engineers interested in mathematics.     

Staircase transportation problems with superadditive rewards and cumulative capacities

A cumulative-capacitated transportation problem is studied. The supply nodes and demand nodes are each chains. Shipments from a supply node to a demand node are possible only if the pair lies in a sublattice, or equivalently, in a staircase disjoint union of rectangles, of the product of the two chains. There are (lattice) superadditive upper bounds on the cumulative flows in all leading subrectangles of each rectangle. It is shown that there is a greatest cumulative flow formed by the natural generalization of the South-West Corner Rule that respects cumulative-flow capacities; it has maximum reward when the rewards are (lattice) superadditive; it is integer if the supplies, demands and capacities are integer; and it can be calculated myopically in linear time. The result is specialized to earlier work of Hoeffding (1940), Fréchet (1951), Lorentz (1953), Hoffman (1963) and Barnes and Hoffman (1985). Applications are given to extreme constrained bivariate distributions, optimal distribution with limited one-way product substitution and, generalizing results of Derman and Klein (1958), optimal sales with age-dependent rewards and capacities.To our friend, Philip Wolfe, with admiration and affection, on the occasion of his 65th birthday.Research was supported respectively by the IBM T.J. Watson and IBM Almaden Research Centers and is a minor revision of the IBM Research Report [6].