大家看看我得到功率谱密度的图的程序对不对？

lovelydeath

function [Y]= sheji(file)
data=wavread(file);%语音文件读取
dat(1,:)=data(:,1);%取单路数据
ave=mean(dat);%计算直流分量
dat=dat-ave;%减去直流分量
w=hamming(256);%hamming 窗函数
dat=conv(dat,w);%卷积滤波
n=length(dat);%求dat长度
X=fft(dat);%进行FFT变换
%subplot(2,1,1);
%plot(abs(X));%画图
Y(:,floor(1:n/5))=X(:,floor(1:n/5));%取能量较集中的低频段
%subplot(2,1,2);
Y=abs(Y);
Y=Y.^2;%Y矩阵元素平方
plot(Y);%画图

画出的就是功率谱密度 对不对？ 不对的话错在哪里