5th-order differentiation

Size: px
Start display at page:

Download "5th-order differentiation"

Transcription

1 ARBITRARY-ORDER REAL-TIME EXACT ROBUST DIFFERENTIATION A. Levant Applied Mathematics Dept., Tel-Aviv University, Israel Homepage: 5th-order differentiation

2 Differentiation problematics: f ( t + τ) f ( t) Division by zero: f (t) = lim τ 0 τ Let f(t) = f 0 (t) + η(t), η(t) - noise f ( t + τ) f ( t) f η η = +, (-, ) τ τ τ τ 1st way: sin(ωt) = ω cos(ωt) Fourier transform, high harmonics neglection linear filter with transfer function Pm ( p) p p, m n. p Qn ( p) (0.01p Advantages: robustness Drawbacks: not exact, asymptotic convergence 2nd way: numeric differentiation Df(t) = F(f(t 1 ),..., f(t n )) + 1) 2 Advantages: exact on certain functions (for example polynomials P m (t)) Drawbacks: actual division by zero with t i t

3 Transfer Functions y (n) + a 1 y (n-1) a n y = b 0 u (m) + b 1 u (m-1) b m u (D n + a 1 D n a n )y = (b 0 D m + b 1 D m b m )u G(p) = b 0 p p n m + + b a 1 1 p p m 1 n b a n m u G( p) y u = e λt a y = G(λ) e λt in particular λ = iω, ω y = G(iω) with n < m Realization: Thus, only n m is feasible y (n) + a 1 y (n-1) a n y = u, G = p n + a 1 p 1 n a n Y = (b 0 p m + b 1 p m b m ) n 1 n 1 p + a p a 1 n U

4 Differentiation error = intrinsic differentiator error + noise differentiation error The goal: to build a differentiator with zero intrinsic error for a large class of functions and with error continuously dependent on the noise magnitude; + finite-time convergence The main trade-of: d D: f 1 a dt f1, D: f 2 a dt noise η = f 2 - f 1 d f2, D: f 1 + η a dt d f1 + dt d η no robust differentiator is exact on all polynomials or on all smooth functions

5 Differentiation Problem Input: f(t) = f 0 (t) + η(t), η < ε η(t) - Lebesgue-measurable function, f 0,η, ε are unknown, known : f (n+1) 0 (t) L (or Lipschitz constant of f (n) 0 LC ) The goal: real-time estimation of f & 0(t), & f 0(t),..., f (n) 0 (t) Theorem: Optimal differentiator existence γ ι (n) 1, D n : f( ) a D n i (f )( ), i = 0, 1,..., n sup f D i n (f )(t) - f (i) 0 (t) C = γ ι (n) L i/(n+1) (n-i+1) / (n+1) ε In particular γ 1 (1) = 2 (Landau, Kolmogorov) sup f D n n (f )(t) - f 0 (n) (t) L n/(n+1) ε 1/(n+1) ε = 10-6 f, f 0 (6) (t) 1: sup t D 5 5 (f )(t) - f 0 (5) (t) 0.1 That asymptotics cannot be improved

6 First-order differentiator idea Auxiliary dynamic system: Tracking problem: &z = u σ = z - f(t) 0 Any 2-sliding controller being applied, σ = σ& = 0 z - f(t) = u - f& (t) = 0 With input noises x - f(t) 0, u - f& (t) 0

7 A proof sketch Consider a smooth noise η(t) = ε sin (L/ε) 1/(n+1) t sup η (i) (t) = L i/(n+1) ε (n-i+1)/ (n+1), i = 0, 1,..., n +1 D differentiates both 0 and 0 + η(t) Lemma. βi(n) 1, β 0 (n) = β n (n +1) = 1: L i/(n+1) ε (n-i+1)/ (n+1) sup f ( i) ( t) ( ) L.const f n L, sup f ( t) ε β i (n) L i/(n+1) ε (n-i+1)/ (n+1). Some evaluations: β 1 (1) = 2 2, β 1 (2) = 7.07, β 2 (2) = W(C,n) = {f: [a, b] R, L.const of f (n) L} Let : g a g, g is closest to g in W(L,n) (not unique). Dg = ( g) (i).

8 Robust first order differentiator Input: f(t) = f 0 (t) + η(t), & f 0 < L η < ε, η(t) - measurable (Lebesgue) function z& 0 = v 0, v 0 = - λ 0 z 0 - f 1/2 sign(z 0 - f ) + z 1, z& 1= - λ 1 sign(z 0 - f ) (= - λ 1 sign(z 1 - v 0 )). Output: z 1 (t), v 0 (t) z 1-0 f& ~ L 1/2 ε 1/2 σ = z 0 - f(t) 0 For example: λ 1 = 1.1L, λ 0 =1.5L 1/2

9 Filippov Definition x& = v(x) x& V(x) x(t) is an absolutely continuous function I I V(x) = conv v( O ( x) \ N ) ε> 0µ N = 0 v +, v - - limit values, &x = pv + + (1-p)v -, p [0,1] Solutions exist for any locally bounded Lebesgue-measurable v(x); or for any upper-semicontinuous, convex, closed, non-empty, locally-bounded V(x) ε

