Optimal Convergence Rates Results for Linear Inverse Problems in Hilbert Spaces

In this article, we prove optimal convergence rates results for regularization methods for solving linear ill-posed operator equations in Hilbert spaces. The results generalizes existing convergence rates results on optimality to general source conditions, such as logarithmic source conditions. Moreover, we also provide optimality results under variational source conditions and show the connection to approximative source conditions.     

直接边界元法中边界积分的解析处理

确立了平面位势和弹性力学问题的边界元直接法中边界积分的解析计算框架系统,从而避免了传统的主似求积分,数值算例表胆它具有较高的精度和效率,特别是在边界量和边界附近区域内点物理量的计算可获得较高的精度。     

乘法半群上的线性方程组与晶体对势的反演

用乘法半群上的线性方程组来求解晶体原子间对势反演的逆问题.这种方法是解决此类问题的一般性方法.本文还给出了一个计算实例.     

Convergence Rates for Linear Inverse Problems in the Presence of an Additive Normal Noise

Abstract

In this work, we examine a finite-dimensional linear inverse problem where the measurements are disturbed by an additive normal noise. The problem is solved both in the frequentist and in the Bayesian frameworks. Convergence of the used methods when the noise tends to zero is studied in the Ky Fan metric. The obtained convergence rate results and parameter choice rules are of a similar structure for both approaches.     

平面弹性力学问题的混合元法的泡函数稳定性及其后验误差估计

该文讨论平面弹性力学问题的混合元法的泡函数稳定性,并导出基于简化的稳定化格式的一种先验误差估计和后验误差估计．这种简化的稳定化格式较通常的格式节省自由度．     

Banach空间中线性算子的Drazin广义逆的表示

1981年,1985年,2000年,乔三正、蔡东汉、魏益民分别给出Drazin广义逆不同形式的表示．本文将对上述结果进行推广,给出Drazin广义逆的统一表示,使上述三个结果均成为本文主要结果的特例．在本文主要结果的基础上,利用算子谱理论,给出Drazin广义逆的一种逼近形式的表示,同时给出逼近解的估计．此结果推广了蔡东汉、魏益民的相应结果．     

求解二维波动方程正演反演问题的半离散方法

本文用半离散方法将二维波动方程离散为一维耦合波动方程组．给出了离散的收敛性及波动方程组的适定性．利用这种方法可以求解波动方程系数及演问题．     

A New Energy Inversion for Parameter Identification in Saddle Point Problems with an Application to the Elasticity Imaging Inverse Problem of Predicting Tumor Location

The primary objective of this work is a detailed theoretical and computational study of the elasticity imaging inverse problem for tumor identification within the human body. Apart from this inverse problem's important and interesting application, it also poses noteworthy mathematical challenges since the underlying mathematical model is a system of elasticity involving incompressibility. This gives rise to the “locking” effect and special treatment is necessary for both the direct and inverse problems. To study the inverse problem in an optimization framework, we introduce a general computational scheme for handling parameter identification in saddle point problems along with the introduction and analysis of a new energy output least-squares objective functionals. We also present a treatment of the identification of discontinuous elasticity coefficients using the total variation regularization method. General formulas for the computation of the coefficient-to-solution map and a complete convergence analysis are given for the continuous problem as well as for its discrete analogue. Discrete formulas and implementation issues are discussed in detail and numerical examples for smooth and discontinuous coefficients are given.     

On the Finite Termination of an Entropy Function Based Non-Interior Continuation Method for Vertical Linear Complementarity Problems

By using a smooth entropy function to approximate the non-smooth max-type function, a vertical linear complementarity problem (VLCP) can be treated as a family of parameterized smooth equations. A Newton-type method with a testing procedure is proposed to solve such a system. We show that under some milder than usual assumptions the proposed algorithm finds an exact solution of VLCP in a finite number of iterations. Some computational results are included to illustrate the potential of this approach. This author’s work was partially supported by the National Natural Science Foundation of China (Grant Nos. 10271002 and 10401038). This author’s work was partially supported by the Scientific Research Foundation of Tianjin University for the Returned Overseas Chinese Scholars and the Scientific Research Foundation of Liu Hui Center for Applied Mathematics, Nankai University-Tianjin University.     

