Source term identification in 1-D IHCP

We propose stable numerical solutions for the simultaneous identification of temperature, temperature gradient, and general source terms in the one-dimensional inverse heat conduction problem (IHCP).

The numerical solution consists of a regularization procedure, based on the mollification method,and a marching scheme for the solution of the stabilized problem. The stability, error analysis and implementation of the algorithm are presented together with a set of numerical results.     

Mathematical programs with complementarity constraints in the context of inverse optimal control for locomotion

In this paper an inverse optimal control problem in the form of a mathematical program with complementarity constraints (MPCC) is considered and numerical experiences are discussed. The inverse optimal control problem arises in the context of human navigation where the body is modelled as a dynamical system and it is assumed that the motions are optimally controlled with respect to an unknown cost function. The goal of the inversion is now to find a cost function within a given parametrized family of candidate cost functions such that the corresponding optimal motion minimizes the deviation from given data. MPCCs are known to be a challenging class of optimization problems typically violating all standard constraint qualifications (CQs). We show that under certain assumptions the resulting MPCC fulfills CQs for MPCCs being the basis for theory on MPCC optimality conditions and consequently for numerical solution techniques. Finally, numerical results are presented for the discretized inverse optimal control problem of locomotion using different solution techniques based on relaxation and lifting.     

Efficient solution of nonlinear,underdetermined inverse problems with a generalized PDE model

Several parameter estimation problems (or “inverse” problems) such as those that occur in hydrology and geophysics are solved using partial differential equation (PDE)-based models of the physical system in question. Likewise, these problems are usually underdetermined due to the lack of enough data to constrain a unique solution. In this paper, we present a framework for the solution of underdetermined inverse problems using COMSOL Multiphysics (formerly FEMLAB) that is applicable to a broad range of physical systems governed by PDEs. We present a general adjoint state formulation which may be used in this framework and allows for faster calculation of sensitivity matrices in a variety of commonly encountered underdetermined problems. The aim of this approach is to provide a platform for the solution of inverse problems that is efficient, flexible, and not restricted to one particular scientific application.We present an example application of this framework on a synthetic underdetermined inverse problem in aquifer characterization, and present numerical results on the accuracy and efficiency of this method. Our results indicate that our COMSOL-based routines provide an accurate, flexible, and scalable method for the solution of PDE-based inverse problems.     

Iterative solvers for Tikhonov regularization of dense inverse problems

According to the special demands arising from the development of science and technology, in the last decades appeared a special class of problems that are inverse to the classical direct ones. Such an inverse problem is concerned with the opposite way, usually followed by a direct one: finding the cause of a given effect or finding the law of evolution given the cause and effect. Very frequently, such inverse problems are modelled by Fredholm first-kind integral equations that give rise after discretization to (very) ill-conditioned linear systems, in classical or least squares formulation. Then, an efficient numerical solution can be obtained by using the Tikhonov regularization technique. In this respect, in the present paper, we propose three Kovarik-like algorithms for numerical solution of the regularized problem. We prove convergence for all three methods and present numerical experiments on a mathematical model of an inverse problem concerned with the determination of charge distribution generating a given electric field.     

全变差正则化在抛物型方程初始条件反问题的应用

当反问题反演的函数不连续时,一般的正则化算法反演效果不令人满意,用全变差正则化方法对抛物型方程初始条件反问题进行求解,并进行了数值分析和数值模拟,结果显示数值解与真解吻合较好,表明该方法对于不连续函数求解具有高效、稳定等优点.     

Solvability of the Nonlocal Inverse Parabolic Problem and Numerical Results