10 1-differentiator explanation Let σ = z 0 - f 0 (t), ε = 0 σ& & = - & f λσ& σ -1/2 - α sign σ. σ& & λσ& σ -1/2 - [α - L, α + L] sign σ. invariant with respect to g ν : (σ, σ&, t) a (ν 2 σ, νσ&, νt)

11

12 Optimal n-differentiator z& 0 = v 0, v 0 = -λ 0 L 1/(n+1) z 0 - f(t) n/(n + 1) sign(z 0 -f(t)) + z 1, z& 1 = v 1, v 1 = -λ 1 L 1/n z 1 - v 0 (n-1)/ n sign(z 0 - v 0 ) + z 2,... z& i = v i, v i = -λ i L 1/(n-i+1) z i -v i-1 (n-i)/(n-i+1) sign(z i -v i-1 )+ z i+1,... z& n-1 = v n-1, v n-1 = -λ n-1 L 1/2 z n-1 -v n-2 1/ 2 sign(z n-1 -v n-2 )+ z n, z& n = -λ n L sign(z n - v n-1 ) D n : z& 0 = v, v = -µ z 0 - f(t) n/( n + 1) sign(z 0 - f(t)) + z 1, z 1 = 0 n 1 D (v( ), L),..., z n = n 1 n 1 D (v( ), L). D 0 : z& = -µ sign(z - f(t)), µ > L. Non-recursive form: z& 0 = -κ 0 z 0 - f(t) n/( n + 1) sign(z 0 - f(t)) + z 1,... z& i = -κ i z 0 - f(t) (n - i)/( n + 1) sign(z 0 - f(t)) + z i+1, i = 1,..., n-1... z& n = -κ n sign(z 0 - f(t))

13 5th-order differentiator, L = 1 z& 0 = v 0, v 0 = -12 z 0 - f(t) 5/6 sign(z 0 - f(t)) + z 1, z& 1 = v 1, v 1 = -8 z 1 - v 0 4/5 sign(z 1 - v 0 ) + z 2, z& 2 = v 2, v 2 = -5 z 2 - v 1 3/4 sign(z 2 - v 1 ) + z 3, z& 3 = v 3, v 3 = -3 z 3 - v 2 2/3 sign(z 3 - v 2 ) + z 4, z& 4 = v 4, v 4 = -1.5 z 4 - v 3 1/2 sign(z 4 - v 3 ) + z 5, z& 5 = -1.1 sign(z 5 - v 4 ); f(t ) = 0.5 sin 0.5t cos t

14 Theorem 1. f(t) = f 0 (t) in finite time z 0 = f 0 (t); z i = v i-1 = f 0 (i) (t), i = 1,..., n. 2-sliding modes: z i = f 0 (i) (t) i = 0,..., n-1. Theorem 2. f(t) - f 0 (t) ε z i - f 0 (i) (t) µ i ε (n - i +1)/(n + 1), i = 0,..., n; v i - f 0 (i+1) (t) ν i ε (n - i )/(n + 1), i = 0,..., n-1. z n - f (n) 1/(n + 1) 0 (t) µ n ε Discrete-sampling case: z 0 (t j ) - f(t j ) is substituted for z 0 - f(t) with t j t < t j+1, t j+1 - t j = τ > 0. Theorem 3. z i - f 0 (i) (t) µ i τ n - i + 1, i = 0,..., n; v i - f 0 (i+1) (t) ν i τ n - i, i = 0,..., n - 1. z n - f 0 (n) (t) µ n τ

15 Software restrictions of higher-order differentiation 5th-order differentiator τ = z 0 - f , z 1 - f& , z 2 - & f& , z 3 -& f& z 4 - f (4) , z 5 - f (5) τ = z 0 - f , z 1 - f& , z 2 - & f , z 3 -& f& z 4 - f (4) , z 5 - f (5)

16 Sensitivity to noises 3rd-order differentiator, ε = , ω 1000 D f& , D f& 0 & 0.14, 3 D3 -& f& th-order differentiator, ε = 0.01, ω 10 (the worst case) z 1 - f & , z 2 - & f , z 3 - & f& z 4 - f (4) , z 5 - f (5)

17 On-line calculation of the angular motor velocity and acceleration (data from Volvo Ltd) Experimental data, τ = nd order differentiation L = 625

18 On-line 2nd order differentiation comparison with optimal spline approximation

19 Image Processing: Crack Elimination

20 Edge Detection

Robust exact differentiation via sliding mode technique

Robust exact differentiation via sliding mode technique Robust exact differentiation via sliding mode technique Arie Levant * Abstract. The main problem in differentiator design is to combine differentiation exactness with robustness in respect to possible

More information

HIGHER ORDER SLIDING MODES AND ARBITRARY-ORDER EXACT ROBUST DIFFERENTIATION

HIGHER ORDER SLIDING MODES AND ARBITRARY-ORDER EXACT ROBUST DIFFERENTIATION HIGHER ORDER SLIDING MODES AND ARBITRARY-ORDER EXACT ROBUST DIFFERENTIATION A. Levant Institute for Industrial Mathematics, 4/24 Yehuda Ha-Nachtom St., Beer-Sheva 843, Israel Fax: +972-7-232 and E-mail:

