我这有个初程序不知对你有没有用 w#|L8VAh
clear o[Jzx2A<
dt=; %时间间隔 }|kFHodo
s= %信号数组 &Vz$0{d5
len=128; TJO$r6&
F=2*pi*(0:len/2)/len/dt; \g39>;iR
fs=fft(s,len); FQDf?d5
p=fs.*conj(fs) _{K mj,q
p=p(1:len/2+1); ,_Z(!|
rW
plot(F,p) kt/,& oKI
ylabel('power spectrum density ') 20;9XJmjl
xlabel('sequence ') @Nek;xJ
title('03051171') N]A# ecm
参数自已改吧