Two-level primal–dual proximal decomposition technique to solve large scale optimization problems

We describe a primal–dual application of the proximal point algorithm to nonconvex minimization problems. Motivated by the work of Spingarn and more recently by the work of Hamdi et al. about the primal resource-directive decomposition scheme to solve nonlinear separable problems. This paper discusses some local results of a primal–dual regularization approach that leads to a decomposition algorithm.     

求解第一类Fredholm积分方程的一种新的正则化算法

应用一种新的正则化方法建立了一类新的求解第一类Fredholm积分方程的正则化算法, 并借助Matlab软件给出了数值算例.数值结果与理论分析基本一致,而且表明文中建立的正则化比通常的Tikhonov正则化更精确.     

啁啾激光场中HF分子的经典解离

应用经典理论并采用辛算法，计算了HF分子在啁啾激光作用下的经典解离.讨论了激光场中HF分子的解离概率随时间的变化以及考虑振动-转动能级跃迁对分子解离概率的影响.计算结果与理论分析相符，说明运用经典理论并采用辛格式计算双原子分子在激光作用下的经典解离是有效和可行的. 关键词： 啁啾激光场 经典解离 辛算法     

Identification of source terms in wave equation with dynamic boundary conditions

This paper studies an inverse hyperbolic problem for the wave equation with dynamic boundary conditions. It consists of determining some forcing terms from the final overdetermination of the displacement. First, the Fréchet differentiability of the Tikhonov functional is studied, and a gradient formula is obtained via the solution of an associated adjoint problem. Then, the Lipschitz continuity of the gradient is proved. Furthermore, the existence and the uniqueness for the minimization problem are discussed. Finally, some numerical experiments for the reconstruction of an internal wave force are implemented via a conjugate gradient algorithm.     

Two-Phase Image Segmentation by Nonconvex Nonsmooth Models with Convergent Alternating Minimization Algorithms

Two-phase image segmentation is a fundamental task to partition an image into foreground and background. In this paper, two types of nonconvex and nonsmooth regularization models are proposed for basic two-phase segmentation. They extend the convex regularization on the characteristic function on the image domain to the nonconvex case, which are able to better obtain piecewise constant regions with neat boundaries. By analyzing the proposed non-Lipschitz model, we combine the proximal alternating minimization framework with support shrinkage and linearization strategies to design our algorithm. This leads to two alternating strongly convex subproblems which can be easily solved. Similarly, we present an algorithm without support shrinkage operation for the nonconvex Lipschitz case. Using the Kurdyka-Łojasiewicz property of the objective function, we prove that the limit point of the generated sequence is a critical point of the original nonconvex nonsmooth problem. Numerical experiments and comparisons illustrate the effectiveness of our method in two-phase image segmentation.     

On a Class of Quasilinear Elliptic Equations

We consider a class of quasilinear elliptic boundary problems, including the following Modified Nonlinear Schrödinger Equation as a special case: $$\begin{cases} ∆u+ \frac{1}{2} u∆(u^2)−V(x)u+|u|^{q−2}u=0 \ \ \ in \ Ω, \\u=0 \ \ \ \ \ \ \ ~ ~ ~ on \ ∂Ω, \end{cases}$$ where $Ω$ is the entire space $\mathbb{R}^N$ or $Ω ⊂ \mathbb{R}^N$ is a bounded domain with smooth boundary, $q∈(2,22^∗]$ with $2^∗=2N/(N−2)$ being the critical Sobolev exponent and $22^∗= 4N/(N−2).$ We review the general methods developed in the last twenty years or so for the studies of existence, multiplicity, nodal property of the solutions within this range of nonlinearity up to the new critical exponent $4N/(N−2),$ which is a unique feature for this class of problems. We also discuss some related and more general problems.     

SYMPLECTIC STRUCTURE OF POISSON SYSTEM

When the Poisson matrix of Poisson system is non-constant, classical symplectic methods, such as symplectic Runge-Kutta method, generating function method, cannot preserve the Poisson structure. The non-constant Poisson structure was transformed into the symplectic structure by the nonlinear transform.Arbitrary order symplectic method was applied to the transformed Poisson system. The Euler equation of the free rigid body problem was transformed into the symplectic structure and computed by the midpoint scheme. Numerical results show the effectiveness of the nonlinear transform.     

The annular twist mapping of the restricted three-body problem

This paper presents a new derivation of the twist mapping in the planar restricted problem. It differs from other treatments in the use of a novel canonical transformation which allows for the utilization of symplectic reduction techniques.     

Reisner板的哈密尔顿体系及其辛正交解析法

基于Ｒｅｉｓｎｅｒ板理论，通过对混合能变分原理的修正，建立了更一般的哈密尔顿型广义变分原理，并给出了Ｒｅｉｓｎｅｒ板问题的哈密尔顿正则方程及其共轭辛正交解析法