CN103745726B - 一种自适应的变采样率音频采样方法 - Google Patents

一种自适应的变采样率音频采样方法 Download PDF

Info

Publication number
CN103745726B
CN103745726B CN201310562870.3A CN201310562870A CN103745726B CN 103745726 B CN103745726 B CN 103745726B CN 201310562870 A CN201310562870 A CN 201310562870A CN 103745726 B CN103745726 B CN 103745726B
Authority
CN
China
Prior art keywords
zero crossing
frequency
sampled
sample
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.)
Active
Application number
CN201310562870.3A
Other languages
English (en)
Other versions
CN103745726A (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.)
CLP Kesiyi Technology Co Ltd
Original Assignee
CETC 41 Institute
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 CETC 41 Institute filed Critical CETC 41 Institute
Priority to CN201310562870.3A priority Critical patent/CN103745726B/zh
Publication of CN103745726A publication Critical patent/CN103745726A/zh
Application granted granted Critical
Publication of CN103745726B publication Critical patent/CN103745726B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

本发明提供一种自适应的变采样率音频采样方法,其根据采样数据测量过零点数目,依据过零点数目调整采样频率,具体包括以下步骤:步骤1,通过采样点的电压值估计过零点的时刻:判断前后2个采样值的符号是否相同,若相同,则两个采样值的两个采样点间没有过零点发生;若不同,则两个采样值的两个采样点的电压穿过零点一次,两个采样点的电压符号发生一次变更;步骤2,统计一帧采样数据内变更的次数得出过零点的数目,具体而言,以采样数据的符号位为计数时钟,利用符号位的上升沿和下降沿,在一帧采样数据传输时间内进行计数,得到计数值就是这帧数据的过零点数目;步骤3,根据过零点数目与频率的对应关系,利用分段采样进行采样频率的选择。

Description

