melbankm.m

function [xmnmx]=melbankm（pnfsflfhw）  
%MELBANKM determine matrix for a mel-spaced filterbank [XMNMX]=（PNFSFLFHW）  
%  
% Inputs:   p   number of filters in filterbank  
%       n   length of fft  
%       fs  sample rate in Hz  
%       fl  low end of the lowest filter as a fraction of fs （default = 0）  
%       fh  high end of highest filter as a fraction of fs （default = 0.5）  
%       w   any sensible combination of the following:  
%             ‘t‘  triangular shaped filters in mel domain （default）  
%             ‘n‘  hanning shaped filters in mel domain  
%             ‘m‘  hamming shaped filters in mel domain  
%  
%             ‘z‘  highest and lowest filters taper down to zero （default）  
%             ‘y‘  lowest filter remains at 1 down to 0 frequency and  
%              highest filter remains at 1 up to nyquist freqency  
%  
%              If ‘ty‘ or ‘ny‘ is specified the total power in the fft is preserved.  
%  
% Outputs:  x     a sparse matrix containing the filterbank amplitudes  
%             If x is the only output argument then size（x）=[p1+floor（n/2）]  
%             otherwise size（x）=[pmx-mn+1]  
%       mn   the lowest fft bin with a non-zero coefficient  
%       mx   the highest fft bin with a non-zero coefficient  
%  
% Usage:    f=fft（s）;           f=fft（s）;  
%       x=melbankm（pnfs）;     [xnanb]=melbankm（pnfs）;  
%       n2=1+floor（n/2）;        z=log（x*（f（na:nb））.*conj（f（na:nb）））;  
%       z=log（x*abs（f（1:n2））.^2）;  
%       c=dct（z）; c（1）=[];  
%  
% To plot filterbanks e.g.  plot（melbankm（202568000）‘）  
%  
  
  
%      Copyright （C） Mike Brookes 1997  
%  
%      Last modified Tue May 12 16:15:28 1998  
%  
%   VOICEBOX home page: http://www.ee.ic.ac.uk/hp/staff/dmb/voicebox/voicebox.html  
%  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%