More information

A. LEVANT School of Mathematical Sciences, Tel-Aviv University, Israel,

A. LEVANT School of Mathematical Sciences, Tel-Aviv University, Israel, Chapter 1 INTRODUCTION TO HIGH-ORDER SLIDING MODES A. LEVANT School of Mathematical Sciences, Tel-Aviv University, Israel, 2002-2003 1.1 Introduction One of the most important control problems is control

More information

An introduction to Mathematical Theory of Control

An introduction to Mathematical Theory of Control An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018

More information

C.-H. Lamarque. University of Lyon/ENTPE/LGCB & LTDS UMR CNRS 5513

C.-H. Lamarque. University of Lyon/ENTPE/LGCB & LTDS UMR CNRS 5513 Nonlinear Dynamics of Smooth and Non-Smooth Systems with Application to Passive Controls 3rd Sperlonga Summer School on Mechanics and Engineering Sciences on Dynamics, Stability and Control of Flexible

More information

HARMONIC ANALYSIS. Date:

HARMONIC ANALYSIS. Date: HARMONIC ANALYSIS Contents. Introduction 2. Hardy-Littlewood maximal function 3. Approximation by convolution 4. Muckenhaupt weights 4.. Calderón-Zygmund decomposition 5. Fourier transform 6. BMO (bounded

More information

The (Fast) Fourier Transform

The (Fast) Fourier Transform The (Fast) Fourier Transform The Fourier transform (FT) is the analog, for non-periodic functions, of the Fourier series for periodic functions can be considered as a Fourier series in the limit that the

More information

Computational Methods CMSC/AMSC/MAPL 460. Fourier transform

Computational Methods CMSC/AMSC/MAPL 460. Fourier transform Computational Methods CMSC/AMSC/MAPL 460 Fourier transform Ramani Duraiswami, Dept. of Computer Science Several slides from Prof. Healy s course at UMD Fourier Methods Fourier analysis ( harmonic analysis

More information

High-Gain Observers in Nonlinear Feedback Control. Lecture # 3 Regulation

High-Gain Observers in Nonlinear Feedback Control. Lecture # 3 Regulation High-Gain Observers in Nonlinear Feedback Control Lecture # 3 Regulation High-Gain ObserversinNonlinear Feedback ControlLecture # 3Regulation p. 1/5 Internal Model Principle d r Servo- Stabilizing u y

More information

Optimal stopping for non-linear expectations Part I

Optimal stopping for non-linear expectations Part I Stochastic Processes and their Applications 121 (2011) 185 211 www.elsevier.com/locate/spa Optimal stopping for non-linear expectations Part I Erhan Bayraktar, Song Yao Department of Mathematics, University

More information

Nonparametric regression with martingale increment errors

Nonparametric regression with martingale increment errors S. Gaïffas (LSTA - Paris 6) joint work with S. Delattre (LPMA - Paris 7) work in progress Motivations Some facts: Theoretical study of statistical algorithms requires stationary and ergodicity. Concentration

More information

State Space Representation of Gaussian Processes

State Space Representation of Gaussian Processes State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)

More information

Time Series Examples Sheet

Time Series Examples Sheet Lent Term 2001 Richard Weber Time Series Examples Sheet This is the examples sheet for the M. Phil. course in Time Series. A copy can be found at: http://www.statslab.cam.ac.uk/~rrw1/timeseries/ Throughout,

More information

Stability and Robustness of Homogeneous Differential Inclusions

Stability and Robustness of Homogeneous Differential Inclusions Stability and Robustness of Homogeneous Differential Inclusions Arie Levant 1, Denis Efimov 23, Andrey Polyakov 23, Wilfrid Perruquetti 2 Abstract The known results on asymptotic stability of homogeneous

More information

and in each case give the range of values of x for which the expansion is valid.

and in each case give the range of values of x for which the expansion is valid. α β γ δ ε ζ η θ ι κ λ µ ν ξ ο π ρ σ τ υ ϕ χ ψ ω Mathematics is indeed dangerous in that it absorbs students to such a degree that it dulls their senses to everything else P Kraft Further Maths A (MFPD)

More information

Measure and Integration: Solutions of CW2

Measure and Integration: Solutions of CW2 Measure and Integration: s of CW2 Fall 206 [G. Holzegel] December 9, 206 Problem of Sheet 5 a) Left (f n ) and (g n ) be sequences of integrable functions with f n (x) f (x) and g n (x) g (x) for almost

More information

ADAPTIVE OUTPUT FEEDBACK CONTROL OF NONLINEAR SYSTEMS YONGLIANG ZHU. Bachelor of Science Zhejiang University Hanzhou, Zhejiang, P.R.

ADAPTIVE OUTPUT FEEDBACK CONTROL OF NONLINEAR SYSTEMS YONGLIANG ZHU. Bachelor of Science Zhejiang University Hanzhou, Zhejiang, P.R. ADAPTIVE OUTPUT FEEDBACK CONTROL OF NONLINEAR SYSTEMS By YONGLIANG ZHU Bachelor of Science Zhejiang University Hanzhou, Zhejiang, P.R. China 1988 Master of Science Oklahoma State University Stillwater,

