我这有个初程序不知对你有没有用 idLWe9gC
clear jyZ (RB
dt=; %时间间隔 = bfJ^]R
s= %信号数组 ^uKwB;@
len=128; uaT!(Y6
F=2*pi*(0:len/2)/len/dt; Ro\8ZXUQa
fs=fft(s,len); BF\XEm?!
p=fs.*conj(fs) q
?|,O;?
p=p(1:len/2+1); ,]|#[ 8
plot(F,p) 94h]~GqNi
ylabel('power spectrum density ') Lt u'W22
xlabel('sequence ') #eoome2Q
title('03051171') kjIAep0rT
参数自已改吧