我曾经也遇到过相同的问题,到网上四处找人帮助,现给你一个matlab编的小程序,不知有没有用
Xi~I<&
clear zFY$^Oz"_
dt=0.017; :c(I-xif
s=;%输入序列对数组或波的表达式 LaL{
^wP
len=256; =tX"aCW~
F=2*pi*(0:len/2)/len/dt; aU#r`D@0
fs=fft(s,len);
8hMy$
p=fs.*conj(fs) @Y ?p-&
p=p(1:len/2+1); &KV$x3
plot(F,p) J:I As:e`
ylabel('power spectrum density ') e5w0}/yW/
xlabel('sequence ')