view arrows/dsp/amelceps.m @ 61:eff6bddf82e3 tip

Finally implemented perceptual brightness thing.

% melceps - arrow from audio to log mel cepstrum (MFCCs)
function o=melceps(N,M,L,fs,varargin)
	opts=options('norm_rate',1e-6,'plot_melspec',1,'lowf',100,varargin{:});
	if opts.plot_melspec, 
		mp=obs_with(imgtrace(256,'fig',9,'clim',[-20,20])); 
		figure(9); colormap(bipolar1);
	else; mp=aid; end

	o = amelspec(N,M,fs,'lowf',opts.lowf) ...
	  * states(azeromean(gaussian,'rate',opts.norm_rate),50) * (aid + vplotter('fig',12)) ... 
	  * mp * adct(1:L); % ...
%	  * obs_with(plotter('fig',11,'ylim',[-20,20],'xlim',[0,L+1],'plotfn',@bar));