我这有个初程序不知对你有没有用 Cj31'
clear pEw &i
dt=; %时间间隔 9Psy$
s= %信号数组 Yhb=^)@))
len=128; 3;>|*(cO
F=2*pi*(0:len/2)/len/dt; l{E+j%
fs=fft(s,len); BJ @tUn
p=fs.*conj(fs) *U)!9DvA
p=p(1:len/2+1); ~NMal]Fwx
plot(F,p) qWfG@hn
ylabel('power spectrum density ') ?sdVd
xlabel('sequence ') &w- QMjM>
title('03051171') x{6/di
参数自已改吧