Last time: Completed solution to the optimum linear filter in real-time operation

Similar documents
Last time: Ground rules for filtering and control system design

EECE 301 Signals & Systems Prof. Mark Fowler

u t u 0 ( 7) Intuitively, the maximum principles can be explained by the following observation. Recall

State space systems analysis

ELEC 372 LECTURE NOTES, WEEK 4 Dr. Amir G. Aghdam Concordia University

ECEN620: Network Theory Broadband Circuit Design Fall 2014

Heat Equation: Maximum Principles

Lecture 30: Frequency Response of Second-Order Systems

Applied Mathematical Sciences, Vol. 9, 2015, no. 3, HIKARI Ltd,

CONTROL ENGINEERING LABORATORY

High-Speed Serial Interface Circuits and Systems. Lect. 4 Phase-Locked Loop (PLL) Type 1 (Chap. 8 in Razavi)

THE CONCEPT OF THE ROOT LOCUS. H(s) THE CONCEPT OF THE ROOT LOCUS

ELEC 372 LECTURE NOTES, WEEK 1 Dr. Amir G. Aghdam Concordia University

LECTURE 13 SIMULTANEOUS EQUATIONS

The Performance of Feedback Control Systems

Lecture 11. Course Review. (The Big Picture) G. Hovland Input-Output Limitations (Skogestad Ch. 3) Discrete. Time Domain

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

We will look for series solutions to (1) around (at most) regular singular points, which without

Introduction to Control Systems

Comments on Discussion Sheet 18 and Worksheet 18 ( ) An Introduction to Hypothesis Testing

Erick L. Oberstar Fall 2001 Project: Sidelobe Canceller & GSC 1. Advanced Digital Signal Processing Sidelobe Canceller (Beam Former)

10-716: Advanced Machine Learning Spring Lecture 13: March 5

Automatic Control Systems

13.4 Scalar Kalman Filter

DISCRETE MELLIN CONVOLUTION AND ITS EXTENSIONS, PERRON FORMULA AND EXPLICIT FORMULAE

Lesson 5. Chapter 7. Wiener Filters. Bengt Mandersson. r x k We assume uncorrelated noise v(n). LTH. September 2010

STA 4032 Final Exam Formula Sheet

ME 410 MECHANICAL ENGINEERING SYSTEMS LABORATORY REGRESSION ANALYSIS

x z Increasing the size of the sample increases the power (reduces the probability of a Type II error) when the significance level remains fixed.

a 1 = 1 a a a a n n s f() s = Σ log a 1 + a a n log n sup log a n+1 + a n+2 + a n+3 log n sup () s = an /n s s = + t i

Brief Review of Linear System Theory

Explicit scheme. Fully implicit scheme Notes. Fully implicit scheme Notes. Fully implicit scheme Notes. Notes

STRONG DEVIATION THEOREMS FOR THE SEQUENCE OF CONTINUOUS RANDOM VARIABLES AND THE APPROACH OF LAPLACE TRANSFORM

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 16 11/04/2013. Ito integral. Properties

Chem Exam 1-9/14/16. Frequency. Grade Average = 72, Median = 72, s = 20

Fig. 1: Streamline coordinates

SOLUTION: The 95% confidence interval for the population mean µ is x ± t 0.025; 49

Statistics Problem Set - modified July 25, _. d Q w. i n

Statistics and Chemical Measurements: Quantifying Uncertainty. Normal or Gaussian Distribution The Bell Curve

Analysis of Analytical and Numerical Methods of Epidemic Models

8.6 Order-Recursive LS s[n]

2015 Fermat Contest (Grade 11)

ON THE SCALE PARAMETER OF EXPONENTIAL DISTRIBUTION

VIII. Interval Estimation A. A Few Important Definitions (Including Some Reminders)

Chapter 4 : Laplace Transform

Chapter 9. Key Ideas Hypothesis Test (Two Populations)

Lesson 5. Chapter 7. Wiener Filters. Bengt Mandersson. r k s r x LTH. September Prediction Error Filter PEF (second order) from chapter 4

Chapter 1 Econometrics

1it is said to be overdamped. When 1, the roots of

ECM Control Engineering Dr Mustafa M Aziz (2013) SYSTEM RESPONSE

TUTORIAL 6. Review of Electrostatic

REVIEW OF SIMPLE LINEAR REGRESSION SIMPLE LINEAR REGRESSION

Chapter #5 EEE Control Systems

Control of a Linear Permanent Magnet Synchronous Motor using Multiple Reference Frame Theory

20. CONFIDENCE INTERVALS FOR THE MEAN, UNKNOWN VARIANCE

M227 Chapter 9 Section 1 Testing Two Parameters: Means, Variances, Proportions

Time Response. First Order Systems. Time Constant, T c We call 1/a the time constant of the response. Chapter 4 Time Response

