我曾经也遇到过相同的问题,到网上四处找人帮助,现给你一个matlab编的小程序,不知有没有用 NG-`ag`s
clear D0M!"c>\
dt=0.017; wiV&xl
s=;%输入序列对数组或波的表达式 d=nh
len=256; sMJ#<w}Q
F=2*pi*(0:len/2)/len/dt; m^0A?jBrR
fs=fft(s,len); z\$;'
p=fs.*conj(fs) jeFN*r_
p=p(1:len/2+1); # ITLz!gE
plot(F,p) 03"#J2b
ylabel('power spectrum density ') WB;J1TpM7
xlabel('sequence ')