请教三维谱阵图画法

spiderhui3300

本人是做故障诊断的，在用信号分析的三维谱阵时，发现图总出错，希望得到高手的指点，程序如下：
x=load(“文件”);
[m,n]=size(x);
a=x(:,1);
b=x(:,2);
m;
k=fix(m/1024);
fs=1000;
for i=1:k
l1=1024*(i-1);
if l1==0;
    l1=1;
end
l2=1024*i;
t=l1:1:l2;
Y=fft(b(t),512);
Pyy=Y.*conj(Y)/512;
f=fs*(0:256)/512;
  z(:,i)=Pyy(1:257);
  j=i;
waterfall(z(:,j))
end

三维谱阵的x轴为频率y轴为各频谱的数z轴为振动的幅值。
谢谢了