CN114255780A - 一种基于深度神经网络的噪声鲁棒盲混响时间估计方法 - Google Patents

一种基于深度神经网络的噪声鲁棒盲混响时间估计方法 Download PDF

Info

Publication number
CN114255780A
CN114255780A CN202210004710.6A CN202210004710A CN114255780A CN 114255780 A CN114255780 A CN 114255780A CN 202210004710 A CN202210004710 A CN 202210004710A CN 114255780 A CN114255780 A CN 114255780A
Authority
CN
China
Prior art keywords
reverberation time
noise
model
estimation
network
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210004710.6A
Other languages
English (en)
Other versions
CN114255780B (zh
Inventor
郑成诗
郑凯桐
桑晋秋
张玉龙
张芳杰
李晓东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Acoustics CAS
Original Assignee
Institute of Acoustics CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Acoustics CAS filed Critical Institute of Acoustics CAS
Priority to CN202210004710.6A priority Critical patent/CN114255780B/zh
Publication of CN114255780A publication Critical patent/CN114255780A/zh
Application granted granted Critical
Publication of CN114255780B publication Critical patent/CN114255780B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208—Noise filtering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/04—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using predictive techniques
    • G10L19/16—Vocoder architecture
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208—Noise filtering
    • G10L21/0216—Noise filtering characterised by the method used for estimating noise
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208—Noise filtering
    • G10L21/0216—Noise filtering characterised by the method used for estimating noise
    • G10L21/0232—Processing in the frequency domain
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/27—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the analysis technique
    • G10L25/30—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the analysis technique using neural networks
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208—Noise filtering
    • G10L2021/02082—Noise filtering the noise being echo, reverberation of the speech
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Computational Linguistics (AREA)
  • Multimedia (AREA)
  • Acoustics & Sound (AREA)
  • Human Computer Interaction (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Signal Processing (AREA)
  • Quality & Reliability (AREA)
  • Artificial Intelligence (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Biophysics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Biomedical Technology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Reverberation, Karaoke And Other Acoustics (AREA)

Abstract

本发明公开了一种基于深度神经网络的噪声鲁棒盲混响时间估计方法，所述方法包括：对采集的语音信号进行预处理；将预处理后的语音信号输入预先建立和训练好的混响时间估计模型，得到混响时间值；所述混响时间估计模型基于深度神经网络和时频掩膜实现。本发明的方法利用混响语音进行房间混响时间估计，省去了传统声学测量所需要的专业的仪器与测量人员，相比于传统声学测量方法更加方便快捷；通过收集大量混响语音数据集，采用数据驱动的深度神经网络，可以显著提高房间混响时间估计的准确性。

Description

一种基于深度神经网络的噪声鲁棒盲混响时间估计方法
技术领域
本发明属于噪声鲁棒盲混响时间估计方法技术领域，尤其涉及一种基于深度神经网络的噪声鲁棒盲混响时间估计方法。
背景技术
当声音从封闭空间中的声源处传播到观测点时，在房间将产生多次反射声。这一物理过程产生了混响，即在关闭声源后，室内产生的声音在一段时间内仍可被听到。声音停止发声时减少60分贝所花费的时间被定义为混响时间(T60)，它主要取决于房间的几何形状和表面的反射率。ISO标准中规定的混响时间测量方法有中断声源法以及脉冲响应反向积分法。中断声源法需要使用十二面体球形声源发出全频带噪声，使得标准传声器处的声压级大于本底噪声至少45dB以上，关闭声源，记录声能量曲线(EDC)，通过EDC算得混响时间，中断声源法需要在房间中取多点进行平均，在每点都需要多次测量取平均。脉冲响应反向积分法需要对房间冲激响应(RIR)进行反向积分得到声能量衰减曲线。房间冲激响应可以通过在声源处输入扫频等宽带信号作为激励信号，在传声器处接收经过房间传播的激励信号，再利用激励信号与传声器处接收的信号进行解卷积求出。传统的测量方法需要专业仪器、专业人员进行测量，不利于大规模测量混响时间。因此，需要提出一种方便快捷的盲混响时间估计方法。
盲混响时间估计方法指的是直接从混响语音录音中估计T60。这就意味着在没有任何关于声环境或说话人先验信息的情况下，输入混响语音录音，估计算法直接输出估计的T60。几十年来，人们提出了许多盲T60估计方法，包括基于传统信号处理和基于机器学习的方法。
在真实的声学环境中，语音不仅包含混响，还包含噪声。噪声和晚期混响具有类似的性质，在存在噪声的情况下，混响声被淹没在噪声中，从而降低T60估计器的估计性能。因此，提出一种噪声鲁棒的盲混响时间估计方法是十分必要的。
以往的相关专利或论文均未采用过此种方法来搭建噪声鲁棒盲混响时间估计模型。
发明内容
本发明的目的在于克服现有技术缺陷，提出了一种基于深度神经网络的噪声鲁棒盲混响时间估计方法。
为了实现上述目的，本发明提出了一种基于深度神经网络的噪声鲁棒盲混响时间估计方法，所述方法包括：
对采集的语音信号进行预处理；
将预处理后的语音信号输入预先建立和训练好的混响时间估计模型，得到混响时间值；
所述混响时间估计模型基于深度神经网络和时频掩膜实现。
作为上述方法的一种改进，所述对采集的语音信号进行预处理；具体包括：
使用设定毫秒的汉宁窗口，根据设定的帧移和傅里叶变换点数，对采集的语音信号进行短时傅里叶变换，得到多维的频率特征。
作为上述方法的一种改进，所述混响时间估计模型为第一模型或第二模型；其中，
所述第一模型包括串联的噪声估计网络和混响时间估计网络；预处理后的语音信号同时输入噪声估计网络和混响时间估计网络，输出为语音信号的混响时间值；
所述第二模型包括串联的去噪网络和混响时间估计网络；预处理后的语音信号输入去噪网络，输出为语音信号的混响时间值。
作为上述方法的一种改进，所述第一模型的噪声估计网络包括编码器、MG-TCNs模块和解码器，其中，
所述编码器，用于对预处理后的语音信号提取隐变量特征；
所述MG-TCNs模块，用于对编码器提取的隐变量特征的时间序列进行时序建模，输出隐变量；
所述解码器包括并联的语音解码器和噪声解码器，分别用于将MG-TCNs模块输出的隐变量还原为与语音信号幅度谱维度一致的无噪声语音估计和噪声估计；
所述第一模型的混响时间估计网络，用于结合无噪声语音估计和噪声估计对预处理后的语音信号估算出混响时间值，所述混响时间估计网络包括编码器和MG-TCNs模块，均与噪声估计网络对应模块的结构相同。
作为上述方法的一种改进，所述第二模型的去噪网络包括编码器、MG-TCNs模块和语音解码器，其中，
所述编码器，用于对预处理后的语音信号提取隐变量特征；
所述MG-TCNs，用于对编码器提取的隐变量特征的时间序列进行时序建模，输出隐变量；
所述语音解码器用于将隐变量还原为与语音信号幅度谱维度一致的无噪声语音估计；
所述第二模型的混响时间估计网络，用于根据无噪声语音估计估算出混响时间值，所述混响时间估计网络包括编码器和MG-TCNs模块，均与去噪网络对应模块的结构相同。
作为上述方法的一种改进，所述方法还包括混响时间估计模型的训练步骤；具体包括：
准备纯净语音数据集、噪声数据集和房间冲激响应数据集，构成样本集；
从样本集中获取纯净语音和房间冲激响应，进行卷积并加上噪声得到训练语音构成训练集；
使用设定毫秒的汉宁窗口，根据设定的帧移和傅里叶变换点数，对训练集中的训练语音依次进行短时傅里叶变换，分别得到对应的多维的频率特征；
使用Adam优化器，设置两个训练阶段，分别设置每个训练阶段的学习率和训练批次及周期，将预处理后的训练集数据分别输入第一模型或第二模型，直至满足训练要求，得到训练好的第一模型或第二模型，从而完成混响时间估计模型的训练。
作为上述方法的一种改进，所述第一模型的第一训练阶段的损失函数LNE满足下式：
Figure BDA0003455104330000031
其中，α表示加权因子，取为0.5，
Figure BDA0003455104330000032
表示估计的无噪语音幅度谱，|X|表示无噪语音幅度谱，
Figure BDA0003455104330000033
表示估计的噪声幅度谱，|N|表示噪声幅度谱，
Figure BDA0003455104330000034
表示二范数；
所述第一模型的第二训练阶段的损失函数L满足下式：
L＝LRE+λLNE
其中，LRE为估计的混响时间与真实混响时间的均方误差，满足下式：
Figure BDA0003455104330000035
其中，λ表示微调加权因子，取为0.1，l为帧索引，T表示总帧数，
Figure BDA0003455104330000036
为第l帧的混响时间估计值，T60,l为第l帧的混响时间真实值。
作为上述方法的一种改进，所述第二模型的第一训练阶段的损失函数LNR满足下式：
Figure BDA0003455104330000041
其中，
Figure BDA0003455104330000042
表示估计的无噪语音幅度谱，|X|表示无噪语音幅度谱，
Figure BDA0003455104330000043
表示二范数；
所述第二模型的第二训练阶段的损失函数L满足下式：
L＝LRE+λLNE
其中，LRE为估计的混响时间与真实混响时间的均方误差，满足下式：
Figure BDA0003455104330000044
其中，λ表示微调加权因子，取为0.1，l为帧索引，T表示总帧数，
Figure BDA0003455104330000045
为第l帧的混响时间估计值，T60,l为第l帧的混响时间真实值。
一种基于深度神经网络的噪声鲁棒盲混响时间估计系统，所述系统包括预处理模块、输出模块和混响时间估计模型，其中，
所述预处理模块，用于对采集的语音信号进行预处理；
所述输出模块，用于将预处理后的语音信号输入预先建立和训练好的混响时间估计模型，得到混响时间值；
所述混响时间估计模型基于深度神经网络和时频掩膜实现。
作为上述系统的一种改进，所述混响时间估计模型为第一模型或第二模型；其中，
所述第一模型包括串联的噪声估计网络和混响时间估计网络；预处理后的语音信号同时输入噪声估计网络和混响时间估计网络，输出为语音信号的混响时间值；
所述第二模型包括串联的去噪网络和混响时间估计网络；预处理后的语音信号输入去噪网络，输出为语音信号的混响时间值。
与现有技术相比，本发明的优势在于：
1、本发明的方法利用混响语音进行房间混响时间估计，省去了传统声学测量所需要的专业的仪器与测量人员，相比于传统声学测量方法更加方便快捷；
2、本发明的方法通过收集大量混响语音数据集，采用数据驱动的深度神经网络，可以显著提高房间混响时间估计的准确性；
3、本发明的模型利用深度学习编码器，可以提取更有利于房间混响时间估计的特征，利用MG-TCNs模块，网络可以更加有效地学习输入特征在时间维度上的信息，通过进行两阶段估计，利用噪声估计和去噪网络，可以显著增加混响时间估计网络对于噪声的鲁棒性，提升网络在噪声情况下进行房间混响时间估计的准确度。
附图说明
图1是实施例1的基于深度神经网络的噪声鲁棒盲混响时间估计方法流程图；
图2是实施例1的混响时间估计模型框图；
图3是实施例2的基于深度神经网络的噪声鲁棒盲混响时间估计方法流程图。
具体实施方式
受深度神经网络在降噪方面的出色性能的启发，本方法提出了一种基于全卷积神经网络(CNNs)的噪声感知框架，以提高噪声环境下T60估计的精度。两阶段噪声感知框架由两个子网络组成，即噪声估计网络和混响时间估计网络。在第一阶段，噪声估计网络利用输入语音估计无噪声混响语音和噪声的幅度谱。在第二阶段，将第一阶段估计出的无噪声混响语音和噪声的幅度谱与输入语音的幅度谱一同输入混响时间估计网络，估计混响时间值。
本发明的方法包括以下步骤：
步骤一：准备纯净语音数据集、噪声数据集、房间冲激响应数据集。
步骤二：将纯净语音与房间冲激响应卷积，加上噪声得到训练语音。
步骤三：对训练语音进行短时傅里叶变换。
步骤四：训练噪声估计网络或去噪网络。
步骤五：训练混响时间估计网络。
步骤六：输入需要估计混响时间的语音信号，由模型输出混响时间值。
下面结合附图和实施例对本发明的技术方案进行详细的说明。
实施例1
本发明的实施例1提出了一种基于深度神经网络的噪声鲁棒盲混响时间估计方法。如图1所示，是本实施例的一种目标检测模型，简称为第一模型。
参照图1来说明本实施例的第一模型模型架构流程图。准备纯净语音数据集、噪声数据集、房间冲激响应数据集。将纯净语音与房间冲激响应卷积，加上噪声得到训练语音。对训练语音进行短时傅里叶变换。在训练阶段，训练噪声估计网络。训练混响时间估计网络。深度神经网络训练完成后，加载训练完成的神经网络参数，输入需要估计混响时间的语音信号，由模型输出混响时间值。
步骤一：准备wav格式的纯净语音数据集、噪声数据集、房间冲激响应数据集。本次采用从WSJ0-SI84数据集中的纯净语音，该数据集包含83名说话者(42名男性和41名女性)的7138个语音。分别选取5428和957个纯净语音进行训练和模型验证。采用虚源法仿真不同尺寸、不同吸收系数的长方体房间，总共生成2457个房间冲激响应。声源与麦克风之间的距离从距离D{0.7m,1m,1.7m,2m,2.5m}中随机选取。声源和麦克风之间的角度从-45°以10°的间隔增加到45°。T60的取值范围为0.2s至1.5s。
使用虚源法模拟两个未见房间，生成234个RIR作为测试集。为了减少生成的RIR带来的干扰，我们将生成的RIR与公共RIR数据库(包括OpenAIR数据库、REVERB Challenge数据集和RWCP数据库)的实测RIR相结合。将多通道RIR拆分为单通道RIR，共得到2432个单通道实测RIR。分别选取2115和317个实测RIR进行训练和验证。
从Interspeech 2020 DNS-Challenge数据集中选取6303个环境噪声作为噪声集，将5358和945个噪声文件随机分离到训练集和验证集。
步骤二：将纯净语音与房间冲激响应卷积，加上噪声得到训练语音。为了得到有噪声的混响语音信号，将纯净语音和噪声分别与随机选择的RIR进行卷积。从20dB、10dB和0dB三个不同的信噪比中随机选择一个来调整噪声的幅度。然后，将噪声添加到混响语音信号中。语音被缩短到四秒。共得到63251个噪声混响语音用于训练，4974个噪声混响语音用于验证。训练集的总长度约为70小时。
步骤三：对训练语音进行短时傅里叶变换。所有音频信号的采样率为16千赫兹。对信号进行短时傅里叶变换，使用20毫秒的汉宁窗口，帧移为10毫秒，傅里叶变换点数为320点，得到161维的频率特征。
使用Adam优化器，第一阶段学习率为0.001。第二阶段对第一阶段的网络进行微调，学习率为0.0001，对混响时间估计网络的学习率设置为0.001。为了使训练更好地收敛，如果连续3个周期的验证损失没有减少，学习率就会减半。网络训练总周期为120个，第一阶段为60个，第二阶段为60个。训练批次大小设置为8。
步骤四：训练噪声估计网络。在第一个训练阶段，噪声估计网络利用输入语音估计无噪声混响语音和噪声的幅度谱。噪声估计网络模型框图如图二所示。编码器对输入语音的幅度谱提取特征，MG-TCNs对时间序列进行时序建模，解码器将时序建模输出的隐变量还原回与幅度谱一致的维度。网络中共享编码器模块和MG-TCNs模块，分别使用两种不同的解码器估计语音掩码和噪声掩码。这样的设计的灵感来自于多任务学习，利用不同相关任务之间的信息来获得比单独处理任务更好的整体表现。由于语音掩码和噪声掩码估计是两个相关的任务，参数共享机制被期望作为一种正则化，可能产生更好的泛化能力。使用编码器和解码器对应的隐藏层使用残差连接提升网络性能。
在第一阶段使用的损失函数如下：
Figure BDA0003455104330000071
其中，α表示加权因子，取为0.5，
Figure BDA0003455104330000072
表示估计的无噪语音幅度谱，|X|表示无噪语音幅度谱，
Figure BDA0003455104330000073
表示估计的噪声幅度谱，|N|表示噪声幅度谱，
Figure BDA0003455104330000074
表示二范数。
步骤五：训练混响时间估计网络。在第二个训练阶段，将第一阶段估计出的幅度谱和输入语音的幅度谱输入混响时间估计网络，估计混响时间值。混响时间估计网络和噪声估计网络具有相似的网络拓扑结构，只是混响时间估计网络没有任何解码器并且输入特征通道的数量为3。因此省略了混响时间估计网络的详细网络框图。第二阶段对第一阶段的网络进行微调，在第二阶段使用的损失函数如下：
L＝LRE+λLNE
其中，λ表示微调加权因子，取为0.1。LRE可以用下面的式子表示：
Figure BDA0003455104330000075
其中，l是帧索引。对于每一帧，输出估计的混响时间
Figure BDA0003455104330000076
T60,l为第l帧的混响时间真实值，使用所有帧的混响时间估计值和真实值计算MSE损失。
步骤六：深度神经网络训练完成后，加载训练完成的神经网络参数，输入需要估计混响时间的语音信号，由模型输出混响时间值。
实施例2
本发明的实施例2提出了一种基于深度神经网络的噪声鲁棒盲混响时间估计方法。如图3所示，是本实施例的一种目标检测模型，简称为第二模型。
参照图3来说明本实施例的第二模型模型架构流程图。准备纯净语音数据集、噪声数据集、房间冲激响应数据集。将纯净语音与房间冲激响应卷积，加上噪声得到训练语音。对训练语音进行短时傅里叶变换。在训练阶段，训练去噪网络，训练混响时间估计网络。深度神经网络训练完成后，加载训练完成的神经网络参数，输入需要估计混响时间的语音信号，由模型输出混响时间值。
步骤一：准备wav格式的纯净语音数据集、噪声数据集、房间冲激响应数据集。本次采用从WSJ0-SI84数据集中的纯净语音，该数据集包含83名说话者(42名男性和41名女性)的7138个语音。分别选取5428和957个纯净语音进行训练和模型验证。采用虚源法仿真不同尺寸、不同吸收系数的长方体房间，总共生成2457个房间冲激响应。声源与麦克风之间的距离从距离设置D{0.7m,1m,1.7m,2m,2.5m}中随机选取。声源和麦克风之间的角度从-45°以10°的间隔增加到45°。T60的取值范围为0.2s至1.5s。
使用虚源法模拟两个未见房间，生成234个RIR作为测试集。为了减少生成的RIR带来的干扰，我们将生成的RIR与公共RIR数据库(包括OpenAIR数据库、REVERB Challenge数据集和RWCP数据库)的实测RIR相结合。将多通道RIR拆分为单通道RIR，共得到2432个单通道实测RIR。分别选取2115和317个实测RIR进行训练和验证。
从Interspeech 2020 DNS-Challenge数据集中选取6303个环境噪声作为噪声集，将5358和945个噪声文件随机分离到训练集和验证集。
步骤二：将纯净语音与房间冲激响应卷积，加上噪声得到训练语音。为了得到有噪声的混响语音信号，将纯净语音和噪声分别与随机选择的RIR进行卷积。从20dB、10dB和0dB三个不同的信噪比中随机选择一个来调整噪声的幅度。然后，将噪声添加到混响语音信号中。语音被缩短到四秒。共得到63251个噪声混响语音用于训练，4974个噪声混响语音用于验证。训练集的总长度约为70小时。
步骤三：对训练语音进行短时傅里叶变换。所有音频信号的采样率为16千赫兹。对信号进行短时傅里叶变换，使用20毫秒的汉宁窗口，帧移为10毫秒，傅里叶变换点数为320点，得到161维的频率特征。
使用Adam优化器，第一阶段学习率为0.001。第二阶段对第一阶段的网络进行微调，学习率为0.0001，对混响时间估计网络的学习率设置为0.001。为了使训练更好地收敛，如果连续3个周期的验证损失没有减少，学习率就会减半。网络训练总周期为120个，第一阶段为60个，第二阶段为60个。训练批次大小设置为8。
步骤四：训练去噪网络。在第一个训练阶段，去噪网络利用输入语音估计无噪声混响语音的幅度谱。编码器对输入语音的幅度谱提取特征，MG-TCNs对时间序列进行时序建模，解码器将时序建模输出的隐变量还原回与幅度谱一致的维度。使用编码器和解码器对应的隐藏层使用残差连接提升网络性能。
在第一阶段使用的损失函数如下：
Figure BDA0003455104330000081
其中，
Figure BDA0003455104330000091
表示估计的无噪语音幅度谱，|X|表示无噪语音幅度谱，
Figure BDA0003455104330000092
表示二范数。
步骤五：训练混响时间估计网络。在第二个训练阶段，将第一阶段估计出的幅度谱和输入语音的幅度谱输入混响时间估计网络，估计混响时间值。混响时间估计网络和去噪网络具有相似的网络拓扑结构，只是混响时间估计网络没有任何解码器并且输入特征通道的数量为1。因此省略了混响时间估计网络的详细网络框图。
第二阶段对第一阶段的网络进行微调，在第二阶段使用的损失函数如下：
L＝LRE+λLNR
其中，λ表示微调加权因子，取为0.1。LRE可以用下面的式子表示：
Figure BDA0003455104330000093
其中，l是帧索引。对于每一帧，输出估计的混响时间
Figure BDA0003455104330000094
T60,l为第l帧的混响时间真实值，使用所有帧的混响时间估计值和真实值计算MSE损失。
步骤六：深度神经网络训练完成后，加载训练完成的神经网络参数，输入需要估计混响时间的语音信号，由模型输出混响时间值。
实施例3
本发明的实施例3提出了一种基于深度神经网络的噪声鲁棒盲混响时间估计系统，基于实施例1的方法实现，该系统包括预处理模块、输出模块和混响时间估计模型，其中，
所述预处理模块，用于对采集的语音信号进行预处理；
所述输出模块，用于将预处理后的语音信号输入预先建立和训练好的混响时间估计模型，得到混响时间值；
所述混响时间估计模型基于深度神经网络和时频掩膜实现。该模型为第一模型，包括串联的噪声估计网络和混响时间估计网络；预处理后的语音信号同时输入噪声估计网络和混响时间估计网络，输出为语音信号的混响时间值；
具体模型结构及训练步骤详见实施例1的说明。
实施例4
本发明的实施例4提出了一种基于深度神经网络的噪声鲁棒盲混响时间估计系统，基于实施例2的方法实现，该系统包括预处理模块、输出模块和混响时间估计模型，其中，
所述预处理模块，用于对采集的语音信号进行预处理；
所述输出模块，用于将预处理后的语音信号输入预先建立和训练好的混响时间估计模型，得到混响时间值；
所述混响时间估计模型基于深度神经网络和时频掩膜实现。该模型为第二模型，包括串联的去噪网络和混响时间估计网络；预处理后的语音信号输入去噪网络，输出为语音信号的混响时间值。
具体模型结构及训练步骤详见实施例2的说明。
最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (10)

1.一种基于深度神经网络的噪声鲁棒盲混响时间估计方法，所述方法包括：
对采集的语音信号进行预处理；
将预处理后的语音信号输入预先建立和训练好的混响时间估计模型，得到混响时间值；
所述混响时间估计模型基于深度神经网络和时频掩膜实现。
2.根据权利要求1所述的基于深度神经网络的噪声鲁棒盲混响时间估计方法，其特征在于，所述对采集的语音信号进行预处理；具体包括：
使用设定毫秒的汉宁窗口，根据设定的帧移和傅里叶变换点数，对采集的语音信号进行短时傅里叶变换，得到多维的频率特征。
3.根据权利要求2所述的基于深度神经网络的噪声鲁棒盲混响时间估计方法，其特征在于，所述混响时间估计模型为第一模型或第二模型；其中，
所述第一模型包括串联的噪声估计网络和混响时间估计网络；预处理后的语音信号同时输入噪声估计网络和混响时间估计网络，输出为语音信号的混响时间值；
所述第二模型包括串联的去噪网络和混响时间估计网络；预处理后的语音信号输入去噪网络，输出为语音信号的混响时间值。
4.根据权利要求3所述的基于深度神经网络的噪声鲁棒盲混响时间估计方法，其特征在于，
所述第一模型的噪声估计网络包括编码器、MG-TCNs模块和解码器，其中，
所述编码器，用于对预处理后的语音信号提取隐变量特征；
所述MG-TCNs模块，用于对编码器提取的隐变量特征的时间序列进行时序建模，输出隐变量；
所述解码器包括并联的语音解码器和噪声解码器，分别用于将MG-TCNs模块输出的隐变量还原为与语音信号幅度谱维度一致的无噪声语音估计和噪声估计；
所述第一模型的混响时间估计网络，用于结合无噪声语音估计和噪声估计对预处理后的语音信号估算出混响时间值，所述混响时间估计网络包括编码器和MG-TCNs模块，均与噪声估计网络对应模块的结构相同。
5.根据权利要求4所述的基于深度神经网络的噪声鲁棒盲混响时间估计方法，其特征在于，
所述第二模型的去噪网络包括编码器、MG-TCNs模块和语音解码器，其中，
所述编码器，用于对预处理后的语音信号提取隐变量特征；
所述MG-TCNs，用于对编码器提取的隐变量特征的时间序列进行时序建模，输出隐变量；
所述语音解码器用于将隐变量还原为与语音信号幅度谱维度一致的无噪声语音估计；
所述第二模型的混响时间估计网络，用于根据无噪声语音估计估算出混响时间值，所述混响时间估计网络包括编码器和MG-TCNs模块，均与去噪网络对应模块的结构相同。
6.根据权利要求5所述的基于深度神经网络的噪声鲁棒盲混响时间估计方法，其特征在于，所述方法还包括混响时间估计模型的训练步骤；具体包括：
准备纯净语音数据集、噪声数据集和房间冲激响应数据集，构成样本集；
从样本集中获取纯净语音和房间冲激响应，进行卷积并加上噪声得到训练语音构成训练集；
使用设定毫秒的汉宁窗口，根据设定的帧移和傅里叶变换点数，对训练集中的训练语音依次进行短时傅里叶变换，分别得到对应的多维的频率特征；
使用Adam优化器，设置两个训练阶段，分别设置每个训练阶段的学习率和训练批次及周期，将预处理后的训练集数据分别输入第一模型或第二模型，直至满足训练要求，得到训练好的第一模型或第二模型，从而完成混响时间估计模型的训练。
7.根据权利要求6所述的基于深度神经网络的噪声鲁棒盲混响时间估计方法，其特征在于，所述第一模型的第一训练阶段的损失函数LNE满足下式：
Figure FDA0003455104320000021
其中，α表示加权因子，取为0.5，
Figure FDA0003455104320000022
表示估计的无噪语音幅度谱，|X|表示无噪语音幅度谱，
Figure FDA0003455104320000023
表示估计的噪声幅度谱，|N|表示噪声幅度谱，
Figure FDA0003455104320000024
表示二范数；
所述第一模型的第二训练阶段的损失函数L满足下式：
L＝LRE+λLNE
其中，LRE为估计的混响时间与真实混响时间的均方误差，满足下式：
Figure FDA0003455104320000025
其中，λ表示微调加权因子，取为0.1，l为帧索引，T表示总帧数，
Figure FDA0003455104320000026
为第l帧的混响时间估计值，T60,l为第l帧的混响时间真实值。
8.根据权利要求6所述的基于深度神经网络的噪声鲁棒盲混响时间估计方法，其特征在于，所述第二模型的第一训练阶段的损失函数LNR满足下式：
Figure FDA0003455104320000031
其中，
Figure FDA0003455104320000032
表示估计的无噪语音幅度谱，|X|表示无噪语音幅度谱，
Figure FDA0003455104320000033
表示二范数；
所述第二模型的第二训练阶段的损失函数L满足下式：
L＝LRE+λLNE
其中，LRE为估计的混响时间与真实混响时间的均方误差，满足下式：
Figure FDA0003455104320000034
其中，λ表示微调加权因子，取为0.1，l为帧索引，T表示总帧数，
Figure FDA0003455104320000035
为第l帧的混响时间估计值，T60,l为第l帧的混响时间真实值。
9.一种基于深度神经网络的噪声鲁棒盲混响时间估计系统，其特征在于，所述系统包括预处理模块、输出模块和混响时间估计模型，其中，
所述预处理模块，用于对采集的语音信号进行预处理；
所述输出模块，用于将预处理后的语音信号输入预先建立和训练好的混响时间估计模型，得到混响时间值；
所述混响时间估计模型基于深度神经网络和时频掩膜实现。
10.根据权利要求9所述的基于深度神经网络的噪声鲁棒盲混响时间估计系统，其特征在于，所述混响时间估计模型为第一模型或第二模型；其中，
所述第一模型包括串联的噪声估计网络和混响时间估计网络；预处理后的语音信号同时输入噪声估计网络和混响时间估计网络，输出为语音信号的混响时间值；
所述第二模型包括串联的去噪网络和混响时间估计网络；预处理后的语音信号输入去噪网络，输出为语音信号的混响时间值。
CN202210004710.6A 2022-01-04 2022-01-04 一种基于深度神经网络的噪声鲁棒盲混响时间估计方法 Active CN114255780B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210004710.6A CN114255780B (zh) 2022-01-04 2022-01-04 一种基于深度神经网络的噪声鲁棒盲混响时间估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210004710.6A CN114255780B (zh) 2022-01-04 2022-01-04 一种基于深度神经网络的噪声鲁棒盲混响时间估计方法

Publications (2)

Publication Number Publication Date
CN114255780A true CN114255780A (zh) 2022-03-29
CN114255780B CN114255780B (zh) 2023-07-25

Family

ID=80799360

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210004710.6A Active CN114255780B (zh) 2022-01-04 2022-01-04 一种基于深度神经网络的噪声鲁棒盲混响时间估计方法

Country Status (1)

Country Link
CN (1) CN114255780B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115331687A (zh) * 2022-08-03 2022-11-11 北京声智科技有限公司 样本音频数据的生成方法、装置、电子设备及存储介质
CN115798497A (zh) * 2022-11-03 2023-03-14 中国科学院声学研究所 一种时延估计系统及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200082843A1 (en) * 2016-12-15 2020-03-12 Industry-University Cooperation Foundation Hanyang University Multichannel microphone-based reverberation time estimation method and device which use deep neural network
CN112687284A (zh) * 2020-12-21 2021-04-20 中国科学院声学研究所 混响语音的混响抑制方法及装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200082843A1 (en) * 2016-12-15 2020-03-12 Industry-University Cooperation Foundation Hanyang University Multichannel microphone-based reverberation time estimation method and device which use deep neural network
CN112687284A (zh) * 2020-12-21 2021-04-20 中国科学院声学研究所 混响语音的混响抑制方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
万新旺，吴镇扬: "基于自适应频率选择的鲁棒时延估计算法", 《东南大学学报（自然科学版）》, vol. 40, no. 5 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115331687A (zh) * 2022-08-03 2022-11-11 北京声智科技有限公司 样本音频数据的生成方法、装置、电子设备及存储介质
CN115798497A (zh) * 2022-11-03 2023-03-14 中国科学院声学研究所 一种时延估计系统及装置
CN115798497B (zh) * 2022-11-03 2023-10-13 中国科学院声学研究所 一种时延估计系统及装置

Also Published As

Publication number Publication date
CN114255780B (zh) 2023-07-25

Similar Documents

Publication Publication Date Title
Lebart et al. A new method based on spectral subtraction for speech dereverberation
US8160273B2 (en) Systems, methods, and apparatus for signal separation using data driven techniques
CN108447495B (zh) 一种基于综合特征集的深度学习语音增强方法
JP5587396B2 (ja) 信号分離のためのシステム、方法、および装置
CN112735456B (zh) 一种基于dnn-clstm网络的语音增强方法
CN107452389A (zh) 一种通用的单声道实时降噪方法
KR20180115984A (ko) 심화신경망 기반의 잡음 및 에코의 통합 제거 방법 및 장치
CN103854662A (zh) 基于多域联合估计的自适应语音检测方法
CN105448302B (zh) 一种环境自适应的语音混响消除方法和系统
Ratnarajah et al. Towards improved room impulse response estimation for speech recognition
CN110767244A (zh) 语音增强方法
CN111243617B (zh) 一种基于深度学习降低mfcc特征失真的语音增强方法
CN110660406A (zh) 近距离交谈场景下双麦克风移动电话的实时语音降噪方法
CN115424627A (zh) 基于卷积循环网络和wpe算法的语音增强混合处理方法
WO2019014890A1 (zh) 一种通用的单声道实时降噪方法
CN110808057A (zh) 一种基于约束朴素生成对抗网络的语音增强方法
Lv et al. A permutation algorithm based on dynamic time warping in speech frequency-domain blind source separation
Looney et al. Joint estimation of acoustic parameters from single-microphone speech observations
Callens et al. Joint blind room acoustic characterization from speech and music signals using convolutional recurrent neural networks
CN116013344A (zh) 一种多种噪声环境下的语音增强方法
CN114255780B (zh) 一种基于深度神经网络的噪声鲁棒盲混响时间估计方法
Zheng et al. Noise-robust blind reverberation time estimation using noise-aware time–frequency masking
Astudillo et al. Uncertainty propagation
CN109215635B (zh) 用于语音清晰度增强的宽带语音频谱倾斜度特征参数重建方法
Raikar et al. Effect of Microphone Position Measurement Error on RIR and its Impact on Speech Intelligibility and Quality.

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant