CN108199715A - 数字域正弦波检测方法及检测装置 - Google Patents

数字域正弦波检测方法及检测装置 Download PDF

Info

Publication number
CN108199715A
CN108199715A CN201810006884.XA CN201810006884A CN108199715A CN 108199715 A CN108199715 A CN 108199715A CN 201810006884 A CN201810006884 A CN 201810006884A CN 108199715 A CN108199715 A CN 108199715A
Authority
CN
China
Prior art keywords
sine wave
sequence
len
list entries
input
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
CN201810006884.XA
Other languages
English (en)
Other versions
CN108199715B (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.)
Hunan Shuangchou Network Technology Co ltd
Original Assignee
Chongqing University of Post and Telecommunications
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 Chongqing University of Post and Telecommunications filed Critical Chongqing University of Post and Telecommunications
Priority to CN201810006884.XA priority Critical patent/CN108199715B/zh
Publication of CN108199715A publication Critical patent/CN108199715A/zh
Application granted granted Critical
Publication of CN108199715B publication Critical patent/CN108199715B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M1/00Analogue/digital conversion; Digital/analogue conversion
    • H03M1/12Analogue/digital converters

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Measuring Frequencies, Analyzing Spectra (AREA)

Abstract

本发明涉及信号处理技术领域,特别涉及数字域正弦波检测方法及检测装置,所述方法包括:设置突发有效标志、计数值为0,设定最小有效长度、最大测试序列长度;在正弦序列判决统计模块的A端输入第一输入序列An,B端输入为第二输入序列Bn;利用实正弦波序列和复正弦波序列逐点计算相角θn;利用正弦序列判决统计模块判断正弦波;判断正弦波属于突发正弦波还是持续正弦波,并计算出正弦波的频率;本发明本发明根据三角函数特性进行正弦波判断,存储空间小,计算复杂度低,对正弦波畸变的敏感度高。

Description

