CN113766073B - 会议系统中的啸叫检测 - Google Patents
会议系统中的啸叫检测 Download PDFInfo
- Publication number
- CN113766073B CN113766073B CN202111051519.9A CN202111051519A CN113766073B CN 113766073 B CN113766073 B CN 113766073B CN 202111051519 A CN202111051519 A CN 202111051519A CN 113766073 B CN113766073 B CN 113766073B
- Authority
- CN
- China
- Prior art keywords
- howling
- teleconference
- state
- audio data
- term
- 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.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 84
- 238000000034 method Methods 0.000 claims abstract description 186
- 230000003595 spectral effect Effects 0.000 claims abstract description 54
- 238000004891 communication Methods 0.000 claims description 46
- 230000002776 aggregation Effects 0.000 claims description 12
- 238000004220 aggregation Methods 0.000 claims description 12
- 230000000116 mitigating effect Effects 0.000 claims description 12
- 230000001131 transforming effect Effects 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 5
- 230000002123 temporal effect Effects 0.000 abstract description 17
- 230000008569 process Effects 0.000 description 70
- 238000001228 spectrum Methods 0.000 description 31
- 238000000605 extraction Methods 0.000 description 17
- 238000010801 machine learning Methods 0.000 description 17
- 230000007774 longterm Effects 0.000 description 15
- 230000000875 corresponding effect Effects 0.000 description 14
- 239000000203 mixture Substances 0.000 description 14
- 238000009499 grossing Methods 0.000 description 11
- 230000005236 sound signal Effects 0.000 description 10
- 230000000694 effects Effects 0.000 description 8
- 230000000737 periodic effect Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000012805 post-processing Methods 0.000 description 5
- 238000005070 sampling Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 238000003066 decision tree Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R3/00—Circuits for transducers, loudspeakers or microphones
- H04R3/02—Circuits for transducers, loudspeakers or microphones for preventing acoustic reaction, i.e. acoustic oscillatory feedback
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04M—TELEPHONIC COMMUNICATION
- H04M3/00—Automatic or semi-automatic exchanges
- H04M3/42—Systems providing special services or facilities to subscribers
- H04M3/56—Arrangements for connecting several subscribers to a common circuit, i.e. affording conference facilities
- H04M3/568—Arrangements for connecting several subscribers to a common circuit, i.e. affording conference facilities audio processing specific to telephonic conferencing, e.g. spatial distribution, mixing of participants
-
- 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/48—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04M—TELEPHONIC COMMUNICATION
- H04M9/00—Arrangements for interconnection not involving centralised switching
- H04M9/08—Two-way loud-speaking telephone systems with means for conditioning the signal, e.g. for suppressing echoes for one or both directions of traffic
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04M—TELEPHONIC COMMUNICATION
- H04M2203/00—Aspects of automatic or semi-automatic exchanges
- H04M2203/20—Aspects of automatic or semi-automatic exchanges related to features of supplementary services
- H04M2203/2094—Proximity
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Acoustics & Sound (AREA)
- Computational Linguistics (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- General Health & Medical Sciences (AREA)
- Otolaryngology (AREA)
- Telephonic Communication Services (AREA)
- Circuit For Audible Band Transducer (AREA)
- Telephone Function (AREA)
Abstract
本公开涉及会议系统中的啸叫检测。所公开的远程会议方法涉及在远程会议期间检测啸叫状态,所述远程会议涉及两个或更多个远程会议客户端位置以及远程会议服务器。所述远程会议服务器被配置用于在所述远程会议客户端位置之间提供全双工音频连接。所述啸叫状态是涉及某一远程会议客户端位置中的两个或更多个远程会议装置的声反馈状态。检测所述啸叫状态涉及对远程会议音频数据的频谱特性和时间特性两者的分析。所公开的远程会议方法涉及确定哪个客户端位置正在导致所述啸叫状态并且涉及减轻所述啸叫状态或者发送啸叫状态检测消息。
Description
本申请是申请日为2018年9月27日、申请号为201880062458.6、发明名称为“会议系统中的啸叫检测”的发明专利申请的分案申请。
相关申请的交叉引用
本申请要求2017年9月29日提交的国际专利申请号PCT/CN2017/104604、2017年10月23日提交的美国临时申请号62/575,701和2017年11月10日提交的欧洲专利申请号17201011.8的优先权,所述专利申请中的每一个都通过引用以其全文并入本文。
技术领域
本公开涉及音频信号的处理。特别地,本公开涉及处理与远程会议或视频会议相关的音频信号。
背景技术
声反馈是当音频输入(例如,麦克风)与音频输出(例如,扬声器)之间存在声音环路时发生的一种类型的反馈。例如,麦克风接收的信号可能被扬声器放大和再现。来自扬声器的再现声音然后可能被麦克风再次接收,进一步放大,并且然后被扬声器再次以更高的幅度或音量水平再现。在这种类型的系统中,声反馈的声音可以是响亮的尖叫或长声尖叫,这在本文中可以称为“啸叫(howl)”。在一些实例中,吉他手和其他表演者可能故意制造啸叫(例如,在吉他拾音器与扬声器之间)以便产生期望的音乐效果。虽然一些方法已经被证明成功地检测和减轻无意中产生的啸叫,但是已经证明在远程会议的情景下要解决这个问题更加困难。
发明内容
本文公开了各种方法。一些远程会议方法可以涉及在涉及两个或更多个客户端位置的远程会议期间检测啸叫状态以及确定哪个客户端位置正在导致所述啸叫状态。所述远程会议可以涉及两个或更多个远程会议客户端位置和一个远程会议服务器。所述远程会议服务器可以被配置用于在所述远程会议客户端位置之间提供全双工音频连接。所述啸叫状态可以是涉及某一远程会议客户端位置中的两个或更多个远程会议装置的声反馈状态。检测所述啸叫状态可能涉及对远程会议音频数据的频谱特性和时间特性两者的分析。在一些示例中,所述方法可以涉及减轻所述啸叫状态和/或发送啸叫状态检测消息。
根据一些示例,所述远程会议服务器可以确定哪个客户端位置正在导致所述啸叫状态。在一些实施方式中,所述远程会议服务器可以检测所述啸叫状态。
在一些这样的示例中,远程会议客户端位置处的远程会议装置可以对所述远程会议音频数据中的至少一些远程会议音频数据执行啸叫特征提取过程。所述啸叫特征提取过程可以产生啸叫特征数据。所述远程会议装置可以将所述啸叫特征数据发送到所述远程会议服务器。所述远程会议服务器可以至少部分地基于所述啸叫特征数据来检测所述啸叫状态。
根据一些实施方式,所述远程会议服务器可以至少部分地基于从所述远程会议客户端位置中的一个或多个远程会议客户端位置接收的远程会议音频数据来检测所述啸叫状态。在一些这样的实施方式中,所述远程会议服务器可以创建包括从所述远程会议客户端位置接收的远程会议音频数据的远程会议音频数据混合。在一些示例中,所述远程会议服务器可以基于所述远程会议音频数据混合来检测所述啸叫状态。
根据一些这样的示例,所述远程会议服务器可以从所述远程会议客户端位置接收远程会议元数据。所述远程会议服务器可以至少部分地基于所述远程会议元数据来确定哪个客户端位置正在导致所述啸叫状态。所述远程会议元数据可以例如包括语音活动检测元数据、电平元数据和/或能量元数据。
在一些示例中,所述方法可以涉及根据所述远程会议音频数据来估计啸叫存在概率。检测所述啸叫状态可以至少部分地基于所述啸叫存在概率。根据一些这样的示例,所述啸叫存在概率估计可以基于分层规则集或机器学习方法。
根据一些实施方式,远程会议客户端位置处的启用啸叫检测的远程会议装置可以检测所述啸叫状态。在一些这样的实施方式中,所述启用啸叫检测的远程会议装置可以减轻所述啸叫状态。
在一些示例中,检测所述啸叫状态可以涉及:根据频带的顺序统计量来计算基于功率的指标;计算频谱共振指标;计算百分位数间距指标;事件聚合;计算周期性指标;检测包络相似度;计算频谱峰度指标;基于机器学习的过程;和/或估计啸叫存在概率。
本文描述的一些或所有方法可以由一个或多个装置根据存储在一个或多个非暂态介质上的指令(例如,软件)来执行。这种非暂态介质可以包括诸如本文所述的存储器装置,包括但不限于随机存取存储器(RAM)装置、只读存储器(ROM)装置等。因此,本公开中描述的主题的各种创新方面可以在其上存储有软件的非暂态介质中实施。例如所述软件可以包括用于控制至少一个装置来处理音频数据的指令。例如所述软件可以由如本文所公开的那些控制系统等控制系统的一个或多个部件来执行。
所述软件可以例如包括用于执行本文公开的远程会议方法中的一种或多种远程会议方法的指令。一些这样的远程会议方法可以涉及在远程会议期间检测啸叫状态以及确定哪个客户端位置正在导致所述啸叫状态。所述远程会议可以涉及两个或更多个远程会议客户端位置和一个远程会议服务器。所述远程会议服务器可以被配置用于在所述远程会议客户端位置之间提供全双工音频连接。所述啸叫状态可以是涉及某一远程会议客户端位置中的两个或更多个远程会议装置的声反馈状态。检测所述啸叫状态可能涉及对远程会议音频数据的频谱特性和时间特性两者的分析。在一些示例中,所述方法可以涉及减轻所述啸叫状态和/或发送啸叫状态检测消息。
本公开的至少一些方面可以经由设备来实施。例如,一个或多个装置可以被配置用于至少部分地执行本文所公开的方法。在一些实施方式中,设备可以包括接口系统和控制系统。所述接口系统可以包括一个或多个网络接口、所述控制系统与存储器系统之间的一个或多个接口、所述控制系统与另一装置之间的一个或多个接口和/或一个或多个外部装置接口。所述控制系统可以包括通用单芯片或多芯片处理器、数字信号处理器(DSP)、专用集成电路(ASIC)、现场可编程门阵列(FPGA)或其他可编程逻辑装置、离散门或晶体管逻辑、或离散硬件部件中的至少一个。
在一些实施方式中,所述设备可以是远程会议服务器,或者可以是远程会议服务器的一部分。在一些这样的实施方式中,所述设备可以包括远程会议服务器的一个或多个线路卡或刀片。根据一些这样的示例,所述设备可以包括接口系统和控制系统。所述接口系统可以被配置用于所述远程会议服务器与所述远程会议客户端位置之间的通信。
例如,所述控制系统可以被配置用于在远程会议期间在两个或更多个远程会议客户端位置之间提供全双工音频连接。所述控制系统可以被配置用于在所述远程会议期间检测啸叫状态。所述啸叫状态可以是涉及远程会议客户端位置中的两个或更多个远程会议装置的声反馈状态。检测所述啸叫状态可能涉及对远程会议音频数据的频谱特性和时间特性两者的分析。在一些示例中,所述控制系统可以被配置用于确定哪个客户端位置正在导致所述啸叫状态。在一些实施方式中,所述控制系统可以被配置用于减轻所述啸叫状态和/或发送啸叫状态检测消息。
根据一些实施方式,所述控制系统可以被配置用于至少部分地基于从所述远程会议客户端位置中的一个或多个远程会议客户端位置接收的远程会议音频数据来检所述测啸叫状态。在一些示例中,所述控制系统可以被配置用于创建远程会议音频数据混合。所述远程会议音频数据混合可以包括从所述远程会议客户端位置接收的远程会议音频数据。在一些示例中,所述控制系统可以被配置用于基于所述远程会议音频数据混合来检测所述啸叫状态。在一些这样的示例中,所述控制系统可以被配置用于:经由所述接口系统从所述远程会议客户端位置接收远程会议元数据并且至少部分地基于所述远程会议元数据来确定哪个客户端位置正在导致所述啸叫状态。在一些实施方式中,所述远程会议元数据可以包括语音活动检测元数据、电平元数据和/或能量元数据。根据一些实施方式,所述控制系统可以被配置用于确定哪个客户端装置正在导致所述啸叫状态。
在一些示例中,所述控制系统可以被配置用于根据从所述远程会议客户端位置中的一个或多个远程会议客户端位置接收的远程会议音频数据来估计啸叫存在概率。根据一些这样的示例,检测所述啸叫状态可以至少部分基于所述啸叫存在概率。在一些实例中,所述啸叫存在概率估计可以基于分层规则集和/或机器学习过程。
根据一些实施方式,检测所述啸叫状态可以涉及:根据频带的顺序统计量来计算基于功率的指标;计算频谱共振指标;计算百分位数间距指标;事件聚合;计算周期性指标;检测包络相似度;计算频谱峰度指标;基于机器学习的过程;和/或估计啸叫存在概率。
在一些示例中,所述控制系统可以被配置用于经由所述接口系统从远程会议客户端位置处的远程会议装置接收啸叫特征数据。所述啸叫特征数据可以是由所述远程会议装置执行的啸叫特征提取过程的结果。所述控制系统可以被配置用于至少部分地基于所述啸叫特征数据来检测所述啸叫状态。
根据一些实施方式,所述设备可以是远程会议客户端装置,或者可以是其一部分。在一些示例中,所述远程会议客户端装置可以是电话、会议电话、膝上型电脑等。在一些示例中,所述远程会议客户端装置可以包括接口系统和控制系统。所述接口系统可以被配置用于所述远程会议客户端装置与远程会议服务器之间的通信。在一些实施方式中,所述远程会议服务器可以被配置用于在远程会议期间在多个远程会议客户端位置之间提供全双工音频连接。
在一些实例中,所述控制系统可以被配置用于在所述远程会议期间检测啸叫状态。所述啸叫状态可以是涉及所述远程会议客户端装置和包括所述远程会议客户端装置的远程会议客户端位置中的至少一个其他装置的声反馈状态。在一些示例中,检测所述啸叫状态可以涉及对远程会议音频数据的频谱特性和时间特性两者的分析。根据一些示例,所述控制系统可以被配置用于减轻所述啸叫状态。
根据一些实施方式,检测所述啸叫状态可以涉及:根据频带的顺序统计量来计算基于功率的指标;计算频谱共振指标;计算百分位数间距指标;事件聚合;计算周期性指标;检测包络相似度;计算频谱峰度指标;基于机器学习的过程;和/或估计啸叫存在概率。
在一些示例中,所述控制系统可以被配置用于在所述远程会议期间执行啸叫状态特征提取过程。所述啸叫状态特征提取过程可以涉及提取与检测啸叫状态相关的啸叫特征数据。在一些实施方式中,所述啸叫状态特征提取过程可以涉及对远程会议音频数据的频谱特性和时间特性两者的分析。根据一些这样的实施方式,所述控制系统可以被配置用于将所述啸叫特征数据发送到所述远程会议服务器。
根据一些实施方式,执行啸叫状态特征提取过程可以涉及:根据频带的顺序统计量来计算基于功率的指标;计算频谱共振指标;计算百分位数间距指标;事件聚合;计算周期性指标;检测包络相似度;计算频谱峰度指标;基于机器学习的过程;和/或估计啸叫存在概率。
在以下附图和说明中阐述了本说明书所描述的主题的一个或多个实施方式的细节。从所述描述和附图中,其他特征、方面和优点将变得显而易见。注意,以下附图的相对尺寸可能不是按比例来绘制的。
附图说明
图1示出了在涉及两个远程会议客户端位置的远程会议期间的远程会议装置的示例。
图2在时域和频域中示出了可能涉及远程会议客户端位置中的两个或更多个远程会议装置的啸叫状态的示例的图。
图3是示出可以被配置成用于执行本文所公开方法中的至少一些方法的设备的部件的示例的框图。
图4是概述根据一个示例的方法的框的流程图。
图5示出了根据一个实施方式的远程会议客户端装置和远程会议服务器配置。
图6示出了根据另一实施方式的远程会议客户端装置和远程会议服务器配置。
图7示出了根据又一实施方式的远程会议客户端装置和远程会议服务器配置。
图8是概述根据一些实施方式的啸叫检测方法的框的流程图。
图9是示出啸叫事件周期性检测方法的示例框的流程图。
图10是示出测量啸叫事件周期性的另一种方法的示例框的流程图。
图11是示出根据一个示例的啸叫存在概率确定方法的框的流程图。
图12是示出确定信号包络相似度的方法的示例框的流程图。
图13是概述根据一些实施方式的啸叫检测方法的框的流程图。
在各附图中相同的附图标记和名称指示相似的元件。
具体实施方式
以下说明涉及出于描述本公开的一些创新性方面的目的的某些实施方式以及这些创新性方面所可以实施于的背景的示例。然而,可以以各种不同的方式来应用本文的教导内容。此外,所描述的实施例可以在各种硬件、软件、固件等中实施。例如,本申请的各方面可以至少部分地体现在设备、包括一个以上装置的系统、方法、计算机程序产品等中。因此,本申请的各方面可以采用硬件实施例、软件实施例(包括固件、常驻软件、微代码等)和/或结合了软件方面和硬件方面两者的实施例的形式。这些实施例在本文中可以被称为“电路”、“模块”或“引擎”。本申请的一些方面可以采用在一个或多个非暂态介质中实施的计算机程序产品的形式,所述非暂态介质具有在其上实施的计算机可读程序代码。这种非暂态介质例如可以包括硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦除可编程只读存储器(EPROM或闪速存储器)、便携式致密盘只读存储器(CD-ROM)、光存储装置、磁存储装置、或前述各项的任何合适的组合。因此,本公开的教导内容不旨在限制附图中示出和/或本文所描述的实施方式,而是具有广泛的适用性。
图1示出了在涉及两个远程会议客户端位置的远程会议期间的远程会议装置的示例。在这个示例中,远程会议服务器被配置用于在远程会议客户端位置1中的远程会议客户端装置100a与远程会议客户端位置2中的远程会议客户端装置100b和100c之间提供音频连接(可以是全双工音频连接)。这里,远程会议客户端装置100a和100c是会议电话,而远程会议客户端装置100b是膝上型电脑。
根据这个示例,远程会议客户端装置100b和100c两者启用了它们的麦克风和扬声器。因此,远程会议客户端装置100b和100c的麦克风各自向远程会议服务器提供对应于近端讲话者2的话音的远程会议音频数据。因此,远程会议服务器向远程会议客户端装置100b提供由远程会议客户端装置100c接收的与近端讲话者2的话音相对应的远程会议音频数据,并且还向远程会议客户端装置100c提供由远程会议客户端装置100b接收的与近端讲话者2的话音相对应的远程会议音频数据。
远程会议客户端装置100b和100c的扬声器在与远程会议音频数据从远程会议客户端位置2传输到远程会议服务器并返回的往返时间相对应的时间延迟之后各自再现近端讲话者2的话音。当远程会议客户端装置100b的扬声器再现近端讲话者2的话音时,远程会议客户端装置100c的麦克风检测到近端讲话者2的再现话音。同样,当远程会议客户端装置100c的扬声器再现近端讲话者2的话音时,远程会议客户端装置100b的(多个)麦克风检测到近端讲话者2的再现话音。
因此,远程会议客户端位置2中的近端讲话者2的话音可以触发啸叫状态。然而,所述啸叫状态具有如果远程会议客户端装置100b和100c没有参加远程会议则不会出现的至少一些特性。
图2在时域和频域中示出了可能涉及远程会议客户端位置中的两个或更多个远程会议装置的啸叫状态的示例的曲线图。曲线图205在时域中,时间由水平轴表示,并且电平以dB为单位,由垂直轴表示。曲线图205示出了啸叫事件210a至210f。啸叫事件210a至210f发生的时间在本文中可以称为“啸叫状态”。
可以观察到啸叫事件210a至210f是周期性的或准周期性的,并且它们的幅度随着时间而增加。在这个示例中,啸叫事件210a至210f的周期性与时间延迟相对应,所述时间延迟与远程会议音频数据从远程会议客户端位置传输到远程会议服务器并返回的往返时间相对应。还可以观察到,话音事件的电平(在啸叫事件210a至210f的左边)低于啸叫事件210a至210f的电平。
曲线图215在频域中表示相同的音频数据,包括啸叫事件210a至210f。从曲线图215可以观察到,啸叫事件210a至210f的大部分能量集中在相对窄的频带中。
基于发明人进行的许多观察,包括曲线图205和215中的数据,已经观察到反馈啸叫不同于话音的以下频谱和时间特性:
·啸叫事件的共振频率包括单个频带或多个频带的频谱峰;
·啸叫事件的幅度通常比话音的幅度大;
·啸叫事件在本质上通常是周期性或准周期性的并且在频率上高度局部化;以及
·啸叫事件的特性具有跨时间的一致性。
下文描述了啸叫状态的其他特性。如上所述,啸叫状态具有如果参与产生啸叫状态的装置没有参加远程会议则不会出现的至少一些特性。例如,当麦克风太靠近扬声器时和/或当放大倍数太大时,公共广播音频系统的该麦克风与该扬声器之间可能会出现声反馈。这种声反馈通常会表现为单一且连续的啸叫事件,而不是多个不连续和周期性发生的啸叫事件。与公共广播系统情景中的声反馈相比,双工通信用例中的声反馈可能跨越不同的网络和不同的装置。往返时间通常大于300ms,并且可能在500ms至一秒或更长的范围内。因此,远程会议情景中的声反馈检测和减轻可能与公共广播系统情景中的声反馈检测和减轻非常不同,在公共广播系统情景中,反馈是在同一麦克风-放大器-扬声器系统中导致的。
鉴于上述观察,本公开的一些方面可以提供用于在远程会议期间检测啸叫状态的改进方法。图3是示出可以被配置成用于执行本文所公开方法中的至少一些方法的设备的部件的示例的框图。在一些示例中,设备305可以是远程会议服务器。在其他示例中,设备305可以是远程会议客户端装置。可以经由硬件、经由存储在非暂态介质上的软件、经由固件和/或通过其组合来实施设备305的部件。图3以及本文所公开的其他附图中所示出的部件的类型和数量仅是通过示例的方式来示出的。替代性实施方式可以包括更多、更少和/或不同的部件。例如,如果设备305是远程会议客户端装置,则设备305可以包括一个或多个麦克风、一个或多个扬声器等。
在这个示例中,设备305包括接口系统310和控制系统315。接口系统310可以包括一个或多个网络接口、控制系统315与存储器系统之间的一个或多个接口和/或一个或多个外部装置接口(诸如一个或多个通用串行总线(USB)接口)。在一些实施方式中,接口系统310可以包括用户接口系统。用户接口系统可以被配置用于接收来自用户的输入。在一些实施方式中,用户接口系统可以被配置用于向用户提供反馈。例如,用户接口系统可以包括具有相应触摸和/或手势检测系统的一个或多个显示器。在一些示例中,用户接口系统可以包括一个或多个扬声器。根据一些示例,用户接口系统可以包括用于提供触觉反馈的设备,诸如马达、振动器等。控制系统315例如可以包括通用单芯片或多芯片处理器、数字信号处理器(DSP)、专用集成电路(ASIC)、现场可编程门阵列(FPGA)或其他可编程逻辑装置、离散门或晶体管逻辑、和/或离散硬件部件。
在一些示例中,可以在单一装置中实施设备305。然而,在一些实施方式中,可以在多于一个装置中实施设备305。在一些这样的实施方式中,控制系统315的功能可以包括在多于一个装置中。在一些示例中,设备305可以是另一设备的部件。例如,在一些实施方式中,设备305可以是服务器的部件,例如线路卡。
图4是概述根据一个示例的方法的框的流程图。在一些实例中,可以通过图3的设备或通过本文所公开的另一种类型的设备来执行所述方法。在一些示例中,可以经由存储在一个或多个非暂态介质上的软件来实施方法400的框。与本文所描述的其他方法一样,不必以所指示的顺序来执行方法400的框。而且,这种方法可以包括比所示出和/或所描述的框更多或更少的框。
在这个实施方式中,框405涉及在远程会议期间检测啸叫状态。根据这个示例,远程会议涉及两个或更多个远程会议客户端位置和被配置用于在远程会议客户端位置之间提供全双工音频连接的远程会议服务器。这里,啸叫状态是涉及某一远程会议客户端位置中的两个或更多个远程会议装置的声反馈状态。
在这个示例中,检测啸叫状态涉及对远程会议音频数据的频谱特性和时间特性两者的分析。根据特定实施方式,检测啸叫状态可以涉及:根据频带的顺序统计量来计算基于功率的指标;计算频谱共振指标;计算百分位数间距指标;事件聚合;计算周期性指标;检测包络相似度;计算频谱峰度指标;基于机器学习的过程;和/或估计啸叫存在概率。术语“指标”和“度量”在本文可以同义使用。在一些示例中,检测啸叫状态可以例如涉及根据远程会议音频数据来估计啸叫存在概率。啸叫存在概率估计可以基于分层规则集或机器学习方法。检测啸叫状态可以至少部分地基于啸叫存在概率。下文详细描述了各种示例。
在一些实例中,远程会议服务器可以在框405中检测啸叫状态。然而,在一些实施方式中,远程会议客户端位置处的启用啸叫检测的远程会议装置可以在框405中检测啸叫状态。
在这个示例中,框410涉及确定哪个客户端位置正在导致啸叫状态。在一些实例中,远程会议服务器可以在框410中确定哪个客户端位置正在导致啸叫状态。
然而,在一些实施方式中,远程会议客户端位置处的启用啸叫检测的远程会议装置可以在框410中确定哪个客户端位置正在导致啸叫状态。在一些实例中,启用啸叫检测的远程会议装置可以确定启用啸叫检测的远程会议装置本身或同一远程会议客户端位置处的另一远程会议装置正在导致啸叫状态。
可选框415涉及减轻啸叫状态或发送啸叫状态检测消息。根据一些示例,远程会议服务器可以在框415中向正在导致啸叫状态的客户端位置发送啸叫状态检测消息。在一些实施方式中,减轻啸叫状态可以涉及检测啸叫状态的一个或多个峰值频带以及应用对应于一个或多个峰值频带的一个或多个陷波滤波器。在一些示例中,减轻啸叫状态可以涉及使正在导致啸叫状态的客户端位置处的一个或多个麦克风静音和/或降低正在导致啸叫状态的客户端位置处的扬声器音量。在一些实施方式中,正在导致啸叫状态的客户端位置处的启用啸叫检测的远程会议装置可以被配置成用于减轻啸叫状态。
图5示出了根据一个实施方式的远程会议客户端装置和远程会议服务器配置。在这个示例中,远程会议服务器被配置成用于至少部分地基于从远程会议客户端位置接收到的远程会议音频数据来检测啸叫状态。这里,远程会议客户端位置1至M正在参与远程会议。
在图5所示的示例中,远程会议服务器包括啸叫检测器1至M,每个啸叫检测器与远程会议客户端位置或客户端位置处的远程会议客户端装置相对应。根据这个实施方式,啸叫检测器1至M中的每一个被配置用于在远程会议期间检测啸叫状态。在这个示例中,远程会议服务器还包括被配置用于确定哪个客户端位置正在导致啸叫状态的啸叫客户端检测器。
这里,啸叫检测器1至M中的每一个被配置用于在检测到啸叫状态之后向啸叫客户端检测器发送啸叫状态信号。在这个实施方式中,啸叫客户端检测器可以根据从其接收到啸叫状态信号的啸叫检测器来确定哪个客户端位置正在导致啸叫状态。远程会议服务器可以被配置用于减轻啸叫状态和/或向导致啸叫状态的远程会议客户端位置发送啸叫状态检测消息。
本文示出和描述的啸叫检测器和啸叫客户端检测器,包括但不限于图5中示出的那些,可以根据特定实施方式以各种方式实施。例如,啸叫检测器和啸叫客户端检测器可以被实施为控制系统(诸如图3所示和上文描述的控制系统)的一部分。在一些示例中,啸叫检测器和啸叫客户端检测器可以经由硬件、经由软件(可以包括固件、驻留软件、微代码等)来实施和/或以结合了软件方面和硬件方面的实施例来实施。这些实施例在本文中可以被称为“电路”、“模块”或“引擎”。在一些实施方式中,啸叫检测器和啸叫客户端检测器可以经由实施在一个或多个非暂态介质中的计算机程序产品来实施,所述非暂态介质具有实施于其上的计算机可读程序代码。这种非暂态介质例如可以包括硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦除可编程只读存储器(EPROM或闪速存储器)、便携式致密盘只读存储器(CD-ROM)、光存储装置、磁存储装置、或前述各项的任何合适的组合。
图6示出了根据另一实施方式的远程会议客户端装置和远程会议服务器配置。在这个示例中,远程会议服务器被配置成用于至少部分地基于从正在参加远程会议的远程会议客户端位置1至M接收的远程会议音频数据来检测啸叫状态。在图6所示的示例中,远程会议服务器包括混合器,所述混合器被配置成用于混合从远程会议客户端位置1至M接收的远程会议音频数据并将远程会议音频数据混合输出到啸叫检测器。在这个实施方式中,啸叫检测器被配置成用于基于远程会议音频数据混合来检测啸叫状态。
与图5所示的类型的实施方式相比,这种类型的实施方式具有潜在的优势,因为不需要多个啸叫检测器。因此,与图5所示的类型的实施方式相比,这种类型的实施方式可以提供相对较低的计算开销、较少的功耗等。
与前一个示例一样,图6中所示的远程会议服务器也包括被配置用于确定哪个客户端位置正在导致啸叫状态的啸叫客户端检测器。远程会议服务器可以被配置用于减轻啸叫状态和/或向导致啸叫状态的远程会议客户端位置发送啸叫状态检测消息。
因为在这个示例中只有一个啸叫检测器,因此啸叫客户端检测器不能根据从其接收到啸叫状态信号的啸叫检测器来确定哪个客户端位置正在导致啸叫状态。因此,根据一些这样的实施方式,远程会议服务器可以从参与远程会议的远程会议客户端位置1至M中的每一个接收远程会议元数据。在一些示例中,远程会议元数据可以包括语音活动检测元数据、电平元数据和/或能量元数据。啸叫客户端检测器可以被配置成用于至少部分地基于远程会议元数据来确定哪个客户端位置正在导致啸叫状态。例如,在一些实例中,啸叫客户端检测器可以被配置成用于至少部分地基于指示哪个远程会议客户端位置具有当前正在讲话的远程会议参与者的语音活动检测元数据来确定哪个客户端位置正在导致啸叫状态。这是一个合理的假设,因为在远程会议或其他谈话的大部分时间里,一次只有一个人在说话。
一些实施方式可以提供啸叫客户端检测算法,其也可以被称为“责备”算法,以基于服务器中的啸叫检测结果和由远程会议服务器接收的每个远程会议客户端的元数据信息来将啸叫的区域与特定远程会议客户端相关联。在一个示例中,“责备”算法可以如下:
·在远程会议音频数据混合上运行啸叫检测器,其可以表示为xsum(n)=sum(xm(n)),其中n表示帧号并且xm(n)表示指示第m个远程会议客户端的电平的信号。
·当啸叫检测器观察到超过阈值概率的啸叫事件的概率时:
o选择VAD=真的一个或多个当前远程会议客户端;
o运行混合信号xsum(n)与第m个远程会议客户端xm(n)之间的互相关Rm(n),诸如Rm(n)=α*Rm(n)+(1-α)*Rm(n-1)。α的值可以根据特定实施方式而变化。例如,α的值可以取决于帧速率和/或取决于保持Rm的过去值的期望时间长度。在一些示例中,α可以是0.95、0.96、0.97、0.98、0.99等。根据一些实施方式,可以例如使用具有一秒或更多秒数量级的时间常数的1极平滑滤波器来平滑Rm(n)。这个平滑过程将在相对较长的时间段内保持相关信息。
·选择所有Rm(n)中的最大值,该值应与其电平信号与混合信号最相关的远程会议客户端m相对应。这第m个客户端成为检测到的正在导致啸叫状态的远程会议客户端(或客户端位置)。
图7示出了根据又一实施方式的远程会议客户端装置和远程会议服务器配置。在这个示例中,远程会议客户端位置中的每个包括啸叫检测器。包括这种啸叫检测器的远程会议装置在本文中可以被称为远程会议客户端位置处的启用啸叫检测的远程会议装置。在一些示例中,启用啸叫检测的远程会议装置可以被配置成用于检测啸叫状态。在一些这样的示例中,启用啸叫检测的远程会议装置(或远程会议客户端位置处的另一装置)可以被配置成用于减轻啸叫状态。根据一些这样的示例,不需要远程会议服务器的参与来检测和/或减轻啸叫状态。
然而,在一些替代实施方式中,仍然可能存在某种程度的远程会议服务器参与。根据一些这样的实施方式,远程会议客户端位置处的远程会议装置(诸如图7所示的客户端侧啸叫检测器之一)可以对由远程会议装置接收的远程会议音频数据中的至少一些远程会议音频数据执行啸叫特征提取过程。
啸叫特征提取过程可以产生啸叫特征数据。啸叫特征数据的特性可以根据特定实施方式的啸叫特征提取过程的细节而变化。例如,啸叫特征数据可以包括与基于功率的指标、频谱共振指标、百分位数间距指标、事件聚合、周期性指标、包络相似度、频谱峰度指标和/或啸叫存在概率相对应的一个或多个特征。
远程会议装置可以将啸叫特征数据中的至少一些啸叫特征数据发送到远程会议服务器。根据一些这样的示例,远程会议服务器可以被配置成用于至少部分地基于啸叫特征数据来检测啸叫状态。远程会议服务器可以被配置成用于至少部分地基于从其接收啸叫特征数据的客户端位置来确定哪个客户端位置正在导致啸叫状态。
本文公开了各种啸叫检测方法。一些啸叫检测方法可以涉及检测啸叫的短期频谱结构特征、啸叫的中期频谱-时间结构特征和/或啸叫的长期时间重复模式。这种啸叫检测方法可以涉及啸叫存在概率估计或“二元”啸叫/无啸叫确定。下文的详细讨论将描述以下啸叫检测方法的方面,所述方法包括:
o可以计算基于对数域频带功率的度量的顺序统计量以便检测啸叫的短期频谱结构;
o可以计算频谱共振度量以便检测啸叫的短期频谱结构;
o可以计算一致的高百分位数间距度量以便检测啸叫的频谱-时间结构特征;
o可以计算一致的高共振度量以便检测啸叫的频谱-时间结构特征;
o事件(诸如潜在的啸叫事件)可以被聚合以便检测啸叫的长期时间重复模式;
o可以分析基于调制频谱的周期性度量以便检测啸叫的长期时间重复模式;
o可以分析包络相似度以便检测啸叫的长期时间重复模式;
o啸叫存在概率估计可以基于包括两个或更多个上述计算的分层规则集(在本文中可以称为分层决策树)。
将在以下小节中描述这些方法。在一些方法中,可以将音频信号变换到频域并且可以基于频域音频信号来计算特征。然而,一些方法至少部分地在时域中操作。
图8是概述根据一些实施方式的啸叫检测方法的框的流程图。方法800的框可以例如由图3和图5至7中所示的装置中的一个或多个装置来实施。在一些示例中,可以经由存储在一个或多个非暂态介质上的软件来实施方法800的框。与本文所描述的其他方法一样,不必以所指示的顺序来执行方法800的框。而且,这种方法可以包括比所示出和/或所描述的框更多或更少的框。
在这个示例中,方法800开始于在框805中将音频数据从时域变换到频域。音频数据可以包括远程会议音频数据。在这个实施方式中,在框810中将经变换的音频数据划分成频带。频带的数量可以根据特定实施方式而变化。一些实施方式(其可以对应于相对窄带的音频数据)可以涉及将经变换的音频数据划分成40个音频带,而其他实施方式可以涉及将经变换的音频数据划分成48个音频带(这可以对应于宽带音频数据)、56个音频带(这可以对应于超宽带音频数据)或者其他数量的音频带。根据一些示例,框810可以涉及将经变换的音频数据划分成对数间隔的带,而其他示例可以涉及将经变换的音频数据划分成线性间隔的带。
在这个示例中,经变换且频带划分的音频数据被提供给三个不同的啸叫检测过程:框815涉及检测频谱结构的短期特征,框820涉及检测中期频谱和时间特征,并且框825涉及检测在相对长的时期内重复的模式。下文提供了这些啸叫检测过程中的每一个的示例。在一些实施方式中,“中期”可以是数百毫秒的时间段,例如150ms、200ms、250ms、300ms、350ms、400ms等。“长期”意味着比中期更长。例如,“长期”可以是秒的数量级,诸如1秒、1.5秒、2秒、2.5秒等。“短期”意味着比中期更短。例如,“短期”可以对应于几十毫秒,例如10ms、20ms、30ms等。在一些示例中,“短期”事件可以在单独的音频数据帧中被检测到。
其他示例可能涉及更多或更少的啸叫检测过程。例如,一些替代方法可能只涉及图8所示的三种啸叫检测过程中的两种。一些方法可以仅涉及三种类型的啸叫检测过程中的一种,例如,仅框815或框820。
在这个示例中,来自框815、820和825的结果被提供给框830中的啸叫存在概率估计过程。在这个示例中,框830输出啸叫存在概率的指示。可选的后处理框835接收啸叫存在概率的指示并输出指示啸叫状态或非啸叫状态的啸叫二元指示符。在一些实施方式中,方法800还可以包括响应于啸叫状态的动作,诸如减轻啸叫状态或发送啸叫状态检测消息。
以下段落提供了啸叫检测过程和指标的示例。
一致高峰度指标
这些指标提供了框815的过程的示例。第一个这样的指标,在本文中被称为频谱峰度度量(SPM),基于这样的观察,即啸叫音频通常具有非常高的频谱峰度,能量集中在窄频率范围内。有许多潜在的频谱相关度量可以用来表征信号频谱的“尖峰”条件。频谱峰度度量的一个实施例可以在下面的等式(1)中看到。
在等式1中,X(k)表示频带索引k的频谱幅度,l表示帧的索引,并且K表示频带的数量。在这个示例中,SPM被计算为一减去带功率谱的几何平均值与带功率谱的算术平均值之间的比率。
在一些示例中,等式1的SPM度量可以被平滑如下:
在等式2中,N表示块大小并且Ts表示平滑时间常数。N和Ts的值可以根据特定实施方式而变化并且可以取决于例如采样率和/或帧长度。在帧长度为20ms的一些示例中,对于16K采样率,N可以是320,并且对于32K采样率,N可以是640。在其他示例中,帧长度可以是10ms、30ms等。根据一些实施方式,Ts可以是0.04秒、0.06秒、0.10秒、0.12秒等。
在一些实施方式中,可以使用其他频谱相关度量来表征信号频谱的尖峰条件,诸如:
·基于峰值带之和与其他带之间的功率比的频谱波峰度量;
·基于峰均带功率比的频谱波峰度量(在一些示例中,可以在没有一个或多个峰带的情况下计算平均功率比);
·相邻谱带功率的方差或标准偏差;
·相邻频带之间谱带功率差的总和或最大值;
·频谱质心(spectral centroid)周围的频谱延展度或频谱方差;和/或
·频谱熵指标。
可用于啸叫状态检测的另一种峰度度量在本文中可称为“基于计数器的一致高峰度度量”。这些实施方式涉及基于计数器来检测相对长的具有高峰度信号的块。在一些这样的实施方式中,如果信号的尖峰没有变得高度峰,则计数器可以缓慢释放。等式3显示了基于计数器的一致高峰度度量的一个示例:
在等式3中,PCounter(l)表示高峰度信号的计数,PMaxCounter表示当观察到具有高峰度信号的长块时的最大计数器值,表示用于确定高峰度条件的阈值,/>表示下降的/>的阈值,并且PCounterDown表示用于递减计数器的步度。SPM_Th、SPM_DownTh和P_CounterDown的值可以根据特定实施方式而变化。在一些示例中,SPM_Th可以是0.85、0.88、0.90等,SPM_DownTh可以是0.05、0.08、0.10等,并且P_CounterDown可以是0.2、0.5、1.0等。例如,P_CounterDown的值可以取决于期望的释放时间。
一致高IPR指标
百分位数间距(IPR)指标可用于检测当前帧中的带功率之间的对比度。与话音和噪声信号相比,啸叫音频通常具有非常高的IPR。
对数域带功率的顺序统计量
对数域带功率可表示如下:
b(k,l)=10log 10(|X(k,l)|2),其中k=1,...,K;l=1,...,N。在这个表达式中,k表示带的索引并且l表示帧的索引。带功率的相应顺序统计量版本可以分别表示为B(1,l)、B(2,l)、……B(K,l),其中:
B(1,l)≥B(2,l)≥...≥B(K,l) (4)
表达IPR的一种方式如下:
IPR(l)=B(r,l)-B(κ,l) (5)
r和κ的值已经由发明人根据各种啸叫条件下的实验确定。在一些示例中,r=1并且κ可以表示为在一些实施方式中,PercentileTh=80。在其他实施方式中,PercentileTh可以是82、85等。
在一些示例中,原始IPR度量可以例如如下被平滑:
在等式6中,αIPR可以表示为其中N表示块大小并且TS表示平滑时间常数。如上所述,N和TS的值可以根据特定实施方式而变化并且可以取决于采样率和/或帧长度。在帧长度为20ms的一些示例中,对于16K采样率,N可以是320,并且对于32K采样率,N可以是640。在其他示例中,帧长度可以是10ms、30ms等。根据一些实施方式,Ts可以是0.04秒、0.06秒、0.10秒、0.12秒等。
IPR的加权函数
在一些实施方式中,一个或多个某种类型的加权函数可以用于IPR指标。使用精心选择的加权函数可以提高原始IPR度量的鉴别能力。在一些示例中,可以在应用自动增益控制(AGC)之前确定输入电平。在一些实施方式中,输入电平和频谱峰度指标可以用于加权函数中。等式7提供了一个这样的示例:
在等式7中,inLeveldB表示信号的输入电平。在一些示例(其可以是在远程会议客户端装置中实施啸叫检测器的示例)中,inLeveldB可以表示应用AGC之前的信号的输入电平。然而,在其他示例(其可以是在远程会议服务器中实施啸叫检测器的示例)中,inLeveldB可以表示应用AGC之后的信号的输入电平。在一个示例中,LowTh=-60dB并且HighTh=-30dB。应当理解,LowTh和HighTh的值可以根据特定实施方式而变化。IPR的一个相应的加权函数可以表示如下:
WIPR(l)=WLevel(l)×SPM(l) (8)
IPRComb度量是上文描述的IPR度量的加权版本并且在一些示例中可以表示如下:
另一个有用的啸叫检测指标可以被称为“基于计数器的一致高IPRComb度量”。这个指标可以表示如下:
在等式10中,IPRCounter(l)表示高IPRComb信号的计数器,IPRMaxCounter表示当观察到具有高IPRComb信号的长块时的最大计数器,表示高IPRComb指标的阈值,并且IPRCounterDown表示下降的计数器的步度。在一些示例中,IPR_Th可以是12、15、18、20等。“长块”的示例是大约10或20帧的块,其持续时间将取决于帧长度。例如,如果帧长度是20ms,则大约10或20帧的块将对应于200到400毫秒的持续时间。因此,这个指标提供了可以在图8的框820中发生的过程的示例。
基于事件的中期计数器
这个指标提供了可以在框820中发生的过程的另一个示例。基于中期事件的计数器可以在中期期间检测啸叫音频的频谱结构的重复模式。在一个示例中,可以基于PCounter(l)和IPRCounter(l)来检测啸叫事件的开始,如下:
EventOnset=(IPRCounter(l)≥IPRCounterTh)&&(PCounter(l)≥PCounterTh) (11)
中期事件计数器可以表示如下:
在等式12中,EventCounter(l)表示中期啸叫事件的计数器,EventMaxCuner表示中期啸叫事件的最大计数器,并且EventCounterDown表示下降的事件计数器的步度。在这个示例中,计数器用于检测中期内啸叫频谱结构的重复模式。这里,如果计数器的值增加,啸叫存在的可能性也增加。在这个示例中,最大计数器表示啸叫存在的最高可能的可能性,例如1。在这个实施方式中,CounterDown表示用于减小计数器的步长。当计数器值减小时,这表明啸叫存在的可能性正在下降。
长期信号可变性度量
这个指标是可以与图8的框825相对应的过程的示例。在一些示例中,经平滑带功率指标可以如下确定:
在一些实施方式中,时域分析窗口期间的频谱熵可以表示如下:
η(k,l)本质上是在第l帧结束的R个连续帧上在带k处计算的归一化短期频谱的熵度量。根据一些示例,长期信号可变性度量LTSV可以如下计算:
在这个示例中,信号可变性度量LTSV是K个带的熵(根据等式14计算)的样本方差。在等式(15)中,
基于事件的长期周期性
这个指标是可以与图8的框825相对应的过程的另一个示例。发明人已经观察到,在远程会议环境中的啸叫经常显示出(例如)可以发生在两秒或更多秒的时间段内的长期重复的模式。
图9是示出啸叫事件周期性检测方法的示例框的流程图。方法900的框可以例如由图3和图5至7中所示的装置中的一个或多个装置来实施。在一些示例中,可以经由存储在一个或多个非暂态介质上的软件来实施方法900的框。与本文所描述的其他方法一样,不必以所指示的顺序来执行方法900的框。而且,这种方法可以包括比所示出和/或所描述的框更多或更少的框。
在这个示例中,方法900开始于在框905中将接收到的音频数据从时域变换到频域。音频数据可以包括远程会议音频数据。在这个实施方式中,还在框905中将经变换的音频数据划分成频带。
为了对啸叫事件进行计数并确定啸叫事件长度,需要某种类型的啸叫事件开始检测。在这个示例中,经变换和频带划分的音频数据被提供给三个不同的过程,用于确定啸叫开始检测指标。其他实施方式可以涉及确定更多或更少的啸叫开始检测指标。在这个示例中,框910涉及用于确定是否存在作为啸叫事件的特性的电平增加的“Δ电平”过程。在一些示例中,这个电平增加阈值可以是20dB、25dB、30dB等。这里,框915涉及评估IPR。例如,框915可以涉及一种或多种类型的IPR计算方法,诸如上文描述的那些方法。根据这个实施方式,框920涉及一种或多种类型的峰度检测,其可以表示啸叫事件开始的频谱结构。框920可以涉及上文描述的峰度检测方法中的一种或多种。
在这个示例中,从框910至920的过程输出的啸叫开始检测指标被提供给框925中的啸叫事件开始检测过程。在一些示例中,在啸叫事件开始检测过程期间,从框910至920的过程输出的啸叫开始检测指标被给予相等的权重,而在其他示例中,一个指标可以被给予比其他指标更大的权重。当已经检测到啸叫事件的开始时,向框930的过程输入啸叫事件开始检测指示,其中确定啸叫事件长度。例如,在框930中,可以调用一种或多种啸叫事件检测方法,例如,本文公开的短期啸叫事件检测方法中的一种或多种。根据这个示例,在框935中已经确定了啸叫事件长度之后,方法900进行到啸叫事件计数器过程(框935)。
在这个实施方式中,啸叫事件计数器过程的结果(其可以包括检测到的啸叫事件数量)被提供给框940的一个或多个事件周期性测量过程。在一些示例中,框940可以涉及例如如下应用一个或多个启发式规则来得到周期性度量:
EventPeriodicity=max(0,min(1,(Eventcount-EventNumTh)×PeriodicityStep)) (17)
在等式17中,EventCount表示检测到的啸叫事件数量;EventNumTh表示啸叫事件存在数量的阈值并且periodicityStep表示周期性的增加步度。根据一个示例,EventNumTh是2并且PeriodicityStep是0.5。然而,在替代示例中,EventNumTh可以是1.5、1.8等。在其他示例中,PeriodicitySetp可以是0.4、0.6等。
根据一些示例,可以例如如下调用递归平均方法以便提供基于事件的周期性度量的进一步平滑:
在等式18中,αsmooth、表示平滑系数。在一些示例中,αsmooth可以是0.85、0.9、0.95等。
基于电平的长期周期性
一些实施方式涉及确定基于音频信号电平的长期周期性指标。发明人已经观察到,在远程会议环境中啸叫事件的幅度调制频谱经常显示出长期重复的模式。
图10是示出测量啸叫事件周期性的另一种方法的示例框的流程图。方法1000的框可以例如由图3和图5至7中所示的装置中的一个或多个装置来实施。在一些示例中,可以经由存储在一个或多个非暂态介质上的软件来实施方法1000的框。与本文所描述的其他方法一样,不必以所指示的顺序来执行方法1000的框。而且,这种方法可以包括比所示出和/或所描述的框更多或更少的框。
在这个示例中,方法1000开始于在框1005中接收音频数据。音频数据可以包括远程会议音频数据。在这个实施方式中,接收到的音频数据被提供给框1010。这里,框1010涉及根据滑动时间窗口来评估音频数据中的电平信息。在一个示例中,滑动时间窗口可以是音频数据的128帧。在其他示例中,滑动时间窗口可以更长或更短。在一些示例中,可以在AGC过程之前确定电平信息。
根据这个示例,框1010的结果被提供给框1015,框1015涉及调制频谱分析的过程。在一些示例中,框1015涉及快速傅立叶变换(FFT)计算。例如,框1015的调制频谱分析可以涉及滑动时间窗口中的长期电平信息的傅立叶变换。框1015的结果被提供给框1020的峰值拾取过程。在框1020中选择调制频谱的峰值,因为峰值表示周期性信息。然后,可以在框1025中基于框1020的结果来确定周期性度量。例如,在框1020中选择的峰值电平可以指示啸叫事件的开始。相应的峰值电平可以在稍后时间被标识并且可以使其与较早的峰值电平相关。
啸叫存在概率
在一些实施方式中,局部啸叫存在概率可以基于IPR值的线性插值,例如,如下:
可替代地或附加地,可以确定全局啸叫存在概率。在一些示例中,全局啸叫存在概率可以根据基于计数器的一致性度量,如下:
在一些示例中,原始啸叫概率可以定义如下:
P(l)=PLocal(l)PGlobal(l) (21)
图11是示出根据一个示例的啸叫存在概率确定方法的框的流程图。方法1100的框可以例如由图3和图5至7中所示的装置中的一个或多个装置来实施。在一些示例中,可以经由存储在一个或多个非暂态介质上的软件来实施方法1100的框。与本文所描述的其他方法一样,不必以所指示的顺序来执行方法1100的框。而且,这种方法可以包括比所示出和/或所描述的框更多或更少的框。
在这个示例中,方法1100的显著特征在框1105的啸叫事件检测过程之后开始,所述啸叫事件检测过程可以根据任何上文描述的啸叫事件检测方法来执行。尽管未在图11中示出,但是根据特定实施方式,框1105可以涉及其他过程,诸如接收音频数据、将音频数据变换到频域、频带划分等。这里,框1110涉及确定啸叫事件的数量是否已经超过第一啸叫事件计数器阈值,该阈值也可以被称为中期事件计数器高阈值。在一些示例中,中期事件计数器高阈值可以是10、12、13等。如果在框1110中确定啸叫事件的数量已经超过第一啸叫事件计数器阈值,则过程继续到框1115。
然而,如果在框1110中确定啸叫事件的数量没有超过第一啸叫事件计数器阈值,则过程继续到框1150,其中确定啸叫事件的数量是否小于第二啸叫事件计数器阈值,第二啸叫事件计数器阈值也可以被称为中期事件计数器低阈值。在一些示例中,中期事件计数器低阈值可以是5、6、7等。如果在框1150中确定啸叫事件的数量小于第二啸叫事件计数器阈值,则过程继续到框1155,框1155涉及降低当前啸叫存在概率(HPP)。在这个示例中,框1155涉及将当前HPP值乘以HPP_Clear,HPP_Clear是如果假警报的概率高则用于快速减小HPP值的常数。在一些示例中,HPP_Clear的值可以是0.1、0.2、0.3等。
然而,如果在框1150中确定啸叫事件的数量不小于第二啸叫事件计数器阈值,则过程继续到框1160。在这个示例中,框1160涉及将当前HPP值设置为零或值(HPP(1)[当前HPP值]-HPP_Down)中的较大者。根据这个示例,HPP_Down是用于减小HPP的值的常数。在一些示例中,HPP_Down的值可以是0.02、0.03、0.04等。
如上所述,如果在框1110中确定啸叫事件的数量已经超过第一啸叫事件计数器阈值,则过程继续到框1115。这里,框1115涉及评估信号可变性。在这个示例中,框1115涉及确定长期信号可变性(LTSV)是否大于LTSV阈值,所述LTSV阈值在图11中表示为LTSV_Th1。在一些示例中,LTSV_Th1的值可以是0.15、0.2、0.25等。如果在框1115中确定LTSV不大于LTSV_Th1,则过程进行到框1150。然而,如果在框1115中确定LTSV大于LTSV_Th1,则过程继续到框1120。
在这个实施方式中,框1120涉及确定音频输入数据的一个或多个电平值是否超过阈值。在一些示例中,可以在AGC过程之前确定一个或多个电平值。根据这个实施方式,框1120涉及确定是否InLeveldB>Level_Th1。在一些示例中,LTSV_Th1的值可以是-60Db、-65Db、-70Db等。如果在框1120中确定InLeveldB不大于LTSV_Th1,则过程进行到框1150。然而,如果在框1120中确定InLeveldB>Level_Th1,则过程在这个示例中继续到框1125。
在这个实施方式中,框1125涉及将当前HPP值设置为1或当前HPP值与“原始”啸叫概率之和中的较小者,在一些示例中,“原始”啸叫概率可以根据等式21来确定。过程然后进行到框1130。
在这个实施方式中,框1130涉及确定啸叫事件周期性指标是否小于第一啸叫事件周期性阈值,所述阈值在图11中表示为Periodicity_Th1。在一些示例中,Periodicity_Th1的值可以是0.2、0.3、0.4、0.5等。如果在框1130中确定啸叫事件周期性指标小于Periodicity_Th1,则过程继续到框1135。这里,框1135涉及减小当前HPP。在这个示例中,框1135涉及将当前HPP值乘以HPP_Clear。
然而,如果在框1130中确定啸叫事件周期性指标不小于第一啸叫事件周期性阈值,则过程在这个示例中进行到框1140。这里,框1140涉及确定啸叫事件周期性指标是否大于第二啸叫事件周期性阈值,该阈值在图11中表示为Periodicity_Th2。在一些示例中,Periodicity_Th2的值可以是0.80、0.85、0.90等。在这个示例中,如果在框1140中确定啸叫事件周期性指标不大于Periodicity_Th2,则将不改变当前HPP值。然而,如果在框1140中确定啸叫事件周期性指标大于Periodicity_Th2,则过程继续到框1145。
在这个实施方式中,框1145涉及将当前HPP值设置为1或当前HPP值与在图11中表示为HPP_Up的常数之和中的较小者。HPP_Up是用于递增HPP的当前值的值。在一些示例中,HPP_Up的值可以是0.1、0.15、0.2等。
在完成上述过程之后,可以输出HPP值。尽管在图11中未示出,但是当检测到额外啸叫事件时,过程可以继续。
在替代实施方式中,可以根据一种或多种其他方法来确定HPP值。例如,可以根据如回归决策树方法、自适应提升算法、高斯混合模型、隐马尔可夫模型或深度神经网络等基于机器学习的方法来确定HPP值。
确定啸叫二元指示符的后处理
如上文参考图8的框835所描述的,一些实施方式可以涉及在确定HPP值或类似的啸叫概率值之后的后处理步骤。一些这样的实施方式可以涉及基于投票的决策平滑过程。如果我们将帧l中啸叫状态的二元指示符定义为I(l)并且将观察窗R中的投票率定义为V(l),则投票率可以表示如下:
二元决策平滑可以如下执行:
I(l)=I(l)&(V(l)≥Vth) (23)
在等式23中,Vth表示投票率阈值,该阈值可以根据特定实施方式而变化。在一些示例中,Vth可以是0.5、0.6、0.7等。
根据一些实施方式,中值滤波器方法可以用于进一步的二元决策平滑。在一些这样的实施方式中,帧l中啸叫状态的二元指示符可以表示如下:
I(l)=median{I(l-M+1),I(l-M+2),...,I(l)} (24)
在等式24中,M表示中值滤波器窗口的长度。在一些示例中,M可以在3、4、5等的范围内。
时域中基于RMS的包络相似度确定
图12是示出确定信号包络相似度的方法的示例框的流程图。方法1200的框可以例如由图3和图5至7中所示的装置中的一个或多个装置来实施。在一些示例中,可以经由存储在一个或多个非暂态介质上的软件来实施方法1200的框。与本文所描述的其他方法一样,不必以所指示的顺序来执行方法1200的框。而且,这种方法可以包括比所示出和/或所描述的框更多或更少的框。
在这个示例中,方法1200开始于在框1205中接收音频数据。音频数据可以包括远程会议音频数据。在这个实施方式中,接收到的音频数据被提供给框1210,而无需变换到频域。在这个示例中,框1210涉及对时域音频信号加窗。在一个这样的示例中,帧长度是10ms并且首先经由窗口长度为20ms的汉明(Hamming)窗来对时域音频信号加窗。然而,加窗过程的类型和窗口长度可以根据特定实施方式而变化。在一些替代示例中,可以应用汉宁(Hanning)窗。在一些示例中,窗口大小可以是25ms、30ms等。
过程然后继续到框1212,其中计算每个帧的RMS能量。在一个实施方式中,每个帧的RMS能量计算如下:
在等式25中,xi(i=1,2,…,N)是帧的样本。
在这个示例中,过程然后进行到框1215的平滑和滤波过程。在一些实施方式中,框1215涉及通过对框1212中计算的RMS能量应用中值滤波来获得经平滑的RMS包络,例如如下:
xrms(i)=Med[xrms(i-N),...,xrms(i),...,xrms(i+N)] (26)
因此,中值滤波操作可以涉及对RMS能量进行排序并选择第50百分位的值。框1215的结果被提供给框1220。根据这个示例,对于经平滑的RMS包络,框1220涉及用于找到RMS峰值的峰值拾取方法。在一些示例中,峰值拾取方法可以选择局部最大RMS值。在这个示例中,提取两个相邻峰值之间的块并且然后计算剩余块的长度。在一些实施方式中,这个长度可以用作包络周期性的约束条件。
然后在框1225中对RMS块进行归一化。在一些示例中,框1225涉及使用最大RMS值,例如如下:
vrms1=vrms1/max(abs(vrms1)) (27)
vrms2=vrms2/max(abs(vrms2)) (28)
在等式27和28中,vrms1和vrms2表示三个峰值中的两个相邻块。
在这个实施方式中,然后在框1230中计算两个块的平滑RMS的均方误差(MSE)。根据一些示例,框1230可以涉及如下计算MSE:
MSE=sum((vrms1-vrms2).^2)/length(vrms1) (29)
在这个示例中,接下来是框1235的平滑过程,以便减轻峰值之间的MSE的波动。根据一些这样的示例,框1235涉及中值滤波过程以获得经平滑的MSE值。
发明人已观察到,远程会议环境中的相应啸叫事件的长度可能非常相似。通过使用对应于啸叫事件的块的长度,可以构造一种约束规则,所述约束规则可以应用于经平滑的MSE值。以这种方式,可以获得最终MSE值。
因此,根据这个示例,然后在框1240中应用约束规则。在一些示例中,约束规则可以至少部分地基于最大长度以及len1和len2的比率,其中len1和len2分别表示上述vrms1和vrms2的长度。如果最大长度太长或比率太小或太大,这可以表明len1与len2之间的差太大。因此,根据一些示例,可以增加MSE值。如果比率几乎等于1,并且几个块(例如,4个或5个块)上的比率几乎等于1,这可以表明包络具有很强的相似度。因此,根据一些这样的示例,可以减小MSE值。
通过选择和应用MSE值的适当阈值,然后在框1245中可以得到啸叫事件的二元指示。根据一些实施方式,MSE值的阈值可以很小,例如0.015、0.02、0.025、0.03等。在一些示例中,对于啸叫事件,输出可以是1,并且对于非啸叫事件,输出可以是0。
频域中基于能量的包络相似度确定
在一些实施方式中,可以基于能量在频域中计算包络相似度指标。根据帕塞瓦尔(Parseval)定理,
在等式30中,E表示能量。
在确定基于能量的包络相似度的一个示例中,可以对音频信号加窗并且然后例如通过使用FFT或改进的离散余弦变换来变换到频域。在一些示例中,每个帧的能量可以计算如下:
在等式31中,X[k]表示频域中的系数,k表示频率索引,并且N表示帧长度。然后可以应用滤波器以获得平滑的包络。例如,可以应用移动平均滤波器:
在等式32中,N表示移动平均滤波器的长度,其值可能根据实施方式而不同。在一些示例中,N的值可以是3、5、7等。根据一些实施方式,然后可以拾取包络的峰值并且可以计算MSE值。在一些示例中,这些过程可以如上文参考图12的框1220至1245所述来执行。
子带域中基于能量的包络相似度确定
啸叫事件可能发生在几个频带中。因此,子带域中的集体决策可以有助于确定包络相似度。
确定子带域包络相似度的一个实施方式可以通过对接收的音频信号加窗并且然后例如通过使用FFT或改进的离散余弦变换将音频信号变换到频域来开始。可以例如根据心理声学原理将得到的宽带信号划分成子带。子带的数量可以例如基于宽带信号的采样率。在涉及16kHz采样率的一个示例中,宽带信号可以被划分成48个子带。在涉及32kHz采样率的其他示例中,宽带信号可以被划分成56或64个子带。
然后,过程可以基本上如前一节中所描述的那样进行。然而,这个方法的一些示例可以涉及不仅输出二元非啸叫或啸叫事件指示(例如,1和0),还输出啸叫事件确定中的不确定性的指示(例如,0.5)。例如,啸叫事件可能与话音信号混合并且这可能导致模糊的啸叫事件确定结果。
图13是概述根据一些实施方式的啸叫检测方法的框的流程图。方法1300的框可以例如由图3和图5至7中所示的装置中的一个或多个装置来实施。由图13的框引用的过程可以根据任何相应的上文描述的方法来执行。在一些示例中,可以经由存储在一个或多个非暂态介质上的软件来实施方法1300的框。与本文所描述的其他方法一样,不必以所指示的顺序来执行方法1300的框。而且,这种方法可以包括比所示出和/或所描述的框更多或更少的框。
在这个示例中,方法1300开始于在框1305中将音频数据从时域变换到频域。音频数据可以包括远程会议音频数据。在这个实施方式中,在框1310中将经变换的音频数据划分成频带。
在这个示例中,经变换和频带划分的音频数据被提供给四个不同的啸叫检测过程:根据这个实施方式,框1315涉及检测频谱峰度,框1320涉及确定IPR,框1325涉及确定将AGC应用于音频数据之前的信号电平,并且框1330涉及频谱熵确定过程。在这个实施方式中,框1315、1320和1325的结果被提供给IPRComb框,IPRComb框又向高IPRComb一致性度量框1345提供输入。在这个示例中,框1315的结果也被提供给高峰度一致性度量框1340。这里,框1330的结果也被提供给长期信号可变性框1350。
根据这个实施方式,框1325、1340、1345、1350和1355的结果被提供给啸叫存在概率确定框1360。在这个示例中,框1360输出啸叫存在概率度量,所述度量在本文中也被称为HPP。这里,框1360还向后处理框1365提供输入,后处理框1365输出二元啸叫/无啸叫指示。在一些示例中,方法1300可以涉及减轻啸叫状态或发送啸叫状态检测消息。
本领域普通技术人员能够很容易地明白对本公开中所描述实施方式的各种修改。在不脱离本公开范围的情况下,可以将本文所限定的一般原理应用于其他实施方式。因此,本权利要求并不旨在局限于本文所示出的实施方式,而是符合与本公开内容、本文所公开的原理和新颖特征一致的最宽范围。
可以从以下枚举的示例实施例(EEE)中理解本公开的各个方面:
1.一种远程会议方法,包括:
在远程会议期间检测啸叫状态,所述远程会议涉及两个或更多个远程会议客户端位置以及被配置用于在所述远程会议客户端位置之间提供全双工音频连接的远程会议服务器,所述啸叫状态是涉及某一远程会议客户端位置中的两个或更多个远程会议装置的声反馈状态,其中,检测所述啸叫状态涉及对远程会议音频数据的频谱特性和时间特性两者的分析;以及
确定哪个客户端位置正在导致所述啸叫状态。
2.如EEE 1所述的方法,进一步包括减轻所述啸叫状态或发送啸叫状态检测消息中的至少一个。
3.如EEE 1或EEE 2所述的方法,其中,所述远程会议服务器确定哪个客户端位置正在导致所述啸叫状态。
4.如EEE 1至3中任一项所述的方法,其中,所述远程会议服务器检测所述啸叫状态。
5.如EEE 4所述的方法,其中:
远程会议客户端位置处的远程会议装置对所述远程会议音频数据中的至少一些远程会议音频数据执行啸叫特征提取过程,所述啸叫特征提取过程产生啸叫特征数据;
所述远程会议装置将所述啸叫特征数据发送到所述远程会议服务器;并且
所述远程会议服务器至少部分地基于所述啸叫特征数据来检测所述啸叫状态。
6.如EEE 4所述的方法,其中,所述远程会议服务器至少部分地基于从所述远程会议客户端位置中的一个或多个远程会议客户端位置接收的远程会议音频数据来检测所述啸叫状态。
7.如EEE 6所述的方法,其中,所述远程会议服务器创建包括从所述远程会议客户端位置接收的远程会议音频数据的远程会议音频数据混合,并且其中,所述远程会议服务器基于所述远程会议音频数据混合来检测所述啸叫状态。
8.如EEE 7所述的方法,其中,所述远程会议服务器从所述远程会议客户端位置接收远程会议元数据,并且其中,所述远程会议服务器至少部分地基于所述远程会议元数据来确定哪个客户端位置正在导致所述啸叫状态。
9.如EEE 8所述的方法,其中,所述远程会议元数据包括语音活动检测元数据、电平元数据或能量元数据中的一个或多个。
10.如EEE 1所述的方法,进一步包括根据所述远程会议音频数据来估计啸叫存在概率,其中,检测所述啸叫状态至少部分基于所述啸叫存在概率。
11.如EEE 10所述的方法,其中,所述啸叫存在概率估计是基于分层规则集或机器学习方法。
12.如EEE 1所述的方法,其中,远程会议客户端位置处的启用啸叫检测的远程会议装置检测所述啸叫状态。
13.如EEE 12所述的方法,其中,所述启用啸叫检测的远程会议装置减轻所述啸叫状态。
14.如EEE 1至13中任一项所述的方法,其中,检测所述啸叫状态涉及从由以下各项构成的过程列表中选择的一个或多个过程:根据频带的顺序统计量来计算基于功率的指标;计算频谱共振指标;计算百分位数间距指标;事件聚合;计算周期性指标;检测包络相似度;计算频谱峰度指标;基于机器学习的过程;以及估计啸叫存在概率。
15.一个或多个非暂态介质,其上存储有软件,所述软件包括用于控制一个或多个装置执行如EEE 1至14中任一项所述的方法的指令。
16.一种远程会议设备,包括:
用于在远程会议期间检测啸叫状态的装置,所述远程会议涉及两个或更多个远程会议客户端位置以及被配置用于在所述远程会议客户端位置之间提供全双工音频连接的远程会议服务器,所述啸叫状态是涉及某一远程会议客户端位置中的两个或更多个远程会议装置的声反馈状态,其中,检测所述啸叫状态涉及对远程会议音频数据的频谱特性和时间特性两者的分析;以及
用于确定哪个客户端位置正在导致所述啸叫状态的装置。
17.如EEE 16所述的远程会议设备,进一步包括用于减轻所述啸叫状态的装置。
18.如EEE 16或EEE 17所述的远程会议设备,进一步包括用于发送啸叫状态检测消息的装置。
19.如EEE 16至18中任一项所述的远程会议设备,其中,检测所述啸叫状态至少部分基于从所述远程会议客户端位置中的一个或多个远程会议客户端位置接收的远程会议音频数据。
20.如EEE 16至19中任一项所述的远程会议设备,进一步包括用于创建包括从所述远程会议客户端位置接收的远程会议音频数据的远程会议音频数据混合的装置,其中,检测所述啸叫状态至少部分基于所述远程会议音频数据混合。
21.如EEE 16至20中任一项所述的远程会议设备,其中,检测所述啸叫状态至少部分基于从所述远程会议客户端位置接收的远程会议元数据。
22.如EEE 21所述的远程会议设备,其中,所述远程会议元数据包括语音活动检测元数据、电平元数据或能量元数据中的一个或多个。
23.如EEE 16至22中任一项所述的远程会议设备,进一步包括用于根据所述远程会议音频数据来估计啸叫存在概率的装置,其中,检测所述啸叫状态至少部分基于所述啸叫存在概率。
24.如EEE 23所述的远程会议设备,其中,所述啸叫存在概率估计是基于分层规则集或机器学习方法。
25.如EEE 16至24中任一项所述的远程会议设备,其中,检测所述啸叫状态涉及从由以下各项构成的过程列表中选择的一个或多个过程:根据频带的顺序统计量来计算基于功率的指标;计算频谱共振指标;计算百分位数间距指标;事件聚合;计算周期性指标;检测包络相似度;计算频谱峰度指标;基于机器学习的过程;以及估计啸叫存在概率。
26.一种远程会议服务器,包括:
接口系统,所述接口系统被配置用于所述远程会议服务器与远程会议客户端位置之间的通信;以及
控制系统,所述控制系统被配置用于:
在远程会议期间在两个或更多个远程会议客户端位置之间提供全双工音频连接;
在所述远程会议期间检测啸叫状态,所述啸叫状态是涉及某一远程会议客户端位置中的两个或更多个远程会议装置的声反馈状态,其中,检测所述啸叫状态涉及对远程会议音频数据的频谱特性和时间特性两者的分析;以及
确定哪个客户端位置正在导致所述啸叫状态。
27.如EEE 26所述的远程会议服务器,其中,所述控制系统被配置用于减轻所述啸叫状态或发送啸叫状态检测消息中的至少一个。
28.如EEE 26或EEE 27所述的远程会议服务器,其中,所述控制系统被配置用于确定哪个客户端装置正在导致所述啸叫状态。
29.如EEE 26至28中任一项所述的远程会议服务器,其中,所述控制系统被配置用于至少部分地基于从所述远程会议客户端位置中的一个或多个远程会议客户端位置接收的远程会议音频数据来检测所述啸叫状态。
30.如EEE 26至29中任一项所述的远程会议服务器,其中,所述控制系统被配置用于:
创建包括从所述远程会议客户端位置接收的远程会议音频数据的远程会议音频数据混合;以及
基于所述远程会议音频数据混合来检测所述啸叫状态。
31.如EEE 30所述的远程会议服务器,其中,所述控制系统被配置用于:
经由所述接口系统从所述远程会议客户端位置接收远程会议元数据;以及
至少部分地基于所述远程会议元数据来确定哪个客户端位置正在导致所述啸叫状态。
32.如EEE 31所述的远程会议服务器,其中,所述远程会议元数据包括语音活动检测元数据、电平元数据或能量元数据中的一个或多个。
33.如EEE 26至32中任一项所述的远程会议服务器,其中,所述控制系统被配置用于根据从所述远程会议客户端位置中的一个或多个远程会议客户端位置接收的远程会议音频数据来估计啸叫存在概率,其中,检测所述啸叫状态至少部分基于所述啸叫存在概率。
34.如EEE 33所述的远程会议服务器,其中,所述啸叫存在概率估计是基于分层规则集或机器学习过程。
35.如EEE 26至34中任一项所述的远程会议服务器,其中,检测所述啸叫状态涉及从由以下各项构成的过程列表中选择的一个或多个过程:
根据频带的顺序统计量来计算基于功率的指标;计算频谱共振指标;计算百分位数间距指标;事件聚合;计算周期性指标;检测包络相似度;计算频谱峰度指标;基于机器学习的过程;以及估计啸叫存在概率。
36.如EEE 26至35中任一项所述的远程会议服务器,其中,所述控制系统被配置用于:
经由所述接口系统从某一远程会议客户端位置处的远程会议装置接收啸叫特征数据,所述啸叫特征数据是由所述远程会议装置执行的啸叫特征提取过程的结果;以及
至少部分地基于所述啸叫特征数据来检测所述啸叫状态。
37.一种远程会议客户端装置,包括:
接口系统,所述接口系统被配置用于所述远程会议客户端装置与远程会议服务器之间的通信,所述远程会议服务器被配置用于在远程会议期间在多个远程会议客户端位置之间提供全双工音频连接;以及
控制系统,所述控制系统被配置用于:
在所述远程会议期间检测啸叫状态,所述啸叫状态是涉及所述远程会议客户端装置和包括所述远程会议客户端装置的远程会议客户端位置中的至少一个其他装置的声反馈状态,其中,检测所述啸叫状态涉及对远程会议音频数据的频谱特性和时间特性两者的分析;以及
减轻所述啸叫状态。
38.如EEE 37所述的远程会议客户端装置,其中,检测所述啸叫状态涉及从由以下各项构成的过程列表中选择的一个或多个过程:根据频带的顺序统计量来计算基于功率的指标;计算频谱共振指标;计算百分位数间距指标;事件聚合;计算周期性指标;检测包络相似度;计算频谱峰度指标;基于机器学习的过程;以及估计啸叫存在概率。
39.一种远程会议客户端装置,包括:
接口系统,所述接口系统被配置用于所述远程会议客户端装置与远程会议服务器之间的通信,所述远程会议服务器被配置用于在远程会议期间在多个远程会议客户端位置之间提供全双工音频连接;以及
控制系统,所述控制系统被配置用于:
在所述远程会议期间执行啸叫状态特征提取过程,所述啸叫状态特征提取过程涉及提取与检测啸叫状态相关的啸叫特征数据,所述啸叫状态是涉及所述远程会议客户端装置和包括所述远程会议客户端装置的远程会议客户端位置中的至少一个其他装置的声反馈状态,其中,所述啸叫状态特征提取过程涉及对远程会议音频数据的频谱特性和时间特性两者的分析;以及
将所述啸叫特征数据发送到所述远程会议服务器。
40.如EEE 39所述的远程会议客户端装置,其中,执行啸叫状态特征提取过程涉及从由以下各项构成的过程列表中选择的一个或多个过程:根据频带的顺序统计量来计算基于功率的指标;计算频谱共振指标;计算百分位数间距指标;事件聚合;计算周期性指标;检测包络相似度;计算频谱峰度指标;基于机器学习的过程;以及估计啸叫存在概率。
Claims (16)
1.一种用于啸叫检测的通信设备,包括:
控制系统,所述控制系统被配置用于:
在位于两个或更多个位置的两个或更多个通信设备之间的通信会话期间检测啸叫状态,所述啸叫状态是涉及两个或更多个通信设备的声反馈状态,其中,检测所述啸叫状态包括:
从所述两个或更多个通信设备中一个或更多个通信设备接收音频数据;
将所述音频数据从时域变换到频域;
将经变换的音频数据划分成频带;
在经变换且经频带划分的音频数据中检测表示啸叫的结构特征的短期频谱特征或中期频谱-时间特征中的至少一个,其中短期是几十毫秒的数量级并且中期是几百毫秒的数量级,并且,其中短期频谱特征包括百分位数间距IPR,其中帧l的所述IPR被定义为IPR(l)=B(r,l)–B(k,l),其中B(r,l)是帧l的带功率的第r个顺序统计量版本并且k是频带数量的预定义部分;
估计啸叫存在概率并输出指示所述啸叫状态的存在的啸叫二元指示符;以及
确定所述两个或更多个通信设备中的哪个通信设备正在导致所述啸叫状态。
2.如权利要求1所述的通信设备,其中,短期频谱特征包括与集中在窄频率范围内的能量有关的频谱峰度度量SPM。
3.如权利要求1所述的通信设备,其中,k是所述频带数量的20%、18%或15%。
4.如权利要求1所述的通信设备,其中,在加权函数中使用输入电平和频谱峰度度量SPM对IPR(l)进行加权。
5.如权利要求1所述的通信设备,其中,所述控制系统被配置用于减轻所述啸叫状态或发送啸叫状态检测消息中的至少一个。
6.如权利要求1所述的通信设备,其中,所述控制系统被配置用于:
创建包括从两个或更多个通信设备接收的所述音频数据的音频数据混合;以及
基于所述音频数据混合来检测所述啸叫状态。
7.一种用于啸叫检测的方法,包括:
提取与检测啸叫状态相关的啸叫特征数据,所述啸叫状态是涉及彼此并置的两个或更多个通信设备的声反馈状态;
将音频数据从时域变换到频域;
将经变换的音频数据划分成频带;
在经变换且经频带划分的音频数据中检测表示啸叫的结构特征的短期频谱特征或中期频谱-时间特征中的至少一个,其中短期是几十毫秒的数量级并且中期是几百毫秒的数量级,并且其中短期频谱特征包括百分位数间距IPR,其中帧l的所述IPR被定义为IPR(l)=B(r,l)–B(k,l),其中B(r,l)是帧l的带功率的第r个顺序统计量版本并且k是频带数量的预定义部分。
8.如权利要求7所述的方法,其中,所述短期频谱特征以及中期频谱-时间特征进一步包括以下项中的至少一项:根据频带的顺序统计量的基于功率的指标;频谱共振指标;事件聚合;周期性指标;包络相似度;或频谱峰度指标。
9.一种用于啸叫检测的方法,包括:
在位于两个或更多个位置的两个或更多个通信设备之间的通信会话期间检测啸叫状态,所述啸叫状态是涉及两个或更多个通信设备的声反馈状态,其中,检测所述啸叫状态包括:
从两个或更多个通信设备接收音频数据;
将所述音频数据从时域变换到频域;
将经变换的音频数据划分成频带;
在经变换且经频带划分的音频数据中检测表示啸叫的结构特征的短期频谱特征或中期频谱-时间特征中的至少一个,其中短期是几十毫秒的数量级并且中期是几百毫秒的数量级,并且,其中短期频谱特征包括百分位数间距IPR,其中帧l的所述IPR被定义为IPR(l)=B(r,l)–B(k,l),其中B(r,l)是帧l的带功率的第r个顺序统计量版本并且k是频带数量的预定义部分;以及
估计啸叫存在概率并输出指示所述啸叫状态的存在的啸叫二元指示符;以及
确定所述两个或更多个通信设备中的哪个通信设备正在导致所述啸叫状态。
10.如权利要求9所述的方法,包括减轻所述啸叫状态。
11.如权利要求10所述的方法,其中减轻所述啸叫状态包括以下项中的至少一项:使正在导致所述啸叫状态的所述通信设备的一个或多个麦克风静音、或降低正在导致所述啸叫状态的所述通信设备的扬声器音量。
12.如权利要求11所述的方法,包括向所述通信设备发送啸叫检测消息。
13.一种用于啸叫检测的通信系统,包括:
通信设备,所述通信设备被配置为通过一个或多个处理器至少部分地基于从位于多个位置的多个通信设备接收到的音频数据检测啸叫状态,所述通信设备包括:
多个啸叫检测器,所述多个啸叫检测器中的每一个啸叫检测器被配置为基于对应的通信设备的音频数据的短期频谱特征或中期频谱-时间特征中的至少一个检测所述对应的通信设备是否存在啸叫状态,其中短期是几十毫秒的数量级并且中期是几百毫秒的数量级,并且,其中短期频谱特征包括百分位数间距IPR,其中帧l的所述IPR被定义为IPR(l)=B(r,l)–B(k,l),其中B(r,l)是帧l的带功率的第r个顺序统计量版本并且k是频带数量的预定义部分;
啸叫通信设备检测器,所述啸叫通信设备检测器被配置为根据检测到啸叫状态信号的啸叫检测器来确定两个或更多个通信设备中的哪个通信设备正在导致所述啸叫状态;以及
硬件和软件,所述硬件和软件被配置为执行以下项中的至少一项:减轻所述啸叫状态、或向正在导致所述啸叫状态的所述通信设备发送啸叫状态检测消息。
14.一种用于啸叫检测的通信系统,包括:
位于多个位置的多个通信设备,所述多个位置中的每一个位置对应于所述多个通信设备中各自的特定通信设备,每个特定通信设备被配置为通过一个或多个处理器基于对应的通信设备的音频数据的短期频谱特征或中期频谱-时间特征中的至少一个检测所述对应的通信设备是否存在啸叫状态,其中短期是几十毫秒的数量级并且中期是几百毫秒的数量级,并且,其中短期频谱特征包括百分位数间距IPR,其中帧l的所述IPR被定义为IPR(l)=B(r,l)–B(k,l),其中B(r,l)是帧l的带功率的第r个顺序统计量版本并且k是频带数量的预定义部分;以及
服务器,所述服务器包括:
啸叫通信设备检测器,所述啸叫通信设备检测器被配置为根据检测到啸叫状态信号的啸叫检测器来确定所述多个通信设备中的哪个通信设备正在导致所述啸叫状态;以及
硬件和软件,所述硬件和软件被配置为执行以下项中的至少一项:减轻所述啸叫状态、或向正在导致所述啸叫状态的所述通信设备发送啸叫状态检测消息。
15.一种计算机可读存储介质,所述计算机可读存储介质存储指令,所述指令在由一个或多个处理器执行时,使得所述一个或多个处理器执行根据权利要求7至12中任一项所述的方法。
16.一种计算机程序产品,包括计算机程序,所述计算机程序在由一个或多个处理器执行时,使得所述一个或多个处理器执行根据权利要求7至12中任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111051519.9A CN113766073B (zh) | 2017-09-29 | 2018-09-27 | 会议系统中的啸叫检测 |
Applications Claiming Priority (9)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNPCT/CN2017/104604 | 2017-09-29 | ||
CN2017104604 | 2017-09-29 | ||
US201762575701P | 2017-10-23 | 2017-10-23 | |
US62/575,701 | 2017-10-23 | ||
EP17201011 | 2017-11-10 | ||
EP17201011.8 | 2017-11-10 | ||
CN201880062458.6A CN111149370B (zh) | 2017-09-29 | 2018-09-27 | 会议系统中的啸叫检测 |
PCT/US2018/053119 WO2019067718A2 (en) | 2017-09-29 | 2018-09-27 | ACOUSTIC REACTION DETECTION IN CONFERENCE SYSTEMS |
CN202111051519.9A CN113766073B (zh) | 2017-09-29 | 2018-09-27 | 会议系统中的啸叫检测 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201880062458.6A Division CN111149370B (zh) | 2017-09-29 | 2018-09-27 | 会议系统中的啸叫检测 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113766073A CN113766073A (zh) | 2021-12-07 |
CN113766073B true CN113766073B (zh) | 2024-04-16 |
Family
ID=63794778
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111051519.9A Active CN113766073B (zh) | 2017-09-29 | 2018-09-27 | 会议系统中的啸叫检测 |
CN201880062458.6A Active CN111149370B (zh) | 2017-09-29 | 2018-09-27 | 会议系统中的啸叫检测 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201880062458.6A Active CN111149370B (zh) | 2017-09-29 | 2018-09-27 | 会议系统中的啸叫检测 |
Country Status (4)
Country | Link |
---|---|
US (2) | US11277518B2 (zh) |
EP (1) | EP3689002A2 (zh) |
CN (2) | CN113766073B (zh) |
WO (1) | WO2019067718A2 (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11277518B2 (en) * | 2017-09-29 | 2022-03-15 | Dolby Laboratories Licensing Corporation | Howl detection in conference systems |
EP4005230A1 (en) * | 2019-07-25 | 2022-06-01 | Unify Patente GmbH & Co. KG | Method and system for avoiding howling disturbance on conferences |
US11057702B1 (en) | 2019-12-20 | 2021-07-06 | Microsoft Technology Licensing, Llc | Method and system for reducing audio feedback |
GB2596169B (en) * | 2020-02-11 | 2022-04-27 | Tymphany Acoustic Tech Ltd | A method and an audio processing unit for detecting a tone |
CN111883163B (zh) * | 2020-04-02 | 2023-11-28 | 珠海市杰理科技股份有限公司 | 音频啸叫抑制方法、设备和系统及神经网络训练方法 |
CN112992190B (zh) * | 2021-02-02 | 2021-12-10 | 北京字跳网络技术有限公司 | 音频信号的处理方法、装置、电子设备和存储介质 |
US20220415336A1 (en) * | 2021-06-28 | 2022-12-29 | Panasonic Intellectual Property Management Co., Ltd. | Voice communication apparatus and howling detection method |
CN113724725B (zh) * | 2021-11-04 | 2022-01-18 | 北京百瑞互联技术有限公司 | 一种蓝牙音频啸叫检测抑制方法、装置、介质及蓝牙设备 |
CN113905310B (zh) * | 2021-12-09 | 2022-03-01 | 北京百瑞互联技术有限公司 | 一种蓝牙音频的啸叫检测和抑制方法、装置及介质 |
CN114420148B (zh) * | 2022-03-30 | 2022-06-14 | 北京百瑞互联技术有限公司 | 啸叫检测和抑制方法、装置、介质及设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1934897A (zh) * | 2004-02-20 | 2007-03-21 | 松下电器产业株式会社 | 啸叫检测方法及装置、以及具有它们的音响装置 |
WO2015160617A1 (en) * | 2014-04-16 | 2015-10-22 | Dolby Laboratories Licensing Corporation | Jitter buffer control based on monitoring of delay jitter and conversational dynamics |
CN105812993A (zh) * | 2014-12-29 | 2016-07-27 | 联芯科技有限公司 | 啸叫检测和抑制方法及其装置 |
CN105979197A (zh) * | 2016-07-19 | 2016-09-28 | 块互动(北京)科技有限公司 | 基于啸叫音自动识别的远程会议控制方法及装置 |
KR20170072783A (ko) * | 2015-12-16 | 2017-06-27 | 한국전자통신연구원 | 다자간 회의 서비스에서의 채널 상태 적응형 오디오 믹싱 방법 |
CN111149370A (zh) * | 2017-09-29 | 2020-05-12 | 杜比实验室特许公司 | 会议系统中的啸叫检测 |
Family Cites Families (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3654470B2 (ja) * | 1996-09-13 | 2005-06-02 | 日本電信電話株式会社 | サブバンド多チャネル音声通信会議用反響消去方法 |
WO2002054744A1 (en) * | 2000-12-29 | 2002-07-11 | Nokia Corporation | Audio signal quality enhancement in a digital network |
JP2004165888A (ja) * | 2002-11-12 | 2004-06-10 | Oki Electric Ind Co Ltd | ハウリング処理装置 |
GB0317158D0 (en) * | 2003-07-23 | 2003-08-27 | Mitel Networks Corp | A method to reduce acoustic coupling in audio conferencing systems |
TWI289020B (en) * | 2004-02-06 | 2007-10-21 | Fortemedia Inc | Apparatus and method of a dual microphone communication device applied for teleconference system |
US20050285935A1 (en) * | 2004-06-29 | 2005-12-29 | Octiv, Inc. | Personal conferencing node |
US20050286443A1 (en) * | 2004-06-29 | 2005-12-29 | Octiv, Inc. | Conferencing system |
US7742608B2 (en) * | 2005-03-31 | 2010-06-22 | Polycom, Inc. | Feedback elimination method and apparatus |
US7991146B2 (en) * | 2005-11-15 | 2011-08-02 | Mitel Networks Corporation | Anti-howling structure |
US7873175B2 (en) * | 2005-12-27 | 2011-01-18 | Polycom, Inc. | Multiplexed microphone signals with multiple signal processing paths |
JP2007300295A (ja) * | 2006-04-28 | 2007-11-15 | Matsushita Electric Ind Co Ltd | 会議マイクロフォン装置および会議マイクロフォン制御方法 |
US20080091415A1 (en) * | 2006-10-12 | 2008-04-17 | Schafer Ronald W | System and method for canceling acoustic echoes in audio-conference communication systems |
US8160239B2 (en) * | 2006-11-10 | 2012-04-17 | Sony Corporation | Echo canceller and speech processing apparatus |
JP4569618B2 (ja) * | 2006-11-10 | 2010-10-27 | ソニー株式会社 | エコーキャンセラ及び通話音声処理装置 |
JP2008263280A (ja) * | 2007-04-10 | 2008-10-30 | Iwatsu Electric Co Ltd | ハウリング防止装置 |
CN101511050A (zh) * | 2009-03-24 | 2009-08-19 | 北京中星微电子有限公司 | 助听装置及消除声反馈的方法 |
US8488745B2 (en) * | 2009-06-17 | 2013-07-16 | Microsoft Corporation | Endpoint echo detection |
US8462193B1 (en) | 2010-01-08 | 2013-06-11 | Polycom, Inc. | Method and system for processing audio signals |
CN102740214B (zh) * | 2011-04-01 | 2014-09-10 | 中国科学院声学研究所 | 一种基于反馈信号频谱估计的啸叫抑制方法 |
EP3917161B1 (en) * | 2012-11-12 | 2024-01-31 | Yamaha Corporation | Signal processing system and signal processing method |
US9521258B2 (en) * | 2012-11-21 | 2016-12-13 | Castel Communications, LLC | Real-time call center call monitoring and analysis |
DK2736271T3 (da) * | 2012-11-27 | 2019-09-16 | Oticon As | Fremgangsmåde til styring af en opdateringsalgoritme for et adaptivt feedbackestimeringssystem og en de-korrelationsenhed |
EP2975831A1 (en) | 2014-07-16 | 2016-01-20 | Dolby Laboratories Licensing Corporation | In-service monitoring of voice quality in teleconferencing |
US10142484B2 (en) * | 2015-02-09 | 2018-11-27 | Dolby Laboratories Licensing Corporation | Nearby talker obscuring, duplicate dialogue amelioration and automatic muting of acoustically proximate participants |
DK3139636T3 (da) * | 2015-09-07 | 2019-12-09 | Bernafon Ag | Høreanordning, der omfatter et tilbagekoblingsundertrykkelsessystem baseret på signalenergirelokation |
CN105516876A (zh) * | 2015-12-09 | 2016-04-20 | 天津大学 | 一种基于谱熵的啸叫检测方法 |
CN105788611A (zh) * | 2016-02-25 | 2016-07-20 | 成都普创通信技术股份有限公司 | 一种音频质量在线监测系统 |
DK3481085T3 (da) * | 2017-11-01 | 2020-11-30 | Oticon As | Tilbagekoblingsdetektor og en høreanordning, der omfatter en tilbagekoblingsdetektor |
-
2018
- 2018-09-27 US US16/640,169 patent/US11277518B2/en active Active
- 2018-09-27 CN CN202111051519.9A patent/CN113766073B/zh active Active
- 2018-09-27 WO PCT/US2018/053119 patent/WO2019067718A2/en unknown
- 2018-09-27 CN CN201880062458.6A patent/CN111149370B/zh active Active
- 2018-09-27 EP EP18783352.0A patent/EP3689002A2/en active Pending
-
2022
- 2022-03-10 US US17/691,966 patent/US11677879B2/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1934897A (zh) * | 2004-02-20 | 2007-03-21 | 松下电器产业株式会社 | 啸叫检测方法及装置、以及具有它们的音响装置 |
WO2015160617A1 (en) * | 2014-04-16 | 2015-10-22 | Dolby Laboratories Licensing Corporation | Jitter buffer control based on monitoring of delay jitter and conversational dynamics |
CN105812993A (zh) * | 2014-12-29 | 2016-07-27 | 联芯科技有限公司 | 啸叫检测和抑制方法及其装置 |
KR20170072783A (ko) * | 2015-12-16 | 2017-06-27 | 한국전자통신연구원 | 다자간 회의 서비스에서의 채널 상태 적응형 오디오 믹싱 방법 |
CN105979197A (zh) * | 2016-07-19 | 2016-09-28 | 块互动(北京)科技有限公司 | 基于啸叫音自动识别的远程会议控制方法及装置 |
CN111149370A (zh) * | 2017-09-29 | 2020-05-12 | 杜比实验室特许公司 | 会议系统中的啸叫检测 |
Also Published As
Publication number | Publication date |
---|---|
WO2019067718A3 (en) | 2019-05-09 |
US20210136230A1 (en) | 2021-05-06 |
US11277518B2 (en) | 2022-03-15 |
CN113766073A (zh) | 2021-12-07 |
WO2019067718A2 (en) | 2019-04-04 |
US20220201125A1 (en) | 2022-06-23 |
CN111149370B (zh) | 2021-10-01 |
EP3689002A2 (en) | 2020-08-05 |
US11677879B2 (en) | 2023-06-13 |
CN111149370A (zh) | 2020-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113766073B (zh) | 会议系统中的啸叫检测 | |
CN111048119B (zh) | 通话音频混音处理方法、装置、存储介质和计算机设备 | |
CN112071328B (zh) | 音频降噪 | |
CN105118522B (zh) | 噪声检测方法及装置 | |
CN110853664B (zh) | 评估语音增强算法性能的方法及装置、电子设备 | |
EP2896126B1 (en) | Long term monitoring of transmission and voice activity patterns for regulating gain control | |
CN101010722A (zh) | 音频信号中话音活动的检测 | |
EP2283484A1 (en) | System and method for dynamic sound delivery | |
CN111292758B (zh) | 语音活动检测方法及装置、可读存储介质 | |
CN115280414B (zh) | 基于期望信号的机器学习水平估计的自动增益控制 | |
US20230317096A1 (en) | Audio signal processing method and apparatus, electronic device, and storage medium | |
CN104580764B (zh) | 电话会议系统中的超声配对信号控制 | |
CN110148421B (zh) | 一种残余回声检测方法、终端和装置 | |
US11380312B1 (en) | Residual echo suppression for keyword detection | |
Diether et al. | Efficient blind estimation of subband reverberation time from speech in non-diffuse environments | |
JPH09311696A (ja) | 自動利得調整装置 | |
JP6942282B2 (ja) | 補助信号を用いたオーディオデバイスの送信制御 | |
CN113593599A (zh) | 一种去除语音信号中噪声信号的方法 | |
CN113593604A (zh) | 检测音频质量方法、装置及存储介质 | |
CN109218917B (zh) | 一种实时通信系统中的自动声反馈监测与消除方法 | |
RU2436173C1 (ru) | Способ обнаружения пауз в речевых сигналах и устройство его реализующее | |
KR20090082700A (ko) | 노이지 음성 신호의 처리 장치 및 그 장치를 포함하는 음성기반 어플리케이션 장치 | |
Zhang et al. | Fundamental frequency estimation combining air-conducted speech with bone-conducted speech in noisy environment | |
Gong et al. | Noise power spectral density matrix estimation based on modified IMCRA | |
WO2022068440A1 (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 | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 40065600 Country of ref document: HK |
|
GR01 | Patent grant | ||
GR01 | Patent grant |