Files for Gravitational Wave Astronomy Summer School 2012
Peter S. Shawhan

Monday

• Lecture slides
• Handout with Matlab exercises
• Matlab solutions to exercises - see commands and comments in script
• Gary LaMotte's Mathematica solutions: Parts 1 and 2, Part 3, Parts 4 and 5

Tuesday

• Lecture slides
• Handout with Matlab exercises
• Files for exercises: day2data.txt, day2template.txt
• If you have Matlab or Octave but don't have the Signal Processing Toolbox, you are probably missing the pwelch, butter, and filter commands.
  • Here is a simple replacement for pwelch: peterwelch.m. It just takes one argument, the time series.
  • Here is a function which applies a Butterworth filter between 35 and 250 Hz: butterbp35_250.m. It takes one input argument, the time series to be filtered, and returns the filtered time series.
• Matlab solutions to exercises - see commands and comments in script
• Gary LaMotte's Mathematica solutions: All parts

• Lecture slides
• Handout with Matlab exercises
• File for exercises: day3asd.txt * This file contains the amplitude spectral density; square it element-by-element to get the PSD.
• If you have Matlab or Octave but don't have the Signal Processing Toolbox:
  • Here is a function which applies a Butterworth low-pass filter with a corner frequency of 512 Hz: butterlp512.m. It takes one input argument, the time series to be filtered, and returns the filtered time series.
• Matlab solutions to exercises
• Gary LaMotte's Mathematica solutions: All parts

• Lecture slides
• Slides from Paul Demorest about GW detection by pulsar timing
• Statistics article by Fred James
• Handout with Matlab exercises
• Solutions to exercises

Friday

• Lecture slides
• Handout with Matlab exercises
• Files for exercises: day5dataA.txt, day5dataB.txt, day5dataC.txt, day5dataD.txt
• If you have Matlab or Octave but don't have the Signal Processing Toolbox, you are probably missing the butter command, although I think (?) you would have the filter command.
  • Here is a function which can design (but not apply) a Butterworth filter with various frequency bands that you need for the exercise today: day5butter.m. It takes two input arguments, the low and high frequencies for the band in Hz (not in fractions of the Nyquist frequency like the regular butter function expects!), and returns the B and A coefficients for the filter. Then you have to call filter with those vectors to do the actual filtering.
• Matlab solutions to exercises
• Gary LaMotte's Mathematica solutions: Part 1 data A, Part 1 data B, Part 2, Part 3