Frequency Estimation from two DFT Bins
|
|
- Evan Ashley Foster
- 5 years ago
- Views:
Transcription
1 Frequency Estimation from two DFT Bins Martin Vicanek 17. June 015, revised 31. May 016 and 11. February 017 In a recent thread in comp.dsp [1] Cedron Dawg presented an exact formula for extracting the frequency of a pure tone from three DFT bins []. The possibility of using fewer bins remained an open question, also whether or not there are other three-bin formulas. Inspired by that thread I was able to derive an exact expression for the frequency based on values from only two bins. It turns out that that expression is not unique: it offers a degree of freedom which may be used to advantage in the presence of noise. Consider a frame of real valued samples x 0, x 1,..., x N 1 which represent a pure tone with some frequency, amplitude, and phase. We may use the center position n c = (N 1)/ to write x n = A cos[(n n c )ω] + B sin[(n n c )ω], (1) where ω denotes a normalized angular frequency in the range 0 ω π with ω = π representing the Nyquist frequency. A and B are real valued constants to allow for an arbitrary amplitude A + B and phase arctan(b/a). The DFT may be written in a form centered around n c, N 1 X k = e πiknc/n n=0 x n e πikn/n, () where the sum on the right hand side represents the familiar DFT and the phase factor in front provides for the symmetry with respect to n c. Inserting (1) into () yields, after some algebra, X k = ( 1) k sin(nω/) cos β k cos ω [ ( ω ) A sin cos ( βk where we have introduced the abbreviation β k intermediate steps at the end of this note. 1 ) ( ω ) ( )] βk + ib cos sin, (3) = πk/n. We outline the
2 Suppose the set of X k is known for all k = 0, 1,..., N 1, and the task is to extract the frequency ω, and possibly A and B. To this end we define the quantities u k = ( 1) k Re(X k )/ cos(β k /), k N/. (4) It is easy to see that the following relation holds for any k j: u j = cos β k cos ω u k cos β j cos ω. (5) Eq.(5) may be solved for ω to yield ( ) uj cos β j u k cos β k ω = arccos u j u k (6) So far so good: equation (6) provides a formula for the frequency ω using information from the two bins j and k. Note that j and k need not denote adjacent bins, although it is reasonable to choose the two bins with maximum modulus, which, in the case of a pure tone spectrum, will be adjacent bins. Also note that these two bins will enclose ω, and that the corresponding u j and u k will have oposite signs. Eq.(5) will break down for A = 0, however. In that case we may define v k = ( 1) k Im(X k )/ sin(β k /), k 0, (7) and obtain similar formulas as in eqs.(5) and (6). Indeed, any linear combination w k = au k + bv k (8) will do. This degree of freedom is not so surprising given that two bins hold two complex valued numbers (i.e. four independent values) to determine three unknowns ω, A and B. By the same token, one would expect Cedron s three-bin formula to be only one among many others. Clearly, we can use this freedom to our advantage and make an optimum choice of e.g. a and b for best signal to noise ratio. Although it is not obvious what the best choice would be, setting a = cos(β m /) Re(X m ) b = sin(β m /) Im(X m ) (9) with m denoting the peak bin optimizes the SNR of w m. Hence that choice is expected to provide nearly optimum accuracy.
3 It would seem that because only two bins are used, much information is actually thrown away - information that could be used if the background characteristics are known (white noise being the simplest case). On the other hand, for unknown background one might argue that a two-bin method provides maximum locality hence minimum interference. I have tested the scheme only briefly in a real-world scenario with noise or another interfering tone. Results with w k according to eq.(9) are indeed better than with either u k or v k plugged into eq.(6). Derivation of Eq.(3) We may write eq.(1) in the form x n = A ib The DFT of the exponentials in eq.(10) is e iω(n nc) + A + ib e iω(n nc). (10) N 1 n=0 e i(β k±ω)(n n c) = sin [N(β k ± ω)] sin[(β k ± ω)/] = ( 1)k sin(nω/) sin[(ω ± β k )/]. (11) Hence we obtain the following expression for the real part: ( ) { } Re(X k ) = ( 1)k A Nω 1 sin sin[(ω + β k )/] + 1 sin[(ω β k )/] (1) The term in curly brackets in eq.(1) may be simplified using trigonometric identities, resulting in 4 sin(ω/) cos(β k/) cos β k cos ω. (13) The imaginary part Im(X k ) may be treated in a similar way. Collecting all bits and pieces, we obtain the result in eq.(3). Derivation of Optimum a and b If we assume additive white noise, then the Fourier components will fluctuate with some variance σ, regardless of k, and the variance will be equal for the real and imaginary parts: Var[Re(X k )] = Var[Im(X k )] = σ. (14) 3
4 Hence the variances of u k and v k will be Var(u k ) = which yields the variance of w k, [ Var(w k ) = σ cos (β k /), Var(v k) = a cos (β k /) + b sin (β k /) σ sin (β k /), (15) ] σ (16) Now the relative mean standard deviation of w k is Var(w k )/ w k. We want the relative error of w k to be small because we are taking ratios. Minimizing the above expression (or easier, its square) results in a condition for a and b, Minimum Variance Estimate a b = cos(β k/)re(x k ) sin(β k /)Im(X k ). (17) The information from two bins leads to an overdetermined system of equations for ω. With the abreviations µ = cos ω, and µ k = cos β k, these equations are: U(µ) := (u j u k )µ u j µ j + u k µ k = 0 V (µ) := (v j v k )µ v j µ j + v k µ k = 0 (18) In general, the two equations cannot be solved simultaneously. The best one can do is to minimize a suitable error. One may square each equation and minimize their weighted sum. The optimum weights are the inverse variances of each term. Observe that for adjacent j and k at the peak, Var(U) σ / cos (β m /) Var(V ) σ / sin (β m /) (19) where the index m denotes the peak. Hence the optimum cost function is cos (β m /) U(µ) + sin (β m /) V (µ) = min. (0) Taking the derivative with respect to µ and equating to zero, we arrive at the expression µ = cos ω = cos (β m /)(u j u k ) µ u + sin (β m /)(v j v k ) µ v cos (β m /)(u j u k ) + sin (β m /)(v j v k ) [Re(X j + X k )] µ u + [Im(X j + X k )] µ v X j + X k (1) 4
5 where µ u and µ v are the solutions of U(µ) = 0 and V (µ) = 0, respectively. This result coincides with the one found in the previous section if we substitute in eq. (9) u m u j u k / and v m v j v k /. Recall that u j and u k have opposite signs, likewise for v j and v k. Formulas in Terms of Ordinary DFT For easier application of the results we provide the main expressions in terms of the standard DFT which we denote by Z k, Z k = N 1 n=0 x n e iβ kn. () The relation between the centered and the standard DFT is ( 1) k Re(X k ) = cos(β k /) Re(Z k ) + sin(β k /) Im(Z k ) ( 1) k Im(X k ) = sin(β k /) Re(Z k ) + cos(β k /) Im(Z k ), (3) The auxiliary quantities u k and v k become u k = Re(Z k ) + tan(β k /) Im(Z k ) and the optimum choice for a and b is Complex Tone v k = Re(Z k ) + cot(β k /) Im(Z k ), (4) a b = cot(β m/) Re(Z m ) + Im(Z m ) tan(β m /) Re(Z m ) + Im(Z m ). (5) As an alternative to eq.(1), consider a complex tone x n = Ce i(n nc)ω, (6) where C represents a complex amplitude. The DFT according to eq.() is X k = ( 1) k C sin(nω/) sin[(β k ω)/]. (7) To derive a two-bin formula for a complex tone, we form the ratio ( 1) j X j = sin[(β k ω)/] ( 1) k X k sin[(β j ω)/] = cos(ω/) sin(β k/) sin(ω/) cos(β k /) cos(ω/) sin(β j /) sin(ω/) cos(β j /) (8) 5
6 which is easily solved for tan(ω/) = ( 1)j X j sin(β j /) ( 1) k X k sin(β k /) ( 1) j X j cos(β j /) ( 1) k X k cos(β k /). (9) For adjacent bins j = k + 1 this expression becomes References tan(ω/) = X k+1 sin(β k+1 /) + X k sin(β k /) X k+1 cos(β k+1 /) + X k cos(β k /). (30) [1] [] 6
Fourier Series and Fourier Transforms
Fourier Series and Fourier Transforms EECS2 (6.082), MIT Fall 2006 Lectures 2 and 3 Fourier Series From your differential equations course, 18.03, you know Fourier s expression representing a T -periodic
More informationSystems. x 1 = ax 1 + bx 2, x 2 = cx 1 + dx 2,
Systems We consider systems of differential equations of the fm x ax + bx, x cx + dx, where a,b,c, and d are real numbers At first glance the system may seem to be first-der; however, it is coupled, and
More information10.1 Complex Arithmetic Argand Diagrams and the Polar Form The Exponential Form of a Complex Number De Moivre s Theorem 29
10 Contents Complex Numbers 10.1 Complex Arithmetic 2 10.2 Argand Diagrams and the Polar Form 12 10.3 The Exponential Form of a Complex Number 20 10.4 De Moivre s Theorem 29 Learning outcomes In this Workbook
More information1 Review of complex numbers
1 Review of complex numbers 1.1 Complex numbers: algebra The set C of complex numbers is formed by adding a square root i of 1 to the set of real numbers: i = 1. Every complex number can be written uniquely
More informationCMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation
CMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University January 23, 2006 1 Exponentials The exponential is
More informationToday s lecture. The Fourier transform. Sampling, aliasing, interpolation The Fast Fourier Transform (FFT) algorithm
Today s lecture The Fourier transform What is it? What is it useful for? What are its properties? Sampling, aliasing, interpolation The Fast Fourier Transform (FFT) algorithm Jean Baptiste Joseph Fourier
More informationOverview of Complex Numbers
Overview of Complex Numbers Definition 1 The complex number z is defined as: z = a+bi, where a, b are real numbers and i = 1. General notes about z = a + bi Engineers typically use j instead of i. Examples
More informationMath 104 Midterm 3 review November 12, 2018
Math 04 Midterm review November, 08 If you want to review in the textbook, here are the relevant sections: 4., 4., 4., 4.4, 4..,.,. 6., 6., 6., 6.4 7., 7., 7., 7.4. Consider a right triangle with base
More informationTable of contents. d 2 y dx 2, As the equation is linear, these quantities can only be involved in the following manner:
M ath 0 1 E S 1 W inter 0 1 0 Last Updated: January, 01 0 Solving Second Order Linear ODEs Disclaimer: This lecture note tries to provide an alternative approach to the material in Sections 4. 4. 7 and
More informationREFLECTION AND REFRACTION OF PLANE EM WAVES
REFLECTION AND REFRACTION OF PLANE EM WAVES When an electromagnetic wave hits a boundary between different materials, some of the wave s energy is reflected back while the rest continues on through the
More informationNotes for the Physics-Based Calculus workshop
Notes for the hysics-based Calculus workshop Adam Coffman June 9, 24 Trigonometric Functions We recall the following identities for trigonometric functions. Theorem.. For all x R, cos( x = cos(x and sin(
More information(x 1, y 1 ) = (x 2, y 2 ) if and only if x 1 = x 2 and y 1 = y 2.
1. Complex numbers A complex number z is defined as an ordered pair z = (x, y), where x and y are a pair of real numbers. In usual notation, we write z = x + iy, where i is a symbol. The operations of
More informationGeneral Physics I. Lecture 14: Sinusoidal Waves. Prof. WAN, Xin ( 万歆 )
General Physics I Lecture 14: Sinusoidal Waves Prof. WAN, Xin ( 万歆 ) xinwan@zju.edu.cn http://zimp.zju.edu.cn/~xinwan/ Motivation When analyzing a linear medium that is, one in which the restoring force
More informationComplex symmetry Signals and Systems Fall 2015
18-90 Signals and Systems Fall 015 Complex symmetry 1. Complex symmetry This section deals with the complex symmetry property. As an example I will use the DTFT for a aperiodic discrete-time signal. The
More informationBackground ODEs (2A) Young Won Lim 3/7/15
Background ODEs (2A) Copyright (c) 2014-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
More informationFinal Exam January 31, Solutions
Final Exam January 31, 014 Signals & Systems (151-0575-01) Prof. R. D Andrea & P. Reist Solutions Exam Duration: Number of Problems: Total Points: Permitted aids: Important: 150 minutes 7 problems 50 points
More informationIB Paper 6: Signal and Data Analysis
IB Paper 6: Signal and Data Analysis Handout 2: Fourier Series S Godsill Signal Processing and Communications Group, Engineering Department, Cambridge, UK Lent 2015 1 / 1 Fourier Series Revision of Basics
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More informationTopic 5 Notes Jeremy Orloff. 5 Homogeneous, linear, constant coefficient differential equations
Topic 5 Notes Jeremy Orloff 5 Homogeneous, linear, constant coefficient differential equations 5.1 Goals 1. Be able to solve homogeneous constant coefficient linear differential equations using the method
More informationHB Coupled Pendulums Lab Coupled Pendulums
HB 04-19-00 Coupled Pendulums Lab 1 1 Coupled Pendulums Equipment Rotary Motion sensors mounted on a horizontal rod, vertical rods to hold horizontal rod, bench clamps to hold the vertical rods, rod clamps
More informationLikewise, any operator, including the most generic Hamiltonian, can be written in this basis as H11 H
Finite Dimensional systems/ilbert space Finite dimensional systems form an important sub-class of degrees of freedom in the physical world To begin with, they describe angular momenta with fixed modulus
More informationComplex 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 informationWe consider systems of differential equations of the form. x 1 = ax 1 + bx 2,
Chapter 5 Systems We consider systems of differential equations of the fm (5) a + b, c + d, where a, b, c, and d are real numbers At first glance the system may seem to be first-der; however, it is coupled,
More informationSolutions 2: Simple Harmonic Oscillator and General Oscillations
Massachusetts Institute of Technology MITES 2017 Physics III Solutions 2: Simple Harmonic Oscillator and General Oscillations Due Wednesday June 21, at 9AM under Rene García s door Preface: This problem
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More informationCHAPTER 4 FOURIER SERIES S A B A R I N A I S M A I L
CHAPTER 4 FOURIER SERIES 1 S A B A R I N A I S M A I L Outline Introduction of the Fourier series. The properties of the Fourier series. Symmetry consideration Application of the Fourier series to circuit
More informationChapter 2 Algorithms for Periodic Functions
Chapter 2 Algorithms for Periodic Functions In this chapter we show how to compute the Discrete Fourier Transform using a Fast Fourier Transform (FFT) algorithm, including not-so special case situations
More informationNotes on Noncoherent Integration Gain
Notes on Noncoherent Integration Gain Mark A. Richards July 2014 1 Noncoherent Integration Gain Consider threshold detection of a radar target in additive complex Gaussian noise. To achieve a certain probability
More informationHongbing Zhang October 2018
Exact and Intuitive Geometry Explanation: Why Does a Half-angle-rotation in Spin or Quaternion Space Correspond to the Whole-angle-rotation in Normal 3D Space? Hongbing Zhang October 2018 Abstract Why
More information3 + 4i 2 + 3i. 3 4i Fig 1b
The introduction of complex numbers in the 16th century was a natural step in a sequence of extensions of the positive integers, starting with the introduction of negative numbers (to solve equations of
More informationSinusoids. Amplitude and Magnitude. Phase and Period. CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation
Sinusoids CMPT 889: Lecture Sinusoids, Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University September 6, 005 Sinusoids are
More informationCMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation
CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University September 26, 2005 1 Sinusoids Sinusoids
More informationDiscrete Simulation of Power Law Noise
Discrete Simulation of Power Law Noise Neil Ashby 1,2 1 University of Colorado, Boulder, CO 80309-0390 USA 2 National Institute of Standards and Technology, Boulder, CO 80305 USA ashby@boulder.nist.gov
More informationω 0 = 2π/T 0 is called the fundamental angular frequency and ω 2 = 2ω 0 is called the
he ime-frequency Concept []. Review of Fourier Series Consider the following set of time functions {3A sin t, A sin t}. We can represent these functions in different ways by plotting the amplitude versus
More informationMulti-Robotic Systems
CHAPTER 9 Multi-Robotic Systems The topic of multi-robotic systems is quite popular now. It is believed that such systems can have the following benefits: Improved performance ( winning by numbers ) Distributed
More informationLOPE3202: Communication Systems 10/18/2017 2
By Lecturer Ahmed Wael Academic Year 2017-2018 LOPE3202: Communication Systems 10/18/2017 We need tools to build any communication system. Mathematics is our premium tool to do work with signals and systems.
More informationMusic 270a: Complex Exponentials and Spectrum Representation
Music 270a: Complex Exponentials and Spectrum Representation Tamara Smyth, trsmyth@ucsd.edu Department of Music, University of California, San Diego (UCSD) October 24, 2016 1 Exponentials The exponential
More information8/19/16. Fourier Analysis. Fourier analysis: the dial tone phone. Fourier analysis: the dial tone phone
Patrice Koehl Department of Biological Sciences National University of Singapore http://www.cs.ucdavis.edu/~koehl/teaching/bl5229 koehl@cs.ucdavis.edu Fourier analysis: the dial tone phone We use Fourier
More informationz = x + iy ; x, y R rectangular or Cartesian form z = re iθ ; r, θ R polar form. (1)
11 Complex numbers Read: Boas Ch. Represent an arb. complex number z C in one of two ways: z = x + iy ; x, y R rectangular or Cartesian form z = re iθ ; r, θ R polar form. (1) Here i is 1, engineers call
More informationSection 3: Complex numbers
Essentially: Section 3: Complex numbers C (set of complex numbers) up to different notation: the same as R 2 (euclidean plane), (i) Write the real 1 instead of the first elementary unit vector e 1 = (1,
More informationFourier s representation for functions on R, T p, Z, and P N
1 Fourier s representation for functions on R, T p, Z, and P N 1.1 Synthesis and analysis equations Introduction In mathematics we often try to synthesize a rather arbitrary function f using a suitable
More informationEach of these functions represents a signal in terms of its spectral components in the frequency domain.
N INTRODUCTION TO SPECTRL FUNCTIONS Revision B By Tom Irvine Email: tomirvine@aol.com March 3, 000 INTRODUCTION This tutorial presents the Fourier transform. It also discusses the power spectral density
More informationMAT187H1F Lec0101 Burbulla
Spring 2017 Second Order Linear Homogeneous Differential Equation DE: A(x) d 2 y dx 2 + B(x)dy dx + C(x)y = 0 This equation is called second order because it includes the second derivative of y; it is
More informationANoteonEuler sformula
ANoteonEuler sformula Dr. Mike Wilkes ASC Chastain Indian River State College 4/1/014 1 Euler s Formula The general complex exponential function e z,wherez is any complex number of the form (a + ib), has
More informationWith this expanded version of what we mean by a solution to an equation we can solve equations that previously had no solution.
M 74 An introduction to Complex Numbers. 1 Solving equations Throughout the calculus sequence we have limited our discussion to real valued solutions to equations. We know the equation x 1 = 0 has distinct
More informationIn Z: x + 3 = 2 3x = 2 x = 1 No solution In Q: 3x = 2 x 2 = 2. x = 2 No solution. In R: x 2 = 2 x = 0 x = ± 2 No solution Z Q.
THE UNIVERSITY OF NEW SOUTH WALES SCHOOL OF MATHEMATICS AND STATISTICS MATH 1141 HIGHER MATHEMATICS 1A ALGEBRA. Section 1: - Complex Numbers. 1. The Number Systems. Let us begin by trying to solve various
More informationNORMAL MODES, WAVE MOTION AND THE WAVE EQUATION. Professor G.G.Ross. Oxford University Hilary Term 2009
NORMAL MODES, WAVE MOTION AND THE WAVE EQUATION Professor G.G.Ross Oxford University Hilary Term 009 This course of twelve lectures covers material for the paper CP4: Differential Equations, Waves and
More informationLecture 04: Discrete Frequency Domain Analysis (z-transform)
Lecture 04: Discrete Frequency Domain Analysis (z-transform) John Chiverton School of Information Technology Mae Fah Luang University 1st Semester 2009/ 2552 Outline Overview Lecture Contents Introduction
More informationEECS 20N: Midterm 2 Solutions
EECS 0N: Midterm Solutions (a) The LTI system is not causal because its impulse response isn t zero for all time less than zero. See Figure. Figure : The system s impulse response in (a). (b) Recall that
More informationMAE 200B Homework #3 Solutions University of California, Irvine Winter 2005
Problem 1 (Haberman 5.3.2): Consider this equation: MAE 200B Homework #3 Solutions University of California, Irvine Winter 2005 a) ρ 2 u t = T 2 u u 2 0 + αu + β x2 t The term αu describes a force that
More informationModes and Roots ... mx + bx + kx = 0. (2)
A solution of the form x(t) = ce rt to the homogeneous constant coefficient linear equation a x (n) + (n 1). n a n 1 x + + a 1 x + a 0 x = 0 (1) is called a modal solution and ce rt is called a mode of
More informationMATH 135: COMPLEX NUMBERS
MATH 135: COMPLEX NUMBERS (WINTER, 010) The complex numbers C are important in just about every branch of mathematics. These notes 1 present some basic facts about them. 1. The Complex Plane A complex
More informationChapter 2: Complex numbers
Chapter 2: Complex numbers Complex numbers are commonplace in physics and engineering. In particular, complex numbers enable us to simplify equations and/or more easily find solutions to equations. We
More informationEuler-Cauchy Using Undetermined Coefficients
Euler-Cauchy Using Undetermined Coefficients Department of Mathematics California State University, Fresno doreendl@csufresno.edu Joint Mathematics Meetings January 14, 2010 Outline 1 2 3 Second Order
More information01 Harmonic Oscillations
Utah State University DigitalCommons@USU Foundations of Wave Phenomena Library Digital Monographs 8-2014 01 Harmonic Oscillations Charles G. Torre Department of Physics, Utah State University, Charles.Torre@usu.edu
More informationHypothesis elimination on a quantum computer
Hypothesis elimination on a quantum computer arxiv:quant-ph/042025v 3 Dec 2004 Andrei N. Soklakov and Rüdiger Schack Department of athematics, Royal Holloway, University of London Egham, Surrey, TW20 0EX,
More informationDavid Abrecht. February 17, 2012
Addendum to An analytical solution method for the unsteady, unbounded, incompressible three-dimensional Navier-Stokes equations in Cartesian coordinates using coordinate axis symmetry degeneracy David
More informationLab10: FM Spectra and VCO
Lab10: FM Spectra and VCO Prepared by: Keyur Desai Dept. of Electrical Engineering Michigan State University ECE458 Lab 10 What is FM? A type of analog modulation Remember a common strategy in analog modulation?
More informationQuick Overview: Complex Numbers
Quick Overview: Complex Numbers February 23, 2012 1 Initial Definitions Definition 1 The complex number z is defined as: z = a + bi (1) where a, b are real numbers and i = 1. Remarks about the definition:
More information2 Fourier Transforms and Sampling
2 Fourier ransforms and Sampling 2.1 he Fourier ransform he Fourier ransform is an integral operator that transforms a continuous function into a continuous function H(ω) =F t ω [h(t)] := h(t)e iωt dt
More information1.1 A Scattering Experiment
1 Transfer Matrix In this chapter we introduce and discuss a mathematical method for the analysis of the wave propagation in one-dimensional systems. The method uses the transfer matrix and is commonly
More informationDamped Harmonic Oscillator
Damped Harmonic Oscillator Wednesday, 23 October 213 A simple harmonic oscillator subject to linear damping may oscillate with exponential decay, or it may decay biexponentially without oscillating, or
More informationMath 112 Section 10 Lecture notes, 1/7/04
Math 11 Section 10 Lecture notes, 1/7/04 Section 7. Integration by parts To integrate the product of two functions, integration by parts is used when simpler methods such as substitution or simplifying
More informationChapter 7: Trigonometric Equations and Identities
Chapter 7: Trigonometric Equations and Identities In the last two chapters we have used basic definitions and relationships to simplify trigonometric expressions and equations. In this chapter we will
More informationPhysics 142 Wave Optics 1 Page 1. Wave Optics 1. For every complex problem there is one solution that is simple, neat, and wrong. H.L.
Physics 142 Wave Optics 1 Page 1 Wave Optics 1 For every complex problem there is one solution that is simple, neat, and wrong. H.L. Mencken Interference and diffraction of waves The essential characteristic
More informationChapter 7: Trigonometric Equations and Identities
Chapter 7: Trigonometric Equations and Identities In the last two chapters we have used basic definitions and relationships to simplify trigonometric expressions and equations. In this chapter we will
More information23.6. The Complex Form. Introduction. Prerequisites. Learning Outcomes
he Complex Form 3.6 Introduction In this Section we show how a Fourier series can be expressed more concisely if we introduce the complex number i where i =. By utilising the Euler relation: e iθ cos θ
More informationBME 50500: Image and Signal Processing in Biomedicine. Lecture 2: Discrete Fourier Transform CCNY
1 Lucas Parra, CCNY BME 50500: Image and Signal Processing in Biomedicine Lecture 2: Discrete Fourier Transform Lucas C. Parra Biomedical Engineering Department CCNY http://bme.ccny.cuny.edu/faculty/parra/teaching/signal-and-image/
More information7.3 Getting on the Right Wavelength
7.3 Getting on the Right Wavelength A Practice Understanding Task The Ferris wheel in the following diagram has a radius of 40 feet, its center is 50 feet from the ground, and it makes one revolution counterclockwise
More information2. Force Systems. 2.1 Introduction. 2.2 Force
2. Force Systems 2.1 Introduction 2.2 Force - A force is an action of one body on another. - A force is an action which tends to cause acceleration of a body (in dynamics). - A force is a vector quantity.
More informationImage Acquisition and Sampling Theory
Image Acquisition and Sampling Theory Electromagnetic Spectrum The wavelength required to see an object must be the same size of smaller than the object 2 Image Sensors 3 Sensor Strips 4 Digital Image
More informationFourier Series and Integrals
Fourier Series and Integrals Fourier Series et f(x) beapiece-wiselinearfunctionon[, ] (Thismeansthatf(x) maypossessa finite number of finite discontinuities on the interval). Then f(x) canbeexpandedina
More informationThe Wave Function. Chapter The Harmonic Wave Function
Chapter 3 The Wave Function On the basis of the assumption that the de Broglie relations give the frequency and wavelength of some kind of wave to be associated with a particle, plus the assumption that
More informationSection 4.7: Variable-Coefficient Equations
Cauchy-Euler Equations Section 4.7: Variable-Coefficient Equations Before concluding our study of second-order linear DE s, let us summarize what we ve done. In Sections 4.2 and 4.3 we showed how to find
More informationProblem Set 4 for MAE280A Linear Systems Theory, Fall 2017: due Tuesday November 28, 2017, in class
Problem Set 4 for MAE8A Linear Systems Theory, Fall 7: due Tuesday November 8, 7, in class Problem Controllability and reachability Consider the C-circuit depicted in class on Slide 3 of the Controllability
More informationNONLINEAR COMPLEX MODULUS IN BITUMENS
NONLINEAR COMPLEX MODULUS IN BITUMENS J. Stastna, K. Jorshari and L. Zanzotto Bituminous Materials Chair University of Calgary 2500 University Drive NW Calgary, Alberta T2N 1N4 Canada Large Amplitude Oscillations
More informationCOMPLEX NOISE-BITS AND LARGE-SCALE INSTANTANEOUS PARALLEL OPERATIONS WITH LOW COMPLEXITY
COMPLEX NOISE-BITS AND LARGE-SCALE INSTANTANEOUS PARALLEL OPERATIONS WITH LOW COMPLEXITY HE WEN (a, b, LASZLO B. KISH (a, AND ANDREAS KLAPPENECKER (c (a Texas A&M University, Department of Electrical and
More informationLecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More informationA523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2011
A523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2011 Lecture 6 PDFs for Lecture 1-5 are on the web page Problem set 2 is on the web page Article on web page A Guided
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 9 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,
More informationJackson 7.6 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jackson 7.6 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell PROBLEM: A plane wave of frequency ω is incident normally from vacuum on a semi-infinite slab of material
More informationContinuous functions. Limits of non-rational functions. Squeeze Theorem. Calculator issues. Applications of limits
Calculus Lia Vas Continuous functions. Limits of non-rational functions. Squeeze Theorem. Calculator issues. Applications of limits Continuous Functions. Recall that we referred to a function f() as a
More informationContents. Digital Signal Processing, Part II: Power Spectrum Estimation
Contents Digital Signal Processing, Part II: Power Spectrum Estimation 5. Application of the FFT for 7. Parametric Spectrum Est. Filtering and Spectrum Estimation 7.1 ARMA-Models 5.1 Fast Convolution 7.2
More informationThe Exponential of a Matrix
The Exponential of a Matrix 5-8- The solution to the exponential growth equation dx dt kx is given by x c e kt It is natural to ask whether you can solve a constant coefficient linear system x A x in a
More information2.161 Signal Processing: Continuous and Discrete
MI OpenCourseWare http://ocw.mit.edu.6 Signal Processing: Continuous and Discrete Fall 8 For information about citing these materials or our erms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSES INSIUE
More information21.3. z-transforms and Difference Equations. Introduction. Prerequisites. Learning Outcomes
-Transforms and Difference Equations 1.3 Introduction In this we apply -transforms to the solution of certain types of difference equation. We shall see that this is done by turning the difference equation
More informationFrequency Response of Linear Time Invariant Systems
ME 328, Spring 203, Prof. Rajamani, University of Minnesota Frequency Response of Linear Time Invariant Systems Complex Numbers: Recall that every complex number has a magnitude and a phase. Example: z
More informationFourier Series and Transforms
Fourier Series and Transforms Website is now online at: http://www.unc.edu/courses/2008fall/comp/665/001/ 9/2/08 Comp 665 Real and Special Signals 1 Discrete Exponen8al Func8on Discrete Convolu?on: Convolu?on
More informationDCSP-2: Fourier Transform
DCSP-2: Fourier Transform Jianfeng Feng Department of Computer Science Warwick Univ., UK Jianfeng.feng@warwick.ac.uk http://www.dcs.warwick.ac.uk/~feng/dcsp.html Data transmission Channel characteristics,
More informationDIFFERENCE AND DIFFERENTIAL EQUATIONS COMPARED
DIFFERENCE AND DIFFERENTIAL EQUATIONS COMPARED Frequency Response Function: Continuous Systems Given the unit impulse response h(t) of a continuous-time system, the mapping from the system s input x(t)
More informationEECS 123 Digital Signal Processing University of California, Berkeley: Fall 2007 Gastpar November 7, Exam 2
EECS 3 Digital Signal Processing University of California, Berkeley: Fall 7 Gastpar November 7, 7 Exam Last name First name SID You have hour and 45 minutes to complete this exam. he exam is closed-book
More informationDigital Signal Processing
Digital Signal Processing Introduction Moslem Amiri, Václav Přenosil Embedded Systems Laboratory Faculty of Informatics, Masaryk University Brno, Czech Republic amiri@mail.muni.cz prenosil@fi.muni.cz February
More information7.3 Inverse Trigonometric Functions
58 transcendental functions 73 Inverse Trigonometric Functions We now turn our attention to the inverse trigonometric functions, their properties and their graphs, focusing on properties and techniques
More informationELEG 305: Digital Signal Processing
ELEG 5: Digital Signal Processing Lecture 6: The Fast Fourier Transform; Radix Decimatation in Time Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 8 K.
More informationA NOVEL METHOD TO DERIVE EXPLICIT KLT KERNEL FOR AR(1) PROCESS. Mustafa U. Torun and Ali N. Akansu
A NOVEL METHOD TO DERIVE EXPLICIT KLT KERNEL FOR AR() PROCESS Mustafa U. Torun and Ali N. Akansu Department of Electrical and Computer Engineering New Jersey Institute of Technology University Heights,
More informationEXAMINATION: MATHEMATICAL TECHNIQUES FOR IMAGE ANALYSIS
EXAMINATION: MATHEMATICAL TECHNIQUES FOR IMAGE ANALYSIS Course code: 8D Date: Thursday April 8 th, Time: 4h 7h Place: AUD 3 Read this first! Write your name and student identification number on each paper
More informationSecond-Order Linear ODEs
C0.tex /4/011 16: 3 Page 13 Chap. Second-Order Linear ODEs Chapter presents different types of second-order ODEs and the specific techniques on how to solve them. The methods are systematic, but it requires
More informationPartial Fractions. June 27, In this section, we will learn to integrate another class of functions: the rational functions.
Partial Fractions June 7, 04 In this section, we will learn to integrate another class of functions: the rational functions. Definition. A rational function is a fraction of two polynomials. For example,
More informationProperties of Fourier Series - GATE Study Material in PDF
Properties of Fourier Series - GAE Study Material in PDF In the previous article, we learnt the Basics of Fourier Series, the different types and all about the different Fourier Series spectrums. Now,
More informationWaves Encountering Barriers
Waves Encountering Barriers Reflection and Refraction: When a wave is incident on a boundary that separates two regions of different wave speed, part of the wave is reflected and part is transmitted. Figure
More information