Han*_*nna 3 python signals signal-processing fft python-2.7
我可以从 .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)
Run Code Online (Sandbox Code Playgroud)
复频谱包含幅度和相位。您可以通过计算绝对值来获得幅度,并通过计算角度来获得相位。您可以使用numpy.angle()来获取阶段:
spectrum = fft(self.signal)
magnitude = np.abs(spectrum)
phase = np.angle(spectrum)
Run Code Online (Sandbox Code Playgroud)
| 归档时间: |
|
| 查看次数: |
12695 次 |
| 最近记录: |