On the statistical estimation of diffusion processes: a survey

Data available on continuos-time diffusions are always sampled discretely in time. In most cases, the likelihood function of the observations is not directly computable. This survey covers a sample of the statistical methods that have been developed to solve this problem. We concentrate on some recent contributions to the literature based on three di§erent approaches to the problem: an improvement of the Euler-Maruyama discretization scheme, the use of Martingale Estimating Functions and the application of Generalized Method of Moments (GMM).

On the statistical estimation of diffusion processes - a partial survey

Data available on continuous-time diffusions are always sampled discretely in time. In most cases, the likelihood function of the observations is not directly computable. This survey covers a sample of the statistical methods that have been developed to solve this problem. We concentrate on some recent contributions to the literature based on three di§erent approaches to the problem: an improvement of the Euler-Maruyama discretization scheme, the employment of Martingale Estimating Functions, and the application of Generalized Method of Moments (GMM).

Contributions to the Theory of Optimal Tests

This paper considers tests which maximize the weighted average power (WAP). The focus is on determining WAP tests subject to an uncountable number of equalities and/or inequalities. The unifying theory allows us to obtain tests with correct size, similar tests, and unbiased tests, among others. A WAP test may be randomized and its characterization is not always possible. We show how to approximate the power of the optimal test by sequences of nonrandomized tests. Two alternative approximations are considered. The rst approach considers a sequence of similar tests for an increasing number of boundary conditions. This discretization allows us to implement the WAP tests in practice. The second method nds a sequence of tests which approximate the WAP test uniformly. This approximation allows us to show that WAP similar tests are admissible. The theoretical framework is readily applicable to several econometric models, including the important class of the curved-exponential family. In this paper, we consider the instrumental variable model with heteroskedastic and autocorrelated errors (HAC-IV) and the nearly integrated regressor model. In both models, we nd WAP similar and (locally) unbiased tests which dominate other available tests.

Comparison of methodologies for degree-day estimation using numerical methods

O desenvolvimento de projetos relacionados ao desempenho de diversas culturas tem recebido aperfeiçoamento cada vez maior, incorporado a modelos matemáticos sendo indispensável à utilização de equações cada vez mais consistentes que possibilitem previsão e maior aproximação do comportamento real, diminuindo o erro na obtenção das estimativas. Entre as operações unitárias que demandam maior estudo estão aquelas relacionadas com o crescimento da cultura, caracterizadas pela temperatura ideal para o acréscimo de matéria seca. Pelo amplo uso dos métodos matemáticos na representação, análise e obtenção de estimativas de graus-dia, juntamente com a grande importância que a cultura da cana-de-açúcar tem para a economia brasileira, foi realizada uma avaliação dos modelos matemáticos comumente usados e dos métodos numéricos de integração na estimativa da disponibilidade de graus-dia para essa cultura, na região de Botucatu, Estado de São Paulo. Os modelos de integração, com discretização de 6 em 6 h, apresentaram resultados satisfatórios na estimativa de graus-dia. As metodologias tradicionais apresentaram desempenhos satisfatórios quanto à estimativa de grausdia com base na curva de temperatura horária para cada dia e para os agrupamentos de três, sete, 15 e 30 dias. Pelo método numérico de integração, a região de Botucatu, Estado de São Paulo, apresentou disponibilidade térmica anual média de 1.070,6 GD para a cultura da cana-de-açúcar.

A continuously differentiable upwinding scheme for the simulation of fluid flow problems

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A combination of implicit and adaptative upwind tools for the numerical solution of incompressible free surface flows

This paper is concerned with the numerical solutions of time dependent two-dimensional incompressible flows. By using the primitive variables of velocity and pressure, the Navier-Stokes and mass conservation equations are solved by a semi-implicit finite difference projection method. A new bounded higher order upwind convection scheme is employed to deal with the non-linear (advective) terms. The procedure is an adaptation of the GENSMAC (J. Comput. Phys. 1994; 110: 171-186) methodology for calculating confined and free surface fluid flows at both low and high Reynolds numbers. The calculations were performed by using the 2D version of the Freeflow simulation system (J. Comp. Visual. Science 2000; 2:199-210). In order to demonstrate the capabilities of the numerical method, various test cases are presented. These are the fully developed flow in a channel, the flow over a backward facing step, the die-swell problem, the broken dam flow, and an impinging jet onto a flat plate. The numerical results compare favourably with the experimental data and the analytical solutions. Copyright (c) 2006 John Wiley & Sons, Ltd.

