结合修正Curry-Altman步长搜索一个新的共轭梯度算法的收敛性

Conjugate gradient optimization algorithms depend on the search directions with different choices for the parameter in the search directions. In this note, conditions are given on the parameter in the conjugate gradient directions to ensure the descent property of the search directions. Global convergence of such a class of methods is discussed. It is shown that, using reverse modulus of continuity function and forcing function, the new method for solving unconstrained optimization can work for a continuously differentiable function with a modification of the Curry-Altman‘s step-size rule and a bounded level set. Combining PR method with our new method, PR method is modified to have global convergence property.Numerical experiments show that the new methods are efficient by comparing with FR conjugate gradient method.     

Geometric interpretation of several classical iterative methods for linear system of equations and diverse relaxation parameter of the SOR method

Two kinds of iterative methods are designed to solve the linear system of equations, we obtain a new interpretation in terms of a geometric concept. Therefore, we have a better insight into the essence of the iterative methods and provide a reference for further study and design. Finally, a new iterative method is designed named as the diverse relaxation parameter of the SOR method which, in particular, demonstrates the geometric characteristics. Many examples prove that the method is quite effective.     

线性互补问题的一种新Lagrange乘子法

A new multiplier method for solving the linear complementarity problem LCP(q, M) is proposed. Based on the Lagrangian of LCP(q,M) introduced here, we construct a new differentiable merit function θ(x,λ) which containing a multiplier vector λ and satisfying θ(x,λ) ≥ 0 and θ(x,λ) = 0 if and if only x solves LCP(q,M). A simple damped Newton-type algorithm which based on the merit function θ(x,λ) is presented. The main feature of the method is that the multiplier self-adjusting step accelerates the local convergence rate without losing global convergence. When M is the P-matrix, the sequence {θ(x^k,λ^k)}where {(x^k,λ^k)} generated by the algorithm is globally linearly convergent to zero and convergent in finite number of iterations if the solution is nondegenerate. Numerical results suggest that the method is high efficient and promising.     

A weak condition for secant method to solve systems of nonlinear equations

In this paper, a new weak condition for the convergence of secant method to solve the systems of nonlinear equations is proposed. A convergence ball with the center x0 is replaced by that with xl, the first approximation generated by the secant method with the initial data x-1 and x0. Under the bounded conditions of the divided difference, a convergence theorem is obtained and two examples to illustrate the weakness of convergence conditions are provided. Moreover, the secant method is applied to a system of nonlinear equations to demonstrate the viability and effectiveness of the results in the paper.     

Restricted summability of Fourier transforms and local Hardy spaces

A general summability method, the so-called θ-summability is considered for multi-dimensional Fourier transforms. Under some conditions on θ, it is proved that the maximal operator of the θ-means defined in a cone is bounded from the amalgam Hardy space W（hp, e∞） to W（Lp,e∞）. This implies the almost everywhere convergence of the θ-means in a cone for all f ∈ W（L1, e∞） velong to L1.     

Composite penalty method of a low order anisotropic nonconforming finite element for the Stokes problem

Composite penalty method of a low order anisotropic nonconforming quadrilateral finite element for the Stokes problem is presented. This method with a large penalty parameter can achieve the same accuracy as the stand method with a small penalty parameter and the convergence rate of this method is two times as that of the standard method under the condition of the same order penalty parameter. The superconvergence for velocity is established as well. The results of this paper are also valid to the most of the known nonconforming finite element methods.     

NONCONFORMING FINITE ELEMENT PENALTY METHOD FOR STOKES EQUATION

