An inverse problem of identifying the coefficient of parabolic equation

This paper studies an inverse problem of identifying the coefficient of parabolic equation when the final observation is given, which has important application in a large fields of applied science. Based on the optimal control framework, the existence and necessary condition of the minimum for the control functional are established. Since the optimal control problem is nonconvex, one may not expect a unique solution. However, in this paper the solution is proved to be locally unique. After the necessary condition is transformed into an elliptic bilateral variational inequality, an algorithm and some numerical experiments are proposed in the paper. The numerical results show that the algorithm designed in this paper is stable and that the coefficient is recovered very well.     

A hybrid method for pricing European options based on multiple assets with transaction costs

The problem of pricing European options based on multiple assets with transaction costs is considered. These options include, for example, quality options and options on the minimum of two or more risky assets. The value of these options is the solution of a nonlinear parabolic partial differential equation subject to a final condition given by the payoff function associated with the option. A computationally efficient method to solve this final-value problem is proposed. This method is based on an asymptotic expansion of the required solution with respect to the parameters related to the transaction costs followed by the numerical solution of the linear partial differential equations obtained at each order in perturbation theory. The numerical solution of these linear problems involves an implicit finite-difference scheme for the parabolic equation and the use of the fast Fourier sine transform to solve the resulting elliptic problems. Numerical results obtained on test problems with the method proposed here are shown and discussed.     

Computational identification of the lowest space-wise dependent coefficient of a parabolic equation

In the present work, we consider a nonlinear inverse problem of identifying the lowest coefficient of a parabolic equation. The desired coefficient depends on spatial variables only. Additional information about the solution is given at the final time moment, i.e., we consider the final redefinition. An iterative process is used to evaluate the lowest coefficient, where at each iteration we solve the standard initial-boundary value problem for the parabolic equation. On the basis of the maximum principle for the solution of the differential problem, the monotonicity of the iterative process is established along with the fact that the coefficient is approached from above. The possibilities of the proposed computational algorithm are illustrated by numerical examples for a model two-dimensional problem.     

Singular solution to the problem of minimizing resource consumption

An iterative method of finding a singular solution to the problem of minimizing resource consumption has been developed. This method is based on the information about the finite control structure. A condition for existence of a singular solution is obtained. The limit value for transferring the time between the normal and the singular solutions is found. A relation between the variations of the control switching instants and the variations of the initial conditions of the adjoint system is found. A system of linear algebraic equations relating the variations of the initial conditions of the adjoint system to the deviations of the phase coordinates from a given final state of the system is obtained. The calculation algorithm and the results of modeling and numerical calculations are presented.     

A regularization method for solving the radially symmetric backward heat conduction problem

This work is devoted to solving the radially symmetric backward heat conduction problem, starting from the final temperature distribution. The problem is ill-posed: the solution (if it exists) does not depend continuously on the given data. A modified Tikhonov regularization method is proposed for solving this inverse problem. A quite sharp estimate of the error between the approximate solution and the exact solution is obtained with a suitable choice of regularization parameter. A numerical example is presented to verify the efficiency and accuracy of the method.     

一类求解分片延迟微分方程的线性多步法的散逸性

本文研究分片延迟微分方程本身及数值方法的散逸性问题．给出了一个关于此类问题本身散逸性的充分条件,同时得到了一类求解此类问题的线性多步法的数值散逸性结果,此结果表明所考虑的数值方法继承了方程本身的散逸性．数值试验进一步验证了理论结果的正确性．     

On the coefficient inverse problem of heat conduction in a degenerating domain

ABSTRACT

In the paper, we consider a coefficient inverse problem for the heat equation in a degenerating angular domain. It has been shown that the inverse problem for the homogeneous heat equation with homogeneous boundary conditions has a nontrivial solution up to a constant factor consistent with the integral condition. Moreover, the solution of the considered inverse problem is found in explicit form. In conclusion, statements of possible generalizations and the results of numerical calculations are given.     

Numerical solution of two‐dimensional parabolic equation subject to nonstandard boundary specifications using the pseudospectral Legendre method

In this research, the problem of solving the two‐dimensional parabolic equation subject to a given initial condition and nonlocal boundary specifications is considered. A technique based on the pseudospectral Legendre method is proposed for the numerical solution of the studied problem. Several examples are given and the numerical results are shown to demonstrate the efficiently of the newly proposed method. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A numerical technique for solving IHCPs using Tikhonov regularization method

This study is intended to provide a numerical algorithm for solving a one-dimensional inverse heat conduction problem. The given heat conduction equation, the boundary conditions, and the initial condition are presented in a dimensionless form. The numerical approach is developed based on the use of the solution to the auxiliary problem as a basis function. To regularize the resultant ill-conditioned linear system of equations, we apply the Tikhonov regularization method to obtain the stable numerical approximation to the solution.     

求多目标优化问题Pareto最优解集的方法

主要讨论了无约束多目标优化问题Pareto最优解集的求解方法,其中问题的目标函数是C1连续函数.给出了Pareto最优解集的一个充要条件,定义了α强有效解,并结合区间分析的方法,建立了求解无约束多目标优化问题Pareto最优解集的区间算法,理论分析和数值结果均表明该算法是可靠和有效的.     

小周期型复合材料稳态热传导问题的一种双尺度渐近展开收敛性分析

