An Introduction to Stochastic Control

Total Page:16

File Type:pdf, Size:1020Kb

An Introduction to Stochastic Control An introduction to stochastic control F. J. Silva Dip. di Matematica Guido Castelnuovo March 2012 Some disperse but useful results This section is based on [24, Chapter 1]. Let Ω be a nonempty set, and let F ⊂ 2Ω. We say that F is a π-system if A; B 2 F ) A \ B 2 F. F is a λ-system if (i) Ω 2 F; (ii) A, B 2 F and A ⊆ B imply B n A 2 F; (iii) Ai 2 F, Ai " A implies that A 2 F. Lemma 1. [σ-field of a π-system, or λ-π lemma] If a π system A is contained in a λ system F then σ(A) ⊆ F. Proof: See any standard book on measure theory (e.g. [2]) Example of application: [Uniqueness of the extension of a measure defined on a π-system] Let P and Q be two probability measures on (Ω; F) which coincide on a 1 π-system A, then they coincide on σ(A). In fact, it is enough to define C = fC 2 F ; P (C) = Q(C)g and two verify that it is a λ-system. Corollary 1. [Measurable function w.r.t. the σ-field of a π system] Let A be a π system. Let H be a linear space of functions from Ω to R such that (i)1 2 H;( ii) IA 2 H for all A 2 A (iii) φi 2 H; 0 ≤ φi " φ, φ is bounded ) φ 2 H: Then H contains all σ(A)-measurable functions from Ω to R. Proof: Let φ be σ(A)-measurable. Clearly we have that n + n X −n φ " φ ; with φ (!) := j2 Ifφ+(!)2[j2−n;(j+1)2−n]g (1) j≥0 with an analogous approximation for φ−. Therefore, by (iii), it is enough to show that φn 2 H. But φ is the sum of indicator of elements in σ(A). Therefore it is natural to consider the set F := fA 2 Ω; IA 2 Hg: 2 and to prove that σ(A) ⊆ F. But this follows from lemma 1, since A is a π-system and A ⊂ F, which is easily shown to be a λ-system. Theorem 1. [Dynkin theorem] Let (Ω; F) and (Ω0; F 0) two measurable spaces, and let (U; d) be a Polish space. Let ξ :Ω ! Ω0 and φ :Ω ! U be r.v.'s. Then φ 2 σ(ξ)- measurable, i.e. φ−1(B(U)) ⊆ ξ−1(F 0), iff there exists a measurable η :Ω0 ! U such that φ(!) = η(ξ(!)) for all ! 2 Ω: Proof: Consider the case U = R (the general case can be obtained using an isomorphism theorem, see [19]) and define the set 0 0 H := η(ξ); for some F -measurable map η :Ω ! R : We have to shown that the set of σ(ξ)=B(R) measurable maps belongs to H. This can be done by checking the assumptions of corollary 1 with A = σ(ξ) and proving that H satisfies (i), (ii) and (iii). Exercise: Do the details of the above proof. 3 Lemma 2. [Borel-Cantelli] Let (Ω; F; P) be a probability space and consider a sequence of events Ai 2 F. We have X P(Ai) < 1 ) P(\i [j≥i Aj) = 0: i Proof: Straightforward. Note that P(\i [j≥i Aj) ≤ P([j≥iAj) for all i and use the convergence of the series. Lemma 3. [Chebyshev inequality] Consider a nonnegative r.v. X. Then, for all p 2 (0; 1) and " > 0 we have p E(X ) P(X ≥ ") ≤ : "p Proof: It suffices to note that Z Z p p p X (X ≥ ") = (X ≥ " ) = p p d (!) ≤ d (!): P P IfX ≥" g P p P Ω Ω " 4 Conditional expectation This section is based on [24, Chapter 1]. Consider a probability space (Ω; F; P). Let 1 X 2 L (Ω) and G be a sub-σ-field of F. Define the signed measure µ : G! R as Z µ(A) := X(!)dP(!) for all A 2 G: A 1 By the Radon-Nikodym theorem there exists a unique PjG a.s. r.v. f 2 LG(Ω) (i..e. in particular G-measurable) such that Z Z fdP = XdP for all A 2 G: (2) A A The function f is called the conditional expectation of X given G, and we write E(XjG) := f: 5 Fundamental properties of E(·|G) All the properties below are simple consequences of (2) (Do them as an exercise!) (i) E(·|G) is a linear bounded functional. (ii) For a constant a 2 R, we have E(ajG) = a. 1 (iii) [Mononoticity] For X; Y 2 LF with X ≥ Y we have E(XjG) ≥ E(Y jG). 2 2 (iv) [Take out the measurable part] For X 2 LF and Y 2 LG, we have E(YXjG) = Y E(XjG). (v) [Characterization of independence] X is independent of G iff for every Borel f 1 such that f(X) 2 L (Ω) we have E(f(X)jG) = E(f(X)). (vi) [Tower or \projection" property] If G1 ⊆ G2 ⊆ F, then E(E(XjG1)jG2) = E(E(XjG2)jG1) = E(XjG1): (vii) [Jensen inequality] Let φ convex such that φ(X) 2 L1(Ω), then φ(E(XjG)) ≤ E(φ(X)jG): 6 [Conditioning one r.v. X w.r.t. another r.v. ξ] 1 Let X 2 LF and ξ : (Ω; F) ! (U; B(U)). Note that we can always define −1 E(Xjξ) := E(Xjξ (B(U))) = η(ξ) for some B(U)=B(R)- measurable function η; by Dynkin theorem. Therefore, it is natural to define E(Xjξ = x) := η(x): Another way to define this is by appealing to the Radon-Nikodym. In fact, let us define in B(U) the measure Z ν(B) := X(!)dP(!); ξ−1(B) −1 which is absolutely continuous with respect to Pξ := P ◦ ξ (the image measure of P dν under ξ). Therefore, there exists a Radon -Nikodym derivative (unique "-a.s.) such dP" P that Z dν Z (x)dP"(x) = X(!)dP(!) for all B 2 B(U): B dP" ξ−1(B) 7 We define dν E(Xjξ = x) := (x): dP" It can be checked that −1 dν η(ξ(!)) = E(Xjξ (B(U)))(!) = (ξ(!)) Pjξ−1(B(U)) − a.s. (3) dP" −1 Integrating w.r.t. a set of the form ξ (A) and using the definition of Pξ, we obtain that dν η(x) = (x) P" − a.s. dP" Note that incidentally, this gives another proof of Dynkin theorem. Let us prove (3). We 8 have R dν R dν (ξ(!))d (!) = (x)d "(x) ξ−1(B) dP" P B dP" P = ν(B) R = ξ−1(B) X(!)dP(!) R −1 = ξ−1(B) E(Xjξ (B(U)))(!)dP(!): which yields the result. Now, we define the conditional probability w.r.t. a σ-field G as P(AjG) := E(IAjG): Nota that for each B 2 F we have that P(BjG) is only defined PjG a.s. Thus, it is not sure that we can find some A 2 G with P(A) = 1 such that if we fix any ! 2 A, we can measure all the sets B 2 F . However, we can give a sense to this using the concept of regular conditional probability (see [2] for more on this). 9 [Characterization of E(·|ξ) in terms of fg(ξ); g bounded continuousg] Let us now prove that E(Xjξ) = 0 iff for all bounded continuous g we have E(g(ξ)X) = 0: The \only if" part is direct. To prove the \if part", first note that E(Xjξ) = 0 iff E(φX) = 0 for all φ that are σ(ξ) measurable. Now, consider the set H := fφ : (Ω; F) ! (R; B(R)) ; E(φX) = 0g: We have to show that H contains the σ(ξ)-measurable functions. It is clear that H satisfies the assumptions (i) and (iii) of corollary 1. We only have to construct a π-system A such that IA 2 H for all A 2 A and σ(A) = σ(ξ). Let us take A := fξ−1([a; b]) ; for some a < bg: 10 If A = ξ−1([a; b]) 2 A we have E(IA(ξ)X) = E(I[a;b](ξ)X): n Now, take any sequence of continuous functions g ! I[a;b](ξ), where the convergence is n pointwise. Since E(g (ξ)X) = 0, by passing to the limit we get that E(IA(ξ)X) = 0 and so assumption (ii) of corollary 1 is verified. Note that the same proof yields that if G = σ(ξ1; :::; ξn) we have n E(XjG) = 0 iff for all bounded continuous g : R ! R we have E(g(ξ1; : : : ; ξn)X) = 0: The interesting fact is that the result is also valid when we condition to a countable set of r.v's. More precisely, using the above result and the technique of its proof we get the following result. Proposition 1. [Checking conditions only on a finite set of variables] Consider a 11 sequence of variables ξ1; ξ2;::: and define G := σ(fξi ; i 2 Ng). Then n E(XjG) = 0 iff for all n 2 N and g 2 Cb(R ) we have E(g(ξ1; : : : ; ξn)X) = 0: 12 Stochastic process: Basic definitions Good references for this part are [7, 14]. Let I be a non-empty index set and (Ω; F; P) a probability space.
Recommended publications
  • STOCHASTIC CONTROL Contents 1. Introduction 1 2. Mathematical Probability 1 3. Stochastic Calculus 3 4. Diffusions 7 5. the Stoc
    STOCHASTIC CONTROL NAFTALI HARRIS Abstract. In this paper we discuss the Stochastic Control Problem, which deals with the best way to control a stochastic system and has applications in fields as diverse as robotics, finance, and industry. After quickly reviewing mathematical probability and stochastic integration, we prove the Hamilton- Jacobi-Bellman equations for stochastic control, which transform the Stochas- tic Control Problem into a partial differential equation. Contents 1. Introduction 1 2. Mathematical Probability 1 3. Stochastic Calculus 3 4. Diffusions 7 5. The Stochastic Control Problem 11 6. The Hamilton-Jacobi-Bellman Equation 12 7. Acknowledgements 18 References 18 1. Introduction The Stochastic Control Problem arises when we wish to control a continuous- time random system optimally. For example, we wish to fly an airplane through turbulence, or manage a stock portfolio, or precisely maintain the temperature and pressure of a system in the presence of random air currents and temperatures. To build up to a formal definition of this problem, we will first review mathematical probability and stochastic calculus. These will give us the tools to talk about It^o processes and diffusions, which will allow us to state the Stochastic Control Problem formally. Finally, we will spend the latter half of this paper proving necessary conditions and sufficient conditions for its solution. 2. Mathematical Probability In this section we briskly define the mathematical model of probability required in the rest of this paper. We define probability spaces, random variables, and conditional expectation. For more detail, see Williams [4]. Definition 2.1 (Probability Space). A probability space is a triple (Ω; F;P ), where Ω is an arbitrary set, F is a sigma algebra over Ω, and P is a measure on sets in F with P (Ω) = 1.
    [Show full text]
  • MARKOV DECISION PROCESSES Abstract: the Theory of Markov
    July 06, 2010, 15:00 Draft, don’t distribute MARKOV DECISION PROCESSES NICOLE BAUERLE¨ ∗ AND ULRICH RIEDER‡ Abstract: The theory of Markov Decision Processes is the theory of controlled Markov chains. Its origins can be traced back to R. Bellman and L. Shapley in the 1950’s. During the decades of the last century this theory has grown dramatically. It has found applications in various areas like e.g. computer science, engineering, operations research, biology and economics. In this article we give a short intro- duction to parts of this theory. We treat Markov Decision Processes with finite and infinite time horizon where we will restrict the presentation to the so-called (gen- eralized) negative case. Solution algorithms like Howard’s policy improvement and linear programming are also explained. Various examples show the application of the theory. We treat stochastic linear-quadratic control problems, bandit problems and dividend pay-out problems. AMS 2010 Classification: 90C40, 60J05, 93E20 Keywords and Phrases: Markov Decision Process, Markov Chain, Bell- man Equation, Policy Improvement, Linear Programming 1. Introduction Do you want to play a card game? Yes? Then I will tell you how it works. We have a well-shuffled standard 32-card deck which is also known as a piquet deck. 16 cards are red and 16 cards are black. Initially the card deck lies on the table face down. Then I start to remove the cards and you are able to see its faces. Once you have to say ”stop”. If the next card is black you win 10 Euro, if it is red you loose 10 Euro.
    [Show full text]
  • On Some Recent Aspects of Stochastic Control and Their Applications 507
    Probability Surveys Vol. 2 (2005) 506–549 ISSN: 1549-5787 DOI: 10.1214/154957805100000195 On some recent aspects of stochastic control and their applications∗ Huyˆen Pham† Laboratoire de Probabilit´es et Mod`eles Al´eatoires, CNRS, UMR 7599 Universit´eParis 7, e-mail: [email protected], and CREST Abstract: This paper is a survey on some recent aspects and developments in stochastic control. We discuss the two main historical approaches, Bell- man’s optimality principle and Pontryagin’s maximum principle, and their modern exposition with viscosity solutions and backward stochastic differ- ential equations. Some original proofs are presented in a unifying context including degenerate singular control problems. We emphasize key results on characterization of optimal control for diffusion processes, with a view towards applications. Some examples in finance are detailed with their ex- plicit solutions. We also discuss numerical issues and open questions. AMS 2000 subject classifications: 93E20, 49J20, 49L20, 60H30. Keywords and phrases: Controlled diffusions, dynamic programming, maximum principle, viscosity solutions, backward stochastic differential equations, finance. Received December 2005. arXiv:math/0509711v3 [math.PR] 12 Jan 2006 ∗This is an original survey paper †I would like to thank Philip Protter for his encouragement to write this paper. 506 H. Pham/On some recent aspects of stochastic control and their applications 507 Contents 1 Introduction 508 2 Theproblemanddiscussiononmethodology 508 2.1 Problemformulation . 508 2.2 Bellman’s optimality principle . 509 2.2.1 The Hamilton-Jacobi-Bellman (in)equation . 511 2.2.2 Theclassicalverificationapproach . 512 2.3 Pontryagin’s maximum principle . 514 2.4 Othercontrolproblems . 515 2.4.1 Randomhorizon .
    [Show full text]
  • Stochastic Control
    1 Stochastic Control Roger Brockett April 3, 2009 2 PREFACE Teaching stochastic processes to students whose primary interests are in applications has long been a problem. On one hand, the subject can quickly become highly technical and if mathematical concerns are allowed to dominate there may be no time available for exploring the many interesting areas of applications. On the other hand, the treatment of stochastic calculus in a cavalier fashion leaves the student with a feeling of great uncertainty when it comes to exploring new material. Moreover, the problem has become more acute as the power of the differential equation point of view has become more widely appreciated. In these notes, an attempt is made to resolve this dilemma with the needs of those interested in building models and designing algorithms for estimation and control in mind. The approach is to start with Poisson counters and to identify the Wiener process with a certain limiting form. We do not attempt to define the Wiener process per se. Instead, everything is done in terms of limits of jump processes. The Poisson counter and differential equations whose right-hand sides include the differential of Poisson counters are developed first. This leads to the construction of a sample path representations of a continuous time jump process using Poisson counters. This point of view leads to an efficient problem solving technique and permits a unified treatment of time varying and nonlinear problems. More importantly, it provides sound intuition for stochastic differential equations and their uses without allowing the technicalities to dominate. In treating estimation theory, the conditional density equa- tion is given a central role.
    [Show full text]
  • Stochastic Bridges of Linear Systems Yongxin Chen and Tryphon Georgiou
    1 Stochastic bridges of linear systems Yongxin Chen and Tryphon Georgiou Abstract We study a generalization of the Brownian bridge as a stochastic process that models the position and velocity of inertial particles between the two end-points of a time interval. The particles experience random acceleration and are assumed to have known states at the boundary. Thus, the movement of the particles can be modeled as an Ornstein-Uhlenbeck process conditioned on position and velocity measurements at the two end-points. It is shown that optimal stochastic control provides a stochastic differential equation (SDE) that generates such a bridge as a degenerate diffusion process. Generalizations to higher order linear diffusions are considered. I. INTRODUCTION The theoretical foundations on how molecular dynamics affect large scale properties of ensembles were layed down more than a hundred years ago. A most prominent place among mathematical models has been occupied by the Brownian motion which provides a basis for studying diffusion and noise [1], [2], [3], [4]. The Brownian motion is captured by the mathematical model of a Wiener process, herein denoted by w(t). It represents the random motion of particles suspended in a fluid where their inertia is negligible compared to viscous forces. Taking into account inertial effects under a “delta-correlated” stationary Gaussian force field η(t) (that is, white noise, loosely thought of as dw=dt [1, p. 46]) d2x(t) dx(t) m = −λ + η(t) dt2 dt represents the Langevin dynamics ; x represents position, m mass, t time, and λ viscous friction parameter. The corresponding SDE dx(t) 0 1 x(t) 0 = dt + dw(t); dv(t) 0 −λ/m v(t) 1=m where w is a Wiener process and v the velocity, is a degenerate diffusion in that the stochastic term does not affect all degrees of freedom.
    [Show full text]
  • Numerical Methods for Mean Field Games and Mean Field Type Control
    Numerical Methods for Mean Field Games and Mean Field Type Control Mathieu Lauri`ere Abstract. Mean Field Games (MFG) have been introduced to tackle games with a large number of competing players. Considering the limit when the number of players is infinite, Nash equilibria are studied by considering the interaction of a typical player with the population's distribution. The situa- tion in which the players cooperate corresponds to Mean Field Control (MFC) problems, which can also be viewed as optimal control problems driven by a McKean-Vlasov dynamics. These two types of problems have found a wide range of potential applications, for which numerical methods play a key role since most models do not have analytical solutions. In these notes, we review several aspects of numerical methods for MFG and MFC. We start by pre- senting some heuristics in a basic linear-quadratic setting. We then discuss numerical schemes for forward-backward systems of partial differential equa- tions (PDEs), optimization techniques for variational problems driven by a Kolmogorov-Fokker-Planck PDE, an approach based on a monotone operator viewpoint, and stochastic methods relying on machine learning tools. Contents 1. Introduction 2 1.1. Outline 3 1.2. Definition of the problems and notation 3 2. Warm-up: Linear-Quadratic problems 5 2.1. Problem definition and MFG solution 5 2.2. Time discretization 6 2.3. Picard (fixed point) iterations 6 arXiv:2106.06231v1 [math.OC] 11 Jun 2021 2.4. Fictitious play 9 2.5. Newton iterations 10 2.6. Mean Field Type Control and the Price of Anarchy 10 3.
    [Show full text]
  • Markov Chains
    Chapter 1 Markov Chains A sequence of random variables X0,X1,...with values in a countable set S is a Markov chain if at any time n, the future states (or values) Xn+1,Xn+2,... depend on the history X0,...,Xn only through the present state Xn.Markov chains are fundamental stochastic processes that have many diverse applica- tions. This is because a Markov chain represents any dynamical system whose states satisfy the recursion Xn = f(Xn−1,Yn), n ≥ 1, where Y1,Y2 ... are independent and identically distributed (i.i.d.) and f is a deterministic func- tion. That is, the new state Xn is simply a function of the last state and an auxiliary random variable. Such system dynamics are typical of those for queue lengths in call centers, stresses on materials, waiting times in produc- tion and service facilities, inventories in supply chains, parallel-processing software, water levels in dams, insurance funds, stock prices, etc. This chapter begins by describing the basic structure of a Markov chain and how its single-step transition probabilities determine its evolution. For in- stance, what is the probability of reaching a certain state, and how long does it take to reach it? The next and main part of the chapter characterizes the stationary or equilibrium distribution of Markov chains. These distributions are the basis of limiting averages of various cost and performance param- eters associated with Markov chains. Considerable discussion is devoted to branching phenomena, stochastic networks, and time-reversible chains. In- cluded are examples of Markov chains that represent queueing, production systems, inventory control, reliability, and Monte Carlo simulations.
    [Show full text]
  • Applications of Stochastic Optimal Control to Economics and Finance
    A Service of Leibniz-Informationszentrum econstor Wirtschaft Leibniz Information Centre Make Your Publications Visible. zbw for Economics Federico, Salvatore (Ed.); Ferrari, Giorgio (Ed.); Regis, Luca (Ed.) Book — Published Version Applications of stochastic optimal control to economics and finance Provided in Cooperation with: MDPI – Multidisciplinary Digital Publishing Institute, Basel Suggested Citation: Federico, Salvatore (Ed.); Ferrari, Giorgio (Ed.); Regis, Luca (Ed.) (2020) : Applications of stochastic optimal control to economics and finance, ISBN 978-3-03936-059-8, MDPI, Basel, http://dx.doi.org/10.3390/books978-3-03936-059-8 This Version is available at: http://hdl.handle.net/10419/230540 Standard-Nutzungsbedingungen: Terms of use: Die Dokumente auf EconStor dürfen zu eigenen wissenschaftlichen Documents in EconStor may be saved and copied for your Zwecken und zum Privatgebrauch gespeichert und kopiert werden. personal and scholarly purposes. Sie dürfen die Dokumente nicht für öffentliche oder kommerzielle You are not to copy documents for public or commercial Zwecke vervielfältigen, öffentlich ausstellen, öffentlich zugänglich purposes, to exhibit the documents publicly, to make them machen, vertreiben oder anderweitig nutzen. publicly available on the internet, or to distribute or otherwise use the documents in public. Sofern die Verfasser die Dokumente unter Open-Content-Lizenzen (insbesondere CC-Lizenzen) zur Verfügung gestellt haben sollten, If the documents have been made available under an Open gelten abweichend von
    [Show full text]
  • Machine Learning, Lecture 5 Kernel Machines
    Outline lecture 5 2(39) Machine Learning, Lecture 5 1. Summary of lecture 4 Kernel Machines 2. Introductory GP example 3. Stochastic processes Thomas Schön 4. Gaussian processes (GP) Division of Automatic Control Construct a GP from a Bayesian linear regression model • GP regression Linköping University • Examples where we have made use of GPs in recent research Linköping, Sweden. • 5. Support vector machines (SVM) Email: [email protected], Phone: 013 - 281373, Chapter 6.4 – 7.2 (Chapter 12 in HTF, GP not covered in HTF) Office: House B, Entrance 27. AUTOMATIC CONTROL Machine Learning REGLERTEKNIK Machine Learning T. Schön LINKÖPINGS UNIVERSITET T. Schön Summary of lecture 4 (I/II) 3(39) Summary of lecture 4 (II/II) 4(39) A kernel function k(x, z) is defined as an inner product A neural network is a nonlinear function (as a function expansion) from a set of input variables to a set of output variables controlled by k(x, z) = φ(x)Tφ(z), adjustable parameters w. where φ(x) is a fixed mapping. This function expansion is found by formulating the problem as usual, which results in a (non-convex) optimization problem. This problem is Introduced the kernel trick (a.k.a. kernel substitution). In an solved using numerical methods. algorithm where the input data x enters only in the form of scalar products we can replace this scalar product with another choice of Backpropagation refers to a way of computing the gradients by kernel. making use of the chain rule, combined with clever reuse of information that is needed for more than one gradient.
    [Show full text]
  • Mean Field Simulation for Monte Carlo Integration MONOGRAPHS on STATISTICS and APPLIED PROBABILITY
    Mean Field Simulation for Monte Carlo Integration MONOGRAPHS ON STATISTICS AND APPLIED PROBABILITY General Editors F. Bunea, V. Isham, N. Keiding, T. Louis, R. L. Smith, and H. Tong 1. Stochastic Population Models in Ecology and Epidemiology M.S. Barlett (1960) 2. Queues D.R. Cox and W.L. Smith (1961) 3. Monte Carlo Methods J.M. Hammersley and D.C. Handscomb (1964) 4. The Statistical Analysis of Series of Events D.R. Cox and P.A.W. Lewis (1966) 5. Population Genetics W.J. Ewens (1969) 6. Probability, Statistics and Time M.S. Barlett (1975) 7. Statistical Inference S.D. Silvey (1975) 8. The Analysis of Contingency Tables B.S. Everitt (1977) 9. Multivariate Analysis in Behavioural Research A.E. Maxwell (1977) 10. Stochastic Abundance Models S. Engen (1978) 11. Some Basic Theory for Statistical Inference E.J.G. Pitman (1979) 12. Point Processes D.R. Cox and V. Isham (1980) 13. Identification of OutliersD.M. Hawkins (1980) 14. Optimal Design S.D. Silvey (1980) 15. Finite Mixture Distributions B.S. Everitt and D.J. Hand (1981) 16. ClassificationA.D. Gordon (1981) 17. Distribution-Free Statistical Methods, 2nd edition J.S. Maritz (1995) 18. Residuals and Influence in RegressionR.D. Cook and S. Weisberg (1982) 19. Applications of Queueing Theory, 2nd edition G.F. Newell (1982) 20. Risk Theory, 3rd edition R.E. Beard, T. Pentikäinen and E. Pesonen (1984) 21. Analysis of Survival Data D.R. Cox and D. Oakes (1984) 22. An Introduction to Latent Variable Models B.S. Everitt (1984) 23. Bandit Problems D.A. Berry and B.
    [Show full text]
  • Hidden Markov Model Multiarm Bandits: a Methodology for Beam Scheduling in Multitarget Tracking Vikram Krishnamurthy, Senior Member, IEEE, and Robin J
    IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 49, NO. 12, DECEMBER 2001 2893 Hidden Markov Model Multiarm Bandits: A Methodology for Beam Scheduling in Multitarget Tracking Vikram Krishnamurthy, Senior Member, IEEE, and Robin J. Evans Abstract—In this paper, we derive optimal and suboptimal beam scheduling algorithms for electronically scanned array tracking systems. We formulate the scheduling problem as a multiarm bandit problem involving hidden Markov models (HMMs). A finite-dimensional optimal solution to this multiarm bandit problem is presented. The key to solving any multiarm bandit problem is to compute the Gittins index. We present a finite-dimensional algorithm that computes the Gittins index. Suboptimal algorithms for computing the Gittins index are also presented. Numerical examples are presented to illustrate the algorithms. Index Terms—Dynamic programming, hidden Markov models, optimal beam steering, scheduling, sequential decision procedures. Fig. 1. Multitarget tracking with one intelligent sensor. The scheduling problem involves which single target to track with the steerable electronic scanned array (ESA). I. INTRODUCTION AND PROBLEM FORMULATION ONSIDER the following hidden Markov model (HMM) work considers one important part of a complete agile beam C target tracking problem: targets (e.g., aircraft) are tracking system. We show that the multiarm bandit formulation being tracked by an electronically scanned array (ESA) with can provide a rigorous and potentially useful approach to the only one steerable beam. The coordinates of each target problem of scheduling an agile beam among a finite number of evolve according to a finite state Markov established target tracks. chain (see [27]), and the conditional mean estimate target state The multiarm bandit problem is an example of a dynamic sto- is computed by an HMM tracker.
    [Show full text]
  • A Framework of Stochastic Power Management Using Hidden Markov Model
    A Framework of Stochastic Power Management Using Hidden Markov Model Ying Tan, Qinru Qiu Department of Electrical and Computer Engineering Binghamton University, State University of New York Binghamton, New York 13902, USA {ying, qqiu}@binghamton.edu Abstract - The effectiveness of stochastic power management correspondence because the workload is not only controlled by the relies on the accurate system and workload model and effective hardware and software but also affected by user and environment. policy optimization. Workload modeling is a machine learning For example, user working style and user mood have significant procedure that finds the intrinsic pattern of the incoming tasks impact of the workload distribution of a computer system. Consider based on the observed workload attributes. Markov Decision a power manager (PM) of a wireless adapter. The observed Process (MDP) based model has been widely adopted for workload attribute is the frequency and the size of the stochastic power management because it delivers provable incoming/outgoing TCP/IP packets. During a teleconference, the optimal policy. Given a sequence of observed workload user may choose to send a video image, send the audio message, attributes, the hidden Markov model (HMM) of the workload share files or share the image on whiteboard. Different operations is trained. If the observed workload attributes and states in the generate different communication requirements. An accurate workload model do not have one-to-one correspondence, the workload model of the wireless adapter must reflect how the user MDP becomes a Partially Observable Markov Decision Process switches from one operation to another. However, this information (POMDP). This paper presents a framework of modeling and is not observable by the PM.
    [Show full text]