% 

global sr;
sr = 10000;
level = 70;			% dB
nchan = 40;			% number of BM channels
flo = 100;			% Hz - lowest BM channel frequency (fhi = sr/2)
range = 0.03; 	% s	- range for ACF
windowsize = 0.025;	% s	- window size for ACF

y = synth_ramp;

% scale
level = level - 30;	% 30 dB corresponds to RMS=1 in M&H's model.
y = rms(y, fromdB(level));

yy = bm_haircell(y, nchan, flo);


subplot(2,1,1);
yyy = acf2(yy, range*sr, windowsize*sr);
plot2(yyy);
set(get(gca, 'title'), 'String',...
	'ACF');

subplot(2,1,2);
yyy = rdf2(yy, range*sr, windowsize*sr);
plot2(yyy);
set(get(gca, 'title'), 'String',...
	'RDF');