一种自适应的变采样率音频采样方法
技术领域
本发明属于信号处理技术领域,尤其涉及一种自适应的变采样率音频采样方法。
背景技术
无线电综合测试仪在测试音频信号时,为了测试音频信号失真度和信纳比,采样率是最高频率的10倍以上,以便采集最高音频的5次谐波。过高的采样率,导致在采集低频音频信号时,采集的数据不满一个周期,影响音频频率、音频电平音频失真度和信纳比等参数的精度。传统的解决办法舍弃低频音频信号或降低最高音频信号频率,该方法虽能解决最高频率音频信号5次谐波可测与最低频率音频信号采样数据满周期的矛盾,但是降低了音频测量范围;另一种方法是采用分频段采样,该方法无须压缩音频信号测试范围,但是测量前必须预先估计输入音频信号频率范围,不方便用户操作。
无线电综合测试仪的音频分析功能,通过时域采样,然后傅里叶变换获取所需要的参量进行相应的计算实现。进行傅立叶变换,在频域实现音频参数的测量。失真度的测量要求输入信号至少5次谐波可测;信号幅度的测量要求时域采样至少满一个周期。在音频信号频率范围较宽的情况下,固定采样频率无法同时满足高频音频信号5次谐波可测与低频音频信号采样数据满周期的要求。自适应的变采样率是根据输入音频信号的频率适时快速调整采样频率,同时满足以上两个条件。
但在现有技术中,存在着各种各样的问题。在音频分析中,高频音频信号5次谐波可测、低频音频信号采样数据满周期是一对矛盾,固定采样频率无法兼顾二者的需求。在工程实现中往往采用压缩音频频率范围,或者调整采样频率。压缩音频频率范围,降低仪器的测试能力,无法满足市场需求。通行的采样频率调整分频域、时域进行。频域调整就是变换采样率,对采样数据进行抽取、内插来实现采样率的调整,但前提是高频5次谐波和低频在数据中存在。时域调整是根据采样数据,通过测量过零点估计频率,然后根据估计频率调整采样频率,过零点检测要求输入信号时严格的正弦波形,工程中受噪声和检测方法影响,误差大。其次,傅立叶变换法估计频率;要求有足够的采样数据,在低频时不适用。
发明内容
为解决上述问题,本发明提供一种自适应的变采样率音频采样方法,其在不降低音频测试范围前提下,无须用户估计输入音频信号音频范围,自动实现采样率切换。本发明采用了检测过零点数目,依据过零点数目调整采样频率,有效兼顾了高频音频信号5次谐波可测、低频音频信号采样数据满周期的要求,同时避免了噪声的影响,降低了对采样频率切换时间的要求,扩展了音频频率测试范围。
本发明的一种自适应的变采样率音频采样方法根据采样数据测量过零点数目,依据所述过零点数目调整采样频率,具体包括以下步骤:
步骤1,通过采样点的电压值估计过零点的时刻:判断前后2个采样值的符号是否相同,若相同,则所述两个采样值的两个采样点间没有过零点发生;若不同,则所述两个采样值的两个采样点的电压穿过零点一次,所述两个采样点的电压符号发生一次变更;
步骤2,统计一帧采样数据内变更的次数得出过零点的数目,具体而言,以采样数据的符号位为计数时钟,利用符号位的上升沿和下降沿,在一帧采样数据传输时间内进行计数,得到计数值就是这帧数据的过零点数目;
步骤3,根据过零点数目与频率的对应关系,利用分段采样进行采样频率的选择;
进一步的,所述步骤3包括:
步骤31,同一频率内频率与过零点数目的对应关系如式(1)所示:
N为过零点数目,M帧数据长度,f为输入音频信号频率,fs为采样频率,不大于x的最大整数;
步骤32,采用分段采样减少采样频率级数,分段原则是:在该采样频率下,此段音频频率范围内,过零点数目范围为:[Nmin,Nmax],其中
步骤33,分段结束后,确定各频率段的最大频率在各采样率下的过零点数目,其中i大于等于1、j小于等于L,Nij表示第i频率范围,j大于等于1、j小于等于L,表示第j个采样频率,L表示所分的频率段个数,利用检测到的过零点数目与[N1j,N2j,……,NLj]进行采样频率的调整。
本发明的有益效果在于:
本发明利用采样数据的符号位进行计数,有效避免了过零点抖动和噪声的影响,减少了运算量。
本发明依据过零点数目与频率的对应关系,确定所需要采样率的级数;根据检测的过零点数目与计算出的过零点数目比较,自动选择采样频率。有效解决了高频音频信号5次谐波可测、低频音频信号采样数据满周期的矛盾,提高了测量精度,扩展了频率测量范围。
总而言之,本发明与目前广泛变采样率系统相比,节省了运算资源,提高采样频率转换速度,同时提高了测量精度,扩展了频率测量范围。
附图说明
图1是本发明的自适应的变采样率音频采样方法的正弦信号电压值采样示意图。
具体实施方式
本发明的基本原理是根据采样数据测量过零点数目,依据过零点的数量调整采样频率。如图1,一般情况下,过零点不是采样点。因此,只能通过采样点的电压值间接地估计过零点的时刻。过零点时刻的确定需要先判断前后2个采样值的符号是否相同,如前后符号相同,如图中m-2、m-1点,表明两采样点间没有过零点发生;如前后符号相反,如图中m-1、m点,表明m-1至m时间段内,电压穿过零点一次,采样点电压符号发生一次变更。
因此,可以用统计一帧采样数据内,采样数据变更的次数,从而得出过零点的数目。本方案中,以采样数据的符号位为计数时钟,利用符号位的上升沿和下降沿,在一帧数据传输时间内进行计数;得到计数值就是这帧数据的过零点数目。
采样频率的选择:
(1)过零点数目不是频率,不能直接确定采样频率。同一频率,不同采样频率,过零点数目也不相同,频率与过零点数目的对应关系如下:
N为过零点数目,M帧数据长度,f为输入音频信号频率,fs为采样频率,不大于x的最大整数。
(2)为减少采样频率级数,本方案采用分段采样。分段的原则是:在该采样频率下,此段音频频率范围内,过零点数目范围:[Nmin,Nmax],其中
(3)采样频率的选择。
分段结束后,确定各频率段的最大频率在各采样率下的过零点数目。以Nij表示,i大于等于1、j小于等于L,表示第i频率范围,j大于等于1、j小于等于L,表示第j个采样频率,L表示所分的频率段个数。利用检测到的过零点数目,与[N1j,N2j,……,NLj],进行采样频率的调整。
本发明的数理结果验证
在实验板上,输入幅度为1Vrms,频率为10Hz、100Hz、500Hz、1KHz、2KHz、10KHz、20KHz的音频信号,不同频率音频信号频率、失真度、信纳比测试结果如表1所示:
表1
输入信号 10Hz 100Hz 500Hz 1Kz 2Kz 10KHz 20KHz
频率(HZ) 10 100 500 1000 2000 10000 20000
失真度(%) 0.043 0.142 0.157 0.163 0.195 0.161 0.217
信纳比(dB) 42 55 46 51 44 55 51
本技术方案在时域上检测过零点数目,避免了过零点抖动与噪声的影响,对采样数据的长度也不像傅立叶变换法对数据长度有要求,能实时快速地自动调整采样频率。
当然,本发明还可有其他多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。

