Scipy的快速傅里叶变换

渣渣小码

已于 2023-07-21 16:14:28 修改

阅读量496

点赞数

分类专栏：数字信号处理和应用文章标签： scipy python 开发语言

于 2023-05-24 17:24:44 首次发布

本文链接：https://blog.csdn.net/qq_40088639/article/details/130851632

版权

数字信号处理和应用专栏收录该内容

28 篇文章 3 订阅 ¥19.90 ¥99.00

订阅专栏

超级会员免费看

该代码示例展示了如何使用Python的librosa、numpy和scipy库对音频文件进行傅立叶变换(FFT)和逆傅立叶变换(IFFT)。首先，读取音频文件，然后计算FFT得到复数序列，接着利用相位和模构建新的复数序列，最后通过IFFT将结果转换回时域并保存为新音频文件。过程中，还绘制了波形图以可视化结果。

摘要由CSDN通过智能技术生成

1. 编程实现

读取音频文件，fft后在频域上使用模和相角的方式构造复数序列，再逆变换回时域，对比两个音频文件。

import numpy as np
import librosa
from scipy.fftpack import fft,ifft
import matplotlib.pyplot as plt
import soundfile as sf
from matplotlib import rcParams

config = {
    "font.family": 'serif',     # 衬线字体
    "font.size": 10,            # 相当于小四大小
    "font.serif": ['SimSun'],   # 宋体
    "mathtext.fontset": 'stix', # matplotlib渲染数学字体时使用的字体，和Times New Roman差别不大
    'axes.unicode_minus': False # 处理负号，即-号
}
#配置当前终端
rcParams.update(config)

audio_data = 'num1.wav'
x , fs = librosa.load(audio_data, sr=None)
N = x.shape[0]
xx = np.linspace(0, 1, N)  # 将0到1平分成N份

# FFT 
fft_y = fft(x)

# 求相角
phase_x= np.angle(fft_y)

# print(len(phase_x))
# 取复数的绝对值，即复数的模
mag_x = np.abs(fft_y)

# 使用模和相角构建复数序列
new_spec = mag_x*np.exp(1j*phase_x)
# IFFT