% fig3a
% Figure 3(a) of "Cancellation model of pitch perception" (JASA 97)

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

y = synth_complex1;

% 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);
yyy = acf2(yy, range*sr, windowsize*sr);
yyy = yyy / max(max(yyy));

% hack to suppress stupid lines
yyy(size(yyy, 1), 1) = -0.05;

plot2(yyy);

%set(get(gca, 'title'), 'String',...
%	'"Cancellation model of pitch perception" - Fig. 3a');
set(gca, 'xgrid', 'off');
set(gca, 'ygrid', 'off');
set(gca, 'zgrid', 'off');
set (get(gca, 'xlabel'), 'fontsize', 16);
set (get(gca, 'ylabel'), 'fontsize', 16);
set (gca, 'fontsize', 16);

set (gca, 'zlim', [-0.05, 1]);