CN109374969A - 一种基于相干积累的dft的下采样正弦信号相位差测量方法 - Google Patents
一种基于相干积累的dft的下采样正弦信号相位差测量方法 Download PDFInfo
- Publication number
- CN109374969A CN109374969A CN201811160519.0A CN201811160519A CN109374969A CN 109374969 A CN109374969 A CN 109374969A CN 201811160519 A CN201811160519 A CN 201811160519A CN 109374969 A CN109374969 A CN 109374969A
- Authority
- CN
- China
- Prior art keywords
- signal
- groups
- phase difference
- sampling
- sinusoidal signal
- 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
- G01R25/00—Arrangements for measuring phase angle between a voltage and a current or between voltages or currents
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公开了信号处理领域的一种基于相干积累的DFT的下采样正弦信号相位差测量方法,包括以下步骤:使用下采样的方式分别对两组频率相同的正弦信号进行采样得到两组采样信号;分别对两组采样信号进行相干积累得到两组积累信号;分别对两组累积信号进行离散傅里叶变换,得到两组频率相同的正弦信号的初始相位;根据两组频率相同的正弦信号之间的初始相位计算两组正弦信号之间的相位差。
Description
技术领域
本发明涉及信号处理领域,具体公开了一种基于相干积累的DFT的下采样正弦信号相位差测量方法。
背景技术
正弦信号的相位差测量是电力***监测、射频通信、雷达定位、激光测距等许多应用中的一个重要研究课题。此外,它还可以应用于相位校准,特别是在单航过星载/机载InSAR***上的相位校准中尤其有用。
长久以来,人们已经讨论了许多不同的相位差测量方法,如最小二乘或离散傅里叶变换(DFT)法、希尔伯特变换法、ICA法和零交叉检测法等。然而,基于DFT的频谱分析方法在计算时需要一个完整的采样周期;而零交叉检测法在有噪声叠加于信号中的情况下则不可靠;此外,当信号频率太高时,缺乏采样频率如何选择以及采用何种相位差测量方法的研究。
发明内容
本发明目的在提供,以解决现有技术中信号频率太高时无法获取采用正弦信号相位差的技术缺陷。
为实现上述目的,本发明提供了一种基于相干积累的DFT的下采样正弦信号相位差测量方法,包括以下步骤:
S1:使用下采样的方式分别对两组频率相同的正弦信号进行采样得到两组采样信号;
S2:分别对两组采样信号进行相干积累得到两组积累信号;
S3:分别对两组累积信号进行离散傅里叶变换,得到两组频率相同的正弦信号的初始相位;
S4:根据两组频率相同的正弦信号之间的初始相位计算两组正弦信号之间的相位差。
优选地,S1包括以下步骤:
S11:对所述正弦信号进行采样获得采样后的信号频谱;
S12:对信号频谱进行滤波得到基带信号频谱;
S13:根据所述基带信号频谱采用逆傅里叶变换重构基带信号;
S14:根据基带信号的频率和初相位重建高频信号。
优选地,S11获得的信号频谱如下:
其中,ω表示角频率,fs为采样频率。
优选地,S12对信号进行滤波的增益为1/fs,带通范围为0~0.5fs。
优选地,S13的基带信号如下:
其中,f0为基带信号的频率,为基带信号的初相位。
优选地,S14的高频信号如下:
其中,T为采样周期。
优选地,S3的初始相位如下:
本发明具有以下有益效果:
1、本发明的测量方法能有效地滤除信号中的噪声,提高正弦信号相位差的测量精度,大大减少测量误差;
2、本发明的测量方法在信号频率太高时,下采样条件下采样频率的选择更加广泛。
下面将参照附图,对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是本发明优选实施例的一种基于相干积累的DFT的下采样正弦信号相位差测量方法流程图。
图2是本发明优选实施例的信号频谱图;
图3是本发明优选实施例的相干累积前后的信号图;
图4是本发明优选实施例的相位差示意图。
具体实施方式
以下结合附图对本发明的实施例进行详细说明,但是本发明可以由权利要求限定和覆盖的多种不同方式实施。
DFT全名Discrete Fourier Transform,即为离散傅里叶变换,是傅里叶变换在时域和频域上都呈离散的形式。
本发明首先提供了一种基于相干积累的DFT的下采样正弦信号相位差测量方法,参见图1,包括以下步骤:
S1:使用下采样的方式分别对两组频率相同的正弦信号进行采样得到两组采样信号。
设正弦信号s(t)的采样频率为fs,其数学表达式为:
其中A是未知的振幅,f是信号频率,t是时间,是未知的初相位。根据Nyquist采样定理,fs必须大于等于2f才可以准确地恢复原始信号。但是,当需要测量两个正弦信号的相位差时,fs必须远大于2f。然而,当信号频率本身就非常高时,随着信号频率的增加,采样频率也会越来越高,甚至在现有设备与技术条件下很难达到,以致难以实现。因此,需要根据带通采样定理来降低采样频率,也即使用下采样的方式来对信号进行采样。
使用下采样的方式对正弦信号进行采样又包括以下步骤:
S11:对所述正弦信号进行采样获得采样后的信号频谱。
信号s(t)的频谱S(ω)如图2(a)所示。其中ω表示角频率,箭头代表振幅谱,实心黑点代表相位谱。采样后信号的频谱表达式如下:
采样后信号的频谱Ss(ω)如图2(b)所示,显然,为避免频谱混叠,采样频率和信号频率之间须满足以下条件:
-f+nfs≠f,n=1,2,3,… (3)
也即:
S12:对信号频谱进行滤波得到基带信号频谱。
采样频谱经过增益为1/fs,带通范围为0~0.5fs的滤波器以后,剩下的即为基带信号的频谱。此时,可能存在两种情况,如图2(c)、(d)所示。其中,标有“1”的部分是信号频谱S(ω)右移n次的结果,标有“2”的部分是信号频谱S(ω)左移n次的结果。
在图2(c)所示情况下,需满足条件:
0<f-nfs<0.5fs,n=1,2,3,… (5)
也即:
此时得到的基带信号频谱为:
在图2(d)所示情况下,则需满足条件:
0<-f+nfs<0.5fs,n=1,2,3,… (10)
也即:
此时得到的基带信号频谱为:
S13:根据所述基带信号频谱采用逆傅里叶变换重构基带信号。
对式(7)进行逆傅里叶变化得到重构基带信号:
其中,f0为基带信号的频率,为基带信号的初相位,由式(8)可知
即基带信号y(t)的初相位与信号s(t)的初相位相同。
对式(11)进行逆傅里叶变化得到重构基带信号:
由式(13)可知
即基带信号y(t)的初相位与信号s(t)的初相位相反。
S14:根据基带信号的频率和初相位重建高频信号。
根据低频基带信号y(t)的频率f0和初始相位可以重建高频信号s(t);并且可以通过选择满足条件(6)或(11)的采样频率来测量两个频率相同的正弦信号的相位差。
假设采样频率满足式(11)的条件,采样长度为N点,而在一个基带信号周期中采样点数为N0,且满足关系N=m·N0(m为正整数),则采样后的信号即为:
其中,T为采样周期(T=1/fs)。
S2:分别对两组采样信号进行相干积累得到两组积累信号。
将采样后的信号s1(kT)和s2(kT)以基带信号y(t)的周期T0(T0=N0/fs)为周期进行相干积累,以滤除高斯噪声。即
采样结果如图3(a)所示。
采样结果如图3(b)所示。
S3:分别对两组累积信号进行离散傅里叶变换,得到两组频率相同的正弦信号的初始相位。
对积累后的信号s1a(kT)和s2a(kT)分别进行离散傅里叶变换,以计算出两正弦信号的初始相位和
采样结果如图3(c)所示。
采样结果如图3(d)所示。
S4:根据两组频率相同的正弦信号之间的初始相位计算两组正弦信号之间的相位差。
根据两正弦信号的初相位即可得到它们之间的相位差:
所得相位差结果如图4(c)所示,相位差误差如图4(d)。而采用传统的DFT方法测量得到的相位差如图4(a)所示,相位差误差如图4(b)所示。比较可见,本方法相较于传统的DFT方法测量的相位差更准确,相位差误差更小。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种基于相干积累的DFT的下采样正弦信号相位差测量方法,其特征在于,包括以下步骤:
S1:使用下采样的方式分别对两组频率相同的正弦信号进行采样得到两组采样信号;
S2:分别对两组采样信号进行相干积累得到两组积累信号;
S3:分别对两组累积信号进行离散傅里叶变换,得到两组频率相同的正弦信号的初始相位;
S4:根据两组频率相同的正弦信号之间的初始相位计算两组正弦信号之间的相位差。
2.根据权利要求1所述的一种基于相干积累的DFT的下采样正弦信号相位差测量方法,其特征在于,所述S1包括以下步骤:
S11:对所述正弦信号进行采样获得采样后的信号频谱;
S12:对信号频谱进行滤波得到基带信号频谱;
S13:根据所述基带信号频谱采用逆傅里叶变换重构基带信号;
S14:根据基带信号的频率和初相位重建高频信号。
3.根据权利要求2所述的一种基于相干积累的DFT的下采样正弦信号相位差测量方法,其特征在于,所述S11获得的信号频谱如下:
其中,ω表示角频率,fs为采样频率。
4.根据权利要求2所述的一种基于相干积累的DFT的下采样正弦信号相位差测量方法,其特征在于,所述S12对信号进行滤波的增益为1/fs,带通范围为0~0.5fs。
5.根据权利要求2所述的一种基于相干积累的DFT的下采样正弦信号相位差测量方法,其特征在于,所述S13的基带信号如下:
其中,f0为基带信号的频率,为基带信号的初相位。
6.根据权利要求2所述的一种基于相干积累的DFT的下采样正弦信号相位差测量方法,其特征在于,所述S14的高频信号如下:
其中,T为采样周期。
7.根据权利要求1所述的一种基于相干积累的DFT的下采样正弦信号相位差测量方法,其特征在于,所述S3的初始相位如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811160519.0A CN109374969B (zh) | 2018-09-30 | 2018-09-30 | 一种基于相干积累的dft的下采样正弦信号相位差测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811160519.0A CN109374969B (zh) | 2018-09-30 | 2018-09-30 | 一种基于相干积累的dft的下采样正弦信号相位差测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109374969A true CN109374969A (zh) | 2019-02-22 |
CN109374969B CN109374969B (zh) | 2021-01-15 |
Family
ID=65402831
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811160519.0A Active CN109374969B (zh) | 2018-09-30 | 2018-09-30 | 一种基于相干积累的dft的下采样正弦信号相位差测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109374969B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114019235A (zh) * | 2021-09-22 | 2022-02-08 | 浙江大学 | 一种基于同步离散傅里叶变换的频率响应测量***及方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0378405A2 (en) * | 1989-01-13 | 1990-07-18 | Hewlett-Packard Company | Method and apparatus for measuring phase accuracy and amplitude profile of a continuous-phase-modulation signal |
US6868129B2 (en) * | 2001-03-12 | 2005-03-15 | Freescale Semiconductor, Inc. | Demodulator for a radio receiver and method of operation |
US20100156690A1 (en) * | 2008-12-22 | 2010-06-24 | Electronics And Telecommunications Research Institute | Digital direct conversion receiving apparatus and method |
CN101825660A (zh) * | 2010-05-05 | 2010-09-08 | 天津大学 | 欠采样下的正弦信号频率的高效测量方法及实施装置 |
CN102200550A (zh) * | 2011-05-31 | 2011-09-28 | 中国航空无线电电子研究所 | 一种用于高精度测量相位差的延迟正交数字中频鉴相方法 |
CN102879784A (zh) * | 2012-10-25 | 2013-01-16 | 西安电子科技大学 | 四种模式的合成孔径雷达统一成像方法 |
CN105652086A (zh) * | 2015-12-29 | 2016-06-08 | 国网重庆市电力公司电力科学研究院 | 一种初始相位角的确定方法和*** |
CN106918741A (zh) * | 2017-03-02 | 2017-07-04 | 浙江大学 | 应用于频率宽范围波动电网的自适应采样的相位差校正法 |
-
2018
- 2018-09-30 CN CN201811160519.0A patent/CN109374969B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0378405A2 (en) * | 1989-01-13 | 1990-07-18 | Hewlett-Packard Company | Method and apparatus for measuring phase accuracy and amplitude profile of a continuous-phase-modulation signal |
US6868129B2 (en) * | 2001-03-12 | 2005-03-15 | Freescale Semiconductor, Inc. | Demodulator for a radio receiver and method of operation |
US20100156690A1 (en) * | 2008-12-22 | 2010-06-24 | Electronics And Telecommunications Research Institute | Digital direct conversion receiving apparatus and method |
CN101825660A (zh) * | 2010-05-05 | 2010-09-08 | 天津大学 | 欠采样下的正弦信号频率的高效测量方法及实施装置 |
CN102200550A (zh) * | 2011-05-31 | 2011-09-28 | 中国航空无线电电子研究所 | 一种用于高精度测量相位差的延迟正交数字中频鉴相方法 |
CN102879784A (zh) * | 2012-10-25 | 2013-01-16 | 西安电子科技大学 | 四种模式的合成孔径雷达统一成像方法 |
CN105652086A (zh) * | 2015-12-29 | 2016-06-08 | 国网重庆市电力公司电力科学研究院 | 一种初始相位角的确定方法和*** |
CN106918741A (zh) * | 2017-03-02 | 2017-07-04 | 浙江大学 | 应用于频率宽范围波动电网的自适应采样的相位差校正法 |
Non-Patent Citations (1)
Title |
---|
殷丕磊 等: "分布式全相参雷达的相位差跟踪技术", 《信号处理》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114019235A (zh) * | 2021-09-22 | 2022-02-08 | 浙江大学 | 一种基于同步离散傅里叶变换的频率响应测量***及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109374969B (zh) | 2021-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ramos et al. | Simulation and experimental results of multiharmonic least-squares fitting algorithms applied to periodic signals | |
CN103454495B (zh) | 自适应高精度快速频谱分析方法 | |
Djurić et al. | Frequency measurement of distorted signals using Fourier and zero crossing techniques | |
CN103823177B (zh) | 基于窗函数设计的滤波器性能检测方法和*** | |
CN109374969A (zh) | 一种基于相干积累的dft的下采样正弦信号相位差测量方法 | |
Hua et al. | Matching and reassignment based time-frequency enhancement for rotating machinery fault diagnosis under nonstationary speed operations | |
Yi et al. | Time-varying fault feature extraction of rolling bearing via time–frequency sparsity | |
CN109239463B (zh) | 一种基于线性修正算法的介质损耗测量方法 | |
Liu et al. | Synchronous fault feature extraction for rolling bearings in a generalized demodulation framework | |
CN102095940A (zh) | 一种介质损耗角的测量方法和装置 | |
Du et al. | Fractional iterative variational mode decomposition and its application in fault diagnosis of rotating machinery | |
CN114050867B (zh) | 信号补偿方法、装置、设备及存储介质 | |
Serov et al. | Reseach of influency of mains frequency deviation on voltage spectrum measurement error by DFT method | |
CN113406386B (zh) | 一种基于数字下变频的信号频率精确估计方法 | |
Kuhlmann et al. | Effects of sampling rate and ADC width on the accuracy of amplitude and phase measurements in power-quality monitoring | |
Su et al. | A frequency estimation algorithm based on cross information fusion | |
Belega et al. | Effect of frequency uncertainty on the sine-wave amplitude estimator returned by linear least-squares fitting | |
Nishioka et al. | Tests of AMiBA Data Integrity | |
Belega et al. | Assessment of influence of systematic errors on the precision with which the normalized frequency of a sinusoidal signal is determined by means of a discrete Fourier transformation with interpolation | |
Haoning et al. | Research on envelope demodulation of reducer shock signals based on peak retention and down-sampling method | |
Gao et al. | A multiscale products technique for denoising of DNA capillary electrophoresis signals | |
Alrubei et al. | A new accurate estimator of the frequency using the three-point interpolation of DFT samples | |
Dubey et al. | Analysis of Chirp Signal with Fractional Order Savitzky-Golay Digital Differentiator | |
Machado et al. | Phase sensitive detection for embedded sensors | |
Xu et al. | A new method for instantaneous frequency rate estimation of polynomial phase signal |
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 |