powspecdf.m

计算声质中的Aures Tonality

function [Yxx,f]=PowSpecdF(y,fs,df)

NFFT=round(ceil(fs/df));

NOVERLAP=round(0.75*NFFT);
[Yxx,f]=psd(y,NFFT,fs,NFFT,NOVERLAP);

Yxx=2*Yxx/NFFT;
Yxx=Yxx';
return