arxiv: v1 [physics.data-an] 14 Dec 2015

Similar documents
1/f noise from the nonlinear transformations of the variables

2 int T. is the Fourier transform of f(t) which is the inverse Fourier transform of f. i t e

Theory of! Partial Differential Equations!

Theory of! Partial Differential Equations-I!

Elements of Stochastic Processes Lecture II Hamid R. Rabiee

Scientific Research of the Institute of Mathematics and Computer Science DIFFERENT VARIANTS OF THE BOUNDARY ELEMENT METHOD FOR PARABOLIC EQUATIONS

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL

A Limit Symmetry of Modified KdV Equation and Its Applications

Characteristics of Linear System

A Note on Fractional Electrodynamics. Abstract

1/f Fluctuations from the Microscopic Herding Model

Differential Equations

CHAPTER 2 Signals And Spectra

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Lecture #6: Continuous-Time Signals

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

Improved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method

EELE Lecture 8 Example of Fourier Series for a Triangle from the Fourier Transform. Homework password is: 14445

Backward stochastic dynamics on a filtered probability space

The motions of the celt on a horizontal plane with viscous friction

Tutorial Sheet #2 discrete vs. continuous functions, periodicity, sampling

Math 10B: Mock Mid II. April 13, 2016

Fractional Method of Characteristics for Fractional Partial Differential Equations

arxiv: v1 [math.pr] 19 Feb 2011

7 The Itô/Stratonovich dilemma

ON THE NUMBER OF FAMILIES OF BRANCHING PROCESSES WITH IMMIGRATION WITH FAMILY SIZES WITHIN RANDOM INTERVAL

Robust estimation based on the first- and third-moment restrictions of the power transformation model

GEM4 Summer School OpenCourseWare

Math 221: Mathematical Notation

On a Fractional Stochastic Landau-Ginzburg Equation

A NEW TECHNOLOGY FOR SOLVING DIFFUSION AND HEAT EQUATIONS

Plasma Astrophysics Chapter 3: Kinetic Theory. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University

arxiv: v1 [math.pr] 21 May 2010

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Random variables. A random variable X is a function that assigns a real number, X(ζ), to each outcome ζ in the sample space of a random experiment.

Problem Set on Differential Equations

Dual Representation as Stochastic Differential Games of Backward Stochastic Differential Equations and Dynamic Evaluations

Lecture 16 (Momentum and Impulse, Collisions and Conservation of Momentum) Physics Spring 2017 Douglas Fields

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

Chapter 2. First Order Scalar Equations

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems

INSTANTANEOUS VELOCITY

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Weyl sequences: Asymptotic distributions of the partition lengths

THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University

ψ(t) = V x (0)V x (t)

arxiv: v1 [cond-mat.stat-mech] 25 Aug 2016

Testing the Random Walk Model. i.i.d. ( ) r

in Engineering Prof. Dr. Michael Havbro Faber ETH Zurich, Switzerland Swiss Federal Institute of Technology

Physics Notes - Ch. 2 Motion in One Dimension

The Paradox of Twins Described in a Three-dimensional Space-time Frame

τ /2π S, ρ < ρ <... generate Dirac distribution δ ( τ ). Spectral density

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

Stochastic Structural Dynamics. Lecture-6

Solving a System of Nonlinear Functional Equations Using Revised New Iterative Method

Sensors, Signals and Noise

Vehicle Arrival Models : Headway

2.4 Cuk converter example

JOY: The Journal of Yoga Summer 2008, Volume 7, Number 2

R.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder#

Chapter 14 Wiener Processes and Itô s Lemma. Options, Futures, and Other Derivatives, 9th Edition, Copyright John C. Hull

Lecture #31, 32: The Ornstein-Uhlenbeck Process as a Model of Volatility

A quantum method to test the existence of consciousness

arxiv: v1 [math.na] 23 Feb 2016

An Introduction to Malliavin calculus and its applications