数字域正弦波检测方法及检测装置
技术领域
本法明涉及信号处理技术领域,特别涉及一种数字域正弦波检测方法及检测装置。
背景技术
正弦波的检测可以分为模拟和数字两大类,数字域正弦波检测的方法比较多,主要有DFT法,相关检测方法等。傅里叶变换通过把时域信号变换到频域,实现对正弦波信号的识别及频率计算,正弦波信号经过傅里叶变换后,会形成单一的谱线。利用这一特性,可以对正弦波进行识别及频率计算。
在数字域中,采用离散傅里叶变换(discrete Fourier transform,DFT),DFT 计算量大,快速傅里叶变换(fast Fourier transform,FFT)是DFT的一种快速的计算方法,使DFT的计算效率更高,速度更快。
但是DFT、FFT的运算量随着序列长度的增加而增加,复杂度变高;并且 DFT运算的频率分辨率为fs/N,若运算区间取短了则正弦波频率计算的精度不够,如果要保证精度,运算区间就必须取得比较大,这样造成对瞬态畸变不敏感,难以满足正弦波信号的检测精度。
发明内容
针对上述问题,本发明提出一种数字域正弦波检测方法及检测装置。
本发明的一种种数字域正弦波检测方法,如图1,包括:
S1、设置突发有效标志Flag_BurstValid为0、计数值Counter_Len为0,设定最小有效长度为Len_Validmin、最大测试序列长度为Len_Testmax;
S2、在正弦序列判决统计模块的A端输入第一输入序列An,B端输入第二输入序列Bn
S3、利用第一输入序列An和第二输入序列Bn逐点计算相角θn
S4、判断相邻两点的相角差是否满足|θn-1n|<εth,若不满足,令计数值 Counter_Len=0,判断Flag_BurstValid==1是否成立,若成立则转到步骤S5,若不成立则令Counter_Len=0并返回步骤S2;
若满足,则令Counter_Len+1,然后判断Counter_Len>Len_Validmin成立是否成立,若成立则令Flag_BurstValid=1,并返回步骤S2;若不成立,则是否 Counter_Len≥Len_Testmax,若是,跳至步骤S5,否则返回步骤S2;
S5、判断正弦波属于突发正弦波还是持续正弦波,并根据根据相角θn和计数值Counter_Len的值,计算出正弦波的频率,返回步骤S1;
其中,εth表示误差门限。
优选的,第一输入序列An和第二输入序列Bn的计算包括:
对于实正弦波序列,第一输入序列为An=X(n),第二输入序列为 Bn=X(n-Dm)+X(n+Dm);
对于复正弦波序列,第一输入序列为An=Xr(n)·Xr(n-Dm)+Xc(n)·Xc(n-Dm),第二输入序列为Bn=Xc(n)·Xr(n-Dm)-Xr(n)·Xc(n-Dm);
其中,X(n)表示实正弦波序列,Dm表示延迟单元的延迟参数,Xr(n)表示复正弦波序列的I路,Xc表示复正弦波序列的Q路。
本发明的一种数字域正弦波检测装置,包括:
数字域正弦波检测装置包括实正弦波序列输入模块、复正弦波序列输入模块和正弦序列判决统计模块,实正弦输入模块包括带通滤波器、延迟单元和加法器,复正弦输入模块包括带通滤波器、延迟单元、乘法器和加法器,正弦序列判决统计模块的输入端包括A端和B端,分别在A端输入第一输入序列An、B端输入第二输入序列Bn
实正弦波序列输入模块用于计算实正弦波序列A端的第一输入序列An和B 端的第二输入序列Bn
复正弦波序列输入模块用于计算复正弦波序列A端的第一输入序列An和B 端的第二输入序列Bn
正弦序列判决统计模块用于判断相邻两点的相角θn的绝对差值与误差门限εth的关系,并计算出正弦波的频率。
本发明实现复杂度低,每一个点的运算量远低于FFT的运算量,前面采样点计算过的关于θn的结果可以重复使用;对畸变的敏感度高,正弦波序列的相位满足严格的数学关系,因此本发明能实时的检测出正弦波序列的畸变,并且本发明存储空间小,在高信噪比条件下相对于目前已有的方法,对正弦波的检测效率和准确性都更高。
附图说明
图1为本发明数字域正弦波检测方法流程示意图;
图2为本发明实正弦波序列检测模块示意图;
图3为本发明复正弦波序列检测模块示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。
本发明的一种数字域正弦波检测方法,如图1,包括:
S1、初始化,设置突发有效标志Flag_BurstValid为0、计数值Counter_Len 为0,设定最小有效长度为Len_Validmin、最大测试序列长度为Len_Testmax;
S2、在正弦序列判决统计模块的A端输入第一输入序列An,B端输入为第二输入序列Bn
S3、利用第一输入序列An和第二输入序列Bn逐点计算相角θn
S4、判断相邻两点的相角差是否满足|θn-1n|<εth,若不满足,令计数值 Counter_Len=0,判断Flag_BurstValid==1是否成立,若成立则转到步骤S5,若不成立则令Counter_Len=0并返回步骤S2;
若满足,则令Counter_Len自增1(Counter_Len++),然后判断 Counter_Len>Len_Validmin成立是否成立,若成立则令Flag_BurstValid=1,并返回步骤S2;若不成立,则是否Counter_Len≥Len_Testmax,若是,进入步骤S5,否则返回步骤S2;
S5、判断正弦波属于突发正弦波还是持续正弦波,并根据相角θn和计数值Counter_Len的值,计算出正弦波的频率,返回步骤S1;
其中,εth表示误差门限,最小有效长度Len_Validmin用于判断为正弦波序列的最小长度门限,优选的Len_Validmin=5;
最大测试序列长度Len_Testmax,对持续正弦波序列进行检测时的最大检测长度。
本发明的实正弦波序列检测模块如图2,包括带通滤波器、延迟单元和加法器,实正弦波序列经过带通滤波器后得到X(n+Dm),X(n+Dm)通过延迟单元得到X(n),将X(n)输入到正弦序列判决统计模块的A端;将X(n)通过延迟单元得到X(n-Dm),将X(n-Dm)与X(n+Dm)输入加法器所得到的结果输入到正弦序列判决统计模块的B端,Dm为延迟单元的延迟参数。
本发明的复正弦波序列检测模块如图3,复正弦波序列分为I、Q两路,复正弦波序列的I路经过带通滤波器后得到Xr(n),Xr(n)通过延迟单元得到 Xr(n-Dm),复正弦波序列的Q路经过带通滤波器后得到Xc(n),Xc(n)通过延迟单元得到Xc(n-Dm),将Xr(n)·Xr(n-Dm)+Xc(n)·Xc(n-Dm)的结果输入到正弦序列判决统计模块的A端,将Xc(n)·Xr(n-Dm)-Xr(n)·Xc(n-Dm)的结果输入到正弦序列判决统计模块的B端。
正弦序列判决统计模块进行逐点运算,首先需要初始化系统,令突发有效标志Flag_BurstValid、计数值Counter_Len为0,设定最小有效长度为 Len_Validmin、最大测试序列长度为Len_Testmax,其中最小有效长度和最大测试序列长度可以根据实际应用场景由用户设定。
在正弦序列判决统计模块的A端输入第一输入序列An,B端输入为第二输入序列Bn,对于实正弦波序列,第一输入序列为An=X(n),第二输入序列为 Bn=X(n-Dm)+X(n+Dm);
对于复正弦波序列,第一输入序列为An=Xr(n)·Xr(n-Dm)+Xc(n)·Xc(n-Dm),第二输入序列为Bn=Xc(n)·Xr(n-Dm)-Xr(n)·Xc(n-Dm)。
对于实正弦波序列,相角为θn=arccosBn2An,对反三角函数的计算可以采用级数展开的方法,这里采用泰勒级数的方法展开,令X=Bn2An,可以得到:
对于复正弦波序列,相角为θn=arctanBn/An对反三角函数的计算可以采用级数展开的方法,这里采用泰勒级数的方法展开,令X=Bn/An,可以得到:
当计算出相角θn后,需要判断相邻两点的相角差是否满足|θn-1n|<εth,若不满足,令计数值Counter_Len=0,判断Flag_BurstValid==1是否成立,若成立则转到步骤S5,若不成立则令Counter_Len=0并返回步骤S2;
若满足,则令Counter_Len+1,然后判断Counter_Len>Len_Validmin成立是否成立,若成立则令Flag_BurstValid=1,并返回步骤S2;若不成立,则是否Counter_Len≥Len_Testmax,若是,跳至步骤S5,否则返回步骤S2。
优选的,误差门限εth的计算为:
其中,ferr表示根据实际场景设置的最大检测频率误差,fs表示采样频率;在实际应用中εth需要计算位宽、定点运算及应用场景带来的计算误差进行修正。
优选的,判断正弦波属于突发正弦波还是持续正弦波时,需要判断是否 Counter_Len大于等于Len_Testmax,若大于等于,则认为该段正弦波是一段持续正弦波序列,计算出频率,频率的计算公式为:
若不大于等于,则认为该段正弦波是一段突发正弦波序列,其频率同样可以由上述频率计算公式得出,起点位置和持续长度由参数Counter_Len确定。
本发明实现复杂度低,每一个点的运算量远低于FFT的运算量,比如对于64点的FFT,需要128次复数乘法,256次复数加法,逐点滑动时,每滑动一次都需要做一次FFT,而用本发明的方法,逐点滑动时,每滑动一次只需要做一次需要1次复数乘法,1次除法,1次反三角函数计算,前面采样点计算过的关于θn的结果可以重复使用。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序可以存储于一计算机可读存储介质中,存储介质可以包括:ROM、RAM、磁盘或光盘等。
此外,术语“第一”、“第二”、“第三”、“第四”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量,由此,限定有“第一”、“第二”、“第三”、“第四”的特征可以明示或者隐含地包括至少一个该特征,不能理解为对本发明的限制。
所应理解的是,本发明一种数字域正弦波检测装置与数字域正弦波检测方法是同一构思的不同实现方式,系统实施例未描述部分可以参考方法实施例相应部分的描述,反之亦然。
以上所举实施例,对本发明的目的、技术方案和优点进行了进一步的详细说明,所应理解的是,以上所举实施例仅为本发明的优选实施方式而已,并不用以限制本发明,凡在本发明的精神和原则之内对本发明所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种数字域正弦波检测方法,其特征在于,包括:
S1、设置检测正弦波畸变的参数,即令突发有效标志Flag_BurstValid为0、计数值Counter_Len为0,并设定正弦波的最小有效长度为Len_Validmin、最大测试序列长度为Len_Testmax;
S2、在正弦序列判决统计模块的A端输入第一输入序列An,B端输入为第二输入序列Bn
S3、利用输入的第一输入序列An和第二输入序列Bn,逐点计算相角θn
S4、判断相邻两点的相角θn的差是否满足|θn-1n|<εth,若不满足,令计数值Counter_Len=0,判断Flag_BurstValid==1是否成立,若成立则转到步骤S5,若不成立则令Counter_Len=0并返回步骤S2;若满足,则令Counter_Len自增1,然后判断Counter_Len>Len_Validmin是否成立,若成立则令Flag_BurstValid=1,并返回步骤S2;若不成立,则是否Counter_Len≥Len_Testmax,若是,跳至步骤S5,否则返回步骤S2;
S5、判断正弦波属于突发正弦波还是持续正弦波,根据相角θn和计数值Counter_Len的值,计算出所述正弦波的频率,得出正弦波的类别,返回步骤S1;
其中,εth表示误差门限。
2.根据权利要求1所述的数字域正弦波检测方法,其特征在于,步骤S2第一输入序列An和第二输入序列Bn的计算为:
对于实正弦波序列,A端输入的第一输入序列为An=X(n),B端输入的第二输入序列为Bn=X(n-Dm)+X(n+Dm);
对于复正弦波序列,A端输入的第一输入序列为An=Xr(n)·Xr(n-Dm)+Xc(n)·Xc(n-Dm),B端输入的第二输入序列为Bn=Xc(n)·Xr(n-Dm)-Xr(n)·Xc(n-Dm);
其中,X(n)表示实正弦波序列,Dm表示延迟单元的延迟参数,Xr(n)表示复正弦波序列的I路,Xc(n)表示复正弦波序列的Q路。
3.根据权利要求1所述的一种数字域正弦波检测方法,其特征在于,步骤S3相角θn的计算为:
对于实正弦波序列,相角为θn=arccosBn/2An
对于复正弦波序列,相角为θn=arctanBn/An
4.根据权利要求1所述的数字域正弦波检测方法,其特征在于,误差门限εth的计算为:
其中,ferr表示根据实际场景设置的最大检测频率误差,fs表示采样频率,Dm表示延迟单元的延迟参数。
5.根据权利要求1所述的数字域正弦波检测方法,其特征在于,所述根据相角θn和计数值Counter_Len的值,计算出正弦波的频率得出正弦波的类别,包括:
若Counter_Len≥Len_Testmax,则为一段持续正弦波序列,频率为f;否则为一段突发正弦波序列,频率为f,起点位置和持续长度由计数值Counter_Len确定,其中正弦波的频率表示为:
其中,fs为采样频率,Dm表示表示延迟单元的延迟参数。
6.一种数字域正弦波检测装置,其特征在于,包括:实正弦波序列输入模块、复正弦波序列输入模块和正弦序列判决统计模块,所述实正弦波序列输入模块包括带通滤波器、延迟单元和加法器,所述复正弦输入模块包括带通滤波器、延迟单元、乘法器和加法器,所述正弦序列判决统计模块的输入端包括A端和B端,分别在A端输入第一输入序列An、B端输入第二输入序列Bn
实正弦波序列输入模块用于计算实正弦波序列A端的第一输入序列An和B端的第二输入序列Bn
复正弦波序列输入模块用于计算复正弦波序列A端的第一输入序列An和B端的第二输入序列Bn
正弦序列判决统计模块用于判断相邻两点的相角θn的绝对差值与误差门限εth的关系,并计算出正弦波的频率。
7.根据权利要求6所述的数字域正弦波检测装置,其特征在于,第一输入序列An和第二输入序列Bn为:
对于实正弦波序列输入模块,A端输入为An=X(n),B端输入为Bn=X(n-Dm)+X(n+Dm);
对于复正弦波序列输入模块,A端输入为An=Xr(n)·Xr(n-Dm)+Xc(n)·Xc(n-Dm),B端输入为Bn=Xc(n)·Xr(n-Dm)-Xr(n)·Xc(n-Dm);
其中,X(n)表示实正弦波序列,Dm表示延迟单元的延迟参数,Xr(n)表示复正弦波序列的I路,Xc表示复正弦波序列的Q路。
8.根据权利要求6所述的数字域正弦波检测装置,其特征在于,相角θn的计算为:
对于实正弦波序列,相角为θn=arccosBn/2An
对于复正弦波序列,相角为θn=arctanBn/An
9.根据权利要求6所述的数字域正弦波检测装置,其特征在于,误差门限εth的计算为:
其中,ferr表示根据实际场景设置的最大检测频率误差,fs表示采样频率,Dm表示延迟单元的延迟参数。
10.根据权利要求6所述的数字域正弦波检测装置,其特征在于,正弦波的频率的计算为:
其中,fs为采样频率,Dm表示表示延迟单元的延迟参数,θn表示相角,Counter_Len表示计数值。
CN201810006884.XA 2018-01-04 2018-01-04 数字域正弦波检测方法及检测装置 Active CN108199715B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810006884.XA CN108199715B (zh) 2018-01-04 2018-01-04 数字域正弦波检测方法及检测装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810006884.XA CN108199715B (zh) 2018-01-04 2018-01-04 数字域正弦波检测方法及检测装置

