CN112748285A - 一种基于智能跟踪相关运算的相位测量方法 - Google Patents
一种基于智能跟踪相关运算的相位测量方法 Download PDFInfo
- Publication number
- CN112748285A CN112748285A CN202011521929.0A CN202011521929A CN112748285A CN 112748285 A CN112748285 A CN 112748285A CN 202011521929 A CN202011521929 A CN 202011521929A CN 112748285 A CN112748285 A CN 112748285A
- Authority
- CN
- China
- Prior art keywords
- path
- intelligent tracking
- bandwidth
- paths
- phase measurement
- 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
Images
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)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于智能跟踪相关运算的相位测量方法,属于电子对抗领域。该方法为:首先构建N路采集数据向量,并测量第1路信号的中心频率和带宽;其次根据第1路信号的带宽智能匹配,生成一组带宽适配的低通滤波器系数,并根据第1路信号的中心频率,构建智能跟踪滤波向量;接着N路采集数据向量分别与智能跟踪滤波向量卷积,得到N路智能跟踪数据向量;然后将第2到N路的智能跟踪数据向量分别与第1路智能跟踪数据向量做相关运算;最后分别提取相关结果的相位信息,即为第2到N路与第1路的相位差。
Description
技术领域
本发明属于高精度相位测量技术,具体涉及一种基于智能跟踪相关运算的相位测量方法。
背景技术
在电子对抗领域,相位作为一种典型的信号特征参数,在目标测向定位、辐射源个体识别等方向发挥了非常重要的作用。目前,国内相位测量系统大多基于模拟电路设计,这类相位测量系统往往精度低、体积大、成本高。而得益于大规模数字集成电路的发展,数字化相位测量技术日益成熟,其精度高、系统稳定性好等优点使其获得了广泛的应用。
数字化相位测量技术主要有快速傅立叶变换(FFT)和相关运算这两种方法。快速傅立叶变换理论基于多路弦波信号的相位差等于各自采样序列的离散频谱在最大谱线处的相位差,方法简单,但对于相位调制、线性调频等宽带信号,损失了大量频谱能量导致精度大幅降低;相关运算法基于多路弦波信号的互相关函数和自相关函数,通过反余弦求解多路信号的相位差值,有效利用了宽带内频谱能量,但宽带内的噪声导致精度下降。
针对快速傅立叶变换法和相关运算法精度不足的问题,本发明提出一种基于智能跟踪相关运算的相位测量方法。该方法能够根据原始中频数据的中心频率和带宽智能匹配跟踪,充分利用通带内有效频谱能量,同时基于跟踪滤波去除通带内的噪声,从而有效提升信噪比,提高相位测量精度。本方案原理简单,相位测量精度高,还可以利用FPGA的高速并行运算能力,实现对逐脉冲信号的实时处理,从而实现单脉冲的高精度相位测量功能。
发明内容
本发明的目的在于提供一种基于智能跟踪相关运算的相位测量方法,针对宽带信号和低信噪比条件,传统的快速傅立叶变换、相关运算等算法相位测量精度较低。本发明能够根据原始中频数据的中心频率和带宽智能匹配跟踪,充分利用通带内有效频谱能量,同时去除带内噪声影响,从而有效提升信噪比,提高相位测量精度,实现对信号的逐脉冲高精度相位测量功能。
实现本发明目的的技术解决方案为:一种基于智能跟踪相关运算的相位测量方法,包括以下步骤:
步骤1、构建N路采集数据向量,并测量第1路信号的中心频率和带宽;
步骤2、根据第1路信号的带宽智能匹配,生成一组带宽适配的低通滤波器系数,并根据第1路信号的中心频率,构建智能跟踪滤波向量;
步骤3、N路采集数据向量分别与智能跟踪滤波向量卷积,得到N路智能跟踪数据向量;
步骤4、将第2到N路的智能跟踪数据向量分别与第1路智能跟踪数据向量做相关运算;
步骤5、分别提取相关结果的相位信息,即为第2到N路与第1路的相位差。
本发明与现有技术相比,其显著优点在于:
(1)本发明原理简单、易于实现,得益于FPGA的高速并行运算能力,相位测量精度高、系统稳定性好,可实现对信号的逐脉冲高精度相位测量功能;
(2)根据检测信号的频率和带宽特性,智能匹配跟踪滤波器的中心频率和通带带宽,充分利用通带内有效频谱能量,同时去除通带内的噪声,从而有效提升信噪比,提高相位测量精度。
附图说明
图1为本发明所述的一种基于智能跟踪相关运算的相位测量方法流程图。
图2为实施例中常规信号相位测量误差对比图。
图3为实施例中相位调制信号相位测量误差对比图。
图4为实施例中线性调频信号相位测量误差对比图。
具体实施方式
下面结合附图对本发明作进一步详细描述。
结合图1,本发明的一种基于智能跟踪相关运算的相位测量方法,包括以下步骤:
步骤1、在采样率fs下,采集一定带宽内的原始中频数据,共计采集N路,每路M点,构建采集数据向量Xn如下:
Xn=[xn(1) xn(2) … xn(m) … xn(M)],1≤n≤N,1≤m≤M (1)
测量采集的第1路原始中频数据X1的中心频率为f,带宽为B。
步骤2、根据第1路原始中频数据X1的带宽B智能匹配,生成一组带宽适配的低通滤波器系数,并根据第1路中心频率f,构建智能跟踪滤波向量HTrack,具体包括以下步骤:
步骤2-1、根据第1路原始中频数据X1的带宽B智能匹配,生成一组带宽适配的K阶低通滤波器系数hLowpass如下:
hLowpass=[h(1) h(2) … h(k) … h(K)],1≤k≤K (2)
步骤2-2、根据测得第1路原始中频数据X1的中心频率f对低通滤波器系数hLowpass进行频率搬移,构建智能跟踪滤波向量HTrack如下:
其中1≤k≤K。
步骤3、N路采集数据向量Xn分别与智能跟踪滤波向量HTrack卷积,得到N路智能跟踪数据向量Yn,具体如下:
Yn=HTrack*Xn,1≤n≤N (4)
步骤4、将第2路到第N路的智能跟踪数据向量分别与第1路智能跟踪数据向量做相关运算,得到相关结果Rn如下:
Rn=Yn·Y′1,2≤n≤N (5)
其中,Y′1为Y1的转置。
实施例:
一种基于智能跟踪相关运算的相位测量方法,方法步骤如下:
步骤1、在采样率fs=1000MHz条件下,分别采集窄带常规、带宽为5MHz的相位调制以及带宽为5MHz的线性调频这三种类型的原始中频信号,共计采集N=2路,每路M=512点,以此构建采集数据向量X1和X2;测量采集的第1路原始中频数据的中心频率f=560MHz,带宽分别为窄带、B=5MHz和B=5MHz。
步骤2、根据第1路原始中频数据X1的带宽B=5MHz智能匹配,生成一组通带带宽为6MHz的K=20阶低通滤波器系数hLowpass,并根据第1路原始中频数据的信号频率f=560MHz构建智能跟踪滤波向量HTrack。
步骤3、将2路采集数据向量X1和X2分别与智能跟踪滤波向量HTrack卷积,得到2路智能跟踪数据向量Y1和Y2。
步骤4、将第2路智能跟踪数据向量Y2与第1路智能跟踪数据向量Y1做相关运算,得到相关结果R2。
每次随机产生一组相位差并以此生成2路原始中频数据,重复以上步骤,采用蒙特卡洛仿真计算相位测量误差的均方根结果。蒙特卡洛法的模拟次数为10000时,常规、相位调制、线性调频这三种类型信号的相位测量误差均方根结果如下:
(1)常规信号
采集常规信号的相位测量误差均方根结果如图2所示。
(2)相位调制信号
采集相位调制信号的相位测量误差均方根结果如图3所示。
(3)线性调频信号
采集线性调频信号的相位测量误差均方根结果如图4示。
试验结果表明,针对相位调制、线性调频等宽带信号时,快速傅立叶变换法由于损失了大量频谱能量导致相位测量精度大幅降低;而在低信噪比条件下,相关运算法由于宽带内噪声大,相位测量精度也受到较大影响。本专利提出的基于跟踪相关运算的相位测量方法,能够根据原始中频数据的中心频率和带宽智能适配跟踪,充分利用通带内有效频谱能量,同时基于跟踪滤波去除通带内的噪声,从而有效提升信噪比,提高相位测量精度。
综上所述,本专利提出的基于跟踪相关运算的相位测量方法,针对宽带信号、低信噪比条件,仍可保持较高的相位测量精度,满足了电子对抗领域对高精度相位测量的需求,具有广泛的应用价值。
Claims (6)
1.一种基于智能跟踪相关运算的相位测量方法,其特征在于,包括以下步骤:
步骤1、构建N路采集数据向量,并测量第1路信号的中心频率和带宽;
步骤2、根据第1路信号的带宽智能匹配,生成一组带宽适配的低通滤波器系数,并根据第1路信号的中心频率,构建智能跟踪滤波向量;
步骤3、N路采集数据向量分别与智能跟踪滤波向量卷积,得到N路智能跟踪数据向量;
步骤4、将第2到N路的智能跟踪数据向量分别与第1路智能跟踪数据向量做相关运算;
步骤5、分别提取相关结果的相位信息,即为第2到N路与第1路的相位差。
2.根据权利要求1所述的基于智能跟踪相关运算的相位测量方法,其特征在于,步骤1中,构建N路采集数据向量,并测量第1路信号的中心频率和带宽,具体如下:
在采样率fs下,采集一定带宽内的原始中频数据,共计采集N路,每路M点,构建采集数据向量Xn如下:
Xn=[xn(1) xn(2) … xn(m) … xn(M)],1≤n≤N,1≤m≤M (1)
测量采集的第1路原始中频数据X1的中心频率为f,带宽为B。
3.根据权利要求1所述的基于智能跟踪相关运算的相位测量方法,其特征在于,步骤2中,根据第1路原始中频数据X1的带宽B智能匹配,生成一组带宽适配的低通滤波器系数,并根据第1路中心频率f,构建智能跟踪滤波向量HTrack,具体包括以下步骤:
步骤2-1、根据第1路原始中频数据X1的带宽B智能匹配,生成一组带宽适配的K阶低通滤波器系数hLowpass如下:
hLowpass=[h(1) h(2) … h(k) … h(K)],1≤k≤K (2)
步骤2-2、根据测得第1路原始中频数据X1的中心频率f对低通滤波器系数hLowpass进行频率搬移,构建智能跟踪滤波向量HTrack如下:
其中1≤k≤K。
4.根据权利要求1所述的基于智能跟踪相关运算的相位测量方法,其特征在于,步骤3中,N路采集数据向量Xn分别与智能跟踪滤波向量HTrack卷积,得到N路智能跟踪数据向量Yn,具体如下:
Yn=HTrack*Xn,1≤n≤N (4)。
5.根据权利要求1所述的基于智能跟踪相关运算的相位测量方法,其特征在于,步骤4中,将第2路到第N路的智能跟踪数据向量分别与第1路智能跟踪数据向量做相关运算,得到相关结果Rn,具体如下:
Rn=Yn·Y′1,2≤n≤N (5)
其中,Y′1为Y1的转置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011521929.0A CN112748285B (zh) | 2020-12-21 | 2020-12-21 | 一种基于智能跟踪相关运算的相位测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011521929.0A CN112748285B (zh) | 2020-12-21 | 2020-12-21 | 一种基于智能跟踪相关运算的相位测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112748285A true CN112748285A (zh) | 2021-05-04 |
CN112748285B CN112748285B (zh) | 2023-04-28 |
Family
ID=75648629
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011521929.0A Active CN112748285B (zh) | 2020-12-21 | 2020-12-21 | 一种基于智能跟踪相关运算的相位测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112748285B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117289202A (zh) * | 2023-11-27 | 2023-12-26 | 中国航天科工集团八五一一研究所 | 一种自适应相位差测量方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6239743B1 (en) * | 1999-07-28 | 2001-05-29 | Trimble Navigation Limited | Integrated split spectrum positioning system receiver |
CN101750198A (zh) * | 2008-11-28 | 2010-06-23 | 上海宝钢工业检测公司 | 旋转机械不同测点振动信号相对相位的测量方法 |
CN103105533A (zh) * | 2011-11-09 | 2013-05-15 | 华北电力科学研究院有限责任公司 | 一种次同步振荡智能分析方法及系统 |
WO2014061878A1 (ko) * | 2012-10-19 | 2014-04-24 | 한국전기연구원 | Lpn 필터를 이용한 전력계통의 위상추종 시스템 |
CN111239781A (zh) * | 2020-02-28 | 2020-06-05 | 同济大学 | 一种集群自主协同中的gnss并行基带信号处理方法 |
-
2020
- 2020-12-21 CN CN202011521929.0A patent/CN112748285B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6239743B1 (en) * | 1999-07-28 | 2001-05-29 | Trimble Navigation Limited | Integrated split spectrum positioning system receiver |
CN101750198A (zh) * | 2008-11-28 | 2010-06-23 | 上海宝钢工业检测公司 | 旋转机械不同测点振动信号相对相位的测量方法 |
CN103105533A (zh) * | 2011-11-09 | 2013-05-15 | 华北电力科学研究院有限责任公司 | 一种次同步振荡智能分析方法及系统 |
WO2014061878A1 (ko) * | 2012-10-19 | 2014-04-24 | 한국전기연구원 | Lpn 필터를 이용한 전력계통의 위상추종 시스템 |
CN111239781A (zh) * | 2020-02-28 | 2020-06-05 | 同济大学 | 一种集群自主协同中的gnss并行基带信号处理方法 |
Non-Patent Citations (4)
Title |
---|
SANG HWAN PARK 等: "Manoeuvre detection and target tracking using state-space optimal FIR filters", 《PROCEEDINGS OF THE 1999 AMERICAN CONTROL CONFERENCE (CAT. NO. 99CH36251)》 * |
王显海 等: "一种基于IWOLA的脉冲星消色散数字终端设计", 《东南大学学报(自然科学版)》 * |
王铁丹 等: "基于二阶锥优化的复系数FIR滤波器设计", 《航天电子对抗》 * |
马帅 等: "基于强跟踪滤波器的单相电压基波相位实时提取", 《中国电机工程学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117289202A (zh) * | 2023-11-27 | 2023-12-26 | 中国航天科工集团八五一一研究所 | 一种自适应相位差测量方法 |
CN117289202B (zh) * | 2023-11-27 | 2024-02-13 | 中国航天科工集团八五一一研究所 | 一种自适应相位差测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112748285B (zh) | 2023-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102739272B (zh) | 一种基于fpga实现的信道化接收机子信道实时频谱合成方法 | |
CN102546499B (zh) | 一种实线性调频信号的分数阶信道化接收方法 | |
CN106771905B (zh) | 一种适用于高频电流局部放电检测的脉冲提取方法 | |
CN109490862B (zh) | 一种基于相位差分统计谱的载频估计方法 | |
CN111693944A (zh) | 雷达有源干扰信号参数提取与干扰样式识别方法及装置 | |
CN109975771B (zh) | 基于信号三阶相位差分的宽带数字信道化方法 | |
CN108710027B (zh) | 通道间相位差、幅度差的高精度测量方法 | |
CN101881821B (zh) | 一种分数阶傅里叶域信道化接收方法 | |
CN112162152B (zh) | 基于相位直线拟合的正弦波相参脉冲串信号频率估计方法 | |
CN111159888A (zh) | 一种基于互相关函数的协方差矩阵稀疏迭代时延估计方法 | |
CN107561357A (zh) | 一种基于信道化的高精度瞬时测频方法和装置 | |
CN111562438B (zh) | 一种基于fft和相位差的正弦信号频率估计方法和装置 | |
CN112748285A (zh) | 一种基于智能跟踪相关运算的相位测量方法 | |
CN105044769B (zh) | 提高地震信号的分辨率的方法 | |
Ibrahim et al. | Signal processing considerations for laser Doppler and phase Doppler applications | |
CN108132383B (zh) | 一种时间交错的多谐波信号欠采样方法 | |
CN102790650B (zh) | 一种宽带lfm信号的分数域奇偶交替信道化接收方法 | |
CN113406386B (zh) | 一种基于数字下变频的信号频率精确估计方法 | |
CN112162153B (zh) | 基于相位直线拟合的正弦波信号频率估计方法 | |
CN115510786A (zh) | 一种高动态冲击信号的频率解算方法与系统 | |
CN105700019A (zh) | 一种基于Born-Jordan时频分布的地震信号时频峰值滤波方法 | |
CN113627398B (zh) | 一种基于自适应重构滤波的信号特征检测方法 | |
CN112883787B (zh) | 一种基于频谱匹配的短样本低频正弦信号参数估计方法 | |
CN108919202A (zh) | 一种基于认知机制的非均匀动态滤波器组的构建方法 | |
Jayarama et al. | Enhanced Wideband Frequency Estimation via FFT: Leveraging Polynomial Interpolation and Array Indexing |
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 |