我这有个初程序不知对你有没有用 OxD\e5r
clear pZH
bj2~
dt=; %时间间隔 -du+iOe?
s= %信号数组 eSvu:euv
len=128; -F*j`
F=2*pi*(0:len/2)/len/dt; 'n=FBu^
fs=fft(s,len); }x%"Oq|2]x
p=fs.*conj(fs) ALKhZFuz
p=p(1:len/2+1); /O8'8 sL5
plot(F,p) Qc3!FW<26
ylabel('power spectrum density ') 9)=as/o
xlabel('sequence ') ]$z~;\ T
title('03051171') sG\=_-"v(
参数自已改吧