More information

University of California, Berkeley Department of Mechanical Engineering ME 104, Fall Midterm Exam 1 Solutions

University of California, Berkeley Department of Mechanical Engineering ME 104, Fall Midterm Exam 1 Solutions University of California, Berkeley Department of Mechanical Engineering ME 104, Fall 2013 Midterm Exam 1 Solutions 1. (20 points) (a) For a particle undergoing a rectilinear motion, the position, velocity,

More information

Phasor Young Won Lim 05/19/2015

Phasor Young Won Lim 05/19/2015 Phasor Copyright (c) 2009-2015 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version

More information

Overview of the Numerics of the ECMWF. Atmospheric Forecast Model

Overview of the Numerics of the ECMWF. Atmospheric Forecast Model Overview of the Numerics of the Atmospheric Forecast Model M. Hortal Seminar 6 Sept 2004 Slide 1 Characteristics of the model Hydrostatic shallow-atmosphere approimation Pressure-based hybrid vertical

More information

Periodic functions: simple harmonic oscillator

Periodic functions: simple harmonic oscillator Periodic functions: simple harmonic oscillator Recall the simple harmonic oscillator (e.g. mass-spring system) d 2 y dt 2 + ω2 0y = 0 Solution can be written in various ways: y(t) = Ae iω 0t y(t) = A cos

More information

Tribology in Industry. Some Aspects Concerning the Behaviour of Friction Materials at Low and Very Low Sliding Speeds

Tribology in Industry. Some Aspects Concerning the Behaviour of Friction Materials at Low and Very Low Sliding Speeds RESEARCH Vol. 37, No. 3 (2015) 374-379 Tribology in Industry www.tribology.fink.rs Some Aspects Concerning the Behaviour of Friction Materials at Low and Very Low Sliding Speeds N.A. Stoica a, A. Tudor

More information

Locally Lipschitzian Guiding Function Method for ODEs.

Locally Lipschitzian Guiding Function Method for ODEs. Locally Lipschitzian Guiding Function Method for ODEs. Marta Lewicka International School for Advanced Studies, SISSA, via Beirut 2-4, 3414 Trieste, Italy. E-mail: lewicka@sissa.it 1 Introduction Let f

More information

Input to state Stability

Input to state Stability Input to state Stability Mini course, Universität Stuttgart, November 2004 Lars Grüne, Mathematisches Institut, Universität Bayreuth Part III: Lyapunov functions and quantitative aspects ISS Consider with

More information

Computational Methods CMSC/AMSC/MAPL 460

Computational Methods CMSC/AMSC/MAPL 460 Computational Methods CMSC/AMSC/MAPL 460 Fourier transform Balaji Vasan Srinivasan Dept of Computer Science Several slides from Prof Healy s course at UMD Last time: Fourier analysis F(t) = A 0 /2 + A

More information

L6: Short-time Fourier analysis and synthesis

L6: Short-time Fourier analysis and synthesis L6: Short-time Fourier analysis and synthesis Overview Analysis: Fourier-transform view Analysis: filtering view Synthesis: filter bank summation (FBS) method Synthesis: overlap-add (OLA) method STFT magnitude

More information

A TWO PARAMETERS AMBROSETTI PRODI PROBLEM*

A TWO PARAMETERS AMBROSETTI PRODI PROBLEM* PORTUGALIAE MATHEMATICA Vol. 53 Fasc. 3 1996 A TWO PARAMETERS AMBROSETTI PRODI PROBLEM* C. De Coster** and P. Habets 1 Introduction The study of the Ambrosetti Prodi problem has started with the paper

More information

Trajectory Smoothing as a Linear Optimal Control Problem

Trajectory Smoothing as a Linear Optimal Control Problem Trajectory Smoothing as a Linear Optimal Control Problem Biswadip Dey & P. S. Krishnaprasad Allerton Conference - 212 Monticello, IL, USA. October 4, 212 Background and Motivation 1 Background and Motivation

More information

NONLINEAR DIFFERENTIAL INEQUALITY. 1. Introduction. In this paper the following nonlinear differential inequality

NONLINEAR DIFFERENTIAL INEQUALITY. 1. Introduction. In this paper the following nonlinear differential inequality M athematical Inequalities & Applications [2407] First Galley Proofs NONLINEAR DIFFERENTIAL INEQUALITY N. S. HOANG AND A. G. RAMM Abstract. A nonlinear differential inequality is formulated in the paper.

More information

Gramians based model reduction for hybrid switched systems

Gramians based model reduction for hybrid switched systems Gramians based model reduction for hybrid switched systems Y. Chahlaoui Younes.Chahlaoui@manchester.ac.uk Centre for Interdisciplinary Computational and Dynamical Analysis (CICADA) School of Mathematics

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department 8.323: Relativistic Quantum Field Theory I PROBLEM SET 2

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department 8.323: Relativistic Quantum Field Theory I PROBLEM SET 2 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department 8.323: Relativistic Quantum Field Theory I PROBLEM SET 2 REFERENCES: Peskin and Schroeder, Chapter 2 Problem 1: Complex scalar fields Peskin and

