主要内容

nufft

非均匀快速傅里叶变换

描述

例子

Y= nufft(Xt返回非均匀离散傅里叶变换X使用样本点t

  • 如果X是矢量吗nufft返回向量的变换。

  • 如果X是一个矩阵,然后nufft处理X作为向量并返回每列的变换。

  • 如果X那是一个多维数组nufft沿第一阵列维度对其大小不等于1作为向量并返回每个向量的变换。

例子

Y= nufft(Xtf在查询点计算NUDFTf使用样本点t.指定f不指定样本点,使用nufft (X, [], f)

Y= nufft(Xtf昏暗的沿着维度返回NUDFT昏暗的.例如,nufft (X, t, f, 2)计算每行矩阵的变换X

Y= nufft(X的离散傅里叶变换X,等于fft (X)

例子

全部折叠

创建一个信号X在不均匀间隔的点采样t.计算非均匀快速傅里叶变换Y

T = [0:300 500.5:700.5];s = 2 * sin(0.1 * pi * t)+ sin(0.02 * pi * t);x = s + rand(尺寸(t));y = nufft(x,t);

绘制变换的绝对值作为默认频率的函数。

n =长度(t);f = (0: n - 1) / n;情节(f, abs (Y))

图中包含一个轴。轴包含一个类型为line的对象。

定义和标记一系列音乐音调的频率。

C3 = 440 / (2^(21/12));nOctaves = 3;= C3 * 2.^((0:(12* noctaves1))/12);toneNames = [“C”“c#”“D”“D”#“E”“F”“F#”“G”“G#”“一种”“#”“b”) +字符串(3:(3 + nOctaves-1));toneNames = categorical(toneNames, toneNames);

定义音频信号的采样频率,以Hz为采样点n,以及包含主要和弦的信号X

fs = 16 e3;n = 1:16000;X = 4 * cos(2 *π* (440 / fs) * n) + 2 * cos(2 *π* (554.37 / fs) * n) + 3 * cos(2 *π* (659.2 / fs) * n);

计算和绘制主要和弦的频率分量。

Y = nufft (X, [], musicalTones / fs);栏(toneNames(:)、abs (Y))

图中包含一个轴。轴包含一个类型为bar的对象。

输入参数

全部折叠

输入数组,指定为矢量、矩阵或多维数组。

数据类型:||int8|int16|int32|uint8|uint16|uint32|逻辑
复数的支持:金宝app是的

采样点,指定为长度向量n,在那里n是输入数组的操作维数的长度吗X.默认情况下,样本点向量为0: (n - 1)

数据类型:|

查询点,指定为向量。默认情况下,查询点向量为(0: (n - 1) / n,在那里n是输入数组的操作维数的长度吗X.指定f不指定样本点,使用nufft (X, [], f)

数据类型:|

要操作的维数,指定为正整数标量。默认值是大小不等于1的第一个数组维数。

数据类型:||int8|int16|int32|int64|uint8|uint16|uint32|uint64|逻辑

更多关于

全部折叠

向量的非一致离散傅里叶变换

为一个向量X的长度n,采样点t,频率f的非一致离散傅里叶变换X被定义为

Y k j 1 n X j e 2 π t j f k

在哪里k= 1、2、….当t= 0, 1,…n1,f=(0,1,…n-1)/n(默认值为nufft的一致离散傅里叶变换FFT.函数。

参考

[1]波特,塞缪尔F., Nail A. Gumerov,和Ramani Duraiswami。带限函数的快速插值2017年IEEE声学、语音和信号处理(ICASSP)国际会议,4516-20。新奥尔良,路易斯安那州:IEEE, 2017。https://doi.org/10.1109/ICASSP.2017.7953011。

Dutt, A.和V. Rokhlin。非等距数据的快速傅里叶变换SIAM科学计算杂志14日,没有。6(1993年11月):1368-93。https://doi.org/10.1137/0914081。

另请参阅

|

介绍了R2020a