理想带通使用如下函数,dim是维度,如果input是一维,则dim=1
% FILTERED = ideal_bandpassing(INPUT,DIM,WL,WH,SAMPLINGRATE)
%
% Apply ideal band pass filter on INPUT along dimension DIM.
%
% WL: lower cutoff frequency of ideal band pass filter
% WH: higher cutoff frequency of ideal band pass filter
% SAMPLINGRATE: sampling rate of INPUT
%
% Copyright (c) 2011-2012 Massachusetts Institute of Technology,
% Quanta Research Cambridge, Inc.
%
% Authors: Hao-yu Wu, Michael Rubinstein, Eugene Shih,
% License: Please refer to the LICENCE file
% Date: June 2012
%
function filtered = ideal_bandpassing(input, dim, wl, wh, samplingRate)
if (dim > size(size(input),2))
error('Exceed maximum dimension');
end
input_shifted = shiftdim(input,dim-1);
Dimensions = size(input_shifted);
n = Dimensions(1);
dn = size(Dimensions,2);
Freq = 1:n;
Fre