I. OBJECTIVE OF THE EXPERIMENT.

Lecture 4. Goals: Be able to determine bandwidth of digital signals. Be able to convert a signal from baseband to passband and back IV-1

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

New effective moduli of isotropic viscoelastic composites. Part I. Theoretical justification

Class Meeting # 10: Introduction to the Wave Equation

THE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES

Application of variational iteration method for solving the nonlinear generalized Ito system

Analytical Solutions of an Economic Model by the Homotopy Analysis Method

Effect of an external periodic potential on pairs of dissipative solitons

Modeling point processes by the stochastic differential equation

ENGINEERING FOR RURAL DEVELOPMENT Jelgava, MECHANISM MOTION STUDIES WITH COLLISIONS AT SEVERAL POINTS

HEAT CONDUCTION IN TECHNOLOGICAL PROCESSES

Speaker Adaptation Techniques For Continuous Speech Using Medium and Small Adaptation Data Sets. Constantinos Boulis

Cointegration in Frequency Domain*

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

Homotopy Perturbation Method for Solving Some Initial Boundary Value Problems with Non Local Conditions

Diffusion & Viscosity: Navier-Stokes Equation

Modeling by the nonlinear stochastic differential equation of the power-law distribution of extreme events in the financial systems

Dirac Particles Emission from Reissner-Nordstrom-Vaidya Black Hole

arxiv: v1 [math.ca] 15 Nov 2016

Stochastic Processes. A. Bassi. RMP 85, April-June Appendix

Stochastic Model for Cancer Cell Growth through Single Forward Mutation

Chapter 2: Principles of steady-state converter analysis

First Order RC and RL Transient Circuits

A New Perturbative Approach in Nonlinear Singularity Analysis

The Maxwell Equations, the Lorentz Field and the Electromagnetic Nanofield with Regard to the Question of Relativity

Homework 10 (Stats 620, Winter 2017) Due Tuesday April 18, in class Questions are derived from problems in Stochastic Processes by S. Ross.

EMS SCM joint meeting. On stochastic partial differential equations of parabolic type

EE 435. Lecture 31. Absolute and Relative Accuracy DAC Design. The String DAC

Mean-square Stability Control for Networked Systems with Stochastic Time Delay

Stochastic Loewner Evolution

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

Transcription:

1/ noise rom he nonlinear ransormaions o he variables Bronislovas Kaulakys, Miglius Alaburda, and Julius Ruseckas Insiue o Theoreical Physics and Asronomy, Vilnius Universiy, A. Gošauo 1, 118 Vilnius, Lihuania arxiv:151.98v1 [physics.daa-an] 1 Dec 15 The origin o he low-requency noise wih power specrum 1/ β also known as 1/ lucuaions or licker noise) remains a challenge. Recenly, he nonlinear sochasic dierenial equaions or modeling 1/ β noise have been proposed and analyzed. Here we use he sel-similariy properies o his model wih respec o he nonlinear ransormaions o he variable o hese equaions and show ha 1/ β noise o he observable may yield rom he power-law ransormaions o well-known sandard processes, like he Brownian moion, Bessel and similar sochasic processes. Analyical and numerical invesigaions o such echniques or modeling processes wih 1/ β lucuaions is presened. I. INTRODUCTION Dieren heories and models have been proposed or eplanaion o he ubiquious 1/ β noise phenomena, observable or abou eighy years in dieren sysems rom physics o inancial markes see, e.g., Res. [1 8] and reerences herein). Recenly, he sochasic model o 1/ β noise, based on he nonlinear sochasic dierenial equaions d = η λ ) η 1 d+ η dw, 1) where is he signal wih 1/ β specrum, η 1 is he nonlineariy eponen, λ is he eponen o he seady-sae disribuion P ss λ, and W is a Wiener process Brownian moion), has been proposed and analyzed.[6, 7] The relaion o he eponen β in he specrum o he parameers o Eq. 1) is given by β = 1+ λ 3 η 1). ) Eq. 1) may be derived rom he poin process model [6], rom scaling properies o he signal [7] or rom he agenbased herding model [9]. Here we employ he sel-similariy propery o Eq. 1) wih respec he power-law ransormaions o he variable. We show ha processes wih 1/ β specrum may yield rom he nonlinear ransormaions o he variable o he widespread processes, e.g., rom he Brownian moion, Bessel or similar amiliar processes. II. TRANSFORMATIONS From he scaling o he power specral densiy PSD) S) β, Sa) a β S), 3) according o he Wiener-Khinchine heorem yields he scaling o he auocorrelaion uncion Ca) a β 1 C) ) in someime inerval 1/ ma 1/ min. Recenly [7] i wasshown ha rom he scaling) i ollows he nonlinear sochasic dierenial equaion SDE) 1). Due o his scaling he nonlinear SDE d = η λ ) η 1 d+ η dw, 5) bronislovas.kaulakys@ai.vu.l miglius.alaburda@ai.vu.l julius.ruseckas@ai.vu.l

