f_fft.m

来自「audio compression matlab code」· M 代码 · 共 21 行

function [f, df, spec] = f_fft(fs, s);
%INPUT:
% - semnalului esantionat, s;
% - frecventa de esantionare, fs;
%OUTPUT: 
% - axa frecventei, f;
% - ecartul in frecventa, df;
% - densitatea spectrala de putere, spec;
%-------------------------------

n = length(s);
s = s.* hann(n);                     % windowing ..;
df=fs/n;						    % ecartul in frecventa;
f=0:df:(n/2-1)*df;					% axa frecventei;
ys=fft(s,n);						% densitatea spectrala de energie;
abs_ys=ys.*conj(ys)./n;				% densitatea spectrala de amplitudine;
nf=length(f);						% numarul de puncte din domeniul frecventa;
spec=abs(abs_ys(1:nf));			    % spectrul de putere unilateral;	
%----------------------------------------------------------------------|