An inverse problem for stochastic differential equations

We discuss the problem of reconstructing the drift coefficient of a diffusion from the knowledge of the transition probabilities outside a given bounded region in ^d ,d>1. We also give an interpretation of the solution of this inverse problem in the framework of stochastic mechanics.This paper is dedicated to the dear memory of Paola Calderoni.     

Application of multi-phase particle swarm optimization technique to inverse radiation problem

The multi-phase particle swarm optimization (MPPSO) technique is applied to the inverse radiation problem in the present paper. The directional radiative intensities are served as the measurement data to estimate the radiative source term, optical thickness, scattering albedo, and phase function in one-dimensional semitransparent plane-parallel media by the inverse simulation. To check the performance and accuracy in retrieval, a comparison is presented between three PSO methods, i.e. the MPPSO, the standard PSO, and the Stochastic PSO. The results confirm the potential of the proposed approach MPPSO and show its effectiveness and superiority over the other two PSO algorithms. Furthermore, the effects of swarm size, searching space, phase change frequency, and velocity-reinitializing frequency on the convergence velocity and computational accuracy of MPPSO are also investigated.     

利用辐射场诊断模式成分的病态特性

 根据辐射场测量并确定已知尺寸和工作频率的波导内可能存在的模式(即模式诊断)是一个电磁逆问题。针对该逆问题,利用电磁分析软件计算过模圆波导馈电圆锥喇叭天线的远区辐射场,通过比较多种模式以不同比例组合馈电产生的总辐射场之间的差异,研究模式诊断逆问题的病态程度。结果表明该逆问题病态严重,模式成分比例变化很大,但某些位置上辐射场变化较小甚至没有区别,即辐射场对各模式比例变化不敏感。如果在这些位置上测量则很难判断出各模式的成分,只有通过增加信息量才能缓解病态的严重程度,得到较准确的结果。通过分析,给出了对各模式组合比例变化敏感的测量位置。     

Manipulation with heterogeneity within a species population formulated as an inverse problem

Dependence of the evolutionary dynamics on the population’s heterogeneity has been reliably recognized and studied within the frame of evolutionary optimization theory. As the causal relation between the heterogeneity and dynamics of environment has been revealed, the possibility to influence convergence rate of evolutionary processes by purposeful manipulation with environment emerges.For the above purposes we formulate the task as the inverse problem meaning that desired population heterogeneity, quantified by Tsallis information entropy, represents the model’s input and dynamics of environment leading to desired population heterogeneity is looked for. Here the presented abstract model of evolutionary motion within the inverse model of replicating species is case-independent and it is relevant for the broad range of phenomena observed at cellular, ecological, economic and social scales. We envision relevance of the model for anticancer therapy, in which the effort is to circumvent heterogeneity as it typically correlates with the therapy efficiency.     

An accurate solver for forward and inverse transport

This paper presents a robust and accurate way to solve steady-state linear transport (radiative transfer) equations numerically. Our main objective is to address the inverse transport problem, in which the optical parameters of a domain of interest are reconstructed from measurements performed at the domain’s boundary. This inverse problem has important applications in medical and geophysical imaging, and more generally in any field involving high frequency waves or particles propagating in scattering environments. Stable solutions of the inverse transport problem require that the singularities of the measurement operator, which maps the optical parameters to the available measurements, be captured with sufficient accuracy. This in turn requires that the free propagation of particles be calculated with care, which is a difficult problem on a Cartesian grid.     

Relevance vector machine technique for the inverse scattering problem

A novel method based on the relevance vector machine(RVM) for the inverse scattering problem is presented in this paper.The nonlinearity and the ill-posedness inherent in this problem are simultaneously considered.The nonlinearity is embodied in the relation between the scattered field and the target property,which can be obtained through the RVM training process.Besides,rather than utilizing regularization,the ill-posed nature of the inversion is naturally accounted for because the RVM can produce a probabilistic output.Simulation results reveal that the proposed RVM-based approach can provide comparative performances in terms of accuracy,convergence,robustness,generalization,and improved performance in terms of sparse property in comparison with the support vector machine(SVM) based approach.     

An inverse algorithm for reconstructing the acoustic parameters of layered medium

I.IntroductionThcinvcrsescattcringprob1cmisbasicinarcassuchasradar,sonar,geophysica1explora-honandnondestructivctesting.Whenthescatteringobjectisaconstitutiveparametersvaryinacontinuousmannerandtheinvcrseproblcmistodctcrmineoneormoreoftheseparametcrsfromtheobserveddata.Muchofthcpreviousworkonthisproblcmforthecaseofp1anewaveatnormalincidencchasconsistcdofderivingaSchr6dingerequationfromthebasicacousticandstrcss-straincquations,andthenreconstruchngthepotentialappearinginthisequationbyusingtheGe…     

