frequency1.m

来自「speechtimefre 语音时域处理」· M 代码 · 共 14 行

load mtlb;
figure;
x=mtlb;
framelen=240;
frameinc=80;
specgram(x,240,8000,80);

y=enframe(x,framelen,frameinc);
z=fft(y);
Pyy = z.* conj(y) / 512;
f = 1000*(0:256)/512;
plot(f,Pyy(1:257))
title('Frequency content of y')
xlabel('frequency (Hz)')