图像的zoomfft变换

 

 算法实现基础原理如下：

 

部分核心代码如下：

%这里开始正式处理
% 步骤一：乘以exp
zoom_fft_xx   = (x_real_zoom+j*x_imag_zoom).*exp(-j*2*pi*(0:N-1)*frequency_shift/Fs);
% 步骤二：数字低通+重采样
zoom_fft_xx   =  zoom_fft_xx.*w;
zoom_fft_xx   =  zoom_fft_xx-sum(zoom_fft_xx)/N;
% 步骤三：FFT变化
zoom_fft_xx   =  fft(zoom_fft_xx);
% 步骤四：频率调整
zoomfft_x     =  zoomfft_x+(abs(zoom_fft_xx(1:N/2))/N*2).^2;
zoomfft_x     =  sqrt(zoomfft_x);
%定位到第二个图上
axes(handles.axes2);
plot((0:N/2-1)*Fs/N,fft_x);hold on;%显示原始的频谱图
plot(f,zoomfft_x,'r');hold off;    %显示zoomfft处理后的频谱图
axis([f1 f2 0 1.4*max(fft_x)]);    %显示坐标范围
grid on;                           %显示方格

A23-02