频域消除音频中高频类似微电流唦唦声（见代码）

云开-秋秋

于 2023-06-27 10:43:50 发布

阅读量178

点赞数

CC 4.0 BY-SA版权

文章标签： python 降噪语音增强电流沙沙声傅里叶变换频域

本文链接：https://blog.csdn.net/qiufeng_xinqing/article/details/131412139

该文通过使用librosa库进行STFT（短时傅里叶变换），结合自定义的加窗函数，针对hello.wav文件中的高频类微电流声进行消除。通过调整窗口数据并设置特定频率范围的系数为0，然后应用ISTFT（逆短时傅里叶变换）恢复信号，生成了处理后的输出文件。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

### 傅里叶加窗消除高频类微电流唦唦声
import matplotlib.pyplot as plt
from scipy import signal
from scipy.io import wavfile
import librosa
import numpy as np
import soundfile as sf

y, sr = librosa.load('hello.wav')
#S = librosa.stft(y, n_fft=2048, window="hann")
#S = librosa.stft(y, n_fft=2048, window=filterXY[:,1]/400) # 使用折线 y 作为窗，得到颤音
#S = librosa.stft(y, n_fft=2048, window=filterXY[:,0]/400) # 使用折线 x 作为窗，抹出高频音频
windowdata = []
sep = 1.05/2048
for i in range(2048):
    #windowdata.append(math.fabs(1024-i)*sep)
    windowdata.append((2048-i)*sep)

S = librosa.stft(y, n_fft=2048, window=np.asarray(windowdata))
#print(type(S), len(S), S, S.shape)

for i in range(350,1025):
    for j in range(len(S[i])):
        S[i][j] = 0.0

for i in range(350):
    S[i] = S[i] * 1.33

# 加高频静音不会消除微沙电流音，添加后音频会变调
#zeros = np.zeros(shape=(1,len(S[0])))
#for i in range(200):
#    S = np.append(S, zeros, axis=0)

print(S, S.shape)
D = S
y_inv = librosa.istft(D)
#y_inv = y_inv * 1.33
sf.write('hello-output_file.wav', y_inv, sr)