Asymptotic Distribution of Eigenfunctions and Eigenvalues of the Basic Boundary-Contact Oscillation Problems of the Classical Theory of Elasticity

The basic boundary-contact oscillation problems are considered for a three-dimensional piecewise-homogeneous isotropic elastic medium bounded by several closed surfaces. Using Carleman's method, the asymptotic formulas for the distribution of eigenfunctions and eigenvalues are obtained.     

黄金分割法在求解线性方程组中的应用

把线性方程组转化成为其等价的变分问题 ,借助黄金分割法的思想对该变分问题进行求解 ,给出的算例结果表明 ,本文方法不仅对良态线性方程组的求解有效 ,而且对于病态线性方程组的求解同样是有效实用的 .     

Avoiding Order Reduction of Runge–Kutta Discretizations for Linear Time-Dependent Parabolic Problems

A technique is developed in this paper to avoid order reduction when discretizing linear parabolic problems with time dependent operator using Runge–Kutta methods in time and standard schemes in space. In an abstract framework, the boundaries of the stages of the Runge–Kutta method which would completely avoid the order reduction are given. Then, the possible practical implementations for the calculus of those boundaries from the given data are studied, and the full discretization is completely analyzed. Some numerical experiments are included. This revised version was published online in July 2006 with corrections to the Cover Date.     

电磁散射问题手性障碍边界重构的线性探测法的理论分析

In this paper,we consider the inverse scattering by chiral obstacle inelectromagnetic fields,and prove that the linear sampling method is also effective todetermine the support of a chiral obstacle from the noisy far field data.     

基于三维弹性问题的多极边界元法核函数分解

多极边界元法已经成功地应用于大规模工程计算中.得到并且证明了基于三维弹性问题的多极边界元法核函数分解的定理(定理1),完善了多击边界元法的数学理论.     

车灯线光源的离散法优化设计

采用连续问题离散化的方法 ,将线光源离散成由 n =l△ l 个点光源组成 ,反射面离散为m =S△ S 个小平面的组合 ,利用空间光线在抛物面上的反射规律 ,讨论了光强度与线光源的长度的关系 ,得出了符合要求的线光源的长度的最值及设计方案 ,该方法计算方便 ,便于实际操作且具有很大的推广性 .     

Equivalence of three kinds of well-posed-ness of discontinuous Riemann-Hilbert problem for elliptic complex equation in multiply connected domains

In this article, we first introduce the general linear elliptic complex equation of first order with certain conditions, and then propose discontinuous Riemann-Hilbert problem and some kinds of modified well-posed-ness for the complex equation. Then we verify the equivalence of three kinds of well-posed-ness. The discontinuous boundary value problem possesses many applications in mechanics and physics etc.     

Numerical Solution of Partial Differential Equations in Arbitrary Shaped Domains Using Cartesian Cut-Stencil Finite Difference Method. Part I: Concepts and Fundamentals

A new finite difference (FD) method, referred to as "Cartesian cut-stencil FD", is introduced to obtain the numerical solution of partial differential equations on any arbitrary irregular shaped domain. The 2^nd-order accurate two-dimensional Cartesian cut-stencil FD method utilizes a 5-point stencil and relies on the construction of a unique mapping of each physical stencil, rather than a cell, in any arbitrary domain to a generic uniform computational stencil. The treatment of boundary conditions and quantification of the solution accuracy using the local truncation error are discussed. Numerical solutions of the steady convection-diffusion equation on sample complex domains have been obtained and the results have been compared to exact solutions for manufactured partial differential equations (PDEs) and other numerical solutions.     

最小平方距离法和隐式线性函数关系的参数估计

本文再次指出:最小平方距离法(LSD),也可叫最小模方法,是从解决多维空间n个点的超平面拟合问题而提取的;通过对p个随机变量的n组观测值,此方法是探求它们之间是否存在隐式线形函数关系的好方法;使用此法,可以得出隐式线形函数关系较好的参数估计。     

The Method of Rothe and Two-Scale Convergence in Nonlinear Problems

Modelling of macroscopic behaviour of materials, consisting of several layers or components, cannot avoid their microstructural properties. This article demonstrates how the method of Rothe, described in the book of K. Rektorys The Method of Discretization in Time, together with the two-scale homogenization technique can be applied to the existence and convergence analysis of some strongly nonlinear time-dependent problems of this type.