我这有个初程序不知对你有没有用 W5TqC
clear _Wq7U1v`
dt=; %时间间隔 }'%$7vL`Ft
s= %信号数组 4
3}qaf[
len=128; imC&pPBB/G
F=2*pi*(0:len/2)/len/dt; kVy"+ZebK
fs=fft(s,len); FW/6{tm
p=fs.*conj(fs) ><MGZ?-N
p=p(1:len/2+1); HW^{ ;'kH~
plot(F,p) Lqg7D\7j
ylabel('power spectrum density ') :p,c%"8
xlabel('sequence ') |.^^|@+
title('03051171') VOD1xWrb
参数自已改吧