CN116405823A - 一种骨传导耳机音频智能去噪增强方法 - Google Patents
一种骨传导耳机音频智能去噪增强方法 Download PDFInfo
- Publication number
- CN116405823A CN116405823A CN202310638865.XA CN202310638865A CN116405823A CN 116405823 A CN116405823 A CN 116405823A CN 202310638865 A CN202310638865 A CN 202310638865A CN 116405823 A CN116405823 A CN 116405823A
- Authority
- CN
- China
- Prior art keywords
- signal
- component
- peak
- complexity
- adaptive sampling
- 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
- 238000000034 method Methods 0.000 title claims abstract description 45
- 210000000988 bone and bone Anatomy 0.000 title claims abstract description 34
- 238000005070 sampling Methods 0.000 claims abstract description 126
- 238000001914 filtration Methods 0.000 claims abstract description 23
- 230000001186 cumulative effect Effects 0.000 claims abstract description 6
- 238000000354 decomposition reaction Methods 0.000 claims description 54
- 230000003044 adaptive effect Effects 0.000 claims description 45
- 230000002708 enhancing effect Effects 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 abstract description 14
- 230000005236 sound signal Effects 0.000 abstract description 8
- 230000000694 effects Effects 0.000 abstract description 7
- 239000004973 liquid crystal related substance Substances 0.000 description 16
- 238000004364 calculation method Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000008054 signal transmission Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 210000000613 ear canal Anatomy 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 210000003625 skull Anatomy 0.000 description 1
- 210000003454 tympanic membrane Anatomy 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R1/00—Details of transducers, loudspeakers or microphones
- H04R1/10—Earpieces; Attachments therefor ; Earphones; Monophonic headphones
- H04R1/1083—Reduction of ambient noise
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/16—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/175—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
- G10K11/178—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
- G10K11/1785—Methods, e.g. algorithms; Devices
- G10K11/17853—Methods, e.g. algorithms; Devices of the filter
- G10K11/17854—Methods, e.g. algorithms; Devices of the filter the filter being an adaptive filter
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
- G10L21/0216—Noise filtering characterised by the method used for estimating noise
- G10L21/0232—Processing in the frequency domain
-
- 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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R2201/00—Details of transducers, loudspeakers or microphones covered by H04R1/00 but not provided for in any of its subgroups
- H04R2201/10—Details of earpieces, attachments therefor, earphones or monophonic headphones covered by H04R1/10 but not provided for in any of its subgroups
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- Multimedia (AREA)
- Computational Linguistics (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
Abstract
本发明涉及音频处理处理领域,具体涉及一种骨传导耳机音频智能去噪增强方法,包括:获取参考信号;对参考信号进行分解,获取若干信号分量及若干单峰信号及单分量信号,进而获取每一单峰信号的分解信号集合,根据分解信号集合获取每一单峰信号的所有级分量,获取所有级分量中相邻两级分量的分解复杂程度,对分解复杂程度进行累乘处理,获取每一单峰信号的成分复杂度;根据成分复杂度获取自适应采样点个数,得到调整步长因子,获取每一自适应采样点处的更新滤波权重;利用ANC算法对参考信号进行去噪处理。本发明旨在解决使用ANC算法对骨传导耳机音频去噪时,音频信号在传输通道中的循环损失进而导致去噪效果较差的问题。
Description
技术领域
本发明涉及音频处理技术领域,具体涉及一种骨传导耳机音频智能去噪增强方法。
背景技术
骨传导耳机原理为声音由空气振动传入耳道再敲撞鼓膜,深入耳朵内里带动耳小骨一起振动;由于骨传感耳机通过皮肤、颅骨传递至耳小骨,其音质与传统耳机相比较差,需要对骨传导耳机音频进行去噪增强预处理。
由于骨传导耳机中的音频噪声大多属于外界噪声,现有通常使用ANC主动降噪技术通过采集外接噪音,并产生与噪音反相的信号用耳机等装置回放,用以抵消噪音;但外界噪声与误差传感器之间存在三条传播通路,信号在通道内循环转化和传输过程中会出现采样点丢失及时延问题,导致出现较高的FM(调频)使降噪效果较差;Morgan在LMS的基础上提出了FXLMS,在参考信号送入LMS算法模块前增加一个对次级通道传递函数的估计,从而解决次级通道所引起的相关延时问题,但并没有解决信号在传输通道中的循环损失问题。
发明内容
本发明提供一种骨传导耳机音频智能去噪增强方法,以解决现有的利用ANC主动降噪技术对骨传导耳机音频进行去噪增强时,降噪效果较差的问题。
本发明的一种骨传导耳机音频智能去噪增强方法采用如下技术方案:
本发明一个实施例提供了一种骨传导耳机音频智能去噪增强方法,该方法包括以下步骤:
获取参考信号;
对参考信号进行分解获取若干信号分量,利用参考信号及信号分量的波谷点进行分割,获取若干单峰信号及单分量信号;根据每一单峰信号与不同单分量信号峰值点的时序值的差异及时宽差异获取每一单峰信号的分解信号集合;
将单峰信号视为第零级分量;在单峰信号的分解信号集合中的第一个信号分量中,将时宽上包含第零级分量的单分量信号记为第一级分量,根据第零级分量和第一级分量获得第一个分解复杂程度;在分解信号集合中的第二个信号分量中,将时宽上包含第一级分量的单分量信号记为第二级分量,根据第一级分量和第二级分量获得第二个分解复杂程度,依次类推,获得所有分解复杂程度;对获取的所有分解复杂程度进行累乘处理,得到每一单峰信号的成分复杂度;
根据每一单峰信号的成分复杂度获取每一单峰信号的自适应采样点个数,根据每一自适应采样点与上一自适应采样点的采样间隔变化、每一自适应采样点上一自适应采样点的调整步长因子获取每一自适应采样点的调整步长因子,根据调整步长因子获取每一自适应采样点处的更新滤波权重;
根据自适应采样点个数与更新滤波权重利用ANC算法对参考信号进行去噪处理。
可选的,所述根据每一单峰信号与不同单分量信号峰值点的时序值的差异及时宽差异获取每一单峰信号的分解信号集合,包括的具体步骤如下:
记参考信号的任意一单峰信号为关注单峰信号,根据每一单分量信号与关注单峰信号的峰值时序值差异与时序宽度差异获取每一单分量信号的差异程度,获取差异程度最小值对应的单分量信号所在的信号分量,记为关注单峰信号的起始信号分量;将信号分量按照获取的先后顺序从小到大进行编号,将编号大于等于起始信号分量编号的所有信号分量记为关注单峰信号的信号分量集合;将关注单峰信号与关注单峰信号的信号分量集合共同构成一个新集合,记为关注单峰信号的分解信号集合。
可选的,所述根据每一单分量信号与关注单峰信号的峰值时序值差异与时序宽度差异获取每一单分量信号的差异程度,包括的具体步骤如下:
记参考信号的任意一单峰信号为关注单峰信号,获取每一单分量信号与关注单峰信号峰值点时序值的差值绝对值,记为每一单分量信号的时序差异,获取每一单分量信号与关注单峰信号在时序上宽度的差值绝对值,记为每一单分量信号的宽度差异,对每一单分量信号的时序差异的平方与宽度差异的平方进行求和,记为每一单分量信号的差异程度。
可选的,所述根据第零级分量和第一级分量获得第一个分解复杂程度,包括的具体步骤如下:
获取属于第一级分量时宽范围内、且在第零级分量所在信号分量上的所有单分量信号的峰值幅值总个数;获取峰值幅值类型个数/>,获取所有峰值幅值类型中每类峰值幅值的个数记为/>,其中/>,根据/>与/>的比值得到第零级分量和第一级分量获得的第一个分解复杂程度。
可选的,所述第零级分量所在信号分量是指单峰信号所在的参考信号。
可选的,所述对获取的所有分解复杂程度进行累乘处理,得到每一单峰信号的成分复杂度,包括的具体步骤如下:
可选的,所述根据每一单峰信号的成分复杂度获取每一单峰信号的自适应采样点个数,包括的具体步骤如下:
可选的,所述根据每一自适应采样点与上一自适应采样点的采样间隔变化、每一自适应采样点上一自适应采样点的调整步长因子获取每一自适应采样点的调整步长因子,包括的具体步骤如下:
其中,表示第/>个自适应采样点的调整步长因子,/>表示第/>个自适应采样点至第/>个自适应采样点的调整步长因子,/>表示参考信号的固定步长因子,/>表示参考信号上第/>个自适应采样点所在时刻,/>、/>分别表示参考信号第/>、/>个自适应采样点所在时刻,/>表示第/>个自适应采样点与第/>个自适应采样点的采样间隔,/>表示第/>个自适应采样点与第/>个自适应采样点的采样间隔。
可选的,所述获取参考信号,包括使用声学传感器获取环境噪声信号,记为参考信号。
可选的,所述对参考信号进行分解获取若干信号分量,包括对参考信号进行EMD分解,获取若干个信号分量。
可选的,所述利用参考信号及信号分量的波谷点进行分割,获取若干单峰信号及单分量信号,包括的具体步骤如下:
使用SG滤波算法对参考信号进行去噪平滑处理;
获取去噪后的参考信号的所有波谷点,对参考信号按照波谷点进行分割,获取若干个单峰信号;获取每一信号分量的所有波谷点,对每一信号分量按照波谷点进行分割,获取得到若干个单分量信号。
本发明的技术方案的有益效果是:本方法考虑使用ANC算法对骨传感耳机音频信号进行去噪的过程中,参考信号多次迭代循环会发生数据丢失进而产生FM调频的问题,首先对参考信号进行分解获取信号分量,进而对参考信号及信号分量进行划分,获取单峰信号及单分量信号,根据每一单峰信号其信号分解的变化情况获取每一单峰信号的成分复杂度,进而根据成分复杂度调节采样点数量、采样频率,使冗余度较低,对频率较大、易丢失的单峰信号增加较多的采样点数量,提高保真率,并考虑新增采样点输入自适应滤波器后无响应的问题,通过调节步长因子,获取新增采样点的滤波权重。本方法相比于传统ANC去噪算法而言,在原有FXLMS算法解决时延问题的基础上,优化了信号在传输通道中迭代循环的损失问题,进一步减少了出现FM的可能性,降低其对去噪精度的影响,保证参考信号的保真率,使误差信号迭代收敛更快,去噪效果更好。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明一个实施例提供的一种骨传导耳机音频智能去噪增强方法的步骤流程图;
图2为本发明一实施例所提供的参考信号;
图3、图4、图5、图6为参考信号的信号分量;
具体实施方式
为了更进一步阐述本发明为达成预定发明目的所采取的技术手段及功效,以下结合附图及较佳实施例,对依据本发明提出的一种骨传导耳机音频智能去噪增强方法,其具体实施方式、结构、特征及其功效,详细说明如下。在下述说明中,不同的“一个实施例”或“另一个实施例”指的不一定是同一实施例。此外,一或多个实施例中的特定特征、结构或特点可由任何合适形式组合。
除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。
下面结合附图具体的说明本发明所提供的一种骨传导耳机音频智能去噪增强方法的具体方案。
请参阅图1,其示出了本发明一个实施例提供的一种骨传导耳机音频智能去噪增强方法的步骤流程图,该方法包括以下步骤:
步骤S001、获取参考信号。
由于对于骨传感耳机而言,音频噪声即为外界任意声音信号,使用声学传感器获取环境噪声信号,记为参考信号;其中,声学传感器获取参考信号的步骤为现有技术,本实施例不再赘述。
步骤S002、获取若干信号分量,对参考信号进行平滑处理,根据参考信号及每一信号分量的波谷点获取若干单峰信号及单分量信号,获取每一单峰信号的分解信号集合,根据分解信号集合获取每一单峰信号的所有级分量,获取所有级分量中相邻两级分量的分解复杂程度,对分解复杂程度进行累乘处理,获取每一单峰信号的成分复杂度。
需要说明的是,由于使用ANC算法对参考信号进行去噪时,参考信号在通路中多次循环,每次传输过程中不可避免的传输通道波动会导致采样点丢失,进而导致信号损耗,使得参考信号频率与实际耳机音频信号中的噪声信号频率存在误差;由于信号的冗余度越高即局部重复信号越多,在传输过程中的传输效率越高、传输速率越快,其由于传输通道波动导致信号损失的可能性越小,但外界噪声信号成分复杂,不同频率信号互相叠加,不易区分和评估,因此首先需要对参考信号进行EMD分解,依次获得获取多个信号分量,越向下分解,信号频率越低,每个信号分量内信号频率几乎一致;由于对参考信号进行EMD分解后,参考信号中的各类信号分解在不同的分量信号中,对于局部信号损失概率进行分析时,需要找到参考信号中的每个各类信号初次分解至完全分解的过程,即对参考信号分割为单峰信号进行分析时,对于某一频率较大的单峰信号,可能在中不存在该单峰信号的信号分量,需要先找到该单峰信号首个分解时所在分量,再开始对向下的信号分量进行分析。
具体的,对参考信号进行EMD分解,获取若干个信号分量,其中,/>表示获取得到的信号分量个数;需要说明的是,对信号进行EMD分解获取信号分量的步骤为现有技术,本实施例不再赘述;根据信号分量获得的先后顺序,将这些信号分量从小到大进行编号:1,2,…,/>。
请参阅图2至图6;其中,图2为一参考信号;图3至图6为图2所示参考信号进行EMD分解获取的信号分量。
进一步需要说明的是,由于参考信号上存在幅值极低的小波动,这些波动不具有较大的噪声干扰,但是会对获取参考信号中的各类信号的过程产生干扰,使某些各类信号的峰值点不突出,因此设置阈值将弱的小噪声滤除。
使用SG滤波算法对参考信号进行去噪平滑处理,设定滤波窗口为,本实施例设定/>,具体实施过程实施者可根据实际情况设定窗口大小,为方便后续叙述,记去噪后的参考信号为参考信号,即后续步骤中出现的参考信号均指的是去噪后的参考信号。
获取参考信号的所有波谷点,对参考信号按照波谷点进行分割,获取得到若干个单峰信号,获取每一单峰信号的峰值点及该单峰信号在时序上的宽度;获取每一信号分量的波谷点,对每一信号分量按照波谷点进行分割,获取得到若干个单分量信号,获取每一单分量信号的峰值点及该单分量信号在时序上的宽度。
记参考信号的任意一单峰信号为关注单峰信号,获取所有单分量信号与关注单峰信号峰值点时序值的差值绝对值,记为每一单分量信号的时序差异,获取所有单分量信号与关注单峰信号在时序上宽度的差值绝对值,记为每一单分量信号的宽度差异,对每一单分量信号的时序差异的平方与宽度差异的平方进行求和,记为每一单分量信号的差异程度,获取差异程度最小值对应的单分量信号所在的信号分量,记为关注单峰信号的起始信号分量;将起始信号分量作为关注单信号分量的下一级信号分量,记自起始信号分量起始的所有信号分量为关注单峰信号的信号分量集合,也即编号大于等于起始信号分量编号的所有信号分量为关注单峰信号的信号分量集合。
上述获得了每个单峰信号的信号分量集合。
将每一单峰信号与每一单峰信号的信号分量集合共同构成一个新集合,记为每一单峰信号的分解信号集合,即每一单峰信号的信号分量集合中包含该单峰信号及该单峰信号的信号分量集合。
进一步需要说明的是,由于熵值可表征数据的混乱程度即复杂程度,可借助熵值的计算思想获取每一单峰信号的成分复杂程度。
以参考信号中的第个单峰信号为例,将第/>个单峰信号视为第/>个单峰信号的第零级分量,获取第/>个单峰信号的分解信号集合中的第1个信号分量,在第1个信号分量上,获取包含第/>个单峰信号时宽的单分量信号/>,记为第/>个单峰信号的下级分量/>,记/>为第/>个单峰信号的第一级分量,/>为/>至/>之间(包括/>和/>),其中,/>、/>分别为/>的左侧边界和右侧边界,即单分量信号/>的时宽为/>。
需要说明的是上述的所有信号分量都是基于参考信号分解得到的,本实施例中将参考信号也视可以为一个信号分量。
获取在至/>之间,第/>个单峰信号的第零级分量所在的参考信号中单峰信号的个数/>及每一单峰信号的峰值点的信号幅值,记为峰值幅值,由于每一单峰信号对应一个峰值点,/>也可表示/>至/>之间的峰值幅值的总个数,记为第/>个单峰信号的第零级分量的峰值幅值总个数;由于/>至/>之间包含若干个峰值幅值,这些峰值幅值中取值相同的归为一类,总共获取/>类,记/>为第/>个单峰信号的第零级分量的峰值幅值类型个数,其中第/>类包含了/>个峰值幅值,进而/>表示/>至/>之间第/>类峰值幅值个数与/>至之间峰值幅值总个数的比值。
的计算过程为熵值计算公式,若/>所在时宽在参考信号中单峰信号个数越多,单峰信号的信号幅值的种类越多,说明单分量信号/>越复杂,即若/>越大,说明单分量信号/>越复杂;/>记为得到的第/>个单峰信号的第一个分解复杂程度。
进一步需要说明的是,由于每一单峰信号首次被分解至完全被分解获取得到的所有级信号为该单峰信号的所有信号组成成分,所有单峰信号的信号组成成分在参考信号上相互叠加、相互影响,因此每一单峰信号的成分复杂度为单峰信号的所有级信号中相邻两个信号分量之间的分解复杂程度的累乘。
记任意一单峰信号为目标单峰信号,将目标单峰信号视为目标单峰信号的第零级分量,目标单峰信号的分解信号集合中的第一个信号分量中时宽上包含第零级分量的单分量信号记为目标单峰信号的第一级分量,分解信号集合中的第二个信号分量中时宽上包含第一级分量的单分量信号记为目标单峰信号的第二级分量,依次类推获得目标单峰信号的所有级分量。
按照上述根据每一单峰信号的第零级分量与第一级分量间的分解复杂程度获取方法,同理获取每一单峰信号的所有级分量中相邻两级分量之间的分解复杂程度。
至此,获取得到每一单峰信号的所有级分量中相邻两级分量之间的分解复杂程度。
请参阅图7,图7为在第个单峰信号的下级分量/>所在信号分量上,/>的下级分量所在时宽范围的单分量信号,其中每一单分量信号所对应的峰值幅值依次为/>、/>;在/>所在信号分量上,/>的下级分量/>所在时宽范围的单分量信号,其中每一单分量信号所对应的峰值幅值依次为/>、/> 、/> 、/> 、 、/> 、/> 、/>。
其中,表示第/>个单峰信号的第/>个分解复杂程度,/>表示第/>个单峰信号的分解复杂程度的个数;由于单峰信号向下分解时,每向下一级分解,所对应的分解复杂程度越小,为了避免出现熵值为0的情况,引入/>函数,即/>,使得熵值/>越大,/>越小,且/>不为0。由于/>函数转换了逻辑关系,因此对于第/>个单峰信号而言,其所对应的成分复杂度/>,应为/>;若/>越大,则对于第/>个单峰信号而言,其所在位置包含的信号分量越复杂,那么在传输通道中循环时,传输率更低,则更容易出现信号丢失;对求得的所有单峰信号的成分复杂度进行线性归一化处理,记为每一单峰信号对应的成分复杂度。
至此,获取得到参考信号中每一单峰信号的成分复杂度。
步骤S003、根据每一单峰信号的成分复杂度获取每一单峰信号的自适应采样点个数,进而获取两两自适应采样点的自适应采样间隔,根据相邻自适应采样间隔差异获取调整步长因子,进而根据调整步长因子获取每一自适应采样点处的更新滤波权重。
需要说明的是,对于单峰信号而言,成分复杂度越大,其在传输通道中循环时的传输率越低、出现信号丢失的可能性越大;由于音频信号均采用采样点的方式进行记录,采样点越多,音频信号曲线拟合的越平滑,保证率越好;但当采样点过多时会导致信号传输时间增加,进而也会增加信号丢失的风险,因此对于参考信号的每一单峰信号而言,根据每一单峰信号的成分复杂度调整采样点数量,获取每一单峰信号的自适应采样点个数进而更新采样频率,规避信号丢失风险。
其中,表示第/>个单峰信号的成分复杂度,/>表示第/>个单峰信号的原始采样点个数,由于采样点个数为整数,/>表示进行向下取整操作;若第/>个单峰信号的成分复杂度越大,第/>个单峰信号所在位置包含的信号分量越复杂,在传输通道中循环时,更容易出现信号丢失,需要增加的采样点越多,使得后续误差信号收敛效果更好,进而提高去噪效率;将获取得到的自适应采样点平均分布在第/>个单峰信号的时宽长度上,调整采样频率,需要说明的是,人耳可以听见的最大外界噪声基本在20000HZ左右,即用20000个采样点记录1秒钟的信号。利用pcm编码对参考信号传输时,默认采样频率也为20000HZ,进而每个单峰信号的采样点个数/>为/>。
至此,获取每一单峰信号的自适应采样点个数,进而调整采样频率。
进一步需要说明的是,由于对于单峰信号自适应更新后的采样频率,降低了丢失风险,但是对于同一个设备而言,其记录音频信号的采样频率是固定的,因此更新采样点后的参考信号在输入自适应滤波器时,新增的采样点在自适应滤波器中没有对应的响应点;参考信号进入滤波器后,利用FXLMS算法对每个单峰信号自适应滤波权重,然后得到次级噪声信号,这一过程中当单峰信号的采样点增加时,其采样间隔变小,而FXLMS算法的步长因子固定;FXLMS算法在迭代滤波权重时,其原本的固定步长因子可能会跳过新增采样点,不发生响应,因此需要结合每一单峰信号的自适应采样点个数对步长因子进行调整,进而根据调整后的步长因子对FXLMS算法的权重迭代公式进行更新,获取每一自适应采样点的更新滤波权重。
进一步的,记任意一自适应采样点为目标采样点,记时序上目标采样点的上一个自适应采样点为第一目标采样点,记时序上第一目标采样点的上一个自适应采样点为第二目标采样点,将第一目标采样点与目标采样点之间的步长因子记为目标采样点的步长因子。
其中,表示参考信号的固定步长因子,/>表示参考信号上第/>个自适应采样点所在时刻,/>、/>分别表示参考信号第/>、/>个自适应采样点所在时刻,/>表示第/>个自适应采样点和第/>个自适应采样点之间的采样间隔除以第/>个自适应采样点和第/>个自适应采样点之间的采样间隔;/>表示第/>个自适应采样点的调整步长因子;需要说明的是,固定步长因子为大于0小于参考信号的相关矩阵最大特征值的倒数,固定步长因子根据实际输入信号获取,获取方法为公知技术,本实施例不再赘述;当/>时,调整步长因子为参考信号的固定步长因子,即对于时序上第一采样点与第二采样点之间的步长因子不做调整,仍使用固定步长因子。
表示第/>个自适应采样点处的更新滤波权重,/>表示第/>个自适应采样点所在时刻所得的误差信号,/>为第/>个自适应采样点所在时刻在参考信号中所对应的信号,为FXLMS算法的初始滤波权重;需要说明的是,/>为FXLMS算法的权重迭代公式,为现有技术,本实施例不再赘述,式中/>、/>、/>均可根据FXLMS算法直接获取,本方法将权重迭代公式中的原始步长因子调整为调整步长因子,滤波权重调整为更新滤波权重,原始其他计算步骤不做改变。
至此,获取得到每一自适应采样点处的更新滤波权重。
步骤S004、根据自适应采样点个数与更新滤波权重利用ANC算法对参考信号进行去噪处理。
需要说明的是,本方法根据参考信号中每一单峰信号的成分复杂度获取自适应采样点个数,对采样点个数进行增加,进而根据自适应采样点对应的采样频率获取调整步长因子,更新滤波权重,使得新增的采样点在滤波器中获取滤波权重,进而使得冗余度较低、易发生丢失的局部信号保真度提高,提升去噪精度;使用调整后的滤波器对参考信号进行去噪处理。
具体的,利用自适应采样点个数与更新滤波权重将参考信号经过滤波器后输出为次级噪声信号,与耳机信号进行抵消,得到误差信号反馈回ANC控制器,调整次级噪声信号重新循环直至误差消失,完成ANC去噪过程,获取去噪后的骨传感耳机音频信号。
至此,完成对骨传导耳机音频的增强去噪。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (12)
1.一种骨传导耳机音频智能去噪增强方法,其特征在于,该方法包括以下步骤:
获取参考信号;
对参考信号进行分解获取若干信号分量,利用参考信号及信号分量的波谷点进行分割,获取若干单峰信号及单分量信号;根据每一单峰信号与不同单分量信号峰值点的时序值的差异及时宽差异获取每一单峰信号的分解信号集合;
将单峰信号视为第零级分量;在单峰信号的分解信号集合中的第一个信号分量中,将时宽上包含第零级分量的单分量信号记为第一级分量,根据第零级分量和第一级分量获得第一个分解复杂程度;在分解信号集合中的第二个信号分量中,将时宽上包含第一级分量的单分量信号记为第二级分量,根据第一级分量和第二级分量获得第二个分解复杂程度,依次类推,获得所有分解复杂程度;对获取的所有分解复杂程度进行累乘处理,得到每一单峰信号的成分复杂度;
根据每一单峰信号的成分复杂度获取每一单峰信号的自适应采样点个数,根据每一自适应采样点与上一自适应采样点的采样间隔变化、每一自适应采样点上一自适应采样点的调整步长因子获取每一自适应采样点的调整步长因子,根据调整步长因子获取每一自适应采样点处的更新滤波权重;
根据自适应采样点个数与更新滤波权重利用ANC算法对参考信号进行去噪处理。
2.根据权利要求1所述一种骨传导耳机音频智能去噪增强方法,其特征在于,所述根据每一单峰信号与不同单分量信号峰值点的时序值的差异及时宽差异获取每一单峰信号的分解信号集合,包括的具体步骤如下:
记参考信号的任意一单峰信号为关注单峰信号,根据每一单分量信号与关注单峰信号的峰值时序值差异与时序宽度差异获取每一单分量信号的差异程度,获取差异程度最小值对应的单分量信号所在的信号分量,记为关注单峰信号的起始信号分量;将信号分量按照获取的先后顺序从小到大进行编号,将编号大于等于起始信号分量编号的所有信号分量记为关注单峰信号的信号分量集合;将关注单峰信号与关注单峰信号的信号分量集合共同构成一个新集合,记为关注单峰信号的分解信号集合。
3.根据权利要求2所述一种骨传导耳机音频智能去噪增强方法,其特征在于,所述根据每一单分量信号与关注单峰信号的峰值时序值差异与时序宽度差异获取每一单分量信号的差异程度,包括的具体步骤如下:
记参考信号的任意一单峰信号为关注单峰信号,获取每一单分量信号与关注单峰信号峰值点时序值的差值绝对值,记为每一单分量信号的时序差异,获取每一单分量信号与关注单峰信号在时序上宽度的差值绝对值,记为每一单分量信号的宽度差异,对每一单分量信号的时序差异的平方与宽度差异的平方进行求和,记为每一单分量信号的差异程度。
6.根据权利要求4所述一种骨传导耳机音频智能去噪增强方法,其特征在于,所述第零级分量所在信号分量是指单峰信号所在的参考信号。
9.根据权利要求1所述一种骨传导耳机音频智能去噪增强方法,其特征在于,所述根据每一自适应采样点与上一自适应采样点的采样间隔变化、每一自适应采样点上一自适应采样点的调整步长因子获取每一自适应采样点的调整步长因子,包括的具体步骤如下:
10.根据权利要求1所述一种骨传导耳机音频智能去噪增强方法,其特征在于,所述获取参考信号,包括使用声学传感器获取环境噪声信号,记为参考信号。
11.根据权利要求1所述一种骨传导耳机音频智能去噪增强方法,其特征在于,所述对参考信号进行分解获取若干信号分量,包括对参考信号进行EMD分解,获取若干个信号分量。
12.根据权利要求1所述一种骨传导耳机音频智能去噪增强方法,其特征在于,所述利用参考信号及信号分量的波谷点进行分割,获取若干单峰信号及单分量信号,包括的具体步骤如下:
使用SG滤波算法对参考信号进行去噪平滑处理;
获取去噪后的参考信号的所有波谷点,对参考信号按照波谷点进行分割,获取若干个单峰信号;获取每一信号分量的所有波谷点,对每一信号分量按照波谷点进行分割,获取得到若干个单分量信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310638865.XA CN116405823B (zh) | 2023-06-01 | 2023-06-01 | 一种骨传导耳机音频智能去噪增强方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310638865.XA CN116405823B (zh) | 2023-06-01 | 2023-06-01 | 一种骨传导耳机音频智能去噪增强方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116405823A true CN116405823A (zh) | 2023-07-07 |
CN116405823B CN116405823B (zh) | 2023-08-29 |
Family
ID=87014530
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310638865.XA Active CN116405823B (zh) | 2023-06-01 | 2023-06-01 | 一种骨传导耳机音频智能去噪增强方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116405823B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117093827A (zh) * | 2023-10-16 | 2023-11-21 | 欣灵电气股份有限公司 | 基于物联网的智慧消防给水数据处理系统 |
CN117191305A (zh) * | 2023-11-06 | 2023-12-08 | 临沂市公路事业发展中心兰陵县中心 | 一种公路桥梁的状态评估方法及系统 |
Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080144853A1 (en) * | 2006-12-06 | 2008-06-19 | Sommerfeldt Scott D | Secondary Path Modeling for Active Noise Control |
US20120308025A1 (en) * | 2011-06-03 | 2012-12-06 | Hendrix Jon D | Adaptive noise canceling architecture for a personal audio device |
KR101442700B1 (ko) * | 2013-09-03 | 2014-09-23 | 서강대학교산학협력단 | 독립성분분석 적응필터를 이용한 헤드폰 유입잡음 제거방법 및 장치 |
JP2015152578A (ja) * | 2014-02-19 | 2015-08-24 | 株式会社島津製作所 | マススペクトルデータ処理装置及びマススペクトルデータ処理方法 |
US20150243271A1 (en) * | 2014-02-22 | 2015-08-27 | Apple Inc. | Active noise control with compensation for acoustic leak in personal listening devices |
CN109741727A (zh) * | 2019-01-07 | 2019-05-10 | 哈尔滨工业大学(深圳) | 基于有源噪声控制算法的主动降噪耳机、降噪方法及存储介质 |
CN110706686A (zh) * | 2019-12-13 | 2020-01-17 | 恒玄科技(北京)有限公司 | 降噪方法、自适应滤波器、入耳式耳机和半入耳式耳机 |
WO2020028280A1 (en) * | 2018-08-02 | 2020-02-06 | Dolby Laboratories Licensing Corporation | Auto-calibration of an active noise control system |
CN111128111A (zh) * | 2019-12-18 | 2020-05-08 | 清华大学苏州汽车研究院(相城) | 发动机主动噪声控制的变步长前馈控制系统及控制方法 |
CN111193977A (zh) * | 2019-12-13 | 2020-05-22 | 恒玄科技(上海)股份有限公司 | 耳机的降噪方法、自适应fir滤波器、除噪滤波器组和耳机 |
CN111627415A (zh) * | 2020-04-28 | 2020-09-04 | 重庆邮电大学 | 一种基于自适应MFxLMS算法的主动降噪装置及FPGA实现 |
US20200342889A1 (en) * | 2019-04-26 | 2020-10-29 | Beijing Didi Infinity Technology And Development Co., Ltd. | Systems and methods for enhancing audio signals |
CN113611279A (zh) * | 2021-09-15 | 2021-11-05 | 深圳沃迪声科技股份有限公司 | 一种主动降噪系统及降噪耳机 |
CN113676803A (zh) * | 2020-05-14 | 2021-11-19 | 华为技术有限公司 | 一种主动降噪方法及装置 |
CN114501211A (zh) * | 2020-11-13 | 2022-05-13 | 炬芯科技股份有限公司 | 具有丽音通透性的主动降噪电路、方法、设备及存储介质 |
CN115798502A (zh) * | 2023-01-29 | 2023-03-14 | 深圳市深羽电子科技有限公司 | 一种用于蓝牙耳机的音频去噪方法 |
CN115884032A (zh) * | 2023-02-20 | 2023-03-31 | 深圳市九音科技有限公司 | 一种后馈式耳机的智慧通话降噪方法及系统 |
CN116095553A (zh) * | 2022-12-23 | 2023-05-09 | 恒玄科技(上海)股份有限公司 | 一种用于耳机的anc系统、降噪方法及存储介质 |
-
2023
- 2023-06-01 CN CN202310638865.XA patent/CN116405823B/zh active Active
Patent Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080144853A1 (en) * | 2006-12-06 | 2008-06-19 | Sommerfeldt Scott D | Secondary Path Modeling for Active Noise Control |
US20120308025A1 (en) * | 2011-06-03 | 2012-12-06 | Hendrix Jon D | Adaptive noise canceling architecture for a personal audio device |
KR101442700B1 (ko) * | 2013-09-03 | 2014-09-23 | 서강대학교산학협력단 | 독립성분분석 적응필터를 이용한 헤드폰 유입잡음 제거방법 및 장치 |
JP2015152578A (ja) * | 2014-02-19 | 2015-08-24 | 株式会社島津製作所 | マススペクトルデータ処理装置及びマススペクトルデータ処理方法 |
US20150243271A1 (en) * | 2014-02-22 | 2015-08-27 | Apple Inc. | Active noise control with compensation for acoustic leak in personal listening devices |
WO2020028280A1 (en) * | 2018-08-02 | 2020-02-06 | Dolby Laboratories Licensing Corporation | Auto-calibration of an active noise control system |
CN109741727A (zh) * | 2019-01-07 | 2019-05-10 | 哈尔滨工业大学(深圳) | 基于有源噪声控制算法的主动降噪耳机、降噪方法及存储介质 |
US20200342889A1 (en) * | 2019-04-26 | 2020-10-29 | Beijing Didi Infinity Technology And Development Co., Ltd. | Systems and methods for enhancing audio signals |
CN110706686A (zh) * | 2019-12-13 | 2020-01-17 | 恒玄科技(北京)有限公司 | 降噪方法、自适应滤波器、入耳式耳机和半入耳式耳机 |
CN111193977A (zh) * | 2019-12-13 | 2020-05-22 | 恒玄科技(上海)股份有限公司 | 耳机的降噪方法、自适应fir滤波器、除噪滤波器组和耳机 |
CN111128111A (zh) * | 2019-12-18 | 2020-05-08 | 清华大学苏州汽车研究院(相城) | 发动机主动噪声控制的变步长前馈控制系统及控制方法 |
CN111627415A (zh) * | 2020-04-28 | 2020-09-04 | 重庆邮电大学 | 一种基于自适应MFxLMS算法的主动降噪装置及FPGA实现 |
CN113676803A (zh) * | 2020-05-14 | 2021-11-19 | 华为技术有限公司 | 一种主动降噪方法及装置 |
CN114501211A (zh) * | 2020-11-13 | 2022-05-13 | 炬芯科技股份有限公司 | 具有丽音通透性的主动降噪电路、方法、设备及存储介质 |
CN113611279A (zh) * | 2021-09-15 | 2021-11-05 | 深圳沃迪声科技股份有限公司 | 一种主动降噪系统及降噪耳机 |
CN116095553A (zh) * | 2022-12-23 | 2023-05-09 | 恒玄科技(上海)股份有限公司 | 一种用于耳机的anc系统、降噪方法及存储介质 |
CN115798502A (zh) * | 2023-01-29 | 2023-03-14 | 深圳市深羽电子科技有限公司 | 一种用于蓝牙耳机的音频去噪方法 |
CN115884032A (zh) * | 2023-02-20 | 2023-03-31 | 深圳市九音科技有限公司 | 一种后馈式耳机的智慧通话降噪方法及系统 |
Non-Patent Citations (1)
Title |
---|
郇战;戴永惠;: "基于预测滤波器的数字降噪耳机研究", 测控技术, no. 10 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117093827A (zh) * | 2023-10-16 | 2023-11-21 | 欣灵电气股份有限公司 | 基于物联网的智慧消防给水数据处理系统 |
CN117093827B (zh) * | 2023-10-16 | 2024-01-30 | 欣灵电气股份有限公司 | 基于物联网的智慧消防给水数据处理系统 |
CN117191305A (zh) * | 2023-11-06 | 2023-12-08 | 临沂市公路事业发展中心兰陵县中心 | 一种公路桥梁的状态评估方法及系统 |
CN117191305B (zh) * | 2023-11-06 | 2024-02-02 | 临沂市公路事业发展中心兰陵县中心 | 一种公路桥梁的状态评估方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN116405823B (zh) | 2023-08-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116405823B (zh) | 一种骨传导耳机音频智能去噪增强方法 | |
JP5718251B2 (ja) | 分解されたオーディオ信号の再構成のためのシステムおよび方法 | |
CN109493877B (zh) | 一种助听装置的语音增强方法和装置 | |
CN112562627B (zh) | 前馈滤波器设计方法、主动降噪方法、系统及电子设备 | |
CN110473567A (zh) | 基于深度神经网络的音频处理方法、装置及存储介质 | |
US10755728B1 (en) | Multichannel noise cancellation using frequency domain spectrum masking | |
CN110611871A (zh) | 一种数字助听器的啸叫抑制方法、系统及专用dsp | |
CN110111769B (zh) | 一种电子耳蜗控制方法、装置、可读存储介质及电子耳蜗 | |
CN112565981B (zh) | 啸叫抑制方法、装置、助听器及存储介质 | |
Pu et al. | A joint auditory attention decoding and adaptive binaural beamforming algorithm for hearing devices | |
US20140067384A1 (en) | Method and apparatus for canceling vocal signal from audio signal | |
CN110115049A (zh) | 基于记录对象声音的声音信号建模 | |
TWI459381B (zh) | 語音增強方法 | |
Esra et al. | Speech Separation Methodology for Hearing Aid. | |
CN112802490B (zh) | 一种基于传声器阵列的波束形成方法和装置 | |
US11849283B2 (en) | Mitigating acoustic feedback in hearing aids with frequency warping by all-pass networks | |
CN111028857B (zh) | 基于深度学习的多通道音视频会议降噪的方法及系统 | |
CN110992969B (zh) | 一种电子耳蜗的滤波器组配置方法及装置 | |
CN108877829B (zh) | 一种信号处理的方法和装置 | |
KR101936242B1 (ko) | 잡음 제거 장치 및 방법, 기록 매체 | |
CN106601249B (zh) | 一种基于听觉感知特性的数字语音实时分解/合成方法 | |
US20230262390A1 (en) | Audio denoising method and system | |
US7808335B2 (en) | Circuit arrangement and signal processing device | |
CN113689869A (zh) | 语音增强方法、电子设备以及计算机可读存储介质 | |
Shrivastav et al. | An optimized frequency response masking reconfigurable filter to enhance the performance of the hearing aid system |
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 |