More information

Damped harmonic motion

Damped harmonic motion Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,

More information

RANDOM FIELDS AND GEOMETRY. Robert Adler and Jonathan Taylor

RANDOM FIELDS AND GEOMETRY. Robert Adler and Jonathan Taylor RANDOM FIELDS AND GEOMETRY from the book of the same name by Robert Adler and Jonathan Taylor IE&M, Technion, Israel, Statistics, Stanford, US. ie.technion.ac.il/adler.phtml www-stat.stanford.edu/ jtaylor

More information

Neighboring feasible trajectories in infinite dimension

Neighboring feasible trajectories in infinite dimension Neighboring feasible trajectories in infinite dimension Marco Mazzola Université Pierre et Marie Curie (Paris 6) H. Frankowska and E. M. Marchini Control of State Constrained Dynamical Systems Padova,

More information

Random Locations of Periodic Stationary Processes

Random Locations of Periodic Stationary Processes Random Locations of Periodic Stationary Processes Yi Shen Department of Statistics and Actuarial Science University of Waterloo Joint work with Jie Shen and Ruodu Wang May 3, 2016 The Fields Institute

More information

UNIVERSITY OF MANITOBA

UNIVERSITY OF MANITOBA Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic

More information

Radiation Damping. 1 Introduction to the Abraham-Lorentz equation

Radiation Damping. 1 Introduction to the Abraham-Lorentz equation Radiation Damping Lecture 18 1 Introduction to the Abraham-Lorentz equation Classically, a charged particle radiates energy if it is accelerated. We have previously obtained the Larmor expression for the

More information

Symmetries 2 - Rotations in Space

Symmetries 2 - Rotations in Space Symmetries 2 - Rotations in Space This symmetry is about the isotropy of space, i.e. space is the same in all orientations. Thus, if we continuously rotated an entire system in space, we expect the system

More information

MA3232 Summary 5. d y1 dy1. MATLAB has a number of built-in functions for solving stiff systems of ODEs. There are ode15s, ode23s, ode23t, ode23tb.

MA3232 Summary 5. d y1 dy1. MATLAB has a number of built-in functions for solving stiff systems of ODEs. There are ode15s, ode23s, ode23t, ode23tb. umerical solutions of higher order ODE We can convert a high order ODE into a system of first order ODEs and then apply RK method to solve it. Stiff ODEs Stiffness is a special problem that can arise in

More information

Fast-slow systems with chaotic noise

Fast-slow systems with chaotic noise Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 1, 216 Statistical properties of dynamical systems, ESI Vienna. David

More information

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points Nonlinear Control Lecture # 3 Stability of Equilibrium Points The Invariance Principle Definitions Let x(t) be a solution of ẋ = f(x) A point p is a positive limit point of x(t) if there is a sequence

More information

E2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2)

E2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2) E.5 Signals & Linear Systems Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & ) 1. Sketch each of the following continuous-time signals, specify if the signal is periodic/non-periodic,

More information

Joint work with Nguyen Hoang (Univ. Concepción, Chile) Padova, Italy, May 2018

Joint work with Nguyen Hoang (Univ. Concepción, Chile) Padova, Italy, May 2018 EXTENDED EULER-LAGRANGE AND HAMILTONIAN CONDITIONS IN OPTIMAL CONTROL OF SWEEPING PROCESSES WITH CONTROLLED MOVING SETS BORIS MORDUKHOVICH Wayne State University Talk given at the conference Optimization,

More information

Robust Adaptive Attitude Control of a Spacecraft

Robust Adaptive Attitude Control of a Spacecraft Robust Adaptive Attitude Control of a Spacecraft AER1503 Spacecraft Dynamics and Controls II April 24, 2015 Christopher Au Agenda Introduction Model Formulation Controller Designs Simulation Results 2

More information

Fourier transforms, Generalised functions and Greens functions

Fourier transforms, Generalised functions and Greens functions Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns

More information

Raktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Norms for Signals and Systems

Raktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Norms for Signals and Systems . AERO 632: Design of Advance Flight Control System Norms for Signals and. Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Norms for Signals ...

More information

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems Time-varying Systems ẋ = f(t,x) f(t,x) is piecewise continuous in t and locally Lipschitz in x for all t 0 and all x D, (0 D). The origin

More information

Lecture I: Asymptotics for large GUE random matrices

Lecture I: Asymptotics for large GUE random matrices Lecture I: Asymptotics for large GUE random matrices Steen Thorbjørnsen, University of Aarhus andom Matrices Definition. Let (Ω, F, P) be a probability space and let n be a positive integer. Then a random

More information

1 Lyapunov theory of stability

1 Lyapunov theory of stability M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability

More information

Letting p shows that {B t } t 0. Definition 0.5. For λ R let δ λ : A (V ) A (V ) be defined by. 1 = g (symmetric), and. 3. g

