Example 3.3 Moving-average filter

Similar documents
Chapter 5: Sequences & Discrete Difference Equa6ons

Par$al Fac$on Decomposi$on. Academic Resource Center

Unit 5. Matrix diagonaliza1on

ECE 314 Signals and Systems Fall 2012

2 (n 1). f (i 1, i 2,..., i k )

Chapter 5. Fourier Analysis for Discrete-Time Signals and Systems Chapter

Lecture 2: Transfer theory. Credit: Kees Dullemond

ECE 8440 Unit 17. Parks- McClellan Algorithm

Smith Chart The quarter-wave transformer

Sec$on Summary. Sequences. Recurrence Relations. Summations. Ex: Geometric Progression, Arithmetic Progression. Ex: Fibonacci Sequence

Other types of errors due to using a finite no. of bits: Round- off error due to rounding of products

EESC 9945 Geodesy with the Global Posi6oning System. Class 2: Satellite orbits

ODEs + Singulari0es + Monodromies + Boundary condi0ons. Kerr BH ScaRering: a systema0c study. Schwarzschild BH ScaRering: Quasi- normal modes

Bellman s Curse of Dimensionality

Lecture 7 Rolling Constraints

Nonlinear ensemble data assimila/on in high- dimensional spaces. Peter Jan van Leeuwen Javier Amezcua, Mengbin Zhu, Melanie Ades

14. HOMOGENEOUS FINITE STRAIN

MCS 256 Discrete Calculus and Probability Exam 5: Final Examination 22 May 2007

Rankine- Hugoniot Equa0ons and Detona0ons. Forman A. Williams UCSD

CSCI1950 Z Computa4onal Methods for Biology Lecture 4. Ben Raphael February 2, hhp://cs.brown.edu/courses/csci1950 z/ Algorithm Summary

Unit 7 Exponential Functions. Mrs. Valen+ne Math III

Two common difficul:es with HW 2: Problem 1c: v = r n ˆr

ECE 8440 Unit 13 Sec0on Effects of Round- Off Noise in Digital Filters

Unit 6 Sequences. Mrs. Valen+ne CCM3

Deriva'on of The Kalman Filter. Fred DePiero CalPoly State University EE 525 Stochas'c Processes

SYSTEMS AND TENSES. English for Scien/sts Maria Cris/na Teodorani

Modelling of Equipment, Processes, and Systems

Class Notes. Examining Repeated Measures Data on Individuals

Lecture 4 Introduc-on to Data Flow Analysis

Numerical Methods in Physics

MA 123 (Calculus I) Lecture 6: September 19, 2016 Section A3. Professor Joana Amorim,

(5.5) Multistep Methods

Sec$on The Use of Exponen$al Weigh$ng Exponen$al weigh$ng of a sequence x(n) is defined by. (equa$on 13.70)

Decision Making in Robots and Autonomous Agents

Lecture 12 The Level Set Approach for Turbulent Premixed Combus=on

Methodological Foundations of Biomedical Informatics (BMSC-GA 4449) Optimization

Homework 2 Solution Section 2.2 and 2.3.

Exponen'al growth and differen'al equa'ons

Molecular Programming Models. Based on notes by Dave Doty

Binary black holes and gravita4onal waves: Introduc4on to SpEC

Alfred University Dept. of Physics and Astronomy

Pseudospectral Methods For Op2mal Control. Jus2n Ruths March 27, 2009

The BCS Model. Sara Changizi. This presenta5on closely follows parts of chapter 6 in Ring & Schuck The nuclear many body problem.

Machine Learning and Data Mining. Linear regression. Prof. Alexander Ihler

CH 302 Unit 3 Exam Review NUCLEAR, KINETICS, INORGANIC

Lesson 76 Introduction to Complex Numbers

Reduced Models for Process Simula2on and Op2miza2on

Z Transform (Part - II)

Least Square Es?ma?on, Filtering, and Predic?on: ECE 5/639 Sta?s?cal Signal Processing II: Linear Es?ma?on

What is the ph of a 1.0 L 1M HCl solu7on that has reacted with 1.0L of a 1M NaOH solu7on?

Tangent lines, cont d. Linear approxima5on and Newton s Method

Least Squares Parameter Es.ma.on

Exponen'al func'ons and exponen'al growth. UBC Math 102

Eirik Endeve. Simula'ons of SASI, turbulence, and magne'c field amplifica'on

