重心法修订FFT代码及效果

Kinn1

t=0:0.1:5;
: [' r) w0 c9 p  Q9 M+ @7 @8 I5 Qy=cos(2*pi*2*t)+sin(2*pi*2*t)*j
y=fftshift(fft(y,64));
# F  {& I" q5 v1 m6 Hf=linspace(-5,5,64);
# F7 v# X/ f; T7 Vfigure(1)
. v- o- t" H4 x) v" K- Tsubplot(2,1,1)
plot(f,abs(y))
%%重心法修订FFT主瓣重心%%
[m,n]=find(abs(y)==max(max(abs(y))));
! u* k+ ~, M, J) j( f) Q, [z=abs(y(n+1))-abs(y(n-1));
$ \7 c  f; D9 t/ R$ m# Kif z>=0
9 O0 w0 x* W' ?    n1=n+abs(y(n+1))/(abs(y(n))+abs(y(n+1)));
! n* e* L; l( t8 k: V0 r6 Mend
if z<0
    n1=n+abs(y(n-1))*(-1)/(abs(y(n))+abs(y(n-1)));
# Z' Z6 T7 C3 F/ N- i3 @end
df=-5+(n1-1)*10/63;
f(n)=df;%修订后主瓣重心
6 [% l" ], y, f1 o! V' y% W! usubplot(2,1,2)
6 D2 g# B+ p; @9 @! Uplot(f,abs(y))
+ w$ r: z' A2 K  }
6 @6 ^& U5 ~+ V+ |1 _3 i1 x

freedom1

重心法修订FFT代码及效果