Letting p shows that {B t } t 0. Definition 0.5. For λ R let δ λ : A (V ) A (V ) be defined by. 1 = g (symmetric), and. 3. g 4 Contents.1 Lie group p variation results Suppose G, d) is a group equipped with a left invariant metric, i.e. Let a := d e, a), then d ca, cb) = d a, b) for all a, b, c G. d a, b) = d e, a 1 b ) = a

More information

2 Statement of the problem and assumptions

2 Statement of the problem and assumptions Mathematical Notes, 25, vol. 78, no. 4, pp. 466 48. Existence Theorem for Optimal Control Problems on an Infinite Time Interval A.V. Dmitruk and N.V. Kuz kina We consider an optimal control problem on

More information

Vehicle Dynamic Control Allocation for Path Following Moritz Gaiser

Vehicle Dynamic Control Allocation for Path Following Moritz Gaiser Vehicle Dynamic Control Allocation for Path Following Moritz Gaiser INSTITUT FÜR THEORETISCHE ELEKTROTECHNIK UND SYSTEMOPTIMIERUNG (ITE) KIT Die Forschungsuniversität in der Helmholtz-Gemeinschaft www.ite.kit.edu

More information

ON BOUNDEDNESS OF MAXIMAL FUNCTIONS IN SOBOLEV SPACES

ON BOUNDEDNESS OF MAXIMAL FUNCTIONS IN SOBOLEV SPACES Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 29, 2004, 167 176 ON BOUNDEDNESS OF MAXIMAL FUNCTIONS IN SOBOLEV SPACES Piotr Haj lasz and Jani Onninen Warsaw University, Institute of Mathematics

More information

Keble College - Hilary 2014 CP3&4: Mathematical methods I&II Tutorial 5 - Waves and normal modes II

Keble College - Hilary 2014 CP3&4: Mathematical methods I&II Tutorial 5 - Waves and normal modes II Tomi Johnson 1 Keble College - Hilary 2014 CP3&4: Mathematical methods I&II Tutorial 5 - Waves and normal modes II Prepare full solutions to the problems with a self assessment of your progress on a cover

More information

EE 230 Lecture 43. Data Converters

EE 230 Lecture 43. Data Converters EE 230 Lecture 43 Data Converters Review from Last Time: Amplitude Quantization Unwanted signals in the output of a system are called noise. Distortion Smooth nonlinearities Frequency attenuation Large

More information

Functional Differential Equations with Causal Operators

Functional Differential Equations with Causal Operators ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.11(211) No.4,pp.499-55 Functional Differential Equations with Causal Operators Vasile Lupulescu Constantin Brancusi

More information

Filtering and Edge Detection

Filtering and Edge Detection Filtering and Edge Detection Local Neighborhoods Hard to tell anything from a single pixel Example: you see a reddish pixel. Is this the object s color? Illumination? Noise? The next step in order of complexity

More information

Time Series Examples Sheet

Time Series Examples Sheet Lent Term 2001 Richard Weber Time Series Examples Sheet This is the examples sheet for the M. Phil. course in Time Series. A copy can be found at: http://www.statslab.cam.ac.uk/~rrw1/timeseries/ Throughout,

More information

Parallel Recursion: Powerlists. Greg Plaxton Theory in Programming Practice, Fall 2005 Department of Computer Science University of Texas at Austin

Parallel Recursion: Powerlists. Greg Plaxton Theory in Programming Practice, Fall 2005 Department of Computer Science University of Texas at Austin Parallel Recursion: Powerlists Greg Plaxton Theory in Programming Practice, Fall 2005 Department of Computer Science University of Texas at Austin Overview Definition of a powerlist Basic operations on

More information

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems Time-varying Systems ẋ = f(t,x) f(t,x) is piecewise continuous in t and locally Lipschitz in x for all t 0 and all x D, (0 D). The origin

More information

L -Bounded Robust Control of Nonlinear Cascade Systems

L -Bounded Robust Control of Nonlinear Cascade Systems L -Bounded Robust Control of Nonlinear Cascade Systems Shoudong Huang M.R. James Z.P. Jiang August 19, 2004 Accepted by Systems & Control Letters Abstract In this paper, we consider the L -bounded robust

More information

An introduction to Birkhoff normal form

An introduction to Birkhoff normal form An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an

More information

Linear Filters. L[e iωt ] = 2π ĥ(ω)eiωt. Proof: Let L[e iωt ] = ẽ ω (t). Because L is time-invariant, we have that. L[e iω(t a) ] = ẽ ω (t a).

Linear Filters. L[e iωt ] = 2π ĥ(ω)eiωt. Proof: Let L[e iωt ] = ẽ ω (t). Because L is time-invariant, we have that. L[e iω(t a) ] = ẽ ω (t a). Linear Filters 1. Convolutions and filters. A filter is a black box that takes an input signal, processes it, and then returns an output signal that in some way modifies the input. For example, if the

More information

Applied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.

Applied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books. Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define

More information

be the set of complex valued 2π-periodic functions f on R such that

be the set of complex valued 2π-periodic functions f on R such that . Fourier series. Definition.. Given a real number P, we say a complex valued function f on R is P -periodic if f(x + P ) f(x) for all x R. We let be the set of complex valued -periodic functions f on

More information

Homogeneous High-Order Sliding Modes