Performance-Based Plastic Design (PBPD) Procedure

STABILITY OF THE ACTIVE VIBRATION CONTROL OF CANTILEVER BEAMS

Generalized Likelihood Functions and Random Measures

Another Look at Estimation for MA(1) Processes With a Unit Root

1. (a) If u (I : R J), there exists c 0 in R such that for all q 0, cu q (I : R J) [q] I [q] : R J [q]. Hence, if j J, for all q 0, j q (cu q ) =

Gain-scheduling of Acceleration Estimator for Low-velocity Measurement with Encoders

CONTROL SYSTEMS. Chapter 7 : Bode Plot. 40dB/dec 1.0. db/dec so resultant slope will be 20 db/dec and this is due to the factor s

Professor: Mihnea UDREA DIGITAL SIGNAL PROCESSING. Grading: Web: MOODLE. 1. Introduction. General information

Chapter 8.2. Interval Estimation

Hidden Markov Model Parameters

Digital Control System

Zeta-reciprocal Extended reciprocal zeta function and an alternate formulation of the Riemann hypothesis By M. Aslam Chaudhry

Dorf, R.C., Wan, Z., Johnson, D.E. Laplace Transform The Electrical Engineering Handbook Ed. Richard C. Dorf Boca Raton: CRC Press LLC, 2000

System Control. Lesson #19a. BME 333 Biomedical Signals and Systems - J.Schesser

TESTS OF SIGNIFICANCE

ME2142/ME2142E Feedback Control Systems

September 2012 C1 Note. C1 Notes (Edexcel) Copyright - For AS, A2 notes and IGCSE / GCSE worksheets 1

Assignment 1 - Solutions. ECSE 420 Parallel Computing Fall November 2, 2014

A Tool for the Spectral Analysis of the Laser Doppler Anemometer Data of the Cambridge Stratified Swirl Burner

Production Scheduling with Genetic Algorithm

Reliable Decentralized PID Stabilization of MIMO Systems

ECE 422 Power System Operations & Planning 6 Small Signal Stability. Spring 2015 Instructor: Kai Sun

A Simplified Derivation of Scalar Kalman Filter using Bayesian Probability Theory

Statistical Inference Procedures

Capacitors and PN Junctions. Lecture 8: Prof. Niknejad. Department of EECS University of California, Berkeley. EECS 105 Fall 2003, Lecture 8

STUDENT S t-distribution AND CONFIDENCE INTERVALS OF THE MEAN ( )

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

Cryptographic Hash Functions and Message Authentication Codes. Reading: Chapter 4 of Katz & Lindell

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

Adaptive control design for a Mimo chemical reactor

1 Routh Array: 15 points

Société de Calcul Mathématique, S. A. Algorithmes et Optimisation

Lecture 8: Period Finding: Simon s Problem over Z N

Oktober Fórum 2005 PPGEQ

18.05 Problem Set 9, Spring 2014 Solutions

THE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.

Abstract Vector Spaces. Abstract Vector Spaces

A Faster Product for π and a New Integral for ln π 2

Further Investigation of alternative Formulation of RP Model with Response Error. Ed Stanek

Topics in MMSE Estimation for Sparse Approximation

arxiv: v1 [math.pr] 19 Aug 2015

M2.The Z-Transform and its Properties

Weak formulation and Lagrange equations of motion

Transcription:

6.3 tochatic Etimatio ad Cotrol, Fall 4 ecture at time: Completed olutio to the oimum liear filter i real-time operatio emi-free cofiguratio: t D( p) F( p) i( p) dte dp e π F( ) F( ) ( ) F( p) ( p) 4444443 pecial cae: ( p) ( p) i ratioal: I thi olutio formula we ca carry out the idicated itegratio i literal form i the cae i which ( p) i ratioal. I our work, we deal i a practical way oly with ratioal F, i, ad, o thi fuctio will e ratioal if D( p ) i ratioal. Thi will e true of every deired operatio exce a predictor. Thu exce i the cae of predictio, the aove fuctio which will e ymolized a [ ] ca e expaded ito ( p) ( p) + ( p) where [ ] ha pole oly i HP ad [ ] ha pole oly i HP. The zeroe may e aywhere. For ratioal [ ], thi expaio i made y expadig ito partial fractio, the addig together the term defg HP pole to form [ ] ad addig together the term defg HP pole to form [ ]. ctually, oly [ ] will e required. t t t dte dp{ ( p) + ( p) } e f() t e dt + f() t e dt where f() t ( p) e dp, t < f() t ( p) e dp, t > Note that f() t i the ivere traform of a fuctio which i aalytic i HP; thu f() t for t > ad Page of 7

