Controlling Systemic Inflammation Using NMPC. Using Nonlinear Model Predictive Control with State Estimation
|
|
- Kory Austin
- 5 years ago
- Views:
Transcription
1 Controlling Systemic Inflammation Using Nonlinear Model Predictive Control with State Estimation Gregory Zitelli, Judy Day July 2013 With generous support from the NSF, Award
2 Motivation We re going to be interested in the acute inflammatory response to biological stress in the form of a bacterial pathogen. The inflammatory response aims to remove the presence of the pathogen. However, an excessive response may lead to collateral tissue damage, organ failure, or worse.
3 Motivation To recover from a high-inflammation state, it is necessary for the inflammatory mechanisms to downregulate itself through the use of some anti-inflammatory mediator. We will consider a mathematical formulation of these ideas in terms of a highly coupled, nonlinear ODE model, as well as control applied to both the pro and anti-inflammatory influences within the system.
4 Interaction Diagram
5 Nonlinear Model P represents the bacterial pathogen population. N represents the concentration of pro-inflammatory mediators, such as activated phagocytes and their produced cytokines. D acts as a marker of tissue damage. C A represents the concentration of anti-inflammatory mediators.
6 Nonlinear Model ( dp dt = k pgp 1 P ) k pms m P P µ m + k mp P k pnf(n )P dn dt = s nrr(p, N, D) µ nr + R(P, N, D) µ nn dd dt = k f(n ) 6 dn x 6 dn + f(n ) 6 µ dd dc A dt = s c + k cn f(n + k cnd D) 1 + f(n + k cnd D) µ cc A
7 Nonlinear Model The system (P, N, D, C A ) has three fixed points corresponding to the three biologically relevant scenarios: P = N = D = 0 and C A = sc µ c, where the patient is healthy. All states elevated, where the patient is septic. P = 0, and N, D, C A > 0, where the patient is aseptic. For values of k pg (the pathogen growth rate) in the interval (0.5137, 1.755), all three states are stable. (Reynolds et al 2006)
8 ??? So what does this model actually look like?
9 Septic Simulation P (state ( )) 6 4 N* (state ( )) Time Time D (state ( )) Time C a (state ( )) Time
10 Aseptic Simulation P (state ( ),est( )) Time N* (state ( ),est( )) Time D (state ( ),est( )) C a (state ( ),est( )) Time Time
11 Healthy Simulation P (state ( )) N* (state ( )) Time Time D (state ( )) C a (state ( )) Time Time
12 Goals We would like to apply some kind of control to this model, to direct septic and aseptic patients towards the healthy fixed state. This is accomplished using Nonlinear Model Predictive Control, or NMPC. The NMPC is applied to the pro and anti-inflammatory mediators, N and C A. We only apply positive control, and there are restrictions on how much can be introduced within certain time frames. We assume that the level of pro and anti-inflammatory mediators, N and C A, are measurable with Gaussian noise.
13 Overview of NMPC
14 Overview of NMPC Our reference trajectory is zero levels of P and D, and minimal amounts of control. J = min AI(t),P I(t) Γ DD Γ P P Γ AI AI Γ P I P I 2 2
15 Note on Observations We take noisy measurements of N and C A at each time interval which update our predictive model. The levels of P and D are not measured directly, and are instead estimated from the predictive model. Being unable to measure P often leads the NMPC to be too aggressive. Many virtual patients were unnecessarily harmed under this scheme. Since there are biologically relevant scenarios for loose measurements of P (indicators like body temperature or blood pressure give an idea whether or not the infection persists), a pathogen update is done every four timesteps.
16 Pathogen Update If the predictive model reads low (< 0.05) but the pathogen levels are high (P > 0.5) then the level in the predictive model is reset to 0.5. If the predictive model reads high (> 0.5) but the pathogen levels are low (P < 0.05) then the level in the predictive model is reset to 0.
17 Virtual Patient Pool 1,000 patients are randomly generated with different parameters, including unique values of the pathogen growth rate k pg. 620 acquire elevated P values to suggest treatment. Of the 620, 251 (40%) will return to a healthy state on their own.
18 Initial Results (Patient Snapshot) P (state ( ),est( )) N* (state ( ),est( )) Time Time D (state ( ),est( )) C a (state ( ),est( )) Time Time
19 Initial Results (Day et al 2010) Therapy Type: Placebo Mismatch k pg = 0.52 Mismatch k pg = 0.6 Mismatch k pg = 0.8 Percentage Healthy: 40% (251) 60% (369) 82% (510) 83% (513) Percentage Aseptic: 37% (228) 19% (120) 8% (49) 17% (107) Percentage Septic: 23% (141) 21% (131) 10% (61) 0% (0) Percentage Harmed: na 0% (0/251) 1% (2/251) 6% (16/251) Percentage Rescued: na 32% (118/369) 71% (261/369) 75% (278/369)
20 Particle Filter One way to go about improving these results is to impliment robust state estimation for the unobserved variables P and D. This is accomplished using a particle filter, which tracks its accuracy by comparing its predictions of N and C A to the observed values.
21 Particle Filter The particle filter is initialized with 1,000 particles randomized near our measurements for N and C A and initial predictions for P and D. Each particle p i = (P i, N i, C A,i, D i ) is simulated for one time step. At the next time step, the particles p i are assigned weights q i depending on how close N i, C A,i are to the new measurements for N and C A.
22 Particle Filter At the end of the time step, each of the 1,000 slots holding a particle is randomly assigned a new particle p i according to their weights q i. This causes bad particles with low weights q i to die off, while good particle with high weights replicate.
23 Particle Filter Results Therapy Type: Original k pg = 0.52 Original k pg = 0.6 Original k pg = 0.8 Percentage Healthy: 60% (369) 80% (497) 84% (519) Percentage Aseptic: 19% (121) 10% (60) 16% (101) Percentage Septic: 21% (130) 10% (63) 0% (0) Therapy Type: Particles k pg = 0.52 Particles k pg = 0.6 Particles k pg = 0.8 Percentage Healthy: 60% (369) 80% (493) 82% (511) Percentage Aseptic: 19% (120) 10% (66) 18% (109) Percentage Septic: 21% (131) 10% (61) 0% (0)
24 Particle Filter Results The particle filter does slightly worse than our original mismatched predictive model. On the other hand, the particle filter is self correcting, so it needs no pathogen update to correct a misled pathogen prediction.
25 Controllability Consider the general nonlinear system with affine control ẏ = F (y) + w = h(y) m g i (x)u i i=1 where f : R n R n and h : R n R k. We typically take u : [0, ) R m locally integrable.
26 Controllability Rank Condition ẏ = F (y) + w = h(y) m g i (x)u i i=1 Definition (Controllability Rank Condition) Given the same nonlinear system with affine control, we say that the system satisfies the controllability rank condition if a finite sub-matrix of the following matrix has rank n [ g1... g m [F, g 1 ]... [F, [F, g 1 ]]... ]
27 Controllability If the system is linear, Controllability Rank Condition = Locally Controllable If the system is nonlinear, Controllability Rank Condition = Locally Accessible If the system is nonlinear and has symmetric control (for each control u there is a u such that F (y, u) = F (y, u ) for every y R n ), then Controllability Rank Condition = Locally Controllable
28 Controllability It can be shown that our system does satisfy the controllability rank condition in the first octant. This means that the system is locally accessible there. Unfortunately, our control does not satisfy the symmetry property, so this does not imply that the system is locally controllable.
29 What s Next? Djouadi and Bara are working on the use of adaptive control for this model, to help overcome the mismatch between the predictive model and patient model parameters like k pg. Although the pathogen update was necessary for the variable P, it appears that the damage variable D is well estimated by our predictive model under many different initial conditions and model parameters. It would be nice to make this notion precise. Are there elements of this model that can be exploited to increase the effectiveness of our filter?
30 The End? Thank you very much!
31 Model ( dp dt = k pgp 1 P ) k pms m P P µ m + k mp P k pnf(n )P dn dt = s nrr(p, N, D) µ nr + R(P, N, D) µ nn dd dt = k f(n ) 6 dn x 6 dn + f(n ) 6 µ dd dc A dt = s c + k cn f(n + k cnd D) 1 + f(n + k cnd D) µ cc A
32 Model Where R(P, N, D) = f(k np P + k nn N + k nd D) f(x) = x 1 + ( CA c ) 2 c is chosen so that when C A is at its highest, f(x) 1 4 x.
33 Pathogen The bacterial pathogen follows logistic growth. Healthy individuals have a baseline, non-specific, local immune response. We assume through the initial analysis of this interaction that the local response reaches a quasi-steady state value of k pm µ m+k mpp. The pathogen are directly attacked by the phagocytic immune cells N, which may be inhibited by the presence of the anti-inflammatory C A. ( dp dt = k pgp 1 P ) k pms m P P µ m + k mp P k pnf(n )P
34 Pro-Inflammatory Mediator The first term comes from the assumed quasi-steady state of a subsystem involving resting and activated phagocytes. The second term represents the gradual decay in the concentration of pro-inflammatory mediators. dn dt = s nrr(p, N, D) µ nr + R(P, N, D) µ nn
35 Damage Collateral tissue damage from activated phagocytes motivated the first term. The damage saturates, with the large Hill coefficient necessary to produce a more realistic basin of attraction for the healthy fixed state. The second term represents tissue regeneration. dd dt = k f(n ) 6 dn x 6 dn + f(n ) 6 µ dd
36 Anti-Inflammatory Mediators The initial biological source term s c is augmented by the Michaelis-Menten term. The second term represents the gradual decay in the concentration of anti-inflammatory mediators. dc A dt = s c + k cn f(n + k cnd D) 1 + f(n + k cnd D) µ cc A
37 Essential Components of NMPC (I) Reference Trajectory Desired trajectory we want the system to approach. Our reference trajectory will be P = D = 0. (II) Prediction of Process Output A predictive model is used to provide an estimation of the states which are unobserved. (III) Objective Function An objective function will describe how are current predictions deviate from the reference trajectory. Our objective function is a weighted sum of squares for the undesireable variables P and D, as well as for the control levels. J = min AI(t),P I(t) Γ DD Γ P P Γ AI AI Γ P I P I 2 2
38 Essential Components of NMPC (IV) Sequence of Control Moves at Each Time Step A discrete control is designed for k future steps to minimize the objective function over each time interval. (V) Error Prediction Update The measured values of N and C A from the patient model are compared to the predictive model after the control has been applied to both. It is at this step that implementations of the particle filter can be refined through this comparison.
39 Controllability Consider the general nonlinear system with affine control ẏ = F (y) + w = h(y) m g i (x)u i i=1 where f : R n R n and h : R n R k. We typically take u : [0, ) R m locally integrable. For us, k pg y 1 (1 y 1 F (y) = P ) kpmsmy 1 µ m+k mpy 1 k pn f(y 2 )y 1 s nrr(y 1,y 2,y 3 ) µ nr+r(y 1,y 2,y 3 ) µ ny 2 k dn f(y 2 ) 6 x 6 dn +f(y 2) 6 µ d y 3 s c + k cn f(y 2 +k cnd y 3 ) 1+f(y 2 +k cnd y 3 ) µ cy 4
40 Controllability Definition We say that a point a R n is attainable from ŷ if there exists an appropriate control u y (locally integrable) such that the trajectory given by m ẏ = F (y) + g i (x)u y i i=1 y(0) = ŷ is such that y(t ) = a for some finite trajectory [0, T ].
41 Controllability Definition If ŷ R n and O is an open set containing ŷ, we let A O (ŷ, T ) denote the set of all points in O which are attainable from ŷ in time T, such that y(t) O for all t [0, T ]. Definition We say that the nonlinear system is locally accessible at a point ŷ R n if the sets A O (ŷ, T ) has nonempty interior for every open set O and every T > 0.
42 Controllability Definition Given the same nonlinear system with affine control, we say that the system satisfies the controllability rank condition if a finite sub-matrix of the following matrix has rank n [ g1... g m [F, g 1 ]... [F, [F, g 1 ]]... ]
43 Controllability If our system was linear, then it turns out that the controllability rank condition implies that the system is controllable. In the nonlinear case, we have the following results. Theorem Given the same nonlinear system with affine control, if the system satisfies the controllability rank condition at a point ŷ R n then it is locally accessible there. Theorem If for a nonlinear system ẏ = F (y, u) we have that for any locally integrable control u there exists another locally integrable control u such that F (y, u) = F (y, u ) for all y R n and the system satisfies the controllability rank condition at a point ŷ R n, then the system is locally controllable at ŷ.
44 Controllability Recall our system is in the following form k pg y 1 (1 y 1 kpmsmy 1 µ m+k mpy 1 k pn f(y 2 )y 1 s nrr(y 1,y 2,y 3 ) µ ẏ = nr+r(y 1,y 2,y 3 ) µ ny 2 + u 1 f(y k 2 ) 6 dn x 6 dn +f(y µ 2) 6 d y 3 f(y s c + k 2 +k cnd y 3 ) cn 1+f(y 2 +k cnd y 3 ) µ cy 4 + u u u 2 P )
Kernel-based Characterization of Dynamics in a Heterogeneous Population of Septic Patients Under Therapy
Kernel-based Characterization of Dynamics in a Heterogeneous Population of Septic Patients Under Therapy Kosta Ristovski a* kosta.ristovski@temple.edu Vladan Radosavljevic a* vladan@temple.edu Zoran Obradovic
More informationModeling the Immune System W9. Ordinary Differential Equations as Macroscopic Modeling Tool
Modeling the Immune System W9 Ordinary Differential Equations as Macroscopic Modeling Tool 1 Lecture Notes for ODE Models We use the lecture notes Theoretical Fysiology 2006 by Rob de Boer, U. Utrecht
More informationA mathematical and computational model of necrotizing enterocolitis
A mathematical and computational model of necrotizing enterocolitis Ivan Yotov Department of Mathematics, University of Pittsburgh McGowan Institute Scientific Retreat March 10-12, 2008 Acknowledgment:
More informationComputational model of NOS/TGF- Beta1/Plasmodia System in humans and mosquitoes
Computational model of NOS/TGF- Beta1/Plasmodia System in humans and mosquitoes Neil Parikh, Rice University Mentor: G. Bard Ermentrout, University of Pittsburgh Malaria One of the world s top 10 deadliest
More information. For each initial condition y(0) = y 0, there exists a. unique solution. In fact, given any point (x, y), there is a unique curve through this point,
1.2. Direction Fields: Graphical Representation of the ODE and its Solution Section Objective(s): Constructing Direction Fields. Interpreting Direction Fields. Definition 1.2.1. A first order ODE of the
More informationM469, Fall 2010, Practice Problems for the Final
M469 Fall 00 Practice Problems for the Final The final exam for M469 will be Friday December 0 3:00-5:00 pm in the usual classroom Blocker 60 The final will cover the following topics from nonlinear systems
More informationParameter Synthesis in Nonlinear Dynamical Systems: Application to Systems Biology
Parameter Synthesis in Nonlinear Dynamical Systems: Application to Systems Biology Alexandre Donzé 1, Gilles Clermont, Axel Legay 1, and Christopher J. Langmead 1,3 1 Computer Science Department, Carnegie
More informationYear 09 Science Learning Cycle 5 Overview
e Year 09 Science Learning Cycle 5 Overview Learning Cycle Overview: Biology How do we keep your body healthy L01 4.3.1.1 Communicable (infectious) disease L02 4.3.1.2 Viral diseases L03 4.3.1.3 Bacterial
More informationPrinciples of Synthetic Biology: Midterm Exam
Principles of Synthetic Biology: Midterm Exam October 28, 2010 1 Conceptual Simple Circuits 1.1 Consider the plots in figure 1. Identify all critical points with an x. Put a circle around the x for each
More information1.2. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy
.. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy = f(x, y). In this section we aim to understand the solution
More informationGaussian Conditional Random Fields for Modeling Patients Response to Acute Inflammation Treatment
Gaussian Conditional Random Fields for Modeling Patients Response to Acute Inflammation Treatment Vladan Radosavljevic a* vladan@temple.edu Kosta Ristovski a* kosta.ristovski@temple.edu Zoran Obradovic
More informationMath 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry
Math 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry Junping Shi College of William and Mary, USA Molecular biology and Biochemical kinetics Molecular biology is one of
More informationA Model-Free Control System Based on the Sliding Mode Control Method with Applications to Multi-Input-Multi-Output Systems
Proceedings of the 4 th International Conference of Control, Dynamic Systems, and Robotics (CDSR'17) Toronto, Canada August 21 23, 2017 Paper No. 119 DOI: 10.11159/cdsr17.119 A Model-Free Control System
More informationCHBE320 LECTURE XI CONTROLLER DESIGN AND PID CONTOLLER TUNING. Professor Dae Ryook Yang
CHBE320 LECTURE XI CONTROLLER DESIGN AND PID CONTOLLER TUNING Professor Dae Ryook Yang Spring 2018 Dept. of Chemical and Biological Engineering 11-1 Road Map of the Lecture XI Controller Design and PID
More informationOptimal State Estimation for Boolean Dynamical Systems using a Boolean Kalman Smoother
Optimal State Estimation for Boolean Dynamical Systems using a Boolean Kalman Smoother Mahdi Imani and Ulisses Braga-Neto Department of Electrical and Computer Engineering Texas A&M University College
More informationDay 3 Lecture 3. Optimizing deep networks
Day 3 Lecture 3 Optimizing deep networks Convex optimization A function is convex if for all α [0,1]: f(x) Tangent line Examples Quadratics 2-norms Properties Local minimum is global minimum x Gradient
More informationGeneral Model of the Innate Immune Response
General Model of the Innate Immune Response Katherine Reed, Kathryn Schalla, Souad Sosa, Jackie Tran, Thuy-My Truong, Alicia Prieto Langarica, Betty Scarbrough, Hristo Kojouharov, James Grover Technical
More informationGlobal Analysis of an Epidemic Model with Nonmonotone Incidence Rate
Global Analysis of an Epidemic Model with Nonmonotone Incidence Rate Dongmei Xiao Department of Mathematics, Shanghai Jiaotong University, Shanghai 00030, China E-mail: xiaodm@sjtu.edu.cn and Shigui Ruan
More informationNetworks in systems biology
Networks in systems biology Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Spring 2017 M. Macauley (Clemson) Networks in systems
More informationLocalization. Howie Choset Adapted from slides by Humphrey Hu, Trevor Decker, and Brad Neuman
Localization Howie Choset Adapted from slides by Humphrey Hu, Trevor Decker, and Brad Neuman Localization General robotic task Where am I? Techniques generalize to many estimation tasks System parameter
More informationLearning parameters in ODEs
Learning parameters in ODEs Application to biological networks Florence d Alché-Buc Joint work with Minh Quach and Nicolas Brunel IBISC FRE 3190 CNRS, Université d Évry-Val d Essonne, France /14 Florence
More informationTHE DYNAMICS OF ACUTE INFLAMMATION
THE DYNAMICS OF ACUTE INFLAMMATION by Rukmini Kumar M.Sc. Physics, Indian Institute of Technology, Madras, 1999 B.Sc. Physics, Stella Maris College, Madras, 1997 Submitted to the Graduate Faculty of the
More informationMATHEMATICAL MODELING OF CHEMICAL SIGNALS IN INFLAMMATORY PATHWAYS. Ian Price. B.A., University of Florida, Submitted to the Graduate Faculty of
MATHEMATICAL MODELING OF CHEMICAL SIGNALS IN INFLAMMATORY PATHWAYS by Ian Price B.A., University of Florida, 2004 Submitted to the Graduate Faculty of The Department of Mathematics in partial fulfillment
More informationAlgebraic Immunity of S-boxes and Augmented Functions
Algebraic Immunity of S-boxes and Augmented Functions Simon Fischer and Willi Meier S. Fischer and W. Meier AI of Sbox and AF 1 / 23 Outline 1 Algebraic Properties of S-boxes 2 Augmented Functions 3 Application
More informationPopulation Games and Evolutionary Dynamics
Population Games and Evolutionary Dynamics William H. Sandholm The MIT Press Cambridge, Massachusetts London, England in Brief Series Foreword Preface xvii xix 1 Introduction 1 1 Population Games 2 Population
More informationA MATHEMATICAL MODEL OF ACUTE INFLAMMATORY RESPONSE TO ENDOTOXIN CHALLENGE. Anirban Roy. Silvia Daun. Gilles Clermont. Jonathan Rubin.
MATHEMATICAL BIOSCIENCES AND ENGINEERING Volume xx, Number xx, xx 2xx doi:1.3934/mbe.29.xx.xx pp. 1 xx A MATHEMATICAL MODEL OF ACUTE INFLAMMATORY RESPONSE TO ENDOTOXIN CHALLENGE Anirban Roy Department
More information/639 Final Solutions, Part a) Equating the electrochemical potentials of H + and X on outside and inside: = RT ln H in
580.439/639 Final Solutions, 2014 Question 1 Part a) Equating the electrochemical potentials of H + and X on outside and inside: RT ln H out + zf 0 + RT ln X out = RT ln H in F 60 + RT ln X in 60 mv =
More informationHeart rate control and variability
Heart rate control and variability Na (Lina) Li (CDS13 ) EE @ SEAS Harvard University CDS @ 20 The persistent mystery Young, fit, healthy more extreme Resting Heart Rate (bpm) 60 0 50 100 150 200 250 300
More informationArtificial Intelligence
Artificial Intelligence Roman Barták Department of Theoretical Computer Science and Mathematical Logic Summary of last lecture We know how to do probabilistic reasoning over time transition model P(X t
More informationSupplemental Figures A Nonlinear Dynamical Theory of Cell Injury
upplementary Figures, A Nonlinear ynamical Theory of Cell Injury - upplemental Figures A Nonlinear ynamical Theory of Cell Injury onald J. egracia, Zhi-Feng Huang, ui Huang A. f ( X ) = a Θ n n Θ + X n
More informationSTUDENT VERSION DISEASE SPREAD
STUDENT VERSION DISEASE SPREAD STATEMENT We offer a modeling opportunity in which the phenomenon of the spread of disease can be described by one differential equation. First, let us consider growth in
More informationData assimilation with and without a model
Data assimilation with and without a model Tyrus Berry George Mason University NJIT Feb. 28, 2017 Postdoc supported by NSF This work is in collaboration with: Tim Sauer, GMU Franz Hamilton, Postdoc, NCSU
More information[Part 2] Model Development for the Prediction of Survival Times using Longitudinal Measurements
[Part 2] Model Development for the Prediction of Survival Times using Longitudinal Measurements Aasthaa Bansal PhD Pharmaceutical Outcomes Research & Policy Program University of Washington 69 Biomarkers
More informationStoichiometry Of Tumor Dynamics: Models And Analysis
Stoichiometry Of Tumor Dynamics: Models And Analysis Yang Kuang Department of Mathematics and Statistics Arizona State University Supported by NSF grant DMS-0077790 Y. Kuang, J. Nagy and J. Elser: Disc.
More informationNumerical approximation for optimal control problems via MPC and HJB. Giulia Fabrini
Numerical approximation for optimal control problems via MPC and HJB Giulia Fabrini Konstanz Women In Mathematics 15 May, 2018 G. Fabrini (University of Konstanz) Numerical approximation for OCP 1 / 33
More informationSTUDY OF THE DYNAMICAL MODEL OF HIV
STUDY OF THE DYNAMICAL MODEL OF HIV M.A. Lapshova, E.A. Shchepakina Samara National Research University, Samara, Russia Abstract. The paper is devoted to the study of the dynamical model of HIV. An application
More informationDynamical Systems and Chaos Part II: Biology Applications. Lecture 10: Coupled Systems. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part II: Biology Applications Lecture 10: Coupled Systems. Ilya Potapov Mathematics Department, TUT Room TD325 Foreword In order to model populations of physical/biological
More informationLS-N-IPS: AN IMPROVEMENT OF PARTICLE FILTERS BY MEANS OF LOCAL SEARCH
LS--IPS: A IMPROVEMET OF PARTICLE FILTERS BY MEAS OF LOCAL SEARCH Péter Torma Csaba Szepesvári Eötvös Loránd University, Budapest Mindmaker Ltd. H-1121 Budapest, Konkoly Thege Miklós út 29-33. Abstract:
More informationDERIVATIVE FREE OUTPUT FEEDBACK ADAPTIVE CONTROL
DERIVATIVE FREE OUTPUT FEEDBACK ADAPTIVE CONTROL Tansel YUCELEN, * Kilsoo KIM, and Anthony J. CALISE Georgia Institute of Technology, Yucelen Atlanta, * GA 30332, USA * tansel@gatech.edu AIAA Guidance,
More information3 Single species models. Reading: Otto & Day (2007) section 4.2.1
3 Single species models 3.1 Exponential growth Reading: Otto & Day (2007) section 4.2.1 We can solve equation 17 to find the population at time t given a starting population N(0) = N 0 as follows. N(t)
More informationContinuous-Time Dynamical Systems: Sketching solution curves and bifurcations
Math 360 Winter 07 Section 0 HW#4 Continuous-Time Dynamical Systems: Sketching solution curves and bifurcations (Due in one week! Thursday Nov 9, 07) Materials contained in this assignment are of great
More informationM.5 Modeling the Effect of Functional Responses
M.5 Modeling the Effect of Functional Responses The functional response is referred to the predation rate as a function of the number of prey per predator. It is recognized that as the number of prey increases,
More informationCONDENSATION Conditional Density Propagation for Visual Tracking
CONDENSATION Conditional Density Propagation for Visual Tracking Michael Isard and Andrew Blake Presented by Neil Alldrin Department of Computer Science & Engineering University of California, San Diego
More informationCONTROLLING WOUND HEALING THROUGH DEBRIDEMENT
CONTROLLING WOUND HEALING THROUGH DEBRIDEMENT MICHAEL A JONES AND DIANA M THOMAS Abstract The formation of slough (dead tissue) on a wound is widely accepted as an inhibitor to natural wound healing In
More informationCSC242: Intro to AI. Lecture 21
CSC242: Intro to AI Lecture 21 Administrivia Project 4 (homeworks 18 & 19) due Mon Apr 16 11:59PM Posters Apr 24 and 26 You need an idea! You need to present it nicely on 2-wide by 4-high landscape pages
More informationCONTROLLING CHAOS. Sudeshna Sinha. The Institute of Mathematical Sciences Chennai
CONTROLLING CHAOS Sudeshna Sinha The Institute of Mathematical Sciences Chennai Sinha, Ramswamy and Subba Rao: Physica D, vol. 43, p. 118 Sinha, Physics Letts. A vol. 156, p. 475 Ramswamy, Sinha, Gupte:
More informationLecture 7: Optimal Smoothing
Department of Biomedical Engineering and Computational Science Aalto University March 17, 2011 Contents 1 What is Optimal Smoothing? 2 Bayesian Optimal Smoothing Equations 3 Rauch-Tung-Striebel Smoother
More informationMulti-Level Fringe Rotation. JONATHAN D. ROMNEY National Radio Astronomy Observatory Charlottesville, Virginia
I VLB A Correlator Memo No. S3 (8705) Multi-Level Fringe Rotation JONATHAN D. ROMNEY National Radio Astronomy Observatory Charlottesville, Virginia 1987 February 3 Station-based fringe rotation using a
More informationNonlinear Model Predictive Control Tools (NMPC Tools)
Nonlinear Model Predictive Control Tools (NMPC Tools) Rishi Amrit, James B. Rawlings April 5, 2008 1 Formulation We consider a control system composed of three parts([2]). Estimator Target calculator Regulator
More informationPenalty and Barrier Methods. So we again build on our unconstrained algorithms, but in a different way.
AMSC 607 / CMSC 878o Advanced Numerical Optimization Fall 2008 UNIT 3: Constrained Optimization PART 3: Penalty and Barrier Methods Dianne P. O Leary c 2008 Reference: N&S Chapter 16 Penalty and Barrier
More informationHypothesis Testing, Power, Sample Size and Confidence Intervals (Part 2)
Hypothesis Testing, Power, Sample Size and Confidence Intervals (Part 2) B.H. Robbins Scholars Series June 23, 2010 1 / 29 Outline Z-test χ 2 -test Confidence Interval Sample size and power Relative effect
More informationP = 1 F m(p ) = IP = P I = f(i) = QI = IQ = 1 F m(p ) = Q, so we are done.
Section 1.6: Invertible Matrices One can show (exercise) that the composition of finitely many invertible functions is invertible. As a result, we have the following: Theorem 6.1: Any admissible row operation
More informationDetection of signal transitions by order statistics filtering
Detection of signal transitions by order statistics filtering A. Raji Images, Signals and Intelligent Systems Laboratory Paris-Est Creteil University, France Abstract In this article, we present a non
More informationGene Expression as a Stochastic Process: From Gene Number Distributions to Protein Statistics and Back
Gene Expression as a Stochastic Process: From Gene Number Distributions to Protein Statistics and Back June 19, 2007 Motivation & Basics A Stochastic Approach to Gene Expression Application to Experimental
More informationIntroduction to Compressed Sensing
Introduction to Compressed Sensing Alejandro Parada, Gonzalo Arce University of Delaware August 25, 2016 Motivation: Classical Sampling 1 Motivation: Classical Sampling Issues Some applications Radar Spectral
More informationState Estimation and Motion Tracking for Spatially Diverse VLC Networks
State Estimation and Motion Tracking for Spatially Diverse VLC Networks GLOBECOM Optical Wireless Communications Workshop December 3, 2012 Anaheim, CA Michael Rahaim mrahaim@bu.edu Gregary Prince gbprince@bu.edu
More informationEL 625 Lecture 10. Pole Placement and Observer Design. ẋ = Ax (1)
EL 625 Lecture 0 EL 625 Lecture 0 Pole Placement and Observer Design Pole Placement Consider the system ẋ Ax () The solution to this system is x(t) e At x(0) (2) If the eigenvalues of A all lie in the
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 3, Part 2: Introduction to Affine Parametrisation School of Electrical Engineering and Computer Science Lecture 3, Part 2: Affine Parametrisation p. 1/29 Outline
More informationComputational Intelligence Lecture 6: Associative Memory
Computational Intelligence Lecture 6: Associative Memory Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2011 Farzaneh Abdollahi Computational Intelligence
More informationEL2520 Control Theory and Practice
EL2520 Control Theory and Practice Lecture 8: Linear quadratic control Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden Linear quadratic control Allows to compute the controller
More informationSimple ODE Solvers - Derivation
Simple ODE Solvers - Derivation These notes provide derivations of some simple algorithms for generating, numerically, approximate solutions to the initial value problem y (t =f ( t, y(t y(t 0 =y 0 Here
More informationDynamic measurement: application of system identification in metrology
1 / 25 Dynamic measurement: application of system identification in metrology Ivan Markovsky Dynamic measurement takes into account the dynamical properties of the sensor 2 / 25 model of sensor as dynamical
More informationProportional, Integral & Derivative Control Design. Raktim Bhattacharya
AERO 422: Active Controls for Aerospace Vehicles Proportional, ntegral & Derivative Control Design Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University
More informationMath 2373: Linear Algebra and Differential Equations
Math 373: Linear Algebra and Differential Equations Paul Cazeaux* Fraser Hall 1, MW 15:35-1:5 September 7 th, 1 December th, 1 Contents 1 First-order differential equations 1 1.1 Dynamical systems and
More informationBiochemistry 3100 Sample Problems Binding proteins, Kinetics & Catalysis
(1) Draw an approximate denaturation curve for a typical blood protein (eg myoglobin) as a function of ph. (2) Myoglobin is a simple, single subunit binding protein that has an oxygen storage function
More informationVerification of analog and mixed-signal circuits using hybrid systems techniques
FMCAD, November 2004, Austin Verification of analog and mixed-signal circuits using hybrid systems techniques Thao Dang, Alexandre Donze, Oded Maler VERIMAG Grenoble, France Plan 1. Introduction 2. Verification
More informationCS-E5880 Modeling biological networks Gene regulatory networks
CS-E5880 Modeling biological networks Gene regulatory networks Jukka Intosalmi (based on slides by Harri Lähdesmäki) Department of Computer Science Aalto University January 12, 2018 Outline Modeling gene
More information2.152 Course Notes Contraction Analysis MIT, 2005
2.152 Course Notes Contraction Analysis MIT, 2005 Jean-Jacques Slotine Contraction Theory ẋ = f(x, t) If Θ(x, t) such that, uniformly x, t 0, F = ( Θ + Θ f x )Θ 1 < 0 Θ(x, t) T Θ(x, t) > 0 then all solutions
More informationOCR Biology Checklist
Topic 1. Cell level systems Video: Eukaryotic and prokaryotic cells Compare the structure of animal and plant cells. Label typical and atypical prokaryotic cells. Compare prokaryotic and eukaryotic cells.
More informationPopulation Games and Evolutionary Dynamics
Population Games and Evolutionary Dynamics (MIT Press, 200x; draft posted on my website) 1. Population games 2. Revision protocols and evolutionary dynamics 3. Potential games and their applications 4.
More informationOCR Biology Checklist
Topic 1. Cell level systems Video: Eukaryotic and prokaryotic cells Compare the structure of animal and plant cells. Label typical and atypical prokaryotic cells. Compare prokaryotic and eukaryotic cells.
More informationKalman Filtering in a Mass-Spring System
Department of Mathematics Center for Quantitative Sciences in Biomedicine North Carolina State University July 30, 2016 Practical Overview Goal: Implement Kalman filter for linear system 1. Solve a linear
More informationA new hybrid numerical scheme for modelling elastodynamics in unbounded media with near-source heterogeneities
A new hybrid numerical scheme for modelling elastodynamics in unbounded media with near-source heterogeneities Setare Hajarolasvadi Ahmed E. Elbanna https://www.youtube.com/watch?v=bltx92tuwha MOTIVATION
More informationCounterfactual Evaluation and Learning
SIGIR 26 Tutorial Counterfactual Evaluation and Learning Adith Swaminathan, Thorsten Joachims Department of Computer Science & Department of Information Science Cornell University Website: http://www.cs.cornell.edu/~adith/cfactsigir26/
More informationA nonparametric test for seasonal unit roots
Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna To be presented in Innsbruck November 7, 2007 Abstract We consider a nonparametric test for the
More informationNonlinear Regulation for Motorcycle Maneuvering
Nonlinear Regulation for Motorcycle Maneuvering John Hauser Univ of Colorado in collaboration with Alessandro Saccon* & Ruggero Frezza, Univ Padova * email asaccon@dei.unipd.it for dissertation aggressive
More informationSupporting Information. Methods. Equations for four regimes
Supporting Information A Methods All analytical expressions were obtained starting from quation 3, the tqssa approximation of the cycle, the derivation of which is discussed in Appendix C. The full mass
More informationClassification. Classification is similar to regression in that the goal is to use covariates to predict on outcome.
Classification Classification is similar to regression in that the goal is to use covariates to predict on outcome. We still have a vector of covariates X. However, the response is binary (or a few classes),
More informationA First Course on Kinetics and Reaction Engineering. Class 9 on Unit 9
A First Course on Kinetics and Reaction Engineering Class 9 on Unit 9 Part I - Chemical Reactions Part II - Chemical Reaction Kinetics Where We re Going A. Rate Expressions - 4. Reaction Rates and Temperature
More informationAdaptive modeling of biochemical pathways
Adaptive modeling of biochemical pathways R.Brause J.W.G.University, Frankfurt, Germany Brause@cs.uni-frankfurt.de Abstract In bioinformatics, biochemical pathways can be modeled by many differential equations.
More informationAdaptive Control of Variable-Speed Variable-Pitch Wind Turbines Using RBF Neural Network
Schulich School of Engineering Department of Mechanical and Manufacturing Engineering Adaptive Control of Variable-Speed Variable-Pitch Wind Turbines Using RBF Neural Network By: Hamidreza Jafarnejadsani,
More informationUsing Theorem Provers to Guarantee Closed-Loop Properties
Using Theorem Provers to Guarantee Closed-Loop Properties Nikos Aréchiga Sarah Loos André Platzer Bruce Krogh Carnegie Mellon University April 27, 2012 Aréchiga, Loos, Platzer, Krogh (CMU) Theorem Provers
More informationHarmonic Regression in the Biological Setting. Michael Gaffney, Ph.D., Pfizer Inc
Harmonic Regression in the Biological Setting Michael Gaffney, Ph.D., Pfizer Inc Two primary aims of harmonic regression 1. To describe the timing (phase) or degree of the diurnal variation (amplitude)
More informationIntroduction to stochastic multiscale modelling of tumour growth
Introduction to stochastic multiscale modelling of tumour growth Tomás Alarcón ICREA & Centre de Recerca Matemàtica T. Alarcón (ICREA & CRM, Barcelona, Spain) Lecture 1 CIMPA, Santiago de Cuba, June 2016
More informationApprentissage, réseaux de neurones et modèles graphiques (RCP209) Neural Networks and Deep Learning
Apprentissage, réseaux de neurones et modèles graphiques (RCP209) Neural Networks and Deep Learning Nicolas Thome Prenom.Nom@cnam.fr http://cedric.cnam.fr/vertigo/cours/ml2/ Département Informatique Conservatoire
More informationName Class. 5. Find the particular solution to given the general solution y C cos x and the. x 2 y
10 Differential Equations Test Form A 1. Find the general solution to the first order differential equation: y 1 yy 0. 1 (a) (b) ln y 1 y ln y 1 C y y C y 1 C y 1 y C. Find the general solution to the
More informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/319/5869/1543/dc1 Supporting Online Material for Synaptic Theory of Working Memory Gianluigi Mongillo, Omri Barak, Misha Tsodyks* *To whom correspondence should be addressed.
More informationTime depending inhibition with atypical kinetics, what is one to do? Ken Korzekwa Temple University School of Pharmacy
Time depending inhibition with atypical kinetics, what is one to do? Ken Korzekwa Temple University School of harmacy Acknowledgements TUS Swati Nagar Jaydeep Yadav harma - Donald Tweedie - Andrea Whitcher-Johnstone
More informationWhat is Experimental Design?
One Factor ANOVA What is Experimental Design? A designed experiment is a test in which purposeful changes are made to the input variables (x) so that we may observe and identify the reasons for change
More informationRELATIVE NAVIGATION FOR SATELLITES IN CLOSE PROXIMITY USING ANGLES-ONLY OBSERVATIONS
(Preprint) AAS 12-202 RELATIVE NAVIGATION FOR SATELLITES IN CLOSE PROXIMITY USING ANGLES-ONLY OBSERVATIONS Hemanshu Patel 1, T. Alan Lovell 2, Ryan Russell 3, Andrew Sinclair 4 "Relative navigation using
More informationStochastic contraction BACS Workshop Chamonix, January 14, 2008
Stochastic contraction BACS Workshop Chamonix, January 14, 2008 Q.-C. Pham N. Tabareau J.-J. Slotine Q.-C. Pham, N. Tabareau, J.-J. Slotine () Stochastic contraction 1 / 19 Why stochastic contraction?
More informationSynchronous state transition graph
Heike Siebert, FU Berlin, Molecular Networks WS10/11 2-1 Synchronous state transition graph (0, 2) (1, 2) vertex set X (state space) edges (x,f(x)) every state has only one successor attractors are fixed
More informationSurvival of Dominated Strategies under Evolutionary Dynamics. Josef Hofbauer University of Vienna. William H. Sandholm University of Wisconsin
Survival of Dominated Strategies under Evolutionary Dynamics Josef Hofbauer University of Vienna William H. Sandholm University of Wisconsin a hypnodisk 1 Survival of Dominated Strategies under Evolutionary
More informationUwe Aickelin and Qi Chen, School of Computer Science and IT, University of Nottingham, NG8 1BB, UK {uxa,
On Affinity Measures for Artificial Immune System Movie Recommenders Proceedings RASC-2004, The 5th International Conference on: Recent Advances in Soft Computing, Nottingham, UK, 2004. Uwe Aickelin and
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationA Multi-scale Extensive Petri Net Model of Bacterialmacrophage
A Multi-scale Extensive Petri Net Model of Bacterialmacrophage Interaction Rafael V. Carvalho Imaging & BioInformatics, Leiden Institute of Advanced Computer Science Introduction - Mycobacterial infection
More informationSMA 208: Ordinary differential equations I
SMA 208: Ordinary differential equations I Modeling with First Order differential equations Lecturer: Dr. Philip Ngare (Contacts: pngare@uonbi.ac.ke, Tue 12-2 PM) School of Mathematics, University of Nairobi
More informationOpen Loop Tuning Rules
Open Loop Tuning Rules Based on approximate process models Process Reaction Curve: The process reaction curve is an approximate model of the process, assuming the process behaves as a first order plus
More informationCHAPTER 6 FAST OUTPUT SAMPLING CONTROL TECHNIQUE
80 CHAPTER 6 FAST OUTPUT SAMPLING CONTROL TECHNIQUE 6.1 GENERAL In this chapter a control strategy for hyperthermia system is developed using fast output sampling feedback control law which is a type of
More informationAutonomous Equations and Stability Sections
A B I L E N E C H R I S T I A N U N I V E R S I T Y Department of Mathematics Autonomous Equations and Stability Sections 2.4-2.5 Dr. John Ehrke Department of Mathematics Fall 2012 Autonomous Differential
More information