如何获得信号的相位 fft &-- 我可以在时域中获得相位吗？

Question

我可以得到来自 .wav 文件的信号幅度，但是如何也得到该信号的相位, 这是我浏览 .wav 文件并提取信号的地方

def browse_wav(self):

    filepath = QtGui.QFileDialog.getOpenFileName(self, 'Single File', "C:\Users\Hanna Nabil\Documents",'*.wav')
    f= str(filepath)
    if f != "":
        spf = wave.open(f, 'r')
    import contextlib

    with contextlib.closing(wave.open(f, 'r')) as f:
        frames = f.getnframes()
        rate = f.getframerate()
        duration = frames / float(rate)
        print "Duration is " , duration

    # Extract Raw Audio from Wav File
    self.signal = spf.readframes(-1)
    self.signal = np.fromstring(self.signal, 'Int16')
    self.fs = spf.getframerate()
    print "Sampling Rate is " ,self.fs

    # If Stereo
    if spf.getnchannels() == 2:
        print 'Just mono files'
        sys.exit(0)

    #self.time = np.linspace(0, len(self.signal) / fs, num=len(self.signal))
    self.time = np.linspace(0, duration, self.fs * duration)

    self.xfourier = fftfreq(self.signal.size, d=self.time[1] - self.time[0])
    self.yfourier = np.abs(fft(self.signal))  # signal magnitude

    self.zico = self.yfourier
    self.cut_signal = ifft(self.zico)

Answer 1

复谱包含幅度和相位。您可以通过计算绝对值获得幅度，通过计算角度获得相位。您可以使用 numpy.angle() 获取相位：

spectrum = fft(self.signal)
magnitude = np.abs(spectrum)
phase = np.angle(spectrum)

Answer 2

对于 Lukas Lalinsky 的回答，我发现出于某种原因我必须添加 pi/2 才能获得正确答案。

还有积分：

N = len(t)
Amplitude = 2/N *  np.abs(magnitude)

如何获得信号的相位 fft &-- 我可以在时域中获得相位吗？

how to get phase fft of a signal &-- can i get phase in time domain?

python

signals

signal-processing

fft

python-2.7