Claims (1)

1.一种自适应的变采样率音频采样方法,其特征在于,根据采样数据测量过零点数目,依据所述过零点数目调整采样频率,具体包括以下步骤:
步骤1,通过采样点的电压值估计过零点的时刻:判断前后2个采样值的符号是否相同,若相同,则所述两个采样值的两个采样点间没有过零点发生;若不同,则所述两个采样值的两个采样点的电压穿过零点一次,所述两个采样点的电压符号发生一次变更;
步骤2,统计一帧采样数据内变更的次数得出过零点的数目,具体而言,以采样数据的符号位为计数时钟,利用符号位的上升沿和下降沿,在一帧采样数据传输时间内进行计数,得到计数值就是这帧数据的过零点数目;
步骤3,根据过零点数目与频率的对应关系,利用分段采样进行采样频率的选择;
进一步的,所述步骤3包括:
步骤31,同一频率内频率与过零点数目的对应关系如式(1)所示:
N为过零点数目,M为帧数据长度,f为输入音频信号频率,fs为采样频率,表达式表示不大于x的最大整数;
步骤32,采用分段采样减少采样频率级数,分段原则是:在该采样频率下,此段音频频率范围内,过零点数目范围为:[Nmin,Nmax],其中
步骤33,分段结束后,确定各频率段的最大频率在各采样率下的过零点数目,其中i大于等于1、j小于等于L,Nij表示第i频率范围,j大于等于1、j小于等于L,表示第j个采样频率,L表示所分的频率段个数,利用检测到的过零点数目与[N1j,N2j,……,NLj]进行采样频率的调整。
CN201310562870.3A 2013-11-07 2013-11-07 一种自适应的变采样率音频采样方法 Active CN103745726B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310562870.3A CN103745726B (zh) 2013-11-07 2013-11-07 一种自适应的变采样率音频采样方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310562870.3A CN103745726B (zh) 2013-11-07 2013-11-07 一种自适应的变采样率音频采样方法

Publications (2)

Publication Number Publication Date
CN103745726A CN103745726A (zh) 2014-04-23
CN103745726B true CN103745726B (zh) 2016-08-17

Family

ID=50502739

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310562870.3A Active CN103745726B (zh) 2013-11-07 2013-11-07 一种自适应的变采样率音频采样方法

Country Status (1)

