Topics. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2004 Lecture 4 2D Fourier Transforms

Size: px
Start display at page:

Download "Topics. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2004 Lecture 4 2D Fourier Transforms"

Transcription

1 Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2004 Lecture 4 2D Fourier Transforms Topics 1. 2D Signal Representations 2. 2D Fourier Transform 3. Transform Pairs 4. FT Properties 1

2 2D Signal a 0 0 b a b = + c d c 0 0 d Image Decomposition 1 0 a b 0 1 a b = + c d c d g[m,n] = aδ[m,n] + bδ[m,n 1] + cδ[m 1,n] + dδ[m 1,n 1] = g[k,l] δ[m k,n l] k= 0 l= 0 k= 0 l= = c k,l b k,l [m,n] 2

3 Orthonormal Basis Functions Discrete b k,l,b k, l = b k,l [m,n] b Continuous m= n= = δ[k k,l l ] b kx,b k x, k y = b kx k, (x, yb k x, = δ(k x k x k y l [m,n] k y (x, y dxdy Are these orthonormal? Example

4 Example Are these orthonormal? 1/2 1/2 1/2 1/2 1/2 1/2-1/2-1/2 1/2-1/2 1/2-1/2 1/2-1/2-1/2 1/2 Discrete Expansion Coefficients The discrete expansion is g[m,n] = c k,l k= l= b k,l [m,n] If the basis functions b k,l [m,n] are orthonormal then c k,l = b k,l,g = = m= n= m= n= b k,l = c b k,l [m,n]g[m,n] [m,n] c k, l b k, l [m,n] k = l = k, l b k,l k = l = m= n= [m,n] b k, l [m,n] = c k, l δ[k k,l l ] k = l = = c k,l 4

5 Continuous Expansion Coefficients The continuous expansion is g(x, y = ( c k x b kx (x,ydk x dk y If the basis functions b kx (x,y are orthonormal then ( = b kx,g = b kx c k x (x,yg(x, ydxdy ( = b kx (x, y c k x, k b y k x, k y (x,yd k x d k y dxdy ( k y k x, k y = c k x, ( ( b kx (x, yb k x, k y (x, ydxdyd k x d k y = c δ(k x k x k y d k x d k y = c k x Separable Basis Functions Discrete b k,l [m,n] = b k [m]b l [n] e.g. δ[m k,n l] = δ[m k]δ[n l] Continuous b kx (x, y = b kx (xb ky (y e.g. δ( x x i,y y i = δ( x x i δ y y i ( 5

6 b 1 [m] = 1 0 Separable Basis Functions b 1 [n] = [1 0] b 2 [n] = [0 1] 1 0 b 1 [m] = b 1 [n] = [1 0] b 2 [m] = 0 1 b 2 [n] = [0 1] 0 0 b 2 [m] = Separable Basis Functions b 1 [n] = [1 1]/ 2 b 2 [n] = [1 1]/ 2 b 1 [m] = 1 1 1/2 1/2 1/2-1/2 2 1 b 1 [m] = 1 1 1/2 1/ /2-1/2 b 1 [n] = [1 1]/ 2 1/2 1/2 b 1 [m] = /2-1/2 b 1 [m] = b 2 [n] = [1 1]/ 2 1/2-1/2-1/2 1/2 b k,l [m,n] = b k [m] b l [n] b k [m] = exp( πmk / 2 b l [n] = exp( πnl / 2 6

7 Example x = Basis Functions Coefficients Sum Object Fourier Basis Functions Recall that for 1D the basis functions are complex exponentials b kx ( x = e j 2πk x x For 2D, we use the separable 2D functions b kx ( x, y = b kx ( xb ky y ( = e j 2πk x x e j 2πk yy = e j 2π (k xx +k y y Are they orthonormal? 7

8 Plane Waves e j 2π (k x x +k y y = cos( 2π(k x x + k y y + j sin( 2π(k x x + k y y 1 k x 2 + k y 2 1/k y 1/k x cos(2πk x x cos(2πk y y cos(2πk x x +2πk y y Plane Waves A θ 1/k y D θ B 1/k x C ΔABC ~ ΔBDC AC BC = AB BD BD = AB BC AC = θ = arctan k y k x 1 1 k x k y 1 k x k y = 1 k x 2 + k y 2 8

9 Fourier Transform 2D Fourier Transform = F[ g( x, y ] = e j 2π k x x +k y y (,g = g(x, y ( dxdy e j 2π k x x +k y y Inverse Fourier Transform g(x, y = ( dk x dk y e j 2π k x x +k y y Separable Functions ( is said to be a separable function if it can be g x, y written as g( x, y = g X ( xg Y ( y The Fourier Transform is then separable as well. = g(x,y e j 2π ( k xx +k y y dxdy = g X ( xe j 2πkxx dx g Y ( y e j 2πkyy dy = G X (k x G Y (k y Example g(x, y = Π(xΠ(y = sinc(k x sinc(k y 9

10 Example g(x, y = Π(xΠ(y = sinc(k x sinc(k y Example (sinc/rect y 1/2-1/2 1/2 x -1/2 Example (sinc/rect 10

11 11

12 12

13 g(x, y = δ(x, y = δ(xδ(y =1 g(x, y = δ(x = δ(k y!!! g(x, y =1+ e j 2πax = δ k x ( + δ(k x + aδ(k y g(x, y =1+ e j 2πay = δ(k x + δ(k x δ(k y a 13

14 g(x, y = cos(2π(ax + by = 1 2 δ(k x aδ(k y b δ(k x + aδ(k y + b = δ(k x + g(x, y =??? δ(k x + cδ(k y + δ(k x δ(k y d δ(k x aδ(k y b δ(k x + aδ(k y + b 14

15 15

16 Linearity F ag(x, y + bh(x, y Basic Properties [ ] = a + bh(k x Scaling [ ] = 1 ab G k x F g(ax,by a, k x b Shift F[ g(x a,y b ] = e j 2π (k xa +k y b Modulation [ ] = a b j 2π (xa +yb F g(x, ye Modulation Example x = * = 16

17 Convolution/Multiplication Convolution [ ] = H(k x F g(x, y h(x,y Multiplication [ ] = H(k x F g(x, yh(x, y Multiplication in one domain translates into convolution in the other domain. Convolution Example 17

18 F[ g (x, y ] = G ( k x, k y Symmetry If g(x, y is real then g(x, y = g (x, y, so G ( k x, k y = k y k x Parseval s Relations Energy is preserved 2 2 g(x, y dxdy = G(kx dkx dk y So is the inner product g,h = G,H g (x, y h(x, ydxdy = G (k x H(k x dk x dk y 18

Spin-Warp Pulse Sequence

Spin-Warp Pulse Sequence Bioengineering 28A Principles of Biomedical Imaging Fall Quarter 25 Linear Systems Lecture Spin-Warp Pulse Sequence RF G x (t G y (t k y k x K-space trajectories EPI Spiral k y k y k x k x Credit: Larry

More information

Convolution Theorem Modulation Transfer Function y(t)=g(t)*h(t)

Convolution Theorem Modulation Transfer Function y(t)=g(t)*h(t) MTF = Fourier Transform of PSF Bioengineering 28A Principles of Biomedical Imaging Fall Quarter 214 CT/Fourier Lecture 4 Bushberg et al 21 Convolution Theorem Modulation Transfer Function y(t=g(t*h(t g(t

More information

Convolution/Modulation Theorem. 2D Convolution/Multiplication. Application of Convolution Thm. Bioengineering 280A Principles of Biomedical Imaging

Convolution/Modulation Theorem. 2D Convolution/Multiplication. Application of Convolution Thm. Bioengineering 280A Principles of Biomedical Imaging Bioengineering 8A Principles of Biomedical Imaging Fall Quarter 15 CT/Fourier Lecture 4 Convolution/Modulation Theorem [ ] e j πx F{ g(x h(x } = g(u h(x udu = g(u h(x u e j πx = g(uh( e j πu du = H( dxdu

More information

Backprojection. Projections. Projections " $ & cosθ & Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2014 CT/Fourier Lecture 2

Backprojection. Projections. Projections  $ & cosθ & Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2014 CT/Fourier Lecture 2 Backprojection Bioengineering 280A Principles of Biomedical Imaging 3 0 0 0 0 3 0 0 0 0 0 3 0 0 3 0 3 Fall Quarter 2014 CT/Fourier Lecture 2 0 0 0 1 1 1 0 0 0 1 0 0 1 2 1 0 0 1 1 1 0 1 3 1 0 1 1 1 1 1

More information

Topics. Example. Modulation. [ ] = G(k x. ] = 1 2 G ( k % k x 0) G ( k + k x 0) ] = 1 2 j G ( k x

Topics. Example. Modulation. [ ] = G(k x. ] = 1 2 G ( k % k x 0) G ( k + k x 0) ] = 1 2 j G ( k x Topics Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2008 CT/Fourier Lecture 3 Modulation Modulation Transfer Function Convolution/Multiplication Revisit Projection-Slice Theorem Filtered

More information

Representation of 1D Function

Representation of 1D Function Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2005 Linear Systems Lecture 2 Representation of 1D Function From the sifting property, we can write a 1D function as g(x) = g(ξ)δ(x ξ)dξ.

More information

Topics. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2006 CT/Fourier Lecture 2

Topics. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2006 CT/Fourier Lecture 2 Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2006 CT/Fourier Lecture 2 Topics Modulation Modulation Transfer Function Convolution/Multiplication Revisit Projection-Slice Theorem Filtered

More information

Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2006 X-Rays Lecture 3. Topics

Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2006 X-Rays Lecture 3. Topics Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2006 X-Rays Lecture 3 Topics Address questions from last lecture Review image equation Linearity Delta Functions Superposition Convolution

More information

The Fourier transform, indicated by the operator F, constructs a spectrum A(k x,k y ) = F {E(x,y)} from a spatial distribution E(x,y):

The Fourier transform, indicated by the operator F, constructs a spectrum A(k x,k y ) = F {E(x,y)} from a spatial distribution E(x,y): 1 4. FOURIER TRANSFORMS (vs. 3.0) DEFINITIONS GENERAL PROPERTIES linearity conjugate functions multiplication and convolution multiplication and correlation Parseval's relation time shifting frequency

More information

Kronecker Delta Function # 1 for m = 0,n = 0. Kronecker Delta Function

Kronecker Delta Function # 1 for m = 0,n = 0. Kronecker Delta Function Bioegieerig 280A Priciples of Biomedical Imagig Fall Quarter 202 Sigals ad Fourier Trasforms Sigals ad Images Discrete-time/space sigałimage: cotiuous valued fuctio with a discrete time/space ide, deoted

More information

MIT 2.71/2.710 Optics 10/31/05 wk9-a-1. The spatial frequency domain

MIT 2.71/2.710 Optics 10/31/05 wk9-a-1. The spatial frequency domain 10/31/05 wk9-a-1 The spatial frequency domain Recall: plane wave propagation x path delay increases linearly with x λ z=0 θ E 0 x exp i2π sinθ + λ z i2π cosθ λ z plane of observation 10/31/05 wk9-a-2 Spatial

More information

BME I5000: Biomedical Imaging

BME I5000: Biomedical Imaging 1 BME I5000: Biomedical Imaging Lecture FT Fourier Transform Review Lucas C. Parra, parra@ccny.cuny.edu Blackboard: http://cityonline.ccny.cuny.edu/ Fourier Transform (1D) The Fourier Transform (FT) is

More information

Fourier Sampling. Fourier Sampling. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2005 Linear Systems Lecture 3.

Fourier Sampling. Fourier Sampling. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2005 Linear Systems Lecture 3. Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2005 Linear Sstems Lecture 3 Fourier Sampling F Instead of sampling the signal, we sample its Fourier Transform Sample??? F -1 Fourier

More information

Advanced Training Course on FPGA Design and VHDL for Hardware Simulation and Synthesis. 26 October - 20 November, 2009

Advanced Training Course on FPGA Design and VHDL for Hardware Simulation and Synthesis. 26 October - 20 November, 2009 2065-33 Advanced Training Course on FPGA Design and VHDL for Hardware Simulation and Synthesis 26 October - 20 November, 2009 Introduction to two-dimensional digital signal processing Fabio Mammano University

More information

Representation of 1D Function

Representation of 1D Function Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2005 Linear Systems Lecture 2 Representation of 1D Function From the sifting property, we can write a 1D function as g( = g(ξδ( ξdξ. To

More information

Module Two: Differential Calculus(continued) synopsis of results and problems (student copy)

Module Two: Differential Calculus(continued) synopsis of results and problems (student copy) Module Two: Differential Calculus(continued) synopsis of results and problems (student copy) Srikanth K S 1 Syllabus Taylor s and Maclaurin s theorems for function of one variable(statement only)- problems.

More information

Fourier Transform. sin(n# x)), where! = 2" / L and

Fourier Transform. sin(n# x)), where! = 2 / L and Fourier Transform Henning Stahlberg Introduction The tools provided by the Fourier transform are helpful for the analysis of 1D signals (time and frequency (or Fourier) domains), as well as 2D/3D signals

More information

Nov : Lecture 18: The Fourier Transform and its Interpretations

Nov : Lecture 18: The Fourier Transform and its Interpretations 3 Nov. 04 2005: Lecture 8: The Fourier Transform and its Interpretations Reading: Kreyszig Sections: 0.5 (pp:547 49), 0.8 (pp:557 63), 0.9 (pp:564 68), 0.0 (pp:569 75) Fourier Transforms Expansion of a

More information

Spin. Nuclear Spin Rules

Spin. Nuclear Spin Rules Spin Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 203 MRI Lecture Intrinsic angular momentum of elementary particles -- electrons, protons, neutrons. Spin is quantized. Key concept

More information

2.6 Logarithmic Functions. Inverse Functions. Question: What is the relationship between f(x) = x 2 and g(x) = x?

2.6 Logarithmic Functions. Inverse Functions. Question: What is the relationship between f(x) = x 2 and g(x) = x? Inverse Functions Question: What is the relationship between f(x) = x 3 and g(x) = 3 x? Question: What is the relationship between f(x) = x 2 and g(x) = x? Definition (One-to-One Function) A function f

More information

worked out from first principles by parameterizing the path, etc. If however C is a A path C is a simple closed path if and only if the starting point

worked out from first principles by parameterizing the path, etc. If however C is a A path C is a simple closed path if and only if the starting point III.c Green s Theorem As mentioned repeatedly, if F is not a gradient field then F dr must be worked out from first principles by parameterizing the path, etc. If however is a simple closed path in the

More information

Fourier Approach to Wave Propagation

Fourier Approach to Wave Propagation Phys 531 Lecture 15 13 October 005 Fourier Approach to Wave Propagation Last time, reviewed Fourier transform Write any function of space/time = sum of harmonic functions e i(k r ωt) Actual waves: harmonic

More information

Conformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder.

Conformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder. Lent 29 COMPLEX METHODS G. Taylor A star means optional and not necessarily harder. Conformal maps. (i) Let f(z) = az + b, with ad bc. Where in C is f conformal? cz + d (ii) Let f(z) = z +. What are the

More information

MTH739U/P: Topics in Scientific Computing Autumn 2016 Week 6

MTH739U/P: Topics in Scientific Computing Autumn 2016 Week 6 MTH739U/P: Topics in Scientific Computing Autumn 16 Week 6 4.5 Generic algorithms for non-uniform variates We have seen that sampling from a uniform distribution in [, 1] is a relatively straightforward

More information

Dr. Abdulla Eid. Section 3.8 Derivative of the inverse function and logarithms 3 Lecture. Dr. Abdulla Eid. MATHS 101: Calculus I. College of Science

Dr. Abdulla Eid. Section 3.8 Derivative of the inverse function and logarithms 3 Lecture. Dr. Abdulla Eid. MATHS 101: Calculus I. College of Science Section 3.8 Derivative of the inverse function and logarithms 3 Lecture College of Science MATHS 101: Calculus I (University of Bahrain) Logarithmic Differentiation 1 / 19 Topics 1 Inverse Functions (1

More information

Math 121A: Homework 6 solutions

Math 121A: Homework 6 solutions Math A: Homework 6 solutions. (a) The coefficients of the Fourier sine series are given by b n = π f (x) sin nx dx = x(π x) sin nx dx π = (π x) cos nx dx nπ nπ [x(π x) cos nx]π = n ( )(sin nx) dx + π n

More information

EE 438 Essential Definitions and Relations

EE 438 Essential Definitions and Relations May 2004 EE 438 Essential Definitions and Relations CT Metrics. Energy E x = x(t) 2 dt 2. Power P x = lim T 2T T / 2 T / 2 x(t) 2 dt 3. root mean squared value x rms = P x 4. Area A x = x(t) dt 5. Average

More information

Discrete Fourier Transform

Discrete Fourier Transform Last lecture I introduced the idea that any function defined on x 0,..., N 1 could be written a sum of sines and cosines. There are two different reasons why this is useful. The first is a general one,

More information

Math 115 ( ) Yum-Tong Siu 1. Derivation of the Poisson Kernel by Fourier Series and Convolution

Math 115 ( ) Yum-Tong Siu 1. Derivation of the Poisson Kernel by Fourier Series and Convolution Math 5 (006-007 Yum-Tong Siu. Derivation of the Poisson Kernel by Fourier Series and Convolution We are going to give a second derivation of the Poisson kernel by using Fourier series and convolution.

More information

Joint Distributions: Part Two 1

Joint Distributions: Part Two 1 Joint Distributions: Part Two 1 STA 256: Fall 2018 1 This slide show is an open-source document. See last slide for copyright information. 1 / 30 Overview 1 Independence 2 Conditional Distributions 3 Transformations

More information

1 Review of di erential calculus

1 Review of di erential calculus Review of di erential calculus This chapter presents the main elements of di erential calculus needed in probability theory. Often, students taking a course on probability theory have problems with concepts

More information

x 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7

x 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7 Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)

More information

Finite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product

Finite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )

More information

Suggested Solution to Assignment 7

Suggested Solution to Assignment 7 MATH 422 (25-6) partial diferential equations Suggested Solution to Assignment 7 Exercise 7.. Suppose there exists one non-constant harmonic function u in, which attains its maximum M at x. Then by the

More information

UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH

UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH College of Informatics and Electronics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MS425 SEMESTER: Autumn 25/6 MODULE TITLE: Applied Analysis DURATION OF EXAMINATION:

More information

Chapter 7: Bounded Operators in Hilbert Spaces

Chapter 7: Bounded Operators in Hilbert Spaces Chapter 7: Bounded Operators in Hilbert Spaces I-Liang Chern Department of Applied Mathematics National Chiao Tung University and Department of Mathematics National Taiwan University Fall, 2013 1 / 84

More information

Intro to harmonic analysis on groups Risi Kondor

Intro to harmonic analysis on groups Risi Kondor Risi Kondor Any (sufficiently smooth) function f on the unit circle (equivalently, any 2π periodic f ) can be decomposed into a sum of sinusoidal waves f(x) = k= c n e ikx c n = 1 2π f(x) e ikx dx 2π 0

More information

Chapter 1: Packing your Suitcase

Chapter 1: Packing your Suitcase Chapter : Packing your Suitcase Lesson.. -. a. Independent variable = distance from end of tube to the wall. Dependent variable = width of field of view. e. The equation depends on the length and diameter

More information

ECG782: Multidimensional Digital Signal Processing

ECG782: Multidimensional Digital Signal Processing Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Spring 2014 TTh 14:30-15:45 CBC C313 Lecture 05 Image Processing Basics 13/02/04 http://www.ee.unlv.edu/~b1morris/ecg782/

More information

Chapter 8 Integral Operators

Chapter 8 Integral Operators Chapter 8 Integral Operators In our development of metrics, norms, inner products, and operator theory in Chapters 1 7 we only tangentially considered topics that involved the use of Lebesgue measure,

More information

conditional cdf, conditional pdf, total probability theorem?

conditional cdf, conditional pdf, total probability theorem? 6 Multiple Random Variables 6.0 INTRODUCTION scalar vs. random variable cdf, pdf transformation of a random variable conditional cdf, conditional pdf, total probability theorem expectation of a random

More information

CS-9645 Introduction to Computer Vision Techniques Winter 2018

CS-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 information

2. A die is rolled 3 times, the probability of getting a number larger than the previous number each time is

2. A die is rolled 3 times, the probability of getting a number larger than the previous number each time is . If P(A) = x, P = 2x, P(A B) = 2, P ( A B) = 2 3, then the value of x is (A) 5 8 5 36 6 36 36 2. A die is rolled 3 times, the probability of getting a number larger than the previous number each time

More information

FFTs in Graphics and Vision. Homogenous Polynomials and Irreducible Representations

FFTs in Graphics and Vision. Homogenous Polynomials and Irreducible Representations FFTs in Graphics and Vision Homogenous Polynomials and Irreducible Representations 1 Outline The 2π Term in Assignment 1 Homogenous Polynomials Representations of Functions on the Unit-Circle Sub-Representations

More information

DERIVATIONS. Introduction to non-associative algebra. Playing havoc with the product rule? BERNARD RUSSO University of California, Irvine

DERIVATIONS. Introduction to non-associative algebra. Playing havoc with the product rule? BERNARD RUSSO University of California, Irvine DERIVATIONS Introduction to non-associative algebra OR Playing havoc with the product rule? PART VI COHOMOLOGY OF LIE ALGEBRAS BERNARD RUSSO University of California, Irvine FULLERTON COLLEGE DEPARTMENT

More information

TC412 Microwave Communications. Lecture 8 Rectangular waveguides and cavity resonator

TC412 Microwave Communications. Lecture 8 Rectangular waveguides and cavity resonator TC412 Microwave Communications Lecture 8 Rectangular waveguides and cavity resonator 1 TM waves in rectangular waveguides Finding E and H components in terms of, WG geometry, and modes. From 2 2 2 xye

More information

Properties of surfaces II: Second moment of area

Properties of surfaces II: Second moment of area Properties of surfaces II: Second moment of area Just as we have discussing first moment of an area and its relation with problems in mechanics, we will now describe second moment and product of area of

More information

Math 265 (Butler) Practice Midterm III B (Solutions)

Math 265 (Butler) Practice Midterm III B (Solutions) Math 265 (Butler) Practice Midterm III B (Solutions). Set up (but do not evaluate) an integral for the surface area of the surface f(x, y) x 2 y y over the region x, y 4. We have that the surface are is

More information

2.12: Derivatives of Exp/Log (cont d) and 2.15: Antiderivatives and Initial Value Problems

2.12: Derivatives of Exp/Log (cont d) and 2.15: Antiderivatives and Initial Value Problems 2.12: Derivatives of Exp/Log (cont d) and 2.15: Antiderivatives and Initial Value Problems Mathematics 3 Lecture 14 Dartmouth College February 03, 2010 Derivatives of the Exponential and Logarithmic Functions

More information

Formulas that must be memorized:

Formulas that must be memorized: Formulas that must be memorized: Position, Velocity, Acceleration Speed is increasing when v(t) and a(t) have the same signs. Speed is decreasing when v(t) and a(t) have different signs. Section I: Limits

More information

Math 10C - Fall Final Exam

Math 10C - Fall Final Exam Math 1C - Fall 217 - Final Exam Problem 1. Consider the function f(x, y) = 1 x 2 (y 1) 2. (i) Draw the level curve through the point P (1, 2). Find the gradient of f at the point P and draw the gradient

More information

Sec. 14.3: Partial Derivatives. All of the following are ways of representing the derivative. y dx

Sec. 14.3: Partial Derivatives. All of the following are ways of representing the derivative. y dx Math 2204 Multivariable Calc Chapter 14: Partial Derivatives I. Review from math 1225 A. First Derivative Sec. 14.3: Partial Derivatives 1. Def n : The derivative of the function f with respect to the

More information

Discrete Fourier Transform

Discrete Fourier Transform Discrete Fourier Transform DD2423 Image Analysis and Computer Vision Mårten Björkman Computational Vision and Active Perception School of Computer Science and Communication November 13, 2013 Mårten Björkman

More information

The Dirac δ-function

The Dirac δ-function The Dirac δ-function Elias Kiritsis Contents 1 Definition 2 2 δ as a limit of functions 3 3 Relation to plane waves 5 4 Fourier integrals 8 5 Fourier series on the half-line 9 6 Gaussian integrals 11 Bibliography

More information

Lecture 4 Filtering in the Frequency Domain. Lin ZHANG, PhD School of Software Engineering Tongji University Spring 2016

Lecture 4 Filtering in the Frequency Domain. Lin ZHANG, PhD School of Software Engineering Tongji University Spring 2016 Lecture 4 Filtering in the Frequency Domain Lin ZHANG, PhD School of Software Engineering Tongji University Spring 2016 Outline Background From Fourier series to Fourier transform Properties of the Fourier

More information

Lecture Notes for Math 1000

Lecture Notes for Math 1000 Lecture Notes for Math 1000 Dr. Xiang-Sheng Wang Memorial University of Newfoundland Office: HH-2016, Phone: 864-4321 Office hours: 13:00-15:00 Wednesday, 12:00-13:00 Friday Email: xswang@mun.ca Course

More information

1,3. f x x f x x. Lim. Lim. Lim. Lim Lim. y 13x b b 10 b So the equation of the tangent line is y 13x

1,3. f x x f x x. Lim. Lim. Lim. Lim Lim. y 13x b b 10 b So the equation of the tangent line is y 13x 1.5 Topics: The Derivative lutions 1. Use the limit definition of derivative (the one with x in it) to find f x given f x 4x 5x 6 4 x x 5 x x 6 4x 5x 6 f x x f x f x x0 x x0 x xx x x x x x 4 5 6 4 5 6

More information

Multiple Random Variables

Multiple Random Variables Multiple Random Variables Joint Probability Density Let X and Y be two random variables. Their joint distribution function is F ( XY x, y) P X x Y y. F XY ( ) 1, < x

More information

Lecture 1. 1, 0 x < 1 2 1, 2. x < 1, 0, elsewhere. 1

Lecture 1. 1, 0 x < 1 2 1, 2. x < 1, 0, elsewhere. 1 0 - - -3 Lecture Introductory mathematical ideas The term wavelet, meaning literally little wave, originated in the early 980s in its French version ondelette in the work of Morlet and some French seismologists.

More information

Conditional distributions. Conditional expectation and conditional variance with respect to a variable.

Conditional distributions. Conditional expectation and conditional variance with respect to a variable. Conditional distributions Conditional expectation and conditional variance with respect to a variable Probability Theory and Stochastic Processes, summer semester 07/08 80408 Conditional distributions

More information

Math 210, Final Exam, Fall 2010 Problem 1 Solution. v cosθ = u. v Since the magnitudes of the vectors are positive, the sign of the dot product will

Math 210, Final Exam, Fall 2010 Problem 1 Solution. v cosθ = u. v Since the magnitudes of the vectors are positive, the sign of the dot product will Math, Final Exam, Fall Problem Solution. Let u,, and v,,3. (a) Is the angle between u and v acute, obtuse, or right? (b) Find an equation for the plane through (,,) containing u and v. Solution: (a) The

More information

MATH 307: Problem Set #3 Solutions

MATH 307: Problem Set #3 Solutions : Problem Set #3 Solutions Due on: May 3, 2015 Problem 1 Autonomous Equations Recall that an equilibrium solution of an autonomous equation is called stable if solutions lying on both sides of it tend

More information

Announcements. Topics: Homework: - sections 4.5 and * Read these sections and study solved examples in your textbook!

Announcements. Topics: Homework: - sections 4.5 and * Read these sections and study solved examples in your textbook! Announcements Topics: - sections 4.5 and 5.1-5.5 * Read these sections and study solved examples in your textbook! Homework: - review lecture notes thoroughly - work on practice problems from the textbook

More information

&& && F( u)! "{ f (x)} = f (x)e # j 2$ u x. f (x)! " #1. F(u,v) = f (x, y) e. f (x, y) = 2D Fourier Transform. Fourier Transform - review.

&& && F( u)! { f (x)} = f (x)e # j 2$ u x. f (x)!  #1. F(u,v) = f (x, y) e. f (x, y) = 2D Fourier Transform. Fourier Transform - review. 2D Fourier Transfor 2-D DFT & Properties 2D Fourier Transfor 1 Fourier Transfor - review 1-D: 2-D: F( u)! "{ f (x)} = f (x)e # j 2$ u x % & #% dx f (x)! " #1 { F(u) } = F(u)e j 2$ u x du F(u,v) = f (x,

More information

Lecture 10. (2) Functions of two variables. Partial derivatives. Dan Nichols February 27, 2018

Lecture 10. (2) Functions of two variables. Partial derivatives. Dan Nichols February 27, 2018 Lecture 10 Partial derivatives Dan Nichols nichols@math.umass.edu MATH 233, Spring 2018 University of Massachusetts February 27, 2018 Last time: functions of two variables f(x, y) x and y are the independent

More information

DISCRETE FOURIER TRANSFORM

DISCRETE FOURIER TRANSFORM DD2423 Image Processing and Computer Vision DISCRETE FOURIER TRANSFORM Mårten Björkman Computer Vision and Active Perception School of Computer Science and Communication November 1, 2012 1 Terminology:

More information

Calculus II/III Summer Packet

Calculus II/III Summer Packet Calculus II/III Summer Packet First of all, have a great summer! Enjoy your time away from school. Come back fired up and ready to learn. I know that I will be ready to have a great year of calculus with

More information

XII. Elementary Function

XII. Elementary Function XII. Elementary Function Yuxi Fu BASICS, Shanghai Jiao Tong University What is the class of arithmetic functions we use in mathematics? Computability Theory, by Y. Fu XII. Elementary Function 1 / 17 Definition

More information

EXAMINATION: MATHEMATICAL TECHNIQUES FOR IMAGE ANALYSIS

EXAMINATION: 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 information

The Continuous-time Fourier

The Continuous-time Fourier The Continuous-time Fourier Transform Rui Wang, Assistant professor Dept. of Information and Communication Tongji University it Email: ruiwang@tongji.edu.cn Outline Representation of Aperiodic signals:

More information

Problem Set 8 Mar 5, 2004 Due Mar 10, 2004 ACM 95b/100b 3pm at Firestone 303 E. Sterl Phinney (2 pts) Include grading section number

Problem Set 8 Mar 5, 2004 Due Mar 10, 2004 ACM 95b/100b 3pm at Firestone 303 E. Sterl Phinney (2 pts) Include grading section number Problem Set 8 Mar 5, 24 Due Mar 1, 24 ACM 95b/1b 3pm at Firestone 33 E. Sterl Phinney (2 pts) Include grading section number Useful Readings: For Green s functions, see class notes and refs on PS7 (esp

More information

Stat 366 A1 (Fall 2006) Midterm Solutions (October 23) page 1

Stat 366 A1 (Fall 2006) Midterm Solutions (October 23) page 1 Stat 366 A1 Fall 6) Midterm Solutions October 3) page 1 1. The opening prices per share Y 1 and Y measured in dollars) of two similar stocks are independent random variables, each with a density function

More information

ALGEBRAIC GEOMETRY HOMEWORK 3

ALGEBRAIC GEOMETRY HOMEWORK 3 ALGEBRAIC GEOMETRY HOMEWORK 3 (1) Consider the curve Y 2 = X 2 (X + 1). (a) Sketch the curve. (b) Determine the singular point P on C. (c) For all lines through P, determine the intersection multiplicity

More information

1 Solution to Problem 2.1

1 Solution to Problem 2.1 Solution to Problem 2. I incorrectly worked this exercise instead of 2.2, so I decided to include the solution anyway. a) We have X Y /3, which is a - function. It maps the interval, ) where X lives) onto

More information

Lecture 1.4: Inner products and orthogonality

Lecture 1.4: Inner products and orthogonality Lecture 1.4: Inner products and orthogonality Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4340, Advanced Engineering Mathematics M.

More information

x n cos 2x dx. dx = nx n 1 and v = 1 2 sin(2x). Andreas Fring (City University London) AS1051 Lecture Autumn / 36

x n cos 2x dx. dx = nx n 1 and v = 1 2 sin(2x). Andreas Fring (City University London) AS1051 Lecture Autumn / 36 We saw in Example 5.4. that we sometimes need to apply integration by parts several times in the course of a single calculation. Example 5.4.4: For n let S n = x n cos x dx. Find an expression for S n

More information

Properties of Transformations

Properties of Transformations 6. - 6.4 Properties of Transformations P. Danziger Transformations from R n R m. General Transformations A general transformation maps vectors in R n to vectors in R m. We write T : R n R m to indicate

More information

Lecture 19: Simple Applications of Fourier Analysis & Convolution. Fourier Analysis

Lecture 19: Simple Applications of Fourier Analysis & Convolution. Fourier Analysis Lecture 19: Simple Applications of & Convolution Recall I Let f : {0, 1} n R be a function Let = 2 n Inner product of two functions is defined as follows f, g := 1 f (x)g(x) For S {0, 1} n, define the

More information

FOURIER SERIES, HAAR WAVELETS AND FAST FOURIER TRANSFORM

FOURIER SERIES, HAAR WAVELETS AND FAST FOURIER TRANSFORM FOURIER SERIES, HAAR WAVELETS AD FAST FOURIER TRASFORM VESA KAARIOJA, JESSE RAILO AD SAMULI SILTAE Abstract. This handout is for the course Applications of matrix computations at the University of Helsinki

More information

ISE I Brief Lecture Notes

ISE I Brief Lecture Notes ISE I Brief Lecture Notes 1 Partial Differentiation 1.1 Definitions Let f(x, y) be a function of two variables. The partial derivative f/ x is the function obtained by differentiating f with respect to

More information

Advanced Eng. Mathematics

Advanced Eng. Mathematics Koya University Faculty of Engineering Petroleum Engineering Department Advanced Eng. Mathematics Lecture 6 Prepared by: Haval Hawez E-mail: haval.hawez@koyauniversity.org 1 Second Order Linear Ordinary

More information

Lecture2A--Model QM Problems with Exact Solutions (1-D)

Lecture2A--Model QM Problems with Exact Solutions (1-D) Lecture2A--Model QM Problems with Exact Solutions (1-D) (Ch 2.2-Levine, 3-3 Atkins, Ch. 2-R&S) 1. Free Particle -- If there is no potential then Schroedinger Equation becomes T (x) = E (x) ==> -(h 2 /2m)

More information

Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2004 Lecture 2 Linear Systems. Topics

Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2004 Lecture 2 Linear Systems. Topics Bioegieerig 280A Priciples of Biomedical Imagig Fall Quarter 2004 Lecture 2 Liear Systems Topics 1. Liearity 2. Impulse Respose ad Delta fuctios 3. Superpositio Itegral 4. Shift Ivariace 5. 1D ad 2D covolutio

More information

Let R be the line parameterized by x. Let f be a complex function on R that is integrable. The Fourier transform ˆf = F f is. e ikx f(x) dx. (1.

Let R be the line parameterized by x. Let f be a complex function on R that is integrable. The Fourier transform ˆf = F f is. e ikx f(x) dx. (1. Chapter 1 Fourier transforms 1.1 Introduction Let R be the line parameterized by x. Let f be a complex function on R that is integrable. The Fourier transform ˆf = F f is ˆf(k) = e ikx f(x) dx. (1.1) It

More information

1. Spaces of Functions

1. Spaces of Functions FUNCTION SPACES Analysis of systems of linear algebraic equations leads naturally to the notion of a linear space in which vectors (i.e., the unknowns in the problem) formation of linear combinations (i.e.,

More information

Math 276, Spring 2007 Additional Notes on Vectors

Math 276, Spring 2007 Additional Notes on Vectors Math 276, Spring 2007 Additional Notes on Vectors 1.1. Real Vectors. 1. Scalar Products If x = (x 1,..., x n ) is a vector in R n then the length of x is x = x 2 1 + + x2 n. We sometimes use the notation

More information

Joint Distributions. (a) Scalar multiplication: k = c d. (b) Product of two matrices: c d. (c) The transpose of a matrix:

Joint Distributions. (a) Scalar multiplication: k = c d. (b) Product of two matrices: c d. (c) The transpose of a matrix: Joint Distributions Joint Distributions A bivariate normal distribution generalizes the concept of normal distribution to bivariate random variables It requires a matrix formulation of quadratic forms,

More information

REAL ANALYSIS II HOMEWORK 3. Conway, Page 49

REAL ANALYSIS II HOMEWORK 3. Conway, Page 49 REAL ANALYSIS II HOMEWORK 3 CİHAN BAHRAN Conway, Page 49 3. Let K and k be as in Proposition 4.7 and suppose that k(x, y) k(y, x). Show that K is self-adjoint and if {µ n } are the eigenvalues of K, each

More information

Math 265H: Calculus III Practice Midterm II: Fall 2014

Math 265H: Calculus III Practice Midterm II: Fall 2014 Name: Section #: Math 65H: alculus III Practice Midterm II: Fall 14 Instructions: This exam has 7 problems. The number of points awarded for each question is indicated in the problem. Answer each question

More information

Space-Frequency Atoms

Space-Frequency Atoms Space-Frequency Atoms FREQUENCY FREQUENCY SPACE SPACE FREQUENCY FREQUENCY SPACE SPACE Figure 1: Space-frequency atoms. Windowed Fourier Transform 1 line 1 0.8 0.6 0.4 0.2 0-0.2-0.4-0.6-0.8-1 0 100 200

More information

Fourier Reconstruction. Polar Version of Inverse FT. Filtered Backprojection [ ] Bioengineering 280A Principles of Biomedical Imaging

Fourier Reconstruction. Polar Version of Inverse FT. Filtered Backprojection [ ] Bioengineering 280A Principles of Biomedical Imaging Fourier Reconstruction Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2014 CT/Fourier Lecture 5 F Interpolate onto Cartesian grid with appropriate density correction then take inverse

More information

Lecture 9 Metric spaces. The contraction fixed point theorem. The implicit function theorem. The existence of solutions to differenti. equations.

Lecture 9 Metric spaces. The contraction fixed point theorem. The implicit function theorem. The existence of solutions to differenti. equations. Lecture 9 Metric spaces. The contraction fixed point theorem. The implicit function theorem. The existence of solutions to differential equations. 1 Metric spaces 2 Completeness and completion. 3 The contraction

More information

Inner products. Theorem (basic properties): Given vectors u, v, w in an inner product space V, and a scalar k, the following properties hold:

Inner products. Theorem (basic properties): Given vectors u, v, w in an inner product space V, and a scalar k, the following properties hold: Inner products Definition: An inner product on a real vector space V is an operation (function) that assigns to each pair of vectors ( u, v) in V a scalar u, v satisfying the following axioms: 1. u, v

More information

A Tilt at TILFs. Rod Nillsen University of Wollongong. This talk is dedicated to Gary H. Meisters

A Tilt at TILFs. Rod Nillsen University of Wollongong. This talk is dedicated to Gary H. Meisters A Tilt at TILFs Rod Nillsen University of Wollongong This talk is dedicated to Gary H. Meisters Abstract In this talk I will endeavour to give an overview of some aspects of the theory of Translation Invariant

More information

Math 131 Exam 2 Spring 2016

Math 131 Exam 2 Spring 2016 Math 3 Exam Spring 06 Name: ID: 7 multiple choice questions worth 4.7 points each. hand graded questions worth 0 points each. 0. free points (so the total will be 00). Exam covers sections.7 through 3.0

More information

Sophus Lie s Approach to Differential Equations

Sophus Lie s Approach to Differential Equations Sophus Lie s Approach to Differential Equations IAP lecture 2006 (S. Helgason) 1 Groups Let X be a set. A one-to-one mapping ϕ of X onto X is called a bijection. Let B(X) denote the set of all bijections

More information

Continuous r.v practice problems

Continuous r.v practice problems Continuous r.v practice problems SDS 321 Intro to Probability and Statistics 1. (2+2+1+1 6 pts) The annual rainfall (in inches) in a certain region is normally distributed with mean 4 and standard deviation

More information

Engg. Math. I. Unit-I. Differential Calculus

Engg. Math. I. Unit-I. Differential Calculus Dr. Satish Shukla 1 of 50 Engg. Math. I Unit-I Differential Calculus Syllabus: Limits of functions, continuous functions, uniform continuity, monotone and inverse functions. Differentiable functions, Rolle

More information

Solutions of Math 53 Midterm Exam I

Solutions of Math 53 Midterm Exam I Solutions of Math 53 Midterm Exam I Problem 1: (1) [8 points] Draw a direction field for the given differential equation y 0 = t + y. (2) [8 points] Based on the direction field, determine the behavior

More information

Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem

Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem L u x f x BC u x g x with the weak problem find u V such that B u,v

More information