求各位大佬解答时间序列的相位随机化(如下)如何实现？可以使用R语言实现吗？

对时间序列进行相位随机化处理的步骤如下：
(1)  对原始序列进行离散 Fourier 变换；
(2)  随机旋转任意相位角；
(3)  进行 Fourier 逆变换。