我这有个初程序不知对你有没有用 *tFHM &a
clear k8yEdi`
dt=; %时间间隔 )6MfRw
s= %信号数组 m,28u3@r
len=128; ;1W6G=m
F=2*pi*(0:len/2)/len/dt; 'c&Ed
fs=fft(s,len); OdbEq?3S/?
p=fs.*conj(fs) "9uKtQS0o
p=p(1:len/2+1); mM~qBrwL
plot(F,p) 8,Z_{R#|
ylabel('power spectrum density ') ' {OgN}'{
xlabel('sequence ') M7\szv\Zc=
title('03051171') TprTWod2]t
参数自已改吧