Stabilization of the motions of mechanical systems in prescribed phase-space manifolds

A method for constructing a mathematical model of the dynamics of a mechanical system is proposed. An algorithm is constructed for determining the expressions for the control forces and the components of the constraint reactions. A modification is made to the dynamic equations which enables one to solve the problem of stabilizing the constraints and which ensures the required accuracy in the numerical solution of the corresponding system of differential-algebraic equations describing the constraints imposed on the system, its kinematics and dynamics. By virtue of well-known dynamic analogies, the proposed method can be used to investigate the dynamics of different physical systems. The problem of modelling the dynamics of an adaptive optical system with two degrees of freedom is considered.     

一类齐次多项式优化问题的一个全局优化算法

对图像与信号处理中遇到的一类齐次多项式优化问题,本文首先借助平移技术将目标函数转化为凸函数,然后结合初始点技术提出了求解该类问题的一个全局优化算法.与求解该类问题的幂方法相比,本文给出的方法不但能在一般情形下保证算法的全局收敛性,而且数值结果表明在多数情况下可以得到问题的一个全局最优值解.     

A practical method for computing the largest M‐eigenvalue of a fourth‐order partially symmetric tensor

In this paper, we consider a bi‐quadratic homogeneous polynomial optimization problem over two unit spheres arising in nonlinear elastic material analysis and in entanglement studies in quantum physics. The problem is equivalent to computing the largest M‐eigenvalue of a fourth‐order tensor. To solve the problem, we propose a practical method whose validity is guaranteed theoretically. To make the sequence generated by the method converge to a good solution of the problem, we also develop an initialization scheme. The given numerical experiments show the effectiveness of the proposed method. Copyright © 2009 John Wiley & Sons, Ltd.     

Initialization in semidefinite programming via a self-dual skew-symmetric embedding

The formulation of interior point algorithms for semidefinite programming has become an active research area, following the success of the methods for large-scale linear programming. Many interior point methods for linear programming have now been extended to the more general semidefinite case, but the initialization problem remained unsolved.In this paper we show that the initialization strategy of embedding the problem in a self-dual skew-symmetric problem can also be extended to the semidefinite case. This method also provides a solution for the initialization of quadratic programs and it is applicable to more general convex problems with conic formulation.     

Performance of the MOSA Method for the Bicriteria Assignment Problem

The classical linear Assignment problem is considered with two objectives. The aim is to generate the set of efficient solutions. An exact method is first developed based on the two-phase approach. In the second phase a new upper bound is proposed so that larger instances can be solved exactly. The so-called MOSA (Multi-Objective Simulated Annealing) is then recalled; its efficiency is improved by initialization with a greedy approach. Its results are compared to those obtained with the exact method. Extensive numerical experiments have been realized to measure the performance of the MOSA method.     

Monotone iterations for numerical solutions of reaction-diffusion-convection equations with time delay

In this article we use the monotone method for the computation of numerical solutions of a nonlinear reaction-diffusion-convection problem with time delay. Three monotone iteration processes for a suitably formulated finite-difference system of the problem are presented. It is shown that the sequence of iteration from each of these iterative schemes converges from either above or below to a unique solution of the finite-difference system without any monotone condition on the nonlinear reaction function. An analytical comparison result among the three processes of iterations is given. Also given is the application of the iterative schemes to some model problems in population dynamics, including numerical results of a model problem with known analytical solution. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 339–351, 1998     

Data-driven optimal tracking control of switched linear systems

This paper addresses the optimal tracking control for switched linear systems with unknown dynamics. We convert the problem into an optimal control problem of the augmented switched systems. In view of the augmented systems, we propose a data-driven switched linear quadratic regular algorithm for obtaining the optimal switching signal under unknown system dynamics. It is proved that the optimal switching signal will not cause Zeno behavior and can make the system stable. Besides, with the proposed algorithm, we just need to identify an autonomous system instead of the original systems, which has fewer parameters to be determined. A numerical example is given to illustrate the validity of the main results.     

Feedback Integrators

