Mercurial > hg > aimmat
annotate aim-mat/tools/@signal/lowpass_2003.m @ 4:537f939baef0 tip
various bug fixes and changed copyright message
author | Stefan Bleeck <bleeck@gmail.com> |
---|---|
date | Tue, 16 Aug 2011 14:37:17 +0100 |
parents | 20ada0af3d7d |
children |
rev | line source |
---|---|
tomwalters@0 | 1 % |
tomwalters@0 | 2 % function |
tomwalters@0 | 3 % |
tomwalters@0 | 4 % INPUT VALUES: |
tomwalters@0 | 5 % |
tomwalters@0 | 6 % |
tomwalters@0 | 7 % RETURN VALUE: |
tomwalters@0 | 8 % |
tomwalters@0 | 9 % |
bleeck@3 | 10 % This external file is included as part of the 'aim-mat' distribution package |
bleeck@3 | 11 % (c) 2011, University of Southampton |
bleeck@3 | 12 % Maintained by Stefan Bleeck (bleeck@gmail.com) |
bleeck@3 | 13 % download of current version is on the soundsoftware site: |
bleeck@3 | 14 % http://code.soundsoftware.ac.uk/projects/aimmat |
bleeck@3 | 15 % documentation and everything is on http://www.acousticscale.org |
bleeck@3 | 16 |
tomwalters@0 | 17 |
tomwalters@0 | 18 |
tomwalters@0 | 19 |
tomwalters@0 | 20 function sig=lowpass_2003(sig,cutoff,order) |
tomwalters@0 | 21 |
tomwalters@0 | 22 |
tomwalters@0 | 23 if nargin<2 |
tomwalters@0 | 24 order=1; |
tomwalters@0 | 25 end |
tomwalters@0 | 26 |
tomwalters@0 | 27 |
tomwalters@0 | 28 sr=getsr(sig); |
tomwalters@0 | 29 vals=getvalues(sig); |
tomwalters@0 | 30 |
tomwalters@0 | 31 [b a] = butter(order,cutoff/(sr/2)); |
tomwalters@0 | 32 |
tomwalters@0 | 33 |
tomwalters@0 | 34 sig_len = length(vals); |
tomwalters@0 | 35 |
tomwalters@0 | 36 |
tomwalters@0 | 37 filt_vals(1:sig_len) = filter(b,a,vals); |
tomwalters@0 | 38 vals=filt_vals; |
tomwalters@0 | 39 |
tomwalters@0 | 40 |
tomwalters@0 | 41 sig=setvalues(sig,vals); |