如何在 fftw 中使用离散余弦和离散正弦变换

计算科学 傅立叶分析 fftw
2021-12-24 17:26:38

我从事与流体相关的模拟。我使用 FFT 进行流体模拟。我想使用离散余弦变换 (DCT) 和离散正弦变换 (DST) 将我的速度场转换为波数。

我分别使用REDFT00andRODFT00作为 DCT 和 DST (所以反之亦然)。

我知道标准化是2(n1)2(n+1), 分别。

对于 DCT 和 DST,我采取

u(x)=cos(3x);n=64

但是当我反转它时,之后的值(n/2)+1都是不正确的。

正确的做法是什么?

更新:我使用的是一个经过大量修改的旧库(尽管 FFT 工作正常)。我使用了原始库,一切正常。但是现在我很难将库用于我的目的(不知道如何处理波数),我已将其作为一个单独的问题发布。所以我想我不需要答案。这是一个有缺陷的图书馆。

0个回答
没有发现任何回复~