Stochastic Process Models for Packet/Analytic-Based Network Simulations

Total Page:16

File Type:pdf, Size:1020Kb

Stochastic Process Models for Packet/Analytic-Based Network Simulations ACM/IEEE PADS 2008 1 Stochastic Process Models for Packet/Analytic-Based Network Simulations Robert G. Cole, George Riley, Derya Cansever and William Yurcik Index Terms|Hybrid Simulation Models, Event Driven Simulations, Network Models, Stochastic Queues and Dif- summed internal load fusion Equations. external arrival external departure Abstract network queue E present our preliminary work that develops a new W approach to hybrid packet/analytic network sim- network ulations for improved network simulation fidelity, scale, Fig. 1. An example network of queues with internal and external and simulation efficiency. Much work in the literature ad- arrival and departure traffic. dresses this topic, including [10] [11] [8] [12] [13] and others. Current approaches rely upon models, which we refer to in this paper as Deterministic Fluid Models [9] [12], to address foreground traffic, which is of primary interest, and back- the analytic modeling aspects of these hybrid simulations. ground traffic, which exists solely to provide competition Instead we draw upon an extensive literature on stochas- to the foreground traffic being measured. In the hybrid tic models of queues and (eventually) networks of queues approach discussed here, analytic models are used is to es- to implement a hybrid stochastic model/packet network timate the impact of the background traffic on the fore- simulation. We will refer to our approach as Stochastic ground traffic, which is handled explicitly via discrete- Fluid Models throughout this paper. We outline our ap- event handling with full packet{level detail. In the re- proach, present test cases, and present simulation results mainder of this paper, we will assume that the analytically comparing the measured queue metrics from our approach modeled traffic represents, in some sense, background traf- for hybrid simulation to those of a deterministic fluid model fic and that the explicitly handled packet traffic is fore- hybrid simulation and a full packet{level simulation. We ground traffic. We will use the terms background traffic also discuss plans for future areas of research on this ap- and foreground traffic to distinguish between the analyt- proach. ically modeled versus the explicitly handled event packet traffic. Of course other ways to divide up the analytically I. Introduction modeled and packet handled traffic are possible. Consider a network model as found in Figure 1. The There are two interesting aspects of the analytic model- network model is comprised of nodes and communica- ing in the context of a network of queues. One relates to tions links. Associated with each communications link is a the allocation and estimation of the intermediate load gen- queue. Traffic can arrive and depart the network model at erated by the traffic at the nodes within the network based each node in the network. Each node also carries internal upon the assumed routing patterns, total estimated exter- traffic which is forwarded throughout the network between nal traffic loads, finite queue sizes and associated network the traffic source and destination nodes. internal packet losses. The other aspect relates to meth- In hybrid simulation models, some of the traffic is han- ods of mixing, at a given queue, the analytically modeled dled via analytic methods and some of the traffic is han- background traffic with the explicitly handled foreground dled via more CPU and memory intensive discrete-event, traffic. Our focus in this paper is the later; hence we con- packet-level handling. The more traffic modelled analyt- centrate on methods to mix the analytically modeled traffic ically, the more efficient the simulation becomes (in the with the event-driven, packet-level traffic at a given queue. general case), and thus the simulation can scale to larger A future paper will concentrate on the network equations. networks. Typically, there is a distinction drawn between The remainder of this paper is organized as follows: In the next section we discuss previous methods for hybrid R. G. Cole is with the Applied Physics Lab and the Department of network simulation which rely upon deterministic fluid Computer Science, Johns Hopkins University, Baltimore, MD. Phone: +1 443 778{6951, e-mail: [email protected] model approximations to queue dynamics. In Section III G. Riley is with the School of Electrical and Computer En- we present our approach based upon methods in stochas- gineering, Georgia Institute of Technology, Atlanta, GA. Phone: tic queue dynamics. In Section IV we develop our hy- +1 404 463{1774, e-mail: [email protected] D. Cansever is the Program Director of Advanced Internet- brid equations for a specific instance of a stochastic queue working, SI International, Inc. Phone: +1 703 234{6960, e-mail: model, i.e., a Brownian Motion model based upon solutions [email protected] to the Fokker-Planck Equation. In Section V we report on W. Yurcik is with the Department of Computer Science, Univer- sity of Texas-Dallas, Dallas, TX. Phone: +1 309 531{1570, e-mail: our initial simulations investigations comparing three test [email protected] cases; one based upon pure event-driven packet-level sim- 2 ACM/IEEE PADS 2008 ulations, one based upon hybrid deterministic fluid/packet ωmax simulation and one based upon our hybrid stochastic pro- ω+ (t i ) slope = cess/packet simulation. In Section VI we discuss conclu- λd − µ sions and future investigations. work ω− (t i ) ωi ( work carried by packet i ) II. Deterministic Models time Most current approaches to performance speedup for i i+1 i+2 i+3 network simulations involving hybrid event-analytic sim- packet arrivals ulation rely on analytic models based upon deterministic Fluid Flow Approximations (FFAs), e.g., [7] [4] [10] [11] [8] Fig. 2. An illustration of the hybrid deterministic modeling approach [13]. We refer to these approaches as Deterministic Mod- for mixing traffic at a queue. els because the evolution of the queue dynamics modeling the background analytic traffic is assumed to be a deter- ministic process captured in the form of differential equa- rate at the queue, i.e., the inverse of the link bandwidth. tions. The fluid dynamics is derived from the integration Let w(t) be the total workload in the queue at time t, and of these equations. In some cases the differential equations w−(ti) and w+(ti) be the total workload in the queue just are solvable, e.g., fixed arrival rate models, and numeri- prior to and just following (respectively) a distinguished cal integration of the differential equations is not neces- time event, e.g., the arrival of packet i at time ti. Let wi sary. Both cases address variable mean arrival rates; one be the work carried by the ith packet to the queue. Finally, through time dependent variables in the differential equa- let wmax be the maximum amount of work (or backlog) in tions [9] [13] and one through discrete event fluid models the queue due to its finite size. where rate changes are propagated throughout the network Figure 2 gives a pictorial representation of the method via events [4] [11] [7]. used to mix analytically modeled background traffic with Two approaches to using deterministic analytic models explicitly handled packet traffic. Individual packet arrival for hybrid simulations are found in the literature. One events are indicated along the bottom axis of the figure. approach, used in [13], divides the network into a hierarchy These packets carry with them a given workload which is a consisting of a high speed core and a lower speed edge. function of their packet size and the link bandwidth serv- In the core, all traffic is modeled analytically and in the ing the queue in question. If the packet is allowed to enter edge all traffic is modeled via packet-level, event-driven the queue, then the workload in queue jumps from w−(ti) simulation. Packet traffic traversing the core is converted to w+(ti) where the difference is the workload, wi, carried to fluid load on the core and is then converted back to by the packet. In between packet arrivals, the workload packet-level traffic at the far edge. evolves deterministically. As mentioned, we assume the Another approach, used in [4] [11] and [8], is to treat queue dynamics evolve at a constant rate between packet some of the traffic throughout the network analytically. arrivals (as long as the workload does not exceed the max- Then develop methods to explicitly mix both analytic and imum wmax or drops below the minimum of zero) given by packet traffic at each multiplexing point throughout the the difference λd −µ which is fixed over the averaging time network. Our interests are in this later approach. interval δ. Then the explicit handling of the packet traffic and the A. Deterministic/Packet Mixing Equations process for updating the impact of the background traffic Here, we develop the simplest, least assuming method are as follows: to mix deterministic traffic with packet traffic. There exist • If, during the current δ time averaging period, λd < µ, methods to improve the fidelity of our example determinis- i.e., the background traffic does not overflow the server tic model, e.g., [4] [11] [8], but these improvements require rate. For each packet arrival bringing wi work to the a priori knowledge of the behavior of the foreground traf- queue, do the following tasks: fic, which in general is not known. This choice is made { Increment the packet arrival counter i at the queue. here for several reasons. First, we wish to make an equiva- { Update the queue backlog just prior to the packet lent comparison to our stochastic models presented below arrival time ti, as in order to access its ability to negate a reliance on a priori w (t ) max[w (t )+(λ −µ)(t −t ); 0] (1) traffic knowledge.
Recommended publications
  • Poisson Processes Stochastic Processes
    Poisson Processes Stochastic Processes UC3M Feb. 2012 Exponential random variables A random variable T has exponential distribution with rate λ > 0 if its probability density function can been written as −λt f (t) = λe 1(0;+1)(t) We summarize the above by T ∼ exp(λ): The cumulative distribution function of a exponential random variable is −λt F (t) = P(T ≤ t) = 1 − e 1(0;+1)(t) And the tail, expectation and variance are P(T > t) = e−λt ; E[T ] = λ−1; and Var(T ) = E[T ] = λ−2 The exponential random variable has the lack of memory property P(T > t + sjT > t) = P(T > s) Exponencial races In what follows, T1;:::; Tn are independent r.v., with Ti ∼ exp(λi ). P1: min(T1;:::; Tn) ∼ exp(λ1 + ··· + λn) . P2 λ1 P(T1 < T2) = λ1 + λ2 P3: λi P(Ti = min(T1;:::; Tn)) = λ1 + ··· + λn P4: If λi = λ and Sn = T1 + ··· + Tn ∼ Γ(n; λ). That is, Sn has probability density function (λs)n−1 f (s) = λe−λs 1 (s) Sn (n − 1)! (0;+1) The Poisson Process as a renewal process Let T1; T2;::: be a sequence of i.i.d. nonnegative r.v. (interarrival times). Define the arrival times Sn = T1 + ··· + Tn if n ≥ 1 and S0 = 0: The process N(t) = maxfn : Sn ≤ tg; is called Renewal Process. If the common distribution of the times is the exponential distribution with rate λ then process is called Poisson Process of with rate λ. Lemma. N(t) ∼ Poisson(λt) and N(t + s) − N(s); t ≥ 0; is a Poisson process independent of N(s); t ≥ 0 The Poisson Process as a L´evy Process A stochastic process fX (t); t ≥ 0g is a L´evyProcess if it verifies the following properties: 1.
    [Show full text]
  • A Stochastic Processes and Martingales
    A Stochastic Processes and Martingales A.1 Stochastic Processes Let I be either IINorIR+.Astochastic process on I with state space E is a family of E-valued random variables X = {Xt : t ∈ I}. We only consider examples where E is a Polish space. Suppose for the moment that I =IR+. A stochastic process is called cadlag if its paths t → Xt are right-continuous (a.s.) and its left limits exist at all points. In this book we assume that every stochastic process is cadlag. We say a process is continuous if its paths are continuous. The above conditions are meant to hold with probability 1 and not to hold pathwise. A.2 Filtration and Stopping Times The information available at time t is expressed by a σ-subalgebra Ft ⊂F.An {F ∈ } increasing family of σ-algebras t : t I is called a filtration.IfI =IR+, F F F we call a filtration right-continuous if t+ := s>t s = t. If not stated otherwise, we assume that all filtrations in this book are right-continuous. In many books it is also assumed that the filtration is complete, i.e., F0 contains all IIP-null sets. We do not assume this here because we want to be able to change the measure in Chapter 4. Because the changed measure and IIP will be singular, it would not be possible to extend the new measure to the whole σ-algebra F. A stochastic process X is called Ft-adapted if Xt is Ft-measurable for all t. If it is clear which filtration is used, we just call the process adapted.The {F X } natural filtration t is the smallest right-continuous filtration such that X is adapted.
    [Show full text]
  • Introduction to Stochastic Processes - Lecture Notes (With 33 Illustrations)
    Introduction to Stochastic Processes - Lecture Notes (with 33 illustrations) Gordan Žitković Department of Mathematics The University of Texas at Austin Contents 1 Probability review 4 1.1 Random variables . 4 1.2 Countable sets . 5 1.3 Discrete random variables . 5 1.4 Expectation . 7 1.5 Events and probability . 8 1.6 Dependence and independence . 9 1.7 Conditional probability . 10 1.8 Examples . 12 2 Mathematica in 15 min 15 2.1 Basic Syntax . 15 2.2 Numerical Approximation . 16 2.3 Expression Manipulation . 16 2.4 Lists and Functions . 17 2.5 Linear Algebra . 19 2.6 Predefined Constants . 20 2.7 Calculus . 20 2.8 Solving Equations . 22 2.9 Graphics . 22 2.10 Probability Distributions and Simulation . 23 2.11 Help Commands . 24 2.12 Common Mistakes . 25 3 Stochastic Processes 26 3.1 The canonical probability space . 27 3.2 Constructing the Random Walk . 28 3.3 Simulation . 29 3.3.1 Random number generation . 29 3.3.2 Simulation of Random Variables . 30 3.4 Monte Carlo Integration . 33 4 The Simple Random Walk 35 4.1 Construction . 35 4.2 The maximum . 36 1 CONTENTS 5 Generating functions 40 5.1 Definition and first properties . 40 5.2 Convolution and moments . 42 5.3 Random sums and Wald’s identity . 44 6 Random walks - advanced methods 48 6.1 Stopping times . 48 6.2 Wald’s identity II . 50 6.3 The distribution of the first hitting time T1 .......................... 52 6.3.1 A recursive formula . 52 6.3.2 Generating-function approach .
    [Show full text]
  • Non-Local Branching Superprocesses and Some Related Models
    Published in: Acta Applicandae Mathematicae 74 (2002), 93–112. Non-local Branching Superprocesses and Some Related Models Donald A. Dawson1 School of Mathematics and Statistics, Carleton University, 1125 Colonel By Drive, Ottawa, Canada K1S 5B6 E-mail: [email protected] Luis G. Gorostiza2 Departamento de Matem´aticas, Centro de Investigaci´ony de Estudios Avanzados, A.P. 14-740, 07000 M´exicoD. F., M´exico E-mail: [email protected] Zenghu Li3 Department of Mathematics, Beijing Normal University, Beijing 100875, P.R. China E-mail: [email protected] Abstract A new formulation of non-local branching superprocesses is given from which we derive as special cases the rebirth, the multitype, the mass- structured, the multilevel and the age-reproduction-structured superpro- cesses and the superprocess-controlled immigration process. This unified treatment simplifies considerably the proof of existence of the old classes of superprocesses and also gives rise to some new ones. AMS Subject Classifications: 60G57, 60J80 Key words and phrases: superprocess, non-local branching, rebirth, mul- titype, mass-structured, multilevel, age-reproduction-structured, superprocess- controlled immigration. 1Supported by an NSERC Research Grant and a Max Planck Award. 2Supported by the CONACYT (Mexico, Grant No. 37130-E). 3Supported by the NNSF (China, Grant No. 10131040). 1 1 Introduction Measure-valued branching processes or superprocesses constitute a rich class of infinite dimensional processes currently under rapid development. Such processes arose in appli- cations as high density limits of branching particle systems; see e.g. Dawson (1992, 1993), Dynkin (1993, 1994), Watanabe (1968). The development of this subject has been stimu- lated from different subjects including branching processes, interacting particle systems, stochastic partial differential equations and non-linear partial differential equations.
    [Show full text]
  • Simulation of Markov Chains
    Copyright c 2007 by Karl Sigman 1 Simulating Markov chains Many stochastic processes used for the modeling of financial assets and other systems in engi- neering are Markovian, and this makes it relatively easy to simulate from them. Here we present a brief introduction to the simulation of Markov chains. Our emphasis is on discrete-state chains both in discrete and continuous time, but some examples with a general state space will be discussed too. 1.1 Definition of a Markov chain We shall assume that the state space S of our Markov chain is S = ZZ= f:::; −2; −1; 0; 1; 2;:::g, the integers, or a proper subset of the integers. Typical examples are S = IN = f0; 1; 2 :::g, the non-negative integers, or S = f0; 1; 2 : : : ; ag, or S = {−b; : : : ; 0; 1; 2 : : : ; ag for some integers a; b > 0, in which case the state space is finite. Definition 1.1 A stochastic process fXn : n ≥ 0g is called a Markov chain if for all times n ≥ 0 and all states i0; : : : ; i; j 2 S, P (Xn+1 = jjXn = i; Xn−1 = in−1;:::;X0 = i0) = P (Xn+1 = jjXn = i) (1) = Pij: Pij denotes the probability that the chain, whenever in state i, moves next (one unit of time later) into state j, and is referred to as a one-step transition probability. The square matrix P = (Pij); i; j 2 S; is called the one-step transition matrix, and since when leaving state i the chain must move to one of the states j 2 S, each row sums to one (e.g., forms a probability distribution): For each i X Pij = 1: j2S We are assuming that the transition probabilities do not depend on the time n, and so, in particular, using n = 0 in (1) yields Pij = P (X1 = jjX0 = i): (Formally we are considering only time homogenous MC's meaning that their transition prob- abilities are time-homogenous (time stationary).) The defining property (1) can be described in words as the future is independent of the past given the present state.
    [Show full text]
  • A Representation for Functionals of Superprocesses Via Particle Picture Raisa E
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Elsevier - Publisher Connector stochastic processes and their applications ELSEVIER Stochastic Processes and their Applications 64 (1996) 173-186 A representation for functionals of superprocesses via particle picture Raisa E. Feldman *,‘, Srikanth K. Iyer ‘,* Department of Statistics and Applied Probability, University oj’ Cull~ornia at Santa Barbara, CA 93/06-31/O, USA, and Department of Industrial Enyineeriny and Management. Technion - Israel Institute of’ Technology, Israel Received October 1995; revised May 1996 Abstract A superprocess is a measure valued process arising as the limiting density of an infinite col- lection of particles undergoing branching and diffusion. It can also be defined as a measure valued Markov process with a specified semigroup. Using the latter definition and explicit mo- ment calculations, Dynkin (1988) built multiple integrals for the superprocess. We show that the multiple integrals of the superprocess defined by Dynkin arise as weak limits of linear additive functionals built on the particle system. Keywords: Superprocesses; Additive functionals; Particle system; Multiple integrals; Intersection local time AMS c.lassijication: primary 60555; 60H05; 60580; secondary 60F05; 6OG57; 60Fl7 1. Introduction We study a class of functionals of a superprocess that can be identified as the multiple Wiener-It6 integrals of this measure valued process. More precisely, our objective is to study these via the underlying branching particle system, thus providing means of thinking about the functionals in terms of more simple, intuitive and visualizable objects. This way of looking at multiple integrals of a stochastic process has its roots in the previous work of Adler and Epstein (=Feldman) (1987), Epstein (1989), Adler ( 1989) Feldman and Rachev (1993), Feldman and Krishnakumar ( 1994).
    [Show full text]
  • Immigration Structures Associated with Dawson-Watanabe Superprocesses
    View metadata, citation and similar papers at core.ac.uk brought to you COREby provided by Elsevier - Publisher Connector stochastic processes and their applications ELSEVIER StochasticProcesses and their Applications 62 (1996) 73 86 Immigration structures associated with Dawson-Watanabe superprocesses Zeng-Hu Li Department of Mathematics, Beijing Normal University, Be!ring 100875, Peoples Republic of China Received March 1995: revised November 1995 Abstract The immigration structure associated with a measure-valued branching process may be described by a skew convolution semigroup. For the special type of measure-valued branching process, the Dawson Watanabe superprocess, we show that a skew convolution semigroup corresponds uniquely to an infinitely divisible probability measure on the space of entrance laws for the underlying process. An immigration process associated with a Borel right superpro- cess does not always have a right continuous realization, but it can always be obtained by transformation from a Borel right one in an enlarged state space. Keywords: Dawson-Watanabe superprocess; Skew convolution semigroup; Entrance law; Immigration process; Borel right process AMS classifications: Primary 60J80; Secondary 60G57, 60J50, 60J40 I. Introduction Let E be a Lusin topological space, i.e., a homeomorphism of a Borel subset of a compact metric space, with the Borel a-algebra .N(E). Let B(E) denote the set of all bounded ~(E)-measurable functions on E and B(E) + the subspace of B(E) compris- ing non-negative elements. Denote by M(E) the space of finite measures on (E, ~(E)) equipped with the topology of weak convergence. Forf e B(E) and # e M(E), write ~l(f) for ~Efd/~.
    [Show full text]
  • LECTURES 3 and 4: MARTINGALES 1. Introduction in an Arbitrage–Free
    LECTURES 3 AND 4: MARTINGALES 1. Introduction In an arbitrage–free market, the share price of any traded asset at time t = 0 is the expected value, under an equilibrium measure, of its discounted price at market termination t = T . We shall see that the discounted share price at any time t ≤ T may also be computed by an expectation; however, this expectation is a conditional expectation, given the information about the market scenario that is revealed up to time t. Thus, the discounted share price process, as a (random) function of time, has the property that its value at any time t is the conditional expectation of the value at a future time given information available at t.A stochastic process with this property is called a martingale. The theory of martingales (initiated by Joseph Doob, following earlier work of Paul Levy´ ) is one of the central themes of modern probability. In discrete-time finance, the importance of martingales stems largely from the Optional Stopping Formula, which (a) places fundamental constraints on (discounted) share price processes, and (b) provides a computational tool of considerable utility. 2. Filtrations of Probability Spaces Filtrations. In a multiperiod market, information about the market scenario is revealed in stages. Some events may be completely determined by the end of the first trading period, others by the end of the second, and others not until the termination of all trading. This suggests the following classification of events: for each t ≤ T , (1) Ft = {all events determined in the first t trading periods}. The finite sequence (Ft)0≤t≤T is a filtration of the space Ω of market scenarios.
    [Show full text]
  • Martingales by D
    Martingales by D. Cox December 2, 2009 1 Stochastic Processes. Definition 1.1 Let T be an arbitrary index set. A stochastic process indexed by T is a family of random variables (Xt : t ∈ T) defined on a common probability space (Ω, F,P ). If T is clear from context, we will write (Xt). If T is one of ZZ, IN, or IN \{0}, we usually call (Xt) a discrete time process. If T is an interval in IR (usually IR or [0, ∞)), then we usually call (Xt) a continuous time process. In a sense, all of probability is about stochastic processes. For instance, if T = {1}, then we are just talking about a single random variable. If T = {1,...,n}, then we have a random vector (X1,...,Xn). We have talked about many results for i.i.d. random variables X1, X2, . Assuming an inifinite sequence of such r.v.s, T = IN \{0} for this example. Given any sequence of r.v.s X1, X2, . , we can define a partial sum process n Sn = Xi, n =1, 2,.... Xi=1 One important question that arises about stochastic processes is whether they exist or not. For example, in the above, can we really claim there exists an infinite sequence of i.i.d. random variables? The product measure theorem tells us that for any valid marginal distribution PX , we can construct any finite sequence of r.v.s with this marginal distribution. If such an infinite sequence of i.i.d. r.v.sr does not exist, we have stated a lot of meaniningless theorems.
    [Show full text]
  • A Course in Interacting Particle Systems
    A Course in Interacting Particle Systems J.M. Swart January 14, 2020 arXiv:1703.10007v2 [math.PR] 13 Jan 2020 2 Contents 1 Introduction 7 1.1 General set-up . .7 1.2 The voter model . .9 1.3 The contact process . 11 1.4 Ising and Potts models . 14 1.5 Phase transitions . 17 1.6 Variations on the voter model . 20 1.7 Further models . 22 2 Continuous-time Markov chains 27 2.1 Poisson point sets . 27 2.2 Transition probabilities and generators . 30 2.3 Poisson construction of Markov processes . 31 2.4 Examples of Poisson representations . 33 3 The mean-field limit 35 3.1 Processes on the complete graph . 35 3.2 The mean-field limit of the Ising model . 36 3.3 Analysis of the mean-field model . 38 3.4 Functions of Markov processes . 42 3.5 The mean-field contact process . 47 3.6 The mean-field voter model . 49 3.7 Exercises . 51 4 Construction and ergodicity 53 4.1 Introduction . 53 4.2 Feller processes . 54 4.3 Poisson construction . 63 4.4 Generator construction . 72 4.5 Ergodicity . 79 4.6 Application to the Ising model . 81 4.7 Further results . 85 5 Monotonicity 89 5.1 The stochastic order . 89 5.2 The upper and lower invariant laws . 94 5.3 The contact process . 97 5.4 Other examples . 100 3 4 CONTENTS 5.5 Exercises . 101 6 Duality 105 6.1 Introduction . 105 6.2 Additive systems duality . 106 6.3 Cancellative systems duality . 113 6.4 Other dualities .
    [Show full text]
  • A Tutorial on Hidden Markov Models and Selected Applications in Speech Recognition
    A Tutorial on Hidden Markov Models and Selected Applications in Speech Recognition LAWRENCE R. RABINER, FELLOW, IEEE Although initially introduced and studied in the late 1960s and In this case, with a good signal model, we can simulate the early 1970s, statistical methods of Markov source or hidden Markov source and learn as much as possible via simulations. modeling have become increasingly popular in the last several years. There are two strong reasons why this has occurred. First the Finally, the most important reason why signal models are models are very rich in mathematical structure and hence can form important is that they often workextremelywell in practice, the theoretical basis for use in a wide range of applications. Sec- and enable us to realize important practical systems-e.g., ond the models, when applied properly, work very well in practice prediction systems, recognition systems, identification sys- for several important applications. In this paper we attempt to care- fully and methodically review the theoretical aspects of this type tems, etc., in a very efficient manner. of statistical modeling and show how they have been applied to These are several possible choices for what type of signal selected problems in machine recognition of speech. model is used for characterizing the properties of a given signal. Broadly one can dichotomize the types of signal models into the class of deterministic models, and the class I. INTRODUCTION of statistical models. Deterministic models generally exploit Real-world processes generally produce observable out- some known specific properties of the signal, e.g., that the puts which can be characterized as signals.
    [Show full text]
  • Stochastic Processes: Theory and Methods. Handbook of Statistics 19 (Ed
    Stochastic Processes: Theory and Methods. Handbook of Statistics 19 (ed. C. R. Rao & D. N. Shanbhag) 171-213, Elsevier, 2001) RANDOM WALK AND FLUCTUATION THEORY N. H. BINGHAM CONTENTS PART I: RANDOM WALKS ON EUCLIDEAN SPACE x1. Introduction x2. Simple random walk on ZZ x3. Recurrence and transience x4. Simple random walk on ZZd; P¶olya's theorem x5. Random walk on IRd x6. Harmonic analysis x7. Potential theory x8. Coupling x9. Renewal theory x10. Limit theorems and Brownian motion x11. Conditioned random walk PART II: RANDOM WALK IN MORE GENERAL CONTEXTS x1. Random walk and electrical networks x2. Random walk on graphs x3. Random walk on groups x4. Brownian motion on Riemannian manifolds x5. Random walk on homogeneous spaces, Gelfand pairs, hypergroups and semigroups x6. Random walk on graphs with special structure x7. Variants on random walk PART III: FLUCTUATION THEORY x1. Spitzer's identity x2. Ladder epochs and heights x3. Spitzer's arc-sine law x4. Ballot theorems x5. Queues x6. Continuous time x7. Barrier problems x8. Fluctuation theory in higher dimensions and algebraic generalisations x9. Distribution-free results and non-parametric statistics x10. Postscript 1 PART I: RANDOM WALK ON EUCLIDEAN SPACE x1. Introduction. The theory and applications of random walks are ubiquitous in the modern prob- ability literature, and random walks form perhaps the simplest and most important exam- ples of stochastic processes - random phenomena unfolding with time. The term `random walk' can be traced back (at least) to P¶olya (1921) (`zufÄallige Irrfahrt' in the German), but the context is much older. If we think of X1;X2;::: as the gains or losses of a gambler on successive plays of a gambling game, and the partial sums Sn := X1 + ¢ ¢ ¢ + Xn (S0 := 0) as his cumulative gain or loss to date (by time n), then the behaviour of the stochastic 1 process (Sn)n=0 - the random walk with steps Xn - describes the evolution of the gambler's fortune.
    [Show full text]