a special penalty method is presented to improve the accuracy of the standard penaltymethod (or solving Stokes equation with nonconforming finite element, It is shown that thismethod with a larger penalty parameter can achieve the same accuracy as the staodaxd methodwith a smaller penalty parameter. The convergence rate of the standard method is just hall order of this penalty method when using the same penalty parameter, while the extrapolationmethod proposed by Faik et al can not yield so high accuracy of convergence. At last, we alsoget the super-convergence estimates for total flux.     

GLOBAL CONVERGENCE RESULTS OF A THREE TERM MEMORY GRADIENT METHOD WITH A NON-MONOTONE LINE SEARCH TECHNIQUE

In this paper, a new class of three term memory gradient method with nonmonotone line search technique for unconstrained optimization is presented. Global convergence properties of the new methods are discussed. Combining the quasi-Newton method with the new method, the former is modified to have global convergence property. Numerical results show that the new algorithm is efficient.     

ON A FAMILY OF CHEBYSHEV-HALLEY TYPE METHODS IN BANACH SPACE UNDER WEAKER SMALE CONDITION

In this paper, we discuss local convergence of a family of Chebychev-Halley type methods with a parameter θ∈[0,1] in Banach space using Smale-type δ criterion under 2-th γ-condition. We will see that the properties of the condition used for local convergence is much more different from that used in [6][15] for the semi-local convergence.     

Convergence analysis for the Secant method based on new recurrence relations

A new convergence theorem for the Secant method in Banach spaces based on new recurrence relations is established for approximating a solution of a nonlinear operator equation. It is assumed that the divided difference of order one of the nonlinear operator is Lipschitz continuous. The convergence conditions differ from some existing ones and are easily satisfied. The results of the paper are justified by numerical examples that cannot be handled by earlier works.     

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).     

RECENT PROGRESS ON SPHERICAL HARMONIC APPROXIMATION MADE BY BNU RESEARCH GROUP ——In memory of Professor Sun Yongsheng

As early as in 1990, Professor Sun Yongsheng, suggested his students at Beijing Normal University to consider research problems on the unit sphere. Under his guidance and encouragement his students started the research on spherical harmonic analysis and approximation. In this paper, we incompletely introduce the main achievements in this area obtained by our group and relative researchers during recent 5 years （2001-2005）. The main topics are： convergence of Cesaro summability, a.e. and strong summability of Fourier-Laplace series; smoothness and K-functionals; Kolmogorov and linear widths.     

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.     

Journal of the Operations Research Society of China Special issue on Operations Research in Healthcare Management

正Guest Editors:Hong Chen,Shanghai Jiao Tong University,Shanghai,China Guohua Wan,Shanghai Jiao Tong University,Shanghai,China David Yao,Columbia University,New York,USA Scope:Healthcare delivery worldwide has been fraught with high cost,low efficiency and poor quality of patient care service.For the field of operations research(OR),healthcare offers some of the biggest challenges as well as best opportunities in     

ANALYSIS IN THEORY AND APPLICATIONS Vol. 30, 2014 CONTENTS

正~~     

An alternative mathematical algorithm for the photo- and videokeratoscope

Due to the resolution of current laser technology, the accuracy of corneal topography as measured by the videokeratoscope is no longer adequate to provide precise enough data for refractive surgery or for the fitting of customized contact lenses. We present an algorithm for recovering corneal topography that makes use of modern differential geometric techniques and numerical descent in Sobolev spaces. We believe this algorithm may be used with the photo- and videokeratoscope to increase the accuracy of the recovered corneal topography.     

A note on Sobolev orthogonality for Laguerre matrix polynomials

Let {Ln(A,λ)(x)}n≥0 be the sequence of monic Laguerre matrix polynomials defined on [0, ∞) by Ln(A,λ)(x)=n!/(-λ)n∑nk=0(-λ)κ/k!(n-1)! (A I)n[(A I)k]-1 xk,where A ∈ Cr×r. It is known that {Ln(A,λ)(x)}n≥0 is orthogonal with respect to a matrix moment functional when A satisfies the spectral condition that Re(z) ＞ - 1 for every z ∈σ(A).In this note we show that forA such that σ(A) does not contain negative integers, the Laguerre matrix polynomials Ln(A,λ) (x) are orthogonal with respect to a non-diagonal SobolevLaguerre matrix moment functional, which extends two cases: the above matrix case and the known scalar case.     

CONTENTS OF VOLUME 30

正~~