重心法修订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主瓣重心%%
=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))
页:
[1]