Spectral regularization method for solving a time-fractional inverse diffusion problem

In this paper, we consider an inverse problem for a time-fractional diffusion equation with one-dimensional semi-infinite domain. The temperature and heat flux are sought from a measured temperature history at a fixed location inside the body. We show that such problem is severely ill-posed and further apply a spectral regularization method to solve it based on the solution given by the Fourier method. Convergence estimates are presented under a priori bound assumptions for the exact solution. Finally, numerical examples are given to show that the proposed numerical method is effective.     

Spectral regularization method for a Cauchy problem of the time fractional advection-dispersion equation

In this paper, a Cauchy problem for the time fractional advection-dispersion equation (TFADE) is investigated. Such a problem is obtained from the classical advection-dispersion equation by replacing the first-order time derivative by the Caputo fractional derivative of order . We show that the Cauchy problem of TFADE is severely ill-posed and further apply a spectral regularization method to solve it based on the solution given by the Fourier method. The convergence estimate is obtained under a priori bound assumptions for the exact solution. Numerical examples are given to show the effectiveness of the proposed numerical method.     

Heat Transfer in a Region with Non-Uniform Boundary Conditions

Heat transfer in a rectangular region with non-uniform conditions on the walls is considered. The temperature is given on both vertical walls and a part of the upper wall. The remainder of the upper wall and the lower horizontal wall are perfectly insulated. This boundary value problem is reduced to dual Fourier series equations. Those equations are simplified under the assumption that the height of the region is greater than the length or comparable to it. An exact solution of the simplified equations is constructed by using the Schwinger transformation, which has been used successfully in analyzing the electro-dynamics of wave guides. Numerical solutions also are found using a commercial finite element solver and a finite difference solver written in FORTRAN. Results for the average temperature and the temperature distribution in the region for a variety of high temperature boundary locations are in very good agreement among the three solution techniques.     

An improved non-local boundary value problem method for a cauchy problem of the Laplace equation

In this paper, we propose an improved non-local boundary value problem method to solve a Cauchy problem for the Laplace equation. It is known that the Cauchy problem for the Laplace equation is severely ill-posed, i.e., the solution does not depend continuously on the given Cauchy data. Convergence estimates for the regularized solutions are obtained under a-priori bound assumptions for the exact solution. Some numerical results are given to show the effectiveness of the proposed method.     

A modified method for a non-standard inverse heat conduction problem

A non-standard inverse heat conduction problem is considered. Data are given along the line x = 1 and the solution at x = 0 is sought. The problem is ill-posed in the sense that the solution (if it exists) does not depend continuously on the data. In order to solve the problem numerically it is necessary to employ some regularization method. In this paper, we study a modification of the equation, where a fourth-order mixed derivative term is added. Error estimates for this equation are given, which show that the solution of the modified equation is an approximation of the heat equation. A numerical implementation is considered and a simple example is given. Some numerical results show the usefulness of the modified method.     

Modified regularization method for the Cauchy problem of the Helmholtz equation

In this paper, the Cauchy problem for the Helmholtz equation is investigated. It is known that such problem is severely ill-posed. We propose a modified regularization method to solve it based on the solution given by the method of separation of variables. Convergence estimates are presented under two different a-priori bounded assumptions for the exact solution. Finally, numerical examples are given to show the effectiveness of the proposed numerical method.     

An aggregate deformation homotopy method for min-max-min problems with max-min constraints

In this paper, the constrained min-max-min problem, which is an essentially nonsmooth and nonconvex problem, is considered. Based on a twice aggregate function with a modification, an aggregate deformation homotopy method is established. Under some suitable assumptions, a smooth path from a randomly given point to a solution of the generalized KKT system is proven to exist. By numerically tracing the smooth path, a globally convergent algorithm for some solution of the problem is given. Some numerical results are given to show the feasibility of the method.     

A computational method for large‐scale differential symmetric Stein equation

We propose a numerical method for solving large‐scale differential symmetric Stein equations having low‐rank right constant term. Our approach is based on projection the given problem onto a Krylov subspace then solving the low dimensional matrix problem by using an integration method, and the original problem solution is built by using obtained low‐rank approximate solution. Using the extended block Arnoldi process and backward differentiation formula (BDF), we give statements of the approximate solution and corresponding residual. Some numerical results are given to show the efficiency of the proposed method.     

Numerical solution of an inverse initial boundary-value problem for the full time-dependent Maxwell's equations in the presence of imperfections of small volume

We consider the numerical solution, in a three-dimensional bounded domain, of the inverse problem for identifying the location of small electromagnetic imperfections in a medium with homogeneous background. Our numerical algorithm is based on the coupling of a discontinuous Galerkin method for the time-dependent Maxwell's equations, on the exact controllability method and on a Fourier inversion. Several numerical results are given with one and two imperfections and the robustness and accuracy of the numerical method used for the dynamic detection problem are shown.     

Recovering the Initial Distribution for a Time-Fractional Diffusion Equation

We consider the inverse problem of reconstructing the initial condition of a one-dimensional time-fractional diffusion equation from measurements collected at a single interior location over a finite time-interval. The method relies on the eigenfunction expansion of the forward solution in conjunction with a Tikhonov regularization scheme to control the instability inherent in the problem. We show that the inverse problem has a unique solution provided exact data is given, and prove stability results regarding the regularized solution. Numerical realization of the method and illustrations using a finite-element discretization are given at the end of this paper.     

A new regularization method for a Cauchy problem of the time fractional diffusion equation