Adsorptive Transport of Noble Gas Tracers in Porous Media

7. Quantum Monte Carlo (QMC)

Introduc)on to Ar)ficial Intelligence

Chapter 11 Liquids and Intermolecular Forces

Noncommuta5ve Black Holes at the LHC

Least Mean Squares Regression. Machine Learning Fall 2017

Nonlinear Shi, Registers: A Survey and Open Problems. Tor Helleseth University of Bergen NORWAY

System Theory- Based Iden2fica2on of Dynamical Models and Applica2ons

Basic Mathema,cs. Rende Steerenberg BE/OP. CERN Accelerator School Basic Accelerator Science & Technology at CERN 3 7 February 2014 Chavannes de Bogis

Modeling radiocarbon in the Earth System. Radiocarbon summer school 2012

Machine Learning & Data Mining CS/CNS/EE 155. Lecture 8: Hidden Markov Models

ASTRO 310: Galac/c & Extragalac/c Astronomy Prof. Jeff Kenney

Unit 3: Ra.onal and Radical Expressions. 3.1 Product Rule M1 5.8, M , M , 6.5,8. Objec.ve. Vocabulary o Base. o Scien.fic Nota.

Local search algorithms

Introduction to Particle Filters for Data Assimilation

Cockcro: Ins<tute Postgraduate Lectures Numerical Methods and LaAce Design

Decision Making in Robots and Autonomous Agents

Streaming - 2. Bloom Filters, Distinct Item counting, Computing moments. credits:

Machine Learning & Data Mining CS/CNS/EE 155. Lecture 11: Hidden Markov Models

Review: Basic Concepts

Do we understand the origin of the IMF?

Fourier Series and Transforms

Accumula'on-rate changes through the Holocene from ice-core records and radar data. Michelle Koutnik Earth and Space Sciences

d ORTHOGONAL POLYNOMIALS AND THE WEYL ALGEBRA

Regression Part II. One- factor ANOVA Another dummy variable coding scheme Contrasts Mul?ple comparisons Interac?ons

Unit 2. Projec.ons and Subspaces

LAB # 3 HANDOUT. y(n)=t[x(n)] Discrete systems have different classifications based on the input and output relation:

Group B Heat Flow Across the San Andreas Fault. Chris Trautner Bill Savran

Lecture 3d: Atmospheric Radia4on & Remote Sensing

Vacuum Systems for Scien1fic Modeling

CS6750: Cryptography and Communica7on Security

Reversible vs. irreversible processes

NETWORK CALCULATIONS updated 11/5/13 1:02 PM

Regression.

Tsybakov noise adap/ve margin- based ac/ve learning

The par/cle-hole picture

Cosmological N-Body Simulations and Galaxy Surveys

From differen+al equa+ons to trigonometric func+ons. Introducing sine and cosine

Density Matrix Description of NMR BCMB/CHEM 8190

Chapter 4. Fourier Analysis for Con5nuous-Time Signals and Systems Chapter Objec5ves

Unit 1 Part 1. Order of Opera*ons Factoring Frac*ons

DEplot(D(y)(x)=2*sin(x*y(x)),y(x),x=-2..2,[[y(1)=1]],y=-5..5)

Cassini Mission. Cassini Spacecra.

Z - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.

CSE446: Linear Regression Regulariza5on Bias / Variance Tradeoff Winter 2015

2.1 Using Models to Multiply Integers (pp )

Transcription:

3.3 Difference Equa/ons for Discrete-Time Systems A Discrete-/me system can be modeled with difference equa3ons involving current, past, or future samples of input and output signals Example 3.3 Moving-average filter A length-n moving average filter is a simple system that produces an output equal to the arithme=c average of the most recent N samples of the input signal.

3.3 Difference Equa/ons for Discrete-Time Systems Example 3.3 Moving-average filter

3.3 Difference Equa/ons for Discrete-Time Systems Example 3.3 Moving-average filter

3.3 Difference Equa/ons for Discrete-Time Systems Example 3.4 Length-2 moving-average filter

3.3 Difference Equa/ons for Discrete-Time Systems Example 3.4 Length-2 moving-average filter

3.3 Difference Equa/ons for Discrete-Time Systems Example 3.6 Exponen=al smoother

3.3 Difference Equa/ons for Discrete-Time Systems Example 3.6 Exponen=al smoother

