我这有个初程序不知对你有没有用 I~c}&'V
clear 2{#=Ygb0
dt=; %时间间隔 Aj=c,]2
s= %信号数组 /M_kJe,%
len=128; C6b(\#g(
F=2*pi*(0:len/2)/len/dt; f P1fm
fs=fft(s,len); w~*"mZaG
p=fs.*conj(fs) 1)
G6
p=p(1:len/2+1); ,u>K##X\
plot(F,p) yAVt[+0
ylabel('power spectrum density ') o/\z4Ri)$
xlabel('sequence ') ^KFwO=I@PV
title('03051171') 7kidPAhY
参数自已改吧