A new method is proposed to numerically integrate a dynamical system on a manifold such that the trajectory stably remains on the manifold and preserves the first integrals of the system. The idea is that given an initial point in the manifold we extend the dynamics from the manifold to its ambient Euclidean space and then modify the dynamics outside the intersection of the manifold and the level sets of the first integrals containing the initial point such that the intersection becomes a unique local attractor of the resultant dynamics. While the modified dynamics theoretically produces the same trajectory as the original dynamics, it yields a numerical trajectory that stably remains on the manifold and preserves the first integrals. The big merit of our method is that the modified dynamics can be integrated with any ordinary numerical integrator such as Euler or Runge–Kutta. We illustrate this method by applying it to three famous problems: the free rigid body, the Kepler problem and a perturbed Kepler problem with rotational symmetry. We also carry out simulation studies to demonstrate the excellence of our method and make comparisons with the standard projection method, a splitting method and Störmer–Verlet schemes.     

Numerical Solution of the Inverse Nonequilibrium Sorption Problem with a Time-Dependent Boundary Condition

We consider the quasi-linear problem of nonequilibrium sorption dynamics with external-diffusion kinetics and a boundary condition that contains the time derivative of a solution component. A numerical method is proposed for describing the inverse problem to recover the nonlinear parameter of the system of differential equations—the inverse of the sorption isotherm. Convergence of the difference scheme for the direct problem is proved. Numerical solutions of both the direct and the inverse problem are obtained for various parameter values.     

Control of a fourth-order linear system with mixed constraints

A controlled fourth-order linear mechanical system, containing a vibrating member, is considered. Geometric constraints are imposed on the control and phase variables. The problem of bringing the system to a given state in a finite time is solved. The solution employs an approach based on Kalman's general scheme for constructing controls as linear combinations of characteristic motions of the uncontrolled system. Results of a numerical simulation of the dynamics of a closed system are presented     

BEM-based numerical study of three-dimensional compressible bubble dynamics in stokes flow

The dynamics of compressible gas bubbles in a viscous shear flow and an acoustic field at low Reynolds numbers is studied. The numerical approach is based on the boundary element method (BEM), which is effective as applied to the three-dimensional simulation of bubble deformation. However, the application of the conventional BEM to compressible bubble dynamics faces difficulties caused by the degeneration of the resulting algebraic system. Additional relations based on the Lorentz reciprocity principle are used to cope with this problem. Test computations of the dynamics of a single bubble and bubble clusters in acoustic fields and shear flows are presented.     

On the conservativeness of two-stage symmetric-symplectic Runge-Kutta methods and the Störmer-Verlet method

In the present paper, we discuss the problem on the total energy conservation for the numerical solution of the Cauchy problem for the equations of classical molecular dynamics by symplectic and symmetric methods. We consider the methods from a one-parameter family of two-stage symmetric-symplectic Runge-Kutta methods and the Störmer-Verlet method. In particular, we show that a numerical algorithm preserving the total energy of the system on the approximate solutions of the model Cauchy problem almost on the entire trajectory can be constructed on the basis of the one-parameter family of two-stage symmetric-symplectic Runge-Kutta methods.     

基于直方图初始化的有效的多区域Mumford-Shah图像分割方法（英文）

The Mumford-Shah energy functional is a successful image segmentation model. It is a non-convex variational problem and lacks of good initialization techniques so far. In this paper, motivated by the fact that image histogram is a combination of several Gaussian distributions, and their centers can be considered as approximations of cluster centers, we introduce a histogram-based initialization method to compute the cluster centers. With this technique, we then devise an effective multi-region Mumford-Shah image segmentation method, and adopt the recent proximal alternating minimization method to solve the minimization problem. Experiments indicate that our histogram initialization method is more robust than existing methods,and our segmentation method is very effective for both gray and color images.     

Analysis of the Structured Total Least Squares Problem for Hankel/Toeplitz Matrices

