band-pass code snippets in matlab