Avaliação de métodos numéricos para precificação de derivativos: aplicação ao mercado brasileiro

Os objetivos deste trabalho foram (i) rever métodos numéricos para precificação de derivativos; e (ii) comparar os métodos assumindo que os preços de mercado refletem àqueles obtidos pela fórmula de Black Scholes para precificação de opções do tipo européia. Aplicamos estes métodos para precificar opções de compra da ações Telebrás. Os critérios de acurácia e de custo computacional foram utilizados para comparar os seguintes modelos binomial, Monte Carlo, e diferenças finitas. Os resultados indicam que o modelo binomial possui boa acurácia e custo baixo, seguido pelo Monte Carlo e diferenças finitas. Entretanto, o método Monte Carlo poderia ser usado quando o derivativo depende de mais de dois ativos-objetos. É recomendável usar o método de diferenças finitas quando se obtém uma equação diferencial parcial cuja solução é o valor do derivativo.

Avaliação de projetos na presença de volatilidade: estudo de caso de empresas no setor de telecomunicações

Quando as empresas decidem se devem ou não investir em determinado projeto de investimentos a longo prazo (horizonte de 5 a 10 anos), algumas metodologias alternativas ao Fluxo de Caixa Descontado (FCD) podem se tornar úteis tanto para confirmar a viabilidade do negócio como para indicar o melhor momento para iniciar o Empreendimento. As análises que levam em conta a incerteza dos fluxos de caixa futuros e flexibilidade na data de início do projeto podem ser construídos com a abordagem estocástica, usando metodologias como a solução de equações diferenciais que descrevem o movimento browniano. Sob determinadas condições, as oportunidades de investimentos em projetos podem ser tratados como se fossem opções reais de compra, sem data de vencimento, como no modelo proposto por McDonald-Siegel (1986), para a tomada de decisões e momento ótimo para o investimento. Este trabalho analisa a viabilidade de investimentos no mercado de telecomunicações usando modelos não determinísticos, onde a variável mais relevante é a dispersão dos retornos, ou seja, que a variância representa o risco associado a determinado empreendimento.

WeIfare costs of inflation the case for interest-bearing money and empiricaI Estimates for Brazil

We provide in this paper a closed fonn for the Welfare Cost of Inflation which we prove to be closer than Bailey's expression to the correct solution of the corresponding non-separable differential equation. Next, we extend this approach to an economy with interest-bearing money, once again presenting a better appoximation than the one given by Bailey's approach. Finally, empirical estimates for Brazil are presented.

Affine processes, arbitrage-Free Term structures of legendre polynomials,and option pricing

Multivariate Affine term structure models have been increasingly used for pricing derivatives in fixed income markets. In these models, uncertainty of the term structure is driven by a state vector, while the short rate is an affine function of this vector. The model is characterized by a specific form for the stochastic differential equation (SDE) for the evolution of the state vector. This SDE presents restrictions on its drift term which rule out arbitrages in the market. In this paper we solve the following inverse problem: Suppose the term structure of interest rates is modeled by a linear combination of Legendre polynomials with random coefficients. Is there any SDE for these coefficients which rules out arbitrages? This problem is of particular empirical interest because the Legendre model is an example of factor model with clear interpretation for each factor, in which regards movements of the term structure. Moreover, the Affine structure of the Legendre model implies knowledge of its conditional characteristic function. From the econometric perspective, we propose arbitrage-free Legendre models to describe the evolution of the term structure. From the pricing perspective, we follow Duffie et al. (2000) in exploring Legendre conditional characteristic functions to obtain a computational tractable method to price fixed income derivatives. Closing the article, the empirical section presents precise evidence on the reward of implementing arbitrage-free parametric term structure models: The ability of obtaining a good approximation for the state vector by simply using cross sectional data.

Controle inteligente de sistemas eletroidráulicos utilizando redes neurais artificiais

This work describes the development of a nonlinear control strategy for an electro-hydraulic actuated system. The system to be controlled is represented by a third order ordinary differential equation subject to a dead-zone input. The control strategy is based on a nonlinear control scheme, combined with an artificial intelligence algorithm, namely, the method of feedback linearization and an artificial neural network. It is shown that, when such a hard nonlinearity and modeling inaccuracies are considered, the nonlinear technique alone is not enough to ensure a good performance of the controller. Therefore, a compensation strategy based on artificial neural networks, which have been notoriously used in systems that require the simulation of the process of human inference, is used. The multilayer perceptron network and the radial basis functions network as well are adopted and mathematically implemented within the control law. On this basis, the compensation ability considering both networks is compared. Furthermore, the application of new intelligent control strategies for nonlinear and uncertain mechanical systems are proposed, showing that the combination of a nonlinear control methodology and artificial neural networks improves the overall control system performance. Numerical results are presented to demonstrate the efficacy of the proposed control system

Avaliacao de dados de secagem de suspensoes de polpas de frutas em leito de jorro com alimentacao intermitente