Publications (2)

Publication Number Publication Date
CN108199715A true CN108199715A (zh) 2018-06-22
CN108199715B CN108199715B (zh) 2021-10-01

Family

ID=62587844

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810006884.XA Active CN108199715B (zh) 2018-01-04 2018-01-04 数字域正弦波检测方法及检测装置

Country Status (1)

Country Link
CN (1) CN108199715B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006292434A (ja) * 2005-04-06 2006-10-26 Sony Corp 角度検出信号処理装置
CN101105510A (zh) * 2006-07-11 2008-01-16 联发科技股份有限公司 相位误差测量电路与其方法
CN102072987A (zh) * 2010-11-13 2011-05-25 天津大学 短区间正弦信号的相位估计法及其实验装置
CN103562735A (zh) * 2011-04-29 2014-02-05 美国亚德诺半导体公司 用于检测电力系统的基频的系统和方法
CN104407198A (zh) * 2014-11-28 2015-03-11 深圳供电局有限公司 一种用于dvr装置中检测sag信号的方法和系统
CN107436619A (zh) * 2017-08-08 2017-12-05 重庆邮电大学 一种高精度低代价数字正弦波发生装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006292434A (ja) * 2005-04-06 2006-10-26 Sony Corp 角度検出信号処理装置
CN101105510A (zh) * 2006-07-11 2008-01-16 联发科技股份有限公司 相位误差测量电路与其方法
CN102072987A (zh) * 2010-11-13 2011-05-25 天津大学 短区间正弦信号的相位估计法及其实验装置
CN103562735A (zh) * 2011-04-29 2014-02-05 美国亚德诺半导体公司 用于检测电力系统的基频的系统和方法
CN104407198A (zh) * 2014-11-28 2015-03-11 深圳供电局有限公司 一种用于dvr装置中检测sag信号的方法和系统
CN107436619A (zh) * 2017-08-08 2017-12-05 重庆邮电大学 一种高精度低代价数字正弦波发生装置

