CN106291101A - 一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法 - Google Patents
一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法 Download PDFInfo
- Publication number
- CN106291101A CN106291101A CN201610898371.5A CN201610898371A CN106291101A CN 106291101 A CN106291101 A CN 106291101A CN 201610898371 A CN201610898371 A CN 201610898371A CN 106291101 A CN106291101 A CN 106291101A
- Authority
- CN
- China
- Prior art keywords
- matrix
- resolution
- frequency signal
- super
- estimation
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法;它包括以下步骤予以实现:1)构造增强矩阵;2)计算矩阵;3)进行特征值分解;4)构造矩阵;5)计算多项式;6)计算频率估计值。
Description
技术领域
本发明涉及信号处理领域,尤其涉及一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法。
背景技术
谐波信号的参数估计问题在多个信号处理领域有着广泛的应用,主要是从被噪声污染的观测信号中估计谐波信号的频率。依据噪声污染的不同情况一般可分为加性噪声和乘性和加性噪声两种情况。
目前,乘性和加性噪声中谐波信号频率的估计方法主要有循环统计量方法(李宏伟,程乾生.“乘性和加性噪声中谐波恢复的循环统计量方法”,电子学报,第26卷,第7期,1998年)和广义协方差矩阵方法(杨世永.“基于广义协方差矩阵的乘性和加性噪声中的谐波恢复”,信号处理,第28卷,第2期,2012年)。循环统计量方法是基于循环统计量,利用快速傅里叶变换和峰值搜索方法实现,由于受瑞利限的影响,循环统计量方法的估计精度和频率分辨率不高。广义协方差矩阵方法利用子空间旋转不变技术估计谐波信号的频率,但其实现过程复杂,频率估计的分辨率低。因此,有必要对上述技术问题予以解决。
发明内容
本发明的目的在于克服现有技术的不足,适应现实需要,提供一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法。
为了实现本发明的目的,本发明所采用的技术方案为:
设计一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法,它包括以下步骤:
1)构造增强矩阵D;
2)计算矩阵H;
3)对矩阵H进行特征值分解,将所得的特征值按从大到小的顺序记做λ1,λ2,…,λK,相
应的特征值向量记做e1,e2,…,eK;
4)构造矩阵F;
5)计算多项式Q(z);
6)计算频率估计值。
进一步的,所述构造增强矩阵D方法为:设谐波信号的N个数据测量值为x(1),x(2),…,x(N),谐波分量个数为P,对于一个取值范围在[2P+2,N/2]内的整数K,构造一个K×(N-K)的增强矩阵D,即:
进一步的,所述计算矩阵H的方法为:
其中(·)H表示共轭转置。
进一步的,所述的构造矩阵F的方法为:F=[eP+2,eP+3,…,eK]。
进一步的,所述计算多项式的方法为:令q(z)=[1,z,…,zK-1]T,计算多项式Q(z)=qT(z-1)FFHq(z),其中(·)T表示转置。
进一步的,计算频率估计值的方法为:求出Q(z)=0的所有根,将位于单位圆内且最靠近单位圆的P+1个根记做r1,r2,…,rP,rP+1。计算k=1,2,…,P+1,其中∠表示取幅角运算。去掉中绝对值最小的值,余下的P个值则为谐波信号频率的估计值。
本发明的有益效果在于:
本发明的方法可以实现对乘性和加性噪声中谐波信号进行超分辨率的频率估计,且本方法谐波信号频率估计的精度高,解决了现有方法频率估计的分辨率低和实现过程复杂的技术问题。
附图说明
图1为本发明方法的流程示意图;
具体实施方式
下面结合附图和实施例对本发明进一步说明:
实施例1:一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法,参见图1,本方法它包括以下步骤予以实现:
1)构造增强矩阵D:设谐波信号的N个数据测量值为x(1),x(2),…,x(N),谐波分量个数为P,对于一个取值范围在[2P+2,N/2]内的整数K,构造一个K×(N-K)的增强矩阵D,即:
2)计算矩阵H:其中(·)H表示共轭转置。
3)对矩阵H进行特征值分解,将所得的特征值按从大到小的顺序记做λ1,λ2,…,λK,相应的特征值向量记做e1,e2,…,eK。
4)构造矩阵F:F=[eP+2,eP+3,…,eK]。
5)计算多项式Q(z):令q(z)=[1,z,…,zK-1]T,计算多项式Q(z)=qT(z-1)FFHq(z),其中(·)T表示转置。
6)计算频率估计值:求出Q(z)=0的所有根,将位于单位圆内且最靠近单位圆的P+1个根记做r1,r2,…,rP,rP+1。计算k=1,2,…,P+1,其中∠表示取幅角运算。去掉中绝对值最小的值,余下的P个值则为谐波信号频率的估计值。
虽然,本发明的实施例公布的是较佳的实施例,但并不局限于此,本领域的普通技术人员,极易根据上述实施例,领会本发明的精神,并做出不同的引申和变化,但只要不脱离本发明的精神,都在本发明的保护范围内。
Claims (6)
1.一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法,其特征在于:它包括以下步骤:
1)构造增强矩阵D;
2)计算矩阵H;
3)对矩阵H进行特征值分解,将所得的特征值按从大到小的顺序记做λ1,λ2,…,λK,相应的特征值向量记做e1,e2,…,eK;
4)构造矩阵F;
5)计算多项式Q(z);
6)计算频率估计值。
2.如权利要求1所述的一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法,其特征在于:所述构造增强矩阵D方法为:设谐波信号的N个数据测量值为x(1),x(2),…,x(N),谐波分量个数为P,对于一个取值范围在[2P+2,N/2]内的整数K,构造一个K×(N-K)的增强矩阵D,即:
3.如权利要求1所述的一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法,其特征在于:所述计算矩阵H的方法为:
其中(·)H表示共轭转置。
4.如权利要求1所述的一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法,其特征在于:所述的构造矩阵F的方法为:F=[eP+2,eP+3,…,eK]。
5.如权利要求1所述的一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法,其特征在于:所述计算多项式的方法为:令q(z)=[1,z,…,zK-1]T,计算多项式Q(z)=qT(z-1)FFHq(z),其中(·)T表示转置。
6.如权利要求1所述的一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法,其特征在于:计算频率估计值的方法为:求出Q(z)=0的所有根,将位于单位圆内且最靠近单位圆的P+1个根记做r1,r2,…,rP,rP+1。计算k=1,2,…,P+1,其中∠表示取幅角运算。去掉中绝对值最小的值,余下的P个值则为谐波信号频率的估计值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610898371.5A CN106291101B (zh) | 2016-10-14 | 2016-10-14 | 一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610898371.5A CN106291101B (zh) | 2016-10-14 | 2016-10-14 | 一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106291101A true CN106291101A (zh) | 2017-01-04 |
CN106291101B CN106291101B (zh) | 2018-12-18 |
Family
ID=57717336
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610898371.5A Expired - Fee Related CN106291101B (zh) | 2016-10-14 | 2016-10-14 | 一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106291101B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107800658A (zh) * | 2017-11-10 | 2018-03-13 | 九江学院 | 一种二维谐波信号频率估计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104330622A (zh) * | 2014-10-29 | 2015-02-04 | 广东电网有限责任公司电力科学研究院 | 电力系统中正弦波信号的频率测量方法及系统 |
CN105044453A (zh) * | 2015-08-11 | 2015-11-11 | 杨世永 | 一种适用于复杂噪声背景的谐波信号频率估计方法 |
CN105301354A (zh) * | 2015-07-12 | 2016-02-03 | 九江学院 | 一种乘性和加性噪声中谐波信号频率估计方法 |
-
2016
- 2016-10-14 CN CN201610898371.5A patent/CN106291101B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104330622A (zh) * | 2014-10-29 | 2015-02-04 | 广东电网有限责任公司电力科学研究院 | 电力系统中正弦波信号的频率测量方法及系统 |
CN105301354A (zh) * | 2015-07-12 | 2016-02-03 | 九江学院 | 一种乘性和加性噪声中谐波信号频率估计方法 |
CN105044453A (zh) * | 2015-08-11 | 2015-11-11 | 杨世永 | 一种适用于复杂噪声背景的谐波信号频率估计方法 |
Non-Patent Citations (1)
Title |
---|
窦慧晶 等: "乘性和加性噪声相关背景下的二维谐波频率估计", 《电子学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107800658A (zh) * | 2017-11-10 | 2018-03-13 | 九江学院 | 一种二维谐波信号频率估计方法 |
CN107800658B (zh) * | 2017-11-10 | 2020-05-08 | 九江学院 | 一种二维谐波信号频率估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106291101B (zh) | 2018-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101833035B (zh) | 线性调频信号参数估计方法及其实施装置 | |
CN105044453B (zh) | 一种适用于复杂噪声背景的谐波信号频率估计方法 | |
CN202339381U (zh) | 一种纳托尔自卷积窗加权傅里叶变换的谐波电能计量系统 | |
CN106772226A (zh) | 基于压缩感知时间调制阵列的doa估计方法 | |
Monovasilis et al. | Two new phase-fitted symplectic partitioned Runge–Kutta methods | |
CN102608565A (zh) | 一种基于均匀圆阵列的波达方向估计方法 | |
CN103971029A (zh) | 一种用于网格失配下的doa估计的交替迭代方法 | |
CN102662158B (zh) | 一种对传感器天线阵列接收信号的快速处理方法 | |
Park | Guaranteed-stable sliding DFT algorithm with minimal computational requirements | |
CN104535987A (zh) | 适用于均匀圆阵列声纳系统的幅相误差自校正方法 | |
CN103901396A (zh) | 相干信号源亚分辨率超分辨到达角估计方法 | |
CN105301354B (zh) | 一种乘性和加性噪声中谐波信号频率估计方法 | |
Koteswara Rao et al. | Accurate phasor and frequency estimation during power system oscillations using least squares | |
CN106291101A (zh) | 一种具有超分辨率的乘性和加性噪声中谐波信号频率估计方法 | |
Peng et al. | A blind calibration model for I/Q imbalances of wideband zero-IF receivers | |
CN104914451B (zh) | 一种块Toeplitz矩阵低复杂度求逆的空时抗干扰方法 | |
CN104731762B (zh) | 基于循环移位的立方相位信号参数估计方法 | |
CN106507952B (zh) | 一种基于圆阵的快速空间谱计算方法 | |
CN106569182A (zh) | 基于最小熵的相位编码信号载频估计方法 | |
CN103257270A (zh) | 基于傅里叶级数的锯齿载波混沌spwm频谱分析方法 | |
CN105675084B (zh) | 一种带有迭代插值的高精度液位测量方法 | |
CN103365827A (zh) | 一种高精度正弦/余弦函数的计算方法 | |
Xu et al. | Inter‐harmonic parameters estimation in power grid based on accelerated PSO and T5R11 window | |
Huang et al. | Nonlinear Blind Compensation for Array Signal Processing Application | |
D’Ambrosio et al. | Collocation–Based Two Step Runge–Kutta Methods for Ordinary Differential Equations |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181218 Termination date: 20201014 |