fftpsd.m

这是一个用于语音信号处理的工具箱

%fftpsd.m

function [sourpsd]=fftpsd(data,dsamp)
% Draw the FFT PSD of glottal volume velocity

	threshold=0.00001;

	NP = 1024;
	LP = NP/2;
	Hp = zeros(1,NP);
	for k=1:dsamp,
		Hp(k) = data(k);
	end
	Hp=realft(Hp, LP, 1);

	temp=Hp(1)*Hp(1);
	if (temp<threshold)
		temp=temp+threshold;
	end
	Hp(1)=10.0*log10(temp);

	for k=2:LP,
		temp=Hp(2*k-1)*Hp(2*k-1)+Hp(2*k)*Hp(2*k);
		if (temp<threshold)
			 temp=temp+threshold;
		end
		Hp(k)=10.0*log10(temp);
	end

	temp=Hp(2)*Hp(2);
	if (temp<threshold)
		temp=temp+threshold;
	end
	Hp(LP+1)=10.0*log10(temp);      % in dB */
	sourpsd=Hp(1:LP);