a) Py) 1-1 -7 1-1 1 3 y b) S) 1 8 1 6 1 1 1-1 -8 1-7 1 - FIG. 1. Seady-saedisribuion, a), andapproimaely 1/ PSD,b), ohevariabley generaedbyhesimple equaion1), red circles, as well as o y = 1/ generaed by Eq. 13), blue riangles, wih he appropriae resricions[7] o he diusion inervals. generaing signal wih PSD S) 1 β, β = 1+ λ 3 η 1) 6) aer he nonlinear ransormaion = 1 y δ, 7) wih δ being he ransormaion eponen, yields SDE or he variable y o he same orm, dy = η y λ ) y ηy 1 d+y ηy dw 8) wih Eq. 8) generaes signal y wih he PSD η y = 1 δη 1), λ y = 1 δλ 1). 9) S) 1 βy, β y = 1+ λ y 3 η y 1) = β + 1+δ δη 1). 1) Thereore, 1/ β noise o he observable may yield rom he nonlinear dependence 7) o his observable on anoher variable y resuled rom simple or more common equaion 8), where β = β y + 1+δ η y 1. 11) III. EXAMPLES The simples relaion beween he variable y and he observable is he inverse ransormaion = 1/y. From he inerrelaions 1) and 11) beween he parameers β and β y we obain ha simple equaion wihou he dri erm dy = 1 y dw 1) resuls in 1/ noise o he observable = 1/y, insead o very nonlinear equaion [1] d = d+ 5/ dw 13)

3 a) P) 1 1-1 -8 1-1 1 3 b) S) 1 3 1 1-1 -8 1-7 1 - FIG.. Seady-sae disribuion, a), and 1/ PSD, b), o he observable = 1/y, where he variable y is generaed by Eq. 1), red circles, as well as o generaed by Eq. 13), blue riangles. a) y) 8 7 6 5 3 1 6 8 1 b) ) 1 8 6 6 8 1 c) ) 1 8 6 6 8 1 FIG. 3. Eamples o he signals o 1D Brownian moion o he variable y, a), o he inverse 1D Brownian moion o he observable = 1/y yielding in 1/ noise, b), and o = 1/ y observable resuling in pure 1/ noise, c). or he observable. I should be noed, ha Eq. 1) coincides wih equaion or he inereven ime obained orm he Brownian moion dτ = 1 τ dw, 1) dτ k = dw k 15) o he inereven ime τ k in he evens space or k-space.[1] Figures1anddemonsraeheappearanceo1/ noiseoheobservable, asaresulohe inverseransormaion = 1/y o he variable y generaed by he simple equaion wihou he dri erm 1) wih approimae 1/ PSD o he variable y. Noe, ha special cases o Eq. 8) are: i) an equaion wih he addiive noise, η y =, and nonlinear dri, i.e., he Bessel process o order N = 1 λ y, ii) he order N = 1 λ y ) squared Bessel process when η y = 1/ and λ y = 1 N/, iii) he special eponenial resricions o he variable y in Eq. 8) yield Consan Elasiciy o Variance CEV) process or iv) Co-Ingersoll-Ross CIR) process when η = 1/ [9, 11]. Here we will demonsrae he possibiliy o obain 1/ β noise rom he Bessel process and rom he Brownian moion. Transormaion 7) wih δ = 1/η 1), i.e., y = 1 η, yields he special orm o Eq. 8), i.e., he Bessel equaion dy = N 1 d y +dw 16) corresponding o N-dimensional Wiener process, or Bessel process wih inde ν = N/ 1. Here N = λ 1 η 1 = 1 λ y. 17)

