% 

global sr;
sr = 10000;
level = 90;			% dB
nchan = 40;			% number of BM channels
flo = 100;			% Hz - lowest BM channel frequency (fhi = sr/2)
range = 0.03; 		% s
windowsize = 0.025;		% s
T = 0.01;			%s delay for rippled noise

y = synth_ripple_minus(T);

% 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);
plot1(sum(yyy));
% start ordinate at 0
ylim = get (gca, 'Ylim');
ylim(1) = 0;
set (gca, 'Ylim', ylim);
set(get(gca, 'title'), 'String',...
	'SACF');

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