% Figure 12b of "Cancellation model of pitch perception" (JASA 97)

global sr;
sr = 16000;


range = 0.0075; 	% s	- range for ACF
windowsize = 0.025;	% s	- window size for ACF

D = 0.1;

alpha = 1.3;		% alpha factor in modified neural cancellation filter

y = synth_ramp;

y = hwrect(y);		% half-wave rectify

yy = rdf4(y, range*sr, windowsize*sr, alpha);
yy = yy / max(yy);
plot1(yy);
dump('a', yy);
hold on;

y = synth_damp;

y = hwrect(y);		% half-wave rectify

yy = rdf4(y, range*sr, windowsize*sr, alpha);
yy = yy / max(yy);
plot1_bis(yy, 'r--');
dump('b', yy);
hold off;

set(get(gca, 'XLabel'), 'String', 'time (ms)');
set(get(gca, 'title'), 'String', '"Cancellation model of pitch perception" - Fig. 12b');