CN110289011A - 一种用于分布式无线声传感器网络的语音增强系统 - Google Patents
一种用于分布式无线声传感器网络的语音增强系统 Download PDFInfo
- Publication number
- CN110289011A CN110289011A CN201910649352.2A CN201910649352A CN110289011A CN 110289011 A CN110289011 A CN 110289011A CN 201910649352 A CN201910649352 A CN 201910649352A CN 110289011 A CN110289011 A CN 110289011A
- Authority
- CN
- China
- Prior art keywords
- signal
- module
- wavelet
- layer
- noise power
- 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
Links
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 28
- 230000009466 transformation Effects 0.000 claims abstract description 21
- 230000003595 spectral effect Effects 0.000 claims description 20
- 238000000354 decomposition reaction Methods 0.000 claims description 17
- 238000012545 processing Methods 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 230000010363 phase shift Effects 0.000 claims description 2
- 238000000034 method Methods 0.000 description 16
- 238000010586 diagram Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000000875 corresponding effect Effects 0.000 description 3
- 125000002015 acyclic group Chemical group 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 102000016917 Complement C1 Human genes 0.000 description 1
- 108010028774 Complement C1 Proteins 0.000 description 1
- 241000712899 Lymphocytic choriomeningitis mammarenavirus Species 0.000 description 1
- 241000209140 Triticum Species 0.000 description 1
- 235000021307 Triticum Nutrition 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000037433 frameshift Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
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
- 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/0264—Noise filtering characterised by the type of parameter measurement, e.g. correlation techniques, zero crossing techniques or predictive techniques
-
- 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
- G10L2021/02161—Number of inputs available containing the signal or the noise to be suppressed
- G10L2021/02166—Microphone arrays; Beamforming
Landscapes
- Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Quality & Reliability (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Circuit For Audible Band Transducer (AREA)
Abstract
本发明公开了一种用于分布式无线声传感器网络的语音增强系统,具体包括:相位对齐模块、解析变换模块、小波变换模块、NPSD估计模块、MVDR滤波器模块、分布式算法迭代模块和小波反变换模块。采用该系统首先将每个节点的输入信号要和距离声源最远的那个节点相位对齐,对齐之后的每一路信号分别进行解析小波变换;然后,在小波域中估计NPSD,同时得到MVDR滤波系数;最终,通过分布式算法迭代模块得到每个节点麦克风的输出信号,对输出信号进行小波反变换并取实部获得时域信号。该系统的实现是由每个麦克风和其附近麦克风通过交换特定数据,使得所有麦克风的初始状态值经过迭代收敛到全局平均值完成的。
Description
技术领域
本发明涉及语音信号处理技术领域,尤其涉及一种用于分布式无线声传感器网络的语音增强系统。
背景技术
在实际应用中,音频处理设备接收的音频信号经常受到各种噪声的干扰。干扰噪声会使设备性能下降,严重情况下,设备甚至无法正常工作。为了克服噪声的影响,人们提出了语音增强技术。由于基于单个麦克风的语音增强方法性能有限,可以利用空间信息的麦克风阵列语音增强方法获得重视。但是麦克风阵列也暴露了需要知道阵列几何信息,要有固定且规则的阵列结构等一些局限性。随着无线传感器技术的飞速发展,无线声传感器网络(WASN)的应用越来越广泛,如在会议室利用智能手机或笔记本电脑通过WiFi(或蓝牙)和音频接口构建出一个WASN。当WASN没有数据融合中心时,可以利用声传感器的无线通信功能进行数据传输实现语音增强。
现有技术中:Zeng等在文献[1]研究了分布式语音增强Gossip算法的使用,提出一种分布式延迟求和波束形成语音增强技术。该波束形成是在任意连接的 WASN下,通过每个麦克风与其相邻的麦克风通信达到估计期望信号的目的。此技术通过异步更新数据和分布式求平均两种方法同时更新每个麦克风的数据,从而提高WASN分布式语音增强的性能。但是该方案的缺陷是由于只注重了分布式算法的实现,致使最终的结果虽然达到了全局最优解,但是该算法中的噪声功率谱密度(NPSD)估计利用的是频域已有算法,保证性能的前提下估计成本较高,当噪声环境恶劣时,估计准确率降低。NPSD的估计误差越大使输出结果的性能下降越多。
现有技术中:Sherson等在文献[2]提出了一种分布式最小方差算法,该算法不需要每个节点分享原始数据,在循环和非循环网络中都可以计算出最优波束形成的输出结果。NPSD由时变噪声场的协方差矩阵的低秩结构估计得到。该算法数据整体传输数量要比集中式算法少。在循环网络的任意拓扑结构下该算法的性能促进了最佳声学波束形成的使用。但是该方法的缺陷是:在快速变化的噪声场中,需要按块计算波束形成输出。在非循环网络中,分布式最小方差算法传输数据量相对较高。
现有技术中:Xiao等在文献[3]为了使NPSD估计更加精确从而提高现有波束形成最小方差无失真响应(MVDR)算法的鲁棒性,在阵列信号波束形成模型中提出了基于估计对角加载的MVDR算法。建立了基于用拉格朗日乘子求对角线加载补偿值的MVDR优化模型,并在矩阵理论的基础上推导出对角线加载补偿值的区间。通过实验方法也确定了区间中的最优对角加载值,并在稳定性和抑制干扰加性噪声方面有了一定提高。该算法是在均匀线性阵列模型中研究的MVDR波束形成,应用在分布式算法中不具有低复杂度特性。
发明内容
根据现有技术存在的问题,本发明公开了一种用于分布式无线声传感器网络的语音增强系统,具体包括:
相位对齐模块、将先后接收到的语音信号yi”(k)进行相位平移形成同一时刻语音信号yi'(k);
接收所述相位对齐模块传送的语音信号yi'(k)对该信号进行分帧加窗处理得到信号yi(k)的解析变换模块,所述解析变换模块对每一帧信号yi(k)进行解析变换、输出信号yi;
接收解析变换模块输出的信号yi的小波变换模块,所述小波变换模块采用四层小波分解方式对信号yi进行滤波和降采样处理得到解析的离散小波域信号 Y;
接收所述小波变换模块输出的含有四层分解的信号Y的NPSD估计模块,所述NPSD估计模块利用每一层数据中的阈值对语音信号的噪声功率谱密度分别进行估计、输出每一层的噪声功率谱密度;
接收所述NPSD估计模块传送的噪声功率谱密度信息、以及小波变换模块传送的离散小波域信号Y的MVDR滤波器模块,所述MVDR滤波器模块根据接收到的数据信息构建MVDR准则下的滤波器系数W、将该滤波器系数W与离散小波域信号Y组合成滤波信号WHY;
分布式算法迭代模块,将MVDR滤波器模块输出的滤波信号WHY处理为求平均值的形式、通过初始状态值根据Metropolis权矩阵多次迭代求得初始状态值的平均值得到输出信号Z;
小波反变换模块:将分布式算法迭代模块传送的输出信号Z进行小波反变换、再取实部得到增强后的时域输出语音信号。
进一步的,所述NPSD估计模块对语音信号的每层的噪声功率谱密度进行估计时采用如下方式:
设定每次分解时的阈值,其中阈值表达式如下
其中,j代表小波分解中的第j层,Tj是第j层小波系数的阈值,Nj是第j层小波系数的个数,σj的表达式如下
其中,median(|dj|)是第j层小波系数的中位数,采用如下公式估计语音信号的每层噪声功率谱密度,当前帧的噪声功率谱密度估计值均由上一帧的噪声功率谱密度估计值和当前帧的阈值更新得到,
其中,L是当前帧数,是第i个麦克风的输入信号计算到L帧时在小波域中第j层的噪声功率谱密度估计值,是第l帧数字信号对应的解析小波变换中第j层的阈值,
由于采用了上述技术方案,本发明提供的一种用于分布式无线声传感器网络的语音增强系统,采用该系统首先将每个节点的输入信号要和距离声源最远的那个节点相位对齐,对齐之后的每一路信号分别进行解析小波变换;然后,在小波域中估计NPSD,同时得到MVDR滤波系数;最终通过分布式算法迭代模块得到每个节点麦克风的输出信号,对输出信号进行小波反变换并取实部获得时域信号。该系统的实现是由每个麦克风和其附近麦克风通过交换特定数据,使得所有麦克风的初始状态值经过迭代收敛到全局平均值完成的。其中,在设置迭代次数足够的情况下,每个麦克风的输出信号是相同的,均得到输出信号的方差最小,并达到全局最优解的效果。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明系统的工作原理图;
图2为本发明无线声传感器网络示意图;
图3为本发明四层小波分解图示意图;
图4为本发明分布式语音增强前每个麦克风输入信号的SNR示意图;
图5为本发明分布式语音增强前每个麦克风输入信号的PESQ示意图;
图6为本发明分布式语音增强后每个麦克风输出信号的SNR和PESQ示意图;
具体实施方式
为使本发明的技术方案和优点更加清楚,下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚完整的描述:
如图1所示的一种用于分布式无线声传感器网络的语音增强系统,具体包括:相位对齐模块、解析变换模块、小波变换模块、NPSD估计模块、MVDR 滤波器模块、分布式算法迭代模块和小波反变换模块。
其中实际的无线声传感器网络WASN如图2所示,麦克风和声源的位置随机分布。麦克风之间的通信受到距离的限制,只能和附近的麦克风进行通信。
所述相位对齐模块在无线声传感器网络WASN中,在距声源某个已知距离为l的位置放一个参考麦克风,利用该麦克风接收到的信号能量与其它麦克风接收的信号能量可以估计出其它麦克风到声源的距离li。距离估计公式如下:
其中,E和Ei分别是参考信号和其它麦克风信号的能量,ε和εi是背景噪声的能量,能量公式如下:
其中,len是信号的总长度,fs是采样频率,即一秒信号对应的点数。该公式利用语音第一秒大部分是无语音段的特点,估计背景噪声的能量。定义max{li}对应节点为参考输入信号y”a(k),每一路待对齐信号为y”b(k),且y”b(k)循环通过一个单位延迟器的同时,与y”a(k)做互相关运算,表达式如下:
Rab(τ)=E[y”a(k)y”b(k-τ)],τ=0,1,...,T (3)
其中,T是平移最大量,可适当选取。当τ的取值使以上两路信号对齐时,互相关函数值最大。令script{}为取τ值运算,则待对齐信号与参考信号对齐的输出信号为
根据无线声传感器网络连接图可依次将各路信号对齐。由于麦克风与声源的距离不同,导致每个麦克风接收到信号的先后顺序不一致,相位对齐模块的作用是选择一路参考信号,将其他每一路信号的相位进行一定数量的平移,使得所有麦克风接收到的信号相当于是同一时刻接收的。
解析变换模块是将相位对齐模块输出的各路信号分别进行分帧加窗处理,对每一帧信号进行解析变换。在验证时用的是汉宁窗,帧移为50%,每帧数据长度为K=256点。汉宁窗表达式如下:
ω(k)=0.5-0.5cos(2πk/K),k=0,1,...,K-1 (5)
加窗后的信号为
y(k)=y'(k)ω(k) (6)
信号变为解析信号表示为
y=y+jH[y] (7)
其中,H[]表示希尔伯特变换。希尔伯特变换系数表如下:
当序列数为偶数时,系数均为0值。且后半部分系数有如下关系:h(128)=0, h(129)=-h(127),...,h(256)=-h(1),即后一半系数关于h(128)奇对称。希尔伯特变换为
其中,H[y(k)]的长度为511点,先取前K=256个数,后256-1个数前位对齐相加,即
使得H[y(k)]的长度变为K个数。然后将H[y(k)]作虚部,y(k)作实部,更新y 为解析信号。
小波变换模块采用Mallat快速算法的原理,离散小波变换可以看成是由一个滤波器组对信号进行滤波和降采样处理。信号的四层小波分解图如图3所示,这里分别给出了‘sym4’母小波的高通滤波器分解系数和低通滤波器分解系数:
H0=[-0.0322 -0.0126 0.0992 0.2979 -0.8037 0.4976 0.0296 -0.0758];
G0=[-0.0758 -0.0296 0.4976 0.8037 0.2979 -0.0992 -0.0126 0.0322]。
令解析变换模块输出的信号y分别与高通滤波器分解系数H0和低通滤波器分解系数G0作卷积运算,结果如下:
得到(K+M-1)个数,其中M=8是滤波器系数的个数。得到的卷积结果取前K个数,后M-1个数前位对齐相加,即
分别得到点数为K的高通滤波系数和低通滤波系数。然后进行降采样运算,在 K个数中,以2为步长抽取K/2个数,具体操作如下:
分别得到高频系数HI_D1和低频系数LO_D1。第二层分解时,对低频系数重复上述操作。图4给出了四层小波分解的结构图。本专利在验证时是做了四层分解。将解析信号的离散小波变换表示为Y=DW[y+jH[y]]。
NPSD估计模块采用了一种新的噪声功率谱密度NPSD估计方法,称为阈值法。该方法利用了小波域中提到的阈值,可以避免使用VAD(语音活动检测) 方式。这里使用自适应的阈值对噪声功率谱密度NPSD进行估计,阈值表达式如下
其中,下标j代表小波分解中的第j层,Tj是第j层小波系数的阈值。Nj是第j 层小波系数的个数,σj的表达式如下
其中,median(|dj|)是第j层小波系数的中位数。每一帧用到的NPSD估计值均由上一帧的噪声功率谱密度NPSD估计值和当前帧的阈值更新得到,那么NPSD 的估计表达式为
其中,L是当前帧数,是第i个麦克风的输入信号计算到L帧时在小波域中第j 层的NPSD估计值,是第l帧数字信号对应的解析小波变换中第j层的阈值,
进一步的,为了保证输出信号方差最小的同时没有失真,MVDR滤波器模块要在解析信号的小波域找到一个滤波器系数W,使输出信号Z有如下表达式
Z=WHY (16)
其中,[]H表示向量或矩阵的共轭转置。通过MVDR可得滤波器系数W,即
其中,a=[a1,a2,…,aN]H,δi 2是第i个麦克风的NPSD,ai=0.61/li是衰减系数估计。
分布式算法迭代模块的工作原理是:将MVDR滤波器模块输出的滤波信号 WHY处理为求平均值的形式、通过初始状态值根据Metropolis权矩阵多次迭代求得初始状态值的平均值得到输出信号。其中在分布式实现该算法时,需要将输出信号Z的表达形式写成求平均值的形式,该形式经定义式可变换为
其中,令和上式重新写为
观察上式可以得知分布式MVDR的结果只需要每个麦克风获得所有麦克风初始状态值的平均值(和),就可以获得和上式同样的输出结果。在分布式算法下,是通过令初始状态值不断迭代来求得初始状态值的平均值,迭代公式如下
其中,M是Metropolis权矩阵,定义如下
式中D代表两个不同麦克风可以互相通信的连接集合,即 (i,j)∈D(i,j=1,…,N,i≠j)。di表示第i个麦克风可以与附近麦克风通信的数量。在上述迭代计算使得每个麦克风的输出信号收敛时,达到全局最优解。本专利在验证时设置迭代次数上限为100次,当迭代到上限时默认收敛。
进一步的,分布式算法迭代模块的输出是解析信号小波域中的解。在小波反变换模块中,我们需要将输入值进行小波反变换,再取实部得到时域输出语音信号。离散小波反变换的具体操作与离散小波变换的做法正好相反,先升采样再进行滤波,最终数据相加实现重构。这里给出了‘sym4’母小波的重构滤波器系数:
H0'=[-0.0758 0.0296 0.4976 -0.8037 0.2979 0.0992 -0.0126 -0.0322];
G0'=[0.0322 -0.0126 -0.0992 0.2979 0.8037 0.4976 -0.0296 -0.0758]。
升采样即先对最后一层的高频系数HI_D4和低频系数LO_D4分别每相邻两个点之间补1个零,具体运算如下:
升采样后的高频系数HI_R4”和低频系数LO_R”分别与重构高通滤波器系数H0’和重构低通滤波器系数G0’做卷积如下;
将两组卷积后的数分别取后M-1个数末位对齐相加如下:
最后将两组数据相加,即
LO_D3=HI_R4(k)+LO_R4(k) (25)
就实现了第四层高频系数HI_D4和低频系数LO_D4的小波重构,重构回第三层的低频系数LO_D3,依次类推,即可实现小波变换重构。
为验证本发明方法的有效性,本专利通过Imgae模型模拟了一个5×5×3的封闭式房间,混响时间为0.1秒,声源位置坐标为(2.5米,2.5米,1米),取10个麦克风,麦克风横纵坐标利用均匀分布随机生成且高度也均为1米。麦克风之间的通信距离是在所有通信线路有45条的情况下,选择了最短的19条线路作为通信情况。声源是从TIMIT数据库[https://download.csdn.net/download/sdhyfxh/4086482]随机选取的一个6 秒纯净语音信号,采样频率是16kHz。在每个麦克风接收到输入语音信号的同时加上一个互不相关的高斯白噪声作为输入噪声信号,致使10个麦克风输入信号的SNR值和PESQ值如图4、图5所示。可以看出,每个麦克风输入信号的SNR值和PESQ值各不相同。
采用本方法所提出的分布式MVDR的语音增强方法对每个麦克风输入信号进行降噪,每个麦克风的输出信号均达到全局最优解。消噪后的每个麦克风输出信号的SNR值和PESQ值如图6所示。可以看出SNR值和PESQ值相对文献[1]中的方法都有所提高,从而达到了分布式语音增强的效果。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
文件中的字母缩写的含义
MVDR:最小方差无失真响应
PESQ:语音质量的感知评估 SNR:信噪比 VAD:语音活动检测
NPSD:噪声功率谱密度 WASN:无线声传感器网络
[1]Y.Zeng and R.C.Hendriks,"Distributed Delay and Sum Beamformer forSpeech Enhancement via Randomized Gossip,"in IEEE/ACM Transactions on Audio,Speech,and Language Processing,vol.22,no.1,pp.260-273,Jan.2014.
[2]T.Sherson,W.B.Kleijn and R.Heusdens,"A distributed algorithm forrobust LCMV beamforming,"2016IEEE International Conference on Acoustics,Speech and Signal Processing(ICASSP),Shanghai,2016,pp.101-105.
[3]Yuteng Xiao,Jihang Yin,Honggang Qi,Hongsheng Yin,and Gang Hua,“MVDR Algorithm Based on Estimated Diagonal Loading for Beamforming,”Mathematical Problems in Engineering,vol.2017,Article ID 7904356,7pages,2017。
Claims (2)
1.一种用于分布式无线声传感器网络的语音增强系统,其特征在于包括:
相位对齐模块、将先后接收到的语音信号y”i(k)进行相位平移形成同一时刻语音信号y'i(k);
接收所述相位对齐模块传送的语音信号y'i(k)对该信号进行分帧加窗处理得到信号yi(k)的解析变换模块,所述解析变换模块对每一帧信号yi(k)进行解析变换、输出信号yi;
接收解析变换模块输出的信号yi的小波变换模块,所述小波变换模块采用四层小波分解方式对信号yi进行滤波和降采样处理得到解析的离散小波域信号Y;
接收所述小波变换模块输出的含有四层分解的信号Y的NPSD估计模块,所述NPSD估计模块利用每一层数据中的阈值对语音信号的噪声功率谱密度分别进行估计、输出每一层的噪声功率谱密度;
接收所述NPSD估计模块传送的噪声功率谱密度信息、以及小波变换模块传送的离散小波域信号Y的MVDR滤波器模块,所述MVDR滤波器模块根据接收到的数据信息构建MVDR准则下的滤波器系数W、将该滤波器系数W与离散小波域信号Y组合成滤波信号WHY;
分布式算法迭代模块,将MVDR滤波器模块输出的滤波信号WHY处理为求平均值的形式、通过初始状态值根据Metropolis权矩阵多次迭代求得初始状态值的平均值得到输出信号Z;
小波反变换模块:将分布式算法迭代模块传送的输出信号Z进行小波反变换、再取实部得到增强后的时域输出语音信号。
2.根据权利要求1所述的一种用于分布式无线声传感器网络的语音增强,其特征还在于:所述NPSD估计模块对语音信号的每层的噪声功率谱密度进行估计时采用如下方式:
设定每次分解时的阈值,其中阈值表达式如下
其中,j代表小波分解中的第j层,Tj是第j层小波系数的阈值,Nj是第j层小波系数的个数,σj的表达式如下
其中,median(|dj|)是第j层小波系数的中位数,采用如下公式估计语音信号的每层噪声功率谱密度,当前帧的噪声功率谱密度估计值均由上一帧的噪声功率谱密度估计值和当前帧的阈值更新得到,
其中,L是当前帧数,是第i个麦克风的输入信号计算到L帧时在小波域中第j层的噪声功率谱密度估计值,是第l帧数字信号对应的解析小波变换中第j层的阈值,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910649352.2A CN110289011B (zh) | 2019-07-18 | 2019-07-18 | 一种用于分布式无线声传感器网络的语音增强系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910649352.2A CN110289011B (zh) | 2019-07-18 | 2019-07-18 | 一种用于分布式无线声传感器网络的语音增强系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110289011A true CN110289011A (zh) | 2019-09-27 |
CN110289011B CN110289011B (zh) | 2021-06-25 |
Family
ID=68023319
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910649352.2A Active CN110289011B (zh) | 2019-07-18 | 2019-07-18 | 一种用于分布式无线声传感器网络的语音增强系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110289011B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110739004A (zh) * | 2019-10-25 | 2020-01-31 | 大连理工大学 | 一种用于wasn的分布式语音噪声消除系统 |
CN111696575A (zh) * | 2020-06-19 | 2020-09-22 | 杭州电子科技大学 | 基于混合神经网络模型的低通气和呼吸暂停检测识别系统 |
WO2021138201A1 (en) * | 2019-12-30 | 2021-07-08 | Texas Instruments Incorporated | Background noise estimation and voice activity detection system |
CN114449410A (zh) * | 2021-11-12 | 2022-05-06 | 山东大学 | 一种多通道声纹信号同步采集系统及方法 |
CN114609493A (zh) * | 2022-05-09 | 2022-06-10 | 杭州兆华电子股份有限公司 | 一种信号数据增强的局部放电信号识别方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102347028A (zh) * | 2011-07-14 | 2012-02-08 | 瑞声声学科技(深圳)有限公司 | 双麦克风语音增强装置及方法 |
CN104952459A (zh) * | 2015-04-29 | 2015-09-30 | 大连理工大学 | 一种基于分布式一致性和mvdr波束形成的分布式语音增强方法 |
CN108172231A (zh) * | 2017-12-07 | 2018-06-15 | 中国科学院声学研究所 | 一种基于卡尔曼滤波的去混响方法及系统 |
CN109308904A (zh) * | 2018-10-22 | 2019-02-05 | 上海声瀚信息科技有限公司 | 一种阵列语音增强算法 |
US20190043515A1 (en) * | 2018-07-12 | 2019-02-07 | Intel IP Corporation | Joint beamforming and echo cancellation for reduction of noise and non-linear echo |
-
2019
- 2019-07-18 CN CN201910649352.2A patent/CN110289011B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102347028A (zh) * | 2011-07-14 | 2012-02-08 | 瑞声声学科技(深圳)有限公司 | 双麦克风语音增强装置及方法 |
CN104952459A (zh) * | 2015-04-29 | 2015-09-30 | 大连理工大学 | 一种基于分布式一致性和mvdr波束形成的分布式语音增强方法 |
CN108172231A (zh) * | 2017-12-07 | 2018-06-15 | 中国科学院声学研究所 | 一种基于卡尔曼滤波的去混响方法及系统 |
US20190043515A1 (en) * | 2018-07-12 | 2019-02-07 | Intel IP Corporation | Joint beamforming and echo cancellation for reduction of noise and non-linear echo |
CN109308904A (zh) * | 2018-10-22 | 2019-02-05 | 上海声瀚信息科技有限公司 | 一种阵列语音增强算法 |
Non-Patent Citations (3)
Title |
---|
LAE-HOON KIM ET AL.: "《Generalized Optimal Multi-Microphone Speech Enhancement Using Sequential Minimum Variance Distortionless Response(MVDR) Beamforming and Postfiltering》", 《 2006 IEEE INTERNATIONAL CONFERENCE ON ACOUSTICS SPEECH AND SIGNAL PROCESSING PROCEEDINGS》 * |
MAHDI PARCHAMI ET AL.: "《Recent Developments in Speech Enhancement in the Short-Time Fourier Transform Domain》", 《 IEEE CIRCUITS AND SYSTEMS MAGAZINE ( VOLUME: 16, ISSUE: 3, THIRDQUARTER 2016)》 * |
李达: "《无线声学传感器网络中分布式语音增强方法研究》", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110739004A (zh) * | 2019-10-25 | 2020-01-31 | 大连理工大学 | 一种用于wasn的分布式语音噪声消除系统 |
CN110739004B (zh) * | 2019-10-25 | 2021-12-03 | 大连理工大学 | 一种用于wasn的分布式语音噪声消除系统 |
WO2021138201A1 (en) * | 2019-12-30 | 2021-07-08 | Texas Instruments Incorporated | Background noise estimation and voice activity detection system |
US11270720B2 (en) | 2019-12-30 | 2022-03-08 | Texas Instruments Incorporated | Background noise estimation and voice activity detection system |
CN111696575A (zh) * | 2020-06-19 | 2020-09-22 | 杭州电子科技大学 | 基于混合神经网络模型的低通气和呼吸暂停检测识别系统 |
CN111696575B (zh) * | 2020-06-19 | 2023-01-03 | 杭州电子科技大学 | 基于混合神经网络模型的低通气和呼吸暂停检测识别系统 |
CN114449410A (zh) * | 2021-11-12 | 2022-05-06 | 山东大学 | 一种多通道声纹信号同步采集系统及方法 |
CN114609493A (zh) * | 2022-05-09 | 2022-06-10 | 杭州兆华电子股份有限公司 | 一种信号数据增强的局部放电信号识别方法 |
CN114609493B (zh) * | 2022-05-09 | 2022-08-12 | 杭州兆华电子股份有限公司 | 一种信号数据增强的局部放电信号识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110289011B (zh) | 2021-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110289011B (zh) | 一种用于分布式无线声传感器网络的语音增强系统 | |
CN107479030B (zh) | 基于分频和改进的广义互相关双耳时延估计方法 | |
CN108172231B (zh) | 一种基于卡尔曼滤波的去混响方法及系统 | |
Yoshioka et al. | Generalization of multi-channel linear prediction methods for blind MIMO impulse response shortening | |
KR20180069879A (ko) | 음성 강화를 위해 전역적으로 최적화된 최소 제곱 포스트 필터링 | |
CN108447498B (zh) | 应用于麦克风阵列的语音增强方法 | |
CN106887239A (zh) | 用于高度相关的混合物的增强型盲源分离算法 | |
Furnon et al. | DNN-based mask estimation for distributed speech enhancement in spatially unconstrained microphone arrays | |
CN111048061B (zh) | 回声消除滤波器的步长获取方法、装置及设备 | |
O'Connor et al. | Diffusion-based distributed MVDR beamformer | |
CN115424627A (zh) | 基于卷积循环网络和wpe算法的语音增强混合处理方法 | |
CN110739004B (zh) | 一种用于wasn的分布式语音噪声消除系统 | |
Lv et al. | A permutation algorithm based on dynamic time warping in speech frequency-domain blind source separation | |
Song et al. | An integrated multi-channel approach for joint noise reduction and dereverberation | |
Asaei et al. | Binary sparse coding of convolutive mixtures for sound localization and separation via spatialization | |
CN114242104A (zh) | 语音降噪的方法、装置、设备及存储介质 | |
CN116403594B (zh) | 基于噪声更新因子的语音增强方法和装置 | |
CN112802490A (zh) | 一种基于传声器阵列的波束形成方法和装置 | |
Yoshioka et al. | Dereverberation by using time-variant nature of speech production system | |
CN116052702A (zh) | 一种基于卡尔曼滤波的低复杂度多通道去混响降噪方法 | |
Wang et al. | Two-stage enhancement of noisy and reverberant microphone array speech for automatic speech recognition systems trained with only clean speech | |
Chang et al. | Distributed parameterized topology-independent noise reduction in acoustic sensor networks | |
CN116106826A (zh) | 声源定位方法、相关装置和介质 | |
CN113591537B (zh) | 一种双迭代非正交联合块对角化卷积盲源分离方法 | |
CN113763984A (zh) | 一种用于分布式多说话人的参数化噪声消除系统 |
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 |