gusucode.com > 声音的处理有：LPC，FFT，共振峰，频谱源码程序 > siganlandsystemusingMatlab/SSUM/library/pitch.m

    function pitch_val = pitch(signal,Fs)

	% http://www.edgar-berdahl.com/Projects/Grundfrequenzanalyse/
	vals = xcorr(signal,'coef');
	vals(1:floor(length(vals)/2)) = [];
	vals

	plot(vals)