CN109545188A - 一种实时语音端点检测方法及装置 - Google Patents

一种实时语音端点检测方法及装置 Download PDF

Info

Publication number
CN109545188A
CN109545188A CN201811491292.8A CN201811491292A CN109545188A CN 109545188 A CN109545188 A CN 109545188A CN 201811491292 A CN201811491292 A CN 201811491292A CN 109545188 A CN109545188 A CN 109545188A
Authority
CN
China
Prior art keywords
frame
entropy
spectrum
threshold
value
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
CN201811491292.8A
Other languages
English (en)
Other versions
CN109545188B (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.)
Shenzhen Youjie Zhixin Technology Co Ltd
Original Assignee
Shenzhen Youjie Zhixin Technology Co Ltd
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 Shenzhen Youjie Zhixin Technology Co Ltd filed Critical Shenzhen Youjie Zhixin Technology Co Ltd
Priority to CN201811491292.8A priority Critical patent/CN109545188B/zh
Publication of CN109545188A publication Critical patent/CN109545188A/zh
Application granted granted Critical
Publication of CN109545188B publication Critical patent/CN109545188B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L15/00Speech recognition
    • G10L15/04Segmentation; Word boundary detection
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • G10L25/21Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being power information
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/48Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use
    • G10L25/51Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use for comparison or discrimination
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/78Detection of presence or absence of voice signals