给出一类具有某种对称性小周期复合材料稳态热传导问题解的渐近表示方法.区别于传统多尺度计算方法,将计算过程中需要求解的关于单胞Q的H_per¹(Q)周期边值问题改为齐次边值问题,这样数值方法求解时协调元空间容易构造;另一方面传统的多尺度渐近解不满足原始问题的边界条件,新构造的渐近形式不仅满足原始问题的物理边界条件,同时保持一定的收敛阶,更能被工程上所接受.     

Block monotone iterative methods for numerical solutions of nonlinear elliptic equations

Summary. Two block monotone iterative schemes for a nonlinear algebraic system, which is a finite difference approximation of a nonlinear elliptic boundary-value problem, are presented and are shown to converge monotonically either from above or from below to a solution of the system. This monotone convergence result yields a computational algorithm for numerical solutions as well as an existence-comparison theorem of the system, including a sufficient condition for the uniqueness of the solution. An advantage of the block iterative schemes is that the Thomas algorithm can be used to compute numerical solutions of the sequence of iterations in the same fashion as for one-dimensional problems. The block iterative schemes are compared with the point monotone iterative schemes of Picard, Jacobi and Gauss-Seidel, and various theoretical comparison results among these monotone iterative schemes are given. These comparison results demonstrate that the sequence of iterations from the block iterative schemes converges faster than the corresponding sequence given by the point iterative schemes. Application of the iterative schemes is given to a logistic model problem in ecology and numerical ressults for a test problem with known analytical solution are given. Received August 1, 1993 / Revised version received November 7, 1994     

The construction of numerical integration rules of degree three for product regions

Numerical integration formulas in n-dimensional Euclidean space of degree three are discussed. In this paper, for the product regions a method is presented to construct numerical integration formulas of degree three with 2n real points and positive weights. The presented problem is a little different from those dealt with by other authors. All the corresponding one-dimensional integrals can be different from each other and they are also nonsymmetrical. In this paper an n-dimensional numerical integration problem is turned into n one-dimensional moment problems, which simplifies the construction process. Some explicit numerical formulas are given. Furthermore, a more generalized numerical integration problem is considered, which will shed light on the final solution to the third degree numerical integration problem.     

DEPENDENCE OF QUALITATIVE BEHAVIOR OF THE NUMERICAL SOLUTIONS ON THE IGNITION TEMPERATURE FOR A COMBUSTION MODEL

We study the dependence of qualitative behavior of the numerical solutions （obtained by a projective and upwind finite difference scheme） on the ignition temperature for a combustion model problem with general initial condition. Convergence to weak solution is proved under the Courant-Friedrichs-Lewy condition. Some condition on the ignition temperature is given to guarantee the solution containing a strong detonation wave or a weak detonation wave. Finally, we give some numerical examples which show that a strong detonation wave can be transformed to a weak detonation wave under some well-chosen ignition temperature.     

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     

A generalized inverse eigenvalue problem in structural dynamic model updating

This paper is concerned with the problem of the best approximation for a given matrix pencil under a given spectral constraint and a submatrix pencil constraint. Such a problem arises in structural dynamic model updating. By using the Moore–Penrose generalized inverse and the singular value decomposition (SVD) matrices, the solvability condition and the expression for the solution of the problem are presented. A numerical algorithm for solving the problem is developed.     

基于非均匀参数化的自由终端时间最优控制问题求解

针对自由终端时间最优控制问题,提出了一种基于非均匀控制向量参数化的数值解法.将控制时域离散化为不同长度的时间段,各时间段长度作为新的控制变量.通过引入标准化的时间变量,原问题转化为均匀参数化的固定终端时间最优控制问题.建立目标和约束函数的Hamilton函数,通过求解伴随方程获得目标和约束函数的梯度,采用序列二次规划(SQP)获得数值解.针对两个经典的化工过程自由终端时间最优控制问题进行仿真研究,验证了所提出算法的可行性和有效性.     

Boundary control of a hyperbolic system with one space variable

We consider a mixed problem in a half-strip for a hyperbolic system with one space variable and with constant coefficients. The control problem is to find boundary conditions ensuring that the system has a given state vector at a given instant of time. We study whether the problem is asymptotically solvable, i.e., whether there exists a sequence of boundary conditions such that the corresponding sequence of final state vectors uniformly converges to the given vector. We reduce the construction of a family of such sequences of boundary conditions with a function parameter to the solution of a Fredholm integral equation of the second kind and prove a sufficient condition for its unique solvability in terms of the problem data.     

带固定比例交易费率的跳扩散欧式期权的定价

考虑市场存在交易费率的跳扩散欧式期权的定价问题.由于交易费的存在使得传统的对冲方法不适用,我们将该问题转化为两元的随机控制问题.证明了带固定比例交易费率的跳扩散欧式期权的价格是对应的积分微分不等方程的约束粘性解,并通过马尔科夫链对变分问题进行离散,证明了在粘性意义下离散方法的收敛性.最后给出了数值结果.     

Accelerated monotone iterative methods for finite difference equations of reaction-diffusion

This paper is concerned with numerical methods for a finite difference system of reaction-diffusion-convection equation under nonlinear boundary condition. Various monotone iterative methods are presented, and each of these methods leads to an existence-comparison theorem as well as a computational algorithm for numerical solutions. The monotone property of the iterations gives improved upper and lower bounds of the solution in each iteration, and the rate of convergence of the iterations is either quadratic or nearly quadratic depending on the property of the nonlinear function. Application is given to a model problem from chemical engineering, and some numerical results, including a test problem with known analytical solution, are presented to illustrate the various rates of convergence of the iterations. Received November 2, 1995 / Revised version received February 10, 1997