发布网友 发布时间:2022-04-22 15:21
共1个回答
热心网友 时间:2023-08-17 02:06
通用的回答Matlab程序实现片断:%%Stolt插值P=8/2;%8点sinc插值Stolt_kxky=zeros(Nx,Ny);%初始化dkx=1/dx/Nx;%频域间隔KX_kxky=(sqrt((kx+kxc).^2'*ones(1,Ny)+ones(Nx,1)*ky.^2)-kxc)-kx'*ones(1,Ny);KX=fix(KX_kxky/dkx);for i=P:Nx-max(max(KX))-Pfor j=1:NyT=KX(i,j)+(-P+1 );f=sxcp_kxky(i+T,j).';Stolt_kxky(i,j)=sum(f.*sinc(KX_kxky(i,j)/dkx-T));endiend