a) P) 1 1-1 -8 1-1 1 3 b) S) 1 6 1 1 1-1 - 1 FIG.. Seady-sae disribuion, a), and PSD, b), o 1D Brownian moion o he variable y, yielding 1/ noise, red circles, he inverse 1D Brownian moion o = 1/y yielding in 1/ noise, blue riangles, and o he observable = 1/ y resuling in pure 1/ noise, green squares. a) P) 1 1-1 - b) S) 1-1 - 1 FIG. 5. Seady-sae disribuion, a), and PSD, b), o he observable as inverse o he disance rom he beginning, = 1/r, o he Brownian moion in D. On he oher hand λ = 1+ N δ, β = 1+ N δ. 18) Thus, he Bessel process or N-dimensional Brownian moion can cause 1/ β lucuaions o he observable as a uncion 7) o his process. So, he pure 1/ noise yields when N = δ, e.g., rom 1D Brownian moion o y or he observable = 1/ y and D Brownian moion o y or he observable = 1/y. Figure 3 demonsraes eamples o he signals o 1D Brownian moion o he variable y, he inverse 1D Brownian moion o = 1/y yielding in 1/ noise and observable = 1/ y resuling in pure 1/ noise. In igure he disribuion densiy and PSD o he corresponding observables are presened. In igure 5 he pure 1/ noise o he observable as inverse o he disance rom he beginning, = 1/r, o he Brownian moion in D, r = B1 )+B ), 19) i.e., o he Bessel process wih he inde ν = is shown. Here B 1 ) and B ) are independen Brownian moions. We see he pure 1/ noise.

5 IV. CONCLUSIONS Thus, he proper uncion o he widespread process, e.g., he inverse ransormaion o he order Bessel process and he inverse square roo o he Brownian moion yield he pure 1/ noise, observable, e.g., in condense maer. Possible relevancies o hese ransormaions o 1/ noise in condensed maer: i) 1/ lucuaions o he volage, U = I/G, a consan curren I as a resul o Brownian lucuaions o he conduciviy G, ii) 1/ lucuaions o resisiviy, ρ = 1/σ, as Brownian lucuaions o he conduciviy σ. [1] M. B. Weissman, Rev. Mod. Phys. 6 1988) 537. [] T. Gisiger, Biol. Rev. 76 1) 161. [3] M. Li and W. Zhao, Mah. Probl. Engin. 1 1) 67368. [] A. A. Balandin,Naure Nanoechnol. 8 13) 59. [5] E. Paladino, Y. M. Galperin, G. Falci, and B. L. Alshuler, Rev. Mod. Phys. 86 1) 361. [6] B. Kaulakys and M. Alaburda, J. Sa. Mech. 9 9) P51. [7] J. Ruseckas and B. Kaulakys, J. Sa. Mech. 1 1) P65. [8] M. A. Rodriguez, Phys. Rev. E 9 1) 1. [9] J. Ruseckas, B. Kaulakys and V. Gonis, EPL 96 11) 67. [1] B. Kaulakys and J. Ruseckas, Phys. Rev. E 7 ) 11. [11] M. Jeanblanc, M. Yor and M. Chesney Mahemaical Mehods or Financial Markes Springer, London, 9).