tomwalters@0: %  
tomwalters@0: % function 
tomwalters@0: %
tomwalters@0: %   INPUT VALUES:
tomwalters@0: % 		
tomwalters@0: %
tomwalters@0: %   RETURN VALUE:
tomwalters@0: %		
tomwalters@0: % 
bleeck@3: % This external file is included as part of the 'aim-mat' distribution package
bleeck@3: % (c) 2011, University of Southampton
bleeck@3: % Maintained by Stefan Bleeck (bleeck@gmail.com)
bleeck@3: % download of current version is on the soundsoftware site: 
bleeck@3: % http://code.soundsoftware.ac.uk/projects/aimmat
bleeck@3: % documentation and everything is on http://www.acousticscale.org
bleeck@3: 
tomwalters@0: 
tomwalters@0: 
tomwalters@0: 
tomwalters@0: function sig=lowpass_2003(sig,cutoff,order)
tomwalters@0: 
tomwalters@0: 
tomwalters@0: if nargin<2
tomwalters@0:     order=1;
tomwalters@0: end
tomwalters@0: 
tomwalters@0: 
tomwalters@0: sr=getsr(sig);
tomwalters@0: vals=getvalues(sig);
tomwalters@0: 
tomwalters@0: [b a] = butter(order,cutoff/(sr/2));
tomwalters@0: 
tomwalters@0: 
tomwalters@0: sig_len = length(vals);
tomwalters@0: 
tomwalters@0: 
tomwalters@0: filt_vals(1:sig_len) = filter(b,a,vals);
tomwalters@0: vals=filt_vals;
tomwalters@0: 
tomwalters@0: 
tomwalters@0: sig=setvalues(sig,vals);