nufftn

nufftn 函数在对非等间距采样点或非等间距查询点执行运算时性能得到了改进。

例如，以下代码构造一个由非均匀采样点构成的 262144×3 矩阵 t，并在一个 64×64×64 数组的每个维度上计算非均匀离散傅里叶变换。代码的执行速度比上一版本大约快 3.3 倍。

function timingSamplePoints
rng default
t = rand(64^3,3);
X = rand(64,64,64);
tic
  Y = nufftn(X,t);
toc
end

大致的执行时间是：

R2023a：0.40 秒

R2023b：0.12 秒

代码是在运行 Windows^® 10 的 Intel^® Xeon^® CPU E5-1650 v4 @ 3.60 GHz 测试系统上通过调用 timingSamplePoints 函数进行计时的。

作为另一个示例，以下代码构造一个由非均匀查询点构成的 262144×3 矩阵 f，并在一个 64×64×64 数组的每个维度上计算非均匀离散傅里叶变换。代码执行速度比上一版本大约快 1.6 倍。

function timingQueryPoints
rng default
f = rand(64^3,3);
X = rand(64,64,64);
tic
  Y = nufftn(X,[],f);
toc
end

大致的执行时间是：

R2023a：0.40 秒

R2023b：0.25 秒

代码是在运行 Windows 10 的 Intel Xeon CPU E5-1650 v4 @ 3.60 GHz 测试系统上通过调用 timingQueryPoints 函数进行计时的。

nufftn 函数在对非等间距采样点或非等间距查询点执行运算时性能得到了改进。

例如，以下代码构造一个由非均匀采样点构成的 32768×3 矩阵 t，并在一个 32×32×32 数组的每个维度上计算非均匀离散傅里叶变换。其运行速度约为上一版本的 14.5 倍：

function timingSamplePoints
rng default
t = rand(32^3,3);
X = rand(32,32,32);
tic;
  Y = nufftn(X,t);
toc
end

大致的执行时间是：

R2021b：2.76 秒

R2022a：0.19 秒

代码是在运行 Windows 10 的 Intel Xeon CPU E5-1650 v4 @ 3.60 GHz 测试系统上通过调用 timingSamplePoints 函数进行计时的。

作为另一个示例，以下代码构造一个由非均匀查询点构成的 65536×3 矩阵 f，并在一个 64×32×32 数组的每个维度上计算非均匀离散傅里叶变换。其运行速度约为上一版本的 42.6 倍：

function timingQueryPoints
rng default
f = rand(64*32*32,3);
X = rand(64,32,32);
tic;
  Y = nufftn(X,[],f);
toc
end

大致的执行时间是：

R2021b：4.26 秒

R2022a：0.10 秒

代码是在运行 Windows 10 的 Intel Xeon CPU E5-1650 v4 @ 3.60 GHz 测试系统上通过调用 timingQueryPoints 函数进行计时的。