我这有个初程序不知对你有没有用 Tv,.
clear ^@lg5d3F
dt=; %时间间隔 m:fouMS
s= %信号数组 [j]J_S9jJ
len=128; wU)5Evp[
F=2*pi*(0:len/2)/len/dt; LiD |4(3
fs=fft(s,len); &Z;8J @
p=fs.*conj(fs) [Pwo,L,)
p=p(1:len/2+1); [2 w<F[
plot(F,p) :#:O(K1PW
ylabel('power spectrum density ') pUMB)(<k
xlabel('sequence ') w+q;dc8
title('03051171') 9'#.>Q>0=j
参数自已改吧