A6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2013 Lecture 4 For next week be sure to have read Chapter 5 of
|
|
- Robert Hawkins
- 5 years ago
- Views:
Transcription
1 A6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2013 Lecture 4 For next week be sure to have read Chapter 5 of Gregory (Frequentist Statistical Inference) Today: DFT of complex sinusoid (continued)» Detection issues: threshold and false-alarm probability; ROC curves» Examples: SETI, Continuous wave detection of gravitational waves with LIGO Sampling theorem and interpolation Problem Set 2: look for it on the web page tomorrow (Friday Feb 1)
2 Main Points from Lecture 3 Random walks are nonstationary processes that are not ergodic. CLT works for the DFT but does not for the power-spectrum (PS) estimate. The PS is chi-square distributed with two degrees of freedom, i.e. the PDF is a onesided exponential. Smoothing of the PS increases the number of degrees of freedom, making the PDF tend toward a Gaussian in accord with the CLT. Detection and false-alarm probabilities defined for phasor + noise case.
3 White noise RW 1 RW 2 RW 3
4
5
6
7 Random Walk Examples Spinning objects: Earth, neutron stars Steps in torque or spin rate Observable = spin phase Scattered photon propagation (diffusion) Step = mean-free path Observable = propagation time Cosmic-ray propagation in the Galaxy Step = scattering off of small-scale magnetic field variations Observable = `grammage of interaction based on isotopic content (typically ~ 5 g cm -2 ) Orbital perturbations Asteroid belt objects! Near Earth Objects Motions of planetesimals in protoplanetary disks Galactic orbits of stars from gravitational potential granularity (molecular clouds, spiral arms)! diffusion of stellar populations
8 Other Random Walk Examples MCMC: random walk in parameter space Brownian motion of a dust particle Molecular diffusion Diffusion of biological populations Options pricing in financial markets Step = transaction Observable = price Black-Scholes equation = Fokker-Planck equation
9 Back to Phasor + Noise
10 DFT + noise N=1024 P = 3.99 F = 256
11 Spectrum
12 SETI Search for Extraterrestrial Intelligence First search of nearby stars: 1960 by Frank Drake Typical assumed signal type (radio): monochromatic carrier signal + noise Other cases: sharp pulses (radio, OIR)
13
14
15 Gravitational Waves LIGO = Laser Interferometer Gravitational Wave Observatory Measure effective distance between pendula using lasers Many possible signal types: continuous wave, chirped, bursts (CWs from rotating neutron stars with mountains ) CW analysis done as part of the Einstein@Home volunteer network
16 Gravitational Waves CW = Continuous Wave See for details and image credits
17 Gravitational Waves Chirped GW
18 Gravitational Waves Burst GW
19 Gravitational Waves Stochastic GW
20 Gravitational Waves
21
22 LIGO detector in Hanford, WA; A second detector is in Livingston, Louisiana
23 PDF of Phasor Magnitude (Rice distribution) s = 0, 3, 5, 10 sigma_n = 1
24 PDF of Intensity s = 0, 1, 3, 5 sigma_n = 1
25 ROC Curves Receiver Operating Characteristics In a so-called detection problem, we try to establish whether a signal of some assumed type is present in data that include noise This is a universal problem that applies to many laboratory and observational contexts. In astronomy, ROC curves apply to finding sources/signals in images, spectra, time series, etc. An ROC curve = P d vs P fa (detection vs falsealarm probability)
26 ROC Curve Examples
27 LIGO noise! white noise
28 A B5-/6,5-CD%E,F## A )GHFII,5-/6,5-JH)7/J1>EJ,31I# A 9K!L M #>%2+3>53,#&%+1-6,,2/#H2%;,//#NO4?#3*6*# A NO4?#3*6*#/,2&,3#<2%E#PJ8J#QPJ#R5/;J#S5+>*1T,,U# A $%>#B5-/6,5-CD%E,#&%+1-6,,2/#6%#H2%;,//#V2,;5@%# 0VNWV#3*6*# A 2*>#3*6*#(%2-,++#XY#V+@,26#B5-/6,5-#O-/Z616,#QD*--%&,2U# A [,35/H,2/,3#ZE,#/,25,/#6%#,*;)#BCD#;+5,-6# A 8H,-3#\!9#)2#H2%;,//5-.#,*;)#ZE,#/,25,/#6%#/,*2;)#<%2#@5-*27# H1+/*2/#>56)#0 %2@# */#/)%26#*/#\!L#E5-# A (*-353*6,#/5.-*+/#/,-6#@*;T#6%#/,2&,2/#*6#D*--%&,2#*-3#(%2-,++]# 3*6*@*/,#*6#(%2-,++# A S,6*X*-*+7/5/#6%#/,+,;6#;*-353*6,/#<%2#2,%@/,2&*Z%-#*6#V2,;5@%#!"#$%&#'!!# ()*++,-.,/#####01+/*2/#*-3#42*&567####80# 9:#
29 Arecibo Pulsar Survey Data Flow Survey analysis + Einstein@Home extended survey analysis 13 Nov '11 Challenges Pulsars and Gravity SP 29
30 Fourier Transforms II This lecture is largely about discreteness effects. DISCRETE IN ONE DOMAIN PERIODIC IN THE CONJUGATE DOMAIN We will consider here periodic extensions and sampled versions of continuous functions f(t). Graphically, we have: f(t) F (f) f p (t) Fp (f) f ps (t) the act of sampling assumes or imposes periodicity in the f domain Fps (f) if the periodic extensions overlap, then there are problems (aliasing) Related Kinds of Fourier Transforms/Series The essence of Fourier analysis is the expansion of a function into a series of sinusoidal functions. The particular form of the expansion depends on the kinds of functions we are expanding. We classify functions according to the DT/CT dichotomy but also according to whether they are periodic/aperiodic. The most general case is the aperiodic/ct signal: We have already motivated the Fourier transform relation for these signals x(t) = 1 2π dw X(ω) e iωt (1) X(ω) = dt x(t) e iωt (2) 1
31 Note that ω is the radian frequency which is related to cyclic frequency f by ω =2πf. Defining the F.T. in terms of f yields x(t) = df X(f) e 2πift (3) X(f) = dt x(t) e 2πift (4) Note disappearance of (2π) 1 in front. Some definitions of the ω transform place a factor of 1 2π in front. Now, from the aperiodic /CT transform relations, we can derive Fourier series for periodic /CT functions and the discrete Fourier transform for DT functions. Kinds of Fourier Transforms: Fourier transform Fourier Series DFT function f(t) f p (t) f ps (t) periodic extention sampled ( t) of f(t), period = and periodic extention transform F (f) ak = 1 F ( k ) a k continuous? continuous, continuous, discrete, (fn, transform) continuous discrete discrete periodic? aperiodic, periodic, periodic, (fn, transform) aperiodic aperiodic periodic 2
32 Useful Quantities: We will use the following FT relations for delta functions: 1 δ(f) δ(t) 1 where implies a Fourier transform relation. It is useful to consider a periodic train of impluses s(t) = δ(t n ) n= Shift theorem If, f(t) F (f) then, f(t t 0 ) e 2πift0 F (f) This follows by inspection of the definition for F : F.T. {f(t t 0 )} = dt f(t t 0 ) e 2πift ; change variable to t = t t 0 t = t + t 0 : F.T. = dt f(t ) e 2πift e 2πift0 F (f) Similarly, if f(t) F (f) then e +2πif0t f(t) F (f f 0 ) This is the basis for heterodyned receiver systems. 3
33 Sampling Function: Using the first form of the shift theorem we have so, δ(t n ) e 2πifn S(f) = n= e 2πifn. How does this behave? As a series of delta functions in f. Note for f = integer S(f) is the integral of a sinusoid (many cycles) 0. However, for f = k =0, ±1,..., S(f) Thus, we can write (guessing) S(f) = 1 k= δ(f k ) We introduce the factor 1 to keep dimensionally correct, note S has units of s time and s has units of (time 1 ) S is dimensionless. Note that s(t) has a Fourier transform S(f) that is of the same form as s(t) (a sequence of delta functions). The only other function I know of that has this property is a Gaussian function: e πt2 e πf 2 In the following we will use the sampling function several times. For convenience, we will label it with the period : S(t, ) n δ(t n ). 4
34 FOURIER SERIES Consider a function f(t) F (f) which is time limited: f(t) =0outside some range (t a,t b ). Now convolve f(t) with the impulse train and find the resultant fourier transform. This is trivial because of the convolution theorem. f p (t) f(t) S(t, ) = f(t) n n= f(t n ) δ(t n ) F (f) S(f) = 1 = 1 k k δ f k F (f) k F δ f k = periodic function with period ; has discrete Fourier transform Thus, or (simplifying further) f p (t) = f p (t) 1 k k F δ f k df F p (f) e 2πift (5) = = 1 f p (t) = k df 1 k k k F δ f k e 2πift (6) F ( k ) e2πi( k )t (7) a k e 2πi( k )t, a k F ( k ) (8) 5
35 Note that the F.S. is a series in orthonormal basis functions g k (t) =e 2πi( k )t and coefficients a k. We can show that : /2 /2 dt g k g k = = = /2 /2 /2 /2 dt g k (t) g k (t) =δ kk dt e 2πi(k k )t/ 2πi(k k ) [e2πi(k k )/2 e 2πi(k k )/2 ] π(k k ) sin π(k k )= 1 k = k 0 otherwise Thus we have a Fourier series representation of the periodic signal f p (t) in terms of the Fourier transform of the aperiodic function f(t). We can write this in the form of a cos and sin series by expanding the exponential and using a k = b i + ic k f p (t) = k a k e 2πi( k )t = k (b k + ic k )(cos 2π k t + i sin 2π k t) Graphically, we have f(t) F (f) f p (t) Fp (f) periodicity in the time domain discreteness in the frequency domain By reciprocity we can show that discreteness in the time domain periodicity in the f domain 6
36 DISCRETE FOURIER TRANSFORM A periodic (in t) function can be described by a set of discrete coefficients in a Fourier series. By discretizing the time domain function (still periodic) we can show that a discrete set of time-domain sample describes the resultant periodic and sampled function: 1. Start with f(t) =CT/aperiodic function. 2. Extend it periodically, as before: and write as a Fourier series: f p (t) =f(t) s(t, ) f p (t) = k now CT/periodic a k e 2πi( k )t 3. Now sample f p (t) to get a DT/periodic function: f ps (t) f p (t) S(t, t) =f p (t) n δ(t n t) = k a k e 2πi k t n δ(t n t) k 2πi( 4. Now multiply f ps by e )t and integrate over one period i.e. apply 0 k dt e 2πi( )t to the equation in 3.: 0 k 2πi dt e t f p (t) δ(t n k 2πi t) = dt e t k 2πi a k e t δ(t n t), n 0 k n We integrate over one period so f p (t) =f(t) in this interval. Therefore, LHS = f(n t) e 2πi k n t n :0 n t<. Now RHS = k a k dt e 2πi (k k )t δ(t n t) = n 0 k =e 2πi (k k )n t e 2πi(k k )n t n a k 7
37 Note: {n :0 n t< } or {n :0 n <N = no. samples in interval (assume an integer)} t thus {n :0 n N 1} 0 k = k sum of integer N 1 But e 2πi (k k )n t cycles of sinusoids = n =0 1=N k = k n Nδ kk Kronecker delta Therefore RHS = k a k N δ kk = Na k Thus (using N t ) Na k = n :0 n t< f(n t) e 2πi k n N N 1 a k = N 1 f(n t) e 2πi k n N, k =0,...,N 1 ( ) n 0 In words: The Fourier series coefficients a k for a periodic and sampled function are represented as discrete sums of samples over one period. Note that a k is manifestly periodic: N 1 a k +N N 1 f(n t) e 2πin (k +N)/N n =0 e 2πin k /N e 2πin 1 8
38 N 1 5. Now find the inverse transform: Apply e 2πi k n k =0 N to eqn. : N 1 k =0 a k e 2πi k n N N 1 = N 1 f(n t) n =0 N 1 k =0 e 2πik (n n ) N N δ nn N 1 n f(n t) N δ nn = f(n t) or f(n t) = N 1 k =0 a k e 2πi k n N Thus the samples of a periodic, sampled function can be written as a discrete finite Fourier series of the F.S. coefficients. Eqn. and constitute the discrete Fourier transform. Notation: can drop t, rename a k so that we have DFT : f n = k F k e 2πink/N F k = N 1 n f n e 2πink/N Note that the function is periodic in both domains: f n+n = f n, F k+n = F k 9
39 SAMPLING THEOREM Consider a function that has finite frequency extent, e.g. a function that is band limited in the interval [ f, f]. Then the function can be determined totally by sampling it in the time domain at the rate 2 f. Samples must occur at intervals t 1 2 f Plausibility: consider a sine wave of frequency f f(t) = A cos 2π ft A 2 (e2πi ft + e 2πi ft ) (9) F (f) = A 2 [δ(f f)+δ(f + f)] (10) The sampling theorem says to sample this a rate of at least This requires at least two samples per period. f s =2 f or at intervals t = 1 2 f. One might be unlucky if the samples occur at the zero crossings, but this is ok because, by definition, there is no energy in the sine wave (A =0)for a bandlimited function. 10
40 Logic for derivation of sampling theorem: I. First impose periodicity in the f domain and then derive interpolation formula for critically sampled data. II. Then sample in the t domain at arbitrary rate and derive the resultant interpolation formula. Demonstration of sampling theorem (applies to bandlimited signals) Consider a bandlimited function f(t) F (f), F (f) =0for f f f(t) is CT/aperiodic f(t) = df F (f) e 2πift = f f df F (f) e 2πift our choice here with W =2 f critical sampling Now define a periodic version of F (f) called F p (f): Since F p (f) is periodic, by the duality theorem [if f(t) F ( f)], then F (t) f( f)] we can write it as a Fourier series: period W =2 f (11) F p (f) = C m e 2πi m W f (12) m= where the Fourier coefficients are (as usual) C m = 1 W = 1 2 f W/2 W/2 f f df F p (f) e +2πi( m W )f (13) df F (f) e +2πi m 2 f f (14) 11
41 By comparison with the expression for f(t) (and from our knowledge of how F.S. coefficients relate to the F.T.), we have C m = 1 2 f f(t) t= m = 1 m 2 f 2 f f 2 f Thus, F p (f) = 1 2 f m= m f 2 f e 2πi m 2 f f Consequently, since F (f) = F p (f) in the range f f this result implies that F (f) can be determined by sampling f(t) at intervals 1 2 f. Then, from the fact that f(t) F (f) it follows that the t domain function itself can be determined at all times from a set of discrete samples. The interpolation formula shows that the continuous frequency-domain function can be calculated from discrete data. 12
42
43
44
45
46
A6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring
A6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2015 http://www.astro.cornell.edu/~cordes/a6523 Lecture 4 See web page later tomorrow Searching for Monochromatic Signals
More informationA6523 Linear, Shift-invariant Systems and Fourier Transforms
A6523 Linear, Shift-invariant Systems and Fourier Transforms Linear systems underly much of what happens in nature and are used in instrumentation to make measurements of various kinds. We will define
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 informationA6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring
A653 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 15 http://www.astro.cornell.edu/~cordes/a653 Lecture 3 Power spectrum issues Frequentist approach Bayesian approach (some
More informationA6523 Modeling, Inference, and Mining Jim Cordes, Cornell University. False Positives in Fourier Spectra. For N = DFT length: Lecture 5 Reading
A6523 Modeling, Inference, and Mining Jim Cordes, Cornell University Lecture 5 Reading Notes on web page Stochas
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 1 Organization:» Syllabus (text, requirements, topics)» Course approach (goals, themes) Book: Gregory, Bayesian
More informationGravitational-Wave Data Analysis: Lecture 2
Gravitational-Wave Data Analysis: Lecture 2 Peter S. Shawhan Gravitational Wave Astronomy Summer School May 29, 2012 Outline for Today Matched filtering in the time domain Matched filtering in the frequency
More informationLIGO Detection of Gravitational Waves. Dr. Stephen Ng
LIGO Detection of Gravitational Waves Dr. Stephen Ng Gravitational Waves Predicted by Einstein s general relativity in 1916 Indirect confirmation with binary pulsar PSR B1913+16 (1993 Nobel prize in physics)
More informationLIGO Status and Advanced LIGO Plans. Barry C Barish OSTP 1-Dec-04
LIGO Status and Advanced LIGO Plans Barry C Barish OSTP 1-Dec-04 Science Goals Physics» Direct verification of the most relativistic prediction of general relativity» Detailed tests of properties of gravitational
More informationSearching for gravitational waves. with LIGO detectors
Werner Berger, ZIB, AEI, CCT Searching for gravitational waves LIGO Hanford with LIGO detectors Gabriela González Louisiana State University On behalf of the LIGO Scientific Collaboration KITP Colloquium,
More informationStatus and Prospects for LIGO
Status and Prospects for LIGO Crab Pulsar St Thomas, Virgin Islands Barry C. Barish Caltech 17-March-06 LIGO Livingston, Louisiana 4 km 17-March-06 Confronting Gravity - St Thomas 2 LIGO Hanford Washington
More informationGravitational Wave Astronomy the sound of spacetime. Marc Favata Kavli Institute for Theoretical Physics
Gravitational Wave Astronomy the sound of spacetime Marc Favata Kavli Institute for Theoretical Physics What are gravitational waves? Oscillations in the gravitational field ripples in the curvature of
More informationSearching for gravitational waves from neutron stars
Searching for gravitational waves from neutron stars Ian Jones D.I.Jones@soton.ac.uk General Relativity Group, Southampton University Ian Jones Searching for gravitational waves from neutron stars 1/23
More informationGravitational wave data analysis
Max Planck Institut für Gravitationsphysik Albert Einstein Institut, Germany Pasadena, June 2011 1 Introduction to gravitational waves 2 3 4 5 6 Gravitational Waves GR can be formulated in terms of a spacetime
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 informationSearching for Ripples in Space-Time with Your Home Computer"
Einstein@Home: Searching for Ripples in Space-Time with Your Home Computer Eric Myers LIGO Hanford Observatory! Hanford, Washington! Amateur Astronomers Association of New York 11 April 2008 Overview q
More informationGravitational Waves & Intermediate Mass Black Holes. Lee Samuel Finn Center for Gravitational Wave Physics
Gravitational Waves & Intermediate Mass Black Holes Lee Samuel Finn Center for Gravitational Wave Physics Outline What are gravitational waves? How are they produced? How are they detected? Gravitational
More informationVariable and Periodic Signals in Astronomy
Lecture 14: Variability and Periodicity Outline 1 Variable and Periodic Signals in Astronomy 2 Lomb-Scarle diagrams 3 Phase dispersion minimisation 4 Kolmogorov-Smirnov tests 5 Fourier Analysis Christoph
More informationReview of Discrete-Time System
Review of Discrete-Time System Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides were based on class notes developed by Profs. K.J. Ray Liu and Min Wu.
More informationGravitational Waves and LIGO
Gravitational Waves and LIGO Ray Frey, University of Oregon 1. GW Physics and Astrophysics 2. How to detect GWs The experimental challenge 3. Prospects June 16, 2004 R. Frey QNet 1 General Relativity Some
More informationGravitational-Wave Data Analysis
Gravitational-Wave Data Analysis Peter Shawhan Physics 798G April 12, 2007 Outline Gravitational-wave data General data analysis principles Specific data analysis methods Classification of signals Methods
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 informationThe search for continuous gravitational waves: analyses from LIGO s second science run
The search for continuous gravitational waves: analyses from LIGO s second science run Michael Landry LIGO Hanford Observatory on behalf of the LIGO Scientific Collaboration http://www.ligo.org April APS
More informationLIGO Observational Results
LIGO Observational Results Patrick Brady University of Wisconsin Milwaukee on behalf of LIGO Scientific Collaboration LIGO Science Goals Direct verification of two dramatic predictions of Einstein s general
More informationLecture 28 Continuous-Time Fourier Transform 2
Lecture 28 Continuous-Time Fourier Transform 2 Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/6/14 1 Limit of the Fourier Series Rewrite (11.9) and (11.10) as As, the fundamental
More informationGravity. Newtonian gravity: F = G M1 M2/r 2
Gravity Einstein s General theory of relativity : Gravity is a manifestation of curvature of 4- dimensional (3 space + 1 time) space-time produced by matter (metric equation? g μν = η μν ) If the curvature
More informationSignal Modeling, Statistical Inference and Data Mining in Astrophysics
ASTRONOMY 6523 Spring 2013 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Course Approach The philosophy of the course reflects that of the instructor, who takes a dualistic view
More informationReview of Analog Signal Analysis
Review of Analog Signal Analysis Chapter Intended Learning Outcomes: (i) Review of Fourier series which is used to analyze continuous-time periodic signals (ii) Review of Fourier transform which is used
More informationVariable and Periodic Signals in Astronomy
Lecture 14: Variability and Periodicity Outline 1 Variable and Periodic Signals in Astronomy 2 Lomb-Scarle diagrams 3 Phase dispersion minimisation 4 Kolmogorov-Smirnov tests 5 Fourier Analysis Christoph
More informationThe LIGO Project: a Status Report
The LIGO Project: a Status Report LIGO Hanford Observatory LIGO Livingston Observatory Laura Cadonati LIGO Laboratory, MIT for the LIGO Scientific Collaboration Conference on Gravitational Wave Sources
More informationGravitational wave detection with Virgo and LIGO experiment - Case of the long bursts
Gravitational wave detection with Virgo and LIGO experiment - Case of the long bursts Samuel Franco Supervisor: Patrice Hello Laboratoire de l Accélérateur Linéaire (Orsay) 06/12/2013 Samuel Franco (LAL)
More informationGW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral Lazzaro Claudia for the LIGO Scientific Collaboration and the Virgo Collaboration 25 October 2017 GW170817 PhysRevLett.119.161101
More informationNew Applications of Sparse Methods in Physics. Ra Inta, Centre for Gravitational Physics, The Australian National University
New Applications of Sparse Methods in Physics Ra Inta, Centre for Gravitational Physics, The Australian National University 2 Sparse methods A vector is S-sparse if it has at most S non-zero coefficients.
More informationOverview Ground-based Interferometers. Barry Barish Caltech Amaldi-6 20-June-05
Overview Ground-based Interferometers Barry Barish Caltech Amaldi-6 20-June-05 TAMA Japan 300m Interferometer Detectors LIGO Louisiana 4000m Virgo Italy 3000m AIGO Australia future GEO Germany 600m LIGO
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 informationA6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2011
A6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2011 Reading: Chapter 10 = linear LSQ with Gaussian errors Chapter 11 = Nonlinear fitting Chapter 12 = Markov Chain Monte
More informationIntroduction to Fourier Transforms. Lecture 7 ELE 301: Signals and Systems. Fourier Series. Rect Example
Introduction to Fourier ransforms Lecture 7 ELE 3: Signals and Systems Fourier transform as a limit of the Fourier series Inverse Fourier transform: he Fourier integral theorem Prof. Paul Cuff Princeton
More informationGravitational waves from the early Universe
Gravitational waves from the early Universe Part 1 Sachiko Kuroyanagi (Nagoya University) 26 Aug 2017 Summer Institute 2017 What is a Gravitational Wave? What is a Gravitational Wave? 11 Feb 2016 We have
More informationProbing the Universe for Gravitational Waves
Probing the Universe for Gravitational Waves Barry C. Barish Caltech Crab Pulsar University of Illinois 16-Feb-06 General Relativity the essential idea G μν = 8πΤ μν Gravity is not a force, but a property
More informationContinuous Wave Data Analysis: Fully Coherent Methods
Continuous Wave Data Analysis: Fully Coherent Methods John T. Whelan School of Gravitational Waves, Warsaw, 3 July 5 Contents Signal Model. GWs from rotating neutron star............ Exercise: JKS decomposition............
More information8 The Discrete Fourier Transform (DFT)
8 The Discrete Fourier Transform (DFT) ² Discrete-Time Fourier Transform and Z-transform are de ned over in niteduration sequence. Both transforms are functions of continuous variables (ω and z). For nite-duration
More informationNewtonian instantaneous action at a distance General Relativity information carried by gravitational radiation at the speed of light
Modern View of Gravitation Newtonian instantaneous action at a distance G µ = 8 µ # General Relativity information carried by gravitational radiation at the speed of light Gravitational Waves GR predicts
More informationProspects for joint transient searches with LOFAR and the LSC/Virgo gravitational wave interferometers
Prospects for joint transient searches with LOFAR and the LSC/Virgo gravitational wave interferometers Ed Daw - University of Sheffield On behalf of the LIGO Scientific Collaboration and the Virgo collaboration
More informationGravity -- Studying the Fabric of the Universe
Gravity -- Studying the Fabric of the Universe Barry C. Barish Caltech "Colliding Black Holes" Credit: National Center for Supercomputing Applications (NCSA) AAAS Annual Meeting Denver, Colorado 17-Feb-03
More informationChapter 8 The Discrete Fourier Transform
Chapter 8 The Discrete Fourier Transform Introduction Representation of periodic sequences: the discrete Fourier series Properties of the DFS The Fourier transform of periodic signals Sampling the Fourier
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 informationGravitational Wave Astronomy using 0.1Hz space laser interferometer. Takashi Nakamura GWDAW-8 Milwaukee 2003/12/17 1
Gravitational Wave Astronomy using 0.1Hz space laser interferometer Takashi Nakamura GWDAW-8 Milwaukee 2003/12/17 1 In 2001 we considered what we can do using 0.1 hertz laser interferometer ( Seto, Kawamura
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 informationA6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring
Lecture 8 A6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2015 http://www.astro.cornell.edu/~cordes/a6523 Applications: Bayesian inference: overview and examples Introduction
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 informationStatus of LIGO. David Shoemaker LISA Symposium 13 July 2004 LIGO-G M
Status of LIGO David Shoemaker LISA Symposium 13 July 2004 Ground-based interferometric gravitational-wave detectors Search for GWs above lower frequency limit imposed by gravity gradients» Might go as
More informationProbing for Gravitational Waves
Probing for Gravitational Waves LIGO Reach with LIGO AdLIGO Initial LIGO Barry C. Barish Caltech YKIS2005 Kyoto University 1-July-05 Einstein s Theory of Gravitation a necessary consequence of Special
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 informationProbing the Universe for Gravitational Waves
Probing the Universe for Gravitational Waves "Colliding Black Holes" Credit: National Center for Supercomputing Applications (NCSA) Barry C. Barish Caltech Argonne National Laboratory 16-Jan-04 LIGO-G030523-00-M
More informationReview: Continuous Fourier Transform
Review: Continuous Fourier Transform Review: convolution x t h t = x τ h(t τ)dτ Convolution in time domain Derivation Convolution Property Interchange the order of integrals Let Convolution Property By
More information+56'0+0)614#8+6#6+10#.#8'5U +056'+0510).+0'5(41/6*'0+8'45' Barry C. Barish
+56'0+0)614#8+6#6+10#.#8'5U +056'+0510).+0'5(41/6*'0+8'45' Barry C. Barish .$'46 +056'+0 +45##%'9610 Perhaps the most important scientist of all time! Invented the scientific method in Principia Greatest
More informationThe direct detection of gravitational waves: The first discovery, and what the future might bring
The direct detection of gravitational waves: The first discovery, and what the future might bring Chris Van Den Broeck Nikhef - National Institute for Subatomic Physics Amsterdam, The Netherlands Physics
More information6.003 Homework #10 Solutions
6.3 Homework # Solutions Problems. DT Fourier Series Determine the Fourier Series coefficients for each of the following DT signals, which are periodic in N = 8. x [n] / n x [n] n x 3 [n] n x 4 [n] / n
More informationLecture 7 ELE 301: Signals and Systems
Lecture 7 ELE 30: Signals and Systems Prof. Paul Cuff Princeton University Fall 20-2 Cuff (Lecture 7) ELE 30: Signals and Systems Fall 20-2 / 22 Introduction to Fourier Transforms Fourier transform as
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 informationLINEAR SYSTEMS. J. Elder PSYC 6256 Principles of Neural Coding
LINEAR SYSTEMS Linear Systems 2 Neural coding and cognitive neuroscience in general concerns input-output relationships. Inputs Light intensity Pre-synaptic action potentials Number of items in display
More informationSinc Functions. Continuous-Time Rectangular Pulse
Sinc Functions The Cooper Union Department of Electrical Engineering ECE114 Digital Signal Processing Lecture Notes: Sinc Functions and Sampling Theory October 7, 2011 A rectangular pulse in time/frequency
More informationA3. Frequency Representation of Continuous Time and Discrete Time Signals
A3. Frequency Representation of Continuous Time and Discrete Time Signals Objectives Define the magnitude and phase plots of continuous time sinusoidal signals Extend the magnitude and phase plots to discrete
More informationGravitational Wave Detection from the Ground Up
Gravitational Wave Detection from the Ground Up Peter Shawhan (University of Maryland) for the LIGO Scientific Collaboration LIGO-G080393-00-Z From Simple Beginnings Joe Weber circa 1969 AIP Emilio Segre
More informationDiscrete-Time Fourier Transform
Discrete-Time Fourier Transform Chapter Intended Learning Outcomes: (i) (ii) (iii) Represent discrete-time signals using discrete-time Fourier transform Understand the properties of discrete-time Fourier
More informationBME 50500: Image and Signal Processing in Biomedicine. Lecture 2: Discrete Fourier Transform CCNY
1 Lucas Parra, CCNY BME 50500: Image and Signal Processing in Biomedicine Lecture 2: Discrete Fourier Transform Lucas C. Parra Biomedical Engineering Department CCNY http://bme.ccny.cuny.edu/faculty/parra/teaching/signal-and-image/
More informationThe Present Gravitational Wave Detection Effort
The Present Gravitational Wave Detection Effort Keith Riles University of Michigan LIGO Scientific Collaboration International Conference on Topics in Astroparticle and Underground Physics Rome July 1,
More informationLIGO Present and Future. Barry Barish Directory of the LIGO Laboratory
LIGO Present and Future Barry Barish Directory of the LIGO Laboratory LIGO I Schedule and Plan LIGO I has been built by LIGO Lab (Caltech & MIT) 1996 Construction Underway (mostly civil) 1997 Facility
More informationInteresting times for low frequency gravitational wave detection
Interesting times for low frequency gravitational wave detection Neil Cornish Montana State University Why Interesting? Why Interesting? Why Interesting? ESA L2 Selection 2013/14 0.1 1000 Tev (20??) Pulsar
More informationQuestion Paper Code : AEC11T02
Hall Ticket No Question Paper Code : AEC11T02 VARDHAMAN COLLEGE OF ENGINEERING (AUTONOMOUS) Affiliated to JNTUH, Hyderabad Four Year B. Tech III Semester Tutorial Question Bank 2013-14 (Regulations: VCE-R11)
More informationMonday, October 21, 2013 Third exam this Friday, October 25. Third Sky Watch watch the weather! Review sheet posted today. Review Thursday, 5 6 PM,
Monday, October 21, 2013 Third exam this Friday, October 25. Third Sky Watch watch the weather! Review sheet posted today. Review Thursday, 5 6 PM, Room WEL 2.256 Lecture 19 posted today Reading: Chapter
More informationContents. Signals as functions (1D, 2D)
Fourier Transform The idea A signal can be interpreted as en electromagnetic wave. This consists of lights of different color, or frequency, that can be split apart usign an optic prism. Each component
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 informationGRAVITATIONAL WAVE SOURCES AND RATES FOR LISA
GRAVITATIONAL WAVE SOURCES AND RATES FOR LISA W. Z. Korth, PHZ6607, Fall 2008 Outline Introduction What is LISA? Gravitational waves Characteristics Detection (LISA design) Sources Stochastic Monochromatic
More informationFourier Series Example
Fourier Series Example Let us compute the Fourier series for the function on the interval [ π,π]. f(x) = x f is an odd function, so the a n are zero, and thus the Fourier series will be of the form f(x)
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 informationWhat can LIGO detect? Abstract
What can LIGO detect? Adam Getchell Physics Department, University of California, Davis, 95616 Abstract This briey reviews the literature on gravitational wave astronomy, including theoretical basis, experimental
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 informationStochastic Processes. A stochastic process is a function of two variables:
Stochastic Processes Stochastic: from Greek stochastikos, proceeding by guesswork, literally, skillful in aiming. A stochastic process is simply a collection of random variables labelled by some parameter:
More information06EC44-Signals and System Chapter Fourier Representation for four Signal Classes
Chapter 5.1 Fourier Representation for four Signal Classes 5.1.1Mathematical Development of Fourier Transform If the period is stretched without limit, the periodic signal no longer remains periodic but
More informationGravity Waves and Black Holes
Gravity Waves and Black Holes Mike Whybray Orwell Astronomical Society (Ipswich) 14 th March 2016 Overview Introduction to Special and General Relativity The nature of Black Holes What to expect when Black
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 informationContents. Signals as functions (1D, 2D)
Fourier Transform The idea A signal can be interpreted as en electromagnetic wave. This consists of lights of different color, or frequency, that can be split apart usign an optic prism. Each component
More informationENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM. Dr. Lim Chee Chin
ENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM Dr. Lim Chee Chin Outline Introduction Discrete Fourier Series Properties of Discrete Fourier Series Time domain aliasing due to frequency
More informationFunctions not limited in time
Functions not limited in time We can segment an arbitrary L 2 function into segments of width T. The mth segment is u m (t) = u(t)rect(t/t m). We then have m 0 u(t) =l.i.m. m0 u m (t) m= m 0 This wors
More informationLIGO: The Laser Interferometer Gravitational Wave Observatory
LIGO: The Laser Interferometer Gravitational Wave Observatory Credit: Werner Benger/ZIB/AEI/CCT-LSU Michael Landry LIGO Hanford Observatory/Caltech for the LIGO Scientific Collaboration (LSC) http://www.ligo.org
More informationTIME SERIES ANALYSIS
2 WE ARE DEALING WITH THE TOUGHEST CASES: TIME SERIES OF UNEQUALLY SPACED AND GAPPED ASTRONOMICAL DATA 3 A PERIODIC SIGNAL Dots: periodic signal with frequency f = 0.123456789 d -1. Dotted line: fit for
More informationVer 3808 E1.10 Fourier Series and Transforms (2014) E1.10 Fourier Series and Transforms. Problem Sheet 1 (Lecture 1)
Ver 88 E. Fourier Series and Transforms 4 Key: [A] easy... [E]hard Questions from RBH textbook: 4., 4.8. E. Fourier Series and Transforms Problem Sheet Lecture. [B] Using the geometric progression formula,
More informationSearches for a Stochastic Gravitational-Wave Background
earches for a tochastic Gravitational-Wave Background John T. Whelan chool of Gravitational Waves, Warsaw, 13 July Contents 1 tochastic GW Backgrounds 1 1.1 Exercise........................ 1. patial Distributions.................
More informationSampling and Discrete Time. Discrete-Time Signal Description. Sinusoids. Sampling and Discrete Time. Sinusoids An Aperiodic Sinusoid.
Sampling and Discrete Time Discrete-Time Signal Description Sampling is the acquisition of the values of a continuous-time signal at discrete points in time. x t discrete-time signal. ( ) is a continuous-time
More informationdiscovers a radio-quiet gamma-ray millisecond Journal Group
Einstein@Home discovers a radio-quiet gamma-ray millisecond pulsar @CHEA Journal Group Contribution of the paper Contribution of the paper Millisecond Pulsars (MSPs) Ver y rapid rotating neutron star
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 informationSearch for compact binary systems in LIGO data
Search for compact binary systems in LIGO data Thomas Cokelaer On behalf of the LIGO Scientific Collaboration Cardiff University, U.K. LIGO-G060630-00-Z Plan 1) Overview What kind of gravitational waves
More informationGravitational Waves modes in Extended Teleparallel Gravity
Gravitational Waves modes in Extended Teleparallel Gravity Salvatore Capozziello based on H. Abedi & S. Capozziello EPJC78(2018)474 Plan of the talk Ø Gravitational waves in General Relativity Ø Extended
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 informationA6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2011
A6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2011 Reading Chapter 5 of Gregory (Frequentist Statistical Inference) Lecture 7 Examples of FT applications Simulating
More informationSensors. Chapter Signal Conditioning
Chapter 2 Sensors his chapter, yet to be written, gives an overview of sensor technology with emphasis on how to model sensors. 2. Signal Conditioning Sensors convert physical measurements into data. Invariably,
More informationGravitational Waves Listening to the Universe. Teviet Creighton LIGO Laboratory California Institute of Technology
Gravitational Waves Listening to the Universe Teviet Creighton LIGO Laboratory California Institute of Technology Summary So far, nearly all our knowledge of the Universe comes from electromagnetic radiation.
More informationWednesday, April 19, 2017
Wednesday, April 19, 2017 Exam 4, Skywatch 4, Friday, April 21. Review Session Tomorrow, 5 6 PM, Welch 2.308. Lectures 21 to 30 Reading for Exam 4: Chapter 8 Neutron Stars - Sections 8.1, 8.2, 8.5, 8.6,
More informationThe LIGO Experiment Present and Future
The LIGO Experiment Present and Future Keith Riles University of Michigan For the LIGO Scientific Collaboration APS Meeting Denver May 1 4, 2004 LIGO-G040239-00-Z What are Gravitational Waves? Gravitational
More information