spec_rolloff.m

非常好的数字处理教程

function sc_value = spec_rolloff(signal,Fs)

	FFT = fft(signal);
	mag = abs(FFT);
	mag = mag(1:length(mag)/2);
	spsum = 0;

	% Sum entire spectrum
	for i=1:length(mag),
		spsum = spsum + mag(i);
	end

	% Find bin where 85% of spectral energy is below
	cut = 0.15*spsum;
	i = 1;
	while (spsum > cut)
		spsum = spsum - mag(i);
		i = i+1;
	end

	Hz_per_bin = Fs/length(signal);
	sc_value = Hz_per_bin*i;