An inverse radiative transfer problem of simultaneously estimating profiles of temperature and radiative parameters from boundary intensity and temperature measurements

A parallel-plane space filled with absorbing, emitting, isotropically scattering, gray medium is studied in this paper. The boundary intensity and boundary temperature profiles are calculated for the inverse analysis. For the simultaneous estimation of temperature, absorption and scattering coefficient profiles in the medium, the sum of residuals of boundary intensity and temperature after being weighted by a balance factor is minimized through using a Newton-type iteration algorithm and the least-squares method. To avoid over-updating for the parameters, the relative updating magnitude during the iteration process is constrained not to be >0.5. It is shown that the boundary intensity measurement alone is not enough to estimate simultaneously the temperature (source) and the radiative properties (both absorption and scattering coefficients) when the measurement data contain sensitive random errors. The boundary temperature measurement can serve as a necessary supplementation to the boundary intensity to make this kind of inverse radiative transfer problem resolvable. It was shown that a compensation relationship between absorption and scattering coefficients makes it difficult to fix them accurately. Parabolic profiles for the three parameters are used to validate the estimation method. When the optical thickness approaches 4.0, the results for the radiative properties are not acceptable, although the result for temperature profile is reasonable. This means the method needs further improvements.     

广义Birkhoff系统动力学的一类逆问题

针对广义Birkhoff系统动力学,提出广义Birkhoff系统动力学的一类逆问题,研究由已知积分流形来建立广义Birkhoff方程. 这类逆问题的解通常不是唯一的,需给出必要的补充要求. 最后举例说明结果的应用. 关键词： 广义Birkhoff系统 动力学逆问题 积分流形     

An in-process method for the inverse estimation of the powder catchment efficiency during laser cladding

An estimation of the heat loss by conduction can be obtained from measurements of the surface temperature and an overall heat balance at the clad laser interaction zone. Through an inverse calculation of the boundary temperature from observed surface temperatures the powder catchment efficiency can be estimated along with the variation in the clad height expected during laser cladding. This method shows a possible way to monitor and control the clad height and profile as required by near net shape manufacturing methods based on laser cladding.     

The Method of Fundamental Solutions for Solving Exterior Axisymmetric Helmholtz Problems with High Wave-Number

In this paper, we investigate the method of fundamental solutions (MFS) for solving exterior Helmholtz problems with high wave-number in axisymmetric domains. Since the coefficient matrix in the linear system resulting from the MFS approximation has a block circulant structure, it can be solved by the matrix decomposition algorithm and fast Fourier transform for the fast computation of large-scale problems and meanwhile saving computer memory space. Several numerical examples are provided to demonstrate its applicability and efficacy in two and three dimensional domains.     

A Moving Pseudo-Boundary MFS for Three-Dimensional Void Detection

We propose a new moving pseudo-boundary method of fundamental solutions (MFS) for the determination of the boundary of a three-dimensional void (rigid inclusion or cavity) within a conducting homogeneous host medium from overdetermined Cauchy data on the accessible exterior boundary. The algorithm for imaging the interior of the medium also makes use of radial spherical parametrization of the unknown star-shaped void and its centre in three dimensions. We also include the contraction and dilation factors in selecting the fictitious surfaces where the MFS sources are to be positioned in the set of unknowns in the resulting regularized nonlinear least-squares minimization. The feasibility of this new method is illustrated in several numerical examples.     

一种基于虚拟传感的无需误差传声器的自适应有源降噪方法*

该文基于虚拟传感技术引入了一种用于耳机的无需误差传声器的自适应有源降噪方法。该算法仅使用一个参考传声器实现了一种前馈和反馈自适应算法结合的有源降噪算法,提高了有源降噪稳定性,简化了耳机硬件结构。利用DSP平台实现了该文提出的方案,并通过实验验证了其良好的降噪性能和实用价值。     

A Practical Algorithm for Determining the Optimal Pseudo-Boundary in the Method of Fundamental Solutions

One of the main difficulties in the application of the method of fundamental solutions (MFS) is the determination of the position of the pseudo-boundary on which are placed the singularities in terms of which the approximation is expressed. In this work, we propose a simple practical algorithm for determining an estimate of the pseudo-boundary which yields the most accurate MFS approximation when the method is applied to certain boundary value problems. Several numerical examples are provided.     