Homogeneous High-Order Sliding Modes Homogeneous High-Order Sliding Modes A. Levant* * Applied Mathematics Department, Tel-Aviv University, Tel-Aviv, Israel, (e-mail: levant@post.tau.ac.il) Abstract: Homogeneity features of dynamic systems

More information

Regularity of flat level sets in phase transitions

Regularity of flat level sets in phase transitions Annals of Mathematics, 69 (2009), 4 78 Regularity of flat level sets in phase transitions By Ovidiu Savin Abstract We consider local minimizers of the Ginzburg-Landau energy functional 2 u 2 + 4 ( u2 )

More information

PHASE TRANSITIONS: REGULARITY OF FLAT LEVEL SETS

PHASE TRANSITIONS: REGULARITY OF FLAT LEVEL SETS PHASE TRANSITIONS: REGULARITY OF FLAT LEVEL SETS OVIDIU SAVIN Abstract. We consider local minimizers of the Ginzburg-Landau energy functional 2 u 2 + 4 ( u2 ) 2 dx and prove that, if the level set is included

More information

Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems

Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t

More information

Robot Manipulator Control. Hesheng Wang Dept. of Automation

Robot Manipulator Control. Hesheng Wang Dept. of Automation Robot Manipulator Control Hesheng Wang Dept. of Automation Introduction Industrial robots work based on the teaching/playback scheme Operators teach the task procedure to a robot he robot plays back eecute

More information

5 Applying the Fokker-Planck equation

5 Applying the Fokker-Planck equation 5 Applying the Fokker-Planck equation We begin with one-dimensional examples, keeping g = constant. Recall: the FPE for the Langevin equation with η(t 1 )η(t ) = κδ(t 1 t ) is = f(x) + g(x)η(t) t = x [f(x)p

More information

PHYSICS 110A : CLASSICAL MECHANICS HW 2 SOLUTIONS. Here is a sketch of the potential with A = 1, R = 1, and S = 1. From the plot we can see

PHYSICS 110A : CLASSICAL MECHANICS HW 2 SOLUTIONS. Here is a sketch of the potential with A = 1, R = 1, and S = 1. From the plot we can see PHYSICS 11A : CLASSICAL MECHANICS HW SOLUTIONS (1) Taylor 5. Here is a sketch of the potential with A = 1, R = 1, and S = 1. From the plot we can see 1.5 1 U(r).5.5 1 4 6 8 1 r Figure 1: Plot for problem

More information

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014 Prof. Dr. Eleni Chatzi System Stability - 26 March, 24 Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can

More information

Chaos suppression of uncertain gyros in a given finite time

Chaos suppression of uncertain gyros in a given finite time Chin. Phys. B Vol. 1, No. 11 1 1155 Chaos suppression of uncertain gyros in a given finite time Mohammad Pourmahmood Aghababa a and Hasan Pourmahmood Aghababa bc a Electrical Engineering Department, Urmia

More information

Ring-like structures of frequency domains of wavelets

Ring-like structures of frequency domains of wavelets Ring-like structures of frequency domains of wavelets Zhihua Zhang and Naoki aito Dept. of Math., Univ. of California, Davis, California, 95616, UA. E-mail: zzh@ucdavis.edu saito@math.ucdavis.edu Abstract.

More information

Regularity estimates for fully non linear elliptic equations which are asymptotically convex

Regularity estimates for fully non linear elliptic equations which are asymptotically convex Regularity estimates for fully non linear elliptic equations which are asymptotically convex Luis Silvestre and Eduardo V. Teixeira Abstract In this paper we deliver improved C 1,α regularity estimates

More information

3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?

3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure? MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due ). Show that the open disk x 2 + y 2 < 1 is a countable union of planar elementary sets. Show that the closed disk x 2 + y 2 1 is a countable

More information

Principle of Mathematical Induction

Principle of Mathematical Induction Advanced Calculus I. Math 451, Fall 2016, Prof. Vershynin Principle of Mathematical Induction 1. Prove that 1 + 2 + + n = 1 n(n + 1) for all n N. 2 2. Prove that 1 2 + 2 2 + + n 2 = 1 n(n + 1)(2n + 1)

More information

Control System. Contents

Control System. Contents Contents Chapter Topic Page Chapter- Chapter- Chapter-3 Chapter-4 Introduction Transfer Function, Block Diagrams and Signal Flow Graphs Mathematical Modeling Control System 35 Time Response Analysis of

More information

Rice method for the maximum of Gaussian fields

Rice method for the maximum of Gaussian fields Rice method for the maximum of Gaussian fields IWAP, July 8, 2010 Jean-Marc AZAÏS Institut de Mathématiques, Université de Toulouse Jean-Marc AZAÏS ( Institut de Mathématiques, Université de Toulouse )

More information

Module 25: Outline Resonance & Resonance Driven & LRC Circuits Circuits 2

Module 25: Outline Resonance & Resonance Driven & LRC Circuits Circuits 2 Module 25: Driven RLC Circuits 1 Module 25: Outline Resonance & Driven LRC Circuits 2 Driven Oscillations: Resonance 3 Mass on a Spring: Simple Harmonic Motion A Second Look 4 Mass on a Spring (1) (2)

