此文件夹中的 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 变换(非均匀-->非均匀),但可以通过结合此处
2021-09-29 11:27:23
45KB
matlab
1