跪求哥哥姐姐 matlab中声音文件的时频域仿真图 时域图出来了 但频域图出不来 小妹跪求好心人帮忙纠正程序！

Question

%%读取原始音频文件 仿真时域波形
block_len=1024;
num_sb=16;
sb_switch=2.5;
p_sitch=20;
no_lp=0;
need_lp=1;
[X,fs,bits]=wavread('F:\2012biyesheji\biyesheji5-9\55.wav');
t=0:length(X)-1;
N=length(X);
figure(1);
plot(X);
title('时域波形');
hold on;
%%仿真频域波形（以下不知道哪里出错）
figure(2);
FFTSIZE=8000;
Y=zeros(FFTSIZE,1);
Y=20*log10(abs(fft(Y,FFTSIZE)));
plot(Y(1:4000));
title('频域波形');
hold on;
我觉得很可能是1 对数函数输入有问题 2如果好心大侠不能纠错，帮我另写一个频域仿真程序，小妹定是感激不尽！为表诚意，小妹在此把分数全部酬谢！

Answer 1

把Y=20*log10(abs(fft(Y,FFTSIZE))); 改为
Y=20*log10(abs(fft(X,FFTSIZE)));
就好了。你求了Y的傅里叶变化，因Y的赋值为0，你得到的对数值是 -Inf.