More information

Alternatives. The D Operator

Alternatives. The D Operator Using Smoothness Alternatives Text: Chapter 5 Some disadvantages of basis expansions Discrete choice of number of basis functions additional variability. Non-hierarchical bases (eg B-splines) make life

More information

Bridging the Gap between Center and Tail for Multiscale Processes

Bridging the Gap between Center and Tail for Multiscale Processes Bridging the Gap between Center and Tail for Multiscale Processes Matthew R. Morse Department of Mathematics and Statistics Boston University BU-Keio 2016, August 16 Matthew R. Morse (BU) Moderate Deviations

More information

Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science : Discrete-Time Signal Processing

Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science : Discrete-Time Signal Processing Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.34: Discrete-Time Signal Processing OpenCourseWare 006 ecture 8 Periodogram Reading: Sections 0.6 and 0.7

More information

Calculating the domain of attraction: Zubov s method and extensions

Calculating the domain of attraction: Zubov s method and extensions Calculating the domain of attraction: Zubov s method and extensions Fabio Camilli 1 Lars Grüne 2 Fabian Wirth 3 1 University of L Aquila, Italy 2 University of Bayreuth, Germany 3 Hamilton Institute, NUI

More information

Problem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension

Problem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension 105A Practice Final Solutions March 13, 01 William Kelly Problem 1: Lagrangians and Conserved Quantities Consider the following action for a particle of mass m moving in one dimension S = dtl = mc dt 1

More information

Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches

Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Noise is often considered as some disturbing component of the system. In particular physical situations, noise becomes

More information

In terms of measures: Exercise 1. Existence of a Gaussian process: Theorem 2. Remark 3.

In terms of measures: Exercise 1. Existence of a Gaussian process: Theorem 2. Remark 3. 1. GAUSSIAN PROCESSES A Gaussian process on a set T is a collection of random variables X =(X t ) t T on a common probability space such that for any n 1 and any t 1,...,t n T, the vector (X(t 1 ),...,X(t

More information

Large Deviations, Linear Statistics, and Scaling Limits for Mahler Ensemble of Complex Random Polynomials

Large Deviations, Linear Statistics, and Scaling Limits for Mahler Ensemble of Complex Random Polynomials Large Deviations, Linear Statistics, and Scaling Limits for Mahler Ensemble of Complex Random Polynomials Maxim L. Yattselev joint work with Christopher D. Sinclair International Conference on Approximation

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

Convex Optimization. Ofer Meshi. Lecture 6: Lower Bounds Constrained Optimization

Convex Optimization. Ofer Meshi. Lecture 6: Lower Bounds Constrained Optimization Convex Optimization Ofer Meshi Lecture 6: Lower Bounds Constrained Optimization Lower Bounds Some upper bounds: #iter μ 2 M #iter 2 M #iter L L μ 2 Oracle/ops GD κ log 1/ε M x # ε L # x # L # ε # με f

More information

On reduction of differential inclusions and Lyapunov stability

On reduction of differential inclusions and Lyapunov stability 1 On reduction of differential inclusions and Lyapunov stability Rushikesh Kamalapurkar, Warren E. Dixon, and Andrew R. Teel arxiv:1703.07071v5 [cs.sy] 25 Oct 2018 Abstract In this paper, locally Lipschitz

More information

Laboratory handout 5 Mode shapes and resonance

Laboratory handout 5 Mode shapes and resonance laboratory handouts, me 34 82 Laboratory handout 5 Mode shapes and resonance In this handout, material and assignments marked as optional can be skipped when preparing for the lab, but may provide a useful

More information

Complex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,

Complex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )

More information

Spaces with Ricci curvature bounded from below

Spaces with Ricci curvature bounded from below Spaces with Ricci curvature bounded from below Nicola Gigli February 23, 2015 Topics 1) On the definition of spaces with Ricci curvature bounded from below 2) Analytic properties of RCD(K, N) spaces 3)

More information

LMI Methods in Optimal and Robust Control

LMI Methods in Optimal and Robust Control LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear

More information

Oscillations and Waves

Oscillations and Waves Oscillations and Waves Somnath Bharadwaj and S. Pratik Khastgir Department of Physics and Meteorology IIT Kharagpur Module : Oscillations Lecture : Oscillations Oscillations are ubiquitous. It would be

More information

1 Functions of normal random variables

1 Functions of normal random variables Tel Aviv University, 200 Gaussian measures : results formulated Functions of normal random variables It is of course impossible to even think the word Gaussian without immediately mentioning the most important

More information

Learning to play partially-specified equilibrium

Learning to play partially-specified equilibrium Learning to play partially-specified equilibrium Ehud Lehrer and Eilon Solan August 29, 2007 August 29, 2007 First draft: June 2007 Abstract: In a partially-specified correlated equilibrium (PSCE ) the

More information

6. Approximation and fitting

6. Approximation and fitting 6. Approximation and fitting Convex Optimization Boyd & Vandenberghe norm approximation least-norm problems regularized approximation robust approximation 6 Norm approximation minimize Ax b (A R m n with

More information