重心法修订FFT代码及效果

[复制链接]

字体大小: 正常放大

快到碗里来

54 主题	9 听众	130 积分

TA的每日心情

	衰 2013-11-8 15:00

签到天数: 15 天

[LV.4]偶尔看看III

自我介绍: 我笑一笑从容

电梯直达

1^#

发表于 2013-11-8 16:18 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

t=0:0.1:5;
y=cos(2*pi*2*t)+sin(2*pi*2*t)*j
y=fftshift(fft(y,64));
f=linspace(-5,5,64);
figure(1)
subplot(2,1,1)
plot(f,abs(y))
%%重心法修订FFT主瓣重心%%
[m,n]=find(abs(y)==max(max(abs(y))));
z=abs(y(n+1))-abs(y(n-1));
if z>=0
n1=n+abs(y(n+1))/(abs(y(n))+abs(y(n+1)));
end
if z<0
n1=n+abs(y(n-1))*(-1)/(abs(y(n))+abs(y(n-1)));
end
df=-5+(n1-1)*10/63;
f(n)=df;%修订后主瓣重心
subplot(2,1,2)
plot(f,abs(y))

zan