角谱计算时的fftshift及其原理

做一个fft运算,第一个事先用fftshift对待变换的E0进行操作,第二个没有用fftshift,第三个没有用fftshift但是进行了相位手动修正:

%%用fft进行角谱传输计算
%对比fft运算与傅里叶变换(黎曼和)的区别以及修正
clear
clc
close all

scale = 1.5e-4;
N = 50;
d0 = 20e-6;%sampling interval
wave_len=632.8e-9;%wavelength in nm 
k0=2*pi/wave_len;% vacuum wave number
fs=1/d0;%sample frequency
z=k0*scale^2;

x=d0*(-0.5*N:1:0.5*N-1); 
E0=abs(x)

角谱计算时的fftshift及其原理_第1张图片

相位校正原理:

角谱计算时的fftshift及其原理_第2张图片

可以参考文献:FFT计算菲涅尔衍射相位的跳变与矫正研究

你可能感兴趣的:(matlab,角谱)