% ACF and SDF global sr; % sampling rate sr = 10000; range = 0.012; % s windowsize = 0.025; % s n = 11; % harmonic number of carrier m = 2.55; % modulation index y = synth_QFM(n, m); 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');