我这有个初程序不知对你有没有用 7+I2"Hy
clear J12ZdC'O
dt=; %时间间隔 @pD']=d}t
s= %信号数组 afMIq Q?
len=128; zSiSZMP"
F=2*pi*(0:len/2)/len/dt; mTzzF9n"Y
fs=fft(s,len); g{cHh(S
p=fs.*conj(fs) '?~k`zK
p=p(1:len/2+1); FT*
o;&_QS
plot(F,p) $v6dB {%Qu
ylabel('power spectrum density ') df7z&{R
xlabel('sequence ') 5)yQrS !{:
title('03051171') TvWhy`RQ
参数自已改吧