Peter S. Shawhan

- 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

- 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.

- Here is a simple replacement for
- 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.

- Here is a function which applies a Butterworth low-pass filter with a corner frequency of 512 Hz:
- 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

- 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.

- Here is a function which can design (but not apply) a Butterworth filter with various frequency bands that you need for the exercise today:
- Matlab solutions to exercises
- Gary LaMotte's Mathematica solutions: Part 1 data A, Part 1 data B, Part 2, Part 3