No abstract

No abstract

A time-dependent three-dimensional acoustic scattering problem is considered. An incoming wave packet is scattered by a bounded, simply connected obstacle with locally Lipschitz boundary. The obstacle is assumed to have a constant boundary acoustic impedance. The limit cases of acoustically soft and acoustically hard obstacles are considered. The scattered acoustic field is the solution of an exterior problem for the wave equation. A new numerical method to compute the scattered acoustic field is proposed. This numerical method obtains the time-dependent scattered field as a superposition of time-harmonic acoustic waves and computes the time-harmonic acoustic waves by a new "operator expansion method." That is, the time-harmonic acoustic waves are solutions of an exterior boundary value problem for the Helmholtz equation. The method used to compute the time-harmonic waves improves on the method proposed by Misici, Pacelli, and Zirilli [J. Acoust. Soc. Am. 103, 106-113 (1998)] and is based on a "perturbative series" of the type of the one proposed in the operator expansion method by Milder [J. Acoust. Soc. Am. 89, 529-541 (1991)]. Computationally, the method is highly parallelizable with respect to time and space variables. Some numerical experiments on test problems obtained with a parallel implementation of the numerical method proposed are shown and discussed from the numerical and the physical point of view. The website: http://www.econ.unian.it/recchioni/w1 shows four animations relative to the numerical experiments.

Let us suppose that the dynamics of the stock prices and of their stochastic variance is described by the Heston model, that is by a system of two stochastic differential equations with a suitable initial condition. The aim of this paper is to estimate the parameters of the Heston model and one component of the initial condition, that is the initial stochastic variance, from the knowledge of the stock and option prices observed at discrete times. The option prices considered refer to an European call on the stock whose prices are described by the Heston model. The method proposed to solve this problem is based on a filtering technique to construct a likelihood function and on the maximization of the likelihood function obtained. The estimated parameters and initial value component are characterized as being a maximizer of the likelihood function subject to some constraints. The solution of the filtering problem, used to construct the likelihood function, is based on an integral representation of the fundamental solution of the Fokker-Planck equation associated to the Heston model, on the use of the wavelet expansions presented in (Fatone et al. in High performance algorithms based on a new wavelet expansion for time dependent acoustic obstacle scattering. Commun. Computat. Phys. (2007), Research Developments in Acoustics, vol. 2, pp. 39-69. Transworld Research Network, Kerala (2005), New wavelet bases made of piecewise polynomial functions: approximation theory, quadrature rules and applications to kernel sparsification and image compression. SIAM J. Sci. Comput. (submitted)) to approximate the integral kernel appearing in the representation formula of the fundamental solution, on a simple truncation procedure to exploit the sparsifying properties of the wavelet expansions and on the use of the fast Fourier transform (FFT). The use of these techniques generates a very efficient and fully parallelizable numerical procedure to solve the filtering problem, this last fact makes possible to evaluate very efficiently the likelihood function and its gradient. As a byproduct of the solution of the filtering problem we have developed a stochastic variance tracking technique that gives very good results in numerical experiments. The maximum likelihood problem used in the estimation procedure is a low dimensional constrained optimization problem, its solution with ad hoc techniques is justified by the computational cost of evaluating the likelihood function and its gradient. We use parallel computing and a variable metric steepest ascent method to solve the maximum likelihood problem. Some numerical examples of the estimation problem using synthetic and real data, that is data relative to an index of the Milano stock exchange (S&PMIB30), obtained with a parallel implementation of the previous numerical method are presented. Very impressive speed up factors are obtained in the numerical examples using the parallel implementation of the numerical method proposed. The website: http://www.econ.univpm.it/pacelli/mariani/finance/...

No abstract