3.3 Difference Equa/ons for Discrete-Time Systems Example 3.7 Loan payments Express the act of taking a bank loan an paying it back over =me as a discrete-=me system problem with monthly payments represen=ng the Input signal and monthly balance represen=ng the output signal. y[n] = (1+c)y[n-1] + x[n] c is the monthly interest rate

3.3 Difference Equa/ons for Discrete-Time Systems Example 3.7 Loan payments Express the act of taking a bank loan an paying it back over =me as a discrete-=me system problem with monthly payments represen=ng the Input signal and monthly balance represen=ng the output signal. y[n] = (1+c)y[n-1] + x[n] c is the monthly interest rate We can find out the monthly payment x[n] by solving the difference equa=on

3.4 Constant-Coefficient Linear Difference Equa/ons Constant-coefficient difference equa3on N a k y[n k] = b k x[n k] k=0 k=0 Ini/al condi/ons: y[n 0-1], y[n 0-2],. y[n 0 -N] It is typical, but not required, to have n 0 = 0. M

3.4 Constant-Coefficient Linear Difference Equa/ons Itera3vely solving difference equa3ons. Consider the difference equa=on for the exponen=al smoother of Ex.3.6 y[n] = (1- α)y[n-1] + α x[n] ² Given the ini=al value y[-1] of the output signal, y[0] is found by y[0] = (1- α)y[-1] + α x[0] ² Knowing y[0], the next output sample y[1] is found by y[1] = (1- α)y[0] + α x[1] ² Knowing y[1], the next output sample y[2] is found by y[2] = (1- α)y[1] + α x[2]

3.5 Solving Difference Equa/ons N a k y[n k] = b k x[n k] k=0 M k=0 Ini=al condi=ons: y[n 0 1], y[n 0 2],,.., y[n 0 N ] General solu=on: y[n] = y h [n]+ y p [n] ² y n ([n]: is the homogeneous solu=on of the differen=al equa=on natural response ² y p [n]: is the par3cular solu3on of the differen=al equa=on ² y[n] = y h [n] + y p [n] is the forced solu3on of the differen=al equa=on forced response

3.5 Solving Difference Equa/ons Example 3.11 Natural response of exponen=al smoother Determine the natural response of the exponen=al smoother defined as: y[n] = (1 α) y[n 1]+αx[n] If y[-1] = 2

3.5.1 Finding the natural response of a discrete- /me system General homogeneous differen=al equa=on: Characteris/c equa/on: N a k y[n k] = 0 k=0 a k z k = 0 k=0 To obtain the characteris=c equa=on, subs=tute: N y[n k] z k

3.5.1 Finding the natural response of a discrete- /me system Write the characteris=c equa=on in open form: a 0 + a 1 z 1 +...+ a N 1 z N+1 + a N z N = 0 Mul=ply both sides by z N to obtain a 0 z N + a 1 z N 1 +...+ a N 1 z 1 + a N = 0 In factored form: a 0 (a z 1 )(a z 2 )...(a z N ) = 0 Homogeneous solu=on (assuming roots are dis=nct): y h [n] = c 1 z 1 n + c 2 z 2 n +...+ c N z N n = N k=1 c k z k n Where unknown coefficients c 1, c 2,., c N are determined by the ini=al condi=on, the terms are called the modes of the system z i n

3.5.1 Finding the natural response of a discrete- /me system Example 3.12 Natural response of second-order system A second-order system is described by the difference equa=on: y[n] 5 6 y[n 1]+ 1 6 y[n 2] = 0 Determine the natural response of this system for n >= 0 subject to ini=al condi=on: y[-1] = 19, and y[-2] = 53

3.5.2 Finding the forced response of a discrete- /me system Choosing a par3cular solu3on for various discrete-=me input signals The coefficient of the par3cular solu3on are determined from the Difference equa=on by assuming all ini=al condi=ons are equal to zero.

3.5.2 Finding the forced response of a discrete- /me system Write the homogeneous difference equa=on Solve homogeneous difference equa=on with undetermined coefficient Find the form of the par=cular solu=on Find the coefficients of the par=cular solu=on Add the homogeneous and par=cular solu=on Together to obtain the total solu=on

3.5.2 Finding the forced response of a discrete- /me system Example 3.14 Forced response of exponen=al smoother Determine the forced response of the exponen=al smoother defined as: y[n] = (1 α) y[n 1]+αx[n] The input signal is a unit-step func=on, and y[-1] = 2.5