NUFFT算法及说明

资源结合李银斌师兄的论文第三章，编写得到NUFFT算法，实现从非均匀采样点得到均匀频谱，自己使用正确有效。希望能帮到需要的同学

function [zk] = nufft1_z（zlxlul）
i = sqrt（-1）;
c = 2;
K = 3;
alpha = pi*（2-1/c）-0.01;
N = length（zl）;
M = N;
zk = zeros（1N）;
for k=-N/2:N/2-1
    tmp = 0;
    for l=1:M
        for m=-K:K
            x = c*xl（l）-ul（l）+m;
            if（（K^2-x^2）==0）
                x = x+0.01;
            end
            tmp = tmp + zl（l）*1/pi*sinh（alpha*sqrt（K^2-x^2））/sqrt（K^2-x^2）*exp（-2*pi*i*k*（ul（l）+m）/c/N）;
        end
    end
    y = 2*pi*k/c/N;
    zk（k+N/2+1） = tmp/besseli（0K*sqrt（alpha^2-y^2））;
end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        525  2018-03-01 10:37  NUFFT\nufft1_z.m

     文件        127  2018-03-28 16:40  NUFFT\说明.txt

     文件    2638636  2017-11-22 10:11  NUFFT\非均匀快速傅里叶变换及其在SAR成像中的应用_李银斌.caj

     目录          0  2018-03-28 16:40  NUFFT

----------- ---------  ---------- -----  ----

              2639288                    4