% ACF and SDF global sr; % sampling rate sr = 10000; range = 0.03; % s windowsize = 0.025; % s T = 0.01; %s delay for rippled noise y = synth_ripple_plus(T); subplot(2,1,1); yy = acf2(y, range*sr, windowsize*sr); plot1(yy); set(get(gca, 'title'), 'String',... 'ACF'); subplot(2,1,2); yy = sdf2(y, range*sr, windowsize*sr); plot1(yy); set(get(gca, 'title'), 'String',... 'SDF');