Information-theoretic model selection for optimal prediction of stochastic dynamical systems from data

Science.gov (United States)

Darmon, David

2018-03-01

In the absence of mechanistic or phenomenological models of real-world systems, data-driven models become necessary. The discovery of various embedding theorems in the 1980s and 1990s motivated a powerful set of tools for analyzing deterministic dynamical systems via delay-coordinate embeddings of observations of their component states. However, in many branches of science, the condition of operational determinism is not satisfied, and stochastic models must be brought to bear. For such stochastic models, the tool set developed for delay-coordinate embedding is no longer appropriate, and a new toolkit must be developed. We present an information-theoretic criterion, the negative log-predictive likelihood, for selecting the embedding dimension for a predictively optimal data-driven model of a stochastic dynamical system. We develop a nonparametric estimator for the negative log-predictive likelihood and compare its performance to a recently proposed criterion based on active information storage. Finally, we show how the output of the model selection procedure can be used to compare candidate predictors for a stochastic system to an information-theoretic lower bound.

Information Dynamics of a Nonlinear Stochastic Nanopore System

Directory of Open Access Journals (Sweden)

Claire Gilpin

2018-03-01

Full Text Available Nanopores have become a subject of interest in the scientific community due to their potential uses in nanometer-scale laboratory and research applications, including infectious disease diagnostics and DNA sequencing. Additionally, they display behavioral similarity to molecular and cellular scale physiological processes. Recent advances in information theory have made it possible to probe the information dynamics of nonlinear stochastic dynamical systems, such as autonomously fluctuating nanopore systems, which has enhanced our understanding of the physical systems they model. We present the results of local (LER and specific entropy rate (SER computations from a simulation study of an autonomously fluctuating nanopore system. We learn that both metrics show increases that correspond to fluctuations in the nanopore current, indicating fundamental changes in information generation surrounding these fluctuations.

Stochastic ice stream dynamics.

Science.gov (United States)

Mantelli, Elisa; Bertagni, Matteo Bernard; Ridolfi, Luca

2016-08-09

Ice streams are narrow corridors of fast-flowing ice that constitute the arterial drainage network of ice sheets. Therefore, changes in ice stream flow are key to understanding paleoclimate, sea level changes, and rapid disintegration of ice sheets during deglaciation. The dynamics of ice flow are tightly coupled to the climate system through atmospheric temperature and snow recharge, which are known exhibit stochastic variability. Here we focus on the interplay between stochastic climate forcing and ice stream temporal dynamics. Our work demonstrates that realistic climate fluctuations are able to (i) induce the coexistence of dynamic behaviors that would be incompatible in a purely deterministic system and (ii) drive ice stream flow away from the regime expected in a steady climate. We conclude that environmental noise appears to be crucial to interpreting the past behavior of ice sheets, as well as to predicting their future evolution.

Hybrid Semantics of Stochastic Programs with Dynamic Reconfiguration

Directory of Open Access Journals (Sweden)

Alberto Policriti

2009-10-01

Full Text Available We begin by reviewing a technique to approximate the dynamics of stochastic programs --written in a stochastic process algebra-- by a hybrid system, suitable to capture a mixed discrete/continuous evolution. In a nutshell, the discrete dynamics is kept stochastic while the continuous evolution is given in terms of ODEs, and the overall technique, therefore, naturally associates a Piecewise Deterministic Markov Process with a stochastic program. The specific contribution in this work consists in an increase of the flexibility of the translation scheme, obtained by allowing a dynamic reconfiguration of the degree of discreteness/continuity of the semantics. We also discuss the relationships of this approach with other hybrid simulation strategies for biochemical systems.

Geometric integrators for stochastic rigid body dynamics

KAUST Repository

Tretyakov, Mikhail

2016-01-05

Geometric integrators play an important role in simulating dynamical systems on long time intervals with high accuracy. We will illustrate geometric integration ideas within the stochastic context, mostly on examples of stochastic thermostats for rigid body dynamics. The talk will be mainly based on joint recent work with Rusland Davidchak and Tom Ouldridge.

Geometric integrators for stochastic rigid body dynamics

KAUST Repository

Tretyakov, Mikhail

2016-01-01

Geometric integrators play an important role in simulating dynamical systems on long time intervals with high accuracy. We will illustrate geometric integration ideas within the stochastic context, mostly on examples of stochastic thermostats for rigid body dynamics. The talk will be mainly based on joint recent work with Rusland Davidchak and Tom Ouldridge.

Dynamic behaviors of spin-1/2 bilayer system within Glauber-type stochastic dynamics based on the effective-field theory

International Nuclear Information System (INIS)

Ertaş, Mehmet; Kantar, Ersin; Keskin, Mustafa

2014-01-01

The dynamic phase transitions (DPTs) and dynamic phase diagrams of the kinetic spin-1/2 bilayer system in the presence of a time-dependent oscillating external magnetic field are studied by using Glauber-type stochastic dynamics based on the effective-field theory with correlations for the ferromagnetic/ferromagnetic (FM/FM), antiferromagnetic/ferromagnetic (AFM/FM) and antiferromagnetic/antiferromagnetic (AFM/AFM) interactions. The time variations of average magnetizations and the temperature dependence of the dynamic magnetizations are investigated. The dynamic phase diagrams for the amplitude of the oscillating field versus temperature were presented. The results are compared with the results of the same system within Glauber-type stochastic dynamics based on the mean-field theory. - Highlights: • The Ising bilayer system is investigated within the Glauber dynamics based on EFT. • The time variations of average order parameters to find phases are studied. • The dynamic phase diagrams are found for the different interaction parameters. • The system displays the critical points as well as a re-entrant behavior

Dynamic behaviors of spin-1/2 bilayer system within Glauber-type stochastic dynamics based on the effective-field theory

Energy Technology Data Exchange (ETDEWEB)

Ertaş, Mehmet; Kantar, Ersin, E-mail: ersinkantar@erciyes.edu.tr; Keskin, Mustafa

2014-05-01

The dynamic phase transitions (DPTs) and dynamic phase diagrams of the kinetic spin-1/2 bilayer system in the presence of a time-dependent oscillating external magnetic field are studied by using Glauber-type stochastic dynamics based on the effective-field theory with correlations for the ferromagnetic/ferromagnetic (FM/FM), antiferromagnetic/ferromagnetic (AFM/FM) and antiferromagnetic/antiferromagnetic (AFM/AFM) interactions. The time variations of average magnetizations and the temperature dependence of the dynamic magnetizations are investigated. The dynamic phase diagrams for the amplitude of the oscillating field versus temperature were presented. The results are compared with the results of the same system within Glauber-type stochastic dynamics based on the mean-field theory. - Highlights: • The Ising bilayer system is investigated within the Glauber dynamics based on EFT. • The time variations of average order parameters to find phases are studied. • The dynamic phase diagrams are found for the different interaction parameters. • The system displays the critical points as well as a re-entrant behavior.

Dynamics of the stochastic low concentration trimolecular oscillatory chemical system with jumps

Science.gov (United States)

Wei, Yongchang; Yang, Qigui

2018-06-01

This paper is devoted to discern long time dynamics through the stochastic low concentration trimolecular oscillatory chemical system with jumps. By Lyapunov technique, this system is proved to have a unique global positive solution, and the asymptotic stability in mean square of such model is further established. Moreover, the existence of random attractor and Lyapunov exponents are obtained for the stochastic homeomorphism flow generated by the corresponding global positive solution. And some numerical simulations are given to illustrate the presented results.

Numerical simulation of stochastic point kinetic equation in the dynamical system of nuclear reactor

International Nuclear Information System (INIS)

Saha Ray, S.

2012-01-01

Highlights: ► In this paper stochastic neutron point kinetic equations have been analyzed. ► Euler–Maruyama method and Strong Taylor 1.5 order method have been discussed. ► These methods are applied for the solution of stochastic point kinetic equations. ► Comparison between the results of these methods and others are presented in tables. ► Graphs for neutron and precursor sample paths are also presented. -- Abstract: In the present paper, the numerical approximation methods, applied to efficiently calculate the solution for stochastic point kinetic equations () in nuclear reactor dynamics, are investigated. A system of Itô stochastic differential equations has been analyzed to model the neutron density and the delayed neutron precursors in a point nuclear reactor. The resulting system of Itô stochastic differential equations are solved over each time-step size. The methods are verified by considering different initial conditions, experimental data and over constant reactivities. The computational results indicate that the methods are simple and suitable for solving stochastic point kinetic equations. In this article, a numerical investigation is made in order to observe the random oscillations in neutron and precursor population dynamics in subcritical and critical reactors.

Safety Analysis of Stochastic Dynamical Systems

DEFF Research Database (Denmark)

Sloth, Christoffer; Wisniewski, Rafael

2015-01-01

This paper presents a method for verifying the safety of a stochastic system. In particular, we show how to compute the largest set of initial conditions such that a given stochastic system is safe with probability p. To compute the set of initial conditions we rely on the moment method that via...... that shows how the p-safe initial set is computed numerically....

Stochastic modeling and simulation of reaction-diffusion system with Hill function dynamics.

Science.gov (United States)

Chen, Minghan; Li, Fei; Wang, Shuo; Cao, Young

2017-03-14

Stochastic simulation of reaction-diffusion systems presents great challenges for spatiotemporal biological modeling and simulation. One widely used framework for stochastic simulation of reaction-diffusion systems is reaction diffusion master equation (RDME). Previous studies have discovered that for the RDME, when discretization size approaches zero, reaction time for bimolecular reactions in high dimensional domains tends to infinity. In this paper, we demonstrate that in the 1D domain, highly nonlinear reaction dynamics given by Hill function may also have dramatic change when discretization size is smaller than a critical value. Moreover, we discuss methods to avoid this problem: smoothing over space, fixed length smoothing over space and a hybrid method. Our analysis reveals that the switch-like Hill dynamics reduces to a linear function of discretization size when the discretization size is small enough. The three proposed methods could correctly (under certain precision) simulate Hill function dynamics in the microscopic RDME system.

Sparse learning of stochastic dynamical equations

Science.gov (United States)

Boninsegna, Lorenzo; Nüske, Feliks; Clementi, Cecilia

2018-06-01

With the rapid increase of available data for complex systems, there is great interest in the extraction of physically relevant information from massive datasets. Recently, a framework called Sparse Identification of Nonlinear Dynamics (SINDy) has been introduced to identify the governing equations of dynamical systems from simulation data. In this study, we extend SINDy to stochastic dynamical systems which are frequently used to model biophysical processes. We prove the asymptotic correctness of stochastic SINDy in the infinite data limit, both in the original and projected variables. We discuss algorithms to solve the sparse regression problem arising from the practical implementation of SINDy and show that cross validation is an essential tool to determine the right level of sparsity. We demonstrate the proposed methodology on two test systems, namely, the diffusion in a one-dimensional potential and the projected dynamics of a two-dimensional diffusion process.

The Stochastic Dynamics for Ecological Tourism System with Visitor Educational Intervention

Directory of Open Access Journals (Sweden)

Dongping Wei

2013-01-01

Full Text Available The ever-increasing visitation in parks and protected areas continues to present a considerable challenge for worldwide land managers with allowing recreational use while preserving natural conditions. In China, the fast expanding visitation in protected areas is quickly damaging the natural resources and precious culture without effective visitor education, while regulation and site management are also gaining very limited efficacy. We propose a differential equation to describe the ecological tourism system. Shown by the theoretical proof and numerical simulation, the ecological tourism system is unstable without any perturbed factors, especially visitor educational intervention, because the solution of the dynamic system explodes in a finite time given any initial value. Supposing that the intrinsic increasing rate of stakeholders in the systems stochastically perturbed by the visitor educational intervention, we discover that the stochastic dynamic model can effectively suppress the explosion of the solution. As such, we demonstrate that the tourism system can develop steadily and safely even under a large amount of visitors in public vacation, when employing continuous visitor education intervention programmes.

Research on nonlinear stochastic dynamical price model

International Nuclear Information System (INIS)

Li Jiaorui; Xu Wei; Xie Wenxian; Ren Zhengzheng

2008-01-01

In consideration of many uncertain factors existing in economic system, nonlinear stochastic dynamical price model which is subjected to Gaussian white noise excitation is proposed based on deterministic model. One-dimensional averaged Ito stochastic differential equation for the model is derived by using the stochastic averaging method, and applied to investigate the stability of the trivial solution and the first-passage failure of the stochastic price model. The stochastic price model and the methods presented in this paper are verified by numerical studies

Potential and flux field landscape theory. II. Non-equilibrium thermodynamics of spatially inhomogeneous stochastic dynamical systems

International Nuclear Information System (INIS)

Wu, Wei; Wang, Jin

2014-01-01

We have established a general non-equilibrium thermodynamic formalism consistently applicable to both spatially homogeneous and, more importantly, spatially inhomogeneous systems, governed by the Langevin and Fokker-Planck stochastic dynamics with multiple state transition mechanisms, using the potential-flux landscape framework as a bridge connecting stochastic dynamics with non-equilibrium thermodynamics. A set of non-equilibrium thermodynamic equations, quantifying the relations of the non-equilibrium entropy, entropy flow, entropy production, and other thermodynamic quantities, together with their specific expressions, is constructed from a set of dynamical decomposition equations associated with the potential-flux landscape framework. The flux velocity plays a pivotal role on both the dynamic and thermodynamic levels. On the dynamic level, it represents a dynamic force breaking detailed balance, entailing the dynamical decomposition equations. On the thermodynamic level, it represents a thermodynamic force generating entropy production, manifested in the non-equilibrium thermodynamic equations. The Ornstein-Uhlenbeck process and more specific examples, the spatial stochastic neuronal model, in particular, are studied to test and illustrate the general theory. This theoretical framework is particularly suitable to study the non-equilibrium (thermo)dynamics of spatially inhomogeneous systems abundant in nature. This paper is the second of a series

Optimal Stochastic Control Problem for General Linear Dynamical Systems in Neuroscience

Directory of Open Access Journals (Sweden)

Yan Chen

2017-01-01

Full Text Available This paper considers a d-dimensional stochastic optimization problem in neuroscience. Suppose the arm’s movement trajectory is modeled by high-order linear stochastic differential dynamic system in d-dimensional space, the optimal trajectory, velocity, and variance are explicitly obtained by using stochastic control method, which allows us to analytically establish exact relationships between various quantities. Moreover, the optimal trajectory is almost a straight line for a reaching movement; the optimal velocity bell-shaped and the optimal variance are consistent with the experimental Fitts law; that is, the longer the time of a reaching movement, the higher the accuracy of arriving at the target position, and the results can be directly applied to designing a reaching movement performed by a robotic arm in a more general environment.

Introduction to stochastic dynamic programming

CERN Document Server

Ross, Sheldon M; Lukacs, E

1983-01-01

Introduction to Stochastic Dynamic Programming presents the basic theory and examines the scope of applications of stochastic dynamic programming. The book begins with a chapter on various finite-stage models, illustrating the wide range of applications of stochastic dynamic programming. Subsequent chapters study infinite-stage models: discounting future returns, minimizing nonnegative costs, maximizing nonnegative returns, and maximizing the long-run average return. Each of these chapters first considers whether an optimal policy need exist-providing counterexamples where appropriate-and the

Reliability-based Dynamic Network Design with Stochastic Networks

NARCIS (Netherlands)

Li, H.

2009-01-01

Transportation systems are stochastic and dynamic systems. The road capacities and the travel demand are fluctuating from time to time within a day and at the same time from day to day. For road users, the travel time and travel costs experienced over time and space are stochastic, thus desire

Stochastic dynamics of a delayed bistable system with multiplicative noise

Energy Technology Data Exchange (ETDEWEB)

Dung, Nguyen Tien, E-mail: dung-nguyentien10@yahoo.com, E-mail: dungnt@fpt.edu.vn [Department of Mathematics, FPT University, No 8 Ton That Thuyet, My Dinh, Tu Liem, Hanoi (Viet Nam)

2014-05-15

In this paper we investigate the properties of a delayed bistable system under the effect of multiplicative noise. We first prove the existence and uniqueness of the positive solution and show that its moments are uniformly bounded. Then, we study stochastic dynamics of the solution in long time, the lower and upper bounds for the paths and an estimate for the average value are provided.

Optimal Linear Responses for Markov Chains and Stochastically Perturbed Dynamical Systems

Science.gov (United States)

Antown, Fadi; Dragičević, Davor; Froyland, Gary

2018-03-01

The linear response of a dynamical system refers to changes to properties of the system when small external perturbations are applied. We consider the little-studied question of selecting an optimal perturbation so as to (i) maximise the linear response of the equilibrium distribution of the system, (ii) maximise the linear response of the expectation of a specified observable, and (iii) maximise the linear response of the rate of convergence of the system to the equilibrium distribution. We also consider the inhomogeneous, sequential, or time-dependent situation where the governing dynamics is not stationary and one wishes to select a sequence of small perturbations so as to maximise the overall linear response at some terminal time. We develop the theory for finite-state Markov chains, provide explicit solutions for some illustrative examples, and numerically apply our theory to stochastically perturbed dynamical systems, where the Markov chain is replaced by a matrix representation of an approximate annealed transfer operator for the random dynamical system.

Stochastic Reachability Analysis of Hybrid Systems

CERN Document Server

Bujorianu, Luminita Manuela

2012-01-01

Stochastic reachability analysis (SRA) is a method of analyzing the behavior of control systems which mix discrete and continuous dynamics. For probabilistic discrete systems it has been shown to be a practical verification method but for stochastic hybrid systems it can be rather more. As a verification technique SRA can assess the safety and performance of, for example, autonomous systems, robot and aircraft path planning and multi-agent coordination but it can also be used for the adaptive control of such systems. Stochastic Reachability Analysis of Hybrid Systems is a self-contained and accessible introduction to this novel topic in the analysis and development of stochastic hybrid systems. Beginning with the relevant aspects of Markov models and introducing stochastic hybrid systems, the book then moves on to coverage of reachability analysis for stochastic hybrid systems. Following this build up, the core of the text first formally defines the concept of reachability in the stochastic framework and then...

Collisionally induced stochastic dynamics of fast ions in solids

International Nuclear Information System (INIS)

Burgdoerfer, J.

1989-01-01

Recent developments in the theory of excited state formation in collisions of fast highly charged ions with solids are reviewed. We discuss a classical transport theory employing Monte-Carlo sampling of solutions of a microscopic Langevin equation. Dynamical screening by the dielectric medium as well as multiple collisions are incorporated through the drift and stochastic forces in the Langevin equation. The close relationship between the extrinsically stochastic dynamics described by the Langevin and the intrinsic stochasticity in chaotic nonlinear dynamical systems is stressed. Comparison with experimental data and possible modification by quantum corrections are discussed. 49 refs., 11 figs

Nonlinear stochastic dynamics of mesoscopic homogeneous biochemical reaction systems—an analytical theory

International Nuclear Information System (INIS)

Qian, Hong

2011-01-01

The nonlinear dynamics of biochemical reactions in a small-sized system on the order of a cell are stochastic. Assuming spatial homogeneity, the populations of n molecular species follow a multi-dimensional birth-and-death process on Z n . We introduce the Delbrück–Gillespie process, a continuous-time Markov jump process, whose Kolmogorov forward equation has been known as the chemical master equation, and whose stochastic trajectories can be computed via the Gillespie algorithm. Using simple models, we illustrate that a system of nonlinear ordinary differential equations on R n emerges in the infinite system size limit. For finite system size, transitions among multiple attractors of the nonlinear dynamical system are rare events with exponentially long transit times. There is a separation of time scales between the deterministic ODEs and the stochastic Markov jumps between attractors. No diffusion process can provide a global representation that is accurate on both short and long time scales for the nonlinear, stochastic population dynamics. On the short time scale and near deterministic stable fixed points, Ornstein–Uhlenbeck Gaussian processes give linear stochastic dynamics that exhibit time-irreversible circular motion for open, driven chemical systems. Extending this individual stochastic behaviour-based nonlinear population theory of molecular species to other biological systems is discussed. (invited article)

Stochastic dynamics of melt ponds and sea ice-albedo climate feedback

Science.gov (United States)

Sudakov, Ivan

Evolution of melt ponds on the Arctic sea surface is a complicated stochastic process. We suggest a low-order model with ice-albedo feedback which describes stochastic dynamics of melt ponds geometrical characteristics. The model is a stochastic dynamical system model of energy balance in the climate system. We describe the equilibria in this model. We conclude the transition in fractal dimension of melt ponds affects the shape of the sea ice albedo curve.

The stochastic system approach for estimating dynamic treatments effect.

Science.gov (United States)

Commenges, Daniel; Gégout-Petit, Anne

2015-10-01

The problem of assessing the effect of a treatment on a marker in observational studies raises the difficulty that attribution of the treatment may depend on the observed marker values. As an example, we focus on the analysis of the effect of a HAART on CD4 counts, where attribution of the treatment may depend on the observed marker values. This problem has been treated using marginal structural models relying on the counterfactual/potential response formalism. Another approach to causality is based on dynamical models, and causal influence has been formalized in the framework of the Doob-Meyer decomposition of stochastic processes. Causal inference however needs assumptions that we detail in this paper and we call this approach to causality the "stochastic system" approach. First we treat this problem in discrete time, then in continuous time. This approach allows incorporating biological knowledge naturally. When working in continuous time, the mechanistic approach involves distinguishing the model for the system and the model for the observations. Indeed, biological systems live in continuous time, and mechanisms can be expressed in the form of a system of differential equations, while observations are taken at discrete times. Inference in mechanistic models is challenging, particularly from a numerical point of view, but these models can yield much richer and reliable results.

Automated Flight Routing Using Stochastic Dynamic Programming

Science.gov (United States)

Ng, Hok K.; Morando, Alex; Grabbe, Shon

2010-01-01

Airspace capacity reduction due to convective weather impedes air traffic flows and causes traffic congestion. This study presents an algorithm that reroutes flights in the presence of winds, enroute convective weather, and congested airspace based on stochastic dynamic programming. A stochastic disturbance model incorporates into the reroute design process the capacity uncertainty. A trajectory-based airspace demand model is employed for calculating current and future airspace demand. The optimal routes minimize the total expected traveling time, weather incursion, and induced congestion costs. They are compared to weather-avoidance routes calculated using deterministic dynamic programming. The stochastic reroutes have smaller deviation probability than the deterministic counterpart when both reroutes have similar total flight distance. The stochastic rerouting algorithm takes into account all convective weather fields with all severity levels while the deterministic algorithm only accounts for convective weather systems exceeding a specified level of severity. When the stochastic reroutes are compared to the actual flight routes, they have similar total flight time, and both have about 1% of travel time crossing congested enroute sectors on average. The actual flight routes induce slightly less traffic congestion than the stochastic reroutes but intercept more severe convective weather.

The Dynamic Programming Method of Stochastic Differential Game for Functional Forward-Backward Stochastic System

Directory of Open Access Journals (Sweden)

Shaolin Ji

2013-01-01

Full Text Available This paper is devoted to a stochastic differential game (SDG of decoupled functional forward-backward stochastic differential equation (FBSDE. For our SDG, the associated upper and lower value functions of the SDG are defined through the solution of controlled functional backward stochastic differential equations (BSDEs. Applying the Girsanov transformation method introduced by Buckdahn and Li (2008, the upper and the lower value functions are shown to be deterministic. We also generalize the Hamilton-Jacobi-Bellman-Isaacs (HJBI equations to the path-dependent ones. By establishing the dynamic programming principal (DPP, we derive that the upper and the lower value functions are the viscosity solutions of the corresponding upper and the lower path-dependent HJBI equations, respectively.

Effects of stochastic time-delayed feedback on a dynamical system modeling a chemical oscillator

Science.gov (United States)

González Ochoa, Héctor O.; Perales, Gualberto Solís; Epstein, Irving R.; Femat, Ricardo

2018-05-01

We examine how stochastic time-delayed negative feedback affects the dynamical behavior of a model oscillatory reaction. We apply constant and stochastic time-delayed negative feedbacks to a point Field-Körös-Noyes photosensitive oscillator and compare their effects. Negative feedback is applied in the form of simulated inhibitory electromagnetic radiation with an intensity proportional to the concentration of oxidized light-sensitive catalyst in the oscillator. We first characterize the system under nondelayed inhibitory feedback; then we explore and compare the effects of constant (deterministic) versus stochastic time-delayed feedback. We find that the oscillatory amplitude, frequency, and waveform are essentially preserved when low-dispersion stochastic delayed feedback is used, whereas small but measurable changes appear when a large dispersion is applied.

Exact solutions to chaotic and stochastic systems

Science.gov (United States)

González, J. A.; Reyes, L. I.; Guerrero, L. E.

2001-03-01

We investigate functions that are exact solutions to chaotic dynamical systems. A generalization of these functions can produce truly random numbers. For the first time, we present solutions to random maps. This allows us to check, analytically, some recent results about the complexity of random dynamical systems. We confirm the result that a negative Lyapunov exponent does not imply predictability in random systems. We test the effectiveness of forecasting methods in distinguishing between chaotic and random time series. Using the explicit random functions, we can give explicit analytical formulas for the output signal in some systems with stochastic resonance. We study the influence of chaos on the stochastic resonance. We show, theoretically, the existence of a new type of solitonic stochastic resonance, where the shape of the kink is crucial. Using our models we can predict specific patterns in the output signal of stochastic resonance systems.

Path integral methods for the dynamics of stochastic and disordered systems

International Nuclear Information System (INIS)

Hertz, John A; Roudi, Yasser; Sollich, Peter

2017-01-01

We review some of the techniques used to study the dynamics of disordered systems subject to both quenched and fast (thermal) noise. Starting from the Martin–Siggia–Rose/Janssen–De Dominicis–Peliti path integral formalism for a single variable stochastic dynamics, we provide a pedagogical survey of the perturbative, i.e. diagrammatic, approach to dynamics and how this formalism can be used for studying soft spin models. We review the supersymmetric formulation of the Langevin dynamics of these models and discuss the physical implications of the supersymmetry. We also describe the key steps involved in studying the disorder-averaged dynamics. Finally, we discuss the path integral approach for the case of hard Ising spins and review some recent developments in the dynamics of such kinetic Ising models. (topical review)

A deterministic and stochastic model for the system dynamics of tumor-immune responses to chemotherapy

Science.gov (United States)

Liu, Xiangdong; Li, Qingze; Pan, Jianxin

2018-06-01

Modern medical studies show that chemotherapy can help most cancer patients, especially for those diagnosed early, to stabilize their disease conditions from months to years, which means the population of tumor cells remained nearly unchanged in quite a long time after fighting against immune system and drugs. In order to better understand the dynamics of tumor-immune responses under chemotherapy, deterministic and stochastic differential equation models are constructed to characterize the dynamical change of tumor cells and immune cells in this paper. The basic dynamical properties, such as boundedness, existence and stability of equilibrium points, are investigated in the deterministic model. Extended stochastic models include stochastic differential equations (SDEs) model and continuous-time Markov chain (CTMC) model, which accounts for the variability in cellular reproduction, growth and death, interspecific competitions, and immune response to chemotherapy. The CTMC model is harnessed to estimate the extinction probability of tumor cells. Numerical simulations are performed, which confirms the obtained theoretical results.

System Entropy Measurement of Stochastic Partial Differential Systems

Directory of Open Access Journals (Sweden)

Bor-Sen Chen

2016-03-01

Full Text Available System entropy describes the dispersal of a system’s energy and is an indication of the disorder of a physical system. Several system entropy measurement methods have been developed for dynamic systems. However, most real physical systems are always modeled using stochastic partial differential dynamic equations in the spatio-temporal domain. No efficient method currently exists that can calculate the system entropy of stochastic partial differential systems (SPDSs in consideration of the effects of intrinsic random fluctuation and compartment diffusion. In this study, a novel indirect measurement method is proposed for calculating of system entropy of SPDSs using a Hamilton–Jacobi integral inequality (HJII-constrained optimization method. In other words, we solve a nonlinear HJII-constrained optimization problem for measuring the system entropy of nonlinear stochastic partial differential systems (NSPDSs. To simplify the system entropy measurement of NSPDSs, the global linearization technique and finite difference scheme were employed to approximate the nonlinear stochastic spatial state space system. This allows the nonlinear HJII-constrained optimization problem for the system entropy measurement to be transformed to an equivalent linear matrix inequalities (LMIs-constrained optimization problem, which can be easily solved using the MATLAB LMI-toolbox (MATLAB R2014a, version 8.3. Finally, several examples are presented to illustrate the system entropy measurement of SPDSs.

Path to Stochastic Stability: Comparative Analysis of Stochastic Learning Dynamics in Games

KAUST Repository

Jaleel, Hassan

2018-04-08

Stochastic stability is a popular solution concept for stochastic learning dynamics in games. However, a critical limitation of this solution concept is its inability to distinguish between different learning rules that lead to the same steady-state behavior. We address this limitation for the first time and develop a framework for the comparative analysis of stochastic learning dynamics with different update rules but same steady-state behavior. We present the framework in the context of two learning dynamics: Log-Linear Learning (LLL) and Metropolis Learning (ML). Although both of these dynamics have the same stochastically stable states, LLL and ML correspond to different behavioral models for decision making. Moreover, we demonstrate through an example setup of sensor coverage game that for each of these dynamics, the paths to stochastically stable states exhibit distinctive behaviors. Therefore, we propose multiple criteria to analyze and quantify the differences in the short and medium run behavior of stochastic learning dynamics. We derive and compare upper bounds on the expected hitting time to the set of Nash equilibria for both LLL and ML. For the medium to long-run behavior, we identify a set of tools from the theory of perturbed Markov chains that result in a hierarchical decomposition of the state space into collections of states called cycles. We compare LLL and ML based on the proposed criteria and develop invaluable insights into the comparative behavior of the two dynamics.

The dynamic behaviour of a non-stationary elevator compensating rope system under harmonic and stochastic excitations

Energy Technology Data Exchange (ETDEWEB)

Kaczmarczyk, S [School of Applied Sciences, University of Northampton, St. George' s Avenue, Northampton NN2 6JD (United Kingdom); Iwankiewicz, R [Institute of Mechanics and Ocean Engineering, Hamburg University of Technology, Eissendorfer Strasse 42 D-21073, Hamburg (Germany); Terumichi, Y, E-mail: stefan.kaczmarczyk@northampton.ac.u [Faculty of Science and Technology, Sophia University, 7-1 KIOI-CHO, CHIYODAKU, Tokyo, 102-8554 Japan (Japan)

2009-08-01

Moving slender elastic elements such as ropes, cables and belts are pivotal components of vertical transportation systems such as traction elevators. Their lengths vary within the host building structure during the elevator operation which results in the change of the mass and stiffness characteristics of the system. The structure of modern high-rise buildings is flexible and when subjected to loads due to strong winds and earthquakes it vibrates at low frequencies. The inertial load induced by the building motion excites the flexible components of the elevator system. The compensating ropes due to their lower tension are particularly affected and undergo large dynamic deformations. The paper focuses on the presentation of the non-stationary model of a building-compensating rope system and on the analysis to predict its dynamic response. The excitation mechanism is represented by a harmonic process and the results of computer simulations to predict transient resonance response are presented. The analysis of the simulation results leads to recommendations concerning the selection of the weight of the compensation assembly to minimize the effects of an adverse dynamic response of the system. The scenario when the excitation is represented as a narrow-band stochastic process with the state vector governed by stochastic equations is then discussed and the stochastic differential equations governing the second-order statistical moments of the state vector are developed.

Emerging of Stochastic Dynamical Equalities and Steady State Thermodynamics from Darwinian Dynamics

International Nuclear Information System (INIS)

Ao, P.

2008-01-01

The evolutionary dynamics first conceived by Darwin and Wallace, referring to as Darwinian dynamics in the present paper, has been found to be universally valid in biology. The statistical mechanics and thermodynamics, while enormous successful in physics, have been in an awkward situation of wanting a consistent dynamical understanding. Here we present from a formal point of view an exploration of the connection between thermodynamics and Darwinian dynamics and a few related topics. We first show that the stochasticity in Darwinian dynamics implies the existence temperature, hence the canonical distribution of Boltzmann-Gibbs type. In term of relative entropy the Second Law of thermodynamics is dynamically demonstrated without detailed balance condition, and is valid regardless of size of the system. In particular, the dynamical component responsible for breaking detailed balance condition does not contribute to the change of the relative entropy. Two types of stochastic dynamical equalities of current interest are explicitly discussed in the present approach: One is based on Feynman-Kac formula and another is a generalization of Einstein relation. Both are directly accessible to experimental tests. Our demonstration indicates that Darwinian dynamics represents logically a simple and straightforward starting point for statistical mechanics and thermodynamics and is complementary to and consistent with conservative dynamics that dominates the physical sciences. Present exploration suggests the existence of a unified stochastic dynamical framework both near and far from equilibrium

Stochastic control theory dynamic programming principle

CERN Document Server

Nisio, Makiko

2015-01-01

This book offers a systematic introduction to the optimal stochastic control theory via the dynamic programming principle, which is a powerful tool to analyze control problems. First we consider completely observable control problems with finite horizons. Using a time discretization we construct a nonlinear semigroup related to the dynamic programming principle (DPP), whose generator provides the Hamilton–Jacobi–Bellman (HJB) equation, and we characterize the value function via the nonlinear semigroup, besides the viscosity solution theory. When we control not only the dynamics of a system but also the terminal time of its evolution, control-stopping problems arise. This problem is treated in the same frameworks, via the nonlinear semigroup. Its results are applicable to the American option price problem. Zero-sum two-player time-homogeneous stochastic differential games and viscosity solutions of the Isaacs equations arising from such games are studied via a nonlinear semigroup related to DPP (the min-ma...

Stochastic beam dynamics in storage rings

International Nuclear Information System (INIS)

Pauluhn, A.

1993-12-01

In this thesis several approaches to stochastic dynamics in storage rings are investigated. In the first part the theory of stochastic differential equations and Fokker-Planck equations is used to describe the processes which have been assumed to be Markov processes. The mathematical theory of Markov processes is well known. Nevertheless, analytical solutions can be found only in special cases and numerical algorithms are required. Several numerical integration schemes for stochastic differential equations will therefore be tested in analytical solvable examples and then applied to examples from accelerator physics. In particular the stochastically perturbed synchrotron motion is treated. For the special case of a double rf system several perturbation theoretical methods for deriving the Fokker-Planck equation in the action variable are used and compared with numerical results. The second part is concerned with the dynamics of electron storage rings. Due to the synchrotron radiation the electron motion is influenced by damping and exciting forces. An algorithm for the computation of the density function in the phase space of such a dissipative stochastically excited system is introduced. The density function contains all information of a process, e.g. it determines the beam dimensions and the lifetime of a stored electron beam. The new algorithm consists in calculating a time propagator for the density function. By means of this propagator the time evolution of the density is modelled very computing time efficient. The method is applied to simple models of the beam-beam interaction (one-dimensional, round beams) and the results of the density calculations are compared with results obtained from multiparticle tracking. Furthermore some modifications of the algorithm are introduced to improve its efficiency concerning computing time and storage requirements. Finally, extensions to two-dimensional beam-beam models are described. (orig.)

Long-Time Dynamic Response and Stochastic Resonance of Subdiffusive Overdamped Bistable Fractional Fokker-Planck Systems

International Nuclear Information System (INIS)

Yan-Mei, Kang; Yao-Lin, Jiang

2008-01-01

To explore the influence of anomalous diffusion on stochastic resonance (SR) more deeply and effectively, the method of moments is extended to subdiffusive overdamped bistable fractional Fokker-Planck systems for calculating the long-time linear dynamic response. It is found that the method of moments attains high accuracy with the truncation order N = 10, and in normal diffusion such obtained spectral amplification factor (SAF) of the first-order harmonic is also confirmed by stochastic simulation. Observing the SAF of the odd-order harmonics we find some interesting results, i.e. for smaller driving frequency the decrease of sub diffusion exponent inhibits the stochastic resonance (SR), while for larger driving frequency the decrease of sub diffusion exponent enhances the second SR peak, but the first one vanishes and a double SR is induced in the third-order harmonic at the same time. These observations suggest that the anomalous diffusion has important influence on the bistable dynamics

Efficient estimators for likelihood ratio sensitivity indices of complex stochastic dynamics

Energy Technology Data Exchange (ETDEWEB)

Arampatzis, Georgios; Katsoulakis, Markos A.; Rey-Bellet, Luc [Department of Mathematics and Statistics, University of Massachusetts, Amherst, Massachusetts 01003 (United States)

2016-03-14

We demonstrate that centered likelihood ratio estimators for the sensitivity indices of complex stochastic dynamics are highly efficient with low, constant in time variance and consequently they are suitable for sensitivity analysis in long-time and steady-state regimes. These estimators rely on a new covariance formulation of the likelihood ratio that includes as a submatrix a Fisher information matrix for stochastic dynamics and can also be used for fast screening of insensitive parameters and parameter combinations. The proposed methods are applicable to broad classes of stochastic dynamics such as chemical reaction networks, Langevin-type equations and stochastic models in finance, including systems with a high dimensional parameter space and/or disparate decorrelation times between different observables. Furthermore, they are simple to implement as a standard observable in any existing simulation algorithm without additional modifications.

The intrinsic stochasticity of near-integrable Hamiltonian systems

Energy Technology Data Exchange (ETDEWEB)

Krlin, L [Ceskoslovenska Akademie Ved, Prague (Czechoslovakia). Ustav Fyziky Plazmatu

1989-09-01

Under certain conditions, the dynamics of near-integrable Hamiltonian systems appears to be stochastic. This stochasticity (intrinsic stochasticity, or deterministic chaos) is closely related to the Kolmogorov-Arnold-Moser (KAM) theorem of the stability of near-integrable multiperiodic Hamiltonian systems. The effect of the intrinsic stochasticity attracts still growing attention both in theory and in various applications in contemporary physics. The paper discusses the relation of the intrinsic stochasticity to the modern ergodic theory and to the KAM theorem, and describes some numerical experiments on related astrophysical and high-temperature plasma problems. Some open questions are mentioned in conclusion. (author).

Stochastic Relational Presheaves and Dynamic Logic for Contextuality

Directory of Open Access Journals (Sweden)

Kohei Kishida

2014-12-01

Full Text Available Presheaf models provide a formulation of labelled transition systems that is useful for, among other things, modelling concurrent computation. This paper aims to extend such models further to represent stochastic dynamics such as shown in quantum systems. After reviewing what presheaf models represent and what certain operations on them mean in terms of notions such as internal and external choices, composition of systems, and so on, I will show how to extend those models and ideas by combining them with ideas from other category-theoretic approaches to relational models and to stochastic processes. It turns out that my extension yields a transitional formulation of sheaf-theoretic structures that Abramsky and Brandenburger proposed to characterize non-locality and contextuality. An alternative characterization of contextuality will then be given in terms of a dynamic modal logic of the models I put forward.

Measures of thermodynamic irreversibility in deterministic and stochastic dynamics

International Nuclear Information System (INIS)

Ford, Ian J

2015-01-01

It is generally observed that if a dynamical system is sufficiently complex, then as time progresses it will share out energy and other properties amongst its component parts to eliminate any initial imbalances, retaining only fluctuations. This is known as energy dissipation and it is closely associated with the concept of thermodynamic irreversibility, measured by the increase in entropy according to the second law. It is of interest to quantify such behaviour from a dynamical rather than a thermodynamic perspective and to this end stochastic entropy production and the time-integrated dissipation function have been introduced as analogous measures of irreversibility, principally for stochastic and deterministic dynamics, respectively. We seek to compare these measures. First we modify the dissipation function to allow it to measure irreversibility in situations where the initial probability density function (pdf) of the system is asymmetric as well as symmetric in velocity. We propose that it tests for failure of what we call the obversibility of the system, to be contrasted with reversibility, the failure of which is assessed by stochastic entropy production. We note that the essential difference between stochastic entropy production and the time-integrated modified dissipation function lies in the sequence of procedures undertaken in the associated tests of irreversibility. We argue that an assumed symmetry of the initial pdf with respect to velocity inversion (within a framework of deterministic dynamics) can be incompatible with the Past Hypothesis, according to which there should be a statistical distinction between the behaviour of certain properties of an isolated system as it evolves into the far future and the remote past. Imposing symmetry on a velocity distribution is acceptable for many applications of statistical physics, but can introduce difficulties when discussing irreversible behaviour. (paper)

Fast stochastic algorithm for simulating evolutionary population dynamics

Science.gov (United States)

Tsimring, Lev; Hasty, Jeff; Mather, William

2012-02-01

Evolution and co-evolution of ecological communities are stochastic processes often characterized by vastly different rates of reproduction and mutation and a coexistence of very large and very small sub-populations of co-evolving species. This creates serious difficulties for accurate statistical modeling of evolutionary dynamics. In this talk, we introduce a new exact algorithm for fast fully stochastic simulations of birth/death/mutation processes. It produces a significant speedup compared to the direct stochastic simulation algorithm in a typical case when the total population size is large and the mutation rates are much smaller than birth/death rates. We illustrate the performance of the algorithm on several representative examples: evolution on a smooth fitness landscape, NK model, and stochastic predator-prey system.

Dynamic and stochastic multi-project planning

CERN Document Server

Melchiors, Philipp

2015-01-01

This book deals with dynamic and stochastic methods for multi-project planning. Based on the idea of using queueing networks for the analysis of dynamic-stochastic multi-project environments this book addresses two problems: detailed scheduling of project activities, and integrated order acceptance and capacity planning. In an extensive simulation study, the book thoroughly investigates existing scheduling policies. To obtain optimal and near optimal scheduling policies new models and algorithms are proposed based on the theory of Markov decision processes and Approximate Dynamic programming.

Stochastic inequalities and applications to dynamics analysis of a novel SIVS epidemic model with jumps

Directory of Open Access Journals (Sweden)

Xiaona Leng

2017-06-01

Full Text Available Abstract This paper proposes a new nonlinear stochastic SIVS epidemic model with double epidemic hypothesis and Lévy jumps. The main purpose of this paper is to investigate the threshold dynamics of the stochastic SIVS epidemic model. By using the technique of a series of stochastic inequalities, we obtain sufficient conditions for the persistence in mean and extinction of the stochastic system and the threshold which governs the extinction and the spread of the epidemic diseases. Finally, this paper describes the results of numerical simulations investigating the dynamical effects of stochastic disturbance. Our results significantly improve and generalize the corresponding results in recent literatures. The developed theoretical methods and stochastic inequalities technique can be used to investigate the high-dimensional nonlinear stochastic differential systems.

Stochastic quasi-gradient based optimization algorithms for dynamic reliability applications

International Nuclear Information System (INIS)

Bourgeois, F.; Labeau, P.E.

2001-01-01

On one hand, PSA results are increasingly used in decision making, system management and optimization of system design. On the other hand, when severe accidental transients are considered, dynamic reliability appears appropriate to account for the complex interaction between the transitions between hardware configurations, the operator behavior and the dynamic evolution of the system. This paper presents an exploratory work in which the estimation of the system unreliability in a dynamic context is coupled with an optimization algorithm to determine the 'best' safety policy. Because some reliability parameters are likely to be distributed, the cost function to be minimized turns out to be a random variable. Stochastic programming techniques are therefore envisioned to determine an optimal strategy. Monte Carlo simulation is used at all stages of the computations, from the estimation of the system unreliability to that of the stochastic quasi-gradient. The optimization algorithm is illustrated on a HNO 3 supply system

Effects of stochastic noise on dynamical decoupling procedures

Energy Technology Data Exchange (ETDEWEB)

Bernad, Jozsef Zsolt; Frydrych, Holger; Alber, Gernot [Institut fuer Angewandte Physik, Technische Universitaet Darmstadt, D-64289 Darmstadt (Germany)

2013-07-01

Dynamical decoupling is a well-established technique to protect quantum systems from unwanted influences of their environment by exercising active control. It has been used experimentally to drastically increase the lifetime of qubit states in various implementations. The efficiency of different dynamical decoupling schemes defines the lifetime. However, errors in control operations always limit this efficiency. We propose a stochastic model as a possible description of imperfect control pulses and discuss the impact of this kind of error on different decoupling schemes. In the limit of continuous control, i.e. if the number of pulses N → ∞, we derive a stochastic differential equation for the evolution of the density operator of the controlled system and its environment. In the context of this modified time evolution we discuss possibilities of protecting qubit states against environmental noise.

Nambu mechanics for stochastic magnetization dynamics

Energy Technology Data Exchange (ETDEWEB)

Thibaudeau, Pascal, E-mail: pascal.thibaudeau@cea.fr [CEA DAM/Le Ripault, BP 16, F-37260 Monts (France); Nussle, Thomas, E-mail: thomas.nussle@cea.fr [CEA DAM/Le Ripault, BP 16, F-37260 Monts (France); CNRS-Laboratoire de Mathématiques et Physique Théorique (UMR 7350), Fédération de Recherche “Denis Poisson” (FR2964), Département de Physique, Université de Tours, Parc de Grandmont, F-37200 Tours (France); Nicolis, Stam, E-mail: stam.nicolis@lmpt.univ-tours.fr [CNRS-Laboratoire de Mathématiques et Physique Théorique (UMR 7350), Fédération de Recherche “Denis Poisson” (FR2964), Département de Physique, Université de Tours, Parc de Grandmont, F-37200 Tours (France)

2017-06-15

Highlights: • The LLG equation can be formulated in the framework of dissipative Nambu mechanics. • A master equation is derived for the spin dynamics for additive/multiplicative noises. • The derived stochastic equations are compared to moment equations obtained by closures. - Abstract: The Landau–Lifshitz–Gilbert (LLG) equation describes the dynamics of a damped magnetization vector that can be understood as a generalization of Larmor spin precession. The LLG equation cannot be deduced from the Hamiltonian framework, by introducing a coupling to a usual bath, but requires the introduction of additional constraints. It is shown that these constraints can be formulated elegantly and consistently in the framework of dissipative Nambu mechanics. This has many consequences for both the variational principle and for topological aspects of hidden symmetries that control conserved quantities. We particularly study how the damping terms of dissipative Nambu mechanics affect the consistent interaction of magnetic systems with stochastic reservoirs and derive a master equation for the magnetization. The proposals are supported by numerical studies using symplectic integrators that preserve the topological structure of Nambu equations. These results are compared to computations performed by direct sampling of the stochastic equations and by using closure assumptions for the moment equations, deduced from the master equation.

Dynamics of a Stochastic Intraguild Predation Model

Directory of Open Access Journals (Sweden)

Zejing Xing

2016-04-01

Full Text Available Intraguild predation (IGP is a widespread ecological phenomenon which occurs when one predator species attacks another predator species with which it competes for a shared prey species. The objective of this paper is to study the dynamical properties of a stochastic intraguild predation model. We analyze stochastic persistence and extinction of the stochastic IGP model containing five cases and establish the sufficient criteria for global asymptotic stability of the positive solutions. This study shows that it is possible for the coexistence of three species under the influence of environmental noise, and that the noise may have a positive effect for IGP species. A stationary distribution of the stochastic IGP model is established and it has the ergodic property, suggesting that the time average of population size with the development of time is equal to the stationary distribution in space. Finally, we show that our results may be extended to two well-known biological systems: food chains and exploitative competition.

The Theory of Dynamic Public Transit Priority with Dynamic Stochastic Park and Ride

Directory of Open Access Journals (Sweden)

Chengming Zhu

2014-01-01

Full Text Available Public transit priority is very important for relieving traffic congestion. The connotation of dynamic public transit priority and dynamic stochastic park and ride is presented. Based on the point that the travel cost of public transit is not higher than the travel cost of car, how to determine the level of dynamic public transit priority is discussed. The traffic organization method of dynamic public transit priority is introduced. For dynamic stochastic park and ride, layout principle, scale, and charging standard are discussed. Traveler acceptability is high through the analysis of questionnaire survey. Dynamic public transit priority with dynamic stochastic park and ride has application feasibility.

Stochastic dynamics of new inflation

International Nuclear Information System (INIS)

Nakao, Ken-ichi; Nambu, Yasusada; Sasaki, Misao.

1988-07-01

We investigate thoroughly the dynamics of an inflation-driving scalar field in terms of an extended version of the stochastic approach proposed by Starobinsky and discuss the spacetime structure of the inflationary universe. To avoid any complications which might arise due to quantum gravity, we concentrate our discussions on the new inflationary universe scenario in which all the energy scales involved are well below the planck mass. The investigation is done both analytically and numerically. In particular, we present a full numerical analysis of the stochastic scalar field dynamics on the phase space. Then implications of the results are discussed. (author)

Setting development goals using stochastic dynamical system models.

Science.gov (United States)

Ranganathan, Shyam; Nicolis, Stamatios C; Bali Swain, Ranjula; Sumpter, David J T

2017-01-01

The Millennium Development Goals (MDG) programme was an ambitious attempt to encourage a globalised solution to important but often-overlooked development problems. The programme led to wide-ranging development but it has also been criticised for unrealistic and arbitrary targets. In this paper, we show how country-specific development targets can be set using stochastic, dynamical system models built from historical data. In particular, we show that the MDG target of two-thirds reduction of child mortality from 1990 levels was infeasible for most countries, especially in sub-Saharan Africa. At the same time, the MDG targets were not ambitious enough for fast-developing countries such as Brazil and China. We suggest that model-based setting of country-specific targets is essential for the success of global development programmes such as the Sustainable Development Goals (SDG). This approach should provide clear, quantifiable targets for policymakers.

A data driven nonlinear stochastic model for blood glucose dynamics.

Science.gov (United States)

Zhang, Yan; Holt, Tim A; Khovanova, Natalia

2016-03-01

The development of adequate mathematical models for blood glucose dynamics may improve early diagnosis and control of diabetes mellitus (DM). We have developed a stochastic nonlinear second order differential equation to describe the response of blood glucose concentration to food intake using continuous glucose monitoring (CGM) data. A variational Bayesian learning scheme was applied to define the number and values of the system's parameters by iterative optimisation of free energy. The model has the minimal order and number of parameters to successfully describe blood glucose dynamics in people with and without DM. The model accounts for the nonlinearity and stochasticity of the underlying glucose-insulin dynamic process. Being data-driven, it takes full advantage of available CGM data and, at the same time, reflects the intrinsic characteristics of the glucose-insulin system without detailed knowledge of the physiological mechanisms. We have shown that the dynamics of some postprandial blood glucose excursions can be described by a reduced (linear) model, previously seen in the literature. A comprehensive analysis demonstrates that deterministic system parameters belong to different ranges for diabetes and controls. Implications for clinical practice are discussed. This is the first study introducing a continuous data-driven nonlinear stochastic model capable of describing both DM and non-DM profiles. Copyright © 2015 The Authors. Published by Elsevier Ireland Ltd.. All rights reserved.

Compositional Modelling of Stochastic Hybrid Systems

NARCIS (Netherlands)

Strubbe, S.N.

2005-01-01

In this thesis we present a modelling framework for compositional modelling of stochastic hybrid systems. Hybrid systems consist of a combination of continuous and discrete dynamics. The state space of a hybrid system is hybrid in the sense that it consists of a continuous component and a discrete

Stochastic dynamic modeling of regular and slow earthquakes

Science.gov (United States)

Aso, N.; Ando, R.; Ide, S.

2017-12-01

Both regular and slow earthquakes are slip phenomena on plate boundaries and are simulated by a (quasi-)dynamic modeling [Liu and Rice, 2005]. In these numerical simulations, spatial heterogeneity is usually considered not only for explaining real physical properties but also for evaluating the stability of the calculations or the sensitivity of the results on the condition. However, even though we discretize the model space with small grids, heterogeneity at smaller scales than the grid size is not considered in the models with deterministic governing equations. To evaluate the effect of heterogeneity at the smaller scales we need to consider stochastic interactions between slip and stress in a dynamic modeling. Tidal stress is known to trigger or affect both regular and slow earthquakes [Yabe et al., 2015; Ide et al., 2016], and such an external force with fluctuation can also be considered as a stochastic external force. A healing process of faults may also be stochastic, so we introduce stochastic friction law. In the present study, we propose a stochastic dynamic model to explain both regular and slow earthquakes. We solve mode III problem, which corresponds to the rupture propagation along the strike direction. We use BIEM (boundary integral equation method) scheme to simulate slip evolution, but we add stochastic perturbations in the governing equations, which is usually written in a deterministic manner. As the simplest type of perturbations, we adopt Gaussian deviations in the formulation of the slip-stress kernel, external force, and friction. By increasing the amplitude of perturbations of the slip-stress kernel, we reproduce complicated rupture process of regular earthquakes including unilateral and bilateral ruptures. By perturbing external force, we reproduce slow rupture propagation at a scale of km/day. The slow propagation generated by a combination of fast interaction at S-wave velocity is analogous to the kinetic theory of gasses: thermal

Nonlinear and stochastic dynamics of coherent structures

DEFF Research Database (Denmark)

Rasmussen, Kim

1997-01-01

This Thesis deals with nonlinear and stochastic dynamics in systems which can be described by nonlinear Schrödinger models. Basically three different models are investigated. The first is the continuum nonlinear Schröndinger model in one and two dimensions generalized by a tunable degree of nonli......This Thesis deals with nonlinear and stochastic dynamics in systems which can be described by nonlinear Schrödinger models. Basically three different models are investigated. The first is the continuum nonlinear Schröndinger model in one and two dimensions generalized by a tunable degree...... introduces the nonlinear Schrödinger model in one and two dimensions, discussing the soliton solutions in one dimension and the collapse phenomenon in two dimensions. Also various analytical methods are described. Then a derivation of the nonlinear Schrödinger equation is given, based on a Davydov like...... system described by a tight-binding Hamiltonian and a harmonic lattice coupled b y a deformation-type potential. This derivation results in a two-dimensional nonline ar Schrödinger model, and considering the harmonic lattice to be in thermal contact with a heat bath w e show that the nonlinear...

Adaptive neural network output feedback control for stochastic nonlinear systems with unknown dead-zone and unmodeled dynamics.

Science.gov (United States)

Tong, Shaocheng; Wang, Tong; Li, Yongming; Zhang, Huaguang

2014-06-01

This paper discusses the problem of adaptive neural network output feedback control for a class of stochastic nonlinear strict-feedback systems. The concerned systems have certain characteristics, such as unknown nonlinear uncertainties, unknown dead-zones, unmodeled dynamics and without the direct measurements of state variables. In this paper, the neural networks (NNs) are employed to approximate the unknown nonlinear uncertainties, and then by representing the dead-zone as a time-varying system with a bounded disturbance. An NN state observer is designed to estimate the unmeasured states. Based on both backstepping design technique and a stochastic small-gain theorem, a robust adaptive NN output feedback control scheme is developed. It is proved that all the variables involved in the closed-loop system are input-state-practically stable in probability, and also have robustness to the unmodeled dynamics. Meanwhile, the observer errors and the output of the system can be regulated to a small neighborhood of the origin by selecting appropriate design parameters. Simulation examples are also provided to illustrate the effectiveness of the proposed approach.

Distributed Consensus of Stochastic Delayed Multi-agent Systems Under Asynchronous Switching.

Science.gov (United States)

Wu, Xiaotai; Tang, Yang; Cao, Jinde; Zhang, Wenbing

2016-08-01

In this paper, the distributed exponential consensus of stochastic delayed multi-agent systems with nonlinear dynamics is investigated under asynchronous switching. The asynchronous switching considered here is to account for the time of identifying the active modes of multi-agent systems. After receipt of confirmation of mode's switching, the matched controller can be applied, which means that the switching time of the matched controller in each node usually lags behind that of system switching. In order to handle the coexistence of switched signals and stochastic disturbances, a comparison principle of stochastic switched delayed systems is first proved. By means of this extended comparison principle, several easy to verified conditions for the existence of an asynchronously switched distributed controller are derived such that stochastic delayed multi-agent systems with asynchronous switching and nonlinear dynamics can achieve global exponential consensus. Two examples are given to illustrate the effectiveness of the proposed method.

Modeling and Properties of Nonlinear Stochastic Dynamical System of Continuous Culture

Science.gov (United States)

Wang, Lei; Feng, Enmin; Ye, Jianxiong; Xiu, Zhilong

The stochastic counterpart to the deterministic description of continuous fermentation with ordinary differential equation is investigated in the process of glycerol bio-dissimilation to 1,3-propanediol by Klebsiella pneumoniae. We briefly discuss the continuous fermentation process driven by three-dimensional Brownian motion and Lipschitz coefficients, which is suitable for the factual fermentation. Subsequently, we study the existence and uniqueness of solutions for the stochastic system as well as the boundedness of the Two-order Moment and the Markov property of the solution. Finally stochastic simulation is carried out under the Stochastic Euler-Maruyama method.

Computational Methods in Stochastic Dynamics Volume 2

CERN Document Server

Stefanou, George; Papadopoulos, Vissarion

2013-01-01

The considerable influence of inherent uncertainties on structural behavior has led the engineering community to recognize the importance of a stochastic approach to structural problems. Issues related to uncertainty quantification and its influence on the reliability of the computational models are continuously gaining in significance. In particular, the problems of dynamic response analysis and reliability assessment of structures with uncertain system and excitation parameters have been the subject of continuous research over the last two decades as a result of the increasing availability of powerful computing resources and technology.   This book is a follow up of a previous book with the same subject (ISBN 978-90-481-9986-0) and focuses on advanced computational methods and software tools which can highly assist in tackling complex problems in stochastic dynamic/seismic analysis and design of structures. The selected chapters are authored by some of the most active scholars in their respective areas and...

Systemic risk in dynamical networks with stochastic failure criterion

Science.gov (United States)

Podobnik, B.; Horvatic, D.; Bertella, M. A.; Feng, L.; Huang, X.; Li, B.

2014-06-01

Complex non-linear interactions between banks and assets we model by two time-dependent Erdős-Renyi network models where each node, representing a bank, can invest either to a single asset (model I) or multiple assets (model II). We use a dynamical network approach to evaluate the collective financial failure —systemic risk— quantified by the fraction of active nodes. The systemic risk can be calculated over any future time period, divided into sub-periods, where within each sub-period banks may contiguously fail due to links to either i) assets or ii) other banks, controlled by two parameters, probability of internal failure p and threshold Th (“solvency” parameter). The systemic risk decreases with the average network degree faster when all assets are equally distributed across banks than if assets are randomly distributed. The more inactive banks each bank can sustain (smaller Th), the smaller the systemic risk —for some Th values in I we report a discontinuity in systemic risk. When contiguous spreading becomes stochastic ii) controlled by probability p2 —a condition for the bank to be solvent (active) is stochastic— the systemic risk decreases with decreasing p2. We analyse the asset allocation for the U.S. banks.

Advanced models of neural networks nonlinear dynamics and stochasticity in biological neurons

CERN Document Server

Rigatos, Gerasimos G

2015-01-01

This book provides a complete study on neural structures exhibiting nonlinear and stochastic dynamics, elaborating on neural dynamics by introducing advanced models of neural networks. It overviews the main findings in the modelling of neural dynamics in terms of electrical circuits and examines their stability properties with the use of dynamical systems theory. It is suitable for researchers and postgraduate students engaged with neural networks and dynamical systems theory.

The Theory of Dynamic Public Transit Priority with Dynamic Stochastic Park and Ride

OpenAIRE

Zhu, Chengming; Chen, Yanyan; Ma, Changxi

2014-01-01

Public transit priority is very important for relieving traffic congestion. The connotation of dynamic public transit priority and dynamic stochastic park and ride is presented. Based on the point that the travel cost of public transit is not higher than the travel cost of car, how to determine the level of dynamic public transit priority is discussed. The traffic organization method of dynamic public transit priority is introduced. For dynamic stochastic park and ride, layout principle, scal...

Dynamic stochastic accumulation model with application to pension savings management

Directory of Open Access Journals (Sweden)

Melicherčik Igor

2010-01-01

Full Text Available We propose a dynamic stochastic accumulation model for determining optimal decision between stock and bond investments during accumulation of pension savings. Stock prices are assumed to be driven by the geometric Brownian motion. Interest rates are modeled by means of the Cox-Ingersoll-Ross model. The optimal decision as a solution to the corresponding dynamic stochastic program is a function of the duration of saving, the level of savings and the short rate. Qualitative and quantitative properties of the optimal solution are analyzed. The model is tested on the funded pillar of the Slovak pension system. The results are calculated for various risk preferences of a saver.

Agent-based financial dynamics model from stochastic interacting epidemic system and complexity analysis

International Nuclear Information System (INIS)

Lu, Yunfan; Wang, Jun; Niu, Hongli

2015-01-01

An agent-based financial stock price model is developed and investigated by a stochastic interacting epidemic system, which is one of the statistical physics systems and has been used to model the spread of an epidemic or a forest fire. Numerical and statistical analysis are performed on the simulated returns of the proposed financial model. Complexity properties of the financial time series are explored by calculating the correlation dimension and using the modified multiscale entropy method. In order to verify the rationality of the financial model, the real stock market indexes, Shanghai Composite Index and Shenzhen Component Index, are studied in comparison with the simulation data of the proposed model for the different infectiousness parameters. The empirical research reveals that this financial model can reproduce some important features of the real stock markets. - Highlights: • A new agent-based financial price model is developed by stochastic interacting epidemic system. • The structure of the proposed model allows to simulate the financial dynamics. • Correlation dimension and MMSE are applied to complexity analysis of financial time series. • Empirical results show the rationality of the proposed financial model

Agent-based financial dynamics model from stochastic interacting epidemic system and complexity analysis

Energy Technology Data Exchange (ETDEWEB)

Lu, Yunfan, E-mail: yunfanlu@yeah.net; Wang, Jun; Niu, Hongli

2015-06-12

An agent-based financial stock price model is developed and investigated by a stochastic interacting epidemic system, which is one of the statistical physics systems and has been used to model the spread of an epidemic or a forest fire. Numerical and statistical analysis are performed on the simulated returns of the proposed financial model. Complexity properties of the financial time series are explored by calculating the correlation dimension and using the modified multiscale entropy method. In order to verify the rationality of the financial model, the real stock market indexes, Shanghai Composite Index and Shenzhen Component Index, are studied in comparison with the simulation data of the proposed model for the different infectiousness parameters. The empirical research reveals that this financial model can reproduce some important features of the real stock markets. - Highlights: • A new agent-based financial price model is developed by stochastic interacting epidemic system. • The structure of the proposed model allows to simulate the financial dynamics. • Correlation dimension and MMSE are applied to complexity analysis of financial time series. • Empirical results show the rationality of the proposed financial model.

A Stochastic Model for Malaria Transmission Dynamics

Directory of Open Access Journals (Sweden)

Rachel Waema Mbogo

2018-01-01

Full Text Available Malaria is one of the three most dangerous infectious diseases worldwide (along with HIV/AIDS and tuberculosis. In this paper we compare the disease dynamics of the deterministic and stochastic models in order to determine the effect of randomness in malaria transmission dynamics. Relationships between the basic reproduction number for malaria transmission dynamics between humans and mosquitoes and the extinction thresholds of corresponding continuous-time Markov chain models are derived under certain assumptions. The stochastic model is formulated using the continuous-time discrete state Galton-Watson branching process (CTDSGWbp. The reproduction number of deterministic models is an essential quantity to predict whether an epidemic will spread or die out. Thresholds for disease extinction from stochastic models contribute crucial knowledge on disease control and elimination and mitigation of infectious diseases. Analytical and numerical results show some significant differences in model predictions between the stochastic and deterministic models. In particular, we find that malaria outbreak is more likely if the disease is introduced by infected mosquitoes as opposed to infected humans. These insights demonstrate the importance of a policy or intervention focusing on controlling the infected mosquito population if the control of malaria is to be realized.

Decentralized adaptive neural control for high-order interconnected stochastic nonlinear time-delay systems with unknown system dynamics.

Science.gov (United States)

Si, Wenjie; Dong, Xunde; Yang, Feifei

2018-03-01

This paper is concerned with the problem of decentralized adaptive backstepping state-feedback control for uncertain high-order large-scale stochastic nonlinear time-delay systems. For the control design of high-order large-scale nonlinear systems, only one adaptive parameter is constructed to overcome the over-parameterization, and neural networks are employed to cope with the difficulties raised by completely unknown system dynamics and stochastic disturbances. And then, the appropriate Lyapunov-Krasovskii functional and the property of hyperbolic tangent functions are used to deal with the unknown unmatched time-delay interactions of high-order large-scale systems for the first time. At last, on the basis of Lyapunov stability theory, the decentralized adaptive neural controller was developed, and it decreases the number of learning parameters. The actual controller can be designed so as to ensure that all the signals in the closed-loop system are semi-globally uniformly ultimately bounded (SGUUB) and the tracking error converges in the small neighborhood of zero. The simulation example is used to further show the validity of the design method. Copyright © 2018 Elsevier Ltd. All rights reserved.

Probabilistic DHP adaptive critic for nonlinear stochastic control systems.

Science.gov (United States)

Herzallah, Randa

2013-06-01

Following the recently developed algorithms for fully probabilistic control design for general dynamic stochastic systems (Herzallah & Káarnáy, 2011; Kárný, 1996), this paper presents the solution to the probabilistic dual heuristic programming (DHP) adaptive critic method (Herzallah & Káarnáy, 2011) and randomized control algorithm for stochastic nonlinear dynamical systems. The purpose of the randomized control input design is to make the joint probability density function of the closed loop system as close as possible to a predetermined ideal joint probability density function. This paper completes the previous work (Herzallah & Káarnáy, 2011; Kárný, 1996) by formulating and solving the fully probabilistic control design problem on the more general case of nonlinear stochastic discrete time systems. A simulated example is used to demonstrate the use of the algorithm and encouraging results have been obtained. Copyright © 2013 Elsevier Ltd. All rights reserved.

A Volterra series approach to the approximation of stochastic nonlinear dynamics

NARCIS (Netherlands)

Wouw, van de N.; Nijmeijer, H.; Campen, van D.H.

2002-01-01

A response approximation method for stochastically excited, nonlinear, dynamic systems is presented. Herein, the output of the nonlinear system isapproximated by a finite-order Volterra series. The original nonlinear system is replaced by a bilinear system in order to determine the kernels of this

Path to Stochastic Stability: Comparative Analysis of Stochastic Learning Dynamics in Games

KAUST Repository

Jaleel, Hassan; Shamma, Jeff S.

2018-01-01

dynamics: Log-Linear Learning (LLL) and Metropolis Learning (ML). Although both of these dynamics have the same stochastically stable states, LLL and ML correspond to different behavioral models for decision making. Moreover, we demonstrate through

Quantization of dynamical systems and stochastic control theory

International Nuclear Information System (INIS)

Guerra, F.; Morato, L.M.

1982-09-01

In the general framework of stochastic control theory we introduce a suitable form of stochastic action associated to the controlled process. Then a variational principle gives all main features of Nelson's stochastic mechanics. In particular we derive the expression of the current velocity field as the gradient of the phase action. Moreover the stochastic corrections to the Hamilton-Jacobi equation are in agreement with the quantum mechanical form of the Madelung fluid (equivalent to the Schroedinger equation). Therefore stochastic control theory can provide a very simple model simulating quantum mechanical behavior

Dynamic optimization deterministic and stochastic models

CERN Document Server

Hinderer, Karl; Stieglitz, Michael

2016-01-01

This book explores discrete-time dynamic optimization and provides a detailed introduction to both deterministic and stochastic models. Covering problems with finite and infinite horizon, as well as Markov renewal programs, Bayesian control models and partially observable processes, the book focuses on the precise modelling of applications in a variety of areas, including operations research, computer science, mathematics, statistics, engineering, economics and finance. Dynamic Optimization is a carefully presented textbook which starts with discrete-time deterministic dynamic optimization problems, providing readers with the tools for sequential decision-making, before proceeding to the more complicated stochastic models. The authors present complete and simple proofs and illustrate the main results with numerous examples and exercises (without solutions). With relevant material covered in four appendices, this book is completely self-contained.

Size and stochasticity in irrigated social-ecological systems

Science.gov (United States)

Puy, Arnald; Muneepeerakul, Rachata; Balbo, Andrea L.

2017-03-01

This paper presents a systematic study of the relation between the size of irrigation systems and the management of uncertainty. We specifically focus on studying, through a stylized theoretical model, how stochasticity in water availability and taxation interacts with the stochastic behavior of the population within irrigation systems. Our results indicate the existence of two key population thresholds for the sustainability of any irrigation system: or the critical population size required to keep the irrigation system operative, and N* or the population threshold at which the incentive to work inside the irrigation system equals the incentives to work elsewhere. Crossing irretrievably leads to system collapse. N* is the population level with a sub-optimal per capita payoff towards which irrigation systems tend to gravitate. When subjected to strong stochasticity in water availability or taxation, irrigation systems might suffer sharp population drops and irreversibly disintegrate into a system collapse, via a mechanism we dub ‘collapse trap’. Our conceptual study establishes the basis for further work aiming at appraising the dynamics between size and stochasticity in irrigation systems, whose understanding is key for devising mitigation and adaptation measures to ensure their sustainability in the face of increasing and inevitable uncertainty.

Sharp asymptotics for stochastic dynamics with parallel updating rule

NARCIS (Netherlands)

Nardi, F.R.; Spitoni, C.

2012-01-01

In this paper we study the metastability problem for a stochastic dynamics with a parallel updating rule; in particular we consider a finite volume Probabilistic Cellular Automaton (PCA) in a small external field at low temperature regime. We are interested in the nucleation of the system, i.e., the

Sharp Asymptotics for Stochastic Dynamics with Parallel Updating Rule

NARCIS (Netherlands)

Nardi, F.R.; Spitoni, C.

2012-01-01

In this paper we study the metastability problem for a stochastic dynamics with a parallel updating rule; in particular we consider a finite volume Probabilistic Cellular Automaton (PCA) in a small external field at low temperature regime. We are interested in the nucleation of the system, i.e.,

Hopf Bifurcation Analysis for a Stochastic Discrete-Time Hyperchaotic System

Directory of Open Access Journals (Sweden)

Jie Ran

2015-01-01

Full Text Available The dynamics of a discrete-time hyperchaotic system and the amplitude control of Hopf bifurcation for a stochastic discrete-time hyperchaotic system are investigated in this paper. Numerical simulations are presented to exhibit the complex dynamical behaviors in the discrete-time hyperchaotic system. Furthermore, the stochastic discrete-time hyperchaotic system with random parameters is transformed into its equivalent deterministic system with the orthogonal polynomial theory of discrete random function. In addition, the dynamical features of the discrete-time hyperchaotic system with random disturbances are obtained through its equivalent deterministic system. By using the Hopf bifurcation conditions of the deterministic discrete-time system, the specific conditions for the existence of Hopf bifurcation in the equivalent deterministic system are derived. And the amplitude control with random intensity is discussed in detail. Finally, the feasibility of the control method is demonstrated by numerical simulations.

On an aggregation in birth-and-death stochastic dynamics

Science.gov (United States)

Finkelshtein, Dmitri; Kondratiev, Yuri; Kutoviy, Oleksandr; Zhizhina, Elena

2014-06-01

We consider birth-and-death stochastic dynamics of particle systems with attractive interaction. The heuristic generator of the dynamics has a constant birth rate and density-dependent decreasing death rate. The corresponding statistical dynamics is constructed. Using the Vlasov-type scaling we derive the limiting mesoscopic evolution and prove that this evolution propagates chaos. We study a nonlinear non-local kinetic equation for the first correlation function (density of population). The existence of uniformly bounded solutions as well as solutions growing inside of a bounded domain and expanding in the space are shown. These solutions describe two regimes in the mesoscopic system: regulation and aggregation.

On an aggregation in birth-and-death stochastic dynamics

International Nuclear Information System (INIS)

Finkelshtein, Dmitri; Kondratiev, Yuri; Kutoviy, Oleksandr; Zhizhina, Elena

2014-01-01

We consider birth-and-death stochastic dynamics of particle systems with attractive interaction. The heuristic generator of the dynamics has a constant birth rate and density-dependent decreasing death rate. The corresponding statistical dynamics is constructed. Using the Vlasov-type scaling we derive the limiting mesoscopic evolution and prove that this evolution propagates chaos. We study a nonlinear non-local kinetic equation for the first correlation function (density of population). The existence of uniformly bounded solutions as well as solutions growing inside of a bounded domain and expanding in the space are shown. These solutions describe two regimes in the mesoscopic system: regulation and aggregation. (paper)

Computational singular perturbation analysis of stochastic chemical systems with stiffness

Science.gov (United States)

Wang, Lijin; Han, Xiaoying; Cao, Yanzhao; Najm, Habib N.

2017-04-01

Computational singular perturbation (CSP) is a useful method for analysis, reduction, and time integration of stiff ordinary differential equation systems. It has found dominant utility, in particular, in chemical reaction systems with a large range of time scales at continuum and deterministic level. On the other hand, CSP is not directly applicable to chemical reaction systems at micro or meso-scale, where stochasticity plays an non-negligible role and thus has to be taken into account. In this work we develop a novel stochastic computational singular perturbation (SCSP) analysis and time integration framework, and associated algorithm, that can be used to not only construct accurately and efficiently the numerical solutions to stiff stochastic chemical reaction systems, but also analyze the dynamics of the reduced stochastic reaction systems. The algorithm is illustrated by an application to a benchmark stochastic differential equation model, and numerical experiments are carried out to demonstrate the effectiveness of the construction.

Sharp asymptotics for stochastic dynamics with parallel updating rule

NARCIS (Netherlands)

Nardi, F.R.; Spitoni, C.

2012-01-01

In this paper we study the metastability problem for a stochastic dynamics with a parallel updating rule; in particular we consider a ¿nite volume Probabilistic Cellular Automaton (PCA) in a small external ¿eld at low temperature regime. We are interested in the nucleation of the system, i.e., the

Stochastic Dynamics through Hierarchically Embedded Markov Chains.

Science.gov (United States)

Vasconcelos, Vítor V; Santos, Fernando P; Santos, Francisco C; Pacheco, Jorge M

2017-02-03

Studying dynamical phenomena in finite populations often involves Markov processes of significant mathematical and/or computational complexity, which rapidly becomes prohibitive with increasing population size or an increasing number of individual configuration states. Here, we develop a framework that allows us to define a hierarchy of approximations to the stationary distribution of general systems that can be described as discrete Markov processes with time invariant transition probabilities and (possibly) a large number of states. This results in an efficient method for studying social and biological communities in the presence of stochastic effects-such as mutations in evolutionary dynamics and a random exploration of choices in social systems-including situations where the dynamics encompasses the existence of stable polymorphic configurations, thus overcoming the limitations of existing methods. The present formalism is shown to be general in scope, widely applicable, and of relevance to a variety of interdisciplinary problems.

Stochastic dynamics of complex systems: from glasses to evolution (series on complexity science)

CERN Document Server

Sibani, Paolo

2013-01-01

Dynamical evolution over long time scales is a prominent feature of all the systems we intuitively think of as complex - for example, ecosystems, the brain or the economy. In physics, the term ageing is used for this type of slow change, occurring over time scales much longer than the patience, or indeed the lifetime, of the observer. The main focus of this book is on the stochastic processes which cause ageing, and the surprising fact that the ageing dynamics of systems which are very different at the microscopic level can be treated in similar ways. The first part of this book provides the necessary mathematical and computational tools and the second part describes the intuition needed to deal with these systems. Some of the first few chapters have been covered in several other books, but the emphasis and selection of the topics reflect both the authors' interests and the overall theme of the book. The second part contains an introduction to the scientific literature and deals in some detail with the desc...

Sustainability from the Occurrence of Critical Dynamic Power System Blackout Determined by Using the Stochastic Event Tree Technique

Directory of Open Access Journals (Sweden)

Muhammad Murtadha Othman

2017-06-01

Full Text Available With the advent of advanced technology in smart grid, the implementation of renewable energy in a stressed and complicated power system operation, aggravated by a competitive electricity market and critical system contingencies, this will inflict higher probabilities of the occurrence of a severe dynamic power system blackout. This paper presents the proposed stochastic event tree technique used to assess the sustainability against the occurrence of dynamic power system blackout emanating from implication of critical system contingencies such as the rapid increase in total loading condition and sensitive initial transmission line tripping. An extensive analysis of dynamic power system blackout has been carried out in a case study of the following power systems: IEEE RTS-79 and IEEE RTS-96. The findings have shown that the total loading conditions and sensitive transmission lines need to be given full attention by the utility to prevent the occurrence of dynamic power system blackout.

Stochastic linearization of turbulent dynamics of dispersive waves in equilibrium and non-equilibrium state

International Nuclear Information System (INIS)

Jiang, Shixiao W; Lu, Haihao; Zhou, Douglas; Cai, David

2016-01-01

Characterizing dispersive wave turbulence in the long time dynamics is central to understanding of many natural phenomena, e.g., in atmosphere ocean dynamics, nonlinear optics, and plasma physics. Using the β -Fermi–Pasta–Ulam nonlinear system as a prototypical example, we show that in thermal equilibrium and non-equilibrium steady state the turbulent state even in the strongly nonlinear regime possesses an effective linear stochastic structure in renormalized normal variables. In this framework, we can well characterize the spatiotemporal dynamics, which are dominated by long-wavelength renormalized waves. We further demonstrate that the energy flux is nearly saturated by the long-wavelength renormalized waves in non-equilibrium steady state. The scenario of such effective linear stochastic dynamics can be extended to study turbulent states in other nonlinear wave systems. (paper)

Stochastic Optimisation of Battery System Operation Strategy under different Utility Tariff Structures

OpenAIRE

Erdal, Jørgen Sørgård

2017-01-01

This master thesis develops a stochastic optimisation software for household grid-connected batteries combined with PV-systems. The objective of the optimisation is to operate the battery system in order to minimise the costs of the consumer, and it was implemented in MATLAB using a self-written stochastic dynamic programming algorithm. Load was considered as a stochastic variable and modelled as a Markov Chain. Transition probabilities between time steps were calculated using historic load p...

Dynamics and Physiological Roles of Stochastic Firing Patterns Near Bifurcation Points

Science.gov (United States)

Jia, Bing; Gu, Huaguang

2017-06-01

Different stochastic neural firing patterns or rhythms that appeared near polarization or depolarization resting states were observed in biological experiments on three nervous systems, and closely matched those simulated near bifurcation points between stable equilibrium point and limit cycle in a theoretical model with noise. The distinct dynamics of spike trains and interspike interval histogram (ISIH) of these stochastic rhythms were identified and found to build a relationship to the coexisting behaviors or fixed firing frequency of four different types of bifurcations. Furthermore, noise evokes coherence resonances near bifurcation points and plays important roles in enhancing information. The stochastic rhythms corresponding to Hopf bifurcation points with fixed firing frequency exhibited stronger coherence degree and a sharper peak in the power spectrum of the spike trains than those corresponding to saddle-node bifurcation points without fixed firing frequency. Moreover, the stochastic firing patterns changed to a depolarization resting state as the extracellular potassium concentration increased for the injured nerve fiber related to pathological pain or static blood pressure level increased for aortic depressor nerve fiber, and firing frequency decreased, which were different from the physiological viewpoint that firing frequency increased with increasing pressure level or potassium concentration. This shows that rhythms or firing patterns can reflect pressure or ion concentration information related to pathological pain information. Our results present the dynamics of stochastic firing patterns near bifurcation points, which are helpful for the identification of both dynamics and physiological roles of complex neural firing patterns or rhythms, and the roles of noise.

On the stochastic dynamics of disordered spin models

International Nuclear Information System (INIS)

Semerjian, G.; Montanari, A.; Cugliandolo, L.F.

2003-09-01

In this article we discuss several aspects of the stochastic dynamics of spin models. The paper has two independent parts. Firstly, we explore a few properties of the multi-point correlations and responses of generic systems evolving in equilibrium with a thermal bath. We propose a fluctuation principle that allows us to derive fluctuation-dissipation relations for many-time correlations and linear responses. We also speculate on how these features will be modified in systems evolving slowly out of equilibrium, as finite-dimensional or dilute spin-glasses. Secondly, we present a formalism that allows one to derive a series of approximated equations that determine the dynamics of disordered spin models on random (hyper) graphs. (author)

Simulation of quantum dynamics based on the quantum stochastic differential equation.

Science.gov (United States)

Li, Ming

2013-01-01

The quantum stochastic differential equation derived from the Lindblad form quantum master equation is investigated. The general formulation in terms of environment operators representing the quantum state diffusion is given. The numerical simulation algorithm of stochastic process of direct photodetection of a driven two-level system for the predictions of the dynamical behavior is proposed. The effectiveness and superiority of the algorithm are verified by the performance analysis of the accuracy and the computational cost in comparison with the classical Runge-Kutta algorithm.

Analysis of dispatching rules in a stochastic dynamic job shop manufacturing system with sequence-dependent setup times

Science.gov (United States)

Sharma, Pankaj; Jain, Ajai

2014-12-01

Stochastic dynamic job shop scheduling problem with consideration of sequence-dependent setup times are among the most difficult classes of scheduling problems. This paper assesses the performance of nine dispatching rules in such shop from makespan, mean flow time, maximum flow time, mean tardiness, maximum tardiness, number of tardy jobs, total setups and mean setup time performance measures viewpoint. A discrete event simulation model of a stochastic dynamic job shop manufacturing system is developed for investigation purpose. Nine dispatching rules identified from literature are incorporated in the simulation model. The simulation experiments are conducted under due date tightness factor of 3, shop utilization percentage of 90% and setup times less than processing times. Results indicate that shortest setup time (SIMSET) rule provides the best performance for mean flow time and number of tardy jobs measures. The job with similar setup and modified earliest due date (JMEDD) rule provides the best performance for makespan, maximum flow time, mean tardiness, maximum tardiness, total setups and mean setup time measures.

Parallel replica dynamics method for bistable stochastic reaction networks: Simulation and sensitivity analysis

Science.gov (United States)

Wang, Ting; Plecháč, Petr

2017-12-01

Stochastic reaction networks that exhibit bistable behavior are common in systems biology, materials science, and catalysis. Sampling of stationary distributions is crucial for understanding and characterizing the long-time dynamics of bistable stochastic dynamical systems. However, simulations are often hindered by the insufficient sampling of rare transitions between the two metastable regions. In this paper, we apply the parallel replica method for a continuous time Markov chain in order to improve sampling of the stationary distribution in bistable stochastic reaction networks. The proposed method uses parallel computing to accelerate the sampling of rare transitions. Furthermore, it can be combined with the path-space information bounds for parametric sensitivity analysis. With the proposed methodology, we study three bistable biological networks: the Schlögl model, the genetic switch network, and the enzymatic futile cycle network. We demonstrate the algorithmic speedup achieved in these numerical benchmarks. More significant acceleration is expected when multi-core or graphics processing unit computer architectures and programming tools such as CUDA are employed.

Parallel replica dynamics method for bistable stochastic reaction networks: Simulation and sensitivity analysis.

Science.gov (United States)

Wang, Ting; Plecháč, Petr

2017-12-21

Stochastic reaction networks that exhibit bistable behavior are common in systems biology, materials science, and catalysis. Sampling of stationary distributions is crucial for understanding and characterizing the long-time dynamics of bistable stochastic dynamical systems. However, simulations are often hindered by the insufficient sampling of rare transitions between the two metastable regions. In this paper, we apply the parallel replica method for a continuous time Markov chain in order to improve sampling of the stationary distribution in bistable stochastic reaction networks. The proposed method uses parallel computing to accelerate the sampling of rare transitions. Furthermore, it can be combined with the path-space information bounds for parametric sensitivity analysis. With the proposed methodology, we study three bistable biological networks: the Schlögl model, the genetic switch network, and the enzymatic futile cycle network. We demonstrate the algorithmic speedup achieved in these numerical benchmarks. More significant acceleration is expected when multi-core or graphics processing unit computer architectures and programming tools such as CUDA are employed.

Nonlinear dynamics and bifurcation characteristics of shape memory alloy thin films subjected to in-plane stochastic excitation

International Nuclear Information System (INIS)

Zhu, Zhi-Wen; Zhang, Qing-Xin; Xu, Jia

2014-01-01

A kind of shape memory alloy (SMA) hysteretic nonlinear model was developed, and the nonlinear dynamics and bifurcation characteristics of the SMA thin film subjected to in-plane stochastic excitation were investigated. Van der Pol difference item was introduced to describe the hysteretic phenomena of the SMA strain–stress curves, and the nonlinear dynamic model of the SMA thin film subjected to in-plane stochastic excitation was developed. The conditions of global stochastic stability of the system were determined in singular boundary theory, and the probability density function of the system response was obtained. Finally, the conditions of stochastic Hopf bifurcation were analyzed. The results of theoretical analysis and numerical simulation indicate that self-excited vibration is induced by the hysteretic nonlinear characteristics of SMA, and stochastic Hopf bifurcation appears when the bifurcation parameter was changed; there are two limit cycles in the stationary probability density of the dynamic response of the system in some cases, which means that there are two vibration amplitudes whose probabilities are both very high, and jumping phenomena between the two vibration amplitudes appear with the change in conditions. The results obtained in this current paper are helpful for the application of the SMA thin film in stochastic vibration fields. - Highlights: • Hysteretic nonlinear model of shape memory alloy was developed. • Van der Pol item was introduced to interpret hysteretic strain–stress curves. • Nonlinear dynamic characteristics of the shape memory alloy film were analyzed. • Jumping phenomena were observed in the change of the parameters

Patterns of Stochastic Behavior in Dynamically Unstable High-Dimensional Biochemical Networks

Directory of Open Access Journals (Sweden)

Simon Rosenfeld

2009-01-01

Full Text Available The question of dynamical stability and stochastic behavior of large biochemical networks is discussed. It is argued that stringent conditions of asymptotic stability have very little chance to materialize in a multidimensional system described by the differential equations of chemical kinetics. The reason is that the criteria of asymptotic stability (Routh- Hurwitz, Lyapunov criteria, Feinberg’s Deficiency Zero theorem would impose the limitations of very high algebraic order on the kinetic rates and stoichiometric coefficients, and there are no natural laws that would guarantee their unconditional validity. Highly nonlinear, dynamically unstable systems, however, are not necessarily doomed to collapse, as a simple Jacobian analysis would suggest. It is possible that their dynamics may assume the form of pseudo-random fluctuations quite similar to a shot noise, and, therefore, their behavior may be described in terms of Langevin and Fokker-Plank equations. We have shown by simulation that the resulting pseudo-stochastic processes obey the heavy-tailed Generalized Pareto Distribution with temporal sequence of pulses forming the set of constituent-specific Poisson processes. Being applied to intracellular dynamics, these properties are naturally associated with burstiness, a well documented phenomenon in the biology of gene expression.

Inter-species competition-facilitation in stochastic riparian vegetation dynamics.

Science.gov (United States)

Tealdi, Stefano; Camporeale, Carlo; Ridolfi, Luca

2013-02-07

Riparian vegetation is a highly dynamic community that lives on river banks and which depends to a great extent on the fluvial hydrology. The stochasticity of the discharge and erosion/deposition processes in fact play a key role in determining the distribution of vegetation along a riparian transect. These abiotic processes interact with biotic competition/facilitation mechanisms, such as plant competition for light, water, and nutrients. In this work, we focus on the dynamics of plants characterized by three components: (1) stochastic forcing due to river discharges, (2) competition for resources, and (3) inter-species facilitation due to the interplay between vegetation and fluid dynamics processes. A minimalist stochastic bio-hydrological model is proposed for the dynamics of the biomass of two vegetation species: one species is assumed dominant and slow-growing, the other is subdominant, but fast-growing. The stochastic model is solved analytically and the probability density function of the plant biomasses is obtained as a function of both the hydrologic and biologic parameters. The impact of the competition/facilitation processes on the distribution of vegetation species along the riparian transect is investigated and remarkable effects are observed. Finally, a good qualitative agreement is found between the model results and field data. Copyright © 2012 Elsevier Ltd. All rights reserved.

Outbreak and Extinction Dynamics in a Stochastic Ebola Model

Science.gov (United States)

Nieddu, Garrett; Bianco, Simone; Billings, Lora; Forgoston, Eric; Kaufman, James

A zoonotic disease is a disease that can be passed between animals and humans. In many cases zoonotic diseases can persist in the animal population even if there are no infections in the human population. In this case we call the infected animal population the reservoir for the disease. Ebola virus disease (EVD) and SARS are both notable examples of such diseases. There is little work devoted to understanding stochastic disease extinction and reintroduction in the presence of a reservoir. Here we build a stochastic model for EVD and explicitly consider the presence of an animal reservoir. Using a master equation approach and a WKB ansatz, we determine the associated Hamiltonian of the system. Hamilton's equations are then used to numerically compute the 12-dimensional optimal path to extinction, which is then used to estimate mean extinction times. We also numerically investigate the behavior of the model for dynamic population size. Our results provide an improved understanding of outbreak and extinction dynamics in diseases like EVD.

Reformulation of a stochastic action principle for irregular dynamics

International Nuclear Information System (INIS)

Wang, Q.A.; Bangoup, S.; Dzangue, F.; Jeatsa, A.; Tsobnang, F.; Le Mehaute, A.

2009-01-01

A stochastic action principle for random dynamics is revisited. Numerical diffusion experiments are carried out to show that the diffusion path probability depends exponentially on the Lagrangian action A=∫ a b Ldt. This result is then used to derive the Shannon measure for path uncertainty. It is shown that the maximum entropy principle and the least action principle of classical mechanics can be unified into δA-bar=0 where the average is calculated over all possible paths of the stochastic motion between two configuration points a and b. It is argued that this action principle and the maximum entropy principle are a consequence of the mechanical equilibrium condition extended to the case of stochastic dynamics.

On Stochastic Finite-Time Control of Discrete-Time Fuzzy Systems with Packet Dropout

Directory of Open Access Journals (Sweden)

Yingqi Zhang

2012-01-01

Full Text Available This paper is concerned with the stochastic finite-time stability and stochastic finite-time boundedness problems for one family of fuzzy discrete-time systems over networks with packet dropout, parametric uncertainties, and time-varying norm-bounded disturbance. Firstly, we present the dynamic model description studied, in which the discrete-time fuzzy T-S systems with packet loss can be described by one class of fuzzy Markovian jump systems. Then, the concepts of stochastic finite-time stability and stochastic finite-time boundedness and problem formulation are given. Based on Lyapunov function approach, sufficient conditions on stochastic finite-time stability and stochastic finite-time boundedness are established for the resulting closed-loop fuzzy discrete-time system with Markovian jumps, and state-feedback controllers are designed to ensure stochastic finite-time stability and stochastic finite-time boundedness of the class of fuzzy systems. The stochastic finite-time stability and stochastic finite-time boundedness criteria can be tackled in the form of linear matrix inequalities with a fixed parameter. As an auxiliary result, we also give sufficient conditions on the stochastic stability of the class of fuzzy T-S systems with packet loss. Finally, two illustrative examples are presented to show the validity of the developed methodology.

风力发电高塔系统风致随机动力响应分析%Stochastic dynamic analysis of wind turbine systems under wind loads

Institute of Scientific and Technical Information of China (English)

贺广零; 李杰

2011-01-01

Wind energy development is an effective way to solve problems such as energy shortage and environmental pollution, so it is very useful to accurately compute the stochastic dynamic response of wind turbine systems under wind loads. Firstly< a new method termed as generalized probability density evolution method (GPDEM) is presented in the paper. The GPDEM has been proved to be of high accuracy and efficiency in most kinds of stochastic analysis of dynamical systems. Associated with the physical model of stochastic wind field and the integrated finite element model consisting of the rotor, the nacelle. The tower and the foundation, it is very hopeful to apply the GPDEM in stochastic dynamic analysis of wind turbine systems. Then, a 1. 25 MW wind turbine system subject to wind loads is investigated in detail, and the dynamic responses and the stochastic ones are also compared. The results demonstrate that the randomness of wind velocities places a great influence on the stochastic dynamic analysis of wind turbine systems.%介绍了一种高精度且高效的随机动力系统分析方法-广义概率密度演化方法.基于广义概率密度演化方法,结合随机脉动风场物理模型和“桨叶-机舱-塔体-基础”一体化有限元模型,分别对1.25 MW风力发电高钢塔和钢筋混凝土风力发电高塔进行了风致随机动力响应分析,并将分析结果同确定性动力响应分析结果进行比较.研究表明,随机性对风力发电高塔系统结构风致动力响应分析的影响非常显著.

Stochasticity and transport in Hamiltonian systems

International Nuclear Information System (INIS)

MacKay, R.S.; Meiss, J.D.; Percival, I.C.

1983-08-01

The theory of transport in nonlinear dynamics is developed in terms of leaky barriers which remain when invariant tori are destroyed. We describe the organization of stochastic motion by these barriers and give an explanation of long-time correlations in the stochastic regime

Stochastic analysis of a novel nonautonomous periodic SIRI epidemic system with random disturbances

Science.gov (United States)

Zhang, Weiwei; Meng, Xinzhu

2018-02-01

In this paper, a new stochastic nonautonomous SIRI epidemic model is formulated. Given that the incidence rates of diseases may change with the environment, we propose a novel type of transmission function. The main aim of this paper is to obtain the thresholds of the stochastic SIRI epidemic model. To this end, we investigate the dynamics of the stochastic system and establish the conditions for extinction and persistence in mean of the disease by constructing some suitable Lyapunov functions and using stochastic analysis technique. Furthermore, we show that the stochastic system has at least one nontrivial positive periodic solution. Finally, numerical simulations are introduced to illustrate our results.

Forecasting financial asset processes: stochastic dynamics via learning neural networks.

Science.gov (United States)

Giebel, S; Rainer, M

2010-01-01

Models for financial asset dynamics usually take into account their inherent unpredictable nature by including a suitable stochastic component into their process. Unknown (forward) values of financial assets (at a given time in the future) are usually estimated as expectations of the stochastic asset under a suitable risk-neutral measure. This estimation requires the stochastic model to be calibrated to some history of sufficient length in the past. Apart from inherent limitations, due to the stochastic nature of the process, the predictive power is also limited by the simplifying assumptions of the common calibration methods, such as maximum likelihood estimation and regression methods, performed often without weights on the historic time series, or with static weights only. Here we propose a novel method of "intelligent" calibration, using learning neural networks in order to dynamically adapt the parameters of the stochastic model. Hence we have a stochastic process with time dependent parameters, the dynamics of the parameters being themselves learned continuously by a neural network. The back propagation in training the previous weights is limited to a certain memory length (in the examples we consider 10 previous business days), which is similar to the maximal time lag of autoregressive processes. We demonstrate the learning efficiency of the new algorithm by tracking the next-day forecasts for the EURTRY and EUR-HUF exchange rates each.

Stochastic dynamics modeling solute transport in porous media modeling solute transport in porous media

CERN Document Server

Kulasiri, Don

2002-01-01

Most of the natural and biological phenomena such as solute transport in porous media exhibit variability which can not be modeled by using deterministic approaches. There is evidence in natural phenomena to suggest that some of the observations can not be explained by using the models which give deterministic solutions. Stochastic processes have a rich repository of objects which can be used to express the randomness inherent in the system and the evolution of the system over time. The attractiveness of the stochastic differential equations (SDE) and stochastic partial differential equations (SPDE) come from the fact that we can integrate the variability of the system along with the scientific knowledge pertaining to the system. One of the aims of this book is to explaim some useufl concepts in stochastic dynamics so that the scientists and engineers with a background in undergraduate differential calculus could appreciate the applicability and appropriateness of these developments in mathematics. The ideas ...

Functional Abstraction of Stochastic Hybrid Systems

NARCIS (Netherlands)

Bujorianu, L.M.; Blom, Henk A.P.; Hermanns, H.

2006-01-01

The verification problem for stochastic hybrid systems is quite difficult. One method to verify these systems is stochastic reachability analysis. Concepts of abstractions for stochastic hybrid systems are needed to ease the stochastic reachability analysis. In this paper, we set up different ways

Nonperturbative stochastic dynamics driven by strongly correlated colored noise

Science.gov (United States)

Jing, Jun; Li, Rui; You, J. Q.; Yu, Ting

2015-02-01

We propose a quantum model consisting of two remote qubits interacting with two correlated colored noises and establish an exact stochastic Schrödinger equation for this open quantum system. It is shown that the quantum dynamics of the qubit system is profoundly modulated by the mutual correlation between baths and the bath memory capability through dissipation and fluctuation. We report a physical effect on generating inner correlation and entanglement of two distant qubits arising from the strong bath-bath correlation.

Dynamics of bound vector solitons induced by stochastic perturbations: Soliton breakup and soliton switching

International Nuclear Information System (INIS)

Sun, Zhi-Yuan; Gao, Yi-Tian; Yu, Xin; Liu, Ying

2013-01-01

We respectively investigate breakup and switching of the Manakov-typed bound vector solitons (BVSs) induced by two types of stochastic perturbations: the homogenous and nonhomogenous. Symmetry-recovering is discovered for the asymmetrical homogenous case, while soliton switching is found to relate with the perturbation amplitude and soliton coherence. Simulations show that soliton switching in the circularly-polarized light system is much weaker than that in the Manakov and linearly-polarized systems. In addition, the homogenous perturbations can enhance the soliton switching in both of the Manakov and non-integrable (linearly- and circularly-polarized) systems. Our results might be helpful in interpreting dynamics of the BVSs with stochastic noises in nonlinear optics or with stochastic quantum fluctuations in Bose–Einstein condensates.

Stochastic integer programming by dynamic programming

NARCIS (Netherlands)

Lageweg, B.J.; Lenstra, J.K.; Rinnooy Kan, A.H.G.; Stougie, L.; Ermoliev, Yu.; Wets, R.J.B.

1988-01-01

Stochastic integer programming is a suitable tool for modeling hierarchical decision situations with combinatorial features. In continuation of our work on the design and analysis of heuristics for such problems, we now try to find optimal solutions. Dynamic programming techniques can be used to

Dynamic stochastic optimization

CERN Document Server

Ermoliev, Yuri; Pflug, Georg

2004-01-01

Uncertainties and changes are pervasive characteristics of modern systems involving interactions between humans, economics, nature and technology. These systems are often too complex to allow for precise evaluations and, as a result, the lack of proper management (control) may create significant risks. In order to develop robust strategies we need approaches which explic­ itly deal with uncertainties, risks and changing conditions. One rather general approach is to characterize (explicitly or implicitly) uncertainties by objec­ tive or subjective probabilities (measures of confidence or belief). This leads us to stochastic optimization problems which can rarely be solved by using the standard deterministic optimization and optimal control methods. In the stochastic optimization the accent is on problems with a large number of deci­ sion and random variables, and consequently the focus ofattention is directed to efficient solution procedures rather than to (analytical) closed-form solu­ tions. Objective an...

Studies to the stochastic theory of coupled reactorkinetic-thermohydraulic systems Pt. 2

International Nuclear Information System (INIS)

Mesko, L.

1983-06-01

The description is given of the noise phenomena taking place in a multivariable coupled system by a comprehensive model based on the theory of stochastic fluctuations. A comparison is made with models using transfer function formalism for systems characterized by deterministic open and closed loop signal transmission properties. The advantages of the stochastic model are illustrated by simple reactor dynamical examples having diagnostical importance. (author)

Molecular dynamics with deterministic and stochastic numerical methods

CERN Document Server

Leimkuhler, Ben

2015-01-01

This book describes the mathematical underpinnings of algorithms used for molecular dynamics simulation, including both deterministic and stochastic numerical methods. Molecular dynamics is one of the most versatile and powerful methods of modern computational science and engineering and is used widely in chemistry, physics, materials science and biology. Understanding the foundations of numerical methods means knowing how to select the best one for a given problem (from the wide range of techniques on offer) and how to create new, efficient methods to address particular challenges as they arise in complex applications.  Aimed at a broad audience, this book presents the basic theory of Hamiltonian mechanics and stochastic differential equations, as well as topics including symplectic numerical methods, the handling of constraints and rigid bodies, the efficient treatment of Langevin dynamics, thermostats to control the molecular ensemble, multiple time-stepping, and the dissipative particle dynamics method...

Renormalization group and instantons in stochastic nonlinear dynamics, from self-organized criticality to thermonuclear reactors

International Nuclear Information System (INIS)

Volchenkov, D.

2009-01-01

Stochastic counterparts of nonlinear dynamics are studied by means of nonperturbative functional methods developed in the framework of quantum field theory (QFT). In particular, we discuss fully developed turbulence, including leading corrections on possible compressibility of fluids, transport through porous media, theory of waterspouts and tsunami waves, stochastic magnetohydrodynamics, turbulent transport in crossed fields, self-organized criticality, and dynamics of accelerated wrinkled flame fronts advancing in a wide canal. This report would be of interest to the broad auditorium of physicists and applied mathematicians, with a background in nonperturbative QFT methods or nonlinear dynamical systems, having an interest in both methodological developments and interdisciplinary applications. (author)

Renormalization group and instantons in stochastic nonlinear dynamics, from self-organized criticality to thermonuclear reactors

Energy Technology Data Exchange (ETDEWEB)

Volchenkov, D. [Bielefeld Univ., Center of Excellence Cognitive Interaction Technology (CITEC) (Germany)

2009-03-15

Stochastic counterparts of nonlinear dynamics are studied by means of nonperturbative functional methods developed in the framework of quantum field theory (QFT). In particular, we discuss fully developed turbulence, including leading corrections on possible compressibility of fluids, transport through porous media, theory of waterspouts and tsunami waves, stochastic magnetohydrodynamics, turbulent transport in crossed fields, self-organized criticality, and dynamics of accelerated wrinkled flame fronts advancing in a wide canal. This report would be of interest to the broad auditorium of physicists and applied mathematicians, with a background in nonperturbative QFT methods or nonlinear dynamical systems, having an interest in both methodological developments and interdisciplinary applications. (author)

Economic analysis of energy system considering the uncertainties of crude oil, natural gas and nuclear utilization employing stochastic dynamic programming

International Nuclear Information System (INIS)

Hasegawa, Keita; Komiyama, Ryoichi; Fujii, Yasumasa

2016-01-01

The paper presents an economic rationality analysis of power generation mix by stochastic dynamic programming considering fuel price uncertainties and supply disruption risks such as import disruption and nuclear power plant shutdown risk. The situation revolving around Japan's energy security adopted the past statistics, it cannot be applied to a quantitative analysis of future uncertainties. Further objective and quantitative evaluation methods are required in order to analyze Japan's energy system and make it more resilient in sight of long time scale. In this paper, the authors firstly develop the cost minimization model considering oil and natural gas price respectively by stochastic dynamic programming. Then, the authors show several premises of model and an example of result with related to crude oil stockpile, liquefied natural gas stockpile and nuclear power plant capacity. (author)

Stabilizing simulations of complex stochastic representations for quantum dynamical systems

Energy Technology Data Exchange (ETDEWEB)

Perret, C; Petersen, W P, E-mail: wpp@math.ethz.ch [Seminar for Applied Mathematics, ETH, Zurich (Switzerland)

2011-03-04

Path integral representations of quantum dynamics can often be formulated as stochastic differential equations (SDEs). In a series of papers, Corney and Drummond (2004 Phys. Rev. Lett. 93 260401), Deuar and Drummond (2001 Comput. Phys. Commun. 142 442-5), Drummond and Gardnier (1980 J. Phys. A: Math. Gen. 13 2353-68), Gardiner and Zoller (2004 Quantum Noise: A Handbook of Markovian and Non-Markovian Quantum Stochastic Methods with Applications to Quantum Optics (Springer Series in Synergetics) 3rd edn (Berlin: Springer)) and Gilchrist et al (1997 Phys. Rev. A 55 3014-32) and their collaborators have derived SDEs from coherent states representations for density matrices. Computationally, these SDEs are attractive because they seem simple to simulate. They can be quite unstable, however. In this paper, we consider some of the instabilities and propose a few remedies. Particularly, because the variances of the simulated paths typically grow exponentially, the processes become de-localized in relatively short times. Hence, the issues of boundary conditions and stable integration methods become important. We use the Bose-Einstein Hamiltonian as an example. Our results reveal that it is possible to significantly extend integration times and show the periodic structure of certain functionals.

Stochastic sensitivity analysis of the variability of dynamics and transition to chaos in the business cycles model

Science.gov (United States)

Bashkirtseva, Irina; Ryashko, Lev; Ryazanova, Tatyana

2018-01-01

A problem of mathematical modeling of complex stochastic processes in macroeconomics is discussed. For the description of dynamics of income and capital stock, the well-known Kaldor model of business cycles is used as a basic example. The aim of the paper is to give an overview of the variety of stochastic phenomena which occur in Kaldor model forced by additive and parametric random noise. We study a generation of small- and large-amplitude stochastic oscillations, and their mixed-mode intermittency. To analyze these phenomena, we suggest a constructive approach combining the study of the peculiarities of deterministic phase portrait, and stochastic sensitivity of attractors. We show how parametric noise can stabilize the unstable equilibrium and transform dynamics of Kaldor system from order to chaos.

'PSA-SPN' - A Parameter Sensitivity Analysis Method Using Stochastic Petri Nets: Application to a Production Line System

International Nuclear Information System (INIS)

Labadi, Karim; Saggadi, Samira; Amodeo, Lionel

2009-01-01

The dynamic behavior of a discrete event dynamic system can be significantly affected for some uncertain changes in its decision parameters. So, parameter sensitivity analysis would be a useful way in studying the effects of these changes on the system performance. In the past, the sensitivity analysis approaches are frequently based on simulation models. In recent years, formal methods based on stochastic process including Markov process are proposed in the literature. In this paper, we are interested in the parameter sensitivity analysis of discrete event dynamic systems by using stochastic Petri nets models as a tool for modelling and performance evaluation. A sensitivity analysis approach based on stochastic Petri nets, called PSA-SPN method, will be proposed with an application to a production line system.

A fuzzy-stochastic power system planning model: Reflection of dual objectives and dual uncertainties

International Nuclear Information System (INIS)

Zhang, X.Y.; Huang, G.H.; Zhu, H.; Li, Y.P.

2017-01-01

In this study, a fuzzy stochastic dynamic fractional programming (FSDFP) method is proposed for supporting sustainable management of electric power system (EPS) under dual uncertainties. As an improvement upon the mixed-integer linear fractional programming, FSDFP can not only tackle multi-objective issues effectively without setting weights, but also can deal with uncertain parameters which have both stochastic and fuzzy characteristics. Thus, the developed method can help provide valuable information for supporting capacity-expansion planning and in-depth policy analysis of EPS management problems. For demonstrating these advantages, FSDFP has been applied to a case study of a typical regional EPS planning, where the decision makers have to deal with conflicts between economic development that maximizes the system profit and environmental protection that minimizes the carbon dioxide emissions. The obtained results can be analyzed to generate several decision alternatives, and can then help decision makers make suitable decisions under different input scenarios. Furthermore, comparisons of the solution from FSDFP method with that from fuzzy stochastic dynamic linear programming, linear fractional programming and dynamic stochastic fractional programming methods are undertaken. The contrastive analysis reveals that FSDFP is a more effective approach that can better characterize the complexities and uncertainties of real EPS management problems. - Highlights: • A fuzzy stochastic dynamic fractional programming (FSDFP) method is proposed. • FSDFP can address multiple conflicting objectives without setting weights. • FSDFP can reflect dual uncertainties with both stochastic and fuzzy characteristics. • Some reasonable solutions for a case of power system sustainable planning are generated. • Comparisons of the solutions from FSDFP with other optimization methods are undertaken.

Electricity Market Stochastic Dynamic Model and Its Mean Stability Analysis

Directory of Open Access Journals (Sweden)

Zhanhui Lu

2014-01-01

Full Text Available Based on the deterministic dynamic model of electricity market proposed by Alvarado, a stochastic electricity market model, considering the random nature of demand sides, is presented in this paper on the assumption that generator cost function and consumer utility function are quadratic functions. The stochastic electricity market model is a generalization of the deterministic dynamic model. Using the theory of stochastic differential equations, stochastic process theory, and eigenvalue techniques, the determining conditions of the mean stability for this electricity market model under small Gauss type random excitation are provided and testified theoretically. That is, if the demand elasticity of suppliers is nonnegative and the demand elasticity of consumers is negative, then the stochastic electricity market model is mean stable. It implies that the stability can be judged directly by initial data without any computation. Taking deterministic electricity market data combined with small Gauss type random excitation as numerical samples to interpret random phenomena from a statistical perspective, the results indicate the conclusions above are correct, valid, and practical.

Frequency-difference-dependent stochastic resonance in neural systems

Science.gov (United States)

Guo, Daqing; Perc, Matjaž; Zhang, Yangsong; Xu, Peng; Yao, Dezhong

2017-08-01

Biological neurons receive multiple noisy oscillatory signals, and their dynamical response to the superposition of these signals is of fundamental importance for information processing in the brain. Here we study the response of neural systems to the weak envelope modulation signal, which is superimposed by two periodic signals with different frequencies. We show that stochastic resonance occurs at the beat frequency in neural systems at the single-neuron as well as the population level. The performance of this frequency-difference-dependent stochastic resonance is influenced by both the beat frequency and the two forcing frequencies. Compared to a single neuron, a population of neurons is more efficient in detecting the information carried by the weak envelope modulation signal at the beat frequency. Furthermore, an appropriate fine-tuning of the excitation-inhibition balance can further optimize the response of a neural ensemble to the superimposed signal. Our results thus introduce and provide insights into the generation and modulation mechanism of the frequency-difference-dependent stochastic resonance in neural systems.

Stochastic dynamics and control

CERN Document Server

Sun, Jian-Qiao; Zaslavsky, George

2006-01-01

This book is a result of many years of author's research and teaching on random vibration and control. It was used as lecture notes for a graduate course. It provides a systematic review of theory of probability, stochastic processes, and stochastic calculus. The feedback control is also reviewed in the book. Random vibration analyses of SDOF, MDOF and continuous structural systems are presented in a pedagogical order. The application of the random vibration theory to reliability and fatigue analysis is also discussed. Recent research results on fatigue analysis of non-Gaussian stress proc

Nonlinear Damping Identification in Nonlinear Dynamic System Based on Stochastic Inverse Approach

Directory of Open Access Journals (Sweden)

S. L. Han

2012-01-01

Full Text Available The nonlinear model is crucial to prepare, supervise, and analyze mechanical system. In this paper, a new nonparametric and output-only identification procedure for nonlinear damping is studied. By introducing the concept of the stochastic state space, we formulate a stochastic inverse problem for a nonlinear damping. The solution of the stochastic inverse problem is designed as probabilistic expression via the hierarchical Bayesian formulation by considering various uncertainties such as the information insufficiency in parameter of interests or errors in measurement. The probability space is estimated using Markov chain Monte Carlo (MCMC. The applicability of the proposed method is demonstrated through numerical experiment and particular application to a realistic problem related to ship roll motion.

Modeling and simulation of a controlled steam generator in the context of dynamic reliability using a Stochastic Hybrid Automaton

International Nuclear Information System (INIS)

Babykina, Génia; Brînzei, Nicolae; Aubry, Jean-François; Deleuze, Gilles

2016-01-01

The paper proposes a modeling framework to support Monte Carlo simulations of the behavior of a complex industrial system. The aim is to analyze the system dependability in the presence of random events, described by any type of probability distributions. Continuous dynamic evolutions of physical parameters are taken into account by a system of differential equations. Dynamic reliability is chosen as theoretical framework. Based on finite state automata theory, the formal model is built by parallel composition of elementary sub-models using a bottom-up approach. Considerations of a stochastic nature lead to a model called the Stochastic Hybrid Automaton. The Scilab/Scicos open source environment is used for implementation. The case study is carried out on an example of a steam generator of a nuclear power plant. The behavior of the system is studied by exploring its trajectories. Possible system trajectories are analyzed both empirically, using the results of Monte Carlo simulations, and analytically, using the formal system model. The obtained results are show to be relevant. The Stochastic Hybrid Automaton appears to be a suitable tool to address the dynamic reliability problem and to model real systems of high complexity; the bottom-up design provides precision and coherency of the system model. - Highlights: • A part of a nuclear power plant is modeled in the context of dynamic reliability. • Stochastic Hybrid Automaton is used as an input model for Monte Carlo simulations. • The model is formally built using a bottom-up approach. • The behavior of the system is analyzed empirically and analytically. • A formally built SHA shows to be a suitable tool to approach dynamic reliability.

Parallel Stochastic discrete event simulation of calcium dynamics in neuron.

Science.gov (United States)

Ishlam Patoary, Mohammad Nazrul; Tropper, Carl; McDougal, Robert A; Zhongwei, Lin; Lytton, William W

2017-09-26

The intra-cellular calcium signaling pathways of a neuron depends on both biochemical reactions and diffusions. Some quasi-isolated compartments (e.g. spines) are so small and calcium concentrations are so low that one extra molecule diffusing in by chance can make a nontrivial difference in its concentration (percentage-wise). These rare events can affect dynamics discretely in such way that they cannot be evaluated by a deterministic simulation. Stochastic models of such a system provide a more detailed understanding of these systems than existing deterministic models because they capture their behavior at a molecular level. Our research focuses on the development of a high performance parallel discrete event simulation environment, Neuron Time Warp (NTW), which is intended for use in the parallel simulation of stochastic reaction-diffusion systems such as intra-calcium signaling. NTW is integrated with NEURON, a simulator which is widely used within the neuroscience community. We simulate two models, a calcium buffer and a calcium wave model. The calcium buffer model is employed in order to verify the correctness and performance of NTW by comparing it to a serial deterministic simulation in NEURON. We also derived a discrete event calcium wave model from a deterministic model using the stochastic IP3R structure.

Discriminating chaotic and stochastic dynamics through the permutation spectrum test

Energy Technology Data Exchange (ETDEWEB)

Kulp, C. W., E-mail: Kulp@lycoming.edu [Department of Astronomy and Physics, Lycoming College, Williamsport, Pennsylvania 17701 (United States); Zunino, L., E-mail: lucianoz@ciop.unlp.edu.ar [Centro de Investigaciones Ópticas (CONICET La Plata—CIC), C.C. 3, 1897 Gonnet (Argentina); Departamento de Ciencias Básicas, Facultad de Ingeniería, Universidad Nacional de La Plata (UNLP), 1900 La Plata (Argentina)

2014-09-01

In this paper, we propose a new heuristic symbolic tool for unveiling chaotic and stochastic dynamics: the permutation spectrum test. Several numerical examples allow us to confirm the usefulness of the introduced methodology. Indeed, we show that it is robust in situations in which other techniques fail (intermittent chaos, hyperchaotic dynamics, stochastic linear and nonlinear correlated dynamics, and deterministic non-chaotic noise-driven dynamics). We illustrate the applicability and reliability of this pragmatic method by examining real complex time series from diverse scientific fields. Taking into account that the proposed test has the advantages of being conceptually simple and computationally fast, we think that it can be of practical utility as an alternative test for determinism.

Stochastic dynamics of dengue epidemics.

Science.gov (United States)

de Souza, David R; Tomé, Tânia; Pinho, Suani T R; Barreto, Florisneide R; de Oliveira, Mário J

2013-01-01

We use a stochastic Markovian dynamics approach to describe the spreading of vector-transmitted diseases, such as dengue, and the threshold of the disease. The coexistence space is composed of two structures representing the human and mosquito populations. The human population follows a susceptible-infected-recovered (SIR) type dynamics and the mosquito population follows a susceptible-infected-susceptible (SIS) type dynamics. The human infection is caused by infected mosquitoes and vice versa, so that the SIS and SIR dynamics are interconnected. We develop a truncation scheme to solve the evolution equations from which we get the threshold of the disease and the reproductive ratio. The threshold of the disease is also obtained by performing numerical simulations. We found that for certain values of the infection rates the spreading of the disease is impossible, for any death rate of infected mosquitoes.

Adaptive control of chaotic systems with stochastic time varying unknown parameters

Energy Technology Data Exchange (ETDEWEB)

Salarieh, Hassan [Center of Excellence in Design, Robotics and Automation, Department of Mechanical Engineering, Sharif University of Technology, P.O. Box 11365-9567, Azadi Avenue, Tehran (Iran, Islamic Republic of)], E-mail: salarieh@mech.sharif.edu; Alasty, Aria [Center of Excellence in Design, Robotics and Automation, Department of Mechanical Engineering, Sharif University of Technology, P.O. Box 11365-9567, Azadi Avenue, Tehran (Iran, Islamic Republic of)], E-mail: aalasti@sharif.edu

2008-10-15

In this paper based on the Lyapunov stability theorem, an adaptive control scheme is proposed for stabilizing the unstable periodic orbits (UPO) of chaotic systems. It is assumed that the chaotic system has some linearly dependent unknown parameters which are stochastically time varying. The stochastic parameters are modeled through the Weiner process derivative. To demonstrate the effectiveness of the proposed technique it has been applied to the Lorenz, Chen and Rossler dynamical systems, as some case studies. Simulation results indicate that the proposed adaptive controller has a high performance in stabilizing the UPO of chaotic systems in noisy environment.

The dynamics of stochastic processes

DEFF Research Database (Denmark)

Basse-O'Connor, Andreas

In the present thesis the dynamics of stochastic processes is studied with a special attention to the semimartingale property. This is mainly motivated by the fact that semimartingales provide the class of the processes for which it is possible to define a reasonable stochastic calculus due...... to the Bichteler-Dellacherie Theorem. The semimartingale property of Gaussian processes is characterized in terms of their covariance function, spectral measure and spectral representation. In addition, representation and expansion of filtration results are provided as well. Special attention is given to moving...... average processes, and when the driving process is a Lévy or a chaos process the semimartingale property is characterized in the filtration spanned by the driving process and in the natural filtration when the latter is a Brownian motion. To obtain some of the above results an integrability of seminorm...

A stochastic spatial model of HIV dynamics with an asymmetric battle between the virus and the immune system

International Nuclear Information System (INIS)

Lin Hai; Shuai, J W

2010-01-01

A stochastic spatial model based on the Monte Carlo approach is developed to study the dynamics of human immunodeficiency virus (HIV) infection. We aim to propose a more detailed and realistic simulation frame by incorporating many important features of HIV dynamics, which include infections, replications and mutations of viruses, antigen recognitions, activations and proliferations of lymphocytes, and diffusions, encounters and interactions of virions and lymphocytes. Our model successfully reproduces the three-phase pattern observed in HIV infection, and the simulation results for the time distribution from infection to AIDS onset are also in good agreement with the clinical data. The interactions of viruses and the immune system in all the three phases are investigated. We assess the relative importance of various immune system components in the acute phase. The dynamics of how the two important factors, namely the viral diversity and the asymmetric battle between HIV and the immune system, result in AIDS are investigated in detail with the model.

Variance decomposition in stochastic simulators.

Science.gov (United States)

Le Maître, O P; Knio, O M; Moraes, A

2015-06-28

This work aims at the development of a mathematical and computational approach that enables quantification of the inherent sources of stochasticity and of the corresponding sensitivities in stochastic simulations of chemical reaction networks. The approach is based on reformulating the system dynamics as being generated by independent standardized Poisson processes. This reformulation affords a straightforward identification of individual realizations for the stochastic dynamics of each reaction channel, and consequently a quantitative characterization of the inherent sources of stochasticity in the system. By relying on the Sobol-Hoeffding decomposition, the reformulation enables us to perform an orthogonal decomposition of the solution variance. Thus, by judiciously exploiting the inherent stochasticity of the system, one is able to quantify the variance-based sensitivities associated with individual reaction channels, as well as the importance of channel interactions. Implementation of the algorithms is illustrated in light of simulations of simplified systems, including the birth-death, Schlögl, and Michaelis-Menten models.

Variance decomposition in stochastic simulators

Science.gov (United States)

Le Maître, O. P.; Knio, O. M.; Moraes, A.

2015-06-01

This work aims at the development of a mathematical and computational approach that enables quantification of the inherent sources of stochasticity and of the corresponding sensitivities in stochastic simulations of chemical reaction networks. The approach is based on reformulating the system dynamics as being generated by independent standardized Poisson processes. This reformulation affords a straightforward identification of individual realizations for the stochastic dynamics of each reaction channel, and consequently a quantitative characterization of the inherent sources of stochasticity in the system. By relying on the Sobol-Hoeffding decomposition, the reformulation enables us to perform an orthogonal decomposition of the solution variance. Thus, by judiciously exploiting the inherent stochasticity of the system, one is able to quantify the variance-based sensitivities associated with individual reaction channels, as well as the importance of channel interactions. Implementation of the algorithms is illustrated in light of simulations of simplified systems, including the birth-death, Schlögl, and Michaelis-Menten models.

Variance decomposition in stochastic simulators

Energy Technology Data Exchange (ETDEWEB)

Le Maître, O. P., E-mail: olm@limsi.fr [LIMSI-CNRS, UPR 3251, Orsay (France); Knio, O. M., E-mail: knio@duke.edu [Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708 (United States); Moraes, A., E-mail: alvaro.moraesgutierrez@kaust.edu.sa [King Abdullah University of Science and Technology, Thuwal (Saudi Arabia)

2015-06-28

This work aims at the development of a mathematical and computational approach that enables quantification of the inherent sources of stochasticity and of the corresponding sensitivities in stochastic simulations of chemical reaction networks. The approach is based on reformulating the system dynamics as being generated by independent standardized Poisson processes. This reformulation affords a straightforward identification of individual realizations for the stochastic dynamics of each reaction channel, and consequently a quantitative characterization of the inherent sources of stochasticity in the system. By relying on the Sobol-Hoeffding decomposition, the reformulation enables us to perform an orthogonal decomposition of the solution variance. Thus, by judiciously exploiting the inherent stochasticity of the system, one is able to quantify the variance-based sensitivities associated with individual reaction channels, as well as the importance of channel interactions. Implementation of the algorithms is illustrated in light of simulations of simplified systems, including the birth-death, Schlögl, and Michaelis-Menten models.

Variance decomposition in stochastic simulators

KAUST Repository

Le Maî tre, O. P.; Knio, O. M.; Moraes, Alvaro

2015-01-01

This work aims at the development of a mathematical and computational approach that enables quantification of the inherent sources of stochasticity and of the corresponding sensitivities in stochastic simulations of chemical reaction networks. The approach is based on reformulating the system dynamics as being generated by independent standardized Poisson processes. This reformulation affords a straightforward identification of individual realizations for the stochastic dynamics of each reaction channel, and consequently a quantitative characterization of the inherent sources of stochasticity in the system. By relying on the Sobol-Hoeffding decomposition, the reformulation enables us to perform an orthogonal decomposition of the solution variance. Thus, by judiciously exploiting the inherent stochasticity of the system, one is able to quantify the variance-based sensitivities associated with individual reaction channels, as well as the importance of channel interactions. Implementation of the algorithms is illustrated in light of simulations of simplified systems, including the birth-death, Schlögl, and Michaelis-Menten models.

Stochastic ℋ∞ Finite-Time Control of Discrete-Time Systems with Packet Loss

Directory of Open Access Journals (Sweden)

Yingqi Zhang

2012-01-01

Full Text Available This paper investigates the stochastic finite-time stabilization and ℋ∞ control problem for one family of linear discrete-time systems over networks with packet loss, parametric uncertainties, and time-varying norm-bounded disturbance. Firstly, the dynamic model description studied is given, which, if the packet dropout is assumed to be a discrete-time homogenous Markov process, the class of discrete-time linear systems with packet loss can be regarded as Markovian jump systems. Based on Lyapunov function approach, sufficient conditions are established for the resulting closed-loop discrete-time system with Markovian jumps to be stochastic ℋ∞ finite-time boundedness and then state feedback controllers are designed to guarantee stochastic ℋ∞ finite-time stabilization of the class of stochastic systems. The stochastic ℋ∞ finite-time boundedness criteria can be tackled in the form of linear matrix inequalities with a fixed parameter. As an auxiliary result, we also give sufficient conditions on the robust stochastic stabilization of the class of linear systems with packet loss. Finally, simulation examples are presented to illustrate the validity of the developed scheme.

Stochastic partial differential fluid equations as a diffusive limit of deterministic Lagrangian multi-time dynamics.

Science.gov (United States)

Cotter, C J; Gottwald, G A; Holm, D D

2017-09-01

In Holm (Holm 2015 Proc. R. Soc. A 471 , 20140963. (doi:10.1098/rspa.2014.0963)), stochastic fluid equations were derived by employing a variational principle with an assumed stochastic Lagrangian particle dynamics. Here we show that the same stochastic Lagrangian dynamics naturally arises in a multi-scale decomposition of the deterministic Lagrangian flow map into a slow large-scale mean and a rapidly fluctuating small-scale map. We employ homogenization theory to derive effective slow stochastic particle dynamics for the resolved mean part, thereby obtaining stochastic fluid partial equations in the Eulerian formulation. To justify the application of rigorous homogenization theory, we assume mildly chaotic fast small-scale dynamics, as well as a centring condition. The latter requires that the mean of the fluctuating deviations is small, when pulled back to the mean flow.

Fundamentals of stochastic nature sciences

CERN Document Server

Klyatskin, Valery I

2017-01-01

This book addresses the processes of stochastic structure formation in two-dimensional geophysical fluid dynamics based on statistical analysis of Gaussian random fields, as well as stochastic structure formation in dynamic systems with parametric excitation of positive random fields f(r,t) described by partial differential equations. Further, the book considers two examples of stochastic structure formation in dynamic systems with parametric excitation in the presence of Gaussian pumping. In dynamic systems with parametric excitation in space and time, this type of structure formation either happens – or doesn’t! However, if it occurs in space, then this almost always happens (exponentially quickly) in individual realizations with a unit probability. In the case considered, clustering of the field f(r,t) of any nature is a general feature of dynamic fields, and one may claim that structure formation is the Law of Nature for arbitrary random fields of such type. The study clarifies the conditions under wh...

Stochastic Systems Uncertainty Quantification and Propagation

CERN Document Server

Grigoriu, Mircea

2012-01-01

Uncertainty is an inherent feature of both properties of physical systems and the inputs to these systems that needs to be quantified for cost effective and reliable designs. The states of these systems satisfy equations with random entries, referred to as stochastic equations, so that they are random functions of time and/or space. The solution of stochastic equations poses notable technical difficulties that are frequently circumvented by heuristic assumptions at the expense of accuracy and rigor. The main objective of Stochastic Systems is to promoting the development of accurate and efficient methods for solving stochastic equations and to foster interactions between engineers, scientists, and mathematicians. To achieve these objectives Stochastic Systems presents: ·         A clear and brief review of essential concepts on probability theory, random functions, stochastic calculus, Monte Carlo simulation, and functional analysis   ·          Probabilistic models for random variables an...

Efficient stochastic thermostatting of path integral molecular dynamics.

Science.gov (United States)

Ceriotti, Michele; Parrinello, Michele; Markland, Thomas E; Manolopoulos, David E

2010-09-28

The path integral molecular dynamics (PIMD) method provides a convenient way to compute the quantum mechanical structural and thermodynamic properties of condensed phase systems at the expense of introducing an additional set of high frequency normal modes on top of the physical vibrations of the system. Efficiently sampling such a wide range of frequencies provides a considerable thermostatting challenge. Here we introduce a simple stochastic path integral Langevin equation (PILE) thermostat which exploits an analytic knowledge of the free path integral normal mode frequencies. We also apply a recently developed colored noise thermostat based on a generalized Langevin equation (GLE), which automatically achieves a similar, frequency-optimized sampling. The sampling efficiencies of these thermostats are compared with that of the more conventional Nosé-Hoover chain (NHC) thermostat for a number of physically relevant properties of the liquid water and hydrogen-in-palladium systems. In nearly every case, the new PILE thermostat is found to perform just as well as the NHC thermostat while allowing for a computationally more efficient implementation. The GLE thermostat also proves to be very robust delivering a near-optimum sampling efficiency in all of the cases considered. We suspect that these simple stochastic thermostats will therefore find useful application in many future PIMD simulations.

Markov stochasticity coordinates

International Nuclear Information System (INIS)

Eliazar, Iddo

2017-01-01

Markov dynamics constitute one of the most fundamental models of random motion between the states of a system of interest. Markov dynamics have diverse applications in many fields of science and engineering, and are particularly applicable in the context of random motion in networks. In this paper we present a two-dimensional gauging method of the randomness of Markov dynamics. The method–termed Markov Stochasticity Coordinates–is established, discussed, and exemplified. Also, the method is tweaked to quantify the stochasticity of the first-passage-times of Markov dynamics, and the socioeconomic equality and mobility in human societies.

Markov stochasticity coordinates

Energy Technology Data Exchange (ETDEWEB)

Eliazar, Iddo, E-mail: iddo.eliazar@intel.com

2017-01-15

Markov dynamics constitute one of the most fundamental models of random motion between the states of a system of interest. Markov dynamics have diverse applications in many fields of science and engineering, and are particularly applicable in the context of random motion in networks. In this paper we present a two-dimensional gauging method of the randomness of Markov dynamics. The method–termed Markov Stochasticity Coordinates–is established, discussed, and exemplified. Also, the method is tweaked to quantify the stochasticity of the first-passage-times of Markov dynamics, and the socioeconomic equality and mobility in human societies.

Noise-sustained fluctuations in stochastic dynamics with a delay.

Science.gov (United States)

D'Odorico, Paolo; Laio, Francesco; Ridolfi, Luca

2012-04-01

Delayed responses to external drivers are ubiquitous in environmental, social, and biological processes. Delays may induce oscillations, Hopf bifurcations, and instabilities in deterministic systems even in the absence of nonlinearities. Despite recent advances in the study of delayed stochastic differential equations, the interaction of random drivers with delays remains poorly understood. In particular, it is unclear whether noise-induced behaviors may emerge from these interactions. Here we show that noise may enhance and sustain transient periodic oscillations inherent to deterministic delayed systems. We investigate the conditions conducive to the emergence and disappearance of these dynamics in a linear system in the presence of both additive and multiplicative noise.

Stochastic climate theory

NARCIS (Netherlands)

Gottwald, G.A.; Crommelin, D.T.; Franzke, C.L.E.; Franzke, C.L.E.; O'Kane, T.J.

2017-01-01

In this chapter we review stochastic modelling methods in climate science. First we provide a conceptual framework for stochastic modelling of deterministic dynamical systems based on the Mori-Zwanzig formalism. The Mori-Zwanzig equations contain a Markov term, a memory term and a term suggestive of

Multivariate moment closure techniques for stochastic kinetic models

International Nuclear Information System (INIS)

Lakatos, Eszter; Ale, Angelique; Kirk, Paul D. W.; Stumpf, Michael P. H.

2015-01-01

Stochastic effects dominate many chemical and biochemical processes. Their analysis, however, can be computationally prohibitively expensive and a range of approximation schemes have been proposed to lighten the computational burden. These, notably the increasingly popular linear noise approximation and the more general moment expansion methods, perform well for many dynamical regimes, especially linear systems. At higher levels of nonlinearity, it comes to an interplay between the nonlinearities and the stochastic dynamics, which is much harder to capture correctly by such approximations to the true stochastic processes. Moment-closure approaches promise to address this problem by capturing higher-order terms of the temporally evolving probability distribution. Here, we develop a set of multivariate moment-closures that allows us to describe the stochastic dynamics of nonlinear systems. Multivariate closure captures the way that correlations between different molecular species, induced by the reaction dynamics, interact with stochastic effects. We use multivariate Gaussian, gamma, and lognormal closure and illustrate their use in the context of two models that have proved challenging to the previous attempts at approximating stochastic dynamics: oscillations in p53 and Hes1. In addition, we consider a larger system, Erk-mediated mitogen-activated protein kinases signalling, where conventional stochastic simulation approaches incur unacceptably high computational costs

Multivariate moment closure techniques for stochastic kinetic models

Energy Technology Data Exchange (ETDEWEB)

Lakatos, Eszter, E-mail: e.lakatos13@imperial.ac.uk; Ale, Angelique; Kirk, Paul D. W.; Stumpf, Michael P. H., E-mail: m.stumpf@imperial.ac.uk [Department of Life Sciences, Centre for Integrative Systems Biology and Bioinformatics, Imperial College London, London SW7 2AZ (United Kingdom)

2015-09-07

Stochastic effects dominate many chemical and biochemical processes. Their analysis, however, can be computationally prohibitively expensive and a range of approximation schemes have been proposed to lighten the computational burden. These, notably the increasingly popular linear noise approximation and the more general moment expansion methods, perform well for many dynamical regimes, especially linear systems. At higher levels of nonlinearity, it comes to an interplay between the nonlinearities and the stochastic dynamics, which is much harder to capture correctly by such approximations to the true stochastic processes. Moment-closure approaches promise to address this problem by capturing higher-order terms of the temporally evolving probability distribution. Here, we develop a set of multivariate moment-closures that allows us to describe the stochastic dynamics of nonlinear systems. Multivariate closure captures the way that correlations between different molecular species, induced by the reaction dynamics, interact with stochastic effects. We use multivariate Gaussian, gamma, and lognormal closure and illustrate their use in the context of two models that have proved challenging to the previous attempts at approximating stochastic dynamics: oscillations in p53 and Hes1. In addition, we consider a larger system, Erk-mediated mitogen-activated protein kinases signalling, where conventional stochastic simulation approaches incur unacceptably high computational costs.

Automated planning through abstractions in dynamic and stochastic environments

OpenAIRE

Martínez Muñoz, Moisés

2016-01-01

Mención Internacional en el título de doctor Generating sequences of actions - plans - for an automatic system, like a robot, using Automated Planning is particularly diflicult in stochastic and/or dynamic environments. These plans are composed of actions whose execution, in certain scenarios, might fail, which in tum prevents the execution of the rest of the actions in the plan. Also, in some environments, plans must he generated fast, hoth at the start of the execution and after every ex...

A Class of Stochastic Hybrid Systems with State-Dependent Switching Noise

DEFF Research Database (Denmark)

Leth, John-Josef; Rasmussen, Jakob Gulddahl; Schiøler, Henrik

2012-01-01

In this paper, we develop theoretical results based on a proposed method for modeling switching noise for a class of hybrid systems with piecewise linear partitioned state space, and state-depending switching. We devise a stochastic model of such systems, whose global dynamics is governed...

Stochastic Analysis : A Series of Lectures

CERN Document Server

Dozzi, Marco; Flandoli, Franco; Russo, Francesco

2015-01-01

This book presents in thirteen refereed survey articles an overview of modern activity in stochastic analysis, written by leading international experts. The topics addressed include stochastic fluid dynamics and regularization by noise of deterministic dynamical systems; stochastic partial differential equations driven by Gaussian or Lévy noise, including the relationship between parabolic equations and particle systems, and wave equations in a geometric framework; Malliavin calculus and applications to stochastic numerics; stochastic integration in Banach spaces; porous media-type equations; stochastic deformations of classical mechanics and Feynman integrals and stochastic differential equations with reflection. The articles are based on short courses given at the Centre Interfacultaire Bernoulli of the Ecole Polytechnique Fédérale de Lausanne, Switzerland, from January to June 2012. They offer a valuable resource not only for specialists, but also for other researchers and Ph.D. students in the fields o...

Age distribution dynamics with stochastic jumps in mortality.

Science.gov (United States)

Calabrese, Salvatore; Porporato, Amilcare; Laio, Francesco; D'Odorico, Paolo; Ridolfi, Luca

2017-11-01

While deterministic age distribution models have been extensively studied and applied in various disciplines, little work has been devoted to understanding the role of stochasticity in birth and mortality terms. In this paper, we analyse a stochastic M'Kendrick-von Foerster equation in which jumps in mortality represent intense losses of population due to external events. We present explicit solutions for the probability density functions of the age distribution and the total population and for the temporal dynamics of their moments. We also derive the dynamics of the mean age of the population and its harmonic mean. The framework is then used to calculate the age distribution of salt in the soil root zone, where the accumulation of salt by atmospheric deposition is counteracted by plant uptake and by jump losses due to percolation events.

A Dynamic Momentum Compaction Factor Lattice for Improvements to Stochastic Cooling in Storage Rings

Energy Technology Data Exchange (ETDEWEB)

Olivieri, David Nicholas [Massachusetts U., Amherst

1996-01-01

A dynamic momentum compaction factor, also referred to as a dynamic $\\Delta \\gamma \\tau$, lattice for the FNAL Antiproton Source Debuncher Storage Ring is studied, both theoretically and experimentally, for the purpose of improving stochastic precooling, and hence, improving the global antiproton production and stacking performance. A dynamic $\\Delta \\gamma \\tau$ lattice is proposed due to the competing requirements inherent within the Debuncher storage ring upon $\\gamma \\tau$· Specifically, the Debuncher storage ring performs two disparate functions, $(i)$ accepting and debunching a large number of $\\overline{p}$s/pulse at the outset of the production cycle, which would perform ideally with a large value of $\\gamma\\tau$, and $(ii)$ subsequently employing stochastic cooling throughout the remainder of the $\\overline{p}$ production cycle for improved transfer and stacking efficiency into the Accumulator, for which a small value $\\gamma \\tau$ is ideal in order to reduce the diffusive heating caused by the mixing factor. In the initial design of the Debuncher optical lattice, an intermediate value of $\\gamma \\tau$ was chosen as a compromise between the two functional requirements. The goal of the thesis is to improve stochastic precooling by changing $\\gamma \\tau$ between two desired values during each p production cycle. In particular, the dynamic $\\Delta \\gamma \\tau$ lattice accomplishes a reduction in $\\gamma \\tau$, and hence the mixing factor, through an uniform increase to the dispersion throughout the arc sections of the storage ring. Experimental measurements of cooling rates and system performance parameters, with the implementation of the dynamic $\\Delta \\gamma \\tau$ lattice, are in agreement with theoretical predictions based upon a detailed integration of the stochastic cooling Fokker Planck equations. Based upon the consistency between theory and experiment, predictions of cooling rates are presented for future operational

Limit theorems for Markov chains and stochastic properties of dynamical systems by quasi-compactness

CERN Document Server

Hervé, Loïc

2001-01-01

This book shows how techniques from the perturbation theory of operators, applied to a quasi-compact positive kernel, may be used to obtain limit theorems for Markov chains or to describe stochastic properties of dynamical systems. A general framework for this method is given and then applied to treat several specific cases. An essential element of this work is the description of the peripheral spectra of a quasi-compact Markov kernel and of its Fourier-Laplace perturbations. This is first done in the ergodic but non-mixing case. This work is extended by the second author to the non-ergodic case. The only prerequisites for this book are a knowledge of the basic techniques of probability theory and of notions of elementary functional analysis.

Roles of dispersal, stochasticity, and nonlinear dynamics in the spatial structuring of seasonal natural enemy-victim populations

Science.gov (United States)

Patrick C. Tobin; Ottar N. Bjornstad

2005-01-01

Natural enemy-victim systems may exhibit a range of dynamic space-time patterns. We used a theoretical framework to study spatiotemporal structuring in a transient natural enemy-victim system subject to differential rates of dispersal, stochastic forcing, and nonlinear dynamics. Highly mobile natural enemies that attacked less mobile victims were locally spatially...

Nonlinear stochastic exclusion financial dynamics modeling and time-dependent intrinsic detrended cross-correlation

Science.gov (United States)

Zhang, Wei; Wang, Jun

2017-09-01

In attempt to reproduce price dynamics of financial markets, a stochastic agent-based financial price model is proposed and investigated by stochastic exclusion process. The exclusion process, one of interacting particle systems, is usually thought of as modeling particle motion (with the conserved number of particles) in a continuous time Markov process. In this work, the process is utilized to imitate the trading interactions among the investing agents, in order to explain some stylized facts found in financial time series dynamics. To better understand the correlation behaviors of the proposed model, a new time-dependent intrinsic detrended cross-correlation (TDI-DCC) is introduced and performed, also, the autocorrelation analyses are applied in the empirical research. Furthermore, to verify the rationality of the financial price model, the actual return series are also considered to be comparatively studied with the simulation ones. The comparison results of return behaviors reveal that this financial price dynamics model can reproduce some correlation features of actual stock markets.

Multi types DG expansion dynamic planning in distribution system under stochastic conditions using Covariance Matrix Adaptation Evolutionary Strategy and Monte-Carlo simulation

International Nuclear Information System (INIS)

Sadeghi, Mahmood; Kalantar, Mohsen

2014-01-01

Highlights: • Defining a DG dynamic planning problem. • Applying a new evolutionary algorithm called “CMAES” in planning process. • Considering electricity price and fuel price variation stochastic conditions. • Scenario generation and reduction with MCS and backward reduction programs. • Considering approximately all of the costs of the distribution system. - Abstract: This paper presents a dynamic DG planning problem considering uncertainties related to the intermittent nature of the DG technologies such as wind turbines and solar units in addition to the stochastic economic conditions. The stochastic economic situation includes the uncertainties related to the fuel and electricity price of each year. The Monte Carlo simulation is used to generate the possible scenarios of uncertain situations and the produced scenarios are reduced through backward reduction program. The aim of this paper is to maximize the revenue of the distribution system through the benefit cost analysis alongside the encouraging and punishment functions. In order to close to reality, the different growth rates for the planning period are selected. In this paper the Covariance Matrix Adaptation Evolutionary Strategy is introduced and is used to find the best planning scheme of the DG units. The different DG types are considered in the planning problem. The main assumption of this paper is that the DISCO is the owner of the distribution system and the DG units. The proposed method is tested on a 9 bus test distribution system and the results are compared with the known genetic algorithm and PSO methods to show the applicability of the CMAES method in this problem

Stochastic pump effect and geometric phases in dissipative and stochastic systems

Energy Technology Data Exchange (ETDEWEB)

Sinitsyn, Nikolai [Los Alamos National Laboratory

2008-01-01

The success of Berry phases in quantum mechanics stimulated the study of similar phenomena in other areas of physics, including the theory of living cell locomotion and motion of patterns in nonlinear media. More recently, geometric phases have been applied to systems operating in a strongly stochastic environment, such as molecular motors. We discuss such geometric effects in purely classical dissipative stochastic systems and their role in the theory of the stochastic pump effect (SPE).

Tangent map intermittency as an approximate analysis of intermittency in a high dimensional fully stochastic dynamical system: The Tangled Nature model.

Science.gov (United States)

Diaz-Ruelas, Alvaro; Jeldtoft Jensen, Henrik; Piovani, Duccio; Robledo, Alberto

2016-12-01

It is well known that low-dimensional nonlinear deterministic maps close to a tangent bifurcation exhibit intermittency and this circumstance has been exploited, e.g., by Procaccia and Schuster [Phys. Rev. A 28, 1210 (1983)], to develop a general theory of 1/f spectra. This suggests it is interesting to study the extent to which the behavior of a high-dimensional stochastic system can be described by such tangent maps. The Tangled Nature (TaNa) Model of evolutionary ecology is an ideal candidate for such a study, a significant model as it is capable of reproducing a broad range of the phenomenology of macroevolution and ecosystems. The TaNa model exhibits strong intermittency reminiscent of punctuated equilibrium and, like the fossil record of mass extinction, the intermittency in the model is found to be non-stationary, a feature typical of many complex systems. We derive a mean-field version for the evolution of the likelihood function controlling the reproduction of species and find a local map close to tangency. This mean-field map, by our own local approximation, is able to describe qualitatively only one episode of the intermittent dynamics of the full TaNa model. To complement this result, we construct a complete nonlinear dynamical system model consisting of successive tangent bifurcations that generates time evolution patterns resembling those of the full TaNa model in macroscopic scales. The switch from one tangent bifurcation to the next in the sequences produced in this model is stochastic in nature, based on criteria obtained from the local mean-field approximation, and capable of imitating the changing set of types of species and total population in the TaNa model. The model combines full deterministic dynamics with instantaneous parameter random jumps at stochastically drawn times. In spite of the limitations of our approach, which entails a drastic collapse of degrees of freedom, the description of a high-dimensional model system in terms of a low

Improved Stochastic Subspace System Identification for Structural Health Monitoring

Science.gov (United States)

Chang, Chia-Ming; Loh, Chin-Hsiung

2015-07-01

Structural health monitoring acquires structural information through numerous sensor measurements. Vibrational measurement data render the dynamic characteristics of structures to be extracted, in particular of the modal properties such as natural frequencies, damping, and mode shapes. The stochastic subspace system identification has been recognized as a power tool which can present a structure in the modal coordinates. To obtain qualitative identified data, this tool needs to spend computational expense on a large set of measurements. In study, a stochastic system identification framework is proposed to improve the efficiency and quality of the conventional stochastic subspace system identification. This framework includes 1) measured signal processing, 2) efficient space projection, 3) system order selection, and 4) modal property derivation. The measured signal processing employs the singular spectrum analysis algorithm to lower the noise components as well as to present a data set in a reduced dimension. The subspace is subsequently derived from the data set presented in a delayed coordinate. With the proposed order selection criteria, the number of structural modes is determined, resulting in the modal properties. This system identification framework is applied to a real-world bridge for exploring the feasibility in real-time applications. The results show that this improved system identification method significantly decreases computational time, while qualitative modal parameters are still attained.

Spreading dynamics on complex networks: a general stochastic approach.

Science.gov (United States)

Noël, Pierre-André; Allard, Antoine; Hébert-Dufresne, Laurent; Marceau, Vincent; Dubé, Louis J

2014-12-01

Dynamics on networks is considered from the perspective of Markov stochastic processes. We partially describe the state of the system through network motifs and infer any missing data using the available information. This versatile approach is especially well adapted for modelling spreading processes and/or population dynamics. In particular, the generality of our framework and the fact that its assumptions are explicitly stated suggests that it could be used as a common ground for comparing existing epidemics models too complex for direct comparison, such as agent-based computer simulations. We provide many examples for the special cases of susceptible-infectious-susceptible and susceptible-infectious-removed dynamics (e.g., epidemics propagation) and we observe multiple situations where accurate results may be obtained at low computational cost. Our perspective reveals a subtle balance between the complex requirements of a realistic model and its basic assumptions.

Deconvolution effect of near-fault earthquake ground motions on stochastic dynamic response of tunnel-soil deposit interaction systems

Directory of Open Access Journals (Sweden)

K. Hacıefendioğlu

2012-04-01

Full Text Available The deconvolution effect of the near-fault earthquake ground motions on the stochastic dynamic response of tunnel-soil deposit interaction systems are investigated by using the finite element method. Two different earthquake input mechanisms are used to consider the deconvolution effects in the analyses: the standard rigid-base input and the deconvolved-base-rock input model. The Bolu tunnel in Turkey is chosen as a numerical example. As near-fault ground motions, 1999 Kocaeli earthquake ground motion is selected. The interface finite elements are used between tunnel and soil deposit. The mean of maximum values of quasi-static, dynamic and total responses obtained from the two input models are compared with each other.

Estimation of Parameters in Mean-Reverting Stochastic Systems

Directory of Open Access Journals (Sweden)

Tianhai Tian

2014-01-01

Full Text Available Stochastic differential equation (SDE is a very important mathematical tool to describe complex systems in which noise plays an important role. SDE models have been widely used to study the dynamic properties of various nonlinear systems in biology, engineering, finance, and economics, as well as physical sciences. Since a SDE can generate unlimited numbers of trajectories, it is difficult to estimate model parameters based on experimental observations which may represent only one trajectory of the stochastic model. Although substantial research efforts have been made to develop effective methods, it is still a challenge to infer unknown parameters in SDE models from observations that may have large variations. Using an interest rate model as a test problem, in this work we use the Bayesian inference and Markov Chain Monte Carlo method to estimate unknown parameters in SDE models.

Stochastic linear hybrid systems: Modeling, estimation, and application

Science.gov (United States)

Seah, Chze Eng

Hybrid systems are dynamical systems which have interacting continuous state and discrete state (or mode). Accurate modeling and state estimation of hybrid systems are important in many applications. We propose a hybrid system model, known as the Stochastic Linear Hybrid System (SLHS), to describe hybrid systems with stochastic linear system dynamics in each mode and stochastic continuous-state-dependent mode transitions. We then develop a hybrid estimation algorithm, called the State-Dependent-Transition Hybrid Estimation (SDTHE) algorithm, to estimate the continuous state and discrete state of the SLHS from noisy measurements. It is shown that the SDTHE algorithm is more accurate or more computationally efficient than existing hybrid estimation algorithms. Next, we develop a performance analysis algorithm to evaluate the performance of the SDTHE algorithm in a given operating scenario. We also investigate sufficient conditions for the stability of the SDTHE algorithm. The proposed SLHS model and SDTHE algorithm are illustrated to be useful in several applications. In Air Traffic Control (ATC), to facilitate implementations of new efficient operational concepts, accurate modeling and estimation of aircraft trajectories are needed. In ATC, an aircraft's trajectory can be divided into a number of flight modes. Furthermore, as the aircraft is required to follow a given flight plan or clearance, its flight mode transitions are dependent of its continuous state. However, the flight mode transitions are also stochastic due to navigation uncertainties or unknown pilot intents. Thus, we develop an aircraft dynamics model in ATC based on the SLHS. The SDTHE algorithm is then used in aircraft tracking applications to estimate the positions/velocities of aircraft and their flight modes accurately. Next, we develop an aircraft conformance monitoring algorithm to detect any deviations of aircraft trajectories in ATC that might compromise safety. In this application, the SLHS

Complexity Dynamic Character Analysis of Retailers Based on the Share of Stochastic Demand and Service

Directory of Open Access Journals (Sweden)

Junhai Ma

2017-01-01

Full Text Available Apart from the price fluctuation, the retailers’ service level becomes another key factor that affects the market demand. This paper depicts a modified price and demand game model based on the stochastic demand and the retailer’s service level which influences the market demand decided by customers’ preference, while the market demand is stochastic in this model. We explore how the price adjustment speed affects the stability of the supply chain system with respect to service level and stochastic demand. The dynamic behavior of the system is researched by simulation and the stability domain and the bifurcation phenomenon are shown clearly. The largest Lyapunov exponent and the chaotic attractor are also given to confirm the chaotic characteristic of the system. The simulation results indicate that relatively small price adjustment speed may maintain the system at stable state. With the price adjustment speed gradually increasing, the price system gets unstable and finally becomes chaotic. This chaotic phenomenon will perturb the product market and this phenomenon should be controlled to keep the system stay in the stable region. So the chaos control is done and the chaos can be controlled completely. The conclusion makes significant contribution to the system referring to the price fluctuation based on the service level and stochastic demand.

Benchmarking the stochastic time-dependent variational approach for excitation dynamics in molecular aggregates

Energy Technology Data Exchange (ETDEWEB)

Chorošajev, Vladimir [Department of Theoretical Physics, Faculty of Physics, Vilnius University, Sauletekio 9-III, 10222 Vilnius (Lithuania); Gelzinis, Andrius; Valkunas, Leonas [Department of Theoretical Physics, Faculty of Physics, Vilnius University, Sauletekio 9-III, 10222 Vilnius (Lithuania); Department of Molecular Compound Physics, Center for Physical Sciences and Technology, Sauletekio 3, 10222 Vilnius (Lithuania); Abramavicius, Darius, E-mail: darius.abramavicius@ff.vu.lt [Department of Theoretical Physics, Faculty of Physics, Vilnius University, Sauletekio 9-III, 10222 Vilnius (Lithuania)

2016-12-20

Highlights: • The Davydov ansatze can be used for finite temperature simulations with an extension. • The accuracy is high if the system is strongly coupled to the environmental phonons. • The approach can simulate time-resolved fluorescence spectra. - Abstract: Time dependent variational approach is a convenient method to characterize the excitation dynamics in molecular aggregates for different strengths of system-bath interaction a, which does not require any additional perturbative schemes. Until recently, however, this method was only applicable in zero temperature case. It has become possible to extend this method for finite temperatures with the introduction of stochastic time dependent variational approach. Here we present a comparison between this approach and the exact hierarchical equations of motion approach for describing excitation dynamics in a broad range of temperatures. We calculate electronic population evolution, absorption and auxiliary time resolved fluorescence spectra in different regimes and find that the stochastic approach shows excellent agreement with the exact approach when the system-bath coupling is sufficiently large and temperatures are high. The differences between the two methods are larger, when temperatures are lower or the system-bath coupling is small.

Achieving control and synchronization merely through a stochastically adaptive feedback coupling

Science.gov (United States)

Lin, Wei; Chen, Xin; Zhou, Shijie

2017-07-01

Techniques of deterministically adaptive feedback couplings have been successfully and extensively applied to realize control or/and synchronization in chaotic dynamical systems and even in complex dynamical networks. In this article, a technique of stochastically adaptive feedback coupling is novelly proposed to not only realize control in chaotic dynamical systems but also achieve synchronization in unidirectionally coupled systems. Compared with those deterministically adaptive couplings, the proposed stochastic technique interestingly shows some advantages from a physical viewpoint of time and energy consumptions. More significantly, the usefulness of the proposed stochastic technique is analytically validated by the theory of stochastic processes. It is anticipated that the proposed stochastic technique will be widely used in achieving system control and network synchronization.

Modeling cytoskeletal flow over adhesion sites: competition between stochastic bond dynamics and intracellular relaxation

International Nuclear Information System (INIS)

Sabass, Benedikt; Schwarz, Ulrich S

2010-01-01

In migrating cells, retrograde flow of the actin cytoskeleton is related to traction at adhesion sites located at the base of the lamellipodium. The coupling between the moving cytoskeleton and the stationary adhesions is mediated by the continuous association and dissociation of molecular bonds. We introduce a simple model for the competition between the stochastic dynamics of elastic bonds at the moving interface and relaxation within the moving actin cytoskeleton represented by an internal viscous friction coefficient. Using exact stochastic simulations and an analytical mean field theory, we show that the stochastic bond dynamics lead to biphasic friction laws as observed experimentally. At low internal dissipation, stochastic bond dynamics lead to a regime of irregular stick-and-slip motion. High internal dissipation effectively suppresses cooperative effects among bonds and hence stabilizes the adhesion.

A hybrid stochastic hierarchy equations of motion approach to treat the low temperature dynamics of non-Markovian open quantum systems

Science.gov (United States)

Moix, Jeremy M.; Cao, Jianshu

2013-10-01

The hierarchical equations of motion technique has found widespread success as a tool to generate the numerically exact dynamics of non-Markovian open quantum systems. However, its application to low temperature environments remains a serious challenge due to the need for a deep hierarchy that arises from the Matsubara expansion of the bath correlation function. Here we present a hybrid stochastic hierarchical equation of motion (sHEOM) approach that alleviates this bottleneck and leads to a numerical cost that is nearly independent of temperature. Additionally, the sHEOM method generally converges with fewer hierarchy tiers allowing for the treatment of larger systems. Benchmark calculations are presented on the dynamics of two level systems at both high and low temperatures to demonstrate the efficacy of the approach. Then the hybrid method is used to generate the exact dynamics of systems that are nearly impossible to treat by the standard hierarchy. First, exact energy transfer rates are calculated across a broad range of temperatures revealing the deviations from the Förster rates. This is followed by computations of the entanglement dynamics in a system of two qubits at low temperature spanning the weak to strong system-bath coupling regimes.

Nonlinear stochastic interacting dynamics and complexity of financial gasket fractal-like lattice percolation

Science.gov (United States)

Zhang, Wei; Wang, Jun

2018-05-01

A novel nonlinear stochastic interacting price dynamics is proposed and investigated by the bond percolation on Sierpinski gasket fractal-like lattice, aim to make a new approach to reproduce and study the complexity dynamics of real security markets. Fractal-like lattices correspond to finite graphs with vertices and edges, which are similar to fractals, and Sierpinski gasket is a well-known example of fractals. Fractional ordinal array entropy and fractional ordinal array complexity are introduced to analyze the complexity behaviors of financial signals. To deeper comprehend the fluctuation characteristics of the stochastic price evolution, the complexity analysis of random logarithmic returns and volatility are preformed, including power-law distribution, fractional sample entropy and fractional ordinal array complexity. For further verifying the rationality and validity of the developed stochastic price evolution, the actual security market dataset are also studied with the same statistical methods for comparison. The empirical results show that this stochastic price dynamics can reconstruct complexity behaviors of the actual security markets to some extent.

Dimension reduction for stochastic dynamical systems forced onto a manifold by large drift: a constructive approach with examples from theoretical biology

International Nuclear Information System (INIS)

Parsons, Todd L; Rogers, Tim

2017-01-01

Systems composed of large numbers of interacting agents often admit an effective coarse-grained description in terms of a multidimensional stochastic dynamical system, driven by small-amplitude intrinsic noise. In applications to biological, ecological, chemical and social dynamics it is common for these models to posses quantities that are approximately conserved on short timescales, in which case system trajectories are observed to remain close to some lower-dimensional subspace. Here, we derive explicit and general formulae for a reduced-dimension description of such processes that is exact in the limit of small noise and well-separated slow and fast dynamics. The Michaelis–Menten law of enzyme-catalysed reactions, and the link between the Lotka–Volterra and Wright–Fisher processes are explored as a simple worked examples. Extensions of the method are presented for infinite dimensional systems and processes coupled to non-Gaussian noise sources. (paper)

Dimension reduction for stochastic dynamical systems forced onto a manifold by large drift: a constructive approach with examples from theoretical biology

Science.gov (United States)

Parsons, Todd L.; Rogers, Tim

2017-10-01

Systems composed of large numbers of interacting agents often admit an effective coarse-grained description in terms of a multidimensional stochastic dynamical system, driven by small-amplitude intrinsic noise. In applications to biological, ecological, chemical and social dynamics it is common for these models to posses quantities that are approximately conserved on short timescales, in which case system trajectories are observed to remain close to some lower-dimensional subspace. Here, we derive explicit and general formulae for a reduced-dimension description of such processes that is exact in the limit of small noise and well-separated slow and fast dynamics. The Michaelis-Menten law of enzyme-catalysed reactions, and the link between the Lotka-Volterra and Wright-Fisher processes are explored as a simple worked examples. Extensions of the method are presented for infinite dimensional systems and processes coupled to non-Gaussian noise sources.

Stochastic Modelling, Analysis, and Simulations of the Solar Cycle Dynamic Process

Science.gov (United States)

Turner, Douglas C.; Ladde, Gangaram S.

2018-03-01

Analytical solutions, discretization schemes and simulation results are presented for the time delay deterministic differential equation model of the solar dynamo presented by Wilmot-Smith et al. In addition, this model is extended under stochastic Gaussian white noise parametric fluctuations. The introduction of stochastic fluctuations incorporates variables affecting the dynamo process in the solar interior, estimation error of parameters, and uncertainty of the α-effect mechanism. Simulation results are presented and analyzed to exhibit the effects of stochastic parametric volatility-dependent perturbations. The results generalize and extend the work of Hazra et al. In fact, some of these results exhibit the oscillatory dynamic behavior generated by the stochastic parametric additative perturbations in the absence of time delay. In addition, the simulation results of the modified stochastic models influence the change in behavior of the very recently developed stochastic model of Hazra et al.

Estimation of Dynamic Panel Data Models with Stochastic Volatility Using Particle Filters

Directory of Open Access Journals (Sweden)

Wen Xu

2016-10-01

Full Text Available Time-varying volatility is common in macroeconomic data and has been incorporated into macroeconomic models in recent work. Dynamic panel data models have become increasingly popular in macroeconomics to study common relationships across countries or regions. This paper estimates dynamic panel data models with stochastic volatility by maximizing an approximate likelihood obtained via Rao-Blackwellized particle filters. Monte Carlo studies reveal the good and stable performance of our particle filter-based estimator. When the volatility of volatility is high, or when regressors are absent but stochastic volatility exists, our approach can be better than the maximum likelihood estimator which neglects stochastic volatility and generalized method of moments (GMM estimators.

Stochastic differential equations for quantum dynamics of spin-boson networks

International Nuclear Information System (INIS)

Mandt, Stephan; Sadri, Darius; Houck, Andrew A; Türeci, Hakan E

2015-01-01

A popular approach in quantum optics is to map a master equation to a stochastic differential equation, where quantum effects manifest themselves through noise terms. We generalize this approach based on the positive-P representation to systems involving spin, in particular networks or lattices of interacting spins and bosons. We test our approach on a driven dimer of spins and photons, compare it to the master equation, and predict a novel dynamic phase transition in this system. Our numerical approach has scaling advantages over existing methods, but typically requires regularization in terms of drive and dissipation. (paper)

Control of Networked Traffic Flow Distribution - A Stochastic Distribution System Perspective

Energy Technology Data Exchange (ETDEWEB)

Wang, Hong [Pacific Northwest National Laboratory (PNNL); Aziz, H M Abdul [ORNL; Young, Stan [National Renewable Energy Laboratory (NREL); Patil, Sagar [Pacific Northwest National Laboratory (PNNL)

2017-10-01

Networked traffic flow is a common scenario for urban transportation, where the distribution of vehicle queues either at controlled intersections or highway segments reflect the smoothness of the traffic flow in the network. At signalized intersections, the traffic queues are controlled by traffic signal control settings and effective traffic lights control would realize both smooth traffic flow and minimize fuel consumption. Funded by the Energy Efficient Mobility Systems (EEMS) program of the Vehicle Technologies Office of the US Department of Energy, we performed a preliminary investigation on the modelling and control framework in context of urban network of signalized intersections. In specific, we developed a recursive input-output traffic queueing models. The queue formation can be modeled as a stochastic process where the number of vehicles entering each intersection is a random number. Further, we proposed a preliminary B-Spline stochastic model for a one-way single-lane corridor traffic system based on theory of stochastic distribution control.. It has been shown that the developed stochastic model would provide the optimal probability density function (PDF) of the traffic queueing length as a dynamic function of the traffic signal setting parameters. Based upon such a stochastic distribution model, we have proposed a preliminary closed loop framework on stochastic distribution control for the traffic queueing system to make the traffic queueing length PDF follow a target PDF that potentially realizes the smooth traffic flow distribution in a concerned corridor.

A constrained approach to multiscale stochastic simulation of chemically reacting systems

KAUST Repository

Cotter, Simon L.

2011-01-01

Stochastic simulation of coupled chemical reactions is often computationally intensive, especially if a chemical system contains reactions occurring on different time scales. In this paper, we introduce a multiscale methodology suitable to address this problem, assuming that the evolution of the slow species in the system is well approximated by a Langevin process. It is based on the conditional stochastic simulation algorithm (CSSA) which samples from the conditional distribution of the suitably defined fast variables, given values for the slow variables. In the constrained multiscale algorithm (CMA) a single realization of the CSSA is then used for each value of the slow variable to approximate the effective drift and diffusion terms, in a similar manner to the constrained mean-force computations in other applications such as molecular dynamics. We then show how using the ensuing Fokker-Planck equation approximation, we can in turn approximate average switching times in stochastic chemical systems. © 2011 American Institute of Physics.

Stochastic approximation methods-Powerful tools for simulation and optimization: A survey of some recent work on multi-agent systems and cyber-physical systems

International Nuclear Information System (INIS)

Yin, George; Wang, Le Yi; Zhang, Hongwei

2014-01-01

Stochastic approximation methods have found extensive and diversified applications. Recent emergence of networked systems and cyber-physical systems has generated renewed interest in advancing stochastic approximation into a general framework to support algorithm development for information processing and decisions in such systems. This paper presents a survey on some recent developments in stochastic approximation methods and their applications. Using connected vehicles in platoon formation and coordination as a platform, we highlight some traditional and new methodologies of stochastic approximation algorithms and explain how they can be used to capture essential features in networked systems. Distinct features of networked systems with randomly switching topologies, dynamically evolving parameters, and unknown delays are presented, and control strategies are provided

Modeling and analysis of stochastic systems

CERN Document Server

Kulkarni, Vidyadhar G

2011-01-01

Based on the author's more than 25 years of teaching experience, Modeling and Analysis of Stochastic Systems, Second Edition covers the most important classes of stochastic processes used in the modeling of diverse systems, from supply chains and inventory systems to genetics and biological systems. For each class of stochastic process, the text includes its definition, characterization, applications, transient and limiting behavior, first passage times, and cost/reward models. Along with reorganizing the material, this edition revises and adds new exercises and examples. New to the second edi

Chaotic transitions in deterministic and stochastic dynamical systems applications of Melnikov processes in engineering, physics, and neuroscience

CERN Document Server

Simiu, Emil

2002-01-01

The classical Melnikov method provides information on the behavior of deterministic planar systems that may exhibit transitions, i.e. escapes from and captures into preferred regions of phase space. This book develops a unified treatment of deterministic and stochastic systems that extends the applicability of the Melnikov method to physically realizable stochastic planar systems with additive, state-dependent, white, colored, or dichotomous noise. The extended Melnikov method yields the novel result that motions with transitions are chaotic regardless of whether the excitation is deterministic or stochastic. It explains the role in the occurrence of transitions of the characteristics of the system and its deterministic or stochastic excitation, and is a powerful modeling and identification tool. The book is designed primarily for readers interested in applications. The level of preparation required corresponds to the equivalent of a first-year graduate course in applied mathematics. No previous exposure to d...

A stochastic phase-field model determined from molecular dynamics

KAUST Repository

von Schwerin, Erik

2010-03-17

The dynamics of dendritic growth of a crystal in an undercooled melt is determined by macroscopic diffusion-convection of heat and by capillary forces acting on the nanometer scale of the solid-liquid interface width. Its modelling is useful for instance in processing techniques based on casting. The phase-field method is widely used to study evolution of such microstructural phase transformations on a continuum level; it couples the energy equation to a phenomenological Allen-Cahn/Ginzburg-Landau equation modelling the dynamics of an order parameter determining the solid and liquid phases, including also stochastic fluctuations to obtain the qualitatively correct result of dendritic side branching. This work presents a method to determine stochastic phase-field models from atomistic formulations by coarse-graining molecular dynamics. It has three steps: (1) a precise quantitative atomistic definition of the phase-field variable, based on the local potential energy; (2) derivation of its coarse-grained dynamics model, from microscopic Smoluchowski molecular dynamics (that is Brownian or over damped Langevin dynamics); and (3) numerical computation of the coarse-grained model functions. The coarse-grained model approximates Gibbs ensemble averages of the atomistic phase-field, by choosing coarse-grained drift and diffusion functions that minimize the approximation error of observables in this ensemble average. © EDP Sciences, SMAI, 2010.

A stochastic phase-field model determined from molecular dynamics

KAUST Repository

von Schwerin, Erik; Szepessy, Anders

2010-01-01

The dynamics of dendritic growth of a crystal in an undercooled melt is determined by macroscopic diffusion-convection of heat and by capillary forces acting on the nanometer scale of the solid-liquid interface width. Its modelling is useful for instance in processing techniques based on casting. The phase-field method is widely used to study evolution of such microstructural phase transformations on a continuum level; it couples the energy equation to a phenomenological Allen-Cahn/Ginzburg-Landau equation modelling the dynamics of an order parameter determining the solid and liquid phases, including also stochastic fluctuations to obtain the qualitatively correct result of dendritic side branching. This work presents a method to determine stochastic phase-field models from atomistic formulations by coarse-graining molecular dynamics. It has three steps: (1) a precise quantitative atomistic definition of the phase-field variable, based on the local potential energy; (2) derivation of its coarse-grained dynamics model, from microscopic Smoluchowski molecular dynamics (that is Brownian or over damped Langevin dynamics); and (3) numerical computation of the coarse-grained model functions. The coarse-grained model approximates Gibbs ensemble averages of the atomistic phase-field, by choosing coarse-grained drift and diffusion functions that minimize the approximation error of observables in this ensemble average. © EDP Sciences, SMAI, 2010.

Permanence and asymptotic behaviors of stochastic predator-prey system with Markovian switching and Lévy noise

Science.gov (United States)

Wang, Sheng; Wang, Linshan; Wei, Tengda

2018-04-01

This paper concerns the dynamics of a stochastic predator-prey system with Markovian switching and Lévy noise. First, the existence and uniqueness of global positive solution to the system is proved. Then, by combining stochastic analytical techniques with M-matrix analysis, sufficient conditions of stochastic permanence and extinction are obtained. Furthermore, for the stochastic permanence case, by means of four constants related to the stationary probability distribution of the Markov chain and the parameters of the subsystems, both the superior limit and the inferior limit of the average in time of the sample path of the solution are estimated. Finally, our conclusions are illustrated through an example.

Lévy stable noise-induced transitions: stochastic resonance, resonant activation and dynamic hysteresis

International Nuclear Information System (INIS)

Dybiec, Bartłomiej; Gudowska-Nowak, Ewa

2009-01-01

A standard approach to analysis of noise-induced effects in stochastic dynamics assumes a Gaussian character of the noise term describing interaction of the analyzed system with its complex surroundings. An additional assumption about the existence of timescale separation between the dynamics of the measured observable and the typical timescale of the noise allows external fluctuations to be modeled as temporally uncorrelated and therefore white. However, in many natural phenomena the assumptions concerning the above mentioned properties of 'Gaussianity' and 'whiteness' of the noise can be violated. In this context, in contrast to the spatiotemporal coupling characterizing general forms of non-Markovian or semi-Markovian Lévy walks, so called Lévy flights correspond to the class of Markov processes which can still be interpreted as white, but distributed according to a more general, infinitely divisible, stable and non-Gaussian law. Lévy noise-driven non-equilibrium systems are known to manifest interesting physical properties and have been addressed in various scenarios of physical transport exhibiting a superdiffusive behavior. Here we present a brief overview of our recent investigations aimed at understanding features of stochastic dynamics under the influence of Lévy white noise perturbations. We find that the archetypal phenomena of noise-induced ordering are robust and can be detected also in systems driven by memoryless, non-Gaussian, heavy-tailed fluctuations with infinite variance

The response analysis of fractional-order stochastic system via generalized cell mapping method.

Science.gov (United States)

Wang, Liang; Xue, Lili; Sun, Chunyan; Yue, Xiaole; Xu, Wei

2018-01-01

This paper is concerned with the response of a fractional-order stochastic system. The short memory principle is introduced to ensure that the response of the system is a Markov process. The generalized cell mapping method is applied to display the global dynamics of the noise-free system, such as attractors, basins of attraction, basin boundary, saddle, and invariant manifolds. The stochastic generalized cell mapping method is employed to obtain the evolutionary process of probability density functions of the response. The fractional-order ϕ 6 oscillator and the fractional-order smooth and discontinuous oscillator are taken as examples to give the implementations of our strategies. Studies have shown that the evolutionary direction of the probability density function of the fractional-order stochastic system is consistent with the unstable manifold. The effectiveness of the method is confirmed using Monte Carlo results.

Stochastic bifurcation and fractal and chaos control of a giant magnetostrictive film-shape memory alloy composite cantilever plate subjected to in-plane harmonic and stochastic excitation

International Nuclear Information System (INIS)

Zhu, Zhiwen; Zhang, Qingxin; Xu, Jia

2014-01-01

Stochastic bifurcation and fractal and chaos control of a giant magnetostrictive film–shape memory alloy (GMF–SMA) composite cantilever plate subjected to in-plane harmonic and stochastic excitation were studied. Van der Pol items were improved to interpret the hysteretic phenomena of both GMF and SMA, and the nonlinear dynamic model of a GMF–SMA composite cantilever plate subjected to in-plane harmonic and stochastic excitation was developed. The probability density function of the dynamic response of the system was obtained, and the conditions of stochastic Hopf bifurcation were analyzed. The conditions of noise-induced chaotic response were obtained in the stochastic Melnikov integral method, and the fractal boundary of the safe basin of the system was provided. Finally, the chaos control strategy was proposed in the stochastic dynamic programming method. Numerical simulation shows that stochastic Hopf bifurcation and chaos appear in the parameter variation process. The boundary of the safe basin of the system has fractal characteristics, and its area decreases when the noise intensifies. The system reliability was improved through stochastic optimal control, and the safe basin area of the system increased

Stochastic population dynamics of a montane ground-dwelling squirrel.

Directory of Open Access Journals (Sweden)

Jeffrey A Hostetler

Full Text Available Understanding the causes and consequences of population fluctuations is a central goal of ecology. We used demographic data from a long-term (1990-2008 study and matrix population models to investigate factors and processes influencing the dynamics and persistence of a golden-mantled ground squirrel (Callospermophilus lateralis population, inhabiting a dynamic subalpine habitat in Colorado, USA. The overall deterministic population growth rate λ was 0.94±SE 0.05 but it varied widely over time, ranging from 0.45±0.09 in 2006 to 1.50±0.12 in 2003, and was below replacement (λ<1 for 9 out of 18 years. The stochastic population growth rate λ(s was 0.92, suggesting a declining population; however, the 95% CI on λ(s included 1.0 (0.52-1.60. Stochastic elasticity analysis showed that survival of adult females, followed by survival of juvenile females and litter size, were potentially the most influential vital rates; analysis of life table response experiments revealed that the same three life history variables made the largest contributions to year-to year changes in λ. Population viability analysis revealed that, when the influences of density dependence and immigration were not considered, the population had a high (close to 1.0 in 50 years probability of extinction. However, probability of extinction declined to as low as zero when density dependence and immigration were considered. Destabilizing effects of stochastic forces were counteracted by regulating effects of density dependence and rescue effects of immigration, which allowed our study population to bounce back from low densities and prevented extinction. These results suggest that dynamics and persistence of our study population are determined synergistically by density-dependence, stochastic forces, and immigration.

A numerical scheme for optimal transition paths of stochastic chemical kinetic systems

International Nuclear Information System (INIS)

Liu Di

2008-01-01

We present a new framework for finding the optimal transition paths of metastable stochastic chemical kinetic systems with large system size. The optimal transition paths are identified to be the most probable paths according to the Large Deviation Theory of stochastic processes. Dynamical equations for the optimal transition paths are derived using the variational principle. A modified Minimum Action Method (MAM) is proposed as a numerical scheme to solve the optimal transition paths. Applications to Gene Regulatory Networks such as the toggle switch model and the Lactose Operon Model in Escherichia coli are presented as numerical examples

A Simulation-Based Dynamic Stochastic Route Choice Model for Evacuation

Directory of Open Access Journals (Sweden)

Xing Zhao

2012-01-01

Full Text Available This paper establishes a dynamic stochastic route choice model for evacuation to simulate the propagation process of traffic flow and estimate the stochastic route choice under evacuation situations. The model contains a lane-group-based cell transmission model (CTM which sets different traffic capacities for links with different turning movements to flow out in an evacuation situation, an actual impedance model which is to obtain the impedance of each route in time units at each time interval and a stochastic route choice model according to the probit-based stochastic user equilibrium. In this model, vehicles loading at each origin at each time interval are assumed to choose an evacuation route under determinate road network, signal design, and OD demand. As a case study, the proposed model is validated on the network nearby Nanjing Olympic Center after the opening ceremony of the 10th National Games of the People's Republic of China. The traffic volumes and clearing time at five exit points of the evacuation zone are calculated by the model to compare with survey data. The results show that this model can appropriately simulate the dynamic route choice and evolution process of the traffic flow on the network in an evacuation situation.

Computing the optimal path in stochastic dynamical systems

International Nuclear Information System (INIS)

Bauver, Martha; Forgoston, Eric; Billings, Lora

2016-01-01

In stochastic systems, one is often interested in finding the optimal path that maximizes the probability of escape from a metastable state or of switching between metastable states. Even for simple systems, it may be impossible to find an analytic form of the optimal path, and in high-dimensional systems, this is almost always the case. In this article, we formulate a constructive methodology that is used to compute the optimal path numerically. The method utilizes finite-time Lyapunov exponents, statistical selection criteria, and a Newton-based iterative minimizing scheme. The method is applied to four examples. The first example is a two-dimensional system that describes a single population with internal noise. This model has an analytical solution for the optimal path. The numerical solution found using our computational method agrees well with the analytical result. The second example is a more complicated four-dimensional system where our numerical method must be used to find the optimal path. The third example, although a seemingly simple two-dimensional system, demonstrates the success of our method in finding the optimal path where other numerical methods are known to fail. In the fourth example, the optimal path lies in six-dimensional space and demonstrates the power of our method in computing paths in higher-dimensional spaces.

Structure Learning in Stochastic Non-linear Dynamical Systems

Science.gov (United States)

Morris, R. D.; Smelyanskiy, V. N.; Luchinsky, D. G.

2005-12-01

A great many systems can be modeled in the non-linear dynamical systems framework, as x˙ = f(x) + ξ(t), where f(x) is the potential function for the system, and ξ(t) is the driving noise. Modeling the potential using a set of basis functions, we derive the posterior for the basis coefficients. A more challenging problem is to determine the set of basis functions that are required to model a particular system. We show that using the Bayesian Information Criteria (BIC) to rank models, and the beam search technique, that we can accurately determine the structure of simple non-linear dynamical system models, and the structure of the coupling between non-linear dynamical systems where the individual systems are known. This last case has important ecological applications, for example in predator-prey systems, where the very structure of the coupling between predator-prey pairs can have great ecological significance.

Convergence of Sample Path Optimal Policies for Stochastic Dynamic Programming

National Research Council Canada - National Science Library

Fu, Michael C; Jin, Xing

2005-01-01

.... These results have practical implications for Monte Carlo simulation-based solution approaches to stochastic dynamic programming problems where it is impractical to extract the explicit transition...

Stochastic Online Learning in Dynamic Networks under Unknown Models

Science.gov (United States)

2016-08-02

The key is to develop online learning strategies at each individual node. Specifically, through local information exchange with its neighbors, each...infinitely repeated game with incomplete information and developed a dynamic pricing strategy referred to as Competitive and Cooperative Demand Learning...Stochastic Online Learning in Dynamic Networks under Unknown Models This research aims to develop fundamental theories and practical algorithms for

An Error-Entropy Minimization Algorithm for Tracking Control of Nonlinear Stochastic Systems with Non-Gaussian Variables

Energy Technology Data Exchange (ETDEWEB)

Liu, Yunlong; Wang, Aiping; Guo, Lei; Wang, Hong

2017-07-09

This paper presents an error-entropy minimization tracking control algorithm for a class of dynamic stochastic system. The system is represented by a set of time-varying discrete nonlinear equations with non-Gaussian stochastic input, where the statistical properties of stochastic input are unknown. By using Parzen windowing with Gaussian kernel to estimate the probability densities of errors, recursive algorithms are then proposed to design the controller such that the tracking error can be minimized. The performance of the error-entropy minimization criterion is compared with the mean-square-error minimization in the simulation results.

Neural dynamics as sampling: a model for stochastic computation in recurrent networks of spiking neurons.

Science.gov (United States)

Buesing, Lars; Bill, Johannes; Nessler, Bernhard; Maass, Wolfgang

2011-11-01

The organization of computations in networks of spiking neurons in the brain is still largely unknown, in particular in view of the inherently stochastic features of their firing activity and the experimentally observed trial-to-trial variability of neural systems in the brain. In principle there exists a powerful computational framework for stochastic computations, probabilistic inference by sampling, which can explain a large number of macroscopic experimental data in neuroscience and cognitive science. But it has turned out to be surprisingly difficult to create a link between these abstract models for stochastic computations and more detailed models of the dynamics of networks of spiking neurons. Here we create such a link and show that under some conditions the stochastic firing activity of networks of spiking neurons can be interpreted as probabilistic inference via Markov chain Monte Carlo (MCMC) sampling. Since common methods for MCMC sampling in distributed systems, such as Gibbs sampling, are inconsistent with the dynamics of spiking neurons, we introduce a different approach based on non-reversible Markov chains that is able to reflect inherent temporal processes of spiking neuronal activity through a suitable choice of random variables. We propose a neural network model and show by a rigorous theoretical analysis that its neural activity implements MCMC sampling of a given distribution, both for the case of discrete and continuous time. This provides a step towards closing the gap between abstract functional models of cortical computation and more detailed models of networks of spiking neurons.

Dynamic Programming and Error Estimates for Stochastic Control Problems with Maximum Cost

International Nuclear Information System (INIS)

Bokanowski, Olivier; Picarelli, Athena; Zidani, Hasnaa

2015-01-01

This work is concerned with stochastic optimal control for a running maximum cost. A direct approach based on dynamic programming techniques is studied leading to the characterization of the value function as the unique viscosity solution of a second order Hamilton–Jacobi–Bellman (HJB) equation with an oblique derivative boundary condition. A general numerical scheme is proposed and a convergence result is provided. Error estimates are obtained for the semi-Lagrangian scheme. These results can apply to the case of lookback options in finance. Moreover, optimal control problems with maximum cost arise in the characterization of the reachable sets for a system of controlled stochastic differential equations. Some numerical simulations on examples of reachable analysis are included to illustrate our approach

Dynamic Programming and Error Estimates for Stochastic Control Problems with Maximum Cost

Energy Technology Data Exchange (ETDEWEB)

Bokanowski, Olivier, E-mail: boka@math.jussieu.fr [Laboratoire Jacques-Louis Lions, Université Paris-Diderot (Paris 7) UFR de Mathématiques - Bât. Sophie Germain (France); Picarelli, Athena, E-mail: athena.picarelli@inria.fr [Projet Commands, INRIA Saclay & ENSTA ParisTech (France); Zidani, Hasnaa, E-mail: hasnaa.zidani@ensta.fr [Unité de Mathématiques appliquées (UMA), ENSTA ParisTech (France)

2015-02-15

This work is concerned with stochastic optimal control for a running maximum cost. A direct approach based on dynamic programming techniques is studied leading to the characterization of the value function as the unique viscosity solution of a second order Hamilton–Jacobi–Bellman (HJB) equation with an oblique derivative boundary condition. A general numerical scheme is proposed and a convergence result is provided. Error estimates are obtained for the semi-Lagrangian scheme. These results can apply to the case of lookback options in finance. Moreover, optimal control problems with maximum cost arise in the characterization of the reachable sets for a system of controlled stochastic differential equations. Some numerical simulations on examples of reachable analysis are included to illustrate our approach.

Two-boundary first exit time of Gauss-Markov processes for stochastic modeling of acto-myosin dynamics.

Science.gov (United States)

D'Onofrio, Giuseppe; Pirozzi, Enrica

2017-05-01

We consider a stochastic differential equation in a strip, with coefficients suitably chosen to describe the acto-myosin interaction subject to time-varying forces. By simulating trajectories of the stochastic dynamics via an Euler discretization-based algorithm, we fit experimental data and determine the values of involved parameters. The steps of the myosin are represented by the exit events from the strip. Motivated by these results, we propose a specific stochastic model based on the corresponding time-inhomogeneous Gauss-Markov and diffusion process evolving between two absorbing boundaries. We specify the mean and covariance functions of the stochastic modeling process taking into account time-dependent forces including the effect of an external load. We accurately determine the probability density function (pdf) of the first exit time (FET) from the strip by solving a system of two non singular second-type Volterra integral equations via a numerical quadrature. We provide numerical estimations of the mean of FET as approximations of the dwell-time of the proteins dynamics. The percentage of backward steps is given in agreement to experimental data. Numerical and simulation results are compared and discussed.

Applications of Nonlinear Dynamics Model and Design of Complex Systems

CERN Document Server

In, Visarath; Palacios, Antonio

2009-01-01

This edited book is aimed at interdisciplinary, device-oriented, applications of nonlinear science theory and methods in complex systems. In particular, applications directed to nonlinear phenomena with space and time characteristics. Examples include: complex networks of magnetic sensor systems, coupled nano-mechanical oscillators, nano-detectors, microscale devices, stochastic resonance in multi-dimensional chaotic systems, biosensors, and stochastic signal quantization. "applications of nonlinear dynamics: model and design of complex systems" brings together the work of scientists and engineers that are applying ideas and methods from nonlinear dynamics to design and fabricate complex systems.

Perturbation expansions of stochastic wavefunctions for open quantum systems

Science.gov (United States)

Ke, Yaling; Zhao, Yi

2017-11-01

Based on the stochastic unravelling of the reduced density operator in the Feynman path integral formalism for an open quantum system in touch with harmonic environments, a new non-Markovian stochastic Schrödinger equation (NMSSE) has been established that allows for the systematic perturbation expansion in the system-bath coupling to arbitrary order. This NMSSE can be transformed in a facile manner into the other two NMSSEs, i.e., non-Markovian quantum state diffusion and time-dependent wavepacket diffusion method. Benchmarked by numerically exact results, we have conducted a comparative study of the proposed method in its lowest order approximation, with perturbative quantum master equations in the symmetric spin-boson model and the realistic Fenna-Matthews-Olson complex. It is found that our method outperforms the second-order time-convolutionless quantum master equation in the whole parameter regime and even far better than the fourth-order in the slow bath and high temperature cases. Besides, the method is applicable on an equal footing for any kind of spectral density function and is expected to be a powerful tool to explore the quantum dynamics of large-scale systems, benefiting from the wavefunction framework and the time-local appearance within a single stochastic trajectory.

Optimising stochastic trajectories in exact quantum jump approaches of interacting systems

International Nuclear Information System (INIS)

Lacroix, D.

2004-11-01

The standard methods used to substitute the quantum dynamics of two interacting systems by a quantum jump approach based on the Stochastic Schroedinger Equation (SSE) are described. It turns out that for a given situation, there exists an infinite number of SSE reformulation. This fact is used to propose general strategies to optimise the stochastic paths in order to reduce the statistical fluctuations. In this procedure, called the 'adaptative noise method', a specific SSE is obtained for which the noise depends explicitly on both the initial state and on the properties of the interaction Hamiltonian. It is also shown that this method can be further improved by the introduction of a mean-field dynamics. The different optimisation procedures are illustrated quantitatively in the case of interacting spins. A significant reduction of the statistical fluctuations is obtained. Consequently, a much smaller number of trajectories is needed to accurately reproduce the exact dynamics as compared to the standard SSE method. (author)

Extension of Nelson's stochastic quantization to finite temperature using thermo field dynamics

International Nuclear Information System (INIS)

Kobayashi, K.; Yamanaka, Y.

2011-01-01

We present an extension of Nelson's stochastic quantum mechanics to finite temperature. Utilizing the formulation of Thermo Field Dynamics (TFD), we can show that Ito's stochastic equations for tilde and non-tilde particle positions reproduce the TFD-type Schroedinger equation which is equivalent to the Liouville-von Neumann equation. In our formalism, the drift terms in the Ito's stochastic equation have the temperature dependence and the thermal fluctuation is induced through the correlation of the non-tilde and tilde particles. We show that our formalism satisfies the position-momentum uncertainty relation at finite temperature. -- Highlights: → Utilizing TFD, we extend Nelson's stochastic method to finite temperature. → We introduce stochastic equations for tilde and non-tilde particles. → Our stochastic equations can reproduce the TFD-type Schroedinger equation. → Our formalism satisfies the uncertainly relation at finite temperature.

On the precision of quasi steady state assumptions in stochastic dynamics

Science.gov (United States)

Agarwal, Animesh; Adams, Rhys; Castellani, Gastone C.; Shouval, Harel Z.

2012-07-01

Many biochemical networks have complex multidimensional dynamics and there is a long history of methods that have been used for dimensionality reduction for such reaction networks. Usually a deterministic mass action approach is used; however, in small volumes, there are significant fluctuations from the mean which the mass action approach cannot capture. In such cases stochastic simulation methods should be used. In this paper, we evaluate the applicability of one such dimensionality reduction method, the quasi-steady state approximation (QSSA) [L. Menten and M. Michaelis, "Die kinetik der invertinwirkung," Biochem. Z 49, 333369 (1913)] for dimensionality reduction in case of stochastic dynamics. First, the applicability of QSSA approach is evaluated for a canonical system of enzyme reactions. Application of QSSA to such a reaction system in a deterministic setting leads to Michaelis-Menten reduced kinetics which can be used to derive the equilibrium concentrations of the reaction species. In the case of stochastic simulations, however, the steady state is characterized by fluctuations around the mean equilibrium concentration. Our analysis shows that a QSSA based approach for dimensionality reduction captures well the mean of the distribution as obtained from a full dimensional simulation but fails to accurately capture the distribution around that mean. Moreover, the QSSA approximation is not unique. We have then extended the analysis to a simple bistable biochemical network model proposed to account for the stability of synaptic efficacies; the substrate of learning and memory [J. E. Lisman, "A mechanism of memory storage insensitive to molecular turnover: A bistable autophosphorylating kinase," Proc. Natl. Acad. Sci. U.S.A. 82, 3055-3057 (1985)], 10.1073/pnas.82.9.3055. Our analysis shows that a QSSA based dimensionality reduction method results in errors as big as two orders of magnitude in predicting the residence times in the two stable states.

Multiobjective Output Feedback Control of a Class of Stochastic Hybrid Systems with State-Dependent Noise

Directory of Open Access Journals (Sweden)

S. Aberkane

2007-01-01

Full Text Available This paper deals with dynamic output feedback control of continuous-time active fault tolerant control systems with Markovian parameters (AFTCSMP and state-dependent noise. The main contribution is to formulate conditions for multiperformance design, related to this class of stochastic hybrid systems, that take into account the problematic resulting from the fact that the controller only depends on the fault detection and isolation (FDI process. The specifications and objectives under consideration include stochastic stability, ℋ2 and ℋ∞ (or more generally, stochastic integral quadratic constraints performances. Results are formulated as matrix inequalities. The theoretical results are illustrated using a classical example from literature.

Stochastic Models of Polymer Systems

Science.gov (United States)

2016-01-01

Distribution Unlimited Final Report: Stochastic Models of Polymer Systems The views, opinions and/or findings contained in this report are those of the...ADDRESS. Princeton University PO Box 0036 87 Prospect Avenue - 2nd floor Princeton, NJ 08544 -2020 14-Mar-2014 ABSTRACT Number of Papers published in...peer-reviewed journals: Number of Papers published in non peer-reviewed journals: Final Report: Stochastic Models of Polymer Systems Report Title

Filtering and control of stochastic jump hybrid systems

CERN Document Server

Yao, Xiuming; Zheng, Wei Xing

2016-01-01

This book presents recent research work on stochastic jump hybrid systems. Specifically, the considered stochastic jump hybrid systems include Markovian jump Ito stochastic systems, Markovian jump linear-parameter-varying (LPV) systems, Markovian jump singular systems, Markovian jump two-dimensional (2-D) systems, and Markovian jump repeated scalar nonlinear systems. Some sufficient conditions are first established respectively for the stability and performances of those kinds of stochastic jump hybrid systems in terms of solution of linear matrix inequalities (LMIs). Based on the derived analysis conditions, the filtering and control problems are addressed. The book presents up-to-date research developments and novel methodologies on stochastic jump hybrid systems. The contents can be divided into two parts: the first part is focused on robust filter design problem, while the second part is put the emphasis on robust control problem. These methodologies provide a framework for stability and performance analy...

Effects of demographic structure on key properties of stochastic density-independent population dynamics.

Science.gov (United States)

Vindenes, Yngvild; Sæther, Bernt-Erik; Engen, Steinar

2012-12-01

The development of stochastic demography has largely been based on age structured populations, although other types of demographic structure, especially permanent and dynamic heterogeneity, are likely common in natural populations. The combination of stochasticity and demographic structure is a challenge for analyses of population dynamics and extinction risk, because the population structure will fluctuate around the stable structure and the population size shows transient fluctuations. However, by using a diffusion approximation for the total reproductive value, density-independent dynamics of structured populations can be described with only three population parameters: the expected population growth rate, the environmental variance and the demographic variance. These parameters depend on population structure via the state-specific vital rates and transition rates. Once they are found, the diffusion approximation represents a substantial reduction in model complexity. Here, we review and compare the key population parameters across a wide range of demographic structure, from the case of no structure to the most general case of dynamic heterogeneity, and for both discrete and continuous types. We focus on the demographic variance, but also show how environmental stochasticity can be included. This study brings together results from recent models, each considering a specific type of population structure, and places them in a general framework for structured populations. Comparison across different types of demographic structure reveals that the reproductive value is an essential concept for understanding how population structure affects stochastic dynamics and extinction risk. Copyright © 2011 Elsevier Inc. All rights reserved.

Metastable states and quasicycles in a stochastic Wilson-Cowan model of neuronal population dynamics

KAUST Repository

Bressloff, Paul C.

2010-11-03

We analyze a stochastic model of neuronal population dynamics with intrinsic noise. In the thermodynamic limit N→∞, where N determines the size of each population, the dynamics is described by deterministic Wilson-Cowan equations. On the other hand, for finite N the dynamics is described by a master equation that determines the probability of spiking activity within each population. We first consider a single excitatory population that exhibits bistability in the deterministic limit. The steady-state probability distribution of the stochastic network has maxima at points corresponding to the stable fixed points of the deterministic network; the relative weighting of the two maxima depends on the system size. For large but finite N, we calculate the exponentially small rate of noise-induced transitions between the resulting metastable states using a Wentzel-Kramers- Brillouin (WKB) approximation and matched asymptotic expansions. We then consider a two-population excitatory or inhibitory network that supports limit cycle oscillations. Using a diffusion approximation, we reduce the dynamics to a neural Langevin equation, and show how the intrinsic noise amplifies subthreshold oscillations (quasicycles). © 2010 The American Physical Society.

Hybrid Differential Dynamic Programming with Stochastic Search

Science.gov (United States)

Aziz, Jonathan; Parker, Jeffrey; Englander, Jacob

2016-01-01

Differential dynamic programming (DDP) has been demonstrated as a viable approach to low-thrust trajectory optimization, namely with the recent success of NASAs Dawn mission. The Dawn trajectory was designed with the DDP-based Static Dynamic Optimal Control algorithm used in the Mystic software. Another recently developed method, Hybrid Differential Dynamic Programming (HDDP) is a variant of the standard DDP formulation that leverages both first-order and second-order state transition matrices in addition to nonlinear programming (NLP) techniques. Areas of improvement over standard DDP include constraint handling, convergence properties, continuous dynamics, and multi-phase capability. DDP is a gradient based method and will converge to a solution nearby an initial guess. In this study, monotonic basin hopping (MBH) is employed as a stochastic search method to overcome this limitation, by augmenting the HDDP algorithm for a wider search of the solution space.

Non-cooperative stochastic differential game theory of generalized Markov jump linear systems

CERN Document Server

Zhang, Cheng-ke; Zhou, Hai-ying; Bin, Ning

2017-01-01

This book systematically studies the stochastic non-cooperative differential game theory of generalized linear Markov jump systems and its application in the field of finance and insurance. The book is an in-depth research book of the continuous time and discrete time linear quadratic stochastic differential game, in order to establish a relatively complete framework of dynamic non-cooperative differential game theory. It uses the method of dynamic programming principle and Riccati equation, and derives it into all kinds of existence conditions and calculating method of the equilibrium strategies of dynamic non-cooperative differential game. Based on the game theory method, this book studies the corresponding robust control problem, especially the existence condition and design method of the optimal robust control strategy. The book discusses the theoretical results and its applications in the risk control, option pricing, and the optimal investment problem in the field of finance and insurance, enriching the...

Emergence of dynamic cooperativity in the stochastic kinetics of fluctuating enzymes

International Nuclear Information System (INIS)

Kumar, Ashutosh; Chatterjee, Sambarta; Nandi, Mintu; Dua, Arti

2016-01-01

Dynamic co-operativity in monomeric enzymes is characterized in terms of a non-Michaelis-Menten kinetic behaviour. The latter is believed to be associated with mechanisms that include multiple reaction pathways due to enzymatic conformational fluctuations. Recent advances in single-molecule fluorescence spectroscopy have provided new fundamental insights on the possible mechanisms underlying reactions catalyzed by fluctuating enzymes. Here, we present a bottom-up approach to understand enzyme turnover kinetics at physiologically relevant mesoscopic concentrations informed by mechanisms extracted from single-molecule stochastic trajectories. The stochastic approach, presented here, shows the emergence of dynamic co-operativity in terms of a slowing down of the Michaelis-Menten (MM) kinetics resulting in negative co-operativity. For fewer enzymes, dynamic co-operativity emerges due to the combined effects of enzymatic conformational fluctuations and molecular discreteness. The increase in the number of enzymes, however, suppresses the effect of enzymatic conformational fluctuations such that dynamic co-operativity emerges solely due to the discrete changes in the number of reacting species. These results confirm that the turnover kinetics of fluctuating enzyme based on the parallel-pathway MM mechanism switches over to the single-pathway MM mechanism with the increase in the number of enzymes. For large enzyme numbers, convergence to the exact MM equation occurs in the limit of very high substrate concentration as the stochastic kinetics approaches the deterministic behaviour.

Emergence of dynamic cooperativity in the stochastic kinetics of fluctuating enzymes

Energy Technology Data Exchange (ETDEWEB)

Kumar, Ashutosh; Chatterjee, Sambarta; Nandi, Mintu; Dua, Arti, E-mail: arti@iitm.ac.in [Department of Chemistry, Indian Institute of Technology, Madras, Chennai 600036 (India)

2016-08-28

Dynamic co-operativity in monomeric enzymes is characterized in terms of a non-Michaelis-Menten kinetic behaviour. The latter is believed to be associated with mechanisms that include multiple reaction pathways due to enzymatic conformational fluctuations. Recent advances in single-molecule fluorescence spectroscopy have provided new fundamental insights on the possible mechanisms underlying reactions catalyzed by fluctuating enzymes. Here, we present a bottom-up approach to understand enzyme turnover kinetics at physiologically relevant mesoscopic concentrations informed by mechanisms extracted from single-molecule stochastic trajectories. The stochastic approach, presented here, shows the emergence of dynamic co-operativity in terms of a slowing down of the Michaelis-Menten (MM) kinetics resulting in negative co-operativity. For fewer enzymes, dynamic co-operativity emerges due to the combined effects of enzymatic conformational fluctuations and molecular discreteness. The increase in the number of enzymes, however, suppresses the effect of enzymatic conformational fluctuations such that dynamic co-operativity emerges solely due to the discrete changes in the number of reacting species. These results confirm that the turnover kinetics of fluctuating enzyme based on the parallel-pathway MM mechanism switches over to the single-pathway MM mechanism with the increase in the number of enzymes. For large enzyme numbers, convergence to the exact MM equation occurs in the limit of very high substrate concentration as the stochastic kinetics approaches the deterministic behaviour.

Emergence of dynamic cooperativity in the stochastic kinetics of fluctuating enzymes

Science.gov (United States)

Kumar, Ashutosh; Chatterjee, Sambarta; Nandi, Mintu; Dua, Arti

2016-08-01

Dynamic co-operativity in monomeric enzymes is characterized in terms of a non-Michaelis-Menten kinetic behaviour. The latter is believed to be associated with mechanisms that include multiple reaction pathways due to enzymatic conformational fluctuations. Recent advances in single-molecule fluorescence spectroscopy have provided new fundamental insights on the possible mechanisms underlying reactions catalyzed by fluctuating enzymes. Here, we present a bottom-up approach to understand enzyme turnover kinetics at physiologically relevant mesoscopic concentrations informed by mechanisms extracted from single-molecule stochastic trajectories. The stochastic approach, presented here, shows the emergence of dynamic co-operativity in terms of a slowing down of the Michaelis-Menten (MM) kinetics resulting in negative co-operativity. For fewer enzymes, dynamic co-operativity emerges due to the combined effects of enzymatic conformational fluctuations and molecular discreteness. The increase in the number of enzymes, however, suppresses the effect of enzymatic conformational fluctuations such that dynamic co-operativity emerges solely due to the discrete changes in the number of reacting species. These results confirm that the turnover kinetics of fluctuating enzyme based on the parallel-pathway MM mechanism switches over to the single-pathway MM mechanism with the increase in the number of enzymes. For large enzyme numbers, convergence to the exact MM equation occurs in the limit of very high substrate concentration as the stochastic kinetics approaches the deterministic behaviour.

Stochastic dynamics of penetrable rods in one dimension: Entangled dynamics and transport properties

Energy Technology Data Exchange (ETDEWEB)

Craven, Galen T.; Popov, Alexander V.; Hernandez, Rigoberto, E-mail: hernandez@chemistry.gatech.edu [Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400 (United States)

2015-04-21

The dynamical properties of a system of soft rods governed by stochastic hard collisions (SHCs) have been determined over a varying range of softness using molecular dynamics simulations in one dimension and analytic theory. The SHC model allows for interpenetration of the system’s constituent particles in the simulations, generating overlapping clustering behavior analogous to the spatial structures observed in systems governed by deterministic bounded potentials. Through variation of an assigned softness parameter δ, the limiting ranges of intermolecular softness are bridged, connecting the limiting ensemble behavior from hard to ideal (completely soft). Various dynamical and structural observables are measured from simulation and compared to developed theoretical values. The spatial properties are found to be well predicted by theories developed for the deterministic penetrable-sphere model with a transformation from energetic to probabilistic arguments. While the overlapping spatial structures are complex, the dynamical properties can be adequately approximated through a theory built on impulsive interactions with Enskog corrections. Our theory suggests that as the softness of interaction is varied toward the ideal limit, correlated collision processes are less important to the energy transfer mechanism, and Markovian processes dominate the evolution of the configuration space ensemble. For interaction softness close to hard limit, collision processes are highly correlated and overlapping spatial configurations give rise to entanglement of single-particle trajectories.

Stochastic optimal control in infinite dimension dynamic programming and HJB equations

CERN Document Server

Fabbri, Giorgio; Święch, Andrzej

2017-01-01

Providing an introduction to stochastic optimal control in infinite dimension, this book gives a complete account of the theory of second-order HJB equations in infinite-dimensional Hilbert spaces, focusing on its applicability to associated stochastic optimal control problems. It features a general introduction to optimal stochastic control, including basic results (e.g. the dynamic programming principle) with proofs, and provides examples of applications. A complete and up-to-date exposition of the existing theory of viscosity solutions and regular solutions of second-order HJB equations in Hilbert spaces is given, together with an extensive survey of other methods, with a full bibliography. In particular, Chapter 6, written by M. Fuhrman and G. Tessitore, surveys the theory of regular solutions of HJB equations arising in infinite-dimensional stochastic control, via BSDEs. The book is of interest to both pure and applied researchers working in the control theory of stochastic PDEs, and in PDEs in infinite ...

Improving filtering and prediction of spatially extended turbulent systems with model errors through stochastic parameter estimation

International Nuclear Information System (INIS)

Gershgorin, B.; Harlim, J.; Majda, A.J.

2010-01-01

The filtering and predictive skill for turbulent signals is often limited by the lack of information about the true dynamics of the system and by our inability to resolve the assumed dynamics with sufficiently high resolution using the current computing power. The standard approach is to use a simple yet rich family of constant parameters to account for model errors through parameterization. This approach can have significant skill by fitting the parameters to some statistical feature of the true signal; however in the context of real-time prediction, such a strategy performs poorly when intermittent transitions to instability occur. Alternatively, we need a set of dynamic parameters. One strategy for estimating parameters on the fly is a stochastic parameter estimation through partial observations of the true signal. In this paper, we extend our newly developed stochastic parameter estimation strategy, the Stochastic Parameterization Extended Kalman Filter (SPEKF), to filtering sparsely observed spatially extended turbulent systems which exhibit abrupt stability transition from time to time despite a stable average behavior. For our primary numerical example, we consider a turbulent system of externally forced barotropic Rossby waves with instability introduced through intermittent negative damping. We find high filtering skill of SPEKF applied to this toy model even in the case of very sparse observations (with only 15 out of the 105 grid points observed) and with unspecified external forcing and damping. Additive and multiplicative bias corrections are used to learn the unknown features of the true dynamics from observations. We also present a comprehensive study of predictive skill in the one-mode context including the robustness toward variation of stochastic parameters, imperfect initial conditions and finite ensemble effect. Furthermore, the proposed stochastic parameter estimation scheme applied to the same spatially extended Rossby wave system demonstrates

Regular and stochastic particle motion in plasma dynamics

International Nuclear Information System (INIS)

Kaufman, A.N.

1979-08-01

A Hamiltonian formalism is presented for the study of charged-particle trajectories in the self-consistent field of the particles. The intention is to develop a general approach to plasma dynamics. Transformations of phase-space variables are used to separate out the regular, adiabatic motion from the irregular, stochastic trajectories. Several new techniques are included in this presentation

Stochastic stability of four-wheel-steering system

International Nuclear Information System (INIS)

Huang Dongwei; Wang Hongli; Zhu Zhiwen; Feng Zhang

2007-01-01

A four-wheel-steering system subjected to white noise excitations was reduced to a two-degree-of-freedom quasi-non-integrable-Hamiltonian system. Subsequently we obtained an one-dimensional Ito stochastic differential equation for the averaged Hamiltonian of the system by using the stochastic averaging method for quasi-non-integrable-Hamiltonian systems. Thus, the stochastic stability of four-wheel-steering system was analyzed by analyzing the sample behaviors of the averaged Hamiltonian at the boundary H = 0 and calculating its Lyapunov exponent. An example given at the end demonstrated that the conclusion obtained is of considerable significance

An inexact multistage fuzzy-stochastic programming for regional electric power system management constrained by environmental quality.

Science.gov (United States)

Fu, Zhenghui; Wang, Han; Lu, Wentao; Guo, Huaicheng; Li, Wei

2017-12-01

Electric power system involves different fields and disciplines which addressed the economic system, energy system, and environment system. Inner uncertainty of this compound system would be an inevitable problem. Therefore, an inexact multistage fuzzy-stochastic programming (IMFSP) was developed for regional electric power system management constrained by environmental quality. A model which concluded interval-parameter programming, multistage stochastic programming, and fuzzy probability distribution was built to reflect the uncertain information and dynamic variation in the case study, and the scenarios under different credibility degrees were considered. For all scenarios under consideration, corrective actions were allowed to be taken dynamically in accordance with the pre-regulated policies and the uncertainties in reality. The results suggest that the methodology is applicable to handle the uncertainty of regional electric power management systems and help the decision makers to establish an effective development plan.

Stochastic feeding dynamics arise from the need for information and energy.

Science.gov (United States)

Scholz, Monika; Dinner, Aaron R; Levine, Erel; Biron, David

2017-08-29

Animals regulate their food intake in response to the available level of food. Recent observations of feeding dynamics in small animals showed feeding patterns of bursts and pauses, but their function is unknown. Here, we present a data-driven decision-theoretical model of feeding in Caenorhabditis elegans Our central assumption is that food intake serves a dual purpose: to gather information about the external food level and to ingest food when the conditions are good. The model recapitulates experimentally observed feeding patterns. It naturally implements trade-offs between speed versus accuracy and exploration versus exploitation in responding to a dynamic environment. We find that the model predicts three distinct regimes in responding to a dynamical environment, with a transition region where animals respond stochastically to periodic signals. This stochastic response accounts for previously unexplained experimental data.

Stochastic Wilson–Cowan models of neuronal network dynamics with memory and delay

International Nuclear Information System (INIS)

Goychuk, Igor; Goychuk, Andriy

2015-01-01

We consider a simple Markovian class of the stochastic Wilson–Cowan type models of neuronal network dynamics, which incorporates stochastic delay caused by the existence of a refractory period of neurons. From the point of view of the dynamics of the individual elements, we are dealing with a network of non-Markovian stochastic two-state oscillators with memory, which are coupled globally in a mean-field fashion. This interrelation of a higher-dimensional Markovian and lower-dimensional non-Markovian dynamics is discussed in its relevance to the general problem of the network dynamics of complex elements possessing memory. The simplest model of this class is provided by a three-state Markovian neuron with one refractory state, which causes firing delay with an exponentially decaying memory within the two-state reduced model. This basic model is used to study critical avalanche dynamics (the noise sustained criticality) in a balanced feedforward network consisting of the excitatory and inhibitory neurons. Such avalanches emerge due to the network size dependent noise (mesoscopic noise). Numerical simulations reveal an intermediate power law in the distribution of avalanche sizes with the critical exponent around −1.16. We show that this power law is robust upon a variation of the refractory time over several orders of magnitude. However, the avalanche time distribution is biexponential. It does not reflect any genuine power law dependence. (paper)

Dynamic analysis of stochastic transcription cycles.

Directory of Open Access Journals (Sweden)

Claire V Harper

2011-04-01

Full Text Available In individual mammalian cells the expression of some genes such as prolactin is highly variable over time and has been suggested to occur in stochastic pulses. To investigate the origins of this behavior and to understand its functional relevance, we quantitatively analyzed this variability using new mathematical tools that allowed us to reconstruct dynamic transcription rates of different reporter genes controlled by identical promoters in the same living cell. Quantitative microscopic analysis of two reporter genes, firefly luciferase and destabilized EGFP, was used to analyze the dynamics of prolactin promoter-directed gene expression in living individual clonal and primary pituitary cells over periods of up to 25 h. We quantified the time-dependence and cyclicity of the transcription pulses and estimated the length and variation of active and inactive transcription phases. We showed an average cycle period of approximately 11 h and demonstrated that while the measured time distribution of active phases agreed with commonly accepted models of transcription, the inactive phases were differently distributed and showed strong memory, with a refractory period of transcriptional inactivation close to 3 h. Cycles in transcription occurred at two distinct prolactin-promoter controlled reporter genes in the same individual clonal or primary cells. However, the timing of the cycles was independent and out-of-phase. For the first time, we have analyzed transcription dynamics from two equivalent loci in real-time in single cells. In unstimulated conditions, cells showed independent transcription dynamics at each locus. A key result from these analyses was the evidence for a minimum refractory period in the inactive-phase of transcription. The response to acute signals and the result of manipulation of histone acetylation was consistent with the hypothesis that this refractory period corresponded to a phase of chromatin remodeling which significantly

Analytical vs. Simulation Solution Techniques for Pulse Problems in Non-linear Stochastic Dynamics

DEFF Research Database (Denmark)

Iwankiewicz, R.; Nielsen, Søren R. K.

Advantages and disadvantages of available analytical and simulation techniques for pulse problems in non-linear stochastic dynamics are discussed. First, random pulse problems, both those which do and do not lead to Markov theory, are presented. Next, the analytical and analytically-numerical tec......Advantages and disadvantages of available analytical and simulation techniques for pulse problems in non-linear stochastic dynamics are discussed. First, random pulse problems, both those which do and do not lead to Markov theory, are presented. Next, the analytical and analytically...

Stochastic Modelling of Hydrologic Systems

DEFF Research Database (Denmark)

Jonsdottir, Harpa

2007-01-01

In this PhD project several stochastic modelling methods are studied and applied on various subjects in hydrology. The research was prepared at Informatics and Mathematical Modelling at the Technical University of Denmark. The thesis is divided into two parts. The first part contains...... an introduction and an overview of the papers published. Then an introduction to basic concepts in hydrology along with a description of hydrological data is given. Finally an introduction to stochastic modelling is given. The second part contains the research papers. In the research papers the stochastic methods...... are described, as at the time of publication these methods represent new contribution to hydrology. The second part also contains additional description of software used and a brief introduction to stiff systems. The system in one of the papers is stiff....

Stochastic Modelling Of The Repairable System

Directory of Open Access Journals (Sweden)

Andrzejczak Karol

2015-11-01

Full Text Available All reliability models consisting of random time factors form stochastic processes. In this paper we recall the definitions of the most common point processes which are used for modelling of repairable systems. Particularly this paper presents stochastic processes as examples of reliability systems for the support of the maintenance related decisions. We consider the simplest one-unit system with a negligible repair or replacement time, i.e., the unit is operating and is repaired or replaced at failure, where the time required for repair and replacement is negligible. When the repair or replacement is completed, the unit becomes as good as new and resumes operation. The stochastic modelling of recoverable systems constitutes an excellent method of supporting maintenance related decision-making processes and enables their more rational use.

Renormalization group equations in the stochastic quantization scheme

International Nuclear Information System (INIS)

Pugnetti, S.

1987-01-01

We show that there exists a remarkable link between the stochastic quantization and the theory of critical phenomena and dynamical statistical systems. In the stochastic quantization of a field theory, the stochastic Green functions coverge to the quantum ones when the frictious time goes to infinity. We therefore use the typical techniques of the Renormalization Group equations developed in the framework of critical phenomena to discuss some features of the convergence of the stochastic theory. We are also able, in this way, to compute some dynamical critical exponents and give new numerical valuations for them. (orig.)

The zeroth law of thermodynamics and volume-preserving conservative system in equilibrium with stochastic damping

International Nuclear Information System (INIS)

Qian, Hong

2014-01-01

We propose a mathematical formulation of the zeroth law of thermodynamics and develop a stochastic dynamical theory, with a consistent irreversible thermodynamics, for systems possessing sustained conservative stationary current in phase space while in equilibrium with a heat bath. The theory generalizes underdamped mechanical equilibrium: dx=gdt+{−D∇ϕdt+√(2D)dB(t)}, with ∇⋅g=0 and {⋯} respectively representing phase-volume preserving dynamics and stochastic damping. The zeroth law implies stationary distribution u ss (x)=e −ϕ(x) . We find an orthogonality ∇ϕ⋅g=0 as a hallmark of the system. Stochastic thermodynamics based on time reversal (t,ϕ,g)→(−t,ϕ,−g) is formulated: entropy production e p # (t)=−dF(t)/dt; generalized “heat” h d # (t)=−dU(t)/dt, U(t)=∫ R n ϕ(x)u(x,t)dx being “internal energy”, and “free energy” F(t)=U(t)+∫ R n u(x,t)lnu(x,t)dx never increases. Entropy follows (dS)/(dt) =e p # −h d # . Our formulation is shown to be consistent with an earlier theory of P. Ao. Its contradistinctions to other theories, potential-flux decomposition, stochastic Hamiltonian system with even and odd variables, Klein–Kramers equation, Freidlin–Wentzell's theory, and GENERIC, are discussed.

Stochastic Effects; Application in Nuclear Physics

International Nuclear Information System (INIS)

Mazonka, O.

2000-04-01

Stochastic effects in nuclear physics refer to the study of the dynamics of nuclear systems evolving under stochastic equations of motion. In this dissertation we restrict our attention to classical scattering models. We begin with introduction of the model of nuclear dynamics and deterministic equations of evolution. We apply a Langevin approach - an additional property of the model, which reflect the statistical nature of low energy nuclear behaviour. We than concentrate our attention on the problem of calculating tails of distribution functions, which actually is the problem of calculating probabilities of rare outcomes. Two general strategies are proposed. Result and discussion follow. Finally in the appendix we consider stochastic effects in nonequilibrium systems. A few exactly solvable models are presented. For one model we show explicitly that stochastic behaviour in a microscopic description can lead to ordered collective effects on the macroscopic scale. Two others are solved to confirm the predictions of the fluctuation theorem. (author)

Matrix product algorithm for stochastic dynamics on networks applied to nonequilibrium Glauber dynamics

Science.gov (United States)

Barthel, Thomas; De Bacco, Caterina; Franz, Silvio

2018-01-01

We introduce and apply an efficient method for the precise simulation of stochastic dynamical processes on locally treelike graphs. Networks with cycles are treated in the framework of the cavity method. Such models correspond, for example, to spin-glass systems, Boolean networks, neural networks, or other technological, biological, and social networks. Building upon ideas from quantum many-body theory, our approach is based on a matrix product approximation of the so-called edge messages—conditional probabilities of vertex variable trajectories. Computation costs and accuracy can be tuned by controlling the matrix dimensions of the matrix product edge messages (MPEM) in truncations. In contrast to Monte Carlo simulations, the algorithm has a better error scaling and works for both single instances as well as the thermodynamic limit. We employ it to examine prototypical nonequilibrium Glauber dynamics in the kinetic Ising model. Because of the absence of cancellation effects, observables with small expectation values can be evaluated accurately, allowing for the study of decay processes and temporal correlations.

Metaheuristics for the dynamic stochastic dial-a-ride problem with expected return transports.

Science.gov (United States)

Schilde, M; Doerner, K F; Hartl, R F

2011-12-01

The problem of transporting patients or elderly people has been widely studied in literature and is usually modeled as a dial-a-ride problem (DARP). In this paper we analyze the corresponding problem arising in the daily operation of the Austrian Red Cross. This nongovernmental organization is the largest organization performing patient transportation in Austria. The aim is to design vehicle routes to serve partially dynamic transportation requests using a fixed vehicle fleet. Each request requires transportation from a patient's home location to a hospital (outbound request) or back home from the hospital (inbound request). Some of these requests are known in advance. Some requests are dynamic in the sense that they appear during the day without any prior information. Finally, some inbound requests are stochastic. More precisely, with a certain probability each outbound request causes a corresponding inbound request on the same day. Some stochastic information about these return transports is available from historical data. The purpose of this study is to investigate, whether using this information in designing the routes has a significant positive effect on the solution quality. The problem is modeled as a dynamic stochastic dial-a-ride problem with expected return transports. We propose four different modifications of metaheuristic solution approaches for this problem. In detail, we test dynamic versions of variable neighborhood search (VNS) and stochastic VNS (S-VNS) as well as modified versions of the multiple plan approach (MPA) and the multiple scenario approach (MSA). Tests are performed using 12 sets of test instances based on a real road network. Various demand scenarios are generated based on the available real data. Results show that using the stochastic information on return transports leads to average improvements of around 15%. Moreover, improvements of up to 41% can be achieved for some test instances.

GPELab, a Matlab toolbox to solve Gross-Pitaevskii equations II: Dynamics and stochastic simulations

Science.gov (United States)

Antoine, Xavier; Duboscq, Romain

2015-08-01

GPELab is a free Matlab toolbox for modeling and numerically solving large classes of systems of Gross-Pitaevskii equations that arise in the physics of Bose-Einstein condensates. The aim of this second paper, which follows (Antoine and Duboscq, 2014), is to first present the various pseudospectral schemes available in GPELab for computing the deterministic and stochastic nonlinear dynamics of Gross-Pitaevskii equations (Antoine, et al., 2013). Next, the corresponding GPELab functions are explained in detail. Finally, some numerical examples are provided to show how the code works for the complex dynamics of BEC problems.

Integrating stochastic time-dependent travel speed in solution methods for the dynamic dial-a-ride problem.

Science.gov (United States)

Schilde, M; Doerner, K F; Hartl, R F

2014-10-01

In urban areas, logistic transportation operations often run into problems because travel speeds change, depending on the current traffic situation. If not accounted for, time-dependent and stochastic travel speeds frequently lead to missed time windows and thus poorer service. Especially in the case of passenger transportation, it often leads to excessive passenger ride times as well. Therefore, time-dependent and stochastic influences on travel speeds are relevant for finding feasible and reliable solutions. This study considers the effect of exploiting statistical information available about historical accidents, using stochastic solution approaches for the dynamic dial-a-ride problem (dynamic DARP). The authors propose two pairs of metaheuristic solution approaches, each consisting of a deterministic method (average time-dependent travel speeds for planning) and its corresponding stochastic version (exploiting stochastic information while planning). The results, using test instances with up to 762 requests based on a real-world road network, show that in certain conditions, exploiting stochastic information about travel speeds leads to significant improvements over deterministic approaches.

A stochastic MILP energy planning model incorporating power market dynamics

International Nuclear Information System (INIS)

Koltsaklis, Nikolaos E.; Nazos, Konstantinos

2017-01-01

Highlights: •Stochastic MILP model for the optimal energy planning of a power system. •Power market dynamics (offers/bids) are incorporated in the proposed model. •Monte Carlo method for capturing the uncertainty of some key parameters. •Analytical supply cost composition per power producer and activity. •Clean dark and spark spreads are calculated for each power unit. -- Abstract: This paper presents an optimization-based methodological approach to address the problem of the optimal planning of a power system at an annual level in competitive and uncertain power markets. More specifically, a stochastic mixed integer linear programming model (MILP) has been developed, combining advanced optimization techniques with Monte Carlo method in order to deal with uncertainty issues. The main focus of the proposed framework is the dynamic formulation of the strategy followed by all market participants in volatile market conditions, as well as detailed economic assessment of the power system’s operation. The applicability of the proposed approach has been tested on a real case study of the interconnected Greek power system, quantifying in detail all the relevant technical and economic aspects of the system’s operation. The proposed work identifies in the form of probability distributions the optimal power generation mix, electricity trade at a regional level, carbon footprint, as well as detailed total supply cost composition, according to the assumed market structure. The paper demonstrates that the proposed optimization approach is able to provide important insights into the appropriate energy strategies designed by market participants, as well as on the strategic long-term decisions to be made by investors and/or policy makers at a national and/or regional level, underscoring potential risks and providing appropriate price signals on critical energy projects under real market operating conditions.

Stochastic network optimization with application to communication and queueing systems

CERN Document Server

Neely, Michael

2010-01-01

This text presents a modern theory of analysis, control, and optimization for dynamic networks. Mathematical techniques of Lyapunov drift and Lyapunov optimization are developed and shown to enable constrained optimization of time averages in general stochastic systems. The focus is on communication and queueing systems, including wireless networks with time-varying channels, mobility, and randomly arriving traffic. A simple drift-plus-penalty framework is used to optimize time averages such as throughput, throughput-utility, power, and distortion. Explicit performance-delay tradeoffs are prov

A principle of fractal-stochastic dualism and Gompertzian dynamics of growth and self-organization.

Science.gov (United States)

Waliszewski, Przemyslaw

2005-10-01

The emergence of Gompertzian dynamics at the macroscopic, tissue level during growth and self-organization is determined by the existence of fractal-stochastic dualism at the microscopic level of supramolecular, cellular system. On one hand, Gompertzian dynamics results from the complex coupling of at least two antagonistic, stochastic processes at the molecular cellular level. It is shown that the Gompertz function is a probability function, its derivative is a probability density function, and the Gompertzian distribution of probability is of non-Gaussian type. On the other hand, the Gompertz function is a contraction mapping and defines fractal dynamics in time-space; a prerequisite condition for the coupling of processes. Furthermore, the Gompertz function is a solution of the operator differential equation with the Morse-like anharmonic potential. This relationship indicates that distribution of intrasystemic forces is both non-linear and asymmetric. The anharmonic potential is a measure of the intrasystemic interactions. It attains a point of the minimum (U(0), t(0)) along with a change of both complexity and connectivity during growth and self-organization. It can also be modified by certain factors, such as retinoids.

Numerical schemes for dynamically orthogonal equations of stochastic fluid and ocean flows

International Nuclear Information System (INIS)

Ueckermann, M.P.; Lermusiaux, P.F.J.; Sapsis, T.P.

2013-01-01

The quantification of uncertainties is critical when systems are nonlinear and have uncertain terms in their governing equations or are constrained by limited knowledge of initial and boundary conditions. Such situations are common in multiscale, intermittent and non-homogeneous fluid and ocean flows. The dynamically orthogonal (DO) field equations provide an adaptive methodology to predict the probability density functions of such flows. The present work derives efficient computational schemes for the DO methodology applied to unsteady stochastic Navier–Stokes and Boussinesq equations, and illustrates and studies the numerical aspects of these schemes. Semi-implicit projection methods are developed for the mean and for the DO modes, and time-marching schemes of first to fourth order are used for the stochastic coefficients. Conservative second-order finite-volumes are employed in physical space with new advection schemes based on total variation diminishing methods. Other results include: (i) the definition of pseudo-stochastic pressures to obtain a number of pressure equations that is linear in the subspace size instead of quadratic; (ii) symmetric advection schemes for the stochastic velocities; (iii) the use of generalized inversion to deal with singular subspace covariances or deterministic modes; and (iv) schemes to maintain orthonormal modes at the numerical level. To verify our implementation and study the properties of our schemes and their variations, a set of stochastic flow benchmarks are defined including asymmetric Dirac and symmetric lock-exchange flows, lid-driven cavity flows, and flows past objects in a confined channel. Different Reynolds number and Grashof number regimes are employed to illustrate robustness. Optimal convergence under both time and space refinements is shown as well as the convergence of the probability density functions with the number of stochastic realizations.

Stochastic Reformulations of Linear Systems: Algorithms and Convergence Theory

KAUST Repository

Richtarik, Peter; Taká č, Martin

2017-01-01

We develop a family of reformulations of an arbitrary consistent linear system into a stochastic problem. The reformulations are governed by two user-defined parameters: a positive definite matrix defining a norm, and an arbitrary discrete or continuous distribution over random matrices. Our reformulation has several equivalent interpretations, allowing for researchers from various communities to leverage their domain specific insights. In particular, our reformulation can be equivalently seen as a stochastic optimization problem, stochastic linear system, stochastic fixed point problem and a probabilistic intersection problem. We prove sufficient, and necessary and sufficient conditions for the reformulation to be exact. Further, we propose and analyze three stochastic algorithms for solving the reformulated problem---basic, parallel and accelerated methods---with global linear convergence rates. The rates can be interpreted as condition numbers of a matrix which depends on the system matrix and on the reformulation parameters. This gives rise to a new phenomenon which we call stochastic preconditioning, and which refers to the problem of finding parameters (matrix and distribution) leading to a sufficiently small condition number. Our basic method can be equivalently interpreted as stochastic gradient descent, stochastic Newton method, stochastic proximal point method, stochastic fixed point method, and stochastic projection method, with fixed stepsize (relaxation parameter), applied to the reformulations.

Stochastic Reformulations of Linear Systems: Algorithms and Convergence Theory

KAUST Repository

Richtarik, Peter

2017-06-04

We develop a family of reformulations of an arbitrary consistent linear system into a stochastic problem. The reformulations are governed by two user-defined parameters: a positive definite matrix defining a norm, and an arbitrary discrete or continuous distribution over random matrices. Our reformulation has several equivalent interpretations, allowing for researchers from various communities to leverage their domain specific insights. In particular, our reformulation can be equivalently seen as a stochastic optimization problem, stochastic linear system, stochastic fixed point problem and a probabilistic intersection problem. We prove sufficient, and necessary and sufficient conditions for the reformulation to be exact. Further, we propose and analyze three stochastic algorithms for solving the reformulated problem---basic, parallel and accelerated methods---with global linear convergence rates. The rates can be interpreted as condition numbers of a matrix which depends on the system matrix and on the reformulation parameters. This gives rise to a new phenomenon which we call stochastic preconditioning, and which refers to the problem of finding parameters (matrix and distribution) leading to a sufficiently small condition number. Our basic method can be equivalently interpreted as stochastic gradient descent, stochastic Newton method, stochastic proximal point method, stochastic fixed point method, and stochastic projection method, with fixed stepsize (relaxation parameter), applied to the reformulations.

Stochastic Dynamics of a Time-Delayed Ecosystem Driven by Poisson White Noise Excitation

Directory of Open Access Journals (Sweden)

Wantao Jia

2018-02-01

Full Text Available We investigate the stochastic dynamics of a prey-predator type ecosystem with time delay and the discrete random environmental fluctuations. In this model, the delay effect is represented by a time delay parameter and the effect of the environmental randomness is modeled as Poisson white noise. The stochastic averaging method and the perturbation method are applied to calculate the approximate stationary probability density functions for both predator and prey populations. The influences of system parameters and the Poisson white noises are investigated in detail based on the approximate stationary probability density functions. It is found that, increasing time delay parameter as well as the mean arrival rate and the variance of the amplitude of the Poisson white noise will enhance the fluctuations of the prey and predator population. While the larger value of self-competition parameter will reduce the fluctuation of the system. Furthermore, the results from Monte Carlo simulation are also obtained to show the effectiveness of the results from averaging method.

Stochastic resonance and coherence resonance in groundwater-dependent plant ecosystems.

Science.gov (United States)

Borgogno, Fabio; D'Odorico, Paolo; Laio, Francesco; Ridolfi, Luca

2012-01-21

Several studies have shown that non-linear deterministic dynamical systems forced by external random components can give rise to unexpectedly regular temporal behaviors. Stochastic resonance and coherence resonance, the two best known processes of this type, have been studied in a number of physical and chemical systems. Here, we explore their possible occurrence in the dynamics of groundwater-dependent plant ecosystems. To this end, we develop two eco-hydrological models, which allow us to demonstrate that stochastic and coherence resonance may emerge in the dynamics of phreatophyte vegetation, depending on their deterministic properties and the intensity of external stochastic drivers. Copyright © 2011 Elsevier Ltd. All rights reserved.

Stochastic resonance in a stochastic bistable system with additive noises and square–wave signal

International Nuclear Information System (INIS)

Feng, Guo; Xiang-Dong, Luo; Shao-Fu, Li; Yu-Rong, Zhou

2010-01-01

This paper considers the stochastic resonance in a stochastic bistable system driven by a periodic square-wave signal and a static force as well as by additive white noise and dichotomous noise from the viewpoint of signal-to-noise ratio. It finds that the signal-to-noise ratio appears as stochastic resonance behaviour when it is plotted as a function of the noise strength of the white noise and dichotomous noise, as a function of the system parameters, or as a function of the static force. Moreover, the influence of the strength of the stochastic potential force and the correlation rate of the dichotomous noise on the signal-to-noise ratio is investigated. (general)

Stochastic evolution of the Universe: A possible dynamical process ...

Indian Academy of Sciences (India)

C Sivakumar

2017-12-11

Dec 11, 2017 ... https://doi.org/10.1007/s12043-017-1491-z. Stochastic evolution of the Universe: A possible dynamical process leading to fractal structures. C SIVAKUMAR. Department of Physics, Maharaja's College, Ernakulam 682 011, India. E-mail: thrisivc@yahoo.com. MS received 6 July 2016; revised 26 June 2017; ...

Stochastic linear dynamical programming in order to apply it in energy modelling

Energy Technology Data Exchange (ETDEWEB)

El Hachem, S

1995-11-01

This thesis contributes to the development of new algorithms for the computation of stochastic dynamic problem and its mini-maxi variant for the case of imperfect knowledge on random data. The proposed algorithms are scenarios aggregation type. It also contributes to integrate these algorithms in a decision support approach and to discuss the stochastic modeling of two energy problems: the refining and the portfolio gas contracts. (author). 112 refs., 5 tabs.

Pricing decisions in an experimental dynamic stochastic general equilibrium economy

NARCIS (Netherlands)

Noussair, C.N.; Pfajfar, D.; Zsiros, J.

We construct experimental economies, populated with human subjects, with a structure based on a nonlinear version of the New Keynesian dynamic stochastic general equilibrium (DSGE) model. We analyze the behavior of firms’ pricing decisions in four different experimental economies. We consider how

Homogenization of the stochastic Navier–Stokes equation with a stochastic slip boundary condition

KAUST Repository

Bessaih, Hakima

2015-11-02

The two-dimensional Navier–Stokes equation in a perforated domain with a dynamical slip boundary condition is considered. We assume that the dynamic is driven by a stochastic perturbation on the interior of the domain and another stochastic perturbation on the boundaries of the holes. We consider a scaling (ᵋ for the viscosity and 1 for the density) that will lead to a time-dependent limit problem. However, the noncritical scaling (ᵋ, β > 1) is considered in front of the nonlinear term. The homogenized system in the limit is obtained as a Darcy’s law with memory with two permeabilities and an extra term that is due to the stochastic perturbation on the boundary of the holes. The nonhomogeneity on the boundary contains a stochastic part that yields in the limit an additional term in the Darcy’s law. We use the two-scale convergence method after extending the solution with 0 inside the holes to pass to the limit. By Itô stochastic calculus, we get uniform estimates on the solution in appropriate spaces. Due to the stochastic integral, the pressure that appears in the variational formulation does not have enough regularity in time. This fact made us rely only on the variational formulation for the passage to the limit on the solution. We obtain a variational formulation for the limit that is solution of a Stokes system with two pressures. This two-scale limit gives rise to three cell problems, two of them give the permeabilities while the third one gives an extra term in the Darcy’s law due to the stochastic perturbation on the boundary of the holes.

Parametric sensitivity analysis for stochastic molecular systems using information theoretic metrics

Energy Technology Data Exchange (ETDEWEB)

Tsourtis, Anastasios, E-mail: tsourtis@uoc.gr [Department of Mathematics and Applied Mathematics, University of Crete, Crete (Greece); Pantazis, Yannis, E-mail: pantazis@math.umass.edu; Katsoulakis, Markos A., E-mail: markos@math.umass.edu [Department of Mathematics and Statistics, University of Massachusetts, Amherst, Massachusetts 01003 (United States); Harmandaris, Vagelis, E-mail: harman@uoc.gr [Department of Mathematics and Applied Mathematics, University of Crete, and Institute of Applied and Computational Mathematics (IACM), Foundation for Research and Technology Hellas (FORTH), GR-70013 Heraklion, Crete (Greece)

2015-07-07

In this paper, we present a parametric sensitivity analysis (SA) methodology for continuous time and continuous space Markov processes represented by stochastic differential equations. Particularly, we focus on stochastic molecular dynamics as described by the Langevin equation. The utilized SA method is based on the computation of the information-theoretic (and thermodynamic) quantity of relative entropy rate (RER) and the associated Fisher information matrix (FIM) between path distributions, and it is an extension of the work proposed by Y. Pantazis and M. A. Katsoulakis [J. Chem. Phys. 138, 054115 (2013)]. A major advantage of the pathwise SA method is that both RER and pathwise FIM depend only on averages of the force field; therefore, they are tractable and computable as ergodic averages from a single run of the molecular dynamics simulation both in equilibrium and in non-equilibrium steady state regimes. We validate the performance of the extended SA method to two different molecular stochastic systems, a standard Lennard-Jones fluid and an all-atom methane liquid, and compare the obtained parameter sensitivities with parameter sensitivities on three popular and well-studied observable functions, namely, the radial distribution function, the mean squared displacement, and the pressure. Results show that the RER-based sensitivities are highly correlated with the observable-based sensitivities.

Rigid multibody system dynamics with uncertain rigid bodies

Energy Technology Data Exchange (ETDEWEB)

Batou, A., E-mail: anas.batou@univ-paris-est.fr; Soize, C., E-mail: christian.soize@univ-paris-est.fr [Universite Paris-Est, Laboratoire Modelisation et Simulation Multi Echelle, MSME UMR 8208 CNRS (France)

2012-03-15

This paper is devoted to the construction of a probabilistic model of uncertain rigid bodies for multibody system dynamics. We first construct a stochastic model of an uncertain rigid body by replacing the mass, the center of mass, and the tensor of inertia by random variables. The prior probability distributions of the stochastic model are constructed using the maximum entropy principle under the constraints defined by the available information. The generators of independent realizations corresponding to the prior probability distribution of these random quantities are further developed. Then several uncertain rigid bodies can be linked to each other in order to calculate the random response of a multibody dynamical system. An application is proposed to illustrate the theoretical development.

Noise-induced temporal dynamics in Turing systems

KAUST Repository

Schumacher, Linus J.

2013-04-25

We examine the ability of intrinsic noise to produce complex temporal dynamics in Turing pattern formation systems, with particular emphasis on the Schnakenberg kinetics. Using power spectral methods, we characterize the behavior of the system using stochastic simulations at a wide range of points in parameter space and compare with analytical approximations. Specifically, we investigate whether polarity switching of stochastic patterns occurs at a defined frequency. We find that it can do so in individual realizations of a stochastic simulation, but that the frequency is not defined consistently across realizations in our samples of parameter space. Further, we examine the effect of noise on deterministically predicted traveling waves and find them increased in amplitude and decreased in speed. © 2013 American Physical Society.

Renormalization group method in the theory of dynamical systems

International Nuclear Information System (INIS)

Sinai, Y.G.; Khanin, K.M.

1988-01-01

One of the most important events in the theory of dynamical systems for the last decade has become a wide penetration of ideas and renormalization group methods (RG) into this traditional field of mathematical physics. RG-method has been one of the main tools in statistical physics and it has proved to be rather effective while solving problems of the theory of dynamical systems referring to new types of bifurcations (see further). As in statistical mechanics the application of the RG-method is of great interest in the neighborhood of the critical point concerning the order-chaos transition. First the RG-method was applied in the pioneering papers dedicated to the appearance of a stochastical regime as a result of infinite sequences of period doubling bifurcations. At present this stochasticity mechanism is the most studied one and many papers deal with it. The study of the so-called intermittency phenomenon was the next example of application of the RG-method, i.e. the study of such a situation where the domains of the stochastical and regular behavior do alternate along a trajectory of the dynamical system

Stochastic Simulation of Cardiac Ventricular Myocyte Calcium Dynamics and Waves

OpenAIRE

Tuan, Hoang-Trong Minh; Williams, George S. B.; Chikando, Aristide C.; Sobie, Eric A.; Lederer, W. Jonathan; Jafri, M. Saleet

2011-01-01

A three dimensional model of calcium dynamics in the rat ventricular myocyte was developed to study the mechanism of calcium homeostasis and pathological calcium dynamics during calcium overload. The model contains 20,000 calcium release units (CRUs) each containing 49 ryanodine receptors. The model simulates calcium sparks with a realistic spontaneous calcium spark rate. It suggests that in addition to the calcium spark-based leak, there is an invisible calcium leak caused by the stochastic ...

Interactive macroeconomics stochastic aggregate dynamics with heterogeneous and interacting agents

CERN Document Server

Di Guilmi, Corrado

2017-01-01

One of the major problems of macroeconomic theory is the way in which the people exchange goods in decentralized market economies. There are major disagreements among macroeconomists regarding tools to influence required outcomes. Since the mainstream efficient market theory fails to provide an internal coherent framework, there is a need for an alternative theory. The book provides an innovative approach for the analysis of agent based models, populated by the heterogeneous and interacting agents in the field of financial fragility. The text is divided in two parts; the first presents analytical developments of stochastic aggregation and macro-dynamics inference methods. The second part introduces macroeconomic models of financial fragility for complex systems populated by heterogeneous and interacting agents. The concepts of financial fragility and macroeconomic dynamics are explained in detail in separate chapters. The statistical physics approach is applied to explain theories of macroeconomic modelling a...

A generic methodology for the optimisation of sewer systems using stochastic programming and self-optimizing control

DEFF Research Database (Denmark)

Maurico-Iglesias, Miguel; Castro, Ignacio Montero; Mollerup, Ane Loft

2015-01-01

. Such controller is aimed at keeping the system close to the optimal performance, thanks to an optimal selection of controlled variables. The definition of an optimal performance was carried out by a two-stage optimisation (stochastic and deterministic) to take into account both the overflow during the current......The design of sewer system control is a complex task given the large size of the sewer networks, the transient dynamics of the water flow and the stochastic nature of rainfall. This contribution presents a generic methodology for the design of a self-optimising controller in sewer systems...

Stochastic dynamic stiffness of surface footing for offshore wind turbines

DEFF Research Database (Denmark)

Vahdatirad, Mohammadjavad; Andersen, Lars Vabbersgaard; Ibsen, Lars Bo

2014-01-01

Highlights •This study concerns the stochastic dynamic stiffness of foundations for large offshore wind turbines. •A simple model of wind turbine structure with equivalent coupled springs at the base is utilized. •The level of uncertainties is quantified through a sensitivity analysis. •Estimation...

Stochasticity in the Josephson map

International Nuclear Information System (INIS)

Nomura, Y.; Ichikawa, Y.H.; Filippov, A.T.

1996-04-01

The Josephson map describes nonlinear dynamics of systems characterized by standard map with the uniform external bias superposed. The intricate structures of the phase space portrait of the Josephson map are examined on the basis of the tangent map associated with the Josephson map. Numerical observation of the stochastic diffusion in the Josephson map is examined in comparison with the renormalized diffusion coefficient calculated by the method of characteristic function. The global stochasticity of the Josephson map occurs at the values of far smaller stochastic parameter than the case of the standard map. (author)

The global dynamics for a stochastic SIS epidemic model with isolation

Science.gov (United States)

Chen, Yiliang; Wen, Buyu; Teng, Zhidong

2018-02-01

In this paper, we investigate the dynamical behavior for a stochastic SIS epidemic model with isolation which is as an important strategy for the elimination of infectious diseases. It is assumed that the stochastic effects manifest themselves mainly as fluctuation in the transmission coefficient, the death rate and the proportional coefficient of the isolation of infective. It is shown that the extinction and persistence in the mean of the model are determined by a threshold value R0S . That is, if R0S 1, then the disease is stochastic persistent in the means with probability one. Furthermore, the existence of a unique stationary distribution is discussed, and the sufficient conditions are established by using the Lyapunov function method. Finally, some numerical examples are carried out to confirm the analytical results.

Threshold Dynamics in Stochastic SIRS Epidemic Models with Nonlinear Incidence and Vaccination

Directory of Open Access Journals (Sweden)

Lei Wang

2017-01-01

Full Text Available In this paper, the dynamical behaviors for a stochastic SIRS epidemic model with nonlinear incidence and vaccination are investigated. In the models, the disease transmission coefficient and the removal rates are all affected by noise. Some new basic properties of the models are found. Applying these properties, we establish a series of new threshold conditions on the stochastically exponential extinction, stochastic persistence, and permanence in the mean of the disease with probability one for the models. Furthermore, we obtain a sufficient condition on the existence of unique stationary distribution for the model. Finally, a series of numerical examples are introduced to illustrate our main theoretical results and some conjectures are further proposed.

Dynamic Asset Allocation with Stochastic Income and Interest Rates

DEFF Research Database (Denmark)

Munk, Claus; Sørensen, Carsten

2010-01-01

We solve for optimal portfolios when interest rates and labor income are stochastic with the expected income growth being affine in the short-term interest rate in order to encompass business cycle variations in wages. Our calibration based on the Panel Study of Income Dynamics (PSID) data supports...

Linear dynamical quantum systems analysis, synthesis, and control

CERN Document Server

Nurdin, Hendra I

2017-01-01

This monograph provides an in-depth treatment of the class of linear-dynamical quantum systems. The monograph presents a detailed account of the mathematical modeling of these systems using linear algebra and quantum stochastic calculus as the main tools for a treatment that emphasizes a system-theoretic point of view and the control-theoretic formulations of quantum versions of familiar problems from the classical (non-quantum) setting, including estimation and filtering, realization theory, and feedback control. Both measurement-based feedback control (i.e., feedback control by a classical system involving a continuous-time measurement process) and coherent feedback control (i.e., feedback control by another quantum system without the intervention of any measurements in the feedback loop) are treated. Researchers and graduates studying systems and control theory, quantum probability and stochastics or stochastic control whether from backgrounds in mechanical or electrical engineering or applied mathematics ...

Time-ordered product expansions for computational stochastic system biology

International Nuclear Information System (INIS)

Mjolsness, Eric

2013-01-01

The time-ordered product framework of quantum field theory can also be used to understand salient phenomena in stochastic biochemical networks. It is used here to derive Gillespie’s stochastic simulation algorithm (SSA) for chemical reaction networks; consequently, the SSA can be interpreted in terms of Feynman diagrams. It is also used here to derive other, more general simulation and parameter-learning algorithms including simulation algorithms for networks of stochastic reaction-like processes operating on parameterized objects, and also hybrid stochastic reaction/differential equation models in which systems of ordinary differential equations evolve the parameters of objects that can also undergo stochastic reactions. Thus, the time-ordered product expansion can be used systematically to derive simulation and parameter-fitting algorithms for stochastic systems. (paper)

Optimal and centralized reservoir management for drought and flood protection via Stochastic Dual Dynamic Programming on the Upper Seine-Aube River system

Science.gov (United States)

Chiavico, Mattia; Raso, Luciano; Dorchies, David; Malaterre, Pierre-Olivier

2015-04-01

Seine river region is an extremely important logistic and economic junction for France and Europe. The hydraulic protection of most part of the region relies on four controlled reservoirs, managed by EPTB Seine-Grands Lacs. Presently, reservoirs operation is not centrally coordinated, and release rules are based on empirical filling curves. In this study, we analyze how a centralized release policy can face flood and drought risks, optimizing water system efficiency. The optimal and centralized decisional problem is solved by Stochastic Dual Dynamic Programming (SDDP) method, minimizing an operational indicator for each planning objective. SDDP allows us to include into the system: 1) the hydrological discharge, specifically a stochastic semi-distributed auto-regressive model, 2) the hydraulic transfer model, represented by a linear lag and route model, and 3) reservoirs and diversions. The novelty of this study lies on the combination of reservoir and hydraulic models in SDDP for flood and drought protection problems. The study case covers the Seine basin until the confluence with Aube River: this system includes two reservoirs, the city of Troyes, and the Nuclear power plant of Nogent-Sur-Seine. The conflict between the interests of flood protection, drought protection, water use and ecology leads to analyze the environmental system in a Multi-Objective perspective.

Stochastic Ocean Predictions with Dynamically-Orthogonal Primitive Equations

Science.gov (United States)

Subramani, D. N.; Haley, P., Jr.; Lermusiaux, P. F. J.

2017-12-01

The coastal ocean is a prime example of multiscale nonlinear fluid dynamics. Ocean fields in such regions are complex and intermittent with unstationary heterogeneous statistics. Due to the limited measurements, there are multiple sources of uncertainties, including the initial conditions, boundary conditions, forcing, parameters, and even the model parameterizations and equations themselves. For efficient and rigorous quantification and prediction of these uncertainities, the stochastic Dynamically Orthogonal (DO) PDEs for a primitive equation ocean modeling system with a nonlinear free-surface are derived and numerical schemes for their space-time integration are obtained. Detailed numerical studies with idealized-to-realistic regional ocean dynamics are completed. These include consistency checks for the numerical schemes and comparisons with ensemble realizations. As an illustrative example, we simulate the 4-d multiscale uncertainty in the Middle Atlantic/New York Bight region during the months of Jan to Mar 2017. To provide intitial conditions for the uncertainty subspace, uncertainties in the region were objectively analyzed using historical data. The DO primitive equations were subsequently integrated in space and time. The probability distribution function (pdf) of the ocean fields is compared to in-situ, remote sensing, and opportunity data collected during the coincident POSYDON experiment. Results show that our probabilistic predictions had skill and are 3- to 4- orders of magnitude faster than classic ensemble schemes.

Pesin’s entropy formula for stochastic flows of diffeomorphisms

Institute of Scientific and Technical Information of China (English)

刘培东

1996-01-01

Pesin’s entropy formula relating entropy and Lyapunov exponents within the context of random dynamical systems generated by (discrete or continuous) stochastic flows of diffeomorphisms (including solution flows of stochastic differential equations on manifolds) is proved.

Fractional noise destroys or induces a stochastic bifurcation

Energy Technology Data Exchange (ETDEWEB)

Yang, Qigui, E-mail: qgyang@scut.edu.cn [School of Sciences, South China University of Technology, Guangzhou 510640 (China); Zeng, Caibin, E-mail: zeng.cb@mail.scut.edu.cn [School of Sciences, South China University of Technology, Guangzhou 510640 (China); School of Automation Science and Engineering, South China University of Technology, Guangzhou 510640 (China); Wang, Cong, E-mail: wangcong@scut.edu.cn [School of Automation Science and Engineering, South China University of Technology, Guangzhou 510640 (China)

2013-12-15

Little seems to be known about the stochastic bifurcation phenomena of non-Markovian systems. Our intention in this paper is to understand such complex dynamics by a simple system, namely, the Black-Scholes model driven by a mixed fractional Brownian motion. The most interesting finding is that the multiplicative fractional noise not only destroys but also induces a stochastic bifurcation under some suitable conditions. So it opens a possible way to explore the theory of stochastic bifurcation in the non-Markovian framework.

Fault Detection for Wireless Networked Control Systems with Stochastic Switching Topology and Time Delay

Directory of Open Access Journals (Sweden)

Pengfei Guo

2014-01-01

Full Text Available This paper deals with the fault detection problem for a class of discrete-time wireless networked control systems described by switching topology with uncertainties and disturbances. System states of each individual node are affected not only by its own measurements, but also by other nodes’ measurements according to a certain network topology. As the topology of system can be switched in a stochastic way, we aim to design H∞ fault detection observers for nodes in the dynamic time-delay systems. By using the Lyapunov method and stochastic analysis techniques, sufficient conditions are acquired to guarantee the existence of the filters satisfying the H∞ performance constraint, and observer gains are derived by solving linear matrix inequalities. Finally, an illustrated example is provided to verify the effectiveness of the theoretical results.

Stochastic Control - External Models

DEFF Research Database (Denmark)

Poulsen, Niels Kjølstad

2005-01-01

This note is devoted to control of stochastic systems described in discrete time. We are concerned with external descriptions or transfer function model, where we have a dynamic model for the input output relation only (i.e.. no direct internal information). The methods are based on LTI systems...

Distributed parallel computing in stochastic modeling of groundwater systems.

Science.gov (United States)

Dong, Yanhui; Li, Guomin; Xu, Haizhen

2013-03-01

Stochastic modeling is a rapidly evolving, popular approach to the study of the uncertainty and heterogeneity of groundwater systems. However, the use of Monte Carlo-type simulations to solve practical groundwater problems often encounters computational bottlenecks that hinder the acquisition of meaningful results. To improve the computational efficiency, a system that combines stochastic model generation with MODFLOW-related programs and distributed parallel processing is investigated. The distributed computing framework, called the Java Parallel Processing Framework, is integrated into the system to allow the batch processing of stochastic models in distributed and parallel systems. As an example, the system is applied to the stochastic delineation of well capture zones in the Pinggu Basin in Beijing. Through the use of 50 processing threads on a cluster with 10 multicore nodes, the execution times of 500 realizations are reduced to 3% compared with those of a serial execution. Through this application, the system demonstrates its potential in solving difficult computational problems in practical stochastic modeling. © 2012, The Author(s). Groundwater © 2012, National Ground Water Association.

Optimal Strategy for Integrated Dynamic Inventory Control and Supplier Selection in Unknown Environment via Stochastic Dynamic Programming

International Nuclear Information System (INIS)

Sutrisno; Widowati; Solikhin

2016-01-01

In this paper, we propose a mathematical model in stochastic dynamic optimization form to determine the optimal strategy for an integrated single product inventory control problem and supplier selection problem where the demand and purchasing cost parameters are random. For each time period, by using the proposed model, we decide the optimal supplier and calculate the optimal product volume purchased from the optimal supplier so that the inventory level will be located at some point as close as possible to the reference point with minimal cost. We use stochastic dynamic programming to solve this problem and give several numerical experiments to evaluate the model. From the results, for each time period, the proposed model was generated the optimal supplier and the inventory level was tracked the reference point well. (paper)

A stochastic agent-based model of pathogen propagation in dynamic multi-relational social networks

Science.gov (United States)

Khan, Bilal; Dombrowski, Kirk; Saad, Mohamed

2015-01-01

We describe a general framework for modeling and stochastic simulation of epidemics in realistic dynamic social networks, which incorporates heterogeneity in the types of individuals, types of interconnecting risk-bearing relationships, and types of pathogens transmitted across them. Dynamism is supported through arrival and departure processes, continuous restructuring of risk relationships, and changes to pathogen infectiousness, as mandated by natural history; dynamism is regulated through constraints on the local agency of individual nodes and their risk behaviors, while simulation trajectories are validated using system-wide metrics. To illustrate its utility, we present a case study that applies the proposed framework towards a simulation of HIV in artificial networks of intravenous drug users (IDUs) modeled using data collected in the Social Factors for HIV Risk survey. PMID:25859056

A study on stochastic resonance of one-dimensional bistable system in the neighborhood of bifurcation point with the moment method

Energy Technology Data Exchange (ETDEWEB)

Zhang Guangjun [State Key Laboratory of Mechanical Structural Strength and Vibration, School of Architectural Engineering and Mechanics, Xi' an Jiaotong University, Xi' an, Shaanxi (China); Xu Jianxue [State Key Laboratory of Mechanical Structural Strength and Vibration, School of Architectural Engineering and Mechanics, Xi' an Jiaotong University, Xi' an, Shaanxi (China)] e-mail: jxxu@mail.xjtu.edu.cn

2006-02-01

This paper analyzes the stochastic resonance induced by a novel transition of one-dimensional bistable system in the neighborhood of bifurcation point with the method of moment, which refer to the transition of system motion among a potential well of stable fixed point before bifurcation of original system and double-well potential of two coexisting stable fixed points after original system bifurcation at the presence of internal noise. The results show: the semi-analytical result of stochastic resonance of one-dimensional bistable system in the neighborhood of bifurcation point may be obtained, and the semi-analytical result is in accord with the one of Monte Carlo simulation qualitatively, the occurrence of stochastic resonance is related to the bifurcation of noisy nonlinear dynamical system moment equations, which induce the transfer of energy of ensemble average (Ex) of system response in each frequency component and make the energy of ensemble average of system response concentrate on the frequency of input signal, stochastic resonance occurs.

A study on stochastic resonance of one-dimensional bistable system in the neighborhood of bifurcation point with the moment method

International Nuclear Information System (INIS)

Zhang Guangjun; Xu Jianxue

2006-01-01

This paper analyzes the stochastic resonance induced by a novel transition of one-dimensional bistable system in the neighborhood of bifurcation point with the method of moment, which refer to the transition of system motion among a potential well of stable fixed point before bifurcation of original system and double-well potential of two coexisting stable fixed points after original system bifurcation at the presence of internal noise. The results show: the semi-analytical result of stochastic resonance of one-dimensional bistable system in the neighborhood of bifurcation point may be obtained, and the semi-analytical result is in accord with the one of Monte Carlo simulation qualitatively, the occurrence of stochastic resonance is related to the bifurcation of noisy nonlinear dynamical system moment equations, which induce the transfer of energy of ensemble average (Ex) of system response in each frequency component and make the energy of ensemble average of system response concentrate on the frequency of input signal, stochastic resonance occurs

Optimization of stochastic discrete systems and control on complex networks computational networks

CERN Document Server

Lozovanu, Dmitrii

2014-01-01

This book presents the latest findings on stochastic dynamic programming models and on solving optimal control problems in networks. It includes the authors' new findings on determining the optimal solution of discrete optimal control problems in networks and on solving game variants of Markov decision problems in the context of computational networks. First, the book studies the finite state space of Markov processes and reviews the existing methods and algorithms for determining the main characteristics in Markov chains, before proposing new approaches based on dynamic programming and combinatorial methods. Chapter two is dedicated to infinite horizon stochastic discrete optimal control models and Markov decision problems with average and expected total discounted optimization criteria, while Chapter three develops a special game-theoretical approach to Markov decision processes and stochastic discrete optimal control problems. In closing, the book's final chapter is devoted to finite horizon stochastic con...

Stochastic equations for complex systems theoretical and computational topics

CERN Document Server

Bessaih, Hakima

2015-01-01

Mathematical analyses and computational predictions of the behavior of complex systems are needed to effectively deal with weather and climate predictions, for example, and the optimal design of technical processes. Given the random nature of such systems and the recognized relevance of randomness, the equations used to describe such systems usually need to involve stochastics.  The basic goal of this book is to introduce the mathematics and application of stochastic equations used for the modeling of complex systems. A first focus is on the introduction to different topics in mathematical analysis. A second focus is on the application of mathematical tools to the analysis of stochastic equations. A third focus is on the development and application of stochastic methods to simulate turbulent flows as seen in reality.  This book is primarily oriented towards mathematics and engineering PhD students, young and experienced researchers, and professionals working in the area of stochastic differential equations ...

Optimal Control and Optimization of Stochastic Supply Chain Systems

CERN Document Server

Song, Dong-Ping

2013-01-01

Optimal Control and Optimization of Stochastic Supply Chain Systems examines its subject in the context of the presence of a variety of uncertainties. Numerous examples with intuitive illustrations and tables are provided, to demonstrate the structural characteristics of the optimal control policies in various stochastic supply chains and to show how to make use of these characteristics to construct easy-to-operate sub-optimal policies.                 In Part I, a general introduction to stochastic supply chain systems is provided. Analytical models for various stochastic supply chain systems are formulated and analysed in Part II. In Part III the structural knowledge of the optimal control policies obtained in Part II is utilized to construct easy-to-operate sub-optimal control policies for various stochastic supply chain systems accordingly. Finally, Part IV discusses the optimisation of threshold-type control policies and their robustness. A key feature of the book is its tying together of ...

Metastable states and quasicycles in a stochastic Wilson-Cowan model of neuronal population dynamics

KAUST Repository

Bressloff, Paul C.

2010-01-01

We analyze a stochastic model of neuronal population dynamics with intrinsic noise. In the thermodynamic limit N→∞, where N determines the size of each population, the dynamics is described by deterministic Wilson-Cowan equations. On the other hand

Comparison of stochastic resonance in static and dynamical nonlinearities

International Nuclear Information System (INIS)

Ma, Yumei; Duan, Fabing

2014-01-01

We compare the stochastic resonance (SR) effects in parallel arrays of static and dynamical nonlinearities via the measure of output signal-to-noise ratio (SNR). For a received noisy periodic signal, parallel arrays of both static and dynamical nonlinearities can enhance the output SNR by optimizing the internal noise level. The static nonlinearity is easily implementable, while the dynamical nonlinearity has more parameters to be tuned, at the risk of not exploiting the beneficial role of internal noise components. It is of interest to note that, for an input signal buried in the external Laplacian noise, we show that the dynamical nonlinearity is superior to the static nonlinearity in obtaining a better output SNR. This characteristic is assumed to be closely associated with the kurtosis of noise distribution. - Highlights: • Comparison of SR effects in arrays of both static and dynamical nonlinearities. • Static nonlinearity is easily implementable for the SNR enhancement. • Dynamical nonlinearity yields a better output SNR for external Laplacian noise

Asymptotic analysis for functional stochastic differential equations

CERN Document Server

Bao, Jianhai; Yuan, Chenggui

2016-01-01

This brief treats dynamical systems that involve delays and random disturbances. The study is motivated by a wide variety of systems in real life in which random noise has to be taken into consideration and the effect of delays cannot be ignored. Concentrating on such systems that are described by functional stochastic differential equations, this work focuses on the study of large time behavior, in particular, ergodicity. This brief is written for probabilists, applied mathematicians, engineers, and scientists who need to use delay systems and functional stochastic differential equations in their work. Selected topics from the brief can also be used in a graduate level topics course in probability and stochastic processes.

Stochastic Neural Field Theory and the System-Size Expansion

KAUST Repository

Bressloff, Paul C.

2010-01-01

We analyze a master equation formulation of stochastic neurodynamics for a network of synaptically coupled homogeneous neuronal populations each consisting of N identical neurons. The state of the network is specified by the fraction of active or spiking neurons in each population, and transition rates are chosen so that in the thermodynamic or deterministic limit (N → ∞) we recover standard activity-based or voltage-based rate models. We derive the lowest order corrections to these rate equations for large but finite N using two different approximation schemes, one based on the Van Kampen system-size expansion and the other based on path integral methods. Both methods yield the same series expansion of the moment equations, which at O(1/N) can be truncated to form a closed system of equations for the first-and second-order moments. Taking a continuum limit of the moment equations while keeping the system size N fixed generates a system of integrodifferential equations for the mean and covariance of the corresponding stochastic neural field model. We also show how the path integral approach can be used to study large deviation or rare event statistics underlying escape from the basin of attraction of a stable fixed point of the mean-field dynamics; such an analysis is not possible using the system-size expansion since the latter cannot accurately determine exponentially small transitions. © by SIAM.

Stochastic light-cone CTMRG: a new DMRG approach to stochastic models 02.50.Ey Stochastic processes; 64.60.Ht Dynamic critical phenomena; 02.70.-c Computational techniques; 05.10.Cc Renormalization group methods;

CERN Document Server

Kemper, A; Nishino, T; Schadschneider, A; Zittartz, J

2003-01-01

We develop a new variant of the recently introduced stochastic transfer matrix DMRG which we call stochastic light-cone corner-transfer-matrix DMRG (LCTMRG). It is a numerical method to compute dynamic properties of one-dimensional stochastic processes. As suggested by its name, the LCTMRG is a modification of the corner-transfer-matrix DMRG, adjusted by an additional causality argument. As an example, two reaction-diffusion models, the diffusion-annihilation process and the branch-fusion process are studied and compared with exact data and Monte Carlo simulations to estimate the capability and accuracy of the new method. The number of possible Trotter steps of more than 10 sup 5 shows a considerable improvement on the old stochastic TMRG algorithm.

Manipulating potential wells in Logical Stochastic Resonance to obtain XOR logic

International Nuclear Information System (INIS)

Storni, Remo; Ando, Hiroyasu; Aihara, Kazuyuki; Murali, K.; Sinha, Sudeshna

2012-01-01

Logical Stochastic Resonance (LSR) is the application of Stochastic Resonance to logic computation, namely the phenomenon where a nonlinear system driven by weak signals representing logic inputs, under optimal noise, can yield logic outputs. We extend the existing results, obtained in the context of bistable systems, to multi-stable dynamical systems, allowing us to obtain XOR logic, in addition to the AND (NAND) and OR (NOR) logic observed in earlier studies. This strategy widens the scope of LSR from the application point of view, as XOR forms the basis of ubiquitous bit-by-bit addition, and conceptually, showing the ability to yield non-monotonic input–output logic associations. -- Highlights: ► We generalize Logical Stochastic Resonance from bistable to multi-stable systems. ► We propose a tristable dynamical system formed of piecewise linear functions. ► The system can correctly reproduce XOR logic behavior using the LSR principle. ► The system yields different logic behavior without the need to change the dynamics.

Bifurcation Regulations Governed by Delay Self-Control Feedback in a Stochastic Birhythmic System

Science.gov (United States)

Ma, Zhidan; Ning, Lijuan

2017-12-01

We aim to investigate bifurcation behaviors in a stochastic birhythmic van der Pol (BVDP) system subjected to delay self-control feedback. First, the harmonic approximation is adopted to drive the delay self-control feedback to state variables without delay. Then, Fokker-Planck-Kolmogorov (FPK) equation and stationary probability density function (SPDF) for amplitude are obtained by applying stochastic averaging method. Finally, dynamical scenarios of the change of delay self-control feedback as well as noise that markedly influence bifurcation performance are observed. It is found that: the big feedback strength and delay will suppress the large amplitude limit cycle (LC) while the relatively big noise strength facilitates the large amplitude LC, which imply the proposed regulation strategies are feasible. Interestingly enough, the inner LC is never destroyed due to noise. Furthermore, the validity of analytical results was verified by Monte Carlo simulation of the dynamics.

Entropy Evolution and Uncertainty Estimation with Dynamical Systems

Directory of Open Access Journals (Sweden)

X. San Liang

2014-06-01

Full Text Available This paper presents a comprehensive introduction and systematic derivation of the evolutionary equations for absolute entropy H and relative entropy D, some of which exist sporadically in the literature in different forms under different subjects, within the framework of dynamical systems. In general, both H and D are dissipated, and the dissipation bears a form reminiscent of the Fisher information; in the absence of stochasticity, dH/dt is connected to the rate of phase space expansion, and D stays invariant, i.e., the separation of two probability density functions is always conserved. These formulas are validated with linear systems, and put to application with the Lorenz system and a large-dimensional stochastic quasi-geostrophic flow problem. In the Lorenz case, H falls at a constant rate with time, implying that H will eventually become negative, a situation beyond the capability of the commonly used computational technique like coarse-graining and bin counting. For the stochastic flow problem, it is first reduced to a computationally tractable low-dimensional system, using a reduced model approach, and then handled through ensemble prediction. Both the Lorenz system and the stochastic flow system are examples of self-organization in the light of uncertainty reduction. The latter particularly shows that, sometimes stochasticity may actually enhance the self-organization process.

Nonlinear and Stochastic Dynamics in the Heart

Science.gov (United States)

Qu, Zhilin; Hu, Gang; Garfinkel, Alan; Weiss, James N.

2014-01-01

In a normal human life span, the heart beats about 2 to 3 billion times. Under diseased conditions, a heart may lose its normal rhythm and degenerate suddenly into much faster and irregular rhythms, called arrhythmias, which may lead to sudden death. The transition from a normal rhythm to an arrhythmia is a transition from regular electrical wave conduction to irregular or turbulent wave conduction in the heart, and thus this medical problem is also a problem of physics and mathematics. In the last century, clinical, experimental, and theoretical studies have shown that dynamical theories play fundamental roles in understanding the mechanisms of the genesis of the normal heart rhythm as well as lethal arrhythmias. In this article, we summarize in detail the nonlinear and stochastic dynamics occurring in the heart and their links to normal cardiac functions and arrhythmias, providing a holistic view through integrating dynamics from the molecular (microscopic) scale, to the organelle (mesoscopic) scale, to the cellular, tissue, and organ (macroscopic) scales. We discuss what existing problems and challenges are waiting to be solved and how multi-scale mathematical modeling and nonlinear dynamics may be helpful for solving these problems. PMID:25267872

Nonlinear and stochastic dynamics in the heart

Energy Technology Data Exchange (ETDEWEB)

Qu, Zhilin, E-mail: zqu@mednet.ucla.edu [Department of Medicine (Cardiology), David Geffen School of Medicine, University of California, Los Angeles, CA 90095 (United States); Hu, Gang [Department of Physics, Beijing Normal University, Beijing 100875 (China); Garfinkel, Alan [Department of Medicine (Cardiology), David Geffen School of Medicine, University of California, Los Angeles, CA 90095 (United States); Department of Integrative Biology and Physiology, University of California, Los Angeles, CA 90095 (United States); Weiss, James N. [Department of Medicine (Cardiology), David Geffen School of Medicine, University of California, Los Angeles, CA 90095 (United States); Department of Physiology, David Geffen School of Medicine, University of California, Los Angeles, CA 90095 (United States)

2014-10-10

In a normal human life span, the heart beats about 2–3 billion times. Under diseased conditions, a heart may lose its normal rhythm and degenerate suddenly into much faster and irregular rhythms, called arrhythmias, which may lead to sudden death. The transition from a normal rhythm to an arrhythmia is a transition from regular electrical wave conduction to irregular or turbulent wave conduction in the heart, and thus this medical problem is also a problem of physics and mathematics. In the last century, clinical, experimental, and theoretical studies have shown that dynamical theories play fundamental roles in understanding the mechanisms of the genesis of the normal heart rhythm as well as lethal arrhythmias. In this article, we summarize in detail the nonlinear and stochastic dynamics occurring in the heart and their links to normal cardiac functions and arrhythmias, providing a holistic view through integrating dynamics from the molecular (microscopic) scale, to the organelle (mesoscopic) scale, to the cellular, tissue, and organ (macroscopic) scales. We discuss what existing problems and challenges are waiting to be solved and how multi-scale mathematical modeling and nonlinear dynamics may be helpful for solving these problems.

Nonlinear and stochastic dynamics in the heart

International Nuclear Information System (INIS)

Qu, Zhilin; Hu, Gang; Garfinkel, Alan; Weiss, James N.

2014-01-01

In a normal human life span, the heart beats about 2–3 billion times. Under diseased conditions, a heart may lose its normal rhythm and degenerate suddenly into much faster and irregular rhythms, called arrhythmias, which may lead to sudden death. The transition from a normal rhythm to an arrhythmia is a transition from regular electrical wave conduction to irregular or turbulent wave conduction in the heart, and thus this medical problem is also a problem of physics and mathematics. In the last century, clinical, experimental, and theoretical studies have shown that dynamical theories play fundamental roles in understanding the mechanisms of the genesis of the normal heart rhythm as well as lethal arrhythmias. In this article, we summarize in detail the nonlinear and stochastic dynamics occurring in the heart and their links to normal cardiac functions and arrhythmias, providing a holistic view through integrating dynamics from the molecular (microscopic) scale, to the organelle (mesoscopic) scale, to the cellular, tissue, and organ (macroscopic) scales. We discuss what existing problems and challenges are waiting to be solved and how multi-scale mathematical modeling and nonlinear dynamics may be helpful for solving these problems

Stochastic population dynamics of a montane ground-dwelling squirrel.

Science.gov (United States)

Hostetler, Jeffrey A; Kneip, Eva; Van Vuren, Dirk H; Oli, Madan K

2012-01-01

Understanding the causes and consequences of population fluctuations is a central goal of ecology. We used demographic data from a long-term (1990-2008) study and matrix population models to investigate factors and processes influencing the dynamics and persistence of a golden-mantled ground squirrel (Callospermophilus lateralis) population, inhabiting a dynamic subalpine habitat in Colorado, USA. The overall deterministic population growth rate λ was 0.94±SE 0.05 but it varied widely over time, ranging from 0.45±0.09 in 2006 to 1.50±0.12 in 2003, and was below replacement (λbounce back from low densities and prevented extinction. These results suggest that dynamics and persistence of our study population are determined synergistically by density-dependence, stochastic forces, and immigration.

Advances in analysis and control of timedelayed dynamical systems

CERN Document Server

Sun, Jianqiao

2013-01-01

Analysis and control of timedelayed systems have been applied in a wide range of applications, ranging from mechanical, control, economic, to biological systems. Over the years, there has been a steady stream of interest in timedelayed dynamic systems, this book takes a snap shot of recent research from the world leading experts in analysis and control of dynamic systems with time delay to provide a bird's eye view of its development. The topics covered in this book include solution methods, stability analysis and control of periodic dynamic systems with time delay, bifurcations, stochastic dy

Volatile decision dynamics: experiments, stochastic description, intermittency control and traffic optimization

Science.gov (United States)

Helbing, Dirk; Schönhof, Martin; Kern, Daniel

2002-06-01

The coordinated and efficient distribution of limited resources by individual decisions is a fundamental, unsolved problem. When individuals compete for road capacities, time, space, money, goods, etc, they normally make decisions based on aggregate rather than complete information, such as TV news or stock market indices. In related experiments, we have observed a volatile decision dynamics and far-from-optimal payoff distributions. We have also identified methods of information presentation that can considerably improve the overall performance of the system. In order to determine optimal strategies of decision guidance by means of user-specific recommendations, a stochastic behavioural description is developed. These strategies manage to increase the adaptibility to changing conditions and to reduce the deviation from the time-dependent user equilibrium, thereby enhancing the average and individual payoffs. Hence, our guidance strategies can increase the performance of all users by reducing overreaction and stabilizing the decision dynamics. These results are highly significant for predicting decision behaviour, for reaching optimal behavioural distributions by decision support systems and for information service providers. One of the promising fields of application is traffic optimization.

Stochastic resonance and vibrational resonance in an excitable system: The golden mean barrier

International Nuclear Information System (INIS)

Stan, Cristina; Cristescu, C.P.; Alexandroaei, D.; Agop, M.

2009-01-01

We report on stochastic resonance and vibrational resonance in an electric charge double layer configuration as usually found in electrical discharges, biological cell membranes, chemical systems and nanostructures. The experiment and numerical computation show the existence of a barrier expressible in terms of the golden mean above which the two phenomena do not take place. We consider this as new evidence for the importance of the golden mean criticality in the oscillatory dynamics, in agreement with El Naschie's E-infinity theory. In our experiment, the dynamics of a charge double layer generated in the inter-anode space of a twin electrical discharge is investigated under noise-harmonic and harmonic-harmonic perturbations. In the first case, a Gaussian noise can enhance the response of the system to a weak injected periodic signal, a clear mark of stochastic resonance. In the second case, similar enhancement can appear if the noise is replaced by a harmonic perturbation with a frequency much higher than the frequency of the weak oscillation. The amplitude of the low frequency oscillation shows a maximum versus the amplitude of the high frequency perturbation demonstrating vibrational resonance. In order to model these dynamics, we derived an excitable system by modifying a biased van der Pol oscillator. The computational study considers the behaviour of this system under the same types of perturbation as in the experimental investigations and is found to give consistent results in both situations.

Stochastic transport processes in discrete biological systems

CERN Document Server

Frehland, Eckart

1982-01-01

These notes are in part based on a course for advanced students in the applications of stochastic processes held in 1978 at the University of Konstanz. These notes contain the results of re­ cent studies on the stochastic description of ion transport through biological membranes. In particular, they serve as an introduction to an unified theory of fluctuations in complex biological transport systems. We emphasize that the subject of this volume is not to introduce the mathematics of stochastic processes but to present a field of theoretical biophysics in which stochastic methods are important. In the last years the study of membrane noise has become an important method in biophysics. Valuable information on the ion transport mechanisms in membranes can be obtained from noise analysis. A number of different processes such as the opening and closing of ion channels have been shown to be sources of the measured current or voltage fluctuations. Bio­ logical 'transport systems can be complex. For example, the tr...

Induced topological pressure for topological dynamical systems

International Nuclear Information System (INIS)

Xing, Zhitao; Chen, Ercai

2015-01-01

In this paper, inspired by the article [J. Jaerisch et al., Stochastics Dyn. 14, 1350016, pp. 1-30 (2014)], we introduce the induced topological pressure for a topological dynamical system. In particular, we prove a variational principle for the induced topological pressure

Stochastic fractional differential equations: Modeling, method and analysis

International Nuclear Information System (INIS)

Pedjeu, Jean-C.; Ladde, Gangaram S.

2012-01-01

By introducing a concept of dynamic process operating under multi-time scales in sciences and engineering, a mathematical model described by a system of multi-time scale stochastic differential equations is formulated. The classical Picard–Lindelöf successive approximations scheme is applied to the model validation problem, namely, existence and uniqueness of solution process. Naturally, this leads to the problem of finding closed form solutions of both linear and nonlinear multi-time scale stochastic differential equations of Itô–Doob type. Finally, to illustrate the scope of ideas and presented results, multi-time scale stochastic models for ecological and epidemiological processes in population dynamic are outlined.

Dynamic-stochastic modeling of snow cover formation on the European territory of Russia

OpenAIRE

A. N. Gelfan; V. M. Moreido

2014-01-01

A dynamic-stochastic model, which combines a deterministic model of snow cover formation with a stochastic weather generator, has been developed. The deterministic snow model describes temporal change of the snow depth, content of ice and liquid water, snow density, snowmelt, sublimation, re-freezing of melt water, and snow metamorphism. The model has been calibrated and validated against the long-term data of snow measurements over the territory of the European Russia. The model showed good ...

Stochastic Dynamic Programming for Three-Echelon Inventory System of Limited Shelf Life Products

Directory of Open Access Journals (Sweden)

Galal Noha M.

2016-01-01

Full Text Available Coordination of inventory decisions within the supply chain is one of the major determinants of its competitiveness in the global market. Products with limited shelf life impose additional challenges in managing the inventory across the supply chain because of the additional wastage costs incurred in case of being stored beyond product’s useful life. This paper presents a stochastic dynamic programming model for inventory replenishment in a serial multi-echelon distribution supply chain. The model considers uncertain stationary discrete demand at the retailer and zero lead time. The objective is to minimize expected total costs across the supply chain echelons, while maintaining a preset service level. The results illustrate that a cost saving of around 17% is achievable due to coordinating inventory decisions across the supply chain.

Modeling and stochastic analysis of dynamic mechanisms of the perception

Science.gov (United States)

Pisarchik, A.; Bashkirtseva, I.; Ryashko, L.

2017-10-01

Modern studies in physiology and cognitive neuroscience consider a noise as an important constructive factor of the brain functionality. Under the adequate noise, the brain can rapidly access different ordered states, and provide decision-making by preventing deadlocks. Bistable dynamic models are often used for the study of the underlying mechanisms of the visual perception. In the present paper, we consider a bistable energy model subject to both additive and parametric noise. Using the catastrophe theory formalism and stochastic sensitivity functions technique, we analyze a response of the equilibria to noise, and study noise-induced transitions between equilibria. We demonstrate and analyse the effect of hysteresis squeezing when the intensity of noise is increased. Stochastic bifurcations connected with the suppression of oscillations by parametric noises are discussed.

Stochastic modeling of lift and drag dynamics to obtain aerodynamic forces with local dynamics on rotor blade under unsteady wind inflow

International Nuclear Information System (INIS)

Luhur, M.R.

2014-01-01

This contribution provides the development of a stochastic lift and drag model for an airfoil FX 79-W-151A under unsteady wind inflow based on wind tunnel measurements. Here we present the integration of the stochastic model into a well-known standard BEM (Blade Element Momentum) model to obtain the corresponding aerodynamic forces on a rotating blade element. The stochastic model is integrated as an alternative to static tabulated data used by classical BEM. The results show that in comparison to classical BEM, the BEM with stochastic approach additionally reflects the local force dynamics and therefore provides more information on aerodynamic forces that can be used by wind turbine simulation codes. (author)

Stochastic Modeling of Lift and Drag Dynamics to Obtain Aerodynamic Forces with Local Dynamics on Rotor Blade under Unsteady Wind Inflow

Directory of Open Access Journals (Sweden)

Muhammad Ramzan Luhur

2014-01-01

Full Text Available This contribution provides the development of a stochastic lift and drag model for an airfoil FX 79-W-151A under unsteady wind inflow based on wind tunnel measurements. Here we present the integration of the stochastic model into a well-known standard BEM (Blade Element Momentum model to obtain the corresponding aerodynamic forces on a rotating blade element. The stochastic model is integrated as an alternative to static tabulated data used by classical BEM. The results show that in comparison to classical BEM, the BEM with stochastic approach additionally reflects the local force dynamics and therefore provides more information on aerodynamic forces that can be used by wind turbine simulation codes

Nonlinear and Complex Dynamics in Real Systems

OpenAIRE

William Barnett; Apostolos Serletis; Demitre Serletis

2005-01-01

This paper was produced for the El-Naschie Symposium on Nonlinear Dynamics in Shanghai in December 2005. In this paper we provide a review of the literature with respect to fluctuations in real systems and chaos. In doing so, we contrast the order and organization hypothesis of real systems with nonlinear chaotic dynamics and discuss some techniques used in distinguishing between stochastic and deterministic behavior. Moreover, we look at the issue of where and when the ideas of chaos could p...

Multiscale Hy3S: Hybrid stochastic simulation for supercomputers

Directory of Open Access Journals (Sweden)

Kaznessis Yiannis N

2006-02-01

Full Text Available Abstract Background Stochastic simulation has become a useful tool to both study natural biological systems and design new synthetic ones. By capturing the intrinsic molecular fluctuations of "small" systems, these simulations produce a more accurate picture of single cell dynamics, including interesting phenomena missed by deterministic methods, such as noise-induced oscillations and transitions between stable states. However, the computational cost of the original stochastic simulation algorithm can be high, motivating the use of hybrid stochastic methods. Hybrid stochastic methods partition the system into multiple subsets and describe each subset as a different representation, such as a jump Markov, Poisson, continuous Markov, or deterministic process. By applying valid approximations and self-consistently merging disparate descriptions, a method can be considerably faster, while retaining accuracy. In this paper, we describe Hy3S, a collection of multiscale simulation programs. Results Building on our previous work on developing novel hybrid stochastic algorithms, we have created the Hy3S software package to enable scientists and engineers to both study and design extremely large well-mixed biological systems with many thousands of reactions and chemical species. We have added adaptive stochastic numerical integrators to permit the robust simulation of dynamically stiff biological systems. In addition, Hy3S has many useful features, including embarrassingly parallelized simulations with MPI; special discrete events, such as transcriptional and translation elongation and cell division; mid-simulation perturbations in both the number of molecules of species and reaction kinetic parameters; combinatorial variation of both initial conditions and kinetic parameters to enable sensitivity analysis; use of NetCDF optimized binary format to quickly read and write large datasets; and a simple graphical user interface, written in Matlab, to help users

A stochastic analysis for a phytoplankton-zooplankton model

International Nuclear Information System (INIS)

Ge, G; Wang, H-L; Xu, J

2008-01-01

A simple phytoplankton-zooplankton nonlinear dynamical model was proposed to study the coexistence of all the species and a Hopf bifurcation was observed. In order to study the effect of environmental robustness on this system, we have stochastically perturbed the system with respect to white noise around its positive interior equilibrium. We have observed that the system remains stochastically stable around the positive equilibrium for same parametric values in the deterministic situation

Extended Plefka expansion for stochastic dynamics

International Nuclear Information System (INIS)

Bravi, B; Sollich, P; Opper, M

2016-01-01

We propose an extension of the Plefka expansion, which is well known for the dynamics of discrete spins, to stochastic differential equations with continuous degrees of freedom and exhibiting generic nonlinearities. The scenario is sufficiently general to allow application to e.g. biochemical networks involved in metabolism and regulation. The main feature of our approach is to constrain in the Plefka expansion not just first moments akin to magnetizations, but also second moments, specifically two-time correlations and responses for each degree of freedom. The end result is an effective equation of motion for each single degree of freedom, where couplings to other variables appear as a self-coupling to the past (i.e. memory term) and a coloured noise. This constitutes a new mean field approximation that should become exact in the thermodynamic limit of a large network, for suitably long-ranged couplings. For the analytically tractable case of linear dynamics we establish this exactness explicitly by appeal to spectral methods of random matrix theory, for Gaussian couplings with arbitrary degree of symmetry. (paper)

Extended Plefka expansion for stochastic dynamics

Science.gov (United States)

Bravi, B.; Sollich, P.; Opper, M.

2016-05-01

We propose an extension of the Plefka expansion, which is well known for the dynamics of discrete spins, to stochastic differential equations with continuous degrees of freedom and exhibiting generic nonlinearities. The scenario is sufficiently general to allow application to e.g. biochemical networks involved in metabolism and regulation. The main feature of our approach is to constrain in the Plefka expansion not just first moments akin to magnetizations, but also second moments, specifically two-time correlations and responses for each degree of freedom. The end result is an effective equation of motion for each single degree of freedom, where couplings to other variables appear as a self-coupling to the past (i.e. memory term) and a coloured noise. This constitutes a new mean field approximation that should become exact in the thermodynamic limit of a large network, for suitably long-ranged couplings. For the analytically tractable case of linear dynamics we establish this exactness explicitly by appeal to spectral methods of random matrix theory, for Gaussian couplings with arbitrary degree of symmetry.

Sustainable infrastructure system modeling under uncertainties and dynamics

Science.gov (United States)

Huang, Yongxi

Infrastructure systems support human activities in transportation, communication, water use, and energy supply. The dissertation research focuses on critical transportation infrastructure and renewable energy infrastructure systems. The goal of the research efforts is to improve the sustainability of the infrastructure systems, with an emphasis on economic viability, system reliability and robustness, and environmental impacts. The research efforts in critical transportation infrastructure concern the development of strategic robust resource allocation strategies in an uncertain decision-making environment, considering both uncertain service availability and accessibility. The study explores the performances of different modeling approaches (i.e., deterministic, stochastic programming, and robust optimization) to reflect various risk preferences. The models are evaluated in a case study of Singapore and results demonstrate that stochastic modeling methods in general offers more robust allocation strategies compared to deterministic approaches in achieving high coverage to critical infrastructures under risks. This general modeling framework can be applied to other emergency service applications, such as, locating medical emergency services. The development of renewable energy infrastructure system development aims to answer the following key research questions: (1) is the renewable energy an economically viable solution? (2) what are the energy distribution and infrastructure system requirements to support such energy supply systems in hedging against potential risks? (3) how does the energy system adapt the dynamics from evolving technology and societal needs in the transition into a renewable energy based society? The study of Renewable Energy System Planning with Risk Management incorporates risk management into its strategic planning of the supply chains. The physical design and operational management are integrated as a whole in seeking mitigations against the

Approaching complexity by stochastic methods: From biological systems to turbulence

Energy Technology Data Exchange (ETDEWEB)

Friedrich, Rudolf [Institute for Theoretical Physics, University of Muenster, D-48149 Muenster (Germany); Peinke, Joachim [Institute of Physics, Carl von Ossietzky University, D-26111 Oldenburg (Germany); Sahimi, Muhammad [Mork Family Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, CA 90089-1211 (United States); Reza Rahimi Tabar, M., E-mail: mohammed.r.rahimi.tabar@uni-oldenburg.de [Department of Physics, Sharif University of Technology, Tehran 11155-9161 (Iran, Islamic Republic of); Institute of Physics, Carl von Ossietzky University, D-26111 Oldenburg (Germany); Fachbereich Physik, Universitaet Osnabrueck, Barbarastrasse 7, 49076 Osnabrueck (Germany)

2011-09-15

This review addresses a central question in the field of complex systems: given a fluctuating (in time or space), sequentially measured set of experimental data, how should one analyze the data, assess their underlying trends, and discover the characteristics of the fluctuations that generate the experimental traces? In recent years, significant progress has been made in addressing this question for a class of stochastic processes that can be modeled by Langevin equations, including additive as well as multiplicative fluctuations or noise. Important results have emerged from the analysis of temporal data for such diverse fields as neuroscience, cardiology, finance, economy, surface science, turbulence, seismic time series and epileptic brain dynamics, to name but a few. Furthermore, it has been recognized that a similar approach can be applied to the data that depend on a length scale, such as velocity increments in fully developed turbulent flow, or height increments that characterize rough surfaces. A basic ingredient of the approach to the analysis of fluctuating data is the presence of a Markovian property, which can be detected in real systems above a certain time or length scale. This scale is referred to as the Markov-Einstein (ME) scale, and has turned out to be a useful characteristic of complex systems. We provide a review of the operational methods that have been developed for analyzing stochastic data in time and scale. We address in detail the following issues: (i) reconstruction of stochastic evolution equations from data in terms of the Langevin equations or the corresponding Fokker-Planck equations and (ii) intermittency, cascades, and multiscale correlation functions.

Approaching complexity by stochastic methods: From biological systems to turbulence

International Nuclear Information System (INIS)

Friedrich, Rudolf; Peinke, Joachim; Sahimi, Muhammad; Reza Rahimi Tabar, M.

2011-01-01

This review addresses a central question in the field of complex systems: given a fluctuating (in time or space), sequentially measured set of experimental data, how should one analyze the data, assess their underlying trends, and discover the characteristics of the fluctuations that generate the experimental traces? In recent years, significant progress has been made in addressing this question for a class of stochastic processes that can be modeled by Langevin equations, including additive as well as multiplicative fluctuations or noise. Important results have emerged from the analysis of temporal data for such diverse fields as neuroscience, cardiology, finance, economy, surface science, turbulence, seismic time series and epileptic brain dynamics, to name but a few. Furthermore, it has been recognized that a similar approach can be applied to the data that depend on a length scale, such as velocity increments in fully developed turbulent flow, or height increments that characterize rough surfaces. A basic ingredient of the approach to the analysis of fluctuating data is the presence of a Markovian property, which can be detected in real systems above a certain time or length scale. This scale is referred to as the Markov-Einstein (ME) scale, and has turned out to be a useful characteristic of complex systems. We provide a review of the operational methods that have been developed for analyzing stochastic data in time and scale. We address in detail the following issues: (i) reconstruction of stochastic evolution equations from data in terms of the Langevin equations or the corresponding Fokker-Planck equations and (ii) intermittency, cascades, and multiscale correlation functions.

Stochastic dynamics for reinfection by transmitted diseases

Science.gov (United States)

Barros, Alessandro S.; Pinho, Suani T. R.

2017-06-01

The use of stochastic models to study the dynamics of infectious diseases is an important tool to understand the epidemiological process. For several directly transmitted diseases, reinfection is a relevant process, which can be expressed by endogenous reactivation of the pathogen or by exogenous reinfection due to direct contact with an infected individual (with smaller reinfection rate σ β than infection rate β ). In this paper, we examine the stochastic susceptible, infected, recovered, infected (SIRI) model simulating the endogenous reactivation by a spontaneous reaction, while exogenous reinfection by a catalytic reaction. Analyzing the mean-field approximations of a site and pairs of sites, and Monte Carlo (MC) simulations for the particular case of exogenous reinfection, we obtained continuous phase transitions involving endemic, epidemic, and no transmission phases for the simple approach; the approach of pairs is better to describe the phase transition from endemic phase (susceptible, infected, susceptible (SIS)-like model) to epidemic phase (susceptible, infected, and removed or recovered (SIR)-like model) considering the comparison with MC results; the reinfection increases the peaks of outbreaks until the system reaches endemic phase. For the particular case of endogenous reactivation, the approach of pairs leads to a continuous phase transition from endemic phase (SIS-like model) to no transmission phase. Finally, there is no phase transition when both effects are taken into account. We hope the results of this study can be generalized for the susceptible, exposed, infected, and removed or recovered (SEIRIE) model, for which the state exposed (infected but not infectious), describing more realistically transmitted diseases such as tuberculosis. In future work, we also intend to investigate the effect of network topology on phase transitions when the SIRI model describes both transmitted diseases (σ social contagions (σ >1 ).

Dynamic electricity pricing for electric vehicles using stochastic programming

International Nuclear Information System (INIS)

Soares, João; Ghazvini, Mohammad Ali Fotouhi; Borges, Nuno; Vale, Zita

2017-01-01

Electric Vehicles (EVs) are an important source of uncertainty, due to their variable demand, departure time and location. In smart grids, the electricity demand can be controlled via Demand Response (DR) programs. Smart charging and vehicle-to-grid seem highly promising methods for EVs control. However, high capital costs remain a barrier to implementation. Meanwhile, incentive and price-based schemes that do not require high level of control can be implemented to influence the EVs' demand. Having effective tools to deal with the increasing level of uncertainty is increasingly important for players, such as energy aggregators. This paper formulates a stochastic model for day-ahead energy resource scheduling, integrated with the dynamic electricity pricing for EVs, to address the challenges brought by the demand and renewable sources uncertainty. The two-stage stochastic programming approach is used to obtain the optimal electricity pricing for EVs. A realistic case study projected for 2030 is presented based on Zaragoza network. The results demonstrate that it is more effective than the deterministic model and that the optimal pricing is preferable. This study indicates that adequate DR schemes like the proposed one are promising to increase the customers' satisfaction in addition to improve the profitability of the energy aggregation business. - Highlights: • A stochastic model for energy scheduling tackling several uncertainty sources. • A two-stage stochastic programming is used to tackle the developed model. • Optimal EV electricity pricing seems to improve the profits. • The propose results suggest to increase the customers' satisfaction.

Winter NH low-frequency variability in a hierarchy of low-order stochastic dynamical models of earth-atmosphere system

Science.gov (United States)

Zhao, Nan

2018-02-01

The origin of winter Northern Hemispheric low-frequency variability (hereafter, LFV) is regarded to be related to the coupled earth-atmosphere system characterized by the interaction of the jet stream with mid-latitude mountain ranges. On the other hand, observed LFV usually appears as transitions among multiple planetary-scale flow regimes of Northern Hemisphere like NAO + , AO +, AO - and NAO - . Moreover, the interaction between synoptic-scale eddies and the planetary-scale disturbance is also inevitable in the origin of LFV. These raise a question regarding how to incorporate all these aspects into just one framework to demonstrate (1) a planetary-scale dynamics of interaction of the jet stream with mid-latitude mountain ranges can really produce LFV, (2) such a dynamics can be responsible for the existence of above multiple flow regimes, and (3) the role of interaction with eddy is also clarified. For this purpose, a hierarchy of low-order stochastic dynamical models of the coupled earth-atmosphere system derived empirically from different timescale ranges of indices of Arctic Oscillation (AO), North Atlantic Oscillation (NAO), Pacific/North American (PNA), and length of day (LOD) and related probability density function (PDF) analysis are employed in this study. The results seem to suggest that the origin of LFV cannot be understood completely within the planetary-scale dynamics of the interaction of the jet stream with mid-latitude mountain ranges, because (1) the existence of multiple flow regimes such as NAO+, AO+, AO- and NAO- resulted from processes with timescales much longer than LFV itself, which may have underlying dynamics other than topography-jet stream interaction, and (2) we find LFV seems not necessarily to come directly from the planetary-scale dynamics of the interaction of the jet stream with mid-latitude mountain, although it can produce similar oscillatory behavior. The feedback/forcing of synoptic-scale eddies on the planetary

Stochastic resonance: noise-enhanced order

International Nuclear Information System (INIS)

Anishchenko, Vadim S; Neiman, Arkady B; Moss, F; Shimansky-Geier, L

1999-01-01

Stochastic resonance (SR) provides a glaring example of a noise-induced transition in a nonlinear system driven by an information signal and noise simultaneously. In the regime of SR some characteristics of the information signal (amplification factor, signal-to-noise ratio, the degrees of coherence and of order, etc.) at the output of the system are significantly improved at a certain optimal noise level. SR is realized only in nonlinear systems for which a noise-intensity-controlled characteristic time becomes available. In the present review the physical mechanism and methods of theoretical description of SR are briefly discussed. SR features determined by the structure of the information signal, noise statistics and properties of particular systems with SR are studied. A nontrivial phenomenon of stochastic synchronization defined as locking of the instantaneous phase and switching frequency of a bistable system by external periodic force is analyzed in detail. Stochastic synchronization is explored in single and coupled bistable oscillators, including ensembles. The effects of SR and stochastic synchronization of ensembles of stochastic resonators are studied both with and without coupling between the elements. SR is considered in dynamical and nondynamical (threshold) systems. The SR effect is analyzed from the viewpoint of information and entropy characteristics of the signal, which determine the degree of order or self-organization in the system. Applications of the SR concept to explaining the results of a series of biological experiments are discussed. (reviews of topical problems)

Stochastic resonance: noise-enhanced order

Energy Technology Data Exchange (ETDEWEB)

Anishchenko, Vadim S; Neiman, Arkady B [N.G. Chernyshevskii Saratov State University, Saratov (Russian Federation); Moss, F [Department of Physics and Astronomy, University of Missouri at St. Louis (United States); Shimansky-Geier, L [Humboldt University at Berlin (Germany)

1999-01-31

Stochastic resonance (SR) provides a glaring example of a noise-induced transition in a nonlinear system driven by an information signal and noise simultaneously. In the regime of SR some characteristics of the information signal (amplification factor, signal-to-noise ratio, the degrees of coherence and of order, etc.) at the output of the system are significantly improved at a certain optimal noise level. SR is realized only in nonlinear systems for which a noise-intensity-controlled characteristic time becomes available. In the present review the physical mechanism and methods of theoretical description of SR are briefly discussed. SR features determined by the structure of the information signal, noise statistics and properties of particular systems with SR are studied. A nontrivial phenomenon of stochastic synchronization defined as locking of the instantaneous phase and switching frequency of a bistable system by external periodic force is analyzed in detail. Stochastic synchronization is explored in single and coupled bistable oscillators, including ensembles. The effects of SR and stochastic synchronization of ensembles of stochastic resonators are studied both with and without coupling between the elements. SR is considered in dynamical and nondynamical (threshold) systems. The SR effect is analyzed from the viewpoint of information and entropy characteristics of the signal, which determine the degree of order or self-organization in the system. Applications of the SR concept to explaining the results of a series of biological experiments are discussed. (reviews of topical problems)

Stochastic cooling system in COSY

International Nuclear Information System (INIS)

Brittner, P.; Hacker, H.U.; Prasuhn, D.; Schug, G.; Singer, H.; Spiess, W.; Stassen, R.

1994-01-01

The stochastic cooler system in COSY is designed for proton kinetic energies between 0.8 and 2.5 GeV. Fabrication of the mechanical parts of the system is going on. Test results of the prototype measurements as well as data of the active RF-compontens are presented. (orig.)

Stochastic cooling system in COSY

Energy Technology Data Exchange (ETDEWEB)

Brittner, P [Forschungszentrum Juelich GmbH (Germany); Hacker, H U [Forschungszentrum Juelich GmbH (Germany); Prasuhn, D [Forschungszentrum Juelich GmbH (Germany); Schug, G [Forschungszentrum Juelich GmbH (Germany); Singer, H [Forschungszentrum Juelich GmbH (Germany); Spiess, W [Forschungszentrum Juelich GmbH (Germany); Stassen, R [Forschungszentrum Juelich GmbH (Germany)

1994-09-01

The stochastic cooler system in COSY is designed for proton kinetic energies between 0.8 and 2.5 GeV. Fabrication of the mechanical parts of the system is going on. Test results of the prototype measurements as well as data of the active RF-compontens are presented. (orig.)

Information Processing Capacity of Dynamical Systems

Science.gov (United States)

Dambre, Joni; Verstraeten, David; Schrauwen, Benjamin; Massar, Serge

2012-07-01

Many dynamical systems, both natural and artificial, are stimulated by time dependent external signals, somehow processing the information contained therein. We demonstrate how to quantify the different modes in which information can be processed by such systems and combine them to define the computational capacity of a dynamical system. This is bounded by the number of linearly independent state variables of the dynamical system, equaling it if the system obeys the fading memory condition. It can be interpreted as the total number of linearly independent functions of its stimuli the system can compute. Our theory combines concepts from machine learning (reservoir computing), system modeling, stochastic processes, and functional analysis. We illustrate our theory by numerical simulations for the logistic map, a recurrent neural network, and a two-dimensional reaction diffusion system, uncovering universal trade-offs between the non-linearity of the computation and the system's short-term memory.

Information Processing Capacity of Dynamical Systems

Science.gov (United States)

Dambre, Joni; Verstraeten, David; Schrauwen, Benjamin; Massar, Serge

2012-01-01

Many dynamical systems, both natural and artificial, are stimulated by time dependent external signals, somehow processing the information contained therein. We demonstrate how to quantify the different modes in which information can be processed by such systems and combine them to define the computational capacity of a dynamical system. This is bounded by the number of linearly independent state variables of the dynamical system, equaling it if the system obeys the fading memory condition. It can be interpreted as the total number of linearly independent functions of its stimuli the system can compute. Our theory combines concepts from machine learning (reservoir computing), system modeling, stochastic processes, and functional analysis. We illustrate our theory by numerical simulations for the logistic map, a recurrent neural network, and a two-dimensional reaction diffusion system, uncovering universal trade-offs between the non-linearity of the computation and the system's short-term memory. PMID:22816038

Nonequilibrium steady state in open quantum systems: Influence action, stochastic equation and power balance

International Nuclear Information System (INIS)

Hsiang, J.-T.; Hu, B.L.

2015-01-01

The existence and uniqueness of a steady state for nonequilibrium systems (NESS) is a fundamental subject and a main theme of research in statistical mechanics for decades. For Gaussian systems, such as a chain of classical harmonic oscillators connected at each end to a heat bath, and for classical anharmonic oscillators under specified conditions, definitive answers exist in the form of proven theorems. Answering this question for quantum many-body systems poses a challenge for the present. In this work we address this issue by deriving the stochastic equations for the reduced system with self-consistent backaction from the two baths, calculating the energy flow from one bath to the chain to the other bath, and exhibiting a power balance relation in the total (chain + baths) system which testifies to the existence of a NESS in this system at late times. Its insensitivity to the initial conditions of the chain corroborates to its uniqueness. The functional method we adopt here entails the use of the influence functional, the coarse-grained and stochastic effective actions, from which one can derive the stochastic equations and calculate the average values of physical variables in open quantum systems. This involves both taking the expectation values of quantum operators of the system and the distributional averages of stochastic variables stemming from the coarse-grained environment. This method though formal in appearance is compact and complete. It can also easily accommodate perturbative techniques and diagrammatic methods from field theory. Taken all together it provides a solid platform for carrying out systematic investigations into the nonequilibrium dynamics of open quantum systems and quantum thermodynamics. -- Highlights: •Nonequilibrium steady state (NESS) for interacting quantum many-body systems. •Derivation of stochastic equations for quantum oscillator chain with two heat baths. •Explicit calculation of the energy flow from one bath to the

Analysis of dynamic regimes in stochastically forced Kaldor model

International Nuclear Information System (INIS)

Bashkirtseva, Irina; Ryazanova, Tatyana; Ryashko, Lev

2015-01-01

We consider the business cycle Kaldor model forced by random noise. Detailed parametric analysis of deterministic system is carried out and zones of coexisting stable equilibrium and stable limit cycle are found. Noise-induced transitions between these attractors are studied using stochastic sensitivity function technique and confidence domains method. Critical values of noise intensity corresponding to noise-induced transitions “equilibrium → cycle” and “cycle → equilibrium” are estimated. Dominants in combined stochastic regimes are discussed.

Stochastic systems with cross-correlated Gaussian white noises

International Nuclear Information System (INIS)

Wang Cheng-Yu; Song Yu-Min; Zhou Peng; Yang Hai; Gao Yun

2010-01-01

This paper theoretically investigates three stochastic systems with cross-correlation Gaussian white noises. Both steady state properties of the stochastic nonlinear systems and the nonequilibrium transitions induced by the cross-correlated noises are studied. The stationary solutions of the Fokker—Planck equation for three specific examples are analysed. It is shown explicitly that the cross-correlation of white noises can induce nonequilibrium transitions

Who Is Afraid of Liquidity Risk? : Dynamic Portfolio Choice with Stochastic Illiquidity

NARCIS (Netherlands)

J.J.A.G. Driessen (Joost); R. Xing (Rang)

2016-01-01

textabstractRecent empirical work documents large liquidity risk premiums in stock markets. We calculate the liquidity risk premiums demanded by large investors by solving a dynamic portfolio choice problem with stochastic price impact of trading, CRRA utility and a time-varying investment

e-Dairy: a dynamic and stochastic whole-farm model that predicts biophysical and economic performance of grazing dairy systems.

Science.gov (United States)

Baudracco, J; Lopez-Villalobos, N; Holmes, C W; Comeron, E A; Macdonald, K A; Barry, T N

2013-05-01

A whole-farm, stochastic and dynamic simulation model was developed to predict biophysical and economic performance of grazing dairy systems. Several whole-farm models simulate grazing dairy systems, but most of them work at a herd level. This model, named e-Dairy, differs from the few models that work at an animal level, because it allows stochastic behaviour of the genetic merit of individual cows for several traits, namely, yields of milk, fat and protein, live weight (LW) and body condition score (BCS) within a whole-farm model. This model accounts for genetic differences between cows, is sensitive to genotype × environment interactions at an animal level and allows pasture growth, milk and supplements price to behave stochastically. The model includes an energy-based animal module that predicts intake at grazing, mammary gland functioning and body lipid change. This whole-farm model simulates a 365-day period for individual cows within a herd, with cow parameters randomly generated on the basis of the mean parameter values, defined as input and variance and co-variances from experimental data sets. The main inputs of e-Dairy are farm area, use of land, type of pasture, type of crops, monthly pasture growth rate, supplements offered, nutritional quality of feeds, herd description including herd size, age structure, calving pattern, BCS and LW at calving, probabilities of pregnancy, average genetic merit and economic values for items of income and costs. The model allows to set management policies to define: dry-off cows (ceasing of lactation), target pre- and post-grazing herbage mass and feed supplementation. The main outputs are herbage dry matter intake, annual pasture utilisation, milk yield, changes in BCS and LW, economic farm profit and return on assets. The model showed satisfactory accuracy of prediction when validated against two data sets from farmlet system experiments. Relative prediction errors were <10% for all variables, and concordance

Generalization of uncertainty relation for quantum and stochastic systems

Science.gov (United States)

Koide, T.; Kodama, T.

2018-06-01

The generalized uncertainty relation applicable to quantum and stochastic systems is derived within the stochastic variational method. This relation not only reproduces the well-known inequality in quantum mechanics but also is applicable to the Gross-Pitaevskii equation and the Navier-Stokes-Fourier equation, showing that the finite minimum uncertainty between the position and the momentum is not an inherent property of quantum mechanics but a common feature of stochastic systems. We further discuss the possible implication of the present study in discussing the application of the hydrodynamic picture to microscopic systems, like relativistic heavy-ion collisions.

Effect of spatially correlated noise on stochastic synchronization in globally coupled FitzHugh-Nagumo neuron systems

Directory of Open Access Journals (Sweden)

Yange Shao

2014-01-01

Full Text Available The phenomenon of stochastic synchronization in globally coupled FitzHugh-Nagumo (FHN neuron system subjected to spatially correlated Gaussian noise is investigated based on dynamical mean-field approximation (DMA and direct simulation (DS. Results from DMA are in good quantitative or qualitative agreement with those from DS for weak noise intensity and larger system size. Whether the consisting single FHN neuron is staying at the resting state, subthreshold oscillatory regime, or the spiking state, our investigation shows that the synchronization ratio of the globally coupled system becomes higher as the noise correlation coefficient increases, and thus we conclude that spatial correlation has an active effect on stochastic synchronization, and the neurons can achieve complete synchronization in the sense of statistics when the noise correlation coefficient tends to one. Our investigation also discloses that the noise spatial correlation plays the same beneficial role as the global coupling strength in enhancing stochastic synchronization in the ensemble. The result might be useful in understanding the information coding mechanism in neural systems.

Dynamic multi-stage dispatch of isolated wind–diesel power systems

International Nuclear Information System (INIS)

Hu, Yu; Morales, Juan M.; Pineda, Salvador; Sánchez, María Jesús; Solana, Pablo

2015-01-01

Highlights: • Optimal decision-making model for isolated hybrid wind–diesel power system is proposed. • Wind power uncertainty and conditional operating cost are considered. • Battery wear cost of the energy storage system is included in the model. • The results are compared with deterministic dispatch strategies. - Abstract: An optimal dispatch strategy is crucial for an isolated wind–diesel power system to save diesel fuel and maintain the system stability. The uncertainty associated with the stochastic character of the wind is, though, a challenging problem for this optimization. In this paper, a dynamic multi-stage decision-making model is proposed to determine the diesel power output that minimizes the cost of running and maintaining the wind–diesel power system. Optimized operational decisions for each time period are generated dynamically considering the path-dependent nature of the optimal dispatch policy, given the plausible future realizations of the wind power production. A numerical case study is analyzed and it is demonstrated that the proposed stochastic dynamic optimization model significantly outperforms the traditional deterministic dispatch strategies

Stochastic analysis of biochemical systems

CERN Document Server

Anderson, David F

2015-01-01

This book focuses on counting processes and continuous-time Markov chains motivated by examples and applications drawn from chemical networks in systems biology.  The book should serve well as a supplement for courses in probability and stochastic processes.  While the material is presented in a manner most suitable for students who have studied stochastic processes up to and including martingales in continuous time, much of the necessary background material is summarized in the Appendix. Students and Researchers with a solid understanding of calculus, differential equations, and elementary probability and who are well-motivated by the applications will find this book of interest.    David F. Anderson is Associate Professor in the Department of Mathematics at the University of Wisconsin and Thomas G. Kurtz is Emeritus Professor in the Departments of Mathematics and Statistics at that university. Their research is focused on probability and stochastic processes with applications in biology and other ar...

Maximum principle for a stochastic delayed system involving terminal state constraints.

Science.gov (United States)

Wen, Jiaqiang; Shi, Yufeng

2017-01-01

We investigate a stochastic optimal control problem where the controlled system is depicted as a stochastic differential delayed equation; however, at the terminal time, the state is constrained in a convex set. We firstly introduce an equivalent backward delayed system depicted as a time-delayed backward stochastic differential equation. Then a stochastic maximum principle is obtained by virtue of Ekeland's variational principle. Finally, applications to a state constrained stochastic delayed linear-quadratic control model and a production-consumption choice problem are studied to illustrate the main obtained result.

Dynamic asset allocation for bank under stochastic interest rates.

OpenAIRE

Chakroun, Fatma; Abid, Fathi

2014-01-01

This paper considers the optimal asset allocation strategy for bank with stochastic interest rates when there are three types of asset: Bank account, loans and securities. The asset allocation problem is to maximize the expected utility from terminal wealth of a bank's shareholders over a finite time horizon. As a consequence, we apply a dynamic programming principle to solve the Hamilton-Jacobi-Bellman (HJB) equation explicitly in the case of the CRRA utility function. A case study is given ...

A stochastic differential equation framework for the timewise dynamics of turbulent velocities

DEFF Research Database (Denmark)

Barndorff-Nielsen, Ole Eiler; Schmiegel, Jürgen

2008-01-01

We discuss a stochastic differential equation as a modeling framework for the timewise dynamics of turbulent velocities. The equation is capable of capturing basic stylized facts of the statistics of temporal velocity increments. In particular, we focus on the evolution of the probability density...

Who Is Afraid of Liquidity Risk? : Dynamic Portfolio Choice with Stochastic Illiquidity

NARCIS (Netherlands)

Driessen, Joost; Xing, R.

Recent empirical work documents large liquidity risk premiums in stock markets. We calculate the liquidity risk premiums demanded by large investors by solving a dynamic portfolio choice problem with stochastic price impact of trading, CRRA utility and a time-varying investment opportunity set. We

Portfolio Optimization with Stochastic Dividends and Stochastic Volatility

Science.gov (United States)

Varga, Katherine Yvonne

2015-01-01

We consider an optimal investment-consumption portfolio optimization model in which an investor receives stochastic dividends. As a first problem, we allow the drift of stock price to be a bounded function. Next, we consider a stochastic volatility model. In each problem, we use the dynamic programming method to derive the Hamilton-Jacobi-Bellman…

Stochastic approach to equilibrium and nonequilibrium thermodynamics.

Science.gov (United States)

Tomé, Tânia; de Oliveira, Mário J

2015-04-01

We develop the stochastic approach to thermodynamics based on stochastic dynamics, which can be discrete (master equation) and continuous (Fokker-Planck equation), and on two assumptions concerning entropy. The first is the definition of entropy itself and the second the definition of entropy production rate, which is non-negative and vanishes in thermodynamic equilibrium. Based on these assumptions, we study interacting systems with many degrees of freedom in equilibrium or out of thermodynamic equilibrium and how the macroscopic laws are derived from the stochastic dynamics. These studies include the quasiequilibrium processes; the convexity of the equilibrium surface; the monotonic time behavior of thermodynamic potentials, including entropy; the bilinear form of the entropy production rate; the Onsager coefficients and reciprocal relations; and the nonequilibrium steady states of chemical reactions.

Output-Feedback Control of Unknown Linear Discrete-Time Systems With Stochastic Measurement and Process Noise via Approximate Dynamic Programming.

Science.gov (United States)

Wang, Jun-Sheng; Yang, Guang-Hong

2017-07-25

This paper studies the optimal output-feedback control problem for unknown linear discrete-time systems with stochastic measurement and process noise. A dithered Bellman equation with the innovation covariance matrix is constructed via the expectation operator given in the form of a finite summation. On this basis, an output-feedback-based approximate dynamic programming method is developed, where the terms depending on the innovation covariance matrix are available with the aid of the innovation covariance matrix identified beforehand. Therefore, by iterating the Bellman equation, the resulting value function can converge to the optimal one in the presence of the aforementioned noise, and the nearly optimal control laws are delivered. To show the effectiveness and the advantages of the proposed approach, a simulation example and a velocity control experiment on a dc machine are employed.

Stochastic population dynamics in populations of western terrestrial garter snakes with divergent life histories.

Science.gov (United States)

Miller, David A; Clark, William R; Arnold, Stevan J; Bronikowski, Anne M

2011-08-01

Comparative evaluations of population dynamics in species with temporal and spatial variation in life-history traits are rare because they require long-term demographic time series from multiple populations. We present such an analysis using demographic data collected during the interval 1978-1996 for six populations of western terrestrial garter snakes (Thamnophis elegans) from two evolutionarily divergent ecotypes. Three replicate populations from a slow-living ecotype, found in mountain meadows of northeastern California, were characterized by individuals that develop slowly, mature late, reproduce infrequently with small reproductive effort, and live longer than individuals of three populations of a fast-living ecotype found at lakeshore locales. We constructed matrix population models for each of the populations based on 8-13 years of data per population and analyzed both deterministic dynamics based on mean annual vital rates and stochastic dynamics incorporating annual variation in vital rates. (1) Contributions of highly variable vital rates to fitness (lambda(s)) were buffered against the negative effects of stochastic variation, and this relationship was consistent with differences between the meadow (M-slow) and lakeshore (L-fast) ecotypes. (2) Annual variation in the proportion of gravid females had the greatest negative effect among all vital rates on lambda(s). The magnitude of variation in the proportion of gravid females and its effect on lambda(s) was greater in M-slow than L-fast populations. (3) Variation in the proportion of gravid females, in turn, depended on annual variation in prey availability, and its effect on lambda(s) was 4 23 times greater in M-slow than L-fast populations. In addition to differences in stochastic dynamics between ecotypes, we also found higher mean mortality rates across all age classes in the L-fast populations. Our results suggest that both deterministic and stochastic selective forces have affected the evolution of

A dynamically adaptive wavelet approach to stochastic computations based on polynomial chaos - capturing all scales of random modes on independent grids

International Nuclear Information System (INIS)

Ren Xiaoan; Wu Wenquan; Xanthis, Leonidas S.

2011-01-01

Highlights: → New approach for stochastic computations based on polynomial chaos. → Development of dynamically adaptive wavelet multiscale solver using space refinement. → Accurate capture of steep gradients and multiscale features in stochastic problems. → All scales of each random mode are captured on independent grids. → Numerical examples demonstrate the need for different space resolutions per mode. - Abstract: In stochastic computations, or uncertainty quantification methods, the spectral approach based on the polynomial chaos expansion in random space leads to a coupled system of deterministic equations for the coefficients of the expansion. The size of this system increases drastically when the number of independent random variables and/or order of polynomial chaos expansions increases. This is invariably the case for large scale simulations and/or problems involving steep gradients and other multiscale features; such features are variously reflected on each solution component or random/uncertainty mode requiring the development of adaptive methods for their accurate resolution. In this paper we propose a new approach for treating such problems based on a dynamically adaptive wavelet methodology involving space-refinement on physical space that allows all scales of each solution component to be refined independently of the rest. We exemplify this using the convection-diffusion model with random input data and present three numerical examples demonstrating the salient features of the proposed method. Thus we establish a new, elegant and flexible approach for stochastic problems with steep gradients and multiscale features based on polynomial chaos expansions.

Stochastic Averaging and Stochastic Extremum Seeking

CERN Document Server

Liu, Shu-Jun

2012-01-01

Stochastic Averaging and Stochastic Extremum Seeking develops methods of mathematical analysis inspired by the interest in reverse engineering  and analysis of bacterial  convergence by chemotaxis and to apply similar stochastic optimization techniques in other environments. The first half of the text presents significant advances in stochastic averaging theory, necessitated by the fact that existing theorems are restricted to systems with linear growth, globally exponentially stable average models, vanishing stochastic perturbations, and prevent analysis over infinite time horizon. The second half of the text introduces stochastic extremum seeking algorithms for model-free optimization of systems in real time using stochastic perturbations for estimation of their gradients. Both gradient- and Newton-based algorithms are presented, offering the user the choice between the simplicity of implementation (gradient) and the ability to achieve a known, arbitrary convergence rate (Newton). The design of algorithms...

A stochastic perturbation theory for non-autonomous systems

Energy Technology Data Exchange (ETDEWEB)

Moon, W., E-mail: wm275@damtp.cam.ac.uk [Yale University, New Haven, Connecticut 06520-8109 (United States); Wettlaufer, J. S., E-mail: wettlaufer@maths.ox.ac.uk [Yale University, New Haven, Connecticut 06520-8109 (United States); Mathematical Institute, University of Oxford, Oxford OX2 6GG (United Kingdom)

2013-12-15

We develop a perturbation theory for a class of first order nonlinear non-autonomous stochastic ordinary differential equations that arise in climate physics. The perturbative procedure produces moments in terms of integral delay equations, whose order by order decay is characterized in a Floquet-like sense. Both additive and multiplicative sources of noise are discussed and the question of how the nature of the noise influences the results is addressed theoretically and numerically. By invoking the Martingale property, we rationalize the transformation of the underlying Stratonovich form of the model to an Ito form, independent of whether the noise is additive or multiplicative. The generality of the analysis is demonstrated by developing it both for a Brownian particle moving in a periodically forced quartic potential, which acts as a simple model of stochastic resonance, as well as for our more complex climate physics model. The validity of the approach is shown by comparison with numerical solutions. The particular climate dynamics problem upon which we focus involves a low-order model for the evolution of Arctic sea ice under the influence of increasing greenhouse gas forcing ΔF{sub 0}. The deterministic model, developed by Eisenman and Wettlaufer [“Nonlinear threshold behavior during the loss of Arctic sea ice,” Proc. Natl. Acad. Sci. U.S.A. 106(1), 28–32 (2009)] exhibits several transitions as ΔF{sub 0} increases and the stochastic analysis is used to understand the manner in which noise influences these transitions and the stability of the system.

QUANTUM STOCHASTIC PROCESSES: BOSON AND FERMION BROWNIAN MOTION

Directory of Open Access Journals (Sweden)

A.E.Kobryn

2003-01-01

Full Text Available Dynamics of quantum systems which are stochastically perturbed by linear coupling to the reservoir can be studied in terms of quantum stochastic differential equations (for example, quantum stochastic Liouville equation and quantum Langevin equation. In order to work it out one needs to define the quantum Brownian motion. As far as only its boson version has been known until recently, in the present paper we present the definition which makes it possible to consider the fermion Brownian motion as well.

Stochastic simulation of time-series models combined with geostatistics to predict water-table scenarios in a Guarani Aquifer System outcrop area, Brazil

Science.gov (United States)

Manzione, Rodrigo L.; Wendland, Edson; Tanikawa, Diego H.

2012-11-01

Stochastic methods based on time-series modeling combined with geostatistics can be useful tools to describe the variability of water-table levels in time and space and to account for uncertainty. Monitoring water-level networks can give information about the dynamic of the aquifer domain in both dimensions. Time-series modeling is an elegant way to treat monitoring data without the complexity of physical mechanistic models. Time-series model predictions can be interpolated spatially, with the spatial differences in water-table dynamics determined by the spatial variation in the system properties and the temporal variation driven by the dynamics of the inputs into the system. An integration of stochastic methods is presented, based on time-series modeling and geostatistics as a framework to predict water levels for decision making in groundwater management and land-use planning. The methodology is applied in a case study in a Guarani Aquifer System (GAS) outcrop area located in the southeastern part of Brazil. Communication of results in a clear and understandable form, via simulated scenarios, is discussed as an alternative, when translating scientific knowledge into applications of stochastic hydrogeology in large aquifers with limited monitoring network coverage like the GAS.

The Matrix exponential, Dynamic Systems and Control

DEFF Research Database (Denmark)

Poulsen, Niels Kjølstad

The matrix exponential can be found in various connections in analysis and control of dynamic systems. In this short note we are going to list a few examples. The matrix exponential usably pops up in connection to the sampling process, whatever it is in a deterministic or a stochastic setting...... or it is a tool for determining a Gramian matrix. This note is intended to be used in connection to the teaching post the course in Stochastic Adaptive Control (02421) given at Informatics and Mathematical Modelling (IMM), The Technical University of Denmark. This work is a result of a study of the litterature....

A model of gene expression based on random dynamical systems reveals modularity properties of gene regulatory networks.

Science.gov (United States)

Antoneli, Fernando; Ferreira, Renata C; Briones, Marcelo R S

2016-06-01

Here we propose a new approach to modeling gene expression based on the theory of random dynamical systems (RDS) that provides a general coupling prescription between the nodes of any given regulatory network given the dynamics of each node is modeled by a RDS. The main virtues of this approach are the following: (i) it provides a natural way to obtain arbitrarily large networks by coupling together simple basic pieces, thus revealing the modularity of regulatory networks; (ii) the assumptions about the stochastic processes used in the modeling are fairly general, in the sense that the only requirement is stationarity; (iii) there is a well developed mathematical theory, which is a blend of smooth dynamical systems theory, ergodic theory and stochastic analysis that allows one to extract relevant dynamical and statistical information without solving the system; (iv) one may obtain the classical rate equations form the corresponding stochastic version by averaging the dynamic random variables (small noise limit). It is important to emphasize that unlike the deterministic case, where coupling two equations is a trivial matter, coupling two RDS is non-trivial, specially in our case, where the coupling is performed between a state variable of one gene and the switching stochastic process of another gene and, hence, it is not a priori true that the resulting coupled system will satisfy the definition of a random dynamical system. We shall provide the necessary arguments that ensure that our coupling prescription does indeed furnish a coupled regulatory network of random dynamical systems. Finally, the fact that classical rate equations are the small noise limit of our stochastic model ensures that any validation or prediction made on the basis of the classical theory is also a validation or prediction of our model. We illustrate our framework with some simple examples of single-gene system and network motifs. Copyright © 2016 Elsevier Inc. All rights reserved.

Spatial effect on stochastic dynamics of bistable evolutionary games

International Nuclear Information System (INIS)

So, Kohaku H Z; Ohtsuki, Hisashi; Kato, Takeo

2014-01-01

We consider the lifetimes of metastable states in bistable evolutionary games (coordination games), and examine how they are affected by spatial structure. A semiclassical approximation based on a path integral method is applied to stochastic evolutionary game dynamics with and without spatial structure, and the lifetimes of the metastable states are evaluated. It is shown that the population dependence of the lifetimes is qualitatively different in these two models. Our result indicates that spatial structure can accelerate the transitions between metastable states. (paper)

Long-term influence of asteroids on planet longitudes and chaotic dynamics of the solar system

Science.gov (United States)

Woillez, E.; Bouchet, F.

2017-11-01

Over timescales much longer than an orbital period, the solar system exhibits large-scale chaotic behavior and can thus be viewed as a stochastic dynamical system. The aim of the present paper is to compare different sources of stochasticity in the solar system. More precisely we studied the importance of the long term influence of asteroids on the chaotic dynamics of the solar system. We show that the effects of asteroids on planets is similar to a white noise process, when those effects are considered on a timescale much larger than the correlation time τϕ ≃ 104 yr of asteroid trajectories. We computed the timescale τe after which the effects of the stochastic evolution of the asteroids lead to a loss of information for the initial conditions of the perturbed Laplace-Lagrange secular dynamics. The order of magnitude of this timescale is precisely determined by theoretical argument, and we find that τe ≃ 104 Myr. Although comparable to the full main-sequence lifetime of the sun, this timescale is considerably longer than the Lyapunov time τI ≃ 10 Myr of the solar system without asteroids. This shows that the external sources of chaos arise as a small perturbation in the stochastic secular behavior of the solar system, rather due to intrinsic chaos.

Price-Dynamics of Shares and Bohmian Mechanics: Deterministic or Stochastic Model?

Science.gov (United States)

Choustova, Olga

2007-02-01

We apply the mathematical formalism of Bohmian mechanics to describe dynamics of shares. The main distinguishing feature of the financial Bohmian model is the possibility to take into account market psychology by describing expectations of traders by the pilot wave. We also discuss some objections (coming from conventional financial mathematics of stochastic processes) against the deterministic Bohmian model. In particular, the objection that such a model contradicts to the efficient market hypothesis which is the cornerstone of the modern market ideology. Another objection is of pure mathematical nature: it is related to the quadratic variation of price trajectories. One possibility to reply to this critique is to consider the stochastic Bohm-Vigier model, instead of the deterministic one. We do this in the present note.

Quantifying the contribution of chromatin dynamics to stochastic gene expression reveals long, locus-dependent periods between transcriptional bursts.

Science.gov (United States)

Viñuelas, José; Kaneko, Gaël; Coulon, Antoine; Vallin, Elodie; Morin, Valérie; Mejia-Pous, Camila; Kupiec, Jean-Jacques; Beslon, Guillaume; Gandrillon, Olivier

2013-02-25

A number of studies have established that stochasticity in gene expression may play an important role in many biological phenomena. This therefore calls for further investigations to identify the molecular mechanisms at stake, in order to understand and manipulate cell-to-cell variability. In this work, we explored the role played by chromatin dynamics in the regulation of stochastic gene expression in higher eukaryotic cells. For this purpose, we generated isogenic chicken-cell populations expressing a fluorescent reporter integrated in one copy per clone. Although the clones differed only in the genetic locus at which the reporter was inserted, they showed markedly different fluorescence distributions, revealing different levels of stochastic gene expression. Use of chromatin-modifying agents showed that direct manipulation of chromatin dynamics had a marked effect on the extent of stochastic gene expression. To better understand the molecular mechanism involved in these phenomena, we fitted these data to a two-state model describing the opening/closing process of the chromatin. We found that the differences between clones seemed to be due mainly to the duration of the closed state, and that the agents we used mainly seem to act on the opening probability. In this study, we report biological experiments combined with computational modeling, highlighting the importance of chromatin dynamics in stochastic gene expression. This work sheds a new light on the mechanisms of gene expression in higher eukaryotic cells, and argues in favor of relatively slow dynamics with long (hours to days) periods of quiet state.

Hopf Bifurcation of Compound Stochastic van der Pol System

Directory of Open Access Journals (Sweden)

Shaojuan Ma

2016-01-01

Full Text Available Hopf bifurcation analysis for compound stochastic van der Pol system with a bound random parameter and Gaussian white noise is investigated in this paper. By the Karhunen-Loeve (K-L expansion and the orthogonal polynomial approximation, the equivalent deterministic van der Pol system can be deduced. Based on the bifurcation theory of nonlinear deterministic system, the critical value of bifurcation parameter is obtained and the influence of random strength δ and noise intensity σ on stochastic Hopf bifurcation in compound stochastic system is discussed. At last we found that increased δ can relocate the critical value of bifurcation parameter forward while increased σ makes it backward and the influence of δ is more sensitive than σ. The results are verified by numerical simulations.

Sharp asymptotics for stochastic dynamics with parallel updating rule with self-interaction

NARCIS (Netherlands)

Bovier, A.; Nardi, F.R.; Spitoni, C.

2011-01-01

In this paper we study metastability for a stochastic dynamics with a parallel updating rule in particular for a probabilistic cellular automata. The problem is addressed in the Freidlin Wentzel regime, i.e., finite volume, small magnetic field, and in the limit when temperature tends to zero. We

Distinguishing signatures of determinism and stochasticity in spiking complex systems

Science.gov (United States)

Aragoneses, Andrés; Rubido, Nicolás; Tiana-Alsina, Jordi; Torrent, M. C.; Masoller, Cristina

2013-01-01

We describe a method to infer signatures of determinism and stochasticity in the sequence of apparently random intensity dropouts emitted by a semiconductor laser with optical feedback. The method uses ordinal time-series analysis to classify experimental data of inter-dropout-intervals (IDIs) in two categories that display statistically significant different features. Despite the apparent randomness of the dropout events, one IDI category is consistent with waiting times in a resting state until noise triggers a dropout, and the other is consistent with dropouts occurring during the return to the resting state, which have a clear deterministic component. The method we describe can be a powerful tool for inferring signatures of determinism in the dynamics of complex systems in noisy environments, at an event-level description of their dynamics.

Stochastic and Chaotic Relaxation Oscillations

NARCIS (Netherlands)

Grasman, J.; Roerdink, J.B.T.M.

1988-01-01

For relaxation oscillators stochastic and chaotic dynamics are investigated. The effect of random perturbations upon the period is computed. For an extended system with additional state variables chaotic behavior can be expected. As an example, the Van der Pol oscillator is changed into a

On quantum chaos, stochastic webs and localization in a quantum mechanical kick system

International Nuclear Information System (INIS)

Engel, U.M.

2007-01-01

In this study quantum chaos is discussed using the kicked harmonic oscillator as a model system. The kicked harmonic oscillator is characterized by an exceptional scenario of weak chaos: In the case of resonance between the frequency of the harmonic oscillator and the frequency of the periodic forcing, stochastic webs in phase space are generated by the classical dynamics. For the quantum dynamics of this system it is shown that the resulting Husimi distributions in quantum phase space exhibit the same web-like structures as the classical webs. The quantum dynamics is characterized by diffusive energy growth - just as the classical dynamics in the channels of the webs. In the case of nonresonance, the classically diffusive dynamics is found to be quantum mechanically suppressed. This bounded energy growth, which corresponds to localization in quantum phase space, is explained analytically by mapping the system onto the Anderson model. In this way, within the context of quantum chaos, the kicked harmonic oscillator is characterized by exhibiting its noteworthy geometrical and dynamical properties both classically and quantum mechanically, while at the same time there are also very distinct quantum deviations from classical properties, the most prominent example being quantum localization. (orig.)

Simulation of the stochastic wave loads using a physical modeling approach

DEFF Research Database (Denmark)

Liu, W.F.; Sichani, Mahdi Teimouri; Nielsen, Søren R.K.

2013-01-01

In analyzing stochastic dynamic systems, analysis of the system uncertainty due to randomness in the loads plays a crucial role. Typically time series of the stochastic loads are simulated using traditional random phase method. This approach combined with fast Fourier transform algorithm makes...... reliability or its uncertainty. Moreover applicability of the probability density evolution method on engineering problems faces critical difficulties when the system embeds too many random variables. Hence it is useful to devise a method which can make realization of the stochastic load processes with low...

On the theory of stochastic dynamics of magnetically confined plasma

Energy Technology Data Exchange (ETDEWEB)

El-Sharif, R.N.; El-Atoy, N.S. [Plasma and Nuclear Fusion Dept., N.R.C, Atomic Energy Authority, Cairo (Egypt)]|[Physics Dept., Girls Colleges, KSA (Saudi Arabia)

2004-07-01

This work is devoted to a study of the motion of plasma electrons in a system of two fields, a magnetic field along z-axis and wave-packet field, which propagates in the x-z plane. The strongest interaction between plasma electrons and both fields is due to their resonance with these fields. The motion of plasma electrons become stochastic when a set of resonance overlapping. Conditions for stochasticity are obtained. (orig.)

On the theory of stochastic dynamics of magnetically confined plasma

International Nuclear Information System (INIS)

El-Sharif, R.N.; El-Atoy, N.S.

2004-01-01

This work is devoted to a study of the motion of plasma electrons in a system of two fields, a magnetic field along z-axis and wave-packet field, which propagates in the x-z plane. The strongest interaction between plasma electrons and both fields is due to their resonance with these fields. The motion of plasma electrons become stochastic when a set of resonance overlapping. Conditions for stochasticity are obtained. (orig.)

Evaluation of Electric Power Procurement Strategies by Stochastic Dynamic Programming

Science.gov (United States)

Saisho, Yuichi; Hayashi, Taketo; Fujii, Yasumasa; Yamaji, Kenji

In deregulated electricity markets, the role of a distribution company is to purchase electricity from the wholesale electricity market at randomly fluctuating prices and to provide it to its customers at a given fixed price. Therefore the company has to take risk stemming from the uncertainties of electricity prices and/or demand fluctuation instead of the customers. The way to avoid the risk is to make a bilateral contact with generating companies or install its own power generation facility. This entails the necessity to develop a certain method to make an optimal strategy for electric power procurement. In such a circumstance, this research has the purpose for proposing a mathematical method based on stochastic dynamic programming and additionally considering the characteristics of the start-up cost of electric power generation facility to evaluate strategies of combination of the bilateral contract and power auto-generation with its own facility for procuring electric power in deregulated electricity market. In the beginning we proposed two approaches to solve the stochastic dynamic programming, and they are a Monte Carlo simulation method and a finite difference method to derive the solution of a partial differential equation of the total procurement cost of electric power. Finally we discussed the influences of the price uncertainty on optimal strategies of power procurement.

A stochastic approach to multi-gene expression dynamics

International Nuclear Information System (INIS)

Ochiai, T.; Nacher, J.C.; Akutsu, T.

2005-01-01

In the last years, tens of thousands gene expression profiles for cells of several organisms have been monitored. Gene expression is a complex transcriptional process where mRNA molecules are translated into proteins, which control most of the cell functions. In this process, the correlation among genes is crucial to determine the specific functions of genes. Here, we propose a novel multi-dimensional stochastic approach to deal with the gene correlation phenomena. Interestingly, our stochastic framework suggests that the study of the gene correlation requires only one theoretical assumption-Markov property-and the experimental transition probability, which characterizes the gene correlation system. Finally, a gene expression experiment is proposed for future applications of the model

Potential and flux field landscape theory. I. Global stability and dynamics of spatially dependent non-equilibrium systems.

Science.gov (United States)

Wu, Wei; Wang, Jin

2013-09-28

We established a potential and flux field landscape theory to quantify the global stability and dynamics of general spatially dependent non-equilibrium deterministic and stochastic systems. We extended our potential and flux landscape theory for spatially independent non-equilibrium stochastic systems described by Fokker-Planck equations to spatially dependent stochastic systems governed by general functional Fokker-Planck equations as well as functional Kramers-Moyal equations derived from master equations. Our general theory is applied to reaction-diffusion systems. For equilibrium spatially dependent systems with detailed balance, the potential field landscape alone, defined in terms of the steady state probability distribution functional, determines the global stability and dynamics of the system. The global stability of the system is closely related to the topography of the potential field landscape in terms of the basins of attraction and barrier heights in the field configuration state space. The effective driving force of the system is generated by the functional gradient of the potential field alone. For non-equilibrium spatially dependent systems, the curl probability flux field is indispensable in breaking detailed balance and creating non-equilibrium condition for the system. A complete characterization of the non-equilibrium dynamics of the spatially dependent system requires both the potential field and the curl probability flux field. While the non-equilibrium potential field landscape attracts the system down along the functional gradient similar to an electron moving in an electric field, the non-equilibrium flux field drives the system in a curly way similar to an electron moving in a magnetic field. In the small fluctuation limit, the intrinsic potential field as the small fluctuation limit of the potential field for spatially dependent non-equilibrium systems, which is closely related to the steady state probability distribution functional, is

Distributed Fault Detection for a Class of Nonlinear Stochastic Systems

Directory of Open Access Journals (Sweden)

Bingyong Yan

2014-01-01

Full Text Available A novel distributed fault detection strategy for a class of nonlinear stochastic systems is presented. Different from the existing design procedures for fault detection, a novel fault detection observer, which consists of a nonlinear fault detection filter and a consensus filter, is proposed to detect the nonlinear stochastic systems faults. Firstly, the outputs of the nonlinear stochastic systems act as inputs of a consensus filter. Secondly, a nonlinear fault detection filter is constructed to provide estimation of unmeasurable system states and residual signals using outputs of the consensus filter. Stability analysis of the consensus filter is rigorously investigated. Meanwhile, the design procedures of the nonlinear fault detection filter are given in terms of linear matrix inequalities (LMIs. Taking the influence of the system stochastic noises into consideration, an outstanding feature of the proposed scheme is that false alarms can be reduced dramatically. Finally, simulation results are provided to show the feasibility and effectiveness of the proposed fault detection approach.

A microscopic derivation of stochastic differential equations

International Nuclear Information System (INIS)

Arimitsu, Toshihico

1996-01-01

With the help of the formulation of Non-Equilibrium Thermo Field Dynamics, a unified canonical operator formalism is constructed for the quantum stochastic differential equations. In the course of its construction, it is found that there are at least two formulations, i.e. one is non-hermitian and the other is hermitian. Having settled which framework should be satisfied by the quantum stochastic differential equations, a microscopic derivation is performed for these stochastic differential equations by extending the projector methods. This investigation may open a new field for quantum systems in order to understand the deeper meaning of dissipation

Multi-scenario modelling of uncertainty in stochastic chemical systems

International Nuclear Information System (INIS)

Evans, R. David; Ricardez-Sandoval, Luis A.

2014-01-01

Uncertainty analysis has not been well studied at the molecular scale, despite extensive knowledge of uncertainty in macroscale systems. The ability to predict the effect of uncertainty allows for robust control of small scale systems such as nanoreactors, surface reactions, and gene toggle switches. However, it is difficult to model uncertainty in such chemical systems as they are stochastic in nature, and require a large computational cost. To address this issue, a new model of uncertainty propagation in stochastic chemical systems, based on the Chemical Master Equation, is proposed in the present study. The uncertain solution is approximated by a composite state comprised of the averaged effect of samples from the uncertain parameter distributions. This model is then used to study the effect of uncertainty on an isomerization system and a two gene regulation network called a repressilator. The results of this model show that uncertainty in stochastic systems is dependent on both the uncertain distribution, and the system under investigation. -- Highlights: •A method to model uncertainty on stochastic systems was developed. •The method is based on the Chemical Master Equation. •Uncertainty in an isomerization reaction and a gene regulation network was modelled. •Effects were significant and dependent on the uncertain input and reaction system. •The model was computationally more efficient than Kinetic Monte Carlo

Dynamics and Stochastics Festschrift Keane (Michael S) on the occasion of his 65th birthday on 2 Jan 2005

CERN Document Server

Denteneer, Dee; Verbitskiy, Evgeny

2006-01-01

Contains 29 contributions by Mike's closest colleagues and friends, covering a broad range of topics in Dynamics and Stochastics. - To celebrate Mike's scientific achievements, a conference entitled "Dynamical Systems, Probability Theory and Statistical Mechanics" was organized in Eindhoven, The Netherlands, during the week of January 3-7, 2005. This conference was hosted jointly by EURANDOM and by Philips Research. It drew over 80 participants from 5 continents, which is a sign of the warm affection and high esteem for Mike felt by the international mathematics community.

Stochastic hybrid systems with renewal transitions

NARCIS (Netherlands)

Guerreiro Tome Antunes, D.J.; Hespanha, J.P.; Silvestre, C.J.

2010-01-01

We consider Stochastic Hybrid Systems (SHSs) for which the lengths of times that the system stays in each mode are independent random variables with given distributions. We propose an analysis framework based on a set of Volterra renewal-type equations, which allows us to compute any statistical

Effect of time delay on the upper bound of the time derivative of information entropy in a stochastic dynamical system

International Nuclear Information System (INIS)

Zhang Min-Min; Mei Dong-Cheng; Wang Can-Jun

2011-01-01

The effects of the time delay on the upper bound of the time derivative of information entropy are investigated in a time-delayed dynamical system driven by correlated noise. Using the Markov approximation of the stochastic delay differential equations and the Schwartz inequality principle, we obtain an analytical expression for the upper bound U B (t) of the time derivative of the information entropy. The results show that there is a critical value of τ (delay time), and U B (t) presents opposite behaviours on difference sides of the critical value. For the case of the weak additive noise, τ can induce a reentrance transition. Delay time τ also causes a reversal behaviour in U B (t)-λ plot, where λ denotes the degree of the correlation between the two noises. (general)

The stochastic dynamics of intermittent porescale particle motion

Science.gov (United States)

Dentz, Marco; Morales, Veronica; Puyguiraud, Alexandre; Gouze, Philippe; Willmann, Matthias; Holzner, Markus

2017-04-01

Numerical and experimental data for porescale particle dynamics show intermittent patterns in Lagrangian velocities and accelerations, which manifest in long time intervals of low and short durations of high velocities [1, 2]. This phenomenon is due to the spatial persistence of particle velocities on characteristic heterogeneity length scales. In order to systematically quantify these behaviors and extract the stochastic dynamics of particle motion, we focus on the analysis of Lagrangian velocities sampled equidistantly along trajectories [3]. This method removes the intermittency observed under isochrone sampling. The space-Lagrangian velocity series can be quantified by a Markov process that is continuous in distance along streamline. It is fully parameterized in terms of the flux-weighted Eulerian velocity PDF and the characteristic pore-length. The resulting stochastic particle motion describes a continuous time random walk (CTRW). This approach allows for the process based interpretation of experimental and numerical porescale velocity, acceleration and displacement data. It provides a framework for the characterization and upscaling of particle transport and dispersion from the pore to the Darcy-scale based on the medium geometry and Eulerian flow attributes. [1] P. De Anna, T. Le Borgne, M. Dentz, A.M. Tartakovsky, D. Bolster, and P. Davy, "Flow intermittency, dispersion, and correlated continuous time random walks in porous media," Phys. Rev. Lett. 110, 184502 (2013). [2] M. Holzner, V. L. Morales, M. Willmann, and M. Dentz, "Intermittent Lagrangian velocities and accelerations in three- dimensional porous medium flow," Phys. Rev. E 92, 013015 (2015). [3] M. Dentz, P. K. Kang, A. Comolli, T. Le Borgne, and D. R. Lester, "Continuous time random walks for the evolution of Lagrangian velocities," Phys. Rev. Fluids (2016).

Economic MPC for a linear stochastic system of energy units

DEFF Research Database (Denmark)

Jørgensen, John Bagterp; Sokoler, Leo Emil; Standardi, Laura

2016-01-01

This paper summarizes comprehensively the work in four recent PhD theses from the Technical University of Denmark related to Economic MPC of future power systems. Future power systems will consist of a large number of decentralized power producers and a large number of controllable power consumers...... in addition to stochastic power producers such as wind turbines and solar power plants. Control of such large scale systems requires new control algorithms. In this paper, we formulate the control of such a system as an Economic Model Predictive Control (MPC) problem. When the power producers and controllable...... power consumers have linear dynamics, the Economic MPC may be expressed as a linear program. We provide linear models for a number of energy units in an energy system, formulate an Economic MPC for coordination of such a system. We indicate how advances in computational MPC makes the solutions...

Relevance of control theory to design and maintenance problems in time-variant reliability: The case of stochastic viability

International Nuclear Information System (INIS)

Rougé, Charles; Mathias, Jean-Denis; Deffuant, Guillaume

2014-01-01

The goal of this paper is twofold: (1) to show that time-variant reliability and a branch of control theory called stochastic viability address similar problems with different points of view, and (2) to demonstrate the relevance of concepts and methods from stochastic viability in reliability problems. On the one hand, reliability aims at evaluating the probability of failure of a system subjected to uncertainty and stochasticity. On the other hand, viability aims at maintaining a controlled dynamical system within a survival set. When the dynamical system is stochastic, this work shows that a viability problem belongs to a specific class of design and maintenance problems in time-variant reliability. Dynamic programming, which is used for solving Markovian stochastic viability problems, then yields the set of design states for which there exists a maintenance strategy which guarantees reliability with a confidence level β for a given period of time T. Besides, it leads to a straightforward computation of the date of the first outcrossing, informing on when the system is most likely to fail. We illustrate this approach with a simple example of population dynamics, including a case where load increases with time. - Highlights: • Time-variant reliability tools cannot devise complex maintenance strategies. • Stochastic viability is a control theory that computes a probability of failure. • Some design and maintenance problems are stochastic viability problems. • Used in viability, dynamic programming can find reliable maintenance actions. • Confronting reliability and control theories such as viability is promising

Integrating stochastic age-structured population dynamics into complex fisheries economic models for management evaluations: the North Sea saithe fishery as a case study

NARCIS (Netherlands)

Simons, S.L.; Bartelings, H.; Hamon, K.G.; Kempf, A.J.; Doring, R.; Temming, A.

2014-01-01

There is growing interest in bioeconomic models as tools for understanding pathways of fishery behaviour in order to assess the impact of alternative policies on natural resources. A model system is presented that combines stochastic age-structured population dynamics with complex fisheries

Distributed Adaptive Neural Control for Stochastic Nonlinear Multiagent Systems.

Science.gov (United States)

Wang, Fang; Chen, Bing; Lin, Chong; Li, Xuehua

2016-11-14

In this paper, a consensus tracking problem of nonlinear multiagent systems is investigated under a directed communication topology. All the followers are modeled by stochastic nonlinear systems in nonstrict feedback form, where nonlinearities and stochastic disturbance terms are totally unknown. Based on the structural characteristic of neural networks (in Lemma 4), a novel distributed adaptive neural control scheme is put forward. The raised control method not only effectively handles unknown nonlinearities in nonstrict feedback systems, but also copes with the interactions among agents and coupling terms. Based on the stochastic Lyapunov functional method, it is indicated that all the signals of the closed-loop system are bounded in probability and all followers' outputs are convergent to a neighborhood of the output of leader. At last, the efficiency of the control method is testified by a numerical example.

Optimizing conjunctive use of surface water and groundwater resources with stochastic dynamic programming

DEFF Research Database (Denmark)

Davidsen, Claus; Liu, Suxia; Mo, Xinguo

2014-01-01

. A stochastic dynamic programming (SDP) approach is used to minimize the basin-wide total costs arising from water allocations and water curtailments. Dynamic allocation problems with inclusion of groundwater resources proved to be more complex to solve with SDP than pure surface water allocation problems due...... to head-dependent pumping costs. These dynamic pumping costs strongly affect the total costs and can lead to non-convexity of the future cost function. The water user groups (agriculture, industry, domestic) are characterized by inelastic demands and fixed water allocation and water supply curtailment...

Mean, covariance, and effective dimension of stochastic distributed delay dynamics

Science.gov (United States)

René, Alexandre; Longtin, André

2017-11-01

Dynamical models are often required to incorporate both delays and noise. However, the inherently infinite-dimensional nature of delay equations makes formal solutions to stochastic delay differential equations (SDDEs) challenging. Here, we present an approach, similar in spirit to the analysis of functional differential equations, but based on finite-dimensional matrix operators. This results in a method for obtaining both transient and stationary solutions that is directly amenable to computation, and applicable to first order differential systems with either discrete or distributed delays. With fewer assumptions on the system's parameters than other current solution methods and no need to be near a bifurcation, we decompose the solution to a linear SDDE with arbitrary distributed delays into natural modes, in effect the eigenfunctions of the differential operator, and show that relatively few modes can suffice to approximate the probability density of solutions. Thus, we are led to conclude that noise makes these SDDEs effectively low dimensional, which opens the possibility of practical definitions of probability densities over their solution space.

A robust bi-orthogonal/dynamically-orthogonal method using the covariance pseudo-inverse with application to stochastic flow problems

Science.gov (United States)

Babaee, Hessam; Choi, Minseok; Sapsis, Themistoklis P.; Karniadakis, George Em

2017-09-01

We develop a new robust methodology for the stochastic Navier-Stokes equations based on the dynamically-orthogonal (DO) and bi-orthogonal (BO) methods [1-3]. Both approaches are variants of a generalized Karhunen-Loève (KL) expansion in which both the stochastic coefficients and the spatial basis evolve according to system dynamics, hence, capturing the low-dimensional structure of the solution. The DO and BO formulations are mathematically equivalent [3], but they exhibit computationally complimentary properties. Specifically, the BO formulation may fail due to crossing of the eigenvalues of the covariance matrix, while both BO and DO become unstable when there is a high condition number of the covariance matrix or zero eigenvalues. To this end, we combine the two methods into a robust hybrid framework and in addition we employ a pseudo-inverse technique to invert the covariance matrix. The robustness of the proposed method stems from addressing the following issues in the DO/BO formulation: (i) eigenvalue crossing: we resolve the issue of eigenvalue crossing in the BO formulation by switching to the DO near eigenvalue crossing using the equivalence theorem and switching back to BO when the distance between eigenvalues is larger than a threshold value; (ii) ill-conditioned covariance matrix: we utilize a pseudo-inverse strategy to invert the covariance matrix; (iii) adaptivity: we utilize an adaptive strategy to add/remove modes to resolve the covariance matrix up to a threshold value. In particular, we introduce a soft-threshold criterion to allow the system to adapt to the newly added/removed mode and therefore avoid repetitive and unnecessary mode addition/removal. When the total variance approaches zero, we show that the DO/BO formulation becomes equivalent to the evolution equation of the Optimally Time-Dependent modes [4]. We demonstrate the capability of the proposed methodology with several numerical examples, namely (i) stochastic Burgers equation: we

Minimal representation of matrix valued white stochastic processes and U–D factorisation of algorithms for optimal control

NARCIS (Netherlands)

Willigenburg, van L.G.; Koning, de W.L.

2013-01-01

Two different descriptions are used in the literature to formulate the optimal dynamic output feedback control problem for linear dynamical systems with white stochastic parameters and quadratic criteria, called the optimal compensation problem. One describes the matrix valued white stochastic

A decoupled approach to filter design for stochastic systems

Science.gov (United States)

Barbata, A.; Zasadzinski, M.; Ali, H. Souley; Messaoud, H.

2016-08-01

This paper presents a new theorem to guarantee the almost sure exponential stability for a class of stochastic triangular systems by studying only the stability of each diagonal subsystems. This result allows to solve the filtering problem of the stochastic systems with multiplicative noises by using the almost sure exponential stability concept. Two kinds of observers are treated: the full-order and reduced-order cases.

Stochastic resonance in bistable systems driven by harmonic noise

International Nuclear Information System (INIS)

Neiman, A.; Schimansky-Geier, L.

1994-01-01

We study stochastic resonance in a bistable system which is excited simultaneously by white and harmonic noise which we understand as the signal. In our case the spectral line of the signal has a finite width as it occurs in many real situations. Using techniques of cumulant analysis as well as computer simulations we find that the effect of stochastic resonance is preserved in the case of harmonic noise excitation. Moreover we show that the width of the spectral line of the signal at the output can be decreased via stochastic resonance. The last could be of importance in the practical using of the stochastic resonance

Quantum noise and stochastic reduction

International Nuclear Information System (INIS)

Brody, Dorje C; Hughston, Lane P

2006-01-01

In standard nonrelativistic quantum mechanics the expectation of the energy is a conserved quantity. It is possible to extend the dynamical law associated with the evolution of a quantum state consistently to include a nonlinear stochastic component, while respecting the conservation law. According to the dynamics thus obtained, referred to as the energy-based stochastic Schroedinger equation, an arbitrary initial state collapses spontaneously to one of the energy eigenstates, thus describing the phenomenon of quantum state reduction. In this paper, two such models are investigated: one that achieves state reduction in infinite time and the other in finite time. The properties of the associated energy expectation process and the energy variance process are worked out in detail. By use of a novel application of a nonlinear filtering method, closed-form solutions-algebraic in character and involving no integration-are obtained of both these models. In each case, the solution is expressed in terms of a random variable representing the terminal energy of the system and an independent noise process. With these solutions at hand it is possible to simulate explicitly the dynamics of the quantum states of complicated physical systems

Stochastic Analysis 2010

CERN Document Server

Crisan, Dan

2011-01-01

"Stochastic Analysis" aims to provide mathematical tools to describe and model high dimensional random systems. Such tools arise in the study of Stochastic Differential Equations and Stochastic Partial Differential Equations, Infinite Dimensional Stochastic Geometry, Random Media and Interacting Particle Systems, Super-processes, Stochastic Filtering, Mathematical Finance, etc. Stochastic Analysis has emerged as a core area of late 20th century Mathematics and is currently undergoing a rapid scientific development. The special volume "Stochastic Analysis 2010" provides a sa

BRS invariant stochastic quantization of Einstein gravity

International Nuclear Information System (INIS)

Nakazawa, Naohito.

1989-11-01

We study stochastic quantization of gravity in terms of a BRS invariant canonical operator formalism. By introducing artificially canonical momentum variables for the original field variables, a canonical formulation of stochastic quantization is proposed in the sense that the Fokker-Planck hamiltonian is the generator of the fictitious time translation. Then we show that there exists a nilpotent BRS symmetry in an enlarged phase space of the first-class constrained systems. The phase space is spanned by the dynamical variables, their canonical conjugate momentum variables, Faddeev-Popov ghost and anti-ghost. We apply the general BRS invariant formulation to stochastic quantization of gravity which is described as a second-class constrained system in terms of a pair of Langevin equations coupled with white noises. It is shown that the stochastic action of gravity includes explicitly the De Witt's type superspace metric which leads to a geometrical interpretation of quantum gravity analogous to nonlinear σ-models. (author)

Dynamical Behaviors of Stochastic Reaction-Diffusion Cohen-Grossberg Neural Networks with Delays

Directory of Open Access Journals (Sweden)

Li Wan

2012-01-01

Full Text Available This paper investigates dynamical behaviors of stochastic Cohen-Grossberg neural network with delays and reaction diffusion. By employing Lyapunov method, Poincaré inequality and matrix technique, some sufficient criteria on ultimate boundedness, weak attractor, and asymptotic stability are obtained. Finally, a numerical example is given to illustrate the correctness and effectiveness of our theoretical results.

Threshold Dynamics of a Stochastic SIR Model with Vertical Transmission and Vaccination

OpenAIRE

Miao, Anqi; Zhang, Jian; Zhang, Tongqian; Pradeep, B. G. Sampath Aruna

2017-01-01

A stochastic SIR model with vertical transmission and vaccination is proposed and investigated in this paper. The threshold dynamics are explored when the noise is small. The conditions for the extinction or persistence of infectious diseases are deduced. Our results show that large noise can lead to the extinction of infectious diseases which is conducive to epidemic diseases control.

A stochastic differential equation analysis of cerebrospinal fluid dynamics.

Science.gov (United States)

Raman, Kalyan

2011-01-18

Clinical measurements of intracranial pressure (ICP) over time show fluctuations around the deterministic time path predicted by a classic mathematical model in hydrocephalus research. Thus an important issue in mathematical research on hydrocephalus remains unaddressed--modeling the effect of noise on CSF dynamics. Our objective is to mathematically model the noise in the data. The classic model relating the temporal evolution of ICP in pressure-volume studies to infusions is a nonlinear differential equation based on natural physical analogies between CSF dynamics and an electrical circuit. Brownian motion was incorporated into the differential equation describing CSF dynamics to obtain a nonlinear stochastic differential equation (SDE) that accommodates the fluctuations in ICP. The SDE is explicitly solved and the dynamic probabilities of exceeding critical levels of ICP under different clinical conditions are computed. A key finding is that the probabilities display strong threshold effects with respect to noise. Above the noise threshold, the probabilities are significantly influenced by the resistance to CSF outflow and the intensity of the noise. Fluctuations in the CSF formation rate increase fluctuations in the ICP and they should be minimized to lower the patient's risk. The nonlinear SDE provides a scientific methodology for dynamic risk management of patients. The dynamic output of the SDE matches the noisy ICP data generated by the actual intracranial dynamics of patients better than the classic model used in prior research.

Stochastic Sizing of Energy Storage Systems for Wind Integration

Directory of Open Access Journals (Sweden)

D. D. Le

2018-06-01

Full Text Available In this paper, we present an optimal capacity decision model for energy storage systems (ESSs in combined operation with wind energy in power systems. We use a two-stage stochastic programming approach to take into account both wind and load uncertainties. The planning problem is formulated as an AC optimal power flow (OPF model with the objective of minimizing ESS installation cost and system operation cost. Stochastic wind and load inputs for the model are generated from historical data using clustering technique. The model is tested on the IEEE 39-bus system.

Frictions, Persistence, and Central Bank Policy in an Experimental Dynamic Stochastic General Equilibrium Economy

NARCIS (Netherlands)

Noussair, C.N.; Pfajfar, D.; Zsiros, J.

2011-01-01

New Keynesian dynamic stochastic general equilibrium models are the principal paradigm currently employed for central bank policymaking. In this paper, we construct experimental economies, populated with human subjects, with the structure of a New Keynesian DSGE model. We give individuals monetary

Control of stochastic resonance in bistable systems by using periodic signals

International Nuclear Information System (INIS)

Min, Lin; Li-Min, Fang; Yong-Jun, Zheng

2009-01-01

According to the characteristic structure of double wells in bistable systems, this paper analyses stochastic fluctuations in the single potential well and probability transitions between the two potential wells and proposes a method of controlling stochastic resonance by using a periodic signal. Results of theoretical analysis and numerical simulation show that the phenomenon of stochastic resonance happens when the time scales of the periodic signal and the noise-induced probability transitions between the two potential wells achieve stochastic synchronization. By adding a bistable system with a controllable periodic signal, fluctuations in the single potential well can be effectively controlled, thus affecting the probability transitions between the two potential wells. In this way, an effective control can be achieved which allows one to either enhance or realize stochastic resonance

Stochastic Evolution Dynamic of the Rock-Scissors-Paper Game Based on a Quasi Birth and Death Process.