有关音频信号的傅立叶变换问题

wdwyhy

我的毕业设计里有个关于已有的音频信号进行傅立叶变换的问题，首先是[y,fs,bits]=wavread('lntu.wav')，如果用dft函数：function [Xk]=dft(xn,N)
%计算离散傅立叶变换
%---------------------------
%[XK]=dft(xn,N)
%XK=在0<=n<=N-1间的DFT系数数组
%xn=N点有限持续时间序列
%N=DFT的长度
%
n=[0:1:N-1];   %n的行向量
k=[0:1:N-1];   %k的行向量
WN=exp(-j*2*pi/N); %Wn因子
nk=n'*k;       %产生一个含nk值的N乘N维矩阵
WNnk=WN.^nk;   %DFT矩阵
Xk=xn*WNnk;    %DFT系数的行向量
将它存为dft,m，如何用这个函数来进行语音信号的傅立叶变换，并画出变换后的频谱图，还有怎么查看y中有多少个数据啊？
要是用fft的画，如何画出啊？
我很着急，希望兄弟姐妹门帮帮我！谢谢了！

toniling310

用 size(y) 查看 我也是新手 互相帮助下~~~