Also Published As

Publication number Publication date
CN108199715B (zh) 2021-10-01

Similar Documents

Publication Publication Date Title
Yang et al. A new method of leak location for the natural gas pipeline based on wavelet analysis
Pedas et al. Piecewise polynomial collocation for linear boundary value problems of fractional differential equations
CN114037219A (zh) 一种数据评估方法、装置及电子设备
CN115388344B (zh) 一种基于分布式协调检测的管道泄漏监测方法
CN108900622A (zh) 基于物联网的数据融合方法、装置与计算机可读存储介质
CN108053424A (zh) 目标跟踪方法、装置、电子设备及存储介质
CN114460527B (zh) 关联度延拓Hilbert移相电子式互感器校验仪溯源方法和系统
Isermann et al. Two-step process identification with correlation analysis and least-squares parameter estimation
CN116205889A (zh) 偏移检测方法、装置、电子设备及存储介质
CN108199715A (zh) 数字域正弦波检测方法及检测装置
CN114137015A (zh) 孔隙度校正方法及装置
Mayer Efficient hierarchical triplet merging for camera pose estimation
CN109388781B (zh) 测量数据的处理方法和装置
CN111948454B (zh) 一种抗衰减直流分量干扰的同步相量测量方法及系统
CN106970265B (zh) 一种采用多时间尺度不完全s变换估计谐波参数的方法
Čelar et al. Methodology for control delay estimation using new algorithm for critical points identification
Ostapkowicz et al. Leak detection in liquid transmission pipelines during transient state related to a change of operating point
CN108731714B (zh) 一种频率扫描数据的解码方法及装置
Owen et al. Recent ship hydrodynamics developments in the parallel two-fluid flow solver Alya
CN114861984A (zh) 用于预测高含co2凝析气藏油环体积的方法及处理器
US8990277B2 (en) Method for searching a lookup table
CN106227988A (zh) 一种相对误差计算优化方法
CN106777913A (zh) 一种近似熵和样本熵共同最优参数m、r确定的新方法
Wang et al. Line detection algorithm based on adaptive gradient threshold and weighted mean shift
CN115291528B (zh) 模型不确定性等级确定方法、装置、系统及存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20241015

Address after: Room 802-2, Shuangchuang Building, No. 661 Taolin Road, Shuanggang Community, Zhangmuqiao Street, Changde Economic and Technological Development Zone, Hunan Province 415000

Patentee after: Hunan Shuangchou Network Technology Co.,Ltd.

Country or region after: China

Address before: 400065 Chongwen Road, Nanshan Street, Nanan District, Chongqing

Patentee before: CHONGQING University OF POSTS AND TELECOMMUNICATIONS

Country or region before: China