In this paper, we consider a Cauchy problem of the time fractional diffusion equation (TFDE). Such problem is obtained from the classical diffusion equation by replacing the first-order time derivative by the Caputo fractional derivative of order α (0 < α ≤ 1). We show that the Cauchy problem of TFDE is severely ill-posed and further apply a new regularization method to solve it based on the solution given by the Fourier method. Convergence estimates in the interior and on the boundary of solution domain are obtained respectively under different a-priori bound assumptions for the exact solution and suitable choices of regularization parameters. Finally, numerical examples are given to show that the proposed numerical method is effective.     

A Fourth-Order Modified Method for the Cauchy Problem of the Modified Helmholtz Equation

This paper is concerned with the Cauchy problem for the modified Helmholtz equation in an infinite strip domain 0 < x≤ 1, y ∈ R. The Cauchy data at x = 0 is given and the solution is then sought for the interval 0 < x ≤1. This problem is highly ill-posed and the solution (if it exists) does not depend continuously on the given data. In this paper, we propose a fourth-order modified method to solve the Cauchy problem. Convergence estimates are presented under the suitable choices of regularization parameters and the a priori assumption on the bounds of the exact solution. Numerical implementation is considered and the numerical examples show that our proposed method is effective and stable.     

Determining surface heat flux in the steady state for the Cauchy problem for the Laplace equation

In this paper, we consider the Cauchy problem for the Laplace equation, in a strip where the Cauchy data is given at x = 0 and the flux is sought in the interval 0<x?1. This problem is typical ill-posed: the solution (if it exists) does not depend continuously on the data. We study a modification of the equation, where a fourth-order mixed derivative term is added. Some error stability estimates for the flux are given, which show that the solution of the modified equation is approximate to the solution of the Cauchy problem for the Laplace equation. Furthermore, numerical examples show that the modified method works effectively.     

反对称次对称矩阵反问题的最小二乘解

The least-square solutions of inverse problem for anti-symmetric and skew-symmetric matrices are studied. In addition, the problem of using anti-symmetric and skew-symmetric matrices to construct the optimal approximation to a given matrix is discussed, the necessary and sufficient conditions for the problem are derived,and the expression of the solution is provided. A numerical example is given to show the effectiveness of the proposed method.     

Least-square Solutions of Inverse Problems for Anti-symmetric and Skew-symmetric Matrices

The least-square solutions of inverse problem for anti-symmetric and skew-symmetric matrices are studied. In addition, the problem of using anti-symmetric and skew-symmetric matrices to construct the optimal approximation to a given matrix is discussed, the necessary and sufficient conditions for the problem are derived, and the expression of the solution is provided. A numerical example is given to show the effectiveness of the proposed method.     

微分变换法求解二维非线性Volterra积分微分方程

为了解决二维非线性Volterra积分微分方程的求解问题,本文给出微分变换法.利用该方法将方程中的微分部分和积分部分进行变换,这样简化了原方程,进而得到非线性代数方程组,从而将原问题转换为求解非线性代数方程组的解,使得计算更简便.文中最后数值算例说明了该方法的可行性和有效性.     

矩阵不等式约束下矩阵方程AX=B的双对称解

本文讨论矩阵不等式CXD≥E 约束下矩阵方程AX=B的双对称解,即给定矩阵A,B,C,D和 E, 求双对称矩阵X, 使得AX=B 和 CXD≥E, 其中CXD≥E表示矩阵CXD-E非负.本文将问题转化为矩阵不等式最小非负偏差问题,利用极分解理论给出了求其解的迭代方法,并结合相关矩阵理论说明算法的收敛性.最后给出数值算例验证算法的有效性.     

The modified regularization method for identifying the unknown source on Poisson equation

This paper discusses the problem of determining an unknown source which depends only on one variable in two-dimensional Poisson equation from one supplementary temperature measurement at an internal point. The problem is ill-posed in the sense that the solution (if it exists) does not depend continuously on the data. The regularization solution is obtained by the modified regularization method. For the regularization solution, the Hölder type stability estimate between the regularization solution and the exact solution is given. Numerical results are presented to illustrate the accuracy and efficiency of this method.     

Determination of a spacewise dependent heat source

This paper investigates the inverse problem of determining a spacewise dependent heat source in the parabolic heat equation using the usual conditions of the direct problem and information from a supplementary temperature measurement at a given single instant of time. The spacewise dependent temperature measurement ensures that the inverse problem has a unique solution, but this solution is unstable, hence the problem is ill-posed. For this inverse problem, we propose an iterative algorithm based on a sequence of well-posed direct problems which are solved at each iteration step using the boundary element method (BEM). The instability is overcome by stopping the iterations at the first iteration for which the discrepancy principle is satisfied. Numerical results are presented for various typical benchmark test examples which have the input measured data perturbed by increasing amounts of random noise.     

The Cross-Entropy Method for Combinatorial and Continuous Optimization

We present a new and fast method, called the cross-entropy method, for finding the optimal solution of combinatorial and continuous nonconvex optimization problems with convex bounded domains. To find the optimal solution we solve a sequence of simple auxiliary smooth optimization problems based on Kullback-Leibler cross-entropy, importance sampling, Markov chain and Boltzmann distribution. We use importance sampling as an important ingredient for adaptive adjustment of the temperature in the Boltzmann distribution and use Kullback-Leibler cross-entropy to find the optimal solution. In fact, we use the mode of a unimodal importance sampling distribution, like the mode of beta distribution, as an estimate of the optimal solution for continuous optimization and Markov chains approach for combinatorial optimization. In the later case we show almost surely convergence of our algorithm to the optimal solution. Supporting numerical results for both continuous and combinatorial optimization problems are given as well. Our empirical studies suggest that the cross-entropy method has polynomial in the size of the problem running time complexity.