Assessment of a high-order finite difference upwind scheme for the simulation of convection-diffusion problems

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A marker-and-cell approach to free surface 2-D multiphase flows

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Simulation results and applications of an advection bounded scheme to practical flows

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Discretizing the deformation parameter in the su(q)(2) quantum algebra

Inspired in recent works of Biedenham [1, 2] on the realization of the q-algebra su(q)(2), We show in this note that the condition [2j + 1](q) = N-q(j) = integer, implies the discretization of the deformation parameter alpha, where q = e(alpha). This discretization replaces the continuum associated to ct by an infinite sequence alpha(1), alpha(2), alpha(3),..., obtained for the values of j, which label the irreps of su(q)(2). The algebraic properties of N-q(j) are discussed in some detail, including its role as a trace, which conducts to the Clebsch-Gordan series for the direct product of irreps. The consequences of this process of discretization are discussed and its possible applications are pointed out. Although not a necessary one, the present prescription is valuable due to its algebraic simplicity especially in the regime of appreciable values of alpha.

Quantum discrete phase space dynamics and its continuous limit

The main aspects of a discrete phase space formalism are presented and the discrete dynamical bracket, suitable for the description of time evolution in finite-dimensional spaces, is discussed. A set of operator bases is defined in such a way that the Weyl-Wigner formalism is shown to be obtained as a limiting case. In the same form, the Moyal bracket is shown to be the limiting case of the discrete dynamical bracket. The dynamics in quantum discrete phase spaces is shown not to be attained from discretization of the continuous case.

High-order upwinding and the hydraulic jump

This work is concerned with the computation of incompressible axisymmetric and fall three-dimensional free-surface flows. In particular, the circular-hydraulic jump is simulated and compared with approximate analytic solutions. However, the principal thrust of this paper is to provide a real problem as a test bed for comparing the many existing convective approximations. Their performance is compared; SMART, HLPA and VONOS emerge as acceptable upwinding methods for this problem. Copyright (C) 2002 John Wiley Sons, Ltd.

Wavelet-Galerkin method for one-dimensional elastoplasticity and damage problems: Constitutive modeling and computational aspects

This work presents an analysis of the wavelet-Galerkin method for one-dimensional elastoplastic-damage problems. Time-stepping algorithm for non-linear dynamics is presented. Numerical treatment of the constitutive models is developed by the use of return-mapping algorithm. For spacial discretization we can use wavelet-Galerkin method instead of standard finite element method. This approach allows to locate singularities. The discrete formulation developed can be applied to the simulation of one-dimensional problems for elastic-plastic-damage models. (C) 2007 Elsevier B.V. All rights reserved.

Discrete approximations for strict convex continuous time problems and duality

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

A radar reflectivity-runoff model for use in flood warning

This model connects directly the radar reflectivity data and hydrological variable runoff. The catchment is discretized in pixels (4 Km × 4 Km) with the same resolution of the CAPPI. Careful discretization is made so that every grid catchment pixel corresponds precisely to CAPPI grid cell. The basin is assumed a linear system and also time invariant. The forecast technique takes advantage of spatial and temporal resolutions obtained by the radar. The method uses only the measurements of the factor reflectivity distribution observed over the catchment area without using the reflectivity - rainfall rate transformation by the conventional Z-R relationships. The reflectivity values in each catchment pixel are translated to a gauging station by using a transfer function. This transfer function represents the travel time of the superficial water flowing through pixels in the drainage direction ending at the gauging station. The parameters used to compute the transfer function are concentration time and the physiographic catchment characteristics. -from Authors