The Structured Total Least Squares (STLS) problem is a natural extension of the Total Least Squares (TLS) approach when structured matrices are involved and a similarly structured rank deficient approximation of that matrix is desired. In many of those cases the STLS approach yields a Maximum Likelihood (ML) estimate as opposed to, e.g., TLS.In this paper we analyze the STLS problem for Hankel matrices (the theory can be extended in a straightforward way to Toeplitz matrices, block Hankel and block Toeplitz matrices). Using a particular parametrisation of rank-deficient Hankel matrices, we show that this STLS problem suffers from multiple local minima, the properties of which depend on the parameters of the new parametrisation. The latter observation makes initial estimates an important issue in STLS problems and a new initialization method is proposed. The new initialization method is applied to a speech compression example and the results confirm the improved performance compared to other previously proposed initialization methods.     

On a New Diffeomorphic Multi-Modality Image Registration Model and Its Convergent Gauss-Newton Solver

In this work, we propose a new variational model for multi-modal image registration and present an efficient numerical implementation. The model minimizes a new functional based on using reformulated normalized gradients of the images as the fidelity term and higher-order derivatives as the regularizer. A key feature of the model is its ability of guaranteeing a diffeomorphic transformation which is achieved by a control term motivated by the quasi-conformal map and Beltrami coefficient. The existence of the solution of this model is established. To solve the model numerically, we design a Gauss-Newton method to solve the resulting discrete optimization problem and prove its convergence; a multilevel technique is employed to speed up the initialization and avoid likely local minima of the underlying functional. Finally, numerical experiments demonstrate that this new model can deliver good performances for multi-modal image registration and simultaneously generate an accurate diffeomorphic transformation.     

The numerical simulation of a typical fluid-structure interaction problem

To determine the dynamic response of a structure under the influence of the fluid flow one must solve a coupled computational fluid dynamics (CFD) and computational structural dynamics (CSD) mathematical problem. This paper presents the comparison of two methods for the calculation of the fluid-structure interaction. The first one is of explicit-implicit type and uses a staggered time advancement of the fluid and structure problems. The second uses a fully implicit discretization in the physical time of the fluid-structure equations and an explicit advancement in the dual-time. The physical fluid-structure problem is accompanied by the equations of the mesh motion, which are written as for a pseudo-structural system with its own dynamics. Representative numerical results are presented for the two degrees of freedom tipical section in unsteady transonic flow. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Asymptotic solutions of Fredholm integro-differential equations with rapidly varying kernels and irreversible limit operator

We consider a system of singularly perturbed integro-differential Fredholm equations with rapidly varying kernel in the case of irreversible operator of differential part. We develop an algorithm for constructing regularized asymptotic solutions. It is shown that in the presence of rapidly decreasing multiplier in the kernel the original problem is not on the spectrum (i.e, it is solvable for any right-hand side). We study the limit transition (with small parameter tending to zero), and solve the problem of initialization, i.e., the problem of extracting of the source data for which an exact solution to the system tends to the limit at all duration (including a zone of boundary layer).     

Numerical solution of optimization problems in vibrational mechanics using the characteristics method

Optimal control problems with a terminal pay-off functional are considered. The dynamics of the control system consists of rapid oscillatory and slow non-linear motions. A numerical method for solving these problems using the characteristics of the Hamilton–Jacobi–Bellman equation is presented. Estimates of the accuracy of the method are obtained. A theorem is proved which enables one to determine the class of functions containing the optimal preset control to be obtained. The results of the numerical solution of a terminal optimization problem for a fast non-linear pendulum are presented.     

Phase Retrieval with One or Two Diffraction Patterns by Alternating Projections with the Null Initialization

Alternating projection (AP) of various forms, including the parallel AP (PAP), real-constrained AP (RAP) and the serial AP (SAP), are proposed to solve phase retrieval with at most two coded diffraction patterns. The proofs of geometric convergence are given with sharp bounds on the rates of convergence in terms of a spectral gap condition. To compensate for the local nature of convergence, the null initialization is proposed to produce good-quality initial guess. Numerical experiments show that the null initialization is more accurate than the spectral initialization and that AP converges faster to the true object than other iterative schemes such as the Wirtinger flow (WF). In numerical experiments AP with the null initialization converges globally to the true object.