数学建模社区-数学中国

标题: 重心法修订FFT代码及效果 [打印本页]

作者: 快到碗里来 时间: 2013-11-8 16:18
标题: 重心法修订FFT代码及效果
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))

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)