The objective of this work was the development and improvement of the mathematical models based on mass and heat balances, representing the drying transient process fruit pulp in spouted bed dryer with intermittent feeding. Mass and energy balance for drying, represented by a system of differential equations, were developed in Fortran language and adapted to the condition of intermittent feeding and mass accumulation. Were used the DASSL routine (Differential Algebraic System Solver) for solving the differential equation system and used a heuristic optimization algorithm in parameter estimation, the Particle Swarm algorithm. From the experimental data food drying, the differential models were used to determine the quantity of water and the drying air temperature at the exit of a spouted bed and accumulated mass of powder in the dryer. The models were validated using the experimental data of drying whose operating conditions, air temperature, flow rate and time intermittency, varied within the limits studied. In reviewing the results predicted, it was found that these models represent the experimental data of the kinetics of production and accumulation of powder and humidity and air temperature at the outlet of the dryer

Evolução top-hat hidrodinâmica em campos de velocidades peculiares primordiais

We investigate the cosmology of the vacuum energy decaying into cold dark matter according to thermodynamics description of Alcaniz & Lima. We apply this model to analyze the evolution of primordial density perturbations in the matter that gave rise to the first generation of structures bounded by gravity in the Universe, called Population III Objects. The analysis of the dynamics of those systems will involve the calculation of a differential equation system governing the evolution of perturbations to the case of two coupled fluids (dark matter and baryonic matter), modeled with a Top-Hat profile based in the perturbation of the hydrodynamics equations, an efficient analytical tool to study the properties of dark energy models such as the behavior of the linear growth factor and the linear growth index, physical quantities closely related to the fields of peculiar velocities at any time, for different models of dark energy. The properties and the dynamics of current Universe are analyzed through the exact analytical form of the linear growth factor of density fluctuations, taking into account the influence of several physical cooling mechanisms acting on the density fluctuations of the baryonic component of matter during the evolution of the clouds of matter, studied from the primordial hydrogen recombination. This study is naturally extended to more general models of dark energy with constant equation of state parameter in a flat Universe

Métodos numéricos para resolução de equações diferenciais ordinárias lineares baseados em interpolação por spline

In this work we have elaborated a spline-based method of solution of inicial value problems involving ordinary differential equations, with emphasis on linear equations. The method can be seen as an alternative for the traditional solvers such as Runge-Kutta, and avoids root calculations in the linear time invariant case. The method is then applied on a central problem of control theory, namely, the step response problem for linear EDOs with possibly varying coefficients, where root calculations do not apply. We have implemented an efficient algorithm which uses exclusively matrix-vector operations. The working interval (till the settling time) was determined through a calculation of the least stable mode using a modified power method. Several variants of the method have been compared by simulation. For general linear problems with fine grid, the proposed method compares favorably with the Euler method. In the time invariant case, where the alternative is root calculation, we have indications that the proposed method is competitive for equations of sifficiently high order.

Scattering and bound states of spinless particles in a mixed vector-scalar smooth step potential

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

Zeros of Jacobi functions of second kind

The number of zeros in (- 1, 1) of the Jacobi function of second kind Q(n)((alpha, beta)) (x), alpha, beta > - 1, i.e. The second solution of the differential equation(1 - x(2))y (x) + (beta - alpha - (alpha + beta + 2)x)y' (x) + n(n + alpha + beta + 1)y(x) = 0,is determined for every n is an element of N and for all values of the parameters alpha > - 1 and beta > - 1. It turns out that this number depends essentially on alpha and beta as well as on the specific normalization of the function Q(n)((alpha, beta)) (x). Interlacing properties of the zeros are also obtained. As a consequence of the main result, we determine the number of zeros of Laguerre's and Hermite's functions of second kind. (c) 2005 Elsevier B.V. All rights reserved.

Conformal Klein-Gordon equations and quasinormal modes

Using conformal coordinates associated with conformal relativity-associated with de Sitter spacetime homeomorphic projection into Minkowski spacetime-we obtain a conformal Klein-Gordon partial differential equation, which is intimately related to the production of quasi-normal modes (QNMs) oscillations, in the context of electromagnetic and/or gravitational perturbations around, e.g., black holes. While QNMs arise as the solution of a wave-like equation with a Poschl-Teller potential, here we deduce and analytically solve a conformal 'radial' d'Alembert-like equation, from which we derive QNMs formal solutions, in a proposed alternative to more completely describe QNMs. As a by-product we show that this 'radial' equation can be identified with a Schrodinger-like equation in which the potential is exactly the second Poschl-Teller potential, and it can shed some new light on the investigations concerning QNMs.

Correlation times in stochastic equations with delayed feedback and multiplicative noise

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

On some classes of exactly-solvable Klein-Gordon equations

In this work we discuss some exactly solvable Klein-Gordon equations. We basically discuss the existence of classes of potentials with different nonrelativistic limits, but which shares the intermediate effective Schroedinger differential equation. We comment about the possible use of relativistic exact solutions as approximations for nonrelativistic inexact potentials. (c) 2005 Elsevier B.V. All rights reserved.

Sliding Vector Fields via Slow-Fast Systems

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Analytical optimization of spread and change of exponential decay of generalized Wannier functions in one dimension

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