我曾经也遇到过相同的问题,到网上四处找人帮助,现给你一个matlab编的小程序,不知有没有用 <G+IbUG:
clear {:peArO
dt=0.017; ]Ak/:pu
s=;%输入序列对数组或波的表达式 $G
$147z
len=256; o3=2`BvJ
F=2*pi*(0:len/2)/len/dt; %yr(i 6L
fs=fft(s,len); 1MVzu7
p=fs.*conj(fs) @nPXu2c?u7
p=p(1:len/2+1); y5oC|v7
plot(F,p) qKL:#ny
ylabel('power spectrum density ') Tj[=E
xlabel('sequence ')