Landscapes

  • Engineering & Computer Science (AREA)
  • Computational Linguistics (AREA)
  • Health & Medical Sciences (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

本发明涉及语音技术领域,具体地说,涉及一种实时语音端点检测方法及装置,方法包括以下步骤:信号分帧、加重;去脉冲处理;去除直流成分;计算每帧信号的短时能量和过零率;加窗处理;减谱处理;计算谱熵;计算变换平滑谱熵;语音帧和噪声帧初步判断;变换平滑谱熵与阈值的处理;语音段起始帧、结束帧判断;本发明根据信号是根据何种条件判断出来以及判断的结果,加权更新减谱阈值、变换平滑谱熵、相应的短时能量、相应的短时平均能量、减谱功率谱等参数的阈值,以使各阈值越来越准确,最终判断出的语音起始帧、结束帧也越来越准确;本发明能高效、准确地实时检测语音。

Description

一种实时语音端点检测方法及装置
技术领域
本发明涉及语音技术领域,具体地说,涉及一种实时语音端点检测方法及装置。
背景技术
语音活性检测(VAD),即在信号中检测出语音段,也即是语音端点检测技术。端点检测在语音信号处理领域中一直有着重要的意义。作为语音识别的前端,准确的端点检测可以提高识别的准确率;用于语音增强系统中,可以进行准确的噪声模型估计;在语音编码领域中可以降低编码的平均比特率并降低功耗。
目前,端点检测大体上可以分为两类:基于模型的检测方法和基于特征的检测方法。
基于模型的方法是以数据统计分析为基础建立一个能够更好刻画语音信号内部联系的模型,常见的有基于隐马尔科夫模型(HMM)、矢量量化(vector quantization,VQ)和支持向量机(support vector machine,SVM)等端点检测方法。该方法过程比较复杂、运算量大,并且实际应用环境复杂多变,噪声多种多样,建立的语音和噪声模型对环境的适应性比较差,检测准确率达不到实际应用的需求。比如:CN20141077954中公开的语音端点检测方法和装置、CN201510587721中公开的基于统计模型的语音端点检测方法及装置、CN201610192489中公开的语音端点检测方法及装置、CN201710076757中公开的一种语音端点检测方法及语音识别方法、CN201610886934中公开的一种语音激活检测方法及装置,这些检测方法都需要大量的训练。
基于特征的方法,一般都是在较高的信噪比下均能给出较高的检测率,但是在低信噪比下却不够理想。传统的能量和过零率特征在低信噪比下已不再适用。许多新的特征被提出,比如:倒谱、线性预测编码系数、基频及谐波特性、频带方差、频域能量、差分能量、幅度差、差分过零率、高阶统计量特征、高频能量和低频能量的特征等等。以上的多种特征在低信噪比噪下检测准确度仍然不够理想。比如:
1、基于单一参数的:CN200710179342中公开的一种孤立词语音端点检测的方法及系统,基于平均能量,对于强噪,能量法无法区分;CN201110071269中公开的一种语音识别的端点检测方法,基于线性预测编码系数,判断参数单一,需要背景噪声和语音的模板,对于变化的噪声和语音难以检测;
2、基于少量一些参数组合,但是参数并不是很好区分噪声和语音,或者不是很好计算的:CN200410083807中公开的基于滑动窗口的端点检测方法、装置和语音识别系统,主要是基于能量和信噪比,属于比较粗略的方法,信噪比难以估计准确,强噪基本都难以区分;CN200410090802中公开的一种应用于语音识别系统的语音端点检测方法,逐帧判断,频带划分不同信噪比的子带,再根据谐波特性来判断起始帧,同样,信噪比难于估计准确,判断参数较少;CN201410221983中公开的基于基频的端点检测系统及其处理方法,基于基频,以及基频的谐波位置辅助,基频位置确定容易出错,对于一些频谱成分及其丰富的强噪,比如强白噪,无法区分,判断参数还是比较少;
3、基于信息熵、谱熵这种容易计算且能很好区分噪声和语音的参数,并且和其他参数结合,参数较多的:CN201410292519中公开的一种利用短时时频值的自适应端点检测方法,基于短时能量、短时信息熵和短时幅度相对值,信息熵是逐频点处理,误差较大,计算量大;CN201710086400中公开的一种低信噪比环境下基于谱熵改进的语音端点检测方法,主要基于子带谱熵和能量比值,但是子带不是自适应划分,未加异常子带处理,并且未加平滑等处理。
所以需要一种有效的方法能对实时语音进行准确的检测。
发明内容
本发明的内容是提供一种实时语音端点检测方法及装置,其能够克服现有技术的某种或某些缺陷。
根据本发明的一种实时语音端点检测方法,其包括以下步骤:
一、信号分帧、加重;
二、去脉冲处理;
三、去除直流成分;
四、计算每帧信号的短时能量和过零率;
五、加窗处理;
六、减谱处理;
七、计算谱熵,将功率谱的一半频带,取一个有效频带,将有效频带的功率谱按照设定数目的频点作为一个子带,平均分成若干个子带,设每个子带频点数为Np,则子带个数为:
floor(Nvald/Np);
其中,floor函数表示向下取整,Nvald/Np表示有效频带频点数除以子带频点数;设结果是Nb个;
计算每个子带的功率值占总功率的概率,设为p;
根据这些子带的功率谱概率,自适应重新划分子带,调整子带的频点个数和频带个数,得到自适应划分子带的谱熵值;谱熵的计算公式为;
其中,SP为谱熵,m为子带个数,p(i)为第i个子带概率;
八、计算变换平滑谱熵,先做一定的平滑处理,平滑处理后,平滑谱熵取相反数后平移,得到一个新的参数,将取值区间映射为比较合适的值,并且使语音帧的参数值比较大,噪声帧的参数值比较小,新的参数为变换平滑谱熵;设平滑谱熵为SPm,平移值设为n,变换平滑谱熵为SPt,则计算公式为:
SPt=-SPm+n;
九、语音帧和噪声帧初步判断:根据当前帧的短时能量是否超过设置的语音帧能量阈值、以及其过零率是否在设置阈值范围内,若是,初步判断当前帧是语音帧,若否,则是噪声帧;
A、若是语音帧,则再分析判断其谐波特性以及高频段的功率谱特性:
1)谐波特性分析:
a、考查当前帧的功率谱,在前述有效频带内搜寻前8个有效峰值位置;
b、根据这些峰值的位置,即对应频点序号,判断这些位置之间是否存在近似谐波倍数关系,第一个峰值位置是基频位置,判断后面每个位置是否对其有近似2倍、3倍、4倍、5倍、6倍、7倍和8倍关系,近似关系可以通过设置一个浮动区间实现,如果有则计有一个谐波,计算出这一帧信号的谐波个数;
c、对谐波个数设置阈值,当谐波个数超过阈值即满足要求,继续分析高频段特性;
2)高频段特性分析:
取当前帧的一定频点以后的功率谱,求其峰值,和基频信号功率值的一定比例做比较,若峰值大,则认为此帧不是语音帧,否则认为是语音帧;
3)当前帧通过上述谐波特性和高频段功率谱特性判断后,如果是语音帧,则用其变换平滑谱熵来更新当前变换平滑谱熵阈值,更新方式是通过加权因子处理,并用当前帧的能量更新当前语音帧能量阈值,同样是加权因子方式得到,利用加权因子更新当前阈值的表达式如下:
设更新因子为a,当前阈值为Th,用来更新的变量值为x,则更新后的阈值为:
Th=(1-a)*Th+a*x;
如果不是语音帧,则也通过当前帧的变换平滑谱熵更新当前变换平滑谱熵阈值,采用不同的加权因子;同时也用当前帧的能量更新当前噪声帧能量阈值,同样采用加权因子更新;
B、若是噪声帧,则同样采用不同的加权因子更新当前变换平滑谱熵阈值,以及更新当前减谱去噪的功率谱阈值、当前噪声帧能量阈值;
十、计算当前帧的变换平滑谱熵和调整当前变换平滑谱熵阈值、当前语音帧能量阈值以及当前噪声帧能量阈值后,根据当前帧的变换平滑谱熵和当前变换平滑谱熵阈值之间的关系,对变换平滑谱熵阈值进行处理,即,若当前帧变换平滑谱熵值大于一定值,且减去当前变换平滑谱熵阈值的差值大于一定值,则判定当前变换平滑谱熵阈值过小,采用当前帧变换平滑谱熵值加权更新当前变换平滑谱熵阈值;反之,若当前帧变换平滑谱熵值小于一定值,且当前变换平滑谱熵阈值减去当前帧变换平滑谱熵值的差值大于一定值,则判定当前变换平滑谱熵阈值过大,采用当前帧变换平滑谱熵值加权更新当前变换平滑谱熵阈值;
十一、语音段起始帧、结束帧判断:
1)以帧号为顺序,设置搜寻标志,交替搜寻起始帧和结束帧;
2)判断起始帧:
a、初步判定语音段起始帧:以当前帧和之前一定数目的帧为一段,将这一段的变换平滑谱熵和当前变换平滑谱熵阈值比较,若大于阈值的帧数达到一定数目,则初步判定当前帧为语音段起始帧,进行步骤b;若没达到一定数目,从下一帧信号继续判断;
b、第二次判断:对当前帧的其他特性进行判断,所述其他特性为谱熵、短时能量、过零率、谐波特性和特定噪声特性;若当前帧不满足所述其他特性阈值要求,则判定当前帧为不是语音起始帧,从下一帧信号继续判断;若当前帧满足所述其他特性阈值要求,则进行步骤c;
c、第三次判断:在当前帧及向前推一定帧的帧段内,统计谱熵值小于阈值的帧数、具备谐波特性的帧数、过零率在阈值范围内的帧数和具备特定噪声特性的帧数,对所有这些帧数设定一些阈值,当任意一个帧数不满足这些阈值时,或者当上述帧段内平均能量均值小于起始帧段平均能量阈值,则排除当前帧是语音起始帧,即判断为非语音帧,从下一帧信号继续判断;反之则为语音段起始帧,并进行步骤d;
d、用当前帧的变换平滑谱熵值再次更新当前变换平滑谱熵阈值,用当前帧的能量再次更新当前语音帧能量阈值,用当前帧之前一定帧段内的平均能量更新当前起始帧段平均能量阈值;
3)若经过上述2)的判断,任何一步判断出当前帧不是语音帧,则采用当前帧的功率谱加权更新当前减谱阈值,用当前帧的能量值加权更新当前噪声帧能量阈值以及语音帧能量阈值;
4)判断结束帧:搜寻到起始帧后,设置间隔一定的帧数之后才开始搜寻结束帧,这个间隔为最小语音段长度;当检测出语音起始帧后,进入的帧数累计超过最小语音段长度的帧数,即开始判断结束帧,方法为:
考查当前帧向前推一定帧数的帧段,统计变换平滑谱熵值小于阈值的帧数、谱熵值大于阈值的帧数、无谐波特性的帧数、过零率不在阈值范围内的帧数和有特定噪声特性的帧数,对这些帧数设置阈值,若任意一个帧数超过阈值,或者这一帧段的平均能量小于结束帧段平均能量阈值,则判定当前帧为结束帧;将搜寻标志设置成搜寻起始帧标志,从下一帧起开始判断起始帧;
如果当前帧对于上个起始帧已经过去一定帧数,可设置一个较长帧数,在这个较长的帧段内,将每一帧的变换平滑谱熵值、短时能量值、每一帧的对应的变换平滑谱熵阈值和语音帧能量阈值都保存,则在这个较长的时间段内进行判断,若变换平滑谱熵值比阈值小的帧数太多或者太少,则认为是谱熵阈值自适应调整方式已经不适合,这时应给出结束帧,将阈值参数设置成初始值;若这一帧段的谱熵和变换平滑谱熵值参数都符合阈值,但是平均能量异常,这时也给出结束帧,将阈值参数设置成初始值;若时间太久,可设置一个最长语音段,超过这个时间直接给出结束帧;
5)若经过4)对当前帧进行结束帧判断,判断出不是结束帧,即当前帧仍是语音帧,则根据当前帧的变换平滑谱熵值加权更新当前变换平滑谱熵阈值,利用当前帧的能量值更新当前语音帧能量阈值,以及利用当前帧之前一定帧段的平均能量更新结束帧段平均能量阈值;若判断出是结束帧,则可以根据不同的情况,加权更新变换平滑谱熵阈值、噪声帧能量阈值、减谱功率谱阈值以及结束帧段平均能量阈值。
作为优选,去脉冲处理的方法如下:对于一个信号点,考查其前后一段区间内信号点的幅值均值,若此信号点幅值绝对值超过上述幅值均值一定比例,则认为当前信号点幅值是瞬间过大或者过小,有窄脉冲的特性,此时根据上述幅值将当前信号点幅值调小,幅值符号不变;若当前信号点幅值绝对值不超过上述均值的一定比例,则当前信号点不做处理。
作为优选,加窗处理中,可以是汉明窗和汉宁窗的一种,在频谱变换前对一帧信号做一个按信号点的幅值加权处理,以使频谱平滑。
作为优选,减谱处理的方法如下:设置开始检测的前一定帧段当成噪声帧,用于估计初始底噪,这个帧段可以是从打开检测功能即开始处理,也可以是之后一定帧再处理,以消除打开检测功能的按键开关以及参数初始化处理带来的干扰,使谱熵计算更准确,这里的减谱,只是为了后面计算谱熵更准确,并不真正对原始语音信号做减谱去噪,也就是减谱后信号不需要恢复到时域,这样就减小了计算量;可取前10帧,作为噪声段,得出平均功率,作为减谱初值;在10帧前,直接计算每帧的功率谱谱熵,10帧后,先减谱处理再计算谱熵。
作为优选,自适应重新划分子带的方法为:
1)若某个子带功率谱概率大于0.9,则在有效频带中去掉这个子带的所有频点,即有效频带点数变少,重新计算各个子带的功率谱概率;若没有子带功率谱概率大于0.9,则直接执行2);
2)计算最小子带概率的熵,用Nmin表示,设:
Nmin=-pmin*lg(pmin),设pmin是最小子带功率谱概率;
根据Nmin取值,重新计算子带个数,设重新计算的子带个数为V,采用以下公式:
得出新的子带个数,其中第二个表达式最终得出结果向下取整数,然后,根据有效频带的频点数对每个子带均匀分频点数,设有效频带点数为N'vald,则重新分到每个子带的频点数为:
floor(N'vald/V);
floor函数表示向下取整,/表示做除法;上述过程即是自适应重新划分子带;
再次计算个子带的概率,最终得出这一帧信号的谱熵值,即是自适应划分子带的谱熵值。
作为优选,步骤八中的平滑处理的方法为:以当前帧向前算起连续5帧,设谱熵值依次为S(i-4),S(i-3),S(i-2),S(i-1),S(i),即S(i)为当前帧的谱熵,平滑的计算为:
当i小于5时,即第5帧以前,不做平滑处理。
作为优选,交替搜寻起始帧和结束帧的方法为:初始状态为搜寻起始帧标志,当搜寻到起始帧时,将标志设置为搜寻结束帧标志,最小语音段之后开始搜索结束帧,当搜寻到结束帧时,将搜寻标志设置为搜寻起始帧标志,从下一帧起搜寻起始帧,如此交替搜寻和设置。
作为优选,特定噪声特性为汽车发动机、汽车喇叭、警报、敲击桌子、文具、餐具、键盘、空调和电扇的噪声特性的一种或多种;其他,特定噪声特性判断,是根据这些噪声本身的时频特性通过机器录制大量这些噪声来分析总结出来的。比如汽车发动机噪声,典型的是在小轿车旁边或者里面采集到的比较平稳的噪声,其频谱峰值在比较低的频点,谱峰很窄,且越高频越衰减,无谐波特性;汽车喇叭,属于宽频噪声,在频谱高频段依然有很高的峰值;警报信号,是指火警,是一个随时间变化的区间扫频信号、频谱主峰值在语音范围但是没有谐波特性;敲击桌子、文具、键盘等,类似于单频或者单频带信号,能根据若信号频谱太过于集中来加以判别;空调、电扇等强噪,和汽车发动机类似,频谱主要集中在一定频段,且主峰值对应频点比汽车发动机噪声更低。这些噪声和特征可以根据应用需要和方法的有效性来扩展或者减少,这里最重要的是提供了一种思路,并且在最后的产品中也考虑了我们产品实际应用的一些场景噪声,显示了良好的适应效果。
作为优选,当判断出语音起始帧时,加入信号延时处理。若不加延时处理,则在语音段起始帧为一些清音段时,容易漏判,若采用延时处理,当判断出语音起始帧时,将此帧之前一定数目的帧当成起始帧,即当前帧的判断是在一定帧之后才做出的,这样能非常有效地检测出语音起始段的清音段。此外,对于一些语音段之间偶尔有一些间隙被判断成非语音帧的情况,由于采用了延时处理,可以在出现非语音段帧数太少时直接将这段非语音帧改成语音帧,即不会漏掉语音段的间隙。
本发明基于变换平滑自适应子带功率谱熵,通过特定变换得到变换平滑谱熵参数作为主参数来判断语音和非语音信号,不受信号能量影响,在强噪声环境下能够很好的防止误判;计算谱熵是按照划分子带来处理,不是逐点处理,减小计算量以及误差;划分子带是自适应的,更加提高准确性和实用性;由于语音信号短时平稳性,对谱熵参数加平滑处理,消除不正常的随机干扰,并且,由于直接计算出的谱熵参数比较小,并且是宽频噪声的谱熵值比较大,语音或者含噪语音信号的比较小,对谱熵参数取相反数,并且做一个线性平移,得到一个新的参数,参数取值区间比较合理,且语音或含噪语音信号的值比较大,噪声的值比较小,这样比较符合判断习惯,并且阈值比较好设置和处理;将变换平滑谱熵作为判断的主参数,对其设置阈值,提高检测的准确性和适应性;辅助以语音信号帧能量阈值、噪声信号帧能量阈值、语音信号帧段平均能量阈值、噪声信号帧段平均能量阈值、过零率阈值以及谐波特性分析、高频段特性以及多种典型噪声的特性分析等特征判断,能够使检测结果更具鲁棒性,防止很多类型的强噪误判;多种典型噪声是指产品应用环境的典型噪声,比如敲击桌子、文具、餐具、键盘等噪声,汽车发动机、汽车喇叭噪、警报信号噪声,空调、电扇等噪声,专门用产品录制大量这些噪声,试验时若基本方法不能区分,则针对性分析这些噪声的时域、频域特征,确定排除方法;采用单帧和帧段结合判断,并加了一些符合实际情况和应用习惯的处理和异常处理,比如可设置最小语音段,最长语音段等;变换平滑谱熵、语音帧能量、噪声帧能量、语音段平均能量、减谱功率谱等阈值都是动态更新的;作为可选,采用一点点延时处理,这基本不影响算法实时性。具体有,为排除语音起始段的清音漏检,处理时采用延时存取一定帧数,判断出起始帧后向前推一定帧数当成起始帧,试验效果非常好,以及,对于一些语音段之间偶尔有一些间隙被判断成非语音帧的情况,由于采用了延时处理,可以在出现非语音段帧数太少时直接将这段非语音帧改成语音帧,即不会漏掉语音段的合理间隙,增强语音信号的连续性;加入一些去脉冲滤波、减谱等处理,还可以扩展更多前处理,使语音检测更加准确,其中减谱阈值也是自适应更新的,能够适应不同类型的非平稳背景噪声;程序空间小,对于一些需要保存当前帧以及之前一定帧段内的参数,可以设置循环移位寄存器,计算量少,远没有模型法复杂。可移植性好,已经在TI的浮点DSP平台上的产品中应用,初步优化后判断一帧20ms、采样率8k、量化位数16位数据的时间最大在4ms左右,完全满足实时处理要求。
本发明根据信号是根据何种条件判断出来以及判断的结果,加权更新谱熵、变换平滑谱熵、相应的短时能量、相应的短时平均能量等参数的阈值,以使各阈值越来越准确,最终判断出的语音起始帧、结束帧也越来越准确。
本发明还提供了一种实时语音端点检测装置,其包括:
信号分帧加重模块,用于接收信号数据,对其进行分帧、加重处理;
去脉冲处理模块,与信号分帧加重模块连接,用于消除一些脉冲噪声;
去除直流模块,与去脉冲处理模块连接,用于去除直流成分;
短时能量和过零率计算模块,与去直流模块连接,用于计算每一帧的短时能量和过零率;
加窗处理模块,与短时能量和过零率计算模块连接,用于在频谱变换前对一帧信号做一个按信号点的幅值加权处理,以使频谱平滑;
减谱处理模块,与加窗处理模块连接,用于使谱熵计算更准确;
谱熵计算模块,与减谱处理模块连接,用于计算自适应划分子带的谱熵值;
变换平滑谱熵计算模块,与谱熵计算模块连接,用于平滑、变换处理谱熵,得到变换平滑谱熵;
语音噪声帧初步判断模块,与变换平滑谱熵计算模块连接,用于初步判断当前帧是语音还是噪声,更新变换平滑谱熵阈值、语音帧和噪声帧能量帧阈值;
语音段处理模块,与语音噪声帧初步判断模块连接,用于起始帧、结束帧的判断;
阈值调整模块,与语音段处理模块连接,用于根据不同情况判断出每帧信号情况后,根据判断结果调整相应阈值。
本实时语音端点检测装置通过上述实时语音端点检测方法能有效准确的检测语音。
附图说明
图1为实施例1中一种实时语音端点检测装置的结构框图。
具体实施方式
为进一步了解本发明的内容,结合附图和实施例对本发明作详细描述。应当理解的是,实施例仅仅是对本发明进行解释而并非限定。
实施例1
本实施例提供了一种实时语音端点检测方法,其包括以下步骤:
一、信号分帧、加重;可以设置帧重叠,先做加重处理,加重方法属于现有技术;
二、去脉冲处理;消除一些脉冲噪声;
三、去除直流成分;即将所有信号点幅值减去这一帧信号所有点的幅值均值;
四、计算每帧信号的短时能量和过零率;计算方法属于现有技术;
五、加窗处理;加汉宁窗或汉明窗,属于现有技术;
六、减谱处理;减谱去噪属于现有技术,但本实施例中并不需要将频域减谱的信号恢复到时域中去,仅仅是为了使语音和噪声的谱熵计算更加准确;
七、计算谱熵,将功率谱的一半频带,取一个有效频带,通常是选择去掉0频点开始的一些低频点和奈奎斯特频点以下的一些高频点,将有效频带的功率谱按照设定数目的频点作为一个子带,平均分成若干个子带,设每个子带频点数为Np,则子带个数为:
floor(Nvald/Np);
其中,floor函数表示向下取整,Nvald/Np表示有效频带频点数除以子带频点数;设结果是Nb个;
计算每个子带的功率值占总功率的概率,设为p;
根据这些子带的功率谱概率,自适应重新划分子带,调整子带的频点个数和频带个数,得到自适应划分子带的谱熵值;谱熵的计算公式为;
其中,SP为谱熵,m为子带个数,p(i)为第i个子带概率;
八、计算变换平滑谱熵,先做一定的平滑处理,由于直接计算的谱熵参数,某些相邻帧之间由于有随机、突发噪声干扰可能造成抖动较大,并且,由于语音的短时平稳性,对相邻一定帧数的谱熵进行平滑处理,可以有效防止这些不正常的抖动,相邻帧的谱熵参数也更有连续性;平滑处理后,平滑谱熵取相反数后平移,得到一个新的参数,将取值区间映射为合适的值,并且使语音帧的参数值比较大,噪声帧的参数值比较小,新的参数为变换平滑谱熵;设平滑谱熵为SPm,平移值设为n,变换平滑谱熵为SPt,则计算公式为:
SPt=-SPm+n;
九、语音帧和噪声帧初步判断:根据当前帧的短时能量是否超过设置的语音帧能量阈值、以及其过零率是否在设置阈值范围内,若是,初步判断当前帧是语音帧,若否,则是噪声帧;
A、若是语音帧,则再分析判断其谐波特性以及高频段的功率谱特性,这是基于语音特别是浊音有谐波特性,即存在基频及各倍数的谐波成分;至于高频段功率谱特性,是指语音帧功率谱一般在高频的幅值比较小,不像一些宽频噪声,在高频段还存在较多的高幅度峰值。
1)谐波特性分析:
a、考查当前帧的功率谱,在前述有效频带内搜寻前8个有效峰值位置;谐波特性中的有效峰值,不仅仅是峰值,还需要大于一定的幅值,这里的幅值取当前帧功率谱最大峰值的一定比例,比例参数根据实验设置;
b、根据这些峰值的位置,即对应频点序号,判断这些位置之间是否存在近似谐波倍数关系,第一个峰值位置是基频位置,判断后面每个位置是否对其有近似2倍、3倍、4倍、5倍、6倍、7倍和8倍关系,近似关系可以通过设置一个浮动区间实现,如果有则计有一个谐波,计算出这一帧信号的谐波个数;
c、对谐波个数设置阈值,当谐波个数超过阈值即满足要求,继续分析高频段特性;
2)高频段特性分析:
取当前帧的一定频点以后的功率谱,求其峰值,和基频信号功率值的一定比例做比较,若峰值大,则认为此帧不是语音帧,否则认为是语音帧;
3)当前帧通过上述谐波特性和高频段功率谱特性判断后,如果是语音帧,则用其变换平滑谱熵来更新当前变换平滑谱熵阈值,更新方式是通过加权因子处理,并用当前帧的能量更新当前语音帧能量阈值,同样是加权因子方式得到,利用加权因子更新当前阈值的表达式如下:
设更新因子为a,当前阈值为Th,用来更新的变量值为x,则更新后的阈值为:
Th=(1-a)*Th+a*x;
如果不是语音帧,则也通过当前帧的变换平滑谱熵更新当前变换平滑谱熵阈值,采用不同的加权因子;同时也用当前帧的能量更新当前噪声帧能量阈值,同样采用加权因子更新;
B、若是噪声帧,则同样采用不同的加权因子更新当前变换平滑谱熵阈值,以及更新当前减谱去噪的功率谱阈值、当前噪声帧能量阈值;
十、计算当前帧的变换平滑谱熵和调整当前变换平滑谱熵阈值、当前语音帧能量阈值以及当前噪声帧能量阈值后,根据当前帧的变换平滑谱熵和当前变换平滑谱熵阈值之间的关系,对变换平滑谱熵阈值进行处理,即,若当前帧变换平滑谱熵值大于一定值,且减去当前变换平滑谱熵阈值的差值大于一定值,则判定当前变换平滑谱熵阈值过小,采用当前帧变换平滑谱熵值加权更新当前变换平滑谱熵阈值;反之,若当前帧变换平滑谱熵值小于一定值,且当前变换平滑谱熵阈值减去当前帧变换平滑谱熵值的差值大于一定值,则判定当前变换平滑谱熵阈值过大,采用当前帧变换平滑谱熵值加权更新当前变换平滑谱熵阈值;
十一、语音段起始帧、结束帧判断:
1)以帧号为顺序,设置搜寻标志,交替搜寻起始帧和结束帧;
2)判断起始帧:
a、初步判定语音段起始帧:以当前帧和之前一定数目的帧为一段,将这一段的变换平滑谱熵和当前变换平滑谱熵阈值比较,若大于阈值的帧数达到一定数目,则初步判定当前帧为语音段起始帧,进行步骤b;若没达到一定数目,从下一帧信号继续判断;
b、第二次判断:对当前帧的其他特性进行判断,所述其他特性为谱熵、短时能量、过零率、谐波特性和特定噪声特性;若当前帧不满足所述其他特性阈值要求,则判定当前帧为不是语音起始帧,从下一帧信号继续判断;若当前帧满足所述其他特性阈值要求,则进行步骤c;
c、第三次判断:在当前帧及向前推一定帧的帧段内,统计谱熵值小于阈值的帧数、具备谐波特性的帧数、过零率在阈值范围内的帧数和具备特定噪声特性的帧数,对所有这些帧数设定一些阈值,当任意一个帧数不满足这些阈值时,或者当上述帧段内平均能量均值小于起始帧段平均能量阈值,则排除当前帧是语音起始帧,即判断为非语音帧,从下一帧信号继续判断;反之则为语音段起始帧,将搜寻标志设置为搜寻结束帧标志并进行步骤d;阈值判断示例为,若谱熵值大于谱熵阈值(谱熵值大,认为是噪声帧)的帧数大于设置阈值,则排除当前帧为语音帧;或者,若具有汽车发动机噪声特性的帧数大于设置阈值,则排除当前帧为语音帧,等等。
d、用当前帧的变换平滑谱熵值再次更新当前变换平滑谱熵阈值,用当前帧的能量再次更新当前语音帧能量阈值,用当前帧之前一定帧段内的平均能量更新当前起始帧段平均能量阈值;
3)若经过上述2)的判断,任何一步判断出当前帧不是语音帧,则采用当前帧的功率谱加权更新当前减谱阈值,用当前帧的能量值加权更新当前噪声帧能量阈值以及语音帧能量阈值;
4)判断结束帧:搜寻到起始帧后,设置间隔一定的帧数之后才开始搜寻结束帧,这个间隔为最小语音段长度;在语音起始帧检测不准时,能避免语音起始帧太快切换到语音结束帧,在语音起始帧检测准确时,又能避免后面连续一定帧不必要的语音检测,因为语音段都有一定长度,不止一两帧;当检测出语音起始帧后,进入的帧数累计超过最小语音段长度的帧数,即开始判断结束帧,方法为:
考查当前帧向前推一定帧数的帧段,统计变换平滑谱熵值小于阈值的帧数、谱熵值大于阈值的帧数、无谐波特性的帧数、过零率不在阈值范围内的帧数和有特定噪声特性的帧数,对这些帧数设置阈值,若任意一个帧数超过阈值,或者这一帧段的平均能量小于结束帧段平均能量阈值,则判定当前帧为结束帧;将搜寻标志设置成搜寻起始帧标志,从下一帧起开始判断起始帧;
如果当前帧对于上个起始帧已经过去一定帧数,可设置一个较长帧数,在这个较长的帧段内,将每一帧的变换平滑谱熵值、短时能量值、每一帧的对应的变换平滑谱熵阈值和语音帧能量阈值都保存,则在这个较长的时间段内进行判断,若变换平滑谱熵值比阈值小的帧数太多或者太少,则可能是谱熵阈值自适应调整方式已经不适合,这时应给出结束帧,将阈值参数设置成初始值;若这一帧段的谱熵和变换平滑谱熵值参数都符合阈值,但是平均能量异常,这时也给出结束帧,将阈值参数设置成初始值;若时间太久,可设置一个最长语音段,超过这个时间直接给出结束帧;这些相当于一些合理的异常处理;总之,这些可以防止语音段异常持续不能正常结束或者满足实际需求地去实现可以设置最大长度的语音段检测;
5)若经过4)对当前帧进行结束帧判断,判断出不是结束帧,即当前帧仍是语音帧,则根据当前帧的变换平滑谱熵值加权更新当前变换平滑谱熵阈值,利用当前帧的能量值更新当前语音帧能量阈值,以及利用当前帧之前一定帧段的平均能量更新结束帧段平均能量阈值;若判断出是结束帧,则可以根据不同的情况,加权更新变换平滑谱熵阈值、噪声帧能量阈值、减谱功率谱阈值以及结束帧段平均能量阈值。具体哪些方式采用怎样的加权因子,可以通过大量试验确定。
之所以判断出结束帧之后下一帧就开始判断起始帧,是因为有些结束帧是误判,此时若下一帧或者几帧能马上判断出起始帧,则这个误判带来的一帧或几帧损失,是没多大影响的,所以这样处理能尽可能减小结束帧误判造成的语音帧缺失。
本实施例中,去脉冲处理的方法如下:对于一个信号点,考查其前后一段区间内信号点的幅值均值,若此信号点幅值绝对值超过上述幅值均值一定比例,则认为当前信号点幅值是瞬间过大或者过小,有窄脉冲的特性,此时根据上述幅值将当前信号点幅值调小,幅值符号不变;若当前信号点幅值绝对值不超过上述均值的一定比例,则当前信号点不做处理。
本实施例中,加窗处理中,可以是汉明窗和汉宁窗的一种,在频谱变换前对一帧信号做一个按信号点的幅值加权处理,以使频谱平滑。
本实施例中,减谱处理的方法如下:设置开始检测的前一定帧段当成噪声帧,用于估计初始底噪,这个帧段可以是从打开检测功能即开始处理,也可以是之后一定帧再处理,以消除打开检测功能的按键开关以及参数初始化处理带来的干扰,使谱熵计算更准确;本实施例取前10帧,作为噪声段,得出平均功率,作为减谱初值;在10帧前,直接计算每帧的功率谱谱熵,10帧后,先减谱处理再计算谱熵。每一帧信号我们做512点的FFT,功率谱取前257点。这里的减谱,只是为了后面计算谱熵更准确,并不真正对原始语音信号做减谱去噪,也就是减谱后信号不需要恢复到时域,这样就减小了计算量。
本实施例中,自适应重新划分子带的方法为:
1)若某个子带功率谱概率大于0.9,则在有效频带中去掉这个子带的所有频点,即有效频带点数变少,重新计算各个子带的功率谱概率;若没有子带功率谱概率大于0.9,则直接执行2);
2)计算最小子带概率的熵,用Nmin表示,设:
Nmin=-pmin*lg(pmin),设pmin是最小子带功率谱概率;
根据Nmin取值,重新计算子带个数,设重新计算的子带个数为V,采用以下公式:
得出新的子带个数,其中第二个表达式最终得出结果向下取整数,然后,根据有效频带的频点数对每个子带均匀分频点数,设有效频带点数为N'vald,则重新分到每个子带的频点数为:
floor(N'vald/V);
floor函数表示向下取整,/表示做除法;上述过程即是自适应重新划分子带;
再次计算个子带的概率,最终得出这一帧信号的谱熵值,即是自适应划分子带的谱熵值。
本实施例中,步骤八中的平滑处理的方法为:以当前帧向前算起连续5帧,设谱熵值依次为S(i-4),S(i-3),S(i-2),S(i-1),S(i),即S(i)为当前帧的谱熵,平滑的计算为:
当i小于5时,即第5帧以前,不做平滑处理。
本实施例中,交替搜寻起始帧和结束帧的方法为:初始状态为搜寻起始帧标志,当搜寻到起始帧时,将标志设置为搜寻结束帧标志,最小语音段之后开始搜索结束帧,当搜寻到结束帧时,将搜寻标志设置为搜寻起始帧标志,从下一帧起搜寻起始帧,如此交替搜寻和设置。
本实施例中,特定噪声特性为汽车发动机、汽车喇叭、警报、敲击桌子、文具、餐具、键盘、空调和电扇的噪声特性的一种或多种;特定噪声特性判断,是通过产品录制大量这些噪声来分析其时频特性总结出来的。比如汽车发动机噪声,典型的是在小轿车旁边或者里面采集到的比较平稳的噪声,其频谱峰值在比较低的频点,谱峰很窄,且越高频越衰减,无谐波特性;汽车喇叭,属于宽频噪声,在频谱高频段依然有很高的峰值;警报信号,是指火警,是一个随时间变化的区间扫频信号、频谱主峰值在语音范围但是没有谐波特性;敲击桌子、文具、键盘等,类似于单频或者单频带信号,能根据若信号频谱太过于集中来加以判别;空调、电扇等强噪,和汽车发动机类似,频谱主要集中在一定频段,且主峰值对应频点比汽车发动机噪声更低。这些噪声和特征可以根据应用需要和方法的有效性来扩展或者减少,这里最重要的是提供了一种思路,并且在最后的产品中也考虑了我们产品实际应用的一些场景噪声,显示了良好的适应效果。
本实施例中,当判断出语音起始帧时,加入信号延时处理。若不加延时处理,则在语音段起始帧为一些清音段时,容易漏判,若采用延时处理,当判断出语音起始帧时,将此帧之前一定数目的帧当成起始帧,即当前帧的判断是在一定帧之后才做出的,这样能非常有效地检测出语音起始段的清音段。此外,对于一些语音段之间偶尔有一些间隙被判断成非语音帧的情况,由于采用了延时处理,可以在出现非语音段帧数太少时直接将这段非语音帧改成语音帧,即不会漏掉语音段的间隙。
本发明基于变换平滑自适应子带功率谱熵,通过特定变换得到变换平滑谱熵参数作为主参数来判断语音和非语音信号,不受信号能量影响,在强噪声环境下能够很好的防止误判;计算谱熵是按照划分子带来处理,不是逐点处理,减小计算量以及误差;划分子带是自适应的,更加提高准确性和实用性;由于语音信号短时平稳性,对谱熵参数加平滑处理,消除不正常的随机干扰,并且,由于直接计算出的谱熵参数比较小,并且是宽频噪声的谱熵值比较大,语音或者含噪语音信号的比较小,对谱熵参数取相反数,并且做一个线性平移,得到一个新的参数,参数取值区间比较合理,且语音或含噪语音信号的值比较大,噪声的值比较小,这样比较符合判断习惯,并且阈值比较好设置和处理;将变换平滑谱熵作为判断的主参数,对其设置阈值,提高检测的准确性和适应性;辅助以语音信号帧能量阈值、噪声信号帧能量阈值、语音信号帧段平均能量阈值、过零率阈值以及谐波特性分析、高频段特性以及多种典型噪声的特性分析等特征判断,能够使检测结果更具鲁棒性,防止很多类型的强噪误判;多种典型噪声是指产品应用环境的典型噪声,比如敲击桌子、文具、餐具、键盘等噪声,汽车发动机、汽车喇叭噪、警报信号噪声,空调、电扇等噪声,专门用产品录制大量这些噪声,试验时若基本方法不能区分,则针对性分析这些噪声的时域、频域特征,确定排除方法;采用单帧和帧段结合判断,并加了一些符合实际情况和应用习惯的处理和异常处理,比如可设置最小语音段,最长语音段等;变换平滑谱熵、语音帧能量、噪声帧能量、语音段平均能量、减谱功率谱等阈值都是动态更新的;作为可选,采用一点点延时处理,这基本不影响算法实时性。具体有,为排除语音起始段的清音漏检,处理时采用延时存取一定帧数,典型如5至10帧,判断出起始帧后向前推一定帧数当成起始帧,试验效果非常好,以及,对于一些语音段之间偶尔有一些间隙被判断成非语音帧的情况,由于采用了延时处理,可以在出现非语音段帧数太少时直接将这段非语音帧改成语音帧,即不会漏掉语音段的合理间隙,增强语音信号的连续性;加入一些去脉冲滤波、减谱等处理,还可以扩展更多前处理,使语音检测更加准确,其中减谱阈值也是自适应更新的,能够适应不同类型的非平稳背景噪声;程序空间小,对于一些需要保存当前帧以及之前一定帧段内的参数,可以设置循环移位寄存器,计算量少,远没有模型法复杂。可移植性好,已经在TI的浮点DSP平台上的产品中应用,初步优化后判断一帧20ms、采样率8k、量化位数16位的数据的时间最大在4ms左右,完全满足实时处理要求。
如图1所示,本实施例还提供了一种实时语音端点检测装置,其包括:
信号分帧加重模块,用于接收信号数据,对其进行分帧、加重处理;
去脉冲处理模块,与信号分帧加重模块连接,用于消除一些脉冲噪声;
去除直流模块,与去脉冲处理模块连接,用于去除直流成分;
短时能量和过零率计算模块,与去直流模块连接,用于计算每一帧的短时能量和过零率;
加窗处理模块,与短时能量和过零率计算模块连接,用于在频谱变换前对一帧信号做一个按信号点的幅值加权处理,以使频谱平滑;
减谱处理模块,与加窗处理模块连接,用于使谱熵计算更准确;
谱熵计算模块,与减谱处理模块连接,用于计算自适应划分子带的谱熵值;
变换平滑谱熵计算模块,与谱熵计算模块连接,用于平滑、变换处理谱熵,得到变换平滑谱熵;
语音噪声帧初步判断模块,与变换平滑谱熵计算模块连接,用于初步判断当前帧是语音还是噪声,更新变换平滑谱熵阈值、语音帧和噪声帧能量帧阈值;
语音段处理模块,与语音噪声帧初步判断模块连接,用于起始帧、结束帧的判断;
阈值调整模块,与语音段处理模块连接,用于根据不同情况判断出每帧信号情况后,根据判断结果调整相应阈值。
本实时语音端点检测装置通过上述实时语音端点检测方法能有效准确的检测语音。
以上示意性的对本发明及其实施方式进行了描述,该描述没有限制性,附图中所示的也只是本发明的实施方式之一,实际的结构并不局限于此。所以,如果本领域的普通技术人员受其启示,在不脱离本发明创造宗旨的情况下,不经创造性的设计出与该技术方案相似的结构方式及实施例,均应属于本发明的保护范围。

Claims (10)

1.一种实时语音端点检测方法,其特征在于:包括以下步骤:
一、信号分帧、加重;
二、去脉冲处理;
三、去除直流成分;
四、计算每帧信号的短时能量和过零率;
五、加窗处理;
六、减谱处理;
七、计算谱熵,将功率谱的一半频带,取一个有效频带,将有效频带的功率谱按照设定数目的频点作为一个子带,平均分成若干个子带,设每个子带频点数为Np,则子带个数为:
floor(Nvald/Np);
其中,floor函数表示向下取整,Nvald/Np表示有效频带频点数除以子带频点数;设结果是Nb个;
计算每个子带的功率值占总功率的概率,设为p
根据这些子带的功率谱概率,自适应重新划分子带,调整子带的频点个数和频带个数,得到自适应划分子带的谱熵值;谱熵的计算公式为;
其中,SP为谱熵,m为子带个数,p(i)为第i个子带概率;
八、计算变换平滑谱熵,先做一定的平滑处理,平滑处理后,平滑谱熵取相反数后平移,得到一个新的参数,将取值区间映射为合适的值,并且使语音帧的参数值比较大,噪声帧的参数值比较小,新的参数为变换平滑谱熵;设平滑谱熵为SPm,平移值设为n,变换平滑谱熵为SPt,则计算公式为:
SPt=-SPm+n;
九、语音帧和噪声帧初步判断:根据当前帧的短时能量是否超过设置的语音帧能量阈值以及其过零率是否在设置阈值范围内,若是,初步判断当前帧是语音帧,若否,则是噪声帧;
A、若是语音帧,则再分析判断其谐波特性以及高频段的功率谱特性:
1)谐波特性分析:
a、考查当前帧的功率谱,在前述有效频带内搜寻前8个有效峰值位置;
b、根据这些峰值的位置,即对应频点序号,判断这些位置之间是否存在近似谐波倍数关系,第一个峰值位置是基频位置,判断后面每个位置是否对其有近似2倍、3倍、4倍、5倍、6倍、7倍和8倍关系,近似关系可以通过设置一个浮动区间实现,如果有则计有一个谐波,计算出这一帧信号的谐波个数;
c、对谐波个数设置阈值,当谐波个数超过阈值即满足要求,继续分析高频段特性;
2)高频段特性分析:
取当前帧的一定频点以后的功率谱,求其峰值,和基频信号功率值的一定比例做比较,若峰值大,则认为此帧不是语音帧,否则认为是语音帧;
3)当前帧通过上述谐波特性和高频段功率谱特性判断后,如果是语音帧,则用其变换平滑谱熵来更新当前变换平滑谱熵阈值,更新方式是通过加权因子处理,并用当前帧的能量更新当前语音帧能量阈值,同样是加权因子方式得到,利用加权因子更新当前阈值的表达式如下:
设更新因子为a,当前阈值为Th,用来更新的变量值为x,则更新后的阈值为:
Th=(1-a)*Th+a*x;
如果不是语音帧,则也通过当前帧的变换平滑谱熵更新当前变换平滑谱熵阈值,采用不同的加权因子;同时也用当前帧的能量更新当前噪声帧能量阈值,同样采用加权因子更新;
B、若是噪声帧,则同样采用不同的加权因子更新当前变换平滑谱熵阈值,以及更新当前减谱去噪的功率谱阈值、当前噪声帧能量阈值;
十、计算当前帧的变换平滑谱熵和调整当前变换平滑谱熵阈值、当前语音帧能量阈值以及当前噪声帧能量阈值后,根据当前帧的变换平滑谱熵和当前变换平滑谱熵阈值之间的关系,对变换平滑谱熵阈值进行处理,即,若当前帧变换平滑谱熵值大于一定值,且减去当前变换平滑谱熵阈值的差值大于一定值,则判定当前变换平滑谱熵阈值过小,采用当前帧变换平滑谱熵值加权更新当前变换平滑谱熵阈值;反之,若当前帧变换平滑谱熵值小于一定值,且当前变换平滑谱熵阈值减去当前帧变换平滑谱熵值的差值大于一定值,则判定当前变换平滑谱熵阈值过大,采用当前帧变换平滑谱熵值加权更新当前变换平滑谱熵阈值;
十一、语音段起始帧、结束帧判断:
1)以帧号为顺序,设置搜寻标志,交替搜寻起始帧和结束帧;
2)判断起始帧:
a、初步判定语音段起始帧:以当前帧和之前一定数目的帧为一段,将这一段的变换平滑谱熵和当前变换平滑谱熵阈值比较,若大于阈值的帧数达到一定数目,则初步判定当前帧为语音段起始帧,进行步骤b;若没达到一定数目,从下一帧信号继续判断;
b、第二次判断:对当前帧的其他特性进行判断,所述其他特性为谱熵、短时能量、过零率、谐波特性和特定噪声特性;若当前帧不满足所述其他特性阈值要求,则判定当前帧为不是语音起始帧,从下一帧信号继续判断;若当前帧满足所述其他特性阈值要求,则进行步骤c;
c、第三次判断:在当前帧及向前推一定帧的帧段内,统计谱熵值小于阈值的帧数、具备谐波特性的帧数、过零率在阈值范围内的帧数和具备特定噪声特性的帧数,对所有这些帧数设定一些阈值,当任意一个帧数不满足这些阈值时,或者当上述帧段内平均能量均值小于起始帧段平均能量阈值,则排除当前帧是语音起始帧,即判断为非语音帧,从下一帧信号继续判断;反之则为语音段起始帧,并进行步骤d;
d、用当前帧的变换平滑谱熵值再次更新当前变换平滑谱熵阈值,用当前帧的能量再次更新当前语音帧能量阈值,用当前帧之前一定帧段内的平均能量更新当前起始帧段平均能量阈值;
3)若经过上述2)的判断,任何一步判断出当前帧不是语音帧,则采用当前帧的功率谱加权更新当前减谱阈值,用当前帧的能量值加权更新当前噪声帧能量阈值以及语音帧能量阈值;
4)判断结束帧:搜寻到起始帧后,设置间隔一定的帧数之后才开始搜寻结束帧,这个间隔为最小语音段长度;当检测出语音起始帧后,进入的帧数累计超过最小语音段长度的帧数,即开始判断结束帧,方法为:
考查当前帧向前推一定帧数的帧段,统计变换平滑谱熵值小于阈值的帧数、谱熵值大于阈值的帧数、无谐波特性的帧数、过零率不在阈值范围内的帧数和有特定噪声特性的帧数,对这些帧数设置阈值,若任意一个帧数超过阈值,或者这一帧段的平均能量小于结束帧段平均能量阈值,则判定当前帧为结束帧;将搜寻标志设置成搜寻起始帧标志,从下一帧起开始判断起始帧;
如果当前帧对于上个起始帧已经过去一定帧数,可设置一个较长帧数,在这个较长的帧段内,将每一帧的变换平滑谱熵值、短时能量值、每一帧的对应的变换平滑谱熵阈值和语音帧能量阈值都保存,则在这个较长的时间段内进行判断,若变换平滑谱熵值比阈值小的帧数太多或者太少,则认为是谱熵阈值自适应调整方式已经不适合,这时应给出结束帧,将阈值参数设置成初始值;若这一帧段的谱熵和变换平滑谱熵值参数都符合阈值,但是平均能量异常,这时也给出结束帧,将阈值参数设置成初始值;若时间太久,可设置一个最长语音段,超过这个时间直接给出结束帧;
5)若经过4)对当前帧进行结束帧判断,判断出不是结束帧,即当前帧仍是语音帧,则根据当前帧的变换平滑谱熵值加权更新当前变换平滑谱熵阈值,利用当前帧的能量值更新当前语音帧能量阈值,以及利用当前帧之前一定帧段的平均能量更新结束帧段平均能量阈值;若判断出是结束帧,则可以根据不同的情况,加权更新变换平滑谱熵阈值、噪声帧能量阈值、减谱功率谱阈值以及结束帧段平均能量阈值。
2.根据权利要求1所述的一种实时语音端点检测方法,其特征在于:去脉冲处理的方法如下:对于一个信号点,考查其前后一段区间内信号点的幅值均值,若此信号点幅值绝对值超过上述幅值均值一定比例,则认为当前信号点幅值是瞬间过大或者过小,有窄脉冲的特性,此时根据上述幅值将当前信号点幅值调小,幅值符号不变;若当前信号点幅值绝对值不超过上述均值的一定比例,则当前信号点不做处理。
3.根据权利要求2所述的一种实时语音端点检测方法,其特征在于:加窗处理中,可以是汉明窗和汉宁窗的一种,在频谱变换前对一帧信号做一个按信号点的幅值加权处理,以使频谱平滑。
4.根据权利要求3所述的一种实时语音端点检测方法,其特征在于:减谱处理的方法如下:设置开始检测的前一定帧段当成噪声帧,用于估计初始底噪,这个帧段可以是从打开检测功能即开始处理,也可以是之后一定帧再处理,以消除打开检测功能的按键开关以及参数初始化处理带来的干扰,使谱熵计算更准确;可取前10帧,作为噪声段,得出平均功率,作为减谱初值;在10帧前,直接计算每帧的功率谱谱熵,10帧后,先减谱处理再计算谱熵。
5.根据权利要求4所述的一种实时语音端点检测方法,其特征在于:自适应重新划分子带的方法为:
1)若某个子带功率谱概率大于0.9,则在有效频带中去掉这个子带的所有频点,即有效频带点数变少,重新计算各个子带的功率谱概率;若没有子带功率谱概率大于0.9,则直接执行2);
2)计算最小子带概率的熵,用Nmin表示,设:
Nmin=-pmin*lg(pmin),设pmin是最小子带功率谱概率;
根据Nmin取值,重新计算子带个数,设重新计算的子带个数为V,采用以下公式:
得出新的子带个数,其中第二个表达式最终得出结果向下取整数,然后,根据有效频带的频点数对每个子带均匀分频点数,设有效频带点数为N'vald,则重新分到每个子带的频点数为:
floor(N'vald/V);
floor函数表示向下取整,/表示做除法;上述过程即是自适应重新划分子带;
再次计算个子带的概率,最终得出这一帧信号的谱熵值,即是自适应划分子带的谱熵值。
6.根据权利要求5所述的一种实时语音端点检测方法,其特征在于:步骤八中的平滑处理的方法为:以当前帧向前算起连续5帧,设谱熵值依次为S(i-4),S(i-3),S(i-2),S(i-1),S(i),即S(i)为当前帧的谱熵,平滑的计算为:
当i小于5时,即第5帧以前,不做平滑处理。
7.根据权利要求6所述的一种实时语音端点检测方法,其特征在于:交替搜寻起始帧和结束帧的方法为:初始状态为搜寻起始帧标志,当搜寻到起始帧时,将标志设置为搜寻结束帧标志,最小语音段之后开始搜索结束帧,当搜寻到结束帧时,将搜寻标志设置为搜寻起始帧标志,从下一帧起搜寻起始帧,如此交替搜寻和设置。
8.根据权利要求7所述的一种实时语音端点检测方法,其特征在于:特定噪声特性为汽车发动机、汽车喇叭、警报、敲击桌子、文具、餐具、键盘、空调和电扇的噪声特性的一种或多种。
9.根据权利要求8所述的一种实时语音端点检测方法,其特征在于:当判断出语音起始帧时,加入信号延时处理。
10.一种实时语音端点检测装置,其特征在于:包括:
信号分帧加重模块,用于接收信号数据,对其进行分帧、加重处理;
去脉冲处理模块,与信号分帧加重模块连接,用于消除一些脉冲噪声;
去除直流模块,与去脉冲处理模块连接,用于去除直流成分;
短时能量和过零率计算模块,与去直流模块连接,用于计算每一帧的短时能量和过零率;
加窗处理模块,与短时能量和过零率计算模块连接,用于在频谱变换前对一帧信号做一个按信号点的幅值加权处理,以使频谱平滑;
减谱处理模块,与加窗处理模块连接,用于使谱熵计算更准确;
谱熵计算模块,与减谱处理模块连接,用于计算自适应划分子带的谱熵值;
变换平滑谱熵计算模块,与谱熵计算模块连接,用于平滑、变换处理谱熵,得到变换平滑谱熵;
语音噪声帧初步判断模块,与变换平滑谱熵处理模块连接,用于初步判断当前帧是语音还是噪声,更新变换平滑谱熵阈值、语音帧和噪声帧能量帧阈值;
语音段处理模块,与语音噪声帧初步判断模块连接,用于起始帧、结束帧的判断;
阈值调整模块,与语音段处理模块连接,用于根据不同情况判断出每帧信号情况后,根据判断结果调整相应阈值。
CN201811491292.8A 2018-12-07 2018-12-07 一种实时语音端点检测方法及装置 Active CN109545188B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811491292.8A CN109545188B (zh) 2018-12-07 2018-12-07 一种实时语音端点检测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811491292.8A CN109545188B (zh) 2018-12-07 2018-12-07 一种实时语音端点检测方法及装置

Publications (2)

Publication Number Publication Date
CN109545188A true CN109545188A (zh) 2019-03-29
CN109545188B CN109545188B (zh) 2021-07-09

Family

ID=65853102

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811491292.8A Active CN109545188B (zh) 2018-12-07 2018-12-07 一种实时语音端点检测方法及装置

Country Status (1)

Country Link
CN (1) CN109545188B (zh)

Cited By (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110189747A (zh) * 2019-05-29 2019-08-30 大众问问(北京)信息科技有限公司 语音信号识别方法、装置及设备
CN110211596A (zh) * 2019-05-29 2019-09-06 哈尔滨工程大学 一种基于Mel子带谱熵鲸目动物哨声信号检测方法
CN110232933A (zh) * 2019-06-03 2019-09-13 Oppo广东移动通信有限公司 音频检测方法、装置、存储介质及电子设备
CN110364187A (zh) * 2019-07-03 2019-10-22 深圳华海尖兵科技有限公司 一种语音信号的端点识别方法和装置
CN110444194A (zh) * 2019-06-28 2019-11-12 北京捷通华声科技股份有限公司 一种语音检测方法和装置
CN110600010A (zh) * 2019-09-20 2019-12-20 上海优扬新媒信息技术有限公司 一种语料提取方法及装置
CN110706717A (zh) * 2019-09-06 2020-01-17 西安合谱声学科技有限公司 一种基于麦克风阵列板的人声检测定向方法
CN110706693A (zh) * 2019-10-18 2020-01-17 浙江大华技术股份有限公司 语音端点的确定方法及装置、存储介质、电子装置
CN110827858A (zh) * 2019-11-26 2020-02-21 苏州思必驰信息科技有限公司 语音端点检测方法及系统
CN110931048A (zh) * 2019-12-12 2020-03-27 广州酷狗计算机科技有限公司 语音端点检测方法、装置、计算机设备及存储介质
CN111276164A (zh) * 2020-02-15 2020-06-12 中国人民解放军空军特色医学中心 飞机上高噪音环境自适应话音激活检测装置及方法
CN111429932A (zh) * 2020-06-10 2020-07-17 浙江远传信息技术股份有限公司 语音降噪方法、装置、设备及介质
CN111755028A (zh) * 2020-07-03 2020-10-09 四川长虹电器股份有限公司 一种基于基音特征的近场遥控器语音端点检测方法及系统
CN111918196A (zh) * 2019-05-08 2020-11-10 腾讯科技(深圳)有限公司 一种音频采集器录音异常的诊断方法、装置、设备及存储介质
CN112189232A (zh) * 2019-07-31 2021-01-05 深圳市大疆创新科技有限公司 音频处理方法及装置
CN112367128A (zh) * 2020-10-30 2021-02-12 电子科技大学 一种基于相隔窗的突发信号频谱检测仪及盲检测方法
CN112420079A (zh) * 2020-11-18 2021-02-26 青岛海尔科技有限公司 语音端点检测方法和装置、存储介质及电子设备
CN112489692A (zh) * 2020-11-03 2021-03-12 北京捷通华声科技股份有限公司 语音端点检测方法和装置
CN112543159A (zh) * 2020-11-30 2021-03-23 中国电子科技集团公司第五十四研究所 一种噪声调频信号识别方法
CN112669877A (zh) * 2020-09-09 2021-04-16 珠海市杰理科技股份有限公司 噪声检测及压制方法、装置、终端设备和系统、芯片
CN112951243A (zh) * 2021-02-07 2021-06-11 深圳市汇顶科技股份有限公司 语音唤醒方法、装置、芯片、电子设备及存储介质
WO2021146857A1 (zh) * 2020-01-20 2021-07-29 深圳市大疆创新科技有限公司 音频处理方法及装置
CN113763985A (zh) * 2021-10-15 2021-12-07 广州蓝仕威克医疗科技有限公司 一种基于声音识别生理现象的装置
CN114242116A (zh) * 2022-01-05 2022-03-25 成都锦江电子系统工程有限公司 一种语音的话音与非话音的综合判决方法
CN110600018B (zh) * 2019-09-05 2022-04-26 腾讯科技(深圳)有限公司 语音识别方法及装置、神经网络训练方法及装置
CN114678037A (zh) * 2022-04-13 2022-06-28 北京远鉴信息技术有限公司 一种重叠语音的检测方法、装置、电子设备及存储介质
CN114900246A (zh) * 2022-05-25 2022-08-12 中国电子科技集团公司第十研究所 噪声基底估计方法、装置、设备及存储介质
CN115376548A (zh) * 2022-07-06 2022-11-22 华南理工大学 一种音频信号有声段端点检测方法及系统
CN116631063A (zh) * 2023-05-31 2023-08-22 武汉星巡智能科技有限公司 基于用药行为识别的老人智能看护方法、装置及设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101599269A (zh) * 2009-07-02 2009-12-09 中国农业大学 语音端点检测方法及装置
CN102522081A (zh) * 2011-12-29 2012-06-27 北京百度网讯科技有限公司 一种检测语音端点的方法及系统
CN106653062A (zh) * 2017-02-17 2017-05-10 重庆邮电大学 一种低信噪比环境下基于谱熵改进的语音端点检测方法
CN107731223A (zh) * 2017-11-22 2018-02-23 腾讯科技(深圳)有限公司 语音活性检测方法、相关装置和设备
GB2554943A (en) * 2016-10-16 2018-04-18 Sentimoto Ltd Voice activity detection method and apparatus
CN108538310A (zh) * 2018-03-28 2018-09-14 天津大学 一种基于长时信号功率谱变化的语音端点检测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101599269A (zh) * 2009-07-02 2009-12-09 中国农业大学 语音端点检测方法及装置
CN102522081A (zh) * 2011-12-29 2012-06-27 北京百度网讯科技有限公司 一种检测语音端点的方法及系统
GB2554943A (en) * 2016-10-16 2018-04-18 Sentimoto Ltd Voice activity detection method and apparatus
CN106653062A (zh) * 2017-02-17 2017-05-10 重庆邮电大学 一种低信噪比环境下基于谱熵改进的语音端点检测方法
CN107731223A (zh) * 2017-11-22 2018-02-23 腾讯科技(深圳)有限公司 语音活性检测方法、相关装置和设备
CN108538310A (zh) * 2018-03-28 2018-09-14 天津大学 一种基于长时信号功率谱变化的语音端点检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘华平等: "一种改进的自适应子带谱熵语音端点检测方法", 《系统仿真学报》 *
王晓华等: "基于时频参数融合的自适应语音端点检测算法", 《计算机工程与应用》 *

Cited By (44)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111918196B (zh) * 2019-05-08 2022-04-19 腾讯科技(深圳)有限公司 一种音频采集器录音异常的诊断方法、装置、设备及存储介质
CN111918196A (zh) * 2019-05-08 2020-11-10 腾讯科技(深圳)有限公司 一种音频采集器录音异常的诊断方法、装置、设备及存储介质
CN110211596B (zh) * 2019-05-29 2021-04-06 哈尔滨工程大学 一种基于Mel子带谱熵鲸目动物哨声信号检测方法
CN110211596A (zh) * 2019-05-29 2019-09-06 哈尔滨工程大学 一种基于Mel子带谱熵鲸目动物哨声信号检测方法
CN110189747A (zh) * 2019-05-29 2019-08-30 大众问问(北京)信息科技有限公司 语音信号识别方法、装置及设备
CN110232933A (zh) * 2019-06-03 2019-09-13 Oppo广东移动通信有限公司 音频检测方法、装置、存储介质及电子设备
CN110232933B (zh) * 2019-06-03 2022-02-22 Oppo广东移动通信有限公司 音频检测方法、装置、存储介质及电子设备
CN110444194A (zh) * 2019-06-28 2019-11-12 北京捷通华声科技股份有限公司 一种语音检测方法和装置
CN110444194B (zh) * 2019-06-28 2021-08-20 北京捷通华声科技股份有限公司 一种语音检测方法和装置
CN110364187A (zh) * 2019-07-03 2019-10-22 深圳华海尖兵科技有限公司 一种语音信号的端点识别方法和装置
CN112189232A (zh) * 2019-07-31 2021-01-05 深圳市大疆创新科技有限公司 音频处理方法及装置
CN110600018B (zh) * 2019-09-05 2022-04-26 腾讯科技(深圳)有限公司 语音识别方法及装置、神经网络训练方法及装置
CN110706717B (zh) * 2019-09-06 2021-11-09 西安合谱声学科技有限公司 一种基于麦克风阵列板的人声检测定向方法
CN110706717A (zh) * 2019-09-06 2020-01-17 西安合谱声学科技有限公司 一种基于麦克风阵列板的人声检测定向方法
CN110600010B (zh) * 2019-09-20 2022-05-17 度小满科技(北京)有限公司 一种语料提取方法及装置
CN110600010A (zh) * 2019-09-20 2019-12-20 上海优扬新媒信息技术有限公司 一种语料提取方法及装置
CN110706693B (zh) * 2019-10-18 2022-04-19 浙江大华技术股份有限公司 语音端点的确定方法及装置、存储介质、电子装置
CN110706693A (zh) * 2019-10-18 2020-01-17 浙江大华技术股份有限公司 语音端点的确定方法及装置、存储介质、电子装置
CN110827858A (zh) * 2019-11-26 2020-02-21 苏州思必驰信息科技有限公司 语音端点检测方法及系统
CN110931048B (zh) * 2019-12-12 2024-04-02 广州酷狗计算机科技有限公司 语音端点检测方法、装置、计算机设备及存储介质
CN110931048A (zh) * 2019-12-12 2020-03-27 广州酷狗计算机科技有限公司 语音端点检测方法、装置、计算机设备及存储介质
WO2021146857A1 (zh) * 2020-01-20 2021-07-29 深圳市大疆创新科技有限公司 音频处理方法及装置
CN111276164B (zh) * 2020-02-15 2021-08-03 中国人民解放军空军特色医学中心 飞机上高噪音环境自适应话音激活检测装置及方法
CN111276164A (zh) * 2020-02-15 2020-06-12 中国人民解放军空军特色医学中心 飞机上高噪音环境自适应话音激活检测装置及方法
CN111429932A (zh) * 2020-06-10 2020-07-17 浙江远传信息技术股份有限公司 语音降噪方法、装置、设备及介质
CN111755028A (zh) * 2020-07-03 2020-10-09 四川长虹电器股份有限公司 一种基于基音特征的近场遥控器语音端点检测方法及系统
CN112669877A (zh) * 2020-09-09 2021-04-16 珠海市杰理科技股份有限公司 噪声检测及压制方法、装置、终端设备和系统、芯片
CN112669877B (zh) * 2020-09-09 2023-09-29 珠海市杰理科技股份有限公司 噪声检测及压制方法、装置、终端设备和系统、芯片
CN112367128A (zh) * 2020-10-30 2021-02-12 电子科技大学 一种基于相隔窗的突发信号频谱检测仪及盲检测方法
CN112489692A (zh) * 2020-11-03 2021-03-12 北京捷通华声科技股份有限公司 语音端点检测方法和装置
CN112420079B (zh) * 2020-11-18 2022-12-06 青岛海尔科技有限公司 语音端点检测方法和装置、存储介质及电子设备
CN112420079A (zh) * 2020-11-18 2021-02-26 青岛海尔科技有限公司 语音端点检测方法和装置、存储介质及电子设备
CN112543159B (zh) * 2020-11-30 2022-04-22 中国电子科技集团公司第五十四研究所 一种噪声调频信号识别方法
CN112543159A (zh) * 2020-11-30 2021-03-23 中国电子科技集团公司第五十四研究所 一种噪声调频信号识别方法
CN112951243A (zh) * 2021-02-07 2021-06-11 深圳市汇顶科技股份有限公司 语音唤醒方法、装置、芯片、电子设备及存储介质
CN113763985B (zh) * 2021-10-15 2024-04-02 广州蓝仕威克医疗科技有限公司 一种基于声音识别生理现象的装置
CN113763985A (zh) * 2021-10-15 2021-12-07 广州蓝仕威克医疗科技有限公司 一种基于声音识别生理现象的装置
CN114242116A (zh) * 2022-01-05 2022-03-25 成都锦江电子系统工程有限公司 一种语音的话音与非话音的综合判决方法
CN114678037B (zh) * 2022-04-13 2022-10-25 北京远鉴信息技术有限公司 一种重叠语音的检测方法、装置、电子设备及存储介质
CN114678037A (zh) * 2022-04-13 2022-06-28 北京远鉴信息技术有限公司 一种重叠语音的检测方法、装置、电子设备及存储介质
CN114900246A (zh) * 2022-05-25 2022-08-12 中国电子科技集团公司第十研究所 噪声基底估计方法、装置、设备及存储介质
CN115376548A (zh) * 2022-07-06 2022-11-22 华南理工大学 一种音频信号有声段端点检测方法及系统
CN116631063A (zh) * 2023-05-31 2023-08-22 武汉星巡智能科技有限公司 基于用药行为识别的老人智能看护方法、装置及设备
CN116631063B (zh) * 2023-05-31 2024-05-07 武汉星巡智能科技有限公司 基于用药行为识别的老人智能看护方法、装置及设备

Also Published As

Publication number Publication date
CN109545188B (zh) 2021-07-09

Similar Documents

Publication Publication Date Title
CN109545188A (zh) 一种实时语音端点检测方法及装置
Moattar et al. A simple but efficient real-time voice activity detection algorithm
CN105118502B (zh) 语音识别系统的端点检测方法及系统
US7756700B2 (en) Perceptual harmonic cepstral coefficients as the front-end for speech recognition
CN103886871B (zh) 语音端点的检测方法和装置
EP1083541B1 (en) A method and apparatus for speech detection
US6721699B2 (en) Method and system of Chinese speech pitch extraction
CN111128213B (zh) 一种分频段进行处理的噪声抑制方法及其系统
Pang Spectrum energy based voice activity detection
Labied et al. An overview of automatic speech recognition preprocessing techniques
Özaydın Examination of energy based voice activity detection algorithms for noisy speech signals
CN110379438B (zh) 一种语音信号基频检测与提取方法及系统
Sorin et al. The ETSI extended distributed speech recognition (DSR) standards: client side processing and tonal language recognition evaluation
JPS60114900A (ja) 有音・無音判定法
JP3091537B2 (ja) 音声パターン作成方法
Li et al. Robust speech endpoint detection based on improved adaptive band-partitioning spectral entropy
Li et al. Sub-band based log-energy and its dynamic range stretching for robust in-car speech recognition
KR100194953B1 (ko) 유성음 구간에서 프레임별 피치 검출 방법
Lee et al. Objective pathological voice quality assessment based on HOS features
Cui Pitch extraction based on weighted autocorrelation function in speech signal processing
Zhang An energy-based adaptive voice detection approach
JPH0795239B2 (ja) 音声フレーム中の基本周波数の存在を検出する装置および方法
CN114974231A (zh) 一种噪声环境下的基音周期提取方法
Touazi et al. A Case Study on Back-End Voice Activity Detection for Distributed Specch Recognition System Using Support Vector Machines
Zhu et al. Speech endpoint detection method based on logarithmic energy entropy product of adaptive sub-bands in low signal-to-noise ratio environments

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