In this paper, we consider the unique solvability of the inverse problem of determining the right-hand side of a parabolic equation whose leading coefficient depends on time variable under nonlocal integral overdetermination condition. We obtain sufficient conditions for the unique solvability of the inverse problem. The existence and uniqueness of the solution of the inverse parabolic problem upon the data are established using the fixed point theorem. This inverse problem appears extensively in the modelling of various phenomena in engineering and physics. For example, seismology, medicine, fusion welding, continuous casting, metallurgy, aircraft, oil and gas production during drilling and operation of wells. In addition, the numerical solution of the inverse problem is studied by using the Crank-Nicolson finite difference method together with the Tikhonov regularization to find a stable and accurate approximate solution of finite differences. The resulting nonlinear system of parabolic equation is solved computationally using the MATLAB subroutine lsqnonlin. Both analytical and numerically simulated noisy input data are inverted. The root mean square error values for various noise levels for both continuous and discontinuous time-dependent heat source term are compared. Numerical results presented for two examples show the efficiency of the computational method and the accuracy and stability of the numerical solution even in the presence of noise in the input data. Furthermore, the choice of the regularization parameter is also discussed based on the trial and error technique.     

A solution to the inverse kinematics of redundant manipulators

A solution to the inverse kinematics is a set of joint coordinates which correspond to a given set of task space coordinates (position and orientation of end effector). For the class of kinematically redundant robots, the solution is generically nonunique such that special methods are required for obtaining a solution. The method addressed in the paper, introduced earlier and termed “generalized inverse,” is based on a certain partitioning of the Jacobian functional corresponding to a nonlinear relationship of the inverse kinematics type. The article presents a new algorithm for solving the inverse kinematics using the method of generalized inverse based on a modified Newton-Raphson iterative technique. The new algorithm is efficient, converges rapidly, and completely generalizes the solution of the inverse kinematics problem for redundant robots. The method is illustrated by numerical examples.     

Fast and stable solution¶of banded-plus-semiseparable linear systems

We describe a stable algorithm, having linear complexity, for the solution of banded-plus-semiseparable linear systems. The algorithm exploits the structural properties of the inverse of a semiseparable matrix. Stability is achieved by combining these properties with partial pivoting techniques. Several numerical experiments are shown to confirm the effectiveness of the proposed approach. Received: April 2002 / Accepted: September 2002     

谱约束下反对称正交反对称矩阵束的最佳逼近

东莞理工学院软件学院,广东东莞523808矩阵逆特征值问题的研究已有一些很好的成果,但对矩阵广义逆特征值问题的研究还不多见,本文讨论了反对称正交反对称矩阵的广义逆特征值问题,丰富了矩阵理论和方法,得到了通解和最佳逼近解,并给出了数值例子.     

A perturbation approach for a type of inverse linear programming problems

We consider an inverse linear programming (LP) problem in which the parameters in both the objective function and the constraint set of a given LP problem need to be adjusted as little as possible so that a known feasible solution becomes the optimal one. We formulate this problem as a linear complementarity constrained minimization problem. With the help of the smoothed Fischer–Burmeister function, we propose a perturbation approach to solve the inverse problem and demonstrate its global convergence. An inexact Newton method is constructed to solve the perturbed problem and numerical results are reported to show the effectiveness of the approach.     

Mathematical models of elastic wave processes in seismology and seismic prospecting: forward and inverse problems

This paper presents a simulation system based on the solution of forward and inverse problems of elastic wave propagation. Forward and inverse modeling have become a useful tool for interpretation in exploration geophysics and seismology. By 3D modeling, field observations can be simulated numerically, and computed results can be compared to field data. It is known that 3D seismic modeling requires the up-to-date high-performance multiprocessor computer systems which are not readily available for many geophysical firms in Russia. This situation makes us focus on the creation of efficient numerical–analytical algorithms which allow the solution of 3D forward and inverse seismic problems without supercomputers. This paper presents some algorithms and numerical experiments for different models of media. A special emphasis is given to “non-ray” waves which play an important role in seismic interpretation theory.     

Hamilton矩阵反问题的最小二乘解

本文讨论了Hamilton矩阵反问题最小二乘解,得到解的通式．考虑了解集合对给定矩 阵的最佳逼近问题,给出了唯一最佳逼近解的表达式．最后,我们给出了相应的数值算法及 数值实例．     

Existence and uniqueness of solutions of the inverse dynamics of multilink flexible arms: Convergence of a numerical scheme

