我这有个初程序不知对你有没有用 #A 7|=E
clear Vs9]Gm
dt=; %时间间隔 rAZ~R PrW
s= %信号数组 PB;j4
len=128; 'Gqo{wl
F=2*pi*(0:len/2)/len/dt; >)edha*W]
fs=fft(s,len); gv7@4G
p=fs.*conj(fs) C#@-uo2
p=p(1:len/2+1); srryVqgS
plot(F,p) ~BC~^D&WD
ylabel('power spectrum density ') M~g{}_0Z
xlabel('sequence ') jP\5bg-}
title('03051171') nk"nSXm3SR
参数自已改吧