On the use of the resting potential and level set methods for identifying ischemic heart disease: An inverse problem

The electrical activity in the heart is modeled by a complex, nonlinear, fully coupled system of differential equations. Several scientists have studied how this model, referred to as the bidomain model, can be modified to incorporate the effect of heart infarctions on simulated ECG (electrocardiogram) recordings.We are concerned with the associated inverse problem; how can we use ECG recordings and mathematical models to identify the position, size and shape of heart infarctions? Due to the extreme CPU efforts needed to solve the bidomain equations, this model, in its full complexity, is not well-suited for this kind of problems. In this paper we show how biological knowledge about the resting potential in the heart and level set techniques can be combined to derive a suitable stationary model, expressed in terms of an elliptic PDE, for such applications. This approach leads to a nonlinear ill-posed minimization problem, which we propose to regularize and solve with a simple iterative scheme.Finally, our theoretical findings are illuminated through a series of computer simulations for an experimental setup involving a realistic heart in torso geometry. More specifically, experiments with synthetic ECG recordings, produced by solving the bidomain model, indicate that our method manages to identify the physical characteristics of the ischemic region(s) in the heart. Furthermore, the ill-posed nature of this inverse problem is explored, i.e. several quantitative issues of our scheme are explored.     

对称分类在非线性偏微分方程组边值问题中的应用

研究了微分方程对称分类在非线性偏微分方程组边值问题中的应用.首先,利用偏微分方程(组)完全对称分类微分特征列集算法确定了给定非线性偏微分方程组边值问题的完全对称分类;其次,利用一个扩充对称将非线性偏微分方程组边值问题约化为常微分方程组初值问题;最后,利用龙格-库塔法求解了常微分方程组初值问题的数值解.     

Mathematical and numerical studies on meshless methods for exterior unbounded domain problems

The method of fundamental solutions (MFS) is an efficient meshless method for solving boundary value problems in an exterior unbounded domain. The numerical solution obtained by the MFS is accurate, while the corresponding matrix equation is ill-conditioned. A modified MFS (MMFS) with the proper basis functions is proposed by the introduction of the modified Trefftz method (MTM). The concrete expressions of the corresponding condition numbers are given in mathematical forms and the solvability by these methods is mathematically proven. Thereby, the optimal parameter minimizing the condition number is also mathematically given. Numerical experiments show that the condition numbers of the matrices corresponding to the MTM and the MMFS are reduced and that the numerical solution by the MMFS is more accurate than the one by the conventional method.     

A Meshless Regularization Method for a Two-Dimensional Two-Phase Linear Inverse Stefan Problem

In this paper, a meshless regularization method of fundamental solutions is proposed for a two-dimensional, two-phase linear inverse Stefan problem. The numerical implementation and analysis are challenging since one needs to handle composite materials in higher dimensions. Furthermore, the inverse Stefan problem is ill-posed since small errors in the input data cause large errors in the desired output solution. Therefore, regularization is necessary in order to obtain a stable solution. Numerical results for several benchmark test examples are presented and discussed.     

An exact solution of the moving boundary problem for the expansion of a plasma cylinder in a magnetic field

An exact analytic solution has been obtained for a uniformly expanding, neutral, infinitely conducting plasma cylinder in an external uniform and constant magnetic field. The electrodynamical aspects related to the emission and transformation of energy have been considered as well. The results obtained can be used in analysing the recent experimental and simulation data.     

Residual stresses deriving from holographic interferometry data on a base of inverse problem solution

A transition model, which is capable of obtaining both membrane and bending residual stress components from initial experimental information, is developed for thin-walled plane structures. The determination of residual stresses is based on the combined implementing of the hole-drilling method and reflection hologram interferometry. Required input data are obtained by simultaneous measurements on through hole distortions in two principal strain directions on opposite sides of thin plane specimen. These sides are faces of the drill entrance and exit. Superimposed residual stresses field, which consists of both membrane and bending components, is a reason for the various deviations of each specific fringe pattern from an ideal form. This fact is a clear experimental indication of the bending stress contribution in a total stress field. Two ways of decomposition of superimposed residual stresses field are proposed and analysed in detail. Emphasis is laid on a careful quantitative formulation of the inverse problem needed for an accurate deriving both membrane and bending residual stress components. It is shown that an availability of two-side initial data is both an essential and necessary condition of such a formulation. Detailed analysis of an accuracy of the results obtained is performed. This analysis is based on a wide set of both actual interferograms and analogous reference fringe patterns related to superimposed residual stress field under study. Comparing residual stress values obtained proceeding from one-side and two-side data are presented for different types of superimposed field of interest.