The inverse dynamics of flexible manipulators consists in finding the joint torques that drive the end effector of a manipulator along a given trajectory. A proof of the existence and uniqueness of solutions of the inverse dynamics of planar open-chain multilink flexible manipulators and the convergence of an algorithm previously proposed by the second author for the numerical solution of the problem are presented in this article. Higher-order strain terms leading to the geometric stiffness matrix in finite element models are included in the formulation and their influence on the inverse dynamics solution is studied.     

An inverse problem for estimation of bending stiffness in Kirchhoff–Love plates

This work deals with the development of a numerical method for solving an inverse problem for bending stiffness estimation in a Kirchhoff–Love plate from overdetermined data. The coefficient is identified using a technique called the Method of Variational Imbedding, where the original inverse problem is replaced by a minimization problem. The Euler–Lagrange equations for minimization comprise higher-order equations for the solution of the displacement and an equation for the bending stiffness. The correctness of the embedded problem is discussed. A difference scheme and a numerical algorithm for solving the parameter identification problem are developed. Numerical results for the obtained values of the bending stiffness as an inverse problem are presented.     

由两个特征对构造正定Jacobi矩阵

132 数值计算与计算机应用2002 it巨1.引言具有如下形状的实对称矩阵:ldl hi \Ibid。。6,l人=I”.”·.’·.I 门)16、_,tL、_fo、_11\ b、;_IG、l称为n阶实对称三对角矩阵.若人还满足:(a)b;>     

Jacobi矩阵逆特征问题存在唯一解的条件

（?）1.引 言设n阶Jacobi矩阵为记Jp,q为Jn的主子矩阵;即 关于Jacobi矩阵逆特征值问题的研究文献很多,类型有由两组谱数据或两个特征对（指特征值及相应的特征向量）构造Jacobi矩阵的元素[1].由主子阵及一组谱数据构造Jacobi     

三对角对称矩阵逆特征问题存在唯一解的条件

n阶实对称矩阵称为n阶三对角对称矩阵,全体记为S0,若：i）b;＞0（i二1,2,…,l－1）称T为J。C0bi矩阵,全体记为SI;11）h＜0（d一1;又…,n－1）称T为负JacoN矩阵,全体记为h;iii）b;一叩一1,2,…,n－1）称T为不可约三对角对称矩阵,全体记为龙Jacobi矩阵的逆特征问题有广泛的应用;近年来有了较大发展［‘,’］．对由二个特征对构造相应的Jacobi矩阵或三对角对称矩阵问题的研究相对地比较成熟I’一义而对由三个特征对构造相应的Jacobi矩阵或三对角对称矩阵问题的研究却进展迟缓．文门对此作了一些尝试,本文具体研究：…     

Efficient closed-form solution of inverse kinematics for a specific six-DOF arm

Inverse kinematics solutions for multi-DOF arms can be classified as analytical or numerical. In general, analytical solutions are preferable to numerical solutions because analytical ones yield complete solutions and are computationally fast and reliable. However, analytical closed-form solutions for inverse kinematics of 6-DOF arms rarely exist for real-time control purposes of fast moving arms. In this paper, we propose a fast inverse kinematics algorithm with a closed-form solution for a specific 6-DOF arm. The proposed algorithm is verified using simulation modules developed by us for demonstrations.     

Numerical procedure for the determination of an unknown parameter in a parabolic equation

A numerical procedure for an inverse problem of determining unknown source parameter of one-dimensional parabolic equation subject to the specification of the solution at internal point along with the usual initial boundary conditions is considered. By using some transformation the problem is reformulated to a nonlocal parabolic problem. Some numerical examples using the proposed numerical procedure are presented.     

A Multigrid Approach for Minimizing a Nonlinear Functional for Digital Image Matching

In this paper, we consider a multigrid application in digital image processing. Here, the problem is to find a map, which transforms an image T into another image R such that the grey level of the different images are nearly equal in every picture-element. This problem arises in the investigation of human brains. The complete inverse problem is ill posed in the sense of Hadamard and nonlinear, so the numerical solution is quite difficult. We solve the inverse problem by a Landweber iteration. In each minimization step an approximate solution for the linearized problem is computed with a multigrid method as an inner iteration. Finally, we present some experimental results for synthetic and real images. Received December 30, 1998; revised August 16, 1999