重心法修订FFT代码及效果

Kinn1

t=0:0.1:5;
# A- G& m- k# Q$ ly=cos(2*pi*2*t)+sin(2*pi*2*t)*j
- q: n1 m& s3 Z; P( ^) z! ny=fftshift(fft(y,64));
f=linspace(-5,5,64);
6 V& Q2 h, ~5 a, Y/ _4 h7 w: {figure(1)
# M  v1 F- w4 e1 `) Q" E3 Ksubplot(2,1,1)
8 A- k& r% I5 o/ w+ yplot(f,abs(y))
+ Y' ]4 B+ E, l/ k7 O%%重心法修订FFT主瓣重心%%
[m,n]=find(abs(y)==max(max(abs(y))));
- B/ p$ {$ F5 g& t- x0 d* Tz=abs(y(n+1))-abs(y(n-1));
4 o5 Z& |$ W! Q$ n2 m4 sif z>=0
    n1=n+abs(y(n+1))/(abs(y(n))+abs(y(n+1)));
8 R4 H# x/ X2 e) O' m2 iend
! A9 a* j& |( E" bif z<0
/ k% z$ \& |# A9 [! k" z( }    n1=n+abs(y(n-1))*(-1)/(abs(y(n))+abs(y(n-1)));
: @: T& M4 D% }* }$ Lend
% J3 l7 e* q: K: P. q" udf=-5+(n1-1)*10/63;
' @- U9 n' E* b' N) w& x7 {9 Lf(n)=df;%修订后主瓣重心
subplot(2,1,2)
1 `2 K. }! Q0 ^5 m; Qplot(f,abs(y))

freedom1

重心法修订FFT代码及效果