Other Optimization Techniques. Similar to steepest descent, but slightly different way of choosing direction of next step: x r +1 = x.
|
|
- Lynne Howard
- 5 years ago
- Views:
Transcription
1 Conjugate Gradient Other Optimization Techniques Similar to steepest descent, but slightly different way of choosing direction of next step: x r +1 = x r + r s r s 0 = g 0 s r +1 = g r +1 + r +1 s r { new term r is chosen to minimize h( x r +1 ). This yields g t r +1 g r = 0 Here we allow a further step in the s r direction. One choice (Fletcher - Reeves) for r +1 is r +1 = g r +1 g r Winter Semester 006/7 Computational Physics I Lecture 1 1
2 Newton-Raphson Assume the function that we want to minimize is twice differentiable. Then, a Taylor expansion gives h( x + ) a + b t + 1 t C where Now a = h( x ), h( x + b = ) b + C h( x ) = g ( x ), C = h x i x j = H Because C is symmetric (check) For an extremum, we have b + C = 0 = C 1 b or x r +1 = x r H( x r ) 1 g ( x r ) Winter Semester 006/7 Computational Physics I Lecture 1
3 Newton-Raphson x r +1 = x r H( x r ) 1 g ( x r ) i.e., the search direction is s = H( x r ) 1 g ( x r ) and =1 This converges quickly (if you start with a good guess), but the penalty is that the Hessian needs to be calculated (usually numerically) Again, convergence is when s is sufficiently small How would we calculate the Hessian numerically? Use Lagrange polynomial in several dimensions and work it out Winter Semester 006/7 Computational Physics I Lecture 1 3
4 Bounded Regions The standard tool for minimization in particle physics is the MINUIT program (CERN library). It has also made its way well outside the particle physics community. Author: Fred James Here is how MINUIT handles bounded search regions - it transforms the parameter to be optimized as follows: a = arcsin b a 1 = a + b a is the exernal (user) parameter is the internal parameter MINUIT is available within PAW, ROOT ( sin +1) Winter Semester 006/7 Computational Physics I Lecture 1 4
5 MINUIT MINUIT uses a (variable metric) conjugate gradient search algorithm (along with others). Basic idea: assume that the function to minimize can be approximated by a quadratic form near the minimum build up iteratively an approximation for the inverse of the Hessian matrix. Recall h( x + ) h( x ) + h( x ) + 1 t H the approximation for the Hessian is updated as follows: H i +1 = H i + ( x i +1 x i ) ( x i +1 x i ) ( x i +1 x i ) ( h i +1 h i ) H [ i ( h i +1 h i )] H i ( h i +1 h i ) ( h i +1 h i ) H i ( h i +1 h i ) [ ] where the symbol represents an outer product of two vectors (a matrix) a b ( ) ij = a i b j Winter Semester 006/7 Computational Physics I Lecture 1 5
6 Fourier Transforms Fourier transforms are very important as a way of summarizing the data with a few parameters because the transform of the data is itself very interesting (e.g., power spectrum, momentumcoordinate space representation, ) H( f ) = h(t) e ift dt H( f ) frequency domain representation h(t) = H( f ) e ift dt h(t) time domain representation Warning: there is no unanimity on factors in front of the integral. Often the angular frequency is used = f H() = h(t) e it dt h(t) = 1 H() e it d Winter Semester 006/7 Computational Physics I Lecture 1 6
7 Fourier Transform Fourier Transform is a linear operation: transform of the sum of two functions is the sum of the transforms the transform of a constant times a function is constant times the transform h(t) real H( f ) = [ H( f )] * h(t) imaginary H( f ) = [ H( f )] * h(t) even H( f ) = H( f ) h(t) odd H( f ) = H( f ) h(t) real,even H( f ) real, even h(t) real,odd H( f ) imaginary, odd h(t) imaginary,even H( f ) imaginary, even h(t) imaginary,odd H( f ) real, odd Winter Semester 006/7 Computational Physics I Lecture 1 7
8 Fourier Transform Further properties: h(at) 1 a H 1 b h t b f a Hbf ( ) h(t t 0 ) H( f )e ift 0 h(t)e if 0 t H( f f 0 ) We are typically interested in the Fourier analysis of a discretely sampled data set. Define the time step (taken to be constant here) as. The sampling rate (frequency) is 1/. Define the samples as h n = h(n) n =,3,,1,0,1,,3, Winter Semester 006/7 Computational Physics I Lecture 1 8
9 Nyquist frequency f c 1 Nyquist frequency This is the highest frequency which can be resolved with a sampling frequency f=1/. If a continuous function h(t) is limited in frequency components to frequencies less than f c, then h(t) is completely determined by its samples h n. It can then be written as follows: h(t) = n = sin[ f h c (t n) ] n (t n) However, if there are frequency components which are higher than f c, then they will be spuriously moved in the range f<f c (aliasing). Winter Semester 006/7 Computational Physics I Lecture 1 9
10 Data Example Fit Aliasing 10 Hz sampling 5 Hz sampling Conditions are: Sine wave with f=4 Hz, phase offset =0.1 i.e., h(t) = sin( + ft) = sin( t) Winter Semester 006/7 Computational Physics I Lecture 1 10
11 H( f ) = sin( t)e ift dt Example = sin(0.1)cos(8t)e ift dt + cos(0.1)sin(8t)e ift dt e 8it + e 8it = sin(0.1) e ift dt + cos(0.1) Recall the relation: eifx df so we have H( f ) = sin(0.1) = sin(0.1) = sin(0.1) e 8it e 8it Winter Semester 006/7 Computational Physics I Lecture 1 11 e ift dt = (x) where (x) is the Dirac Delta function e 8it + e 8it e ift dt eit( f + 4) + e it( f 4) dt + cos(0.1) + cos(0.1) e 8it e 8it e ift dt eit( f + 4) e it( f 4) dt [ ( f + 4) + ( f 4) ] + cos(0.1) [ ( f + 4) ( f 4) ]
12 Discrete Fourier Transform Suppose we have N consecutive sampled points h k h(t k ), t k k, k = 0,1,,,N 1 We can extract the amplitude for N frequency components since we have N data points. Define the frequency components as f n n 1 n = N N,..., N (take N even) Note: there are N+1 frequencies, but we will find that the two at the ends are equal, so only N independent. Negative frequencies allows us to include sine and cosine terms. So N 1 H( f n ) = h(t)e ifnt dt h k e if nt k = h k e if nk = h k e ikn / N H n N 1 k = 0 h k e ikn / N k = 0 Winter Semester 006/7 Computational Physics I Lecture 1 1 N 1 k = 0 Discrete Fourier Transform N 1 k = 0
13 Discrete Fourier Transform The discrete fourier transform does not depend on any dimensional parameters. Note H n = H N n [ e ik(n n)/n = e ik e ikn / N = e ikn / N ] In particular H N / = H N / We can therefore rewrite the sum as follows N 1 H n h k e ikn / N n = 0,,N 1 k = 0 Discrete inverse Fourier transform h k = 1 N N 1 n = 0 H n e ikn / N k = 0,,N 1 Winter Semester 006/7 Computational Physics I Lecture 1 13
14 Discrete Fourier Transform * * Get the discrete Fourier components * Do n=0,63 Hn(n,1)=0.D0 Hn(n,)=0.D0 Do k=0,63 Hn(n,1)=Hn(n,1) & +amplitude(k,1)*dcos(twopi*k*n/64.) & -amplitude(k,)*dsin(twopi*k*n/64.) Hn(n,)=Hn(n,) & +amplitude(k,1)*dsin(twopi*k*n/64.) & +amplitude(k,)*dcos(twopi*k*n/64.) Enddo Write (11,*) N,Hn(N,1),Hn(N,) Enddo amplitude(k,1) real components amplitude(k,) imaginary components Winter Semester 006/7 Computational Physics I Lecture 1 14
15 Discrete Fourier Transform Let s try it out on our sine wave data: Recall, signal f=4 Hz Here s the result (64 points): Large components are: f 5 = = 3.9, f 6 = = 4.06 f 38 = f 64 6 = f 6 f 39 = f 64 5 = f 5 cos and sin needed because of phase offset. Winter Semester 006/7 Computational Physics I Lecture 1 15
16 Discrete Fourier Transform Here s the inverse transform * * Now we try the inverse transform * Do n=0,63 amplitude(n,1)=0.d0 amplitude(n,)=0.d0 Do k=0,63 amplitude(n,1)=amplitude(n,1) & +Hn(k,1)*dcos(twopi*k*n/64.) & +Hn(k,)*dsin(twopi*k*n/64.) amplitude(n,)=amplitude(n,) & -Hn(k,1)*dsin(twopi*k*n/64.) & +Hn(k,)*dcos(twopi*k*n/64.) Enddo Write (1,*) N,amplitude(n,1)/64.,amplitude(n,)/64. Enddo * Winter Semester 006/7 Computational Physics I Lecture 1 16
17 Fast Fourier Transform The discrete Fourier transforms as we described it requires a sum over N terms for each of the N components. I.e., the number of operations scales as N. A large part of the success of Fourier transforms for analysis of electronic signals, optical images, x-ray tomography,, results from the fact that a numerical algorithm was found which requires of order Nlog N operations - the socalled Fast Fourier Transform (FFT). Here is how it works: F k = N 1 e ijk / N f j = e ik( j )/N f j + e ik( j +1)/ N f j +1 N /1 N /1 N /1 Winter Semester 006/7 Computational Physics I Lecture 1 17 = e ikj /(N /) f j + W k e ikj /(N /) f j +1 = F k e + W k F k o N /1 where W e i / N What is won? The sums in the individual terms in the last line only have 1/ as many terms, and the same factor appears
18 Fast Fourier Transform To see how this works in detail, we take an explicit example of having 8 data points (taking a power of is important! If you don t have enough data, pad with zeroes). F k = F k e + W k F k o where W e i / N, F k e 3 W kj f j, F o k W kj f j +1 Now use a binary representation for the index k=4k +k 1 +k 0 where the k i s are 0,1. Then, W kj = ( i /8 e ) (4k + k + k ) j 1 0 = e i(k + k 1 /+ k /4) 0 = e i(k /+ k /4) 1 0 = W j(k + k 1 0 ) i.e., the k bit is irrelevant. So, e F k = F (k1,k 0 ) o + W k F (k1,k 0 ) 3 e F (k1,k 0 ) 3 W j(k + k 1 ) 0 f j F (k1,k 0 ) o 3 W j(k + k 1 ) 0 f j +1 Winter Semester 006/7 Computational Physics I Lecture 1 18
19 Let s try again: F k e = 1 3 Fast Fourier Transform W j(k + k 1 ) 0 f j = W ( j )(k + k 1 ) 0 f ( j ) + W ( j +1)(k + k 1 ) 0 f ( j +1) 1 = W ( j )(k + k 1 ) 0 f ( j ) + W (k + k 1 0 ) F k e = F k ee + W (k 1 + k 0 ) F k eo 1 1 W ( j )(k 1 + k 0 ) f ( j +1) and W 4 j(k + k ) 1 0 = ( i /8 e ) (8k + 4 k 1 0 ) j = W jk 0 F k o = F k oe + W (k 1 + k 0 ) F k oo Can perform one more step: F ee k = F eee + W 4k 0 F eeo F eee = f 0 F eeo = f 4 The sums have disappeared! Winter Semester 006/7 Computational Physics I Lecture 1 19
20 Fast Fourier Transform The final pieces are: F ee k = F eee + W 4k 0 F eeo = f 0 + W 4k 0 f 4 F eo k = F eoe + W 4k 0 F eoo = f + W 4k 0 f 6 F oe k = F oee + W 4k 0 F oeo = f 1 + W 4k 0 f 5 F oo k = F ooe + W 4k 0 F ooo = f 3 + W 4k 0 f 7 Then, F e k = F ee + W (k + k ) 1 0 F eo F o k = F oe + W (k + k ) 1 0 F oo Note k 0 =0,1 So, need 8 multiplications and 8 additions for this step Here k 0 =0,1 k 1 =0,1 So, again 8 multiplications and 8 additions for this step Finally F k = F e + W k F o again 8 multiplications and 8 additions for this step Winter Semester 006/7 Computational Physics I Lecture 1 0
21 Fast Fourier Transform So we need N operations per level, and there are log N levels. The scaling of the computational time is therefore Nlog N rather than N. E.g., N=1000 Nlog N1000*10=10 4 N =10 6 How to implement in practice. Note that the trick is to find out which value of n corresponds to which pattern of e,o in F eoeooeoe... = f n Answer: reverse pattern of e,o. Assign e=0, o=1, and the binary value gives n. eee eee examples oeo eoe 010 Winter Semester 006/7 Computational Physics I Lecture 1 1
22 Some examples. Agric. Food Chem., 5 (0), , /jf04940e S (04) Web Release Date: September 9, 004 Copyright 004 American Chemical Society Discrimination of Olives According to Fruit Quality Using Fourier Transform Raman Spectroscopy and Pattern Recognition Techniques Barbara Muik, Bernhard Lendl, Antonio Molina-Díaz, Domingo Ortega-Calderón,# and María José Ayora-Cañada* Department of Physical and Analytical Chemistry, University of Jaén, Paraje las Lagunillas s/n, E-3071 Jaén, Spain; Institute of Chemical Technologies and Analytics, Vienna University of Technology, Getreidemarkt 9/164, A-1060 Wien, Austria; and CIFA Venta del Llano, IFAPA, Ctra. Bailén-Motril km 18.5, E-360 Mengíbar, Jaén, Spain Multiplication of large integers The fastest known algorithms for the multiplication of large integers or polynomials are based on the discrete Fourier transform: the sequences of digits or coefficients are interpreted as vectors whose convolution needs to be computed; in order to do this, they are first Fourier-transformed, then multiplied component-wise, then transformed back. Winter Semester 006/7 Computational Physics I Lecture 1
23 Power Spectrum The autocorrelation of a function is Corr[y]() = y(t) * y(t + ) dt and the power spectrum is defined as the Fourier transform of the autocorrelation PS[y]( f ) = Corr[y]()e if d For a periodic function, the correlation is often defined as the expectation value. There is no convention on the normalization, so be careful about the values. Best to see which frequencies dominate a given spectrum. Here is a practical approach for a discretely sampled function: C k = N 1 c j e ijk / N k = 0,1,,N 1 Winter Semester 006/7 Computational Physics I Lecture 1 3
24 P(0) = P( f 0 ) = 1 N C 0 Power Spectrum P( f k ) = 1 [ N C k + C ] N N k k = 1,,, 1 P( f c ) = P( f N / ) = 1 N C N / where only positive frequencies are considered: f k k N = f c k N k = 0,1,, N Let s try it out on our example: The 4Hz frequency is picked out. Winter Semester 006/7 Computational Physics I Lecture 1 4
25 Noise We now add some noise to our spectrum (Gaussian smearing with =0.5) and see what happens: Winter Semester 006/7 Computational Physics I Lecture 1 5
26 Exercizes 1. Solve the minimization problem from last lecture with MINUIT.. Generate 64 data points using f (t) = cos( /4 + f 1 t) + cos(f t) f 1 = 0.5, f = 1 = 0. and fit with a discrete Fourier transform. Extract the power spectrum. Winter Semester 006/7 Computational Physics I Lecture 1 6
Lecture 7. Fourier Analysis
Lecture 7 Fourier Analysis Summary Lecture 6 Minima and maxima 1 dimension : Bracket minima : 3 values of f(x) : f(2) < f(1) and f(2)
More informationProject 2 The Threaded Two Dimensional Discrete Fourier Transform
ECE6122/6122 Spring Semester, 2017 Project 2 The Threaded Two Dimensional Discrete Fourier Transform Assigned: Feb 1, 2017 Due: Feb 15, 2017, 11:59pm Given a one dimensional array of complex or real input
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 informationE The Fast Fourier Transform
Fourier Transform Methods in Finance By Umberto Cherubini Giovanni Della Lunga Sabrina Mulinacci Pietro Rossi Copyright 2010 John Wiley & Sons Ltd E The Fast Fourier Transform E.1 DISCRETE FOURIER TRASFORM
More informationESS Finite Impulse Response Filters and the Z-transform
9. Finite Impulse Response Filters and the Z-transform We are going to have two lectures on filters you can find much more material in Bob Crosson s notes. In the first lecture we will focus on some of
More informationInformation and Communications Security: Encryption and Information Hiding
Short Course on Information and Communications Security: Encryption and Information Hiding Tuesday, 10 March Friday, 13 March, 2015 Lecture 5: Signal Analysis Contents The complex exponential The complex
More informationMathematics for Chemists 2 Lecture 14: Fourier analysis. Fourier series, Fourier transform, DFT/FFT
Mathematics for Chemists 2 Lecture 14: Fourier analysis Fourier series, Fourier transform, DFT/FFT Johannes Kepler University Summer semester 2012 Lecturer: David Sevilla Fourier analysis 1/25 Remembering
More informationRoot Finding and Optimization
Root Finding and Optimization Root finding, solving equations, and optimization are very closely related subjects, which occur often in practical applications. Root finding : f (x) = 0 Solve for x Equation
More informationSignal Processing COS 323
Signal Processing COS 323 Digital Signals D: functions of space or time e.g., sound 2D: often functions of 2 spatial dimensions e.g. images 3D: functions of 3 spatial dimensions CAT, MRI scans or 2 space,
More informationTime and Spatial Series and Transforms
Time and Spatial Series and Transforms Z- and Fourier transforms Gibbs' phenomenon Transforms and linear algebra Wavelet transforms Reading: Sheriff and Geldart, Chapter 15 Z-Transform Consider a digitized
More informationFourier Transforms For additional information, see the classic book The Fourier Transform and its Applications by Ronald N. Bracewell (which is on the shelves of most radio astronomers) and the Wikipedia
More informationNotes 07 largely plagiarized by %khc
Notes 07 largely plagiarized by %khc Warning This set of notes covers the Fourier transform. However, i probably won t talk about everything here in section; instead i will highlight important properties
More informationTutorial Sheet #2 discrete vs. continuous functions, periodicity, sampling
2.39 utorial Sheet #2 discrete vs. continuous functions, periodicity, sampling We will encounter two classes of signals in this class, continuous-signals and discrete-signals. he distinct mathematical
More information12.1 Fourier Transform of Discretely Sampled Data
494 Chapter 12. Fast Fourier Transform now is: The PSD-per-unit-time converges to finite values at all frequencies except those where h(t) has a discrete sine-wave (or cosine-wave) component of finite
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 informationANALOG AND DIGITAL SIGNAL PROCESSING ADSP - Chapter 8
ANALOG AND DIGITAL SIGNAL PROCESSING ADSP - Chapter 8 Fm n N fnt ( ) e j2mn N X() X() 2 X() X() 3 W Chap. 8 Discrete Fourier Transform (DFT), FFT Prof. J.-P. Sandoz, 2-2 W W 3 W W x () x () x () 2 x ()
More informationFourier analysis of discrete-time signals. (Lathi Chapt. 10 and these slides)
Fourier analysis of discrete-time signals (Lathi Chapt. 10 and these slides) Towards the discrete-time Fourier transform How we will get there? Periodic discrete-time signal representation by Discrete-time
More informationWe consider the problem of finding a polynomial that interpolates a given set of values:
Chapter 5 Interpolation 5. Polynomial Interpolation We consider the problem of finding a polynomial that interpolates a given set of values: x x 0 x... x n y y 0 y... y n where the x i are all distinct.
More informationNumerical Optimization
Numerical Optimization Unit 2: Multivariable optimization problems Che-Rung Lee Scribe: February 28, 2011 (UNIT 2) Numerical Optimization February 28, 2011 1 / 17 Partial derivative of a two variable function
More informationIn many diverse fields physical data is collected or analysed as Fourier components.
1. Fourier Methods In many diverse ields physical data is collected or analysed as Fourier components. In this section we briely discuss the mathematics o Fourier series and Fourier transorms. 1. Fourier
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 information6. Advanced Numerical Methods. Monte Carlo Methods
6. Advanced Numerical Methods Part 1: Part : Monte Carlo Methods Fourier Methods Part 1: Monte Carlo Methods 1. Uniform random numbers Generating uniform random numbers, drawn from the pdf U[0,1], is fairly
More informationPART 1. Review of DSP. f (t)e iωt dt. F(ω) = f (t) = 1 2π. F(ω)e iωt dω. f (t) F (ω) The Fourier Transform. Fourier Transform.
PART 1 Review of DSP Mauricio Sacchi University of Alberta, Edmonton, AB, Canada The Fourier Transform F() = f (t) = 1 2π f (t)e it dt F()e it d Fourier Transform Inverse Transform f (t) F () Part 1 Review
More informationCorrelator I. Basics. Chapter Introduction. 8.2 Digitization Sampling. D. Anish Roshi
Chapter 8 Correlator I. Basics D. Anish Roshi 8.1 Introduction A radio interferometer measures the mutual coherence function of the electric field due to a given source brightness distribution in the sky.
More informationSummary of lecture 1. E x = E x =T. X T (e i!t ) which motivates us to define the energy spectrum Φ xx (!) = jx (i!)j 2 Z 1 Z =T. 2 d!
Summary of lecture I Continuous time: FS X FS [n] for periodic signals, FT X (i!) for non-periodic signals. II Discrete time: DTFT X T (e i!t ) III Poisson s summation formula: describes the relationship
More informationOptimization. Escuela de Ingeniería Informática de Oviedo. (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30
Optimization Escuela de Ingeniería Informática de Oviedo (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30 Unconstrained optimization Outline 1 Unconstrained optimization 2 Constrained
More informationNumerical solution of Least Squares Problems 1/32
Numerical solution of Least Squares Problems 1/32 Linear Least Squares Problems Suppose that we have a matrix A of the size m n and the vector b of the size m 1. The linear least square problem is to find
More informationFall 2011, EE123 Digital Signal Processing
Lecture 6 Miki Lustig, UCB September 11, 2012 Miki Lustig, UCB DFT and Sampling the DTFT X (e jω ) = e j4ω sin2 (5ω/2) sin 2 (ω/2) 5 x[n] 25 X(e jω ) 4 20 3 15 2 1 0 10 5 1 0 5 10 15 n 0 0 2 4 6 ω 5 reconstructed
More informationPower Series Solutions We use power series to solve second order differential equations
Objectives Power Series Solutions We use power series to solve second order differential equations We use power series expansions to find solutions to second order, linear, variable coefficient equations
More informationDiscrete 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 informationModeling Data with Linear Combinations of Basis Functions. Read Chapter 3 in the text by Bishop
Modeling Data with Linear Combinations of Basis Functions Read Chapter 3 in the text by Bishop A Type of Supervised Learning Problem We want to model data (x 1, t 1 ),..., (x N, t N ), where x i is a vector
More informationOptimization II: Unconstrained Multivariable
Optimization II: Unconstrained Multivariable CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Doug James (and Justin Solomon) CS 205A: Mathematical Methods Optimization II: Unconstrained
More informationSECTION FOR DIGITAL SIGNAL PROCESSING DEPARTMENT OF MATHEMATICAL MODELLING TECHNICAL UNIVERSITY OF DENMARK Course 04362 Digital Signal Processing: Solutions to Problems in Proakis and Manolakis, Digital
More informationImages have structure at various scales
Images have structure at various scales Frequency Frequency of a signal is how fast it changes Reflects scale of structure A combination of frequencies 0.1 X + 0.3 X + 0.5 X = Fourier transform Can we
More informationECS455: Chapter 5 OFDM. ECS455: Chapter 5 OFDM. OFDM: Overview. OFDM Applications. Dr.Prapun Suksompong prapun.com/ecs455
ECS455: Chapter 5 OFDM OFDM: Overview Let S = (S 1, S 2,, S ) contains the information symbols. S IFFT FFT Inverse fast Fourier transform Fast Fourier transform 1 Dr.Prapun Suksompong prapun.com/ecs455
More informationTutorial 4. Fast Fourier Transforms Phase factors runs around on the unit circle contains Wave forms
Tutorial 4. Fast Fourier Transforms Phase factors There are functions that produce roots-of-one as a function of time (t) or place (x). A good example is a Bloch function φ(x) = exp(i kx) or the phase
More informationNeural Network Training
Neural Network Training Sargur Srihari Topics in Network Training 0. Neural network parameters Probabilistic problem formulation Specifying the activation and error functions for Regression Binary classification
More information1 Understanding Sampling
1 Understanding Sampling Summary. In Part I, we consider the analysis of discrete-time signals. In Chapter 1, we consider how discretizing a signal affects the signal s Fourier transform. We derive the
More informationSecondary Math 2H Unit 3 Notes: Factoring and Solving Quadratics
Secondary Math H Unit 3 Notes: Factoring and Solving Quadratics 3.1 Factoring out the Greatest Common Factor (GCF) Factoring: The reverse of multiplying. It means figuring out what you would multiply together
More informationCorrelation, discrete Fourier transforms and the power spectral density
Correlation, discrete Fourier transforms and the power spectral density visuals to accompany lectures, notes and m-files by Tak Igusa tigusa@jhu.edu Department of Civil Engineering Johns Hopkins University
More informationQuantum Mechanics-I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras. Lecture - 21 Square-Integrable Functions
Quantum Mechanics-I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras Lecture - 21 Square-Integrable Functions (Refer Slide Time: 00:06) (Refer Slide Time: 00:14) We
More informationNonlinear Optimization
Nonlinear Optimization (Com S 477/577 Notes) Yan-Bin Jia Nov 7, 2017 1 Introduction Given a single function f that depends on one or more independent variable, we want to find the values of those variables
More informationConstrained optimization. Unconstrained optimization. One-dimensional. Multi-dimensional. Newton with equality constraints. Active-set method.
Optimization Unconstrained optimization One-dimensional Multi-dimensional Newton s method Basic Newton Gauss- Newton Quasi- Newton Descent methods Gradient descent Conjugate gradient Constrained optimization
More informationSection 7.4: Inverse Laplace Transform
Section 74: Inverse Laplace Transform A natural question to ask about any function is whether it has an inverse function We now ask this question about the Laplace transform: given a function F (s), will
More informationFunctions of Several Variables
Functions of Several Variables The Unconstrained Minimization Problem where In n dimensions the unconstrained problem is stated as f() x variables. minimize f()x x, is a scalar objective function of vector
More informationHow to Multiply. 5.5 Integer Multiplication. Complex Multiplication. Integer Arithmetic. Complex multiplication. (a + bi) (c + di) = x + yi.
How to ultiply Slides by Kevin Wayne. Copyright 5 Pearson-Addison Wesley. All rights reserved. integers, matrices, and polynomials Complex ultiplication Complex multiplication. a + bi) c + di) = x + yi.
More informationThe Fourier Transform
The Fourier Transorm Fourier Series Fourier Transorm The Basic Theorems and Applications Sampling Bracewell, R. The Fourier Transorm and Its Applications, 3rd ed. New York: McGraw-Hill, 2. Eric W. Weisstein.
More informationGeneration of bandlimited sync transitions for sine waveforms
Generation of bandlimited sync transitions for sine waveforms Vadim Zavalishin May 4, 9 Abstract A common way to generate bandlimited sync transitions for the classical analog waveforms is the BLEP method.
More informationFourier Transform for Continuous Functions
Fourier Transform for Continuous Functions Central goal: representing a signal by a set of orthogonal bases that are corresponding to frequencies or spectrum. Fourier series allows to find the spectrum
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 informationHermite Interpolation
Jim Lambers MAT 77 Fall Semester 010-11 Lecture Notes These notes correspond to Sections 4 and 5 in the text Hermite Interpolation Suppose that the interpolation points are perturbed so that two neighboring
More information!Sketch f(t) over one period. Show that the Fourier Series for f(t) is as given below. What is θ 1?
Second Year Engineering Mathematics Laboratory Michaelmas Term 998 -M L G Oldfield 30 September, 999 Exercise : Fourier Series & Transforms Revision 4 Answer all parts of Section A and B which are marked
More informationMath 56 Homework 5 Michael Downs
1. (a) Since f(x) = cos(6x) = ei6x 2 + e i6x 2, due to the orthogonality of each e inx, n Z, the only nonzero (complex) fourier coefficients are ˆf 6 and ˆf 6 and they re both 1 2 (which is also seen from
More information17 Solution of Nonlinear Systems
17 Solution of Nonlinear Systems We now discuss the solution of systems of nonlinear equations. An important ingredient will be the multivariate Taylor theorem. Theorem 17.1 Let D = {x 1, x 2,..., x m
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 2 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization 1 Outline 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization Summary 2 Outline 3.2
More information3 rd class Mech. Eng. Dept. hamdiahmed.weebly.com Fourier Series
Definition 1 Fourier Series A function f is said to be piecewise continuous on [a, b] if there exists finitely many points a = x 1 < x 2
More informationA6523 Modeling, Inference, and Mining Jim Cordes, Cornell University
A6523 Modeling, Inference, and Mining Jim Cordes, Cornell University Lecture 19 Modeling Topics plan: Modeling (linear/non- linear least squares) Bayesian inference Bayesian approaches to spectral esbmabon;
More informationFig. 1: Fourier series Examples
FOURIER SERIES AND ITS SOME APPLICATIONS P. Sathyabama Assistant Professor, Department of Mathematics, Bharath collage of Science and Management, Thanjavur, Tamilnadu Abstract: The Fourier series, the
More informationA SuperDARN Tutorial. Kile Baker
A SuperDARN Tutorial Kile Baker Introduction n Multi-pulse Technique and ACFs Why ACFs? Bad and missing lags n Basic Theory of FITACF n Basic Fitting Technique Getting the Power and Width Getting the Velocity
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 informationA6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2013
A6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2013 Lecture 26 Localization/Matched Filtering (continued) Prewhitening Lectures next week: Reading Bases, principal
More informationGATE EE Topic wise Questions SIGNALS & SYSTEMS
www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)
More informationLecture 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 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 information18.03 LECTURE NOTES, SPRING 2014
18.03 LECTURE NOTES, SPRING 2014 BJORN POONEN 7. Complex numbers Complex numbers are expressions of the form x + yi, where x and y are real numbers, and i is a new symbol. Multiplication of complex numbers
More informationPart 3.3 Differentiation Taylor Polynomials
Part 3.3 Differentiation 3..3.1 Taylor Polynomials Definition 3.3.1 Taylor 1715 and Maclaurin 1742) If a is a fixed number, and f is a function whose first n derivatives exist at a then the Taylor polynomial
More informationImage Enhancement in the frequency domain. GZ Chapter 4
Image Enhancement in the frequency domain GZ Chapter 4 Contents In this lecture we will look at image enhancement in the frequency domain The Fourier series & the Fourier transform Image Processing in
More informationNewton-Raphson. Relies on the Taylor expansion f(x + δ) = f(x) + δ f (x) + ½ δ 2 f (x) +..
2008 Lecture 7 starts here Newton-Raphson When the derivative of f(x) is known, and when f(x) is well behaved, the celebrated (and ancient) Newton- Raphson method gives the fastest convergence of all (
More informationMachine Learning Lecture 5
Machine Learning Lecture 5 Linear Discriminant Functions 26.10.2017 Bastian Leibe RWTH Aachen http://www.vision.rwth-aachen.de leibe@vision.rwth-aachen.de Course Outline Fundamentals Bayes Decision Theory
More informationEE123 Digital Signal Processing
Announcements EE Digital Signal Processing otes posted HW due Friday SDR give away Today! Read Ch 9 $$$ give me your names Lecture based on slides by JM Kahn M Lustig, EECS UC Berkeley M Lustig, EECS UC
More informationFourier transform. Stefano Ferrari. Università degli Studi di Milano Methods for Image Processing. academic year
Fourier transform Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Methods for Image Processing academic year 27 28 Function transforms Sometimes, operating on a class of functions
More informationFOURIER ANALYSIS using Python
(version September 5) FOURIER ANALYSIS using Python his practical introduces the following: Fourier analysis of both periodic and non-periodic signals (Fourier series, Fourier transform, discrete Fourier
More informationLecture V : Oscillatory motion and spectral analysis
Lecture V : Oscillatory motion and spectral analysis I. IDEAL PENDULUM AND STABILITY ANALYSIS Let us remind ourselves of the equation of motion for the pendulum. Remembering that the external torque applied
More informationDeterminants of 2 2 Matrices
Determinants In section 4, we discussed inverses of matrices, and in particular asked an important question: How can we tell whether or not a particular square matrix A has an inverse? We will be able
More informationLecture 10: October 27, 2016
Mathematical Toolkit Autumn 206 Lecturer: Madhur Tulsiani Lecture 0: October 27, 206 The conjugate gradient method In the last lecture we saw the steepest descent or gradient descent method for finding
More informationLecture 5. Complex Numbers and Euler s Formula
Lecture 5. Complex Numbers and Euler s Formula University of British Columbia, Vancouver Yue-Xian Li March 017 1 Main purpose: To introduce some basic knowledge of complex numbers to students so that they
More informationHandout 1 Systems of linear equations Gaussian elimination
06-21254 Mathematical Techniques for Computer Science The University of Birmingham Autumn Semester 2017 School of Computer Science c Achim Jung September 27, 2017 Handout 1 Systems of linear equations
More informationSignals, Instruments, and Systems W5. Introduction to Signal Processing Sampling, Reconstruction, and Filters
Signals, Instruments, and Systems W5 Introduction to Signal Processing Sampling, Reconstruction, and Filters Acknowledgments Recapitulation of Key Concepts from the Last Lecture Dirac delta function (
More informationChapter 2: The Fourier Transform
EEE, EEE Part A : Digital Signal Processing Chapter Chapter : he Fourier ransform he Fourier ransform. Introduction he sampled Fourier transform of a periodic, discrete-time signal is nown as the discrete
More informationLecture 7. Logistic Regression. Luigi Freda. ALCOR Lab DIAG University of Rome La Sapienza. December 11, 2016
Lecture 7 Logistic Regression Luigi Freda ALCOR Lab DIAG University of Rome La Sapienza December 11, 2016 Luigi Freda ( La Sapienza University) Lecture 7 December 11, 2016 1 / 39 Outline 1 Intro Logistic
More informationSquare Root Raised Cosine Filter
Wireless Information Transmission System Lab. Square Root Raised Cosine Filter Institute of Communications Engineering National Sun Yat-sen University Introduction We consider the problem of signal design
More informationExtended Introduction to Computer Science CS1001.py. Lecture 8 part A: Finding Zeroes of Real Functions: Newton Raphson Iteration
Extended Introduction to Computer Science CS1001.py Lecture 8 part A: Finding Zeroes of Real Functions: Newton Raphson Iteration Instructors: Benny Chor, Amir Rubinstein Teaching Assistants: Yael Baran,
More informationQueens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 2018 3 Lecture 3 3.1 General remarks March 4, 2018 This
More informationFourier Transform and Frequency Domain
Fourier Transform and Frequency Domain http://graphics.cs.cmu.edu/courses/15-463 15-463, 15-663, 15-862 Computational Photography Fall 2017, Lecture 6 Course announcements Last call for responses to Doodle
More informationWarm-Up. Simplify the following terms:
Warm-Up Simplify the following terms: 81 40 20 i 3 i 16 i 82 TEST Our Ch. 9 Test will be on 5/29/14 Complex Number Operations Learning Targets Adding Complex Numbers Multiplying Complex Numbers Rules for
More informationISM206 Lecture Optimization of Nonlinear Objective with Linear Constraints
ISM206 Lecture Optimization of Nonlinear Objective with Linear Constraints Instructor: Prof. Kevin Ross Scribe: Nitish John October 18, 2011 1 The Basic Goal The main idea is to transform a given constrained
More informationESS Dirac Comb and Flavors of Fourier Transforms
6. Dirac Comb and Flavors of Fourier ransforms Consider a periodic function that comprises pulses of amplitude A and duration τ spaced a time apart. We can define it over one period as y(t) = A, τ / 2
More informationAPPENDIX A. The Fourier integral theorem
APPENDIX A The Fourier integral theorem In equation (1.7) of Section 1.3 we gave a description of a signal defined on an infinite range in the form of a double integral, with no explanation as to how that
More informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Filtering in the Frequency Domain http://www.ee.unlv.edu/~b1morris/ecg782/ 2 Outline Background
More informationConjugate gradient algorithm for training neural networks
. Introduction Recall that in the steepest-descent neural network training algorithm, consecutive line-search directions are orthogonal, such that, where, gwt [ ( + ) ] denotes E[ w( t + ) ], the gradient
More information2. If the values for f(x) can be made as close as we like to L by choosing arbitrarily large. lim
Limits at Infinity and Horizontal Asymptotes As we prepare to practice graphing functions, we should consider one last piece of information about a function that will be helpful in drawing its graph the
More informationL29: Fourier analysis
L29: Fourier analysis Introduction The discrete Fourier Transform (DFT) The DFT matrix The Fast Fourier Transform (FFT) The Short-time Fourier Transform (STFT) Fourier Descriptors CSCE 666 Pattern Analysis
More information19. The Fourier Transform in optics
19. 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 informationE2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2)
E.5 Signals & Linear Systems Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & ) 1. Sketch each of the following continuous-time signals, specify if the signal is periodic/non-periodic,
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 informationTopic 3: Fourier Series (FS)
ELEC264: Signals And Systems Topic 3: Fourier Series (FS) o o o o Introduction to frequency analysis of signals CT FS Fourier series of CT periodic signals Signal Symmetry and CT Fourier Series Properties
More informationLecture 3 September 1
STAT 383C: Statistical Modeling I Fall 2016 Lecture 3 September 1 Lecturer: Purnamrita Sarkar Scribe: Giorgio Paulon, Carlos Zanini Disclaimer: These scribe notes have been slightly proofread and may have
More informationNumerical optimization
Numerical optimization Lecture 4 Alexander & Michael Bronstein tosca.cs.technion.ac.il/book Numerical geometry of non-rigid shapes Stanford University, Winter 2009 2 Longest Slowest Shortest Minimal Maximal
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 information