CN104020450B - 一种基于邻波相关和相位分层的雷达回波重构方法 - Google Patents
一种基于邻波相关和相位分层的雷达回波重构方法 Download PDFInfo
- Publication number
- CN104020450B CN104020450B CN201410121222.9A CN201410121222A CN104020450B CN 104020450 B CN104020450 B CN 104020450B CN 201410121222 A CN201410121222 A CN 201410121222A CN 104020450 B CN104020450 B CN 104020450B
- Authority
- CN
- China
- Prior art keywords
- phase
- echo
- sequence
- signal
- phase difference
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明为一种基于邻波相关和相位分层的雷达回波重构方法,该方法是通过邻波相关进行相位对齐,对相邻的回波脉冲进行采样,将得到的采样信号进行相关运算,根据相关结果的峰值实现回波相位预对齐;建立回波采样序列的相位差序列并进行线性拟合,进而将相位差序列分层为具有相同斜率的多个组合;根据分层结果对雷达回波进行筛选并重构,利用最终的重构雷达回波开展频率估计及速度解算。本方法用于解决传统技术中发射信号与回波相关而产生的相位对不齐以及由于目标运动状态或信号转发机制而产生的回波相位跳跃问题,降低了目标机动大时的测速误差。
Description
技术领域
本发明涉及雷达测量领域,具体是一种合作式脉冲雷达测速机的测量回波信号处理方法。
背景技术
传统的合作式脉冲雷达测速机的测速原理如图1所示,脉冲雷达测速机在处理接收回波获取速度值时,以发射信号St(n)为参考信号,将回波Sr(n)与参考信号St(n)进行相关匹配滤波,匹配滤波的结果作为速度估计的回波数据,通过相参积累一段时间的回波数据后进行多普勒频率估计,进而使用αβ滤波器进行速度估计。传统测速机对机动小的目标展示了良好的测速性能,但是传统测速机参考信号固定为发射信号St(n),通过运动目标反射之后的回波Sr(n)与St(n)相比,相当于进行了运动调制,当目标机动较大时,如图3所示,Sr(n)与St(n)进行相关后的相参性会有所损失,而同时受目标运动状态或信号转发机制的影响,回波Sr(n)的相位在目标机动较大时会出现阶梯式跳跃的现象,如图5所示,最终导致测速结果不准确,其测速误差率如图7所示。
发明内容
本发明的目的即是给出一种解决相位对不齐和相位跳跃的雷达回波重构方法,采用的技术方案如下:
步骤一、通过邻波相关进行相位对齐,对相邻的回波脉冲进行采样,将得到的采样信号进行相关运算,根据相关值的模值序列的最大值所对应的索引实现采样序列的相位对齐;
步骤二、相位分层的实现,积累m个相位对齐后的回波采样序列,根据加速度在短时间内可近似为恒定值的特性,求解这m个回波脉冲采样序列的相位差序列并进行线性拟合,将相位差序列分层为具有相同斜率的多个组合,求解各组合的拟合系数,实现最终的相位分层;其中m=1,2...26;
步骤三、雷达回波重构的实现
(1)根据得到的所含元素最多的组合对应的拟合系数对相位差序列和回波脉冲采样序列进行重构;
(2)对于m+j个回波信号,找到第j个回波序列信号所属的相位差集合,去除第j个回波信号的采样序列,保持循环序列中回波信号的个数始终为m个,将第m+j个相位差的值与各组合根据拟合系数预报至m+j的值之差与门限值进行比较,判断其所属的分层,根据分层结果对相位差序列和回波脉冲采样序列进行重构;其中j为正整数,且j≥1;
(3)如果(1)中选定的回波信号所属线性拟合集合的中所含相位差序列的个数减少至零,则回到步骤二重新启动此流程。
发明的有益效果:
1.使用邻波相关的方法,解决了相位预对齐的问题,初步保证了后续处理过程中的相参性;
2.使用相位差分层的方法,解决了回波信号相位中出现的阶梯式跳跃问题,进一步保证了相参的准确性;
3.综合运用上述两种方法,解决了回波信号相位未对齐及相位跳变引起的测速不准确问题,降低了目标机动大时的测速误差。
附图说明
图1是传统测速机原理框图;
图2是改进后测速机原理框图;
图3是相位对齐前回波脉冲信号;
图4是相位对齐后回波脉冲信号;
图5是原始回波采样信号相位差分层示例;
图6是重构的回波采样信号相位差示例;
图7是改进前测速误差;
图8是改进后测速误差;
图9相位分层及重构流程图。
具体实施方式
本发明的目的即是给出一种解决相位对不齐和相位跳跃的雷达回波重构方法,具体的改进流程如图9所示,下面就结合附图进行进一步的阐述。
本发明的测速原理如图2所示,是以回波Sr(n-1)为参考信号,将回波信号Sr(n)与Sr(n-1)进行相关运算,根据相关结果的峰值实现Sr(n)序列的相位预对齐,如图4所示,建立Sr(n)的相位差序列θ(n),分析相位差信息θ(n),根据加速度在短时间内可近似为恒定值的特性,对θ(n)进行线性拟合,进而将θ(n)分层为具有相同斜率的多个组合,根据分层结果对雷达回波进行筛选并重构,如图6所示,利用最终的重构雷达回波开展频率估计及速度计算,最终达到提高测速精度的目的,改进后的测速误差率如图8所示。具体实现步骤如下:
步骤一:通过邻波相关进行相位对齐
对第n和n-1个回波脉冲Sr(n)和Sr(n-1)进行采样,得到采样信号xn与yn,其中n是正整数,n=1...150,将xn与yn进行相关运算得到相关值Rxy(t)(t=-149,...,0,...149),表达式为:
由于相关值Rxy(t)为复数,令其模值为abs(Rxy(t)),根据模值序列abs(Rxy(t))(t=-149,...,0,...149)的最大值对应的索引实现xn与yn采样序列的相位对齐。
步骤二:相位分层的实现
(1)求解回波脉冲采样序列的相位差序列
针对相位对齐后的回波脉冲采样序列Sr(m,p),首先积累m个采样,对该m个采样进行分层与重构,之后对第m+1及以后的每个采样,进行逐点分层和重构。
假设相位对齐后的回波脉冲采样序列为Sr(m,p),其中m=1,2...26;p=1,2,3...150,为了使输出信噪比尽量高,取p=150/2,则Sr(m,150/2)(m=1,2,3...26),为表述方便简化表示为Sr(m),令θ(m)表示Sr(m)的相位差序列,单位为度,则Sr(m)的相位差序列可表示为
θ(m)=deg(Sr(m)*Sr*(m-1)) (3)
其中deg表示求Sr(m)与Sr(m-1)共轭乘积的角度值。
(2)对θ(m)进行线性拟合
θ(mi)=k*mi+bi (4)
其中(mi∈[1,2,3...26];i=1,2,3,4),k和bi分别是线性拟合的拟合系数,k为斜率和bi为斜距值。
令Cbi(i)表示满足θ(mi)=k*mi+bi的θ(m)集合,则选择包含θ(m)最多的Cbi(i)对应的k和bi作为回波信号的最终拟合系数k和blast,由此实现最终的相位分层。
步骤三、雷达回波重构的实现
(1)根据得到的所含元素最多的组合对应的拟合系数对相位差序列θ(m)和回波脉冲采样序列Sr(m)进行重构
其中,为进行重构后的相位差序列,为重构后的回波脉冲采样序列,m=1,2,3...26。
(2)对于下一个回波信号Sr(m+j),判断其所属的分层组合:
对于m+j个回波信号,找到第j个回波序列信号所属的相位差集合,去除第j个回波信号的采样序列,保持循环序列中回波信号的个数始终为m个,将第m+j个相位差的值与各组合根据拟合系数预报至m+j的值之差与门限值进行比较,判断其所属的分层,根据分层结果对相位差序列和回波脉冲采样序列进行重构;
即k*(m+j)+bi-Sr(m+j)<门限值,则Sr(m+1)∈Cbi(i)
其中j为正整数,且j≥1;
为保证测试结果,循环序列中回波信号的个数保持不变,对回波信号Sr(m+j),则找到Sr(j)所属的Cbi(i),将Sr(j)去掉,重新对该组进行系数拟合及重构。根据步骤四中的第(1)步重构Sr(m+j)所在组的如果选定的回波信号所属的Cbi(i)中所含θ(m)的个数为零,则回到步骤二重新启动此流程。
Claims (2)
1.一种基于邻波相关和相位分层的雷达回波重构方法,其特征包括以下步骤:
步骤一、通过邻波相关进行相位对齐
对相邻的回波脉冲进行采样,将得到的采样信号进行相关运算,根据相关值的模值序列的最大值所对应的索引实现采样序列的相位对齐;
对第n和n-1个回波脉冲Sr(n)和Sr(n-1)进行采样,得到采样信号xn与yn,其中n是正整数,n=1...150,将xn与yn进行相关运算得到相关值Rxy(t)(t=-149,...,0,...149),表达式为:
由于相关值Rxy(t)为复数,令其模值为abs(Rxy(t)),根据模值序列abs(Rxy(t))(t=-149,...,0,...149)的最大值对应的索引实现xn与yn采样序列的相位对齐;
步骤二、相位分层的实现
积累m个相位对齐后的回波采样序列,根据加速度在短时间内可近似为恒定值的特性,求解这m个回波脉冲采样序列的相位差序列并进行线性拟合,将相位差序列分层为具有相同斜率的多个组合,求解各组合的拟合系数,实现最终的相位分层;
步骤三、雷达回波重构的实现
(1)选择步骤二中得到的所含元素最多的组合对应的拟合系数,对相位差序列和回波脉冲采样序列进行重构;
(2)对于m+j个回波信号,找到第j个回波序列信号所属的相位差集合,去除第j个回波信号的采样序列,保持循环序列中回波信号的个数始终为m个,将第m+j个相位差的值与各组合根据拟合系数预报至m+j的值之差与门限值进行比较,判断其所属的分层,根据分层结果对相位差序列和回波脉冲采样序列进行重构;
(3)如果(1)中选定的回波信号所属线性拟合集合的中所含相位差序列的个数减少至零,则回到步骤二重新启动此流程。
2.如权利要求1所述的一种基于邻波相关和相位分层的雷达回波重构方法,其特征是:m=1,2...26;j为正整数,且j≥1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410121222.9A CN104020450B (zh) | 2014-03-26 | 2014-03-26 | 一种基于邻波相关和相位分层的雷达回波重构方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410121222.9A CN104020450B (zh) | 2014-03-26 | 2014-03-26 | 一种基于邻波相关和相位分层的雷达回波重构方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104020450A CN104020450A (zh) | 2014-09-03 |
CN104020450B true CN104020450B (zh) | 2019-01-08 |
Family
ID=51437308
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410121222.9A Expired - Fee Related CN104020450B (zh) | 2014-03-26 | 2014-03-26 | 一种基于邻波相关和相位分层的雷达回波重构方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104020450B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107783077B (zh) * | 2016-08-25 | 2020-11-17 | 大连楼兰科技股份有限公司 | 过门限峰值点的处理的方法 |
CN112098927B (zh) * | 2020-08-19 | 2022-06-21 | 中国电子科技集团公司第二十九研究所 | 一种基于信号合成的多信道测向方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10241463A1 (de) * | 2002-09-06 | 2004-03-18 | Robert Bosch Gmbh | Radarmessvorrichtung, insbesondere für ein Kraftfahrzeug, und Verfahren zum Betreiben einer Radarmessvorrichtung |
CN102608603A (zh) * | 2012-03-13 | 2012-07-25 | 北京航空航天大学 | 一种基于完全互补序列的多通道合成孔径雷达成像方法 |
CN103389492A (zh) * | 2013-07-25 | 2013-11-13 | 西安电子科技大学 | 多通道随机谐波调制采样雷达接收机及其方法 |
CN103413296A (zh) * | 2013-07-13 | 2013-11-27 | 西安电子科技大学 | 在特征域内对极化sar目标的检测方法 |
-
2014
- 2014-03-26 CN CN201410121222.9A patent/CN104020450B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10241463A1 (de) * | 2002-09-06 | 2004-03-18 | Robert Bosch Gmbh | Radarmessvorrichtung, insbesondere für ein Kraftfahrzeug, und Verfahren zum Betreiben einer Radarmessvorrichtung |
CN102608603A (zh) * | 2012-03-13 | 2012-07-25 | 北京航空航天大学 | 一种基于完全互补序列的多通道合成孔径雷达成像方法 |
CN103413296A (zh) * | 2013-07-13 | 2013-11-27 | 西安电子科技大学 | 在特征域内对极化sar目标的检测方法 |
CN103389492A (zh) * | 2013-07-25 | 2013-11-13 | 西安电子科技大学 | 多通道随机谐波调制采样雷达接收机及其方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104020450A (zh) | 2014-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109061589B (zh) | 随机跳频雷达的目标运动参数估计方法 | |
CN103630888B (zh) | 基于对称三角lfmcw雷达的高精度实时微波测速测距装置 | |
CN110412558A (zh) | 基于tdm mimo的解车载fmcw雷达速度模糊方法 | |
CN107132534B (zh) | 一种高速雷达目标频域检测的优化方法 | |
CN106597429A (zh) | 一种基于雷达信号处理系统的速度跟踪子系统 | |
CN106814353A (zh) | 一种雷达信号处理系统 | |
CN107678003A (zh) | 一种地波雷达海杂波背景下的目标检测方法和装置 | |
CN104502898B (zh) | 将修正rft和修正mdcft相结合的机动目标参数估计方法 | |
CN106154265B (zh) | 一种基于频域移位的去斜体制雷达isar包络对齐方法 | |
CN102680948A (zh) | 一种线性调频信号调频率和起始频率估计方法 | |
CN104239683B (zh) | 一种基于改变信号调频率的解距离—速度模糊方法 | |
CN113093120B (zh) | 基于capon算法的PRI捷变雷达目标参数估计方法 | |
CN106646451A (zh) | 一种基于雷达信号处理系统的角度跟踪子系统 | |
CN110488241A (zh) | 一种激光雷达多次回波信息提取装置及方法 | |
CN102121990B (zh) | 基于空时分析的逆合成孔径雷达的目标转速估计方法 | |
CN104020450B (zh) | 一种基于邻波相关和相位分层的雷达回波重构方法 | |
CN114553260B (zh) | 一种ds/fh扩频信号载波频率的高精度测量系统 | |
CN108809354B (zh) | 一种高动态多普勒捕获的长时间累积方法 | |
CN107390213A (zh) | 一种基于滑动时窗的探地雷达记录剖面的时延曲线提取方法 | |
CN108333568A (zh) | 冲击噪声环境下基于Sigmoid变换的宽带回波Doppler和时延估计方法 | |
CN105445767A (zh) | 基于平均模糊函数的boc信号参数盲估计方法 | |
CN202661624U (zh) | 基于fpga的激光多普勒雷达信号处理器 | |
CN104765040A (zh) | 单脉冲波形识别提取方法 | |
CN105849584A (zh) | 雷达装置以及距离速度测量方法 | |
CN106027454A (zh) | 基于频率补偿的qam信号的频偏估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 | ||
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: 20190108 Termination date: 20200326 |