NUFFT、NFFT、USFFT：通过“快速高斯网格化”实现 1D-3D NUFFT。-matlab开发

此文件夹中的 Matlab 代码通过快速高斯网格实现 1D-3D NUFFT。 卷积循环被编写为 C 程序，可从 Matlab 命令提示符编译为 mex 文件。 关于代码的更多数学细节可以在 L. Greengard 和 J. Lee 的“加速非均匀快速傅立叶变换”，SIAM 评论，卷。 46，第 3 期，第 443-454 页。 我已经包含了三个示例脚本（fgg_1D_experiment.m 等），它们创建了一个简单的 1D/2D/3D 图像并将 DFT 与 Type-I NUFFT（从非均匀数据到均匀图像网格的 DFT）进行了比较。 图像数据通过伴随算子转换回数据域（Type-II NUFFT——均匀网格 DFT 到非均匀数据位置——用 IFFT 而不是 FFT 实现）并再次返回图像域以证明数值精度. 此代码不包括 Type-III 变换（非均匀-->非均匀），但可以通过结合此处