6.3 tochatic Etimatio ad Cotrol, Fall 4 t f () te dt lo f() t i the ivere traform of a fuctio which i aalytic i HP; thu f () t for t <. Thu ( ) t t f () te dt f() te dt Thu fially, DF () ( ) i() F () () F () F( ) () I the uual cae, Fi () a tale, mmum phae fuctio. I that cae, F () F (), F () ; that i, all the pole ad zeroe of F () are i the HP. imilarly, F( ). The D () i() F () () Thu i thi cae the oimum trafer fuctio from iput to output i D () i() FH () () ad the oimum fuctio to e cacaded with the fixed part i otaied from thi y divo y F, () o that the fixed part i compeated out y cacellatio. Free cofiguratio prolem: D () i() Oimum free cofiguratio filter: Page of 7

6.3 tochatic Etimatio ad Cotrol, Fall 4 We tarted with a cloed-loop cofiguratio: C () H() + CFB ( ) ( ) ( ) H() C () H( ) F( ) B( ) The loop will e tale, ut C () may e utale. pecial commet aout the applicatio of thee formulae: a) Utale F () caot e treated ecaue the Fourier traform of wf () t doe ot coverge i that cae. To treat thi ytem, firt cloe a feedack loop aroud F () to create a tale fixed part ad work with thi tale feedack ytem a F. () Whe the oimum compeatio i foud, it ca e collected with the origial compeatio if deired. ) F () which ha pole o the ω axi i the limitig cae of fuctio for which the Fourier traform coverge. You ca move the pole ut ito the HP y addig a real part + ε to the pole locatio. olve the prolem with thi ε ad at the ed et it to zero. Zeroe of F () o ω axi ca e icluded i either factor ad the reult will e the ame. Thi will permit cacellatio compeatio of pole of F () o the ω axi, icludig pole at the origi. c) I factorig ito, ay cotat factor i ca e divided etwee ad i ay coveiet way. The ame i true of F () ad F( ). d) Prolem hould e well-poed i the firt place. void comiatio of D () ad ( ω ) which imply ifte dt () ecaue that may aume ifte e for ay realizale filter. uch a a differetiator o a igal which fall off a ω. e) The poit at t wa left hagig i everal tep of the derivatio of the olutio formula. Do t other checkig the idividual tep; ut check the fial olutio to ee if it atifie the eceary coditio. Page 3 of 7

6.3 tochatic Etimatio ad Cotrol, Fall 4 The Wieer-Hopf equatio require l( τ ) for τ. Thu () hould e aalytic i HP ad go to zero at leat a fat a for large. () F( ) H () F() () F( ) D() () i We have olved the prolem of the oimum liear filter uder the leat mea quared error criterio. Further aalyi how that if the iput, igal ad oie, are Gauia, the reult we have i the oimum filter. Thi i, there i o filter, liear or oliear which will yield maller mea quared error. If the iput are ot oth Gauia, it i almot ure that ome oliear filter ca do etter tha the Wieer filter. But theory for thi i oly egg to e developed o a approximate ai. Note that if we oly kow the ecod order tatitic of the iput, the oimum liear filter i the et we ca do. To take advatage of oliear filterig we mut kow the ditriutio of the iput. Example: Free cofiguratio predictor (real time) () a () The, are ucorrelated. D () e T Page 4 of 7

6.3 tochatic Etimatio ad Cotrol, Fall 4 Ue the olutio form t D( p) i( p) dte dp e π ( ) ( p) t i p p( t+ T ) ( ) dte dp e π ( ) ( p) () () + () + a + + a a a + a + a where a +. i( p) ( p) ( a+ )( a ) i( p) ( a ) ( p) ( a+ )( a )( ) ( a+ )( ) + a + a+ a+ Uig the itegral form, a( t T) ( T) e + a+, t T e + > dp a + p+ a, otherwie Page 5 of 7

6.3 tochatic Etimatio ad Cotrol, Fall 4 t ( T) ( T) e + a+, t T e + < dp a + p, otherwie t a( t+ T ) at ( + a) t e e dt e e dt a+ a+ at e a + + a at at ( + a) e e a+ ( + a) ( + ) ( a+ ) + where a +. Note the adwidth of thi filter ad the gai ~ etwee t () ad t ( + T) e at which i the correlatio Example: emi-free prolem with o-mmum phae F. Oimum compeator ( c ) ( + ) F () K d () a () The, are ucorrelated. ervo example where we d like the output to track the iput, o the deired operator, D (). Page 6 of 7

6.3 tochatic Etimatio ad Cotrol, Fall 4 a + a + a K F () ( + ε )( d + ) F () c c+ F( ) K ε d F( ) c+ ( )( ) ( ) DF () ( ) i() F () () ( ) ( ) ( )( )( )( ) ( ) ( )( )( ) c+ a c+ a+ a c a+ c Page 7 of 7