重组正弦波的相位和幅度？

Question

我一直在研究音频信号的去相关。我正在关注一篇研究论文，其中的技术是以幅度恒定而相位随机变化的方式生成序列号。
现在，问题是我生成了一个正弦波，我对它进行了 fft 计算，使用 NumPy 计算了相位和幅度，如下面的代码所示。

import numpy as np
from scipy.fftpack import fft, rfft, ifft, irfft

Farendsig = ((np.sin(2 * np.pi * np.arange(48000 * 0.02133) * 2000 / 48000)) * 32767).astype(np.int16)

freq = fft(Farendsig, 1024)  # fft output complex

phase = np.angle(freq)  # phase

magnitude = np.abs(freq)  # magnitude

combined = np.multiply(magnitude, np.exp(1j * phase))

FarendSyth = (ifft(combined, 1024) * (32767 + 0j)).astype(np.int16)

现在，当我按照 Whosebug 的回答重新合成相同的信号时：如何使用 python 将一幅图像的相位和另一幅图像的幅度组合成一幅图像。结果信号在某些指标上与原始信号不同。
我在这里做什么错了？

谢谢。

Answer 1

没关系，我找到了答案：

FarendSyth = (ifft(combined, 1024)).astype(np.int1

我将 ifft 结果与 32767 相乘，这导致了结果中的问题。替换上面代码中的这一行给出了准确的结果。

谢谢

重组正弦波的相位和幅度？

Recombining phase and magnitude of sin wave?

signal-processing

numpy

fft

scipy