Recalculation of the derivation of the image source model for a rectangular room with rigid boundaries
|
|
- Evangeline Watts
- 5 years ago
- Views:
Transcription
1 Recalculation of the derivation of the image source model for a rectangular room with rigid boundaries Vera Erbes Universität Rostock, Institute for Communications Engineering, R.-Wagner-Straße, D Rostock, Germany vera.erbes@uni-rostock.de Fabian Brinkmann TU Berlin, Audio Communication Group, Einsteinufer 17c, D Berlin, Germany fabian.brinkmann@tu-berlin.de March 1, 018
2 Remarks This document is meant to be a helpful resource to comprehend the derivation of the image source model for a rectangular room according to [AB79, Appendix A] in more detail. If you find any mistakes, please contact us! We are happy to improve this document. The following conventions for the Fourier transform are used in this document: The Fourier transform of the sound pressure px, t with respect to time t reads P x, ω = F t {px, t} = and the Fourier transform with respect to space x reads P k x, t = F x {px, t} = px, te jωt dt 1 px, te jkxx dx, where ω denotes the angular frequency and k x denotes the wavenumber in x-direction. A monochromatic plane wave in 3D space is then written as px, t = e jkx e jωt 3 with field point x = x, y, z T, and propagates in direction of the wavenumber vector k = k x, k y, k z T, k = k. Please note that these conventions differ from the ones used in the original paper [AB79] where for both spatial and temporal Fourier transform the opposite sign in the exponential function is used.
3 Derivation The inhomogeneous Helmholtz equation with an exciting source in x 0 = x 0, y 0, z 0 T is ω P x, x 0, ω + P x, x0, ω = δx x 0, 4 c where c denotes the speed of sound, exhibits for a rectangular room with rigid walls Neumann boundary conditions and the volume V = L x L y L z the solution P x, x 0, ω = 1 V m= ψ m xψ m x 0 k m k 5 with mx π my π mz π ψ m x = cos x cos y cos z, 6 L x L y L z m x m = m y 7 m z and the modal wave vector k m = m xπ L x m yπ L y m zπ L z = k x k y k z, k m = k m. 8 Please note that source and receiver point can be interchanged in eq. 4. The inverse formulation can be found in [AB79, eq. A1] but the formulation here is preferred for clarity. Eq. 5 shall now be converted to an expression containing exponential functions to arrive at [AB79, eq. A5], which is done with the exponential expansion of the cosine cosx = 1 e jx + e jx. 9 3
4 Expanding the numerator product of eq. 5 and reordering gives mx π my π mz π ψ m xψ m x 0 = cos x cos y cos z... L x L y L z mx π my π mz π cos x 0 cos y 0 cos z 0 L x L y L z = 1 e jkxx + e jkxx 1 e jkxx 0 + e jkxx e jkyy + e jkyy 1 e jkyy 0 + e jkyy e jkzz + e jkzz 1 e jkzz 0 + e jkzz = 1 e jkxx+x0 + e jkxx x0 + e jkxx x0 + e jkxx+x e jkyy+y0 + e jkyy y0 + e jkyy y0 + e jkyy+y e jkzz+z 0 + e jkzz z 0 + e jkzz z 0 + e jkzz+z 0 1 = cosk xx + x 0 + cosk x x x 0... cosk y y + y 0 + cosk y y y 0... cosk z z + z 0 + cosk z z z 0 13 = cosk xx + x 0 + cosk x x x 0... cosk y y + y 0 cosk z z + z 0 + cosk y y + y 0 cosk z z z cosk y y y 0 cosk z z + z 0 + cosk y y y 0 cosk z z z 0 14 = 1 8 cosk xx + x 0 cosk y y + y 0 cosk z z + z
5 Looking only at the first out of eight completely expanded triple cosine terms yields = 1 cosk x x + x 0 cosk y y + y 0 cosk z z + z 0 16 e jkxx+x0 + e jkxx+x 0 1 e jkyy+y0 + e jkyy+y 0 1 e jkzz+z0 + e jkzz+z 0 17 = 1 e jkxx+x0 e jkyy+y0 e jkzz+z0 + e jkxx+x0 e jkyy+y0 e jkzz+z e jkxx+x0 e jkyy+y0 e jkzz+z0 + e jkxx+x0 e jkyy+y0 e jkzz+z e jkxx+x 0 e jkyy+y 0 e jkzz+z 0 + e jkxx+x 0 e jkyy+y 0 e jkzz+z e jkxx+x 0 e jkyy+y 0 e jkzz+z 0 + e jkxx+x 0 e jkyy+y 0 e jkzz+z Each product of three exponential functions can be written as one exponential function with a longer exponent. This exponent then contains a scalar product of a certain vector k m with all three vector components running from to, cf. eq. 5 together with 8 and a vector for field and source points R p = x ± x 0, y ± y 0, z ± z 0, e.g. e jkxx+x 0+k yy+y 0 +k zz+z 0 = e km Rp. 19 There obviously exist 8 different vectors R p if considering all ±-combinations and the 8 combinations of exponential functions from eq. 18 arise for each R p. Considering now all terms arising from the 8 triple cosines in eq. 15, there are 8 vectors k m resulting from all possible combinations of the components and their opposite numbers because the indices in eq. 5 run from to, and from each vector k m, the same! exponential functions as in eq. 18 result. Thus, the sum in eq. 5 contains each exponential function 8 times, which cancels the 1 8 from eq. 15 the factor 1 8 from eq. 18 remains, however and this results in [AB79, eq. A5] P x, x 0, ω = 1 8V m= p=1 e jkm Rp k m k. 0 The exponent of the exponential function in eq. 0 might as well be negative, as this only changes the order in which the terms are summed up: P x, x 0, ω = 1 8V m= p=1 e jkm Rp k m k. 1 This was done to be able to directly apply the Fourier conventions used throughout this derivation, cf. eq. 1 and. Because of the sifting property of the Dirac function for time instant t 0 [GRS01] ftδt t 0 dt = ft 0, 5
6 eq. 0 can be rewritten as [AB79, eq. A8] P x, x 0, ω = 1 8V p=1 e jξ Rp ξ k m= δξ k m d 3 ξ 3 with ξ = ξ x, ξ y, ξ z T. Due to the Fourier series of the Dirac comb with sampling period T [Wik18a] δt nt = 1 π jn e T t, 4 T n= n= the following expression can be used for each component of m in eq. 3, here exemplarily shown for m x : δξ x m xπ = L x e jmxlxξx. 5 L x π m x= m x= When considering all 3 dimensions of m, this leads to m= δξ k m = L xl y L z π 3 = V π 3 m= = V π 3 m= m= e jmxlxξx e jmylyξy e jmzlzξz 6 e jmxlxξx+mylyξy+mzlzξz 7 e jrm ξ, 8 where the exponent holds the scalar product of the vector ξ with m x L x R m = m y L y. 9 m z L z Substituting eq. 8 in 3 cancels the room volume V and leads to P x, x 0, ω = 1 p=1 m= e jξrp Rm ξ k d 3 ξ. 30 The sign in the exponent can be inverted again, because adding instead of subtracting R m inside the brackets only changes the order of summation cf. [AB79, eq. A10] P x, x 0, ω = 1 p=1 m= e jξrp+rm ξ k d 3 ξ. 31 Again, this step was done to better agree with the Fourier convention used here. It can be shown, that the integral in eq. 31 is a plane wave decomposition PWD of a point source at the point R p + R m. 6
7 Plane wave decomposition of a point source Cf. e.g. [Wik18b], here rewritten according to the Fourier conventions used here see eq. 1 and, [Wik18b] uses the same convention as [AB79]. The solution of the three-dimensional inhomogeneous Helmholtz equation for the field ϕx ϕx + k 0ϕx = δx 3 is given by a point source ϕr = e jk 0r 33 4πr with wavenumber k 0 and distance to the origin r. Assuming that the Fourier transform of ϕx, denoted as ˆϕk, exists, the inverse transform is given by ϕx = 1 ˆϕke jkx dk. 34 When substituting the complete inverse Fourier transform in eq. 3, the prefactor jk x,y,z appears for each derivation. Note further that the Dirac impulse is the inverse Fourier transform of 1: Thus, eq. 3 can be rewritten as δx = 1 1 e jkx dk k x k y k z + k 0 ˆϕke jkx dk = 1 e jkx dk. 36 The integrands can be considered in isolation or the Fourier transforms are performed, which yields Substituting eq. 38 in 34 gives kx ky kz +k0 ˆϕk = 1 37 }{{} = k 1 ˆϕk = k k0. 38 ϕx = 1 e jkx k k 0 dk = e jk 0r 4πr. 39 7
8 Eq. 39 expressed with the nomenclature used so far reads e jk R 4π R = 1 e jξr ξ k d ξ 40 where R denotes the distance between source and field point. Consequently, the integrals in equation eq. 31 disappear and it can be rewritten as [AB79, eq. A13] P x, x 0, ω = p=1 m= e jk Rp+Rm 4π R p + R m. 41 To obtain the desired time domain solution, the inverse Fourier transform with respect to time is applied px, x 0, t = 1 π = p=1 m= p=1 m= e jk Rp+Rm 4π R p + R m e jωt dω π R p + R m π With the Fourier transform of a Dirac shifted by time period τ [GRS01] Rp+Rm jω e c e jωt dω. 43 δt τe jωt dt = e jωτ, 44 the result now not depending anymore on the choice of Fourier transform conventions is [AB79, eq. A14] px, x 0, t = p=1 m= δt Rp+Rm c 4π R p + R m. 45 Obviously, the sound field inside a rectangular room with rigid walls is a superposition of delayed and weighted Dirac pulses. The delay and weighting factors for each Dirac depend on the speed of sound and the distance between the image sources and the field point, which is given by R p + R m. The index p can be reformulated to directly account for the 8 possible R p versions if summing over the vector u = u, v, w T with every component being 0 or 1 and using R u instead of R p x x 0 + ux 0 R u = y y 0 + vy z z 0 + wz 0 8
9 Now, the result is [LJ08, eq. 5] px, x 0, t = u=0 m= δt Ru+Rm c 4π R u + R m where the sum over the vector u equals a triple summation: = u=0 u=0 v=0 w=0, References [AB79] [GRS01] [LJ08] Jont B. Allen and David A. Berkley. Image method for efficiently simulating small-room acoustics. J. Acoust. Soc. Am., 654: , Bernd Girod, Rudolf Rabenstein, and Alexander Stenger. Signals and Systems. Wiley, Chichester, 001. Eric A Lehmann and Anders M Johansson. Prediction of energy decay in room impulse responses simulated with an image-source model. J. Acoust. Soc. Am., 141:69 77, 008. [Wik18a] Wikipedia.org. Dirac comb, read [Wik18b] Wikiversity.org. Waves in composites and metamaterials/point sources and EM vector potentials, in_composites_and_metamaterials/point_sources_and_em_vector_ potentials, read
LECTURE 4 WAVE PACKETS
LECTURE 4 WAVE PACKETS. Comparison between QM and Classical Electrons Classical physics (particle) Quantum mechanics (wave) electron is a point particle electron is wavelike motion described by * * F ma
More informationBENG 221 Mathematical Methods in Bioengineering. Fall 2017 Midterm
BENG Mathematical Methods in Bioengineering Fall 07 Midterm NAME: Open book, open notes. 80 minutes limit (end of class). No communication other than with instructor and TAs. No computers or internet,
More informationSound Propagation through Media. Nachiketa Tiwari Indian Institute of Technology Kanpur
Sound Propagation through Media Nachiketa Tiwari Indian Institute of Technology Kanpur LECTURE-17 HARMONIC PLANE WAVES Introduction In this lecture, we discuss propagation of 1-D planar waves as they travel
More informationCartesian Coordinates
Cartesian Coordinates Daniel S. Weile Department of Electrical and Computer Engineering University of Delaware ELEG 648 Cartesian Coordinates Outline Outline Separation of Variables Away from sources,
More informationWave Phenomena Physics 15c. Lecture 11 Dispersion
Wave Phenomena Physics 15c Lecture 11 Dispersion What We Did Last Time Defined Fourier transform f (t) = F(ω)e iωt dω F(ω) = 1 2π f(t) and F(w) represent a function in time and frequency domains Analyzed
More informationData-based Binaural Synthesis Including Rotational and Translatory Head-Movements
Data-based Binaural Synthesis Including Rotational and Translatory Head-Movements Frank Schultz 1 and Sascha Spors 1 1 Institute of Communications Engineering, Universität Rostock, R.-Wagner-Str. 31 (H8),
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu.6 Signal Processing: Continuous and Discrete Fall 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS
More informationENSC327 Communications Systems 2: Fourier Representations. School of Engineering Science Simon Fraser University
ENSC37 Communications Systems : Fourier Representations School o Engineering Science Simon Fraser University Outline Chap..5: Signal Classiications Fourier Transorm Dirac Delta Function Unit Impulse Fourier
More informationOPAC102. The Acoustic Wave Equation
OPAC102 The Acoustic Wave Equation Acoustic waves in fluid Acoustic waves constitute one kind of pressure fluctuation that can exist in a compressible fluid. The restoring forces responsible for propagating
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 2.6 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS
More informationLecture 3: Acoustics
CSC 83060: Speech & Audio Understanding Lecture 3: Acoustics Michael Mandel mim@sci.brooklyn.cuny.edu CUNY Graduate Center, Computer Science Program http://mr-pc.org/t/csc83060 With much content from Dan
More informationJun-ichi Takada
Wave Theory II Numerical Simulation of Waves 2) Representation of Wave Function by Using Green Function I) Green s Theorem and Free Space Green Function Jun-ichi Takada takada@ide.titech.ac.jp)) This lecture
More informationA.1 THE SAMPLED TIME DOMAIN AND THE Z TRANSFORM. 0 δ(t)dt = 1, (A.1) δ(t)dt =
APPENDIX A THE Z TRANSFORM One of the most useful techniques in engineering or scientific analysis is transforming a problem from the time domain to the frequency domain ( 3). Using a Fourier or Laplace
More informationWave Phenomena Physics 15c. Lecture 10 Fourier Transform
Wave Phenomena Physics 15c Lecture 10 Fourier ransform What We Did Last ime Reflection of mechanical waves Similar to reflection of electromagnetic waves Mechanical impedance is defined by For transverse/longitudinal
More informationWave Phenomena Physics 15c
Wave Phenomena Physics 5c Lecture Fourier Analysis (H&L Sections 3. 4) (Georgi Chapter ) What We Did Last ime Studied reflection of mechanical waves Similar to reflection of electromagnetic waves Mechanical
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 2.161 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Massachusetts
More informationIntroduction to Seismology
1.510 Introduction to Seismology Lecture 5 Feb., 005 1 Introduction At previous lectures, we derived the equation of motion (λ + µ) ( u(x, t)) µ ( u(x, t)) = ρ u(x, t) (1) t This equation of motion can
More informationIntroduction to Acoustics Exercises
. 361-1-3291 Introduction to Acoustics Exercises 1 Fundamentals of acoustics 1. Show the effect of temperature on acoustic pressure. Hint: use the equation of state and the equation of state at equilibrium.
More informationSEISMIC WAVE PROPAGATION. Lecture 2: Fourier Analysis
SEISMIC WAVE PROPAGATION Lecture 2: Fourier Analysis Fourier Series & Fourier Transforms Fourier Series Review of trigonometric identities Analysing the square wave Fourier Transform Transforms of some
More informationDSP-I DSP-I DSP-I DSP-I
NOTES FOR 8-79 LECTURES 3 and 4 Introduction to Discrete-Time Fourier Transforms (DTFTs Distributed: September 8, 2005 Notes: This handout contains in brief outline form the lecture notes used for 8-79
More informationHorizontal Local Sound Field Propagation Based on Sound Source Dimension Mismatch
Journal of Information Hiding and Multimedia Signal Processing c 7 ISSN 73-4 Ubiquitous International Volume 8, Number 5, September 7 Horizontal Local Sound Field Propagation Based on Sound Source Dimension
More informationSignals & Systems. Lecture 5 Continuous-Time Fourier Transform. Alp Ertürk
Signals & Systems Lecture 5 Continuous-Time Fourier Transform Alp Ertürk alp.erturk@kocaeli.edu.tr Fourier Series Representation of Continuous-Time Periodic Signals Synthesis equation: x t = a k e jkω
More informationPROPAGATION PHASE REPRESENTATION IN 3D SPACE USING POLES AND ZEROS IN COMPLEX FREQUENCY PLANE. ( address of lead author)
ICSV14 Cairns Australia 9-12 July, 2007 PROPAGATION PHASE REPRESENTATION IN 3D SPACE USING POLES AND ZEROS IN COMPLEX FREQUENCY PLANE Yoshinori Takahashi 1, Mikio Tohyama 2, and Kazunori Miyoshi 1 1 Kogakuin
More informationSupporting Information
Supporting Information A: Calculation of radial distribution functions To get an effective propagator in one dimension, we first transform 1) into spherical coordinates: x a = ρ sin θ cos φ, y = ρ sin
More informationHomework 4. May An LTI system has an input, x(t) and output y(t) related through the equation y(t) = t e (t t ) x(t 2)dt
Homework 4 May 2017 1. An LTI system has an input, x(t) and output y(t) related through the equation y(t) = t e (t t ) x(t 2)dt Determine the impulse response of the system. Rewriting as y(t) = t e (t
More informationIntroduction to PML in time domain
Introduction to PML in time domain Alexander Thomann Introduction to PML in time domain - Alexander Thomann p.1 Overview 1 Introduction 2 PML in one dimension Classical absorbing layers One-dimensional
More informationLinear second-order differential equations with constant coefficients and nonzero right-hand side
Linear second-order differential equations with constant coefficients and nonzero right-hand side We return to the damped, driven simple harmonic oscillator d 2 y dy + 2b dt2 dt + ω2 0y = F sin ωt We note
More informationX. Chen More on Sampling
X. Chen More on Sampling 9 More on Sampling 9.1 Notations denotes the sampling time in second. Ω s = 2π/ and Ω s /2 are, respectively, the sampling frequency and Nyquist frequency in rad/sec. Ω and ω denote,
More informationECE Spring Prof. David R. Jackson ECE Dept. Notes 26
ECE 6345 Spring 05 Prof. David R. Jackson ECE Dept. Notes 6 Overview In this set of notes we use the spectral-domain method to find the mutual impedance between two rectangular patch ennas. z Geometry
More informationThe Helmholtz theorem at last!
Problem. The Helmholtz theorem at last! Recall in class the Helmholtz theorem that says that if if E =0 then E can be written as E = φ () if B =0 then B can be written as B = A (2) (a) Let n be a unit
More informationECE 598: The Speech Chain. Lecture 5: Room Acoustics; Filters
ECE 598: The Speech Chain Lecture 5: Room Acoustics; Filters Today Room = A Source of Echoes Echo = Delayed, Scaled Copy Addition and Subtraction of Scaled Cosines Frequency Response Impulse Response Filter
More informationPhysics 250 Green s functions for ordinary differential equations
Physics 25 Green s functions for ordinary differential equations Peter Young November 25, 27 Homogeneous Equations We have already discussed second order linear homogeneous differential equations, which
More informationENSC327 Communications Systems 2: Fourier Representations. Jie Liang School of Engineering Science Simon Fraser University
ENSC327 Communications Systems 2: Fourier Representations Jie Liang School of Engineering Science Simon Fraser University 1 Outline Chap 2.1 2.5: Signal Classifications Fourier Transform Dirac Delta Function
More informationSeries FOURIER SERIES. Graham S McDonald. A self-contained Tutorial Module for learning the technique of Fourier series analysis
Series FOURIER SERIES Graham S McDonald A self-contained Tutorial Module for learning the technique of Fourier series analysis Table of contents Begin Tutorial c 24 g.s.mcdonald@salford.ac.uk 1. Theory
More informationTHE WAVE EQUATION. The Wave Equation in the Time Domain
THE WAVE EQUATION Disturbances of materials can be described as a waveform propagating through a medium (the material itself) as particles either push or shear neighboring particles from their equilibrium
More informationMicrophone-Array Signal Processing
Microphone-Array Signal Processing, c Apolinárioi & Campos p. 1/27 Microphone-Array Signal Processing José A. Apolinário Jr. and Marcello L. R. de Campos {apolin},{mcampos}@ieee.org IME Lab. Processamento
More informationImproved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method
Center for Turbulence Research Annual Research Briefs 2006 313 Improved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method By Y. Khalighi AND D. J. Bodony 1. Motivation
More information6 Non-homogeneous Heat Problems
6 Non-homogeneous Heat Problems Up to this point all the problems we have considered for the heat or wave equation we what we call homogeneous problems. This means that for an interval < x < l the problems
More informationDifferential equations. Background mathematics review
Differential equations Background mathematics review David Miller Differential equations First-order differential equations Background mathematics review David Miller Differential equations in one variable
More informationPhysics 506 Winter 2004
Physics 506 Winter 004 G. Raithel January 6, 004 Disclaimer: The purpose of these notes is to provide you with a general list of topics that were covered in class. The notes are not a substitute for reading
More informationReview of Fundamentals of Digital Signal Processing
Chapter 2 Review of Fundamentals of Digital Signal Processing 2.1 (a) This system is not linear (the constant term makes it non linear) but is shift-invariant (b) This system is linear but not shift-invariant
More information1 Signals and systems
978--52-5688-4 - Introduction to Orthogonal Transforms: With Applications in Data Processing and Analysis Signals and systems In the first two chapters we will consider some basic concepts and ideas as
More informationECE Spring Prof. David R. Jackson ECE Dept. Notes 15
ECE 6341 Spring 216 Prof. David R. Jackson ECE Dept. Notes 15 1 Arbitrary Line Current TM : A (, ) Introduce Fourier Transform: I I + ( k ) jk = I e d x y 1 I = I ( k ) jk e dk 2π 2 Arbitrary Line Current
More informationECE Spring Prof. David R. Jackson ECE Dept. Notes 32
ECE 6345 Spring 215 Prof. David R. Jackson ECE Dept. Notes 32 1 Overview In this set of notes we extend the spectral-domain method to analyze infinite periodic structures. Two typical examples of infinite
More informationCHAPTER-5 WAVEPACKETS DESCRIBING the MOTION of a FREE PARTICLE in the context of the WAVE-PARTICLE DUALITY hypothesis
Lecture Notes PH 4/5 ECE 598 A. La Rosa INTRODUCTION TO QUANTUM MECHANICS CHAPTER-5 WAVEPACKETS DESCRIBING the MOTION of a FREE PARTICLE in the context of the WAVE-PARTICLE DUALITY hypothesis 5. Spectral
More informationRadio Propagation Channels Exercise 2 with solutions. Polarization / Wave Vector
/8 Polarization / Wave Vector Assume the following three magnetic fields of homogeneous, plane waves H (t) H A cos (ωt kz) e x H A sin (ωt kz) e y () H 2 (t) H A cos (ωt kz) e x + H A sin (ωt kz) e y (2)
More informationSpectra (2A) Young Won Lim 11/8/12
Spectra (A) /8/ Copyright (c) 0 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version. or any later version
More informationQuantum Mechanics for Scientists and Engineers. David Miller
Quantum Mechanics for Scientists and Engineers David Miller Background mathematics 7 Differential equations Differential equations in one variable A differential equation involves derivatives of some function
More informationINFLUENCE OF ROOM GEOMETRY ON MODAL DENSITY, SPATIAL DISTRIBUTION OF MODES AND THEIR DAMPING M. MEISSNER
ARCHIVES OF ACOUSTICS 30, 4 (Supplement), 203 208 (2005) INFLUENCE OF ROOM GEOMETRY ON MODAL DENSITY, SPATIAL DISTRIBUTION OF MODES AND THEIR DAMPING M. MEISSNER Institute of Fundamental Technological
More informationCourse Secretary: Christine Berber O3.095, phone x-6351,
IMPRS: Ultrafast Source Technologies Franz X. Kärtner (Umit Demirbas) & Thorsten Uphues, Bldg. 99, O3.097 & Room 6/3 Email & phone: franz.kaertner@cfel.de, 040 8998 6350 thorsten.uphues@cfel.de, 040 8998
More information26. The Fourier Transform in optics
26. The Fourier Transform in optics What is the Fourier Transform? Anharmonic waves The spectrum of a light wave Fourier transform of an exponential The Dirac delta function The Fourier transform of e
More informationEvidence that x-rays are wave-like
Evidence that x-rays are wave-like After their discovery in 1895 by Roentgen, their spectrum (including characteristic x-rays) was probed and their penetrating ability was exploited, but it was difficult
More informationSignal and systems. Linear Systems. Luigi Palopoli. Signal and systems p. 1/5
Signal and systems p. 1/5 Signal and systems Linear Systems Luigi Palopoli palopoli@dit.unitn.it Wrap-Up Signal and systems p. 2/5 Signal and systems p. 3/5 Fourier Series We have see that is a signal
More informationWavenumber-Frequency Space
Wavenumber-Frequency Space ECE 6279: Spatial Array Processing Spring 2011 Lecture 3 Prof. Aaron D. Lanterman School of Electrical & Computer Engineering Georgia Institute of Technology AL: 404-385-2548
More informationMandl and Shaw reading assignments
Mandl and Shaw reading assignments Chapter 2 Lagrangian Field Theory 2.1 Relativistic notation 2.2 Classical Lagrangian field theory 2.3 Quantized Lagrangian field theory 2.4 Symmetries and conservation
More informationHomework 7 Solution EE235, Spring Find the Fourier transform of the following signals using tables: te t u(t) h(t) = sin(2πt)e t u(t) (2)
Homework 7 Solution EE35, Spring. Find the Fourier transform of the following signals using tables: (a) te t u(t) h(t) H(jω) te t u(t) ( + jω) (b) sin(πt)e t u(t) h(t) sin(πt)e t u(t) () h(t) ( ejπt e
More informationThe Casimir Effect. Notes by Asaf Szulc
The Casimir Effect Notes by Asaf Szulc Introduction The Casimir effect describes the attraction between two perfectly conducting and uncharged parallel plates confined in vacuum. It was the Dutch physicist
More informationThe (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 informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationIB Paper 6: Signal and Data Analysis
IB Paper 6: Signal and Data Analysis Handout 5: Sampling Theory S Godsill Signal Processing and Communications Group, Engineering Department, Cambridge, UK Lent 2015 1 / 85 Sampling and Aliasing All of
More informationENGIN 211, Engineering Math. Fourier Series and Transform
ENGIN 11, Engineering Math Fourier Series and ransform 1 Periodic Functions and Harmonics f(t) Period: a a+ t Frequency: f = 1 Angular velocity (or angular frequency): ω = ππ = π Such a periodic function
More informationPropagating Wave (1B)
Wave (1B) 3-D Wave Copyright (c) 2011 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
More informationLecture 1 Signals. notation and meaning. common signals. size of a signal. qualitative properties of signals. impulsive signals
EE102 spring 2002 Handout #2 Lecture 1 Signals notation and meaning common signals size of a signal qualitative properties of signals impulsive signals 1 1 Signals a signal is a function of time, e.g.,
More informationChapter 6 THE SAMPLING PROCESS 6.1 Introduction 6.2 Fourier Transform Revisited
Chapter 6 THE SAMPLING PROCESS 6.1 Introduction 6.2 Fourier Transform Revisited Copyright c 2005 Andreas Antoniou Victoria, BC, Canada Email: aantoniou@ieee.org July 14, 2018 Frame # 1 Slide # 1 A. Antoniou
More informationGreen s functions for planarly layered media
Green s functions for planarly layered media Massachusetts Institute of Technology 6.635 lecture notes Introduction: Green s functions The Green s functions is the solution of the wave equation for a point
More informationLecture 34. Fourier Transforms
Lecture 34 Fourier Transforms In this section, we introduce the Fourier transform, a method of analyzing the frequency content of functions that are no longer τ-periodic, but which are defined over the
More informationA second look at waves
A second loo at waves ravelling waves A first loo at Amplitude Modulation (AM) Stationary and reflected waves Lossy waves: dispersion & evanescence I thin this is the MOS IMPORAN of my eight lectures,
More informationA wavenumber approach to characterizing the diffuse field conditions in reverberation rooms
PROCEEDINGS of the 22 nd International Congress on Acoustics Isotropy and Diffuseness in Room Acoustics: Paper ICA2016-578 A wavenumber approach to characterizing the diffuse field conditions in reverberation
More information(TRAVELLING) 1D WAVES. 1. Transversal & Longitudinal Waves
(TRAVELLING) 1D WAVES 1. Transversal & Longitudinal Waves Objectives After studying this chapter you should be able to: Derive 1D wave equation for transversal and longitudinal Relate propagation speed
More informationVibrations of string. Henna Tahvanainen. November 8, ELEC-E5610 Acoustics and the Physics of Sound, Lecture 4
Vibrations of string EEC-E5610 Acoustics and the Physics of Sound, ecture 4 Henna Tahvanainen Department of Signal Processing and Acoustics Aalto University School of Electrical Engineering November 8,
More informationEE 224 Signals and Systems I Review 1/10
EE 224 Signals and Systems I Review 1/10 Class Contents Signals and Systems Continuous-Time and Discrete-Time Time-Domain and Frequency Domain (all these dimensions are tightly coupled) SIGNALS SYSTEMS
More informationUnderstand the basic principles of spectroscopy using selection rules and the energy levels. Derive Hund s Rule from the symmetrization postulate.
CHEM 5314: Advanced Physical Chemistry Overall Goals: Use quantum mechanics to understand that molecules have quantized translational, rotational, vibrational, and electronic energy levels. In a large
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 12.510 Introduction to Seismology Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 12.510 Introduction
More informationECE-700 Review. Phil Schniter. January 5, x c (t)e jωt dt, x[n]z n, Denoting a transform pair by x[n] X(z), some useful properties are
ECE-7 Review Phil Schniter January 5, 7 ransforms Using x c (t) to denote a continuous-time signal at time t R, Laplace ransform: X c (s) x c (t)e st dt, s C Continuous-ime Fourier ransform (CF): ote that:
More informationEE 3054: Signals, Systems, and Transforms Summer It is observed of some continuous-time LTI system that the input signal.
EE 34: Signals, Systems, and Transforms Summer 7 Test No notes, closed book. Show your work. Simplify your answers. 3. It is observed of some continuous-time LTI system that the input signal = 3 u(t) produces
More information2 z = ±b and the bottom is the plane y = h. The free surface is located at y = (x; z; t), which is unknown. We assume irrotational flow and incompress
1 I-campus project School-wide Program on Fluid Mechanics Modules on Waves in fluids T. R. Akylas & C. C. Mei CHAPTER SIX FORCED DISPERSIVE WAVES ALONG A NARROW CHANNEL Linear surface gravity w aves propagating
More informationELEN 4810 Midterm Exam
ELEN 4810 Midterm Exam Wednesday, October 26, 2016, 10:10-11:25 AM. One sheet of handwritten notes is allowed. No electronics of any kind are allowed. Please record your answers in the exam booklet. Raise
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 2.161 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Massachusetts
More informationLine Spectra and their Applications
In [ ]: cd matlab pwd Line Spectra and their Applications Scope and Background Reading This session concludes our introduction to Fourier Series. Last time (http://nbviewer.jupyter.org/github/cpjobling/eg-47-
More informationPartial differential equations (ACM30220)
(ACM3. A pot on a stove has a handle of length that can be modelled as a rod with diffusion constant D. The equation for the temperature in the rod is u t Du xx < x
More informationLinear Theory of Evolution to an Unstable State
Chapter 2 Linear Theory of Evolution to an Unstable State c 2012 by William Klein, Harvey Gould, and Jan Tobochnik 1 October 2012 2.1 Introduction The simple theory of nucleation that we introduced in
More informationCS-9645 Introduction to Computer Vision Techniques Winter 2018
Table of Contents Spectral Analysis...1 Special Functions... 1 Properties of Dirac-delta Functions...1 Derivatives of the Dirac-delta Function... 2 General Dirac-delta Functions...2 Harmonic Analysis...
More informationFourier transform representation of CT aperiodic signals Section 4.1
Fourier transform representation of CT aperiodic signals Section 4. A large class of aperiodic CT signals can be represented by the CT Fourier transform (CTFT). The (CT) Fourier transform (or spectrum)
More information1 Separation of Variables
Jim ambers ENERGY 281 Spring Quarter 27-8 ecture 2 Notes 1 Separation of Variables In the previous lecture, we learned how to derive a PDE that describes fluid flow. Now, we will learn a number of analytical
More informationFourier 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 informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
More informationTransfer-matrix-based approach for an eigenvalue problem of a coupled rectangular cavity
ransfer-matrix-based approach for an eigenvalue problem of a coupled rectangular cavity Hiroyuki IWAMOO ; Nobuo ANAKA Seikei University, Japan okyo Metropolitan University, Japan ABSRAC his study concerns
More informationChapter 5 Wave-Optics Analysis of Coherent Optical Systems
Chapter 5 Wave-Optics Analysis of Coherent Optical Systems January 5, 2016 Chapter 5 Wave-Optics Analysis of Coherent Optical Systems Contents: 5.1 A thin lens as a phase transformation 5.2 Fourier transforming
More informationFrequency Domain Representations of Sampled and Wrapped Signals
Frequency Domain Representations of Sampled and Wrapped Signals Peter Kabal Department of Electrical & Computer Engineering McGill University Montreal, Canada v1.5 March 2011 c 2011 Peter Kabal 2011/03/11
More informationCS711008Z Algorithm Design and Analysis
CS711008Z Algorithm Design and Analysis Lecture 5 FFT and Divide and Conquer Dongbo Bu Institute of Computing Technology Chinese Academy of Sciences, Beijing, China 1 / 56 Outline DFT: evaluate a polynomial
More informationSmoluchowski Diffusion Equation
Chapter 4 Smoluchowski Diffusion Equation Contents 4. Derivation of the Smoluchoswki Diffusion Equation for Potential Fields 64 4.2 One-DimensionalDiffusoninaLinearPotential... 67 4.2. Diffusion in an
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 informationSolution of ODEs using Laplace Transforms. Process Dynamics and Control
Solution of ODEs using Laplace Transforms Process Dynamics and Control 1 Linear ODEs For linear ODEs, we can solve without integrating by using Laplace transforms Integrate out time and transform to Laplace
More informationFunctional Quantization
Functional Quantization In quantum mechanics of one or several particles, we may use path integrals to calculate the transition matrix elements as out Ût out t in in D[allx i t] expis[allx i t] Ψ out allx
More informatione iωt dt and explained why δ(ω) = 0 for ω 0 but δ(0) =. A crucial property of the delta function, however, is that
Phys 531 Fourier Transforms In this handout, I will go through the derivations of some of the results I gave in class (Lecture 14, 1/11). I won t reintroduce the concepts though, so you ll want to refer
More informationS=z y. p(x, t) v. = ρ x (x,t)
LTUR 2: One-imensional Traveling Waves Main Points - xponential and sine-wave solutions to the one-dimensional wave equation. - The distributed compressibility and mass in acoustic plane waves are analogous
More informationSpatial Smoothing and Broadband Beamforming. Bhaskar D Rao University of California, San Diego
Spatial Smoothing and Broadband Beamforming Bhaskar D Rao University of California, San Diego Email: brao@ucsd.edu Reference Books and Papers 1. Optimum Array Processing, H. L. Van Trees 2. Stoica, P.,
More informationDevelopment of a Generalized Corcos model for the prediction of turbulent boundary layer induced noise
Development of a Generalized Corcos model for the prediction of turbulent boundary layer induced noise Anna Caiazzo Roberto D Amico, Wim Desmet KU Leuven - Noise & Vibration Research Group Departement
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationCHAPTER-5 WAVEPACKETS DESCRIBING the MOTION of a FREE PARTICLE in the context of the WAVE-PARTICLE DUALITY hypothesis
Lecture Notes PH 4/5 ECE 598 A. La Rosa INTRODUCTION TO QUANTUM MECHANICS CHAPTER-5 WAVEPACKETS DESCRIBING the MOTION of a FREE PARTICLE in the context of the WAVE-PARTICLE DUALITY hypothesis 5. Spectral
More information