Country Link
CN (1) CN103745726B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107436244B (zh) * 2016-05-25 2021-03-19 上海金艺检测技术有限公司 基于频率分段振动数据采集的设备故障报警方法
CN109586997B (zh) * 2018-11-08 2021-02-02 北京车和家信息技术有限公司 采样方法、采样装置、计算机设备与可读存储介质
CN111341302B (zh) * 2020-03-02 2023-10-31 苏宁云计算有限公司 一种语音流采样率确定方法及装置
CN111354365B (zh) * 2020-03-10 2023-10-31 苏宁云计算有限公司 一种纯语音数据采样率识别方法、装置、系统
CN111462764B (zh) * 2020-06-22 2020-09-25 腾讯科技(深圳)有限公司 音频编码方法、装置、计算机可读存储介质及设备

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6163614A (en) * 1997-10-08 2000-12-19 Winbond Electronics Corp. Pitch shift apparatus and method
CN101878504A (zh) * 2007-08-27 2010-11-03 爱立信电话股份有限公司 使用时间分辨率能选择的低复杂性频谱分析/合成
CN102435844A (zh) * 2011-11-01 2012-05-02 南京磐能电力科技股份有限公司 一种频率无关的正弦信号相量计算方法
CN103067014A (zh) * 2012-12-25 2013-04-24 上海贝岭股份有限公司 用于谐波检测的频率跟随数据采集电路
CN202998069U (zh) * 2012-12-25 2013-06-12 上海贝岭股份有限公司 用于谐波检测的频率跟随数据采集电路

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6163614A (en) * 1997-10-08 2000-12-19 Winbond Electronics Corp. Pitch shift apparatus and method
CN101878504A (zh) * 2007-08-27 2010-11-03 爱立信电话股份有限公司 使用时间分辨率能选择的低复杂性频谱分析/合成
CN102435844A (zh) * 2011-11-01 2012-05-02 南京磐能电力科技股份有限公司 一种频率无关的正弦信号相量计算方法
CN103067014A (zh) * 2012-12-25 2013-04-24 上海贝岭股份有限公司 用于谐波检测的频率跟随数据采集电路
CN202998069U (zh) * 2012-12-25 2013-06-12 上海贝岭股份有限公司 用于谐波检测的频率跟随数据采集电路

Also Published As

Publication number Publication date
CN103745726A (zh) 2014-04-23

Similar Documents

Publication Publication Date Title
CN103745726B (zh) 一种自适应的变采样率音频采样方法
CN108831499B (zh) 利用语音存在概率的语音增强方法
CN104143341B (zh) 爆音检测方法和装置
AU2014386442B9 (en) Method for detecting audio signal and apparatus
CN105785324B (zh) 基于mgcstft的线性调频信号参数估计方法
CN104464722A (zh) 基于时域和频域的语音活性检测方法和设备
CN102606891A (zh) 漏水检测仪及检测系统、及检测方法
CN105508276A (zh) 一种城区变电站主变压器风扇噪声评价方法
CN105118511A (zh) 一种雷声识别方法
CN106504760B (zh) 宽频带背景噪声与语音分离检测系统及方法
CN103632682A (zh) 一种音频特征检测的方法
CN102176667B (zh) 一种基于频率跟踪的工频数字陷波装置与方法
CN103905656B (zh) 残留回声的检测方法及装置
CN109474355A (zh) 基于频谱监测数据的自适应噪声门限估计与信号提取方法
CN102998528B (zh) 动态正弦畸变信号中频率分量的提取方法及其实现装置
JP2013074305A (ja) 電磁干渉検出装置、電磁干渉検出方法および電磁干渉検出プログラム
CN116828354A (zh) 一种用于领夹无线麦克风的收音质量优化方法及系统
CN103280225B (zh) 一种低复杂度的静音检测方法
CN109903775A (zh) 一种音频爆音检测方法和装置
CN105791530B (zh) 输出音量调节方法和装置
CN104424954B (zh) 噪声估计方法与装置
CN105738698B (zh) 一种基于中心频移的谐波参数估计算法
CN102208190B (zh) 抑制非平稳噪声设备收敛时间的测量方法及装置
CN107561359A (zh) 一种基于布莱克曼窗的间谐波提取方法
CN203616430U (zh) 具备自适应抗干扰功能的低成本局部放电在线监测装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20190313

Address after: 266000 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province

Patentee after: China Electronics Technology Instrument and Meter Co., Ltd.

Address before: 266555 No. 98 Xiangjiang Road, Qingdao economic and Technological Development Zone, Shandong

Patentee before: The 41st Institute of CETC

CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: Huangdao Xiangjiang Road 266555 Shandong city of Qingdao Province, No. 98

Patentee after: CLP kesiyi Technology Co.,Ltd.

Address before: 266000 No. 98 Xiangjiang Road, Huangdao District, Shandong, Qingdao

Patentee before: CHINA ELECTRONIC TECHNOLOGY INSTRUMENTS Co.,Ltd.