CN108318737A - 一种高精度数字鉴相方法 - Google Patents
一种高精度数字鉴相方法 Download PDFInfo
- Publication number
- CN108318737A CN108318737A CN201810060448.0A CN201810060448A CN108318737A CN 108318737 A CN108318737 A CN 108318737A CN 201810060448 A CN201810060448 A CN 201810060448A CN 108318737 A CN108318737 A CN 108318737A
- Authority
- CN
- China
- Prior art keywords
- signal
- zero
- signals
- straight line
- ways
- 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.)
- Pending
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 Phase Differences (AREA)
Abstract
一种高精度数字鉴相方法,用于确定两路模拟信号的相位差,其鉴相过程为使用双通道同步采样单元对两路信号进行同步采样并经过模数转换得到两路数字信号,对两路数字信号分别求平均值,初步确定其交流信号零点位置对应幅值,在两路信号零点幅值位置附近取一定长度的线性段为过零段,并采用直线拟合算法由信号的过零段得到信号的过零点,采用过零点的位置为参照,得到信号波峰位置,由两路数字信号各自的相邻波峰位置分别确定信号周期,并结合由两路数字信号波峰位置得到的时间差,计算得到两路信号的相位差。本发明消除了零位误差导致的过零点位置误差,因此具有极高的鉴相精度,同时该方法也适用于正弦信号之外其他类型周期信号的鉴相测量。
Description
技术领域
本发明涉及信号处理与分析技术领域,尤其涉及一种数字鉴相方法。
背景技术
鉴相器用于测量两路信号的相位差,广泛应用于激光相位测距、干涉测向、电机控制、锁相环、激光位移干涉仪等领域。鉴相器分为模拟鉴相器和数字鉴相器,其中模拟鉴相器采用分立元件,线性鉴相范围窄,测量精度低。而数字鉴相器又分为数字电路鉴相器和数字算法鉴相器,其中数字算法鉴相器通过采样得到离散的数字信号,采用数字信号处理算法测量相位差。数字算法鉴相器使用CPU、MPU或者DSP等处理器做为核心处理单元,电路简单、抗干扰能力强、测量精度高,在高精度鉴相应用中具有明显的优势。
典型的数字鉴相算法有数字相关鉴相法、FFT谱分析鉴相法、过零检测鉴相法等。数字相关鉴相法的原理为两个同频率的正弦信号的零延时互相关函数值与其相位差的余弦值成正比。由于噪声与有效信号的相关性很小,所以该方法有很好的噪声抑制能力。数字相关鉴相法的缺陷为根据反余弦函数计算相位差,在相位差接近0时,灵敏度很低;由于余弦函数为偶函数,该算法无法分辨相位差的正负号;同时该方法只能测量正弦信号相位差,对其它周期信号或者非周期信号无法测量。FFT谱分析鉴相法对信号进行傅里叶级数分析,得到信号的基波参数。步骤为:通过FFT得到信号的幅度谱和相位谱,幅度谱的峰值频率做为信号的基波频率,再由相位谱中查找基波频率对应的相位,做为基波相位。谱分析鉴相法的缺点是只能对正弦波形的周期信号进行初相位测量,无法对非周期信号进行测量,同时FFT采用有限序列计算信号频谱,不可避免会产生频谱泄露,从而带来相位测量误差。要减小相位差测量误差,就必须提高频率分辨率。当采样频率一定时,只能通过增加数据长度来提高FFT的频率分辨率,计算量大。过零检测鉴相法是通过测量两信号过零点的时间差计算相位差,过零检测鉴相法原理简单可以采用电路来实现,也可采用算法来实现。信号在过零点附近具有较大的斜率,所以过零点检测具有较高的灵敏度。但是,受AD偏置误差、运放输入失调等因素的影响,信号的零位通常会出现误差,零位误差引起过零检测的误差。另外,受AD采样间隔、AD量化误差以及信号噪声的影响,直接由AD采样得到的过零点时,过零点的时间分辨率和精度都较差。
发明内容
本发明提供一种高精度数字鉴相方法,以消除传统过零检测鉴相方法中零位误差和AD量化误差引起的鉴相误差。
为了实现上述目的,本发明提供的技术方案是:使用双通道同步采样单元对两路模拟信号同步采样,并通过数模转换得到两路数字信号,通过对信号整体求平均值初步确定交流信号零点位置,并在初步确定的零点位置前后取一段线性段的信号数据作为过零段,采用直线拟合算法由该过零段得到信号的过零点,并以所得过零点的位置做为参照,以两个相邻过零点的中心位置为信号的波峰位置,由信号两相邻波峰位置的间隔为信号周期,并以两路信号波峰位置的时间差,计算得到两路信号的相位差。
上述方法具体包含以下步骤:
1、采用双通道同步采样单元对两路模拟信号进行同步采样,并经过模数转换得到两路数字信号;
2、对两路数字信号分别求平均值,初步确定交流信号零点位置对应幅值,在两路信号零点幅值位置前后取一定长度的线性段为过零段;
3、对过零段进行直线拟合,并计算拟合直线与时间轴的交点,得到信号的过零点;
4、由过零段直线拟合得到的两个相邻过零点的中点做为信号波峰位置,测量两个相邻波峰位置的时间差作为信号周期T;
5、测量两路数字信号波峰位置的时间差,得到两路信号的延时ΔT,再根据信号周期T,由公式:
计算得到两路信号的相位差。
进一步的,步骤3中的直线拟合方法是:采用方法为一元线性回归法,回归方程为:
y=ax+b,
根据一元线性回归理论,可以得道拟合直线的系数:
其中,
式中,n为线性段采样点数量。
进一步的,步骤3中的直线拟合方法是建立超定线性方程组并求解得到拟合直线的最小二乘解,从而确定拟合直线的两个系数,拟合直线用y=ax+b表示,对信号过零段中的任意一点都可以列出一个线性方程:
yi=axi+b,
其中yi为信号电压,xi为采样时刻。
由此得到一个超定线性方程组,用矩阵形式表示为AX=B,对于由n个点组成的过零段,其中,系数矩阵A和常数项向量B分别为:
X为拟合系数a和b构成的列矩阵,即
通过变换ATAX=ATB得到正定线性方程组,求解该正定方程组,可以得到拟合直线的最小二乘解。
与现有技术相比,本发明的有益效果如下:
1.本算法以信号的峰值为信号特征点,测量两路信号的周期与延时。峰值是模拟信号最基本的特征,对非正弦波形式的周期信号也适用,和其他鉴相算法相比,具有更强的适用性。
2.采用过零点来定位信号波峰点,信号在过零点的斜率最大,所以算法具有较高的灵敏度。
3.过零点定位采用信号过零附近的一段信号进行直线拟合,得到过零直线的最小二乘解,再由过零直线得到信号过零点。最小二乘解对噪声的容忍能力很强,同时该算法的时间分辨率不受AD的采样率限制,因此在低信噪比和低采样率的条件下,过零点定位和峰值点定位的精度也较高,保证了鉴相的精度。
4.通过峰值点定位,可以补偿***的偏置误差。不需要对电路的偏置以及AD的偏置误差进行标定。
5.该算法的鉴相范围是[-π,π],并在整个鉴相区间具有相同的灵敏度。对比相关鉴相法,相关鉴相法只有在相位差为±π/2附近时是线性的,灵敏度高;当相位差接近0和π时灵敏度很低。
6.本发明方法对于正弦信号以外的周期信号也适用。
附图说明
图1是以本发明所述鉴相方法构建的数字算法鉴相器的原理框图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
本发明提供的一种高精度数字鉴相方法,使用双通道同步采样单元对两路模拟信号同步采样,并通过数模转换得到两路数字信号,通过对信号整体求平均值初步确定交流信号零点位置,并在初步确定的零点位置前后取一段线性段的信号数据作为过零段,采用直线拟合算法由该过零段得到信号的过零点,并以所得过零点的位置做为参照,以两个相邻过零点的中心位置为信号的波峰位置,由信号两相邻波峰位置的间隔为信号周期,并以两路信号波峰位置的时间差,计算两路信号的相位差。
参见图1,本发明所述高精度数字鉴相方法构建的数字算法鉴相器,由双通道同步采样单元1、过零检测单元2、波峰定位单元3、周期计算单元4、延时计算单元5以及鉴相输出单元6构成,双通道同步采样单元1、过零检测单元2和波峰定位单元3依次相接,波峰定位单元3产生的信号波峰位置同时输出到周期计算单元4和延时计算单元5,周期计算单元4和延时计算单元5分别输出信号周期T和两路信号之间的延时ΔT至鉴相输出单元6。所说的双通道同步采样单元使用独立的AD芯片实现,也可以使用ARM、MCU、DSP等芯片内置的AD模块实现。
其基本工作过程是,使用双通道同步采样单元1对两路模拟信号进行同步采样,并通过模数转换得到两路数字信号。
过零检测单元2,采用信号整体求平均的方法,初步确定零点位置,并在初步确定的零点位置前后一定范围取线性段,对线性段做线性拟合,精确确定零点位置。
波峰定位单元3,由精确确定的相邻零点位置取中间值,作为信号波峰位置。信号波峰位置输出到周期计算单元4和延时计算单元5。
周期计算单元4,以两相邻波峰位置的时间间隔为信号周期T。
延时计算单元5,由两路信号的波峰位置对应时间差确定两路信号之间的延时ΔT。
鉴相输出单元6,根据周期计算单元4和延时计算单元5分别输出的信号周期T和两路信号之间的延时ΔT,由公式:
计算得到两路信号的相位差,并输出。
一种高精度数字鉴相方法具体包含以下步骤:
1、采用双通道同步采样单元对两路模拟信号进行同步采样,并经过模数转换得到两路数字信号;
2、对两路数字信号分别求平均值,初步确定交流信号零点位置对应幅值,在两路信号零点幅值位置前后取一定长度的线性段为过零段;
3、对过零段进行直线拟合,并计算拟合直线与时间轴的交点,得到信号的过零点;
直线拟合方法之一是:采用方法为一元线性回归法,回归方程为:
y=ax+b,
根据一元线性回归理论,可以得道拟合直线的系数:
其中,
式中,n为线性段采样点数量。
直线拟合方法之二是建立超定线性方程组并求解得到拟合直线的最小二乘解,从而确定拟合直线的两个系数,拟合直线用y=ax+b表示,对信号过零段中的任意一点都可以列出一个线性方程:
yi=axi+b,
其中yi为信号电压,xi为采样时刻。
由此得到一个超定线性方程组,用矩阵形式表示为AX=B,对于由n个点组成的过零段,系数矩阵A和常数项向量B分别为:
X为拟合系数a和b构成的列矩阵,表示为:
通过变换ATAX=ATB得到正定线性方程组,求解该正定方程组,可以得到拟合直线的最小二乘解。
4、由过零段直线拟合得到的两个相邻过零点的中点做为信号波峰位置,测量两个相邻波峰位置的时间差作为信号周期T;
5、测量两路数字信号波峰位置的时间差,得到两路信号的延时ΔT,再根据信号周期T,由公式:
计算得到两路信号的相位差。
以上应用了具体个例对本发明进行阐述,只是用于帮助理解本发明,并不用以限制本发明。对于本发明所属技术领域的技术人员,依据本发明的思想,还可以做出若干简单推演、变形或替换。
Claims (4)
1.一种高精度数字鉴相方法,使用双通道同步采样单元对两路模拟信号进行同步采样,并经过模数转换得到两路数字信号,对两路数字信号分别求平均值,初步确定信号零点位置对应幅值,在两路信号零点幅值位置附近取一定长度的线性段为过零段,并采用直线拟合算法由信号的过零段得到信号的过零点,采用过零点的位置为参照,得到信号波峰位置,由两路数字信号各自的相邻波峰位置分别确定信号周期,并结合由两路数字信号波峰位置得到的时间差,计算得到两路信号的相位差。
2.如权利要求1所述的高精度数字鉴相方法,其特征在于,包含以下步骤:
步骤1:对两路模拟信号进行同步采样并经过模数转换得到两路数字信号;
步骤2:对两路数字信号分别求平均值,初步确定信号零点位置对应幅值,在两路信号零点幅值位置前后取一定长度的线性段为过零段;
步骤3:对过零段进行直线拟合,并计算拟合直线与时间轴的交点,得到信号的过零点;
步骤4:由过零段直线拟合得到的两个相邻过零点的中点做为信号波峰位置,测量两个相邻波峰位置的时间差作为信号周期T;
步骤5:测量两路数字信号波峰位置的时间差,得到两路信号的延时ΔT,再根据信号周期T,由公式:
计算得到两路信号的相位差。
3.如权利要求2所述高精度数字鉴相方法,其特征在于:所述步骤3中过零段的直线拟合,采用一元线性回归法,回归方程为:
y=ax+b,
根据一元线性回归理论,可以得到拟合直线的系数:
其中,
式中,n为线性段采样点数量。
4.如权利要求2所述一种高精度数字鉴相方法,其特征在于:所述步骤3中过零段的直线拟合,采用建立超定线性方程组并求解得到拟合直线的最小二乘解,从而确定拟合直线的两个系数,拟合直线用y=ax+b表示,对信号过零段中的任意一点都可以列出一个线性方程:
yi=axi+b,
其中yi为信号电压,xi为采样时刻。
由此得到一个超定线性方程组,用矩阵形式表示为AX=B,对于由n个点组成的过零段,系数矩阵A和常数项向量B分别为:
X为拟合系数a和b构成的列矩阵,即:
通过变换ATAX=ATB得到正定线性方程组,求解该正定方程组,可以得到拟合直线的最小二乘解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810060448.0A CN108318737A (zh) | 2018-01-22 | 2018-01-22 | 一种高精度数字鉴相方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810060448.0A CN108318737A (zh) | 2018-01-22 | 2018-01-22 | 一种高精度数字鉴相方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108318737A true CN108318737A (zh) | 2018-07-24 |
Family
ID=62887607
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810060448.0A Pending CN108318737A (zh) | 2018-01-22 | 2018-01-22 | 一种高精度数字鉴相方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108318737A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109030936A (zh) * | 2018-07-27 | 2018-12-18 | 中国空间技术研究院 | 小型化相频测试仪 |
CN109115109A (zh) * | 2018-07-26 | 2019-01-01 | 西安工业大学 | 数字式调频连续波干涉激光驱动与信号处理方法及电路 |
CN110190634A (zh) * | 2019-05-30 | 2019-08-30 | 神驰机电股份有限公司 | 一种多台发电机组并联运行相位同步方法 |
CN111458558A (zh) * | 2020-04-03 | 2020-07-28 | 宁波奥克斯电气股份有限公司 | 一种参数获取方法、装置、电子设备及存储介质 |
CN111722013A (zh) * | 2020-06-24 | 2020-09-29 | 合肥工大高科信息科技股份有限公司 | 一种用于25hz相敏轨道电路信号检测的方法及*** |
CN112542823A (zh) * | 2020-11-05 | 2021-03-23 | 上海合凯电气科技有限公司 | 一种重合闸控制方法、***及一种重合闸控制设备 |
CN113504406A (zh) * | 2021-06-22 | 2021-10-15 | 昇辉控股有限公司 | 一种过零检测方法及具有过零检测功能的智能开关 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813725A (zh) * | 2010-04-15 | 2010-08-25 | 南京邮电大学 | 一种低频率信号的相位差测量方法 |
CN101871965A (zh) * | 2010-06-11 | 2010-10-27 | 威胜集团有限公司 | 电力正弦信号过零时间、频率、相位差的检测方法 |
CN201716368U (zh) * | 2009-06-11 | 2011-01-19 | 严海林 | 一种电缆相位分析识别装置 |
CN102315622A (zh) * | 2010-06-30 | 2012-01-11 | 日立汽车***株式会社 | 电能变换***以及电能变换装置 |
-
2018
- 2018-01-22 CN CN201810060448.0A patent/CN108318737A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201716368U (zh) * | 2009-06-11 | 2011-01-19 | 严海林 | 一种电缆相位分析识别装置 |
CN101813725A (zh) * | 2010-04-15 | 2010-08-25 | 南京邮电大学 | 一种低频率信号的相位差测量方法 |
CN101871965A (zh) * | 2010-06-11 | 2010-10-27 | 威胜集团有限公司 | 电力正弦信号过零时间、频率、相位差的检测方法 |
CN102315622A (zh) * | 2010-06-30 | 2012-01-11 | 日立汽车***株式会社 | 电能变换***以及电能变换装置 |
Non-Patent Citations (2)
Title |
---|
李敏等: "《信号分析与处理的软硬件实现(第2版)》", 28 February 2009, 大连海事大学出版社 * |
马昌凤等: "《最优化计算方法及其MATLAB程序实现》", 30 June 2015, 国防工业出版社 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109115109A (zh) * | 2018-07-26 | 2019-01-01 | 西安工业大学 | 数字式调频连续波干涉激光驱动与信号处理方法及电路 |
CN109030936A (zh) * | 2018-07-27 | 2018-12-18 | 中国空间技术研究院 | 小型化相频测试仪 |
CN109030936B (zh) * | 2018-07-27 | 2020-12-18 | 中国空间技术研究院 | 小型化相频测试仪 |
CN110190634A (zh) * | 2019-05-30 | 2019-08-30 | 神驰机电股份有限公司 | 一种多台发电机组并联运行相位同步方法 |
CN111458558A (zh) * | 2020-04-03 | 2020-07-28 | 宁波奥克斯电气股份有限公司 | 一种参数获取方法、装置、电子设备及存储介质 |
CN111722013A (zh) * | 2020-06-24 | 2020-09-29 | 合肥工大高科信息科技股份有限公司 | 一种用于25hz相敏轨道电路信号检测的方法及*** |
CN112542823A (zh) * | 2020-11-05 | 2021-03-23 | 上海合凯电气科技有限公司 | 一种重合闸控制方法、***及一种重合闸控制设备 |
CN112542823B (zh) * | 2020-11-05 | 2023-03-10 | 上海合凯电气科技有限公司 | 一种重合闸控制方法、***及一种重合闸控制设备 |
CN113504406A (zh) * | 2021-06-22 | 2021-10-15 | 昇辉控股有限公司 | 一种过零检测方法及具有过零检测功能的智能开关 |
CN113504406B (zh) * | 2021-06-22 | 2024-01-02 | 昇辉控股有限公司 | 一种过零检测方法及具有过零检测功能的智能开关 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108318737A (zh) | 一种高精度数字鉴相方法 | |
CN1996986B (zh) | 全相位时移相位差频谱校正法 | |
Brennan et al. | On the relationship between time and frequency domain methods in time delay estimation for leak detection in water distribution pipes | |
Kang et al. | Phase difference correction method for phase and frequency in spectral analysis | |
CN101806832A (zh) | 一种低频率信号的频率测量方法 | |
CN101907656B (zh) | 利用固定相移来测量同频信号相位差的方法 | |
CN104777375A (zh) | 一种互相关相位噪声检测装置 | |
Chen et al. | Robust precise time difference estimation based on digital zero-crossing detection algorithm | |
CN106154257B (zh) | 基于FFT与apFFT的精密测量雷达二次测频方法 | |
CN109990713A (zh) | 一种基于平面光栅激光干涉仪的高分辨率相位检测方法 | |
CN106093905B (zh) | 基于apFFT相位差的精密测量雷达频偏分类测频方法 | |
CN109030944A (zh) | 一种基于频移和dft的实正弦信号相位差估计方法 | |
CN105606893B (zh) | 基于空间平滑修正music的电力间谐波检测方法 | |
CN109521269A (zh) | 一种幅度调制信号数字化测频方法 | |
CN112541157B (zh) | 一种信号频率精确估计方法 | |
US9759751B1 (en) | Line cycle correlated spectral analysis for power measurement systems | |
CN106885941A (zh) | 基于频谱极值点的电网基波频率检测方法 | |
CN103575991B (zh) | 一种交流电的瞬时相位的精确测量方法 | |
RU2363005C1 (ru) | Способ спектрального анализа полигармонических сигналов и устройство для его реализации | |
TWI438445B (zh) | Measurement method of signal delay time | |
CN112595889B (zh) | 非理想多指数衰减正弦信号的欠Nyquist采样与参数测量方法 | |
CN110865577B (zh) | 一种用于交流电阻校准的数字采样方法及装置 | |
Kostiukov | Decreasing of Bias of Modified Cross-Correlation Method for Phase Shift Measurement | |
CN104407212A (zh) | 一种实时高精度相参脉冲信号测频方法 | |
CN114624513B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180724 |
|
WD01 | Invention patent application deemed withdrawn after publication |