CN103901419B - 一种基于频域相位校正的外辐射源雷达距离徙动补偿方法 - Google Patents
一种基于频域相位校正的外辐射源雷达距离徙动补偿方法 Download PDFInfo
- Publication number
- CN103901419B CN103901419B CN201410161409.1A CN201410161409A CN103901419B CN 103901419 B CN103901419 B CN 103901419B CN 201410161409 A CN201410161409 A CN 201410161409A CN 103901419 B CN103901419 B CN 103901419B
- Authority
- CN
- China
- Prior art keywords
- frequency domain
- phase correction
- echoed signal
- signal
- pulse
- 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/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
本发明涉及一种基于频域相位校正的外辐射源雷达距离徙动补偿方法,属于雷达目标探测技术领域。首先,先将回波信号进行分段处理;其次将回波变换到频域;然后在频域对回波进行相位校正,使得信号在带宽内可以同相叠加,消除距离徙动的影响;最后将相位校正后的回波再经过傅里叶逆变换到时域,与直达波信号进行互模糊函数计算,得到互模糊函数图。在跟踪模式下(目标速度已知)或目标大致速度已知情况下,本发明方法积累增益和Keystone变换相当,但运算量要小的多。
Description
技术领域
本发明涉及一种基于频域相位校正的外辐射源雷达距离徙动补偿方法,属于雷达目标探测技术领域。
背景技术
外辐射源雷达是一种双(多)基地雷达,它利用已有的非合作辐射源(如广播、电视信号等),通过被动接收目标散射的回波信号,实现对目标的检测,其配置如附图1所示。由于目标回波通常比较微弱,因此需要进行长时间相参积累以提高信噪比。为保证积累增益,要求在相参积累时间内,目标的距离徙动不能超过一个距离单元。然而由于外辐射源雷达目标回波十分微弱,需要的积累时间较长。在较长的积累时间内,目标运动距离会超过多个距离单元,即发生距离徙动现象,如附图2所示。距离徙动会造成目标峰值能量的分散和信噪比下降,进而影响***作用距离,因此需要寻找方法解决。
对于脉冲体制雷达来说,Keystone变换是一种解决距离徙动的有效方法。在赵永波的“一种用于弱信号检测的广义Keystone变换算法”一文中提到了脉冲雷达Keystone变换的原理。以点目标为例,目标的基带回波信号表示为:
其中tm=mT分别表示快时间和慢时间,T为脉冲重复周期,A为回波目标的幅度,R(tm)为tm时刻目标的距离,为归一化的回波包络,fc为载波频率。
对(1)式在快时间维进行傅里叶变换可得:
经过距离向的脉压后得到:
假设目标在tm里匀速运动,即速度v保持不变,则
上式中第一个指数项表示目标0时刻的位置,第二个指数项表示相位随时间tm的变化,由于f的存在,相位的变化率是不同的。
Keystone变换采用对慢时间进行尺度变换的方法:
将式(5)代入式(4)得到:
由式(6)可知,以t′m为新的时间来度量,相位的变化率是恒定的,即不同段的目标回波信号之间的包络不再有平移,从而解决了距离徙动问题。
Keystone变换的原理如图3所示。变换前的数据如图3(a)所示,图3(b)中的空心点表示变换后的数据点。由于变换后的数据需要进行IFFT运算变换回时域,所以需要采用SINC内插的方法调整变换后的数据的采样间隔,内插后的数据如图3(b)中的实心点所示。
Keystone变换应用于外辐射源雷达时能很好的解决距离徙动的问题,但是运算量非常大,需要寻找快速的计算方法。
发明内容
本发明的目的是为解决检测目标微弱回波时进行长时间积累而出现距离徙动,不利于目标检测的问题,提出一种基于频域相位校正的外辐射源雷达距离徙动补偿方法,在跟踪情况下(目标速度已知)或目标大致速度已知情况下,通过外辐射源雷达长时间相参积累消除目标距离徙动、提高目标积累增益。
本发明通过在频域对目标回波进行相位校正,使得信号在带宽内可以进行同相叠加,以达到提高目标检测性能的目的,与此同时还降低了运算量。
本发明的目的是通过以下技术方案实现的,具体包括以下步骤:
步骤一,数据分段
对回波信号进行数据分段处理,把要处理的回波信号分成N个脉冲,每个脉冲的数据点数为L,即要处理的数据总点数为NL。
步骤二,频域变换
把步骤一分段后的回波信号进行快速傅里叶变换,变换到频域:
E=FFT{e(n)}(7)
其中e(n)为分段后的回波信号,E为回波信号在频域的表示。
步骤三,相位校正
为了能够对信号进行同相叠加,需要使相频响应在信号带宽内为一条直线,因此对步骤二得到的回波信号在频域进行相位校正,相邻脉冲间相位校正量φ(k)为:
其中
式中各变量的定义为:
m:回波信号的第m个脉冲,m=1,2,...,N;
L:一个脉冲内的数据点数;
fs:回波信号的基带采样率;
fc:回波信号的载波频率;
θc:回波信号的相邻脉冲间对应信号中心频率处的相位变化量;
fd:目标的多普勒频率;
N:回波信号的脉冲个数。
在频域对回波进行相位校正,具体公式为:
E'=E·exp(jφ(k))(10)
其中E'为在频域进行相位校正补偿后的回波信号。
步骤四,频域变换到时域
通过快速傅里叶逆变换,把步骤三进行相位校正后的频域回波信号变换到时域:
e'(n)=IFFT{E'}(11)
其中e'(n)为经过相位校正后的回波信号。
有益效果
与传统的Keystone方法进行距离徙动补偿相比,相位校正方法的计算量小。具体运算量为:Keystone方法总的复乘次数为3N(Llog2(L)/2)+N(N+1)L,总的复加次数3N(Llog2(L))+N(N-1)L。而相位校正算法总的复乘次数为2N(Llog2(L)/2),总的复加次数为2N(Llog2(L))。本发明的相位校正方法明显比Keystone方法的运算量小得多。
附图说明
图1为背景技术中外辐射源雷达配置示意图;
图2为背景技术中长时间积累距离徙动现象;
图3为背景技术中Keystone变换原理示意图,其中(a)为Keystone变换前的数据,(b)为Keystone变换后的数据;
图4为具体实施方式中积累时间为1.0s的情况下,互模糊函数侧视图,其中(a)为频域相位校正算法,(b)为Keystone算法;
图5为具体实施方式中不同积累时间下,频域相位校正算法和Keystone算法运算量的比较。
具体实施方式
下面将结合附图及具体实施方法对本发明作进一步说明。
一种基于频域相位校正的外辐射源雷达距离徙动补偿方法,集成于外辐射源雷达信号处理模块中。外辐射源雷达一般有两副天线,参考天线接收直达波信号,回波天线接收目标反射的电视信号,所测目标为民航飞机,如附图1所示。
本实施例利用实验的手段来验证本发明的有效性。
实验参数:中央电视广播塔发射信号带宽为7.56MHz,载波频率为674MHz,基带采样率fs=10MHz,检测目标为民航飞机,其多普勒频率fd=-483Hz。
具体实施步骤为:
1)对回波信号进行数据分段处理。设积累时间为Ti,则回波的数据总点数为Ti·fs,每个脉冲的数据点数L=5000,即把回波分成N=Ti·fs/5000个脉冲。由此可得当积累时间为1.0s时,回波总点数为10M,每个脉冲数据点数L=5000,把回波分成N=2000个脉冲。
2)对回波信号进行快速傅里叶变换到频域。对每段回波信号进行5000点的FFT运算,得到回波信号的频域表示E=FFT{e(n)}。
3)相位校正。当积累时间为1.0s时,目标的多普勒频率fd=-483Hz,相邻脉冲间对应信号中心频率处的相位变化量θc=-0.483π。相邻脉冲间相位校正量为φ(k)=-(m-1)·k/337000·(-0.483π),其中k=-2500,-2499,...,2498,2499,m=1,2,...,2000。求得θc后频域对回波信号进行相位校正补偿,即E'=E·exp[-(m-1)·k·j/337000·(-0.483π)]。
4)把已经进行相位校正的频域回波信号变换到时域,即对频域的每段回波信号进行5000点的IFFT得到校正后的时域回波信号:e'(n)=IFFT{E'}。
在时域对已经进行相位校正的回波信号和直达波信号做互模糊函数得到互模糊函数图。当积累时间为1.0s时,互模糊函数图如附图4所示。采用频域相位校正算法目标的信噪比为17.1dB,总的复乘次数为1.2288×1010,总的复加次数为2.4575×108;采用Keystone算法目标的信噪比为16.9dB,总的复乘次数为2.0194×1010,总的复加次数为2.0359×1010。由此可见在相同的积累时间下,采用频域相位校正算法和采用Keystone算法目标的检测性能相当,但是采用频域相位校正算法运算量大大降低了。
当积累时间为1.5s、2.0s、2.5s、3.0s时,采用频域相位校正算法和采用Keystone算法目标的检测性能相当,但是采用频域相位校正算法运算量大大降低了。不同积累时间下,频域相位校正算法和Keystone算法运算量的比较如附图5所示。
以上所述仅为本发明的具体实施例,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种基于频域相位校正的外辐射源雷达距离徙动补偿方法,其特征在于:具体包括以下步骤:
步骤一,数据分段:
对回波信号进行数据分段处理,把要处理的回波信号分成N个脉冲,每个脉冲的数据点数为L,要处理的数据总点数为NL;
步骤二,频域变换:
把步骤一分段后的回波信号进行快速傅里叶变换,变换到频域:
E=FFT{e(n)}
其中e(n)为分段后的回波信号,E为回波信号在频域的表示;
步骤三,相位校正:
对步骤二得到的回波信号在频域进行相位校正,使相频响应在信号带宽内为一条直线,相邻脉冲间相位校正量φ(k)为:
其中
式中各变量的定义为:
m:回波信号的第m个脉冲,m=1,2,...,N;
L:一个脉冲内的数据点数;
fs:回波信号的基带采样率;
fc:回波信号的载波频率;
θc:回波信号的相邻脉冲间对应信号中心频率处的相位变化量;
fd:目标的多普勒频率;
N:回波信号的脉冲个数;
在频域对回波进行相位校正,具体公式为:
E'=E·exp(jφ(k))
其中E'为在频域进行相位校正后的回波信号;
步骤四,频域变换到时域:
通过快速傅里叶逆变换,把步骤三进行相位校正后的频域回波信号变换到时域:
e'(n)=IFFT{E'}
其中e'(n)为经过相位校正后的回波信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410161409.1A CN103901419B (zh) | 2014-04-22 | 2014-04-22 | 一种基于频域相位校正的外辐射源雷达距离徙动补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410161409.1A CN103901419B (zh) | 2014-04-22 | 2014-04-22 | 一种基于频域相位校正的外辐射源雷达距离徙动补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103901419A CN103901419A (zh) | 2014-07-02 |
CN103901419B true CN103901419B (zh) | 2016-03-30 |
Family
ID=50992869
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410161409.1A Expired - Fee Related CN103901419B (zh) | 2014-04-22 | 2014-04-22 | 一种基于频域相位校正的外辐射源雷达距离徙动补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103901419B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104535999A (zh) * | 2015-01-02 | 2015-04-22 | 中国人民解放军国防科学技术大学 | 一种校正天线方向图影响的雷达成像数据预处理方法 |
CN104714213A (zh) * | 2015-03-30 | 2015-06-17 | 宁波高新区宁源科技服务有限公司 | 一种基于多抽样率和速度通道划分的徒动补偿方法 |
CN107144825B (zh) * | 2017-04-25 | 2019-10-01 | 西安电子科技大学 | 一种非参数化数字电视外辐射源雷达多普勒扩散补偿方法 |
CN108398676B (zh) * | 2018-05-04 | 2021-10-26 | 电子科技大学 | 一种外辐射源雷达微弱运动目标探测方法 |
CN108896976A (zh) * | 2018-07-05 | 2018-11-27 | 电子科技大学 | 一种用于多频道外辐射源雷达的相参处理方法 |
CN108931768B (zh) * | 2018-07-05 | 2022-07-01 | 电子科技大学 | 一种基于信号融合的外辐射源微弱目标定位的方法 |
CN110632573B (zh) * | 2019-10-17 | 2023-01-31 | 南京航空航天大学 | 一种机载宽带雷达空时二维keystone变换方法 |
CN110943951B (zh) * | 2019-11-29 | 2022-03-29 | 慧众行知科技(北京)有限公司 | 一种用于数据分段解调的相位补偿方法 |
CN113030886B (zh) * | 2021-03-02 | 2022-10-18 | 上海无线电设备研究所 | 一种高速目标距离徙动校正方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5865794B2 (ja) * | 2012-07-13 | 2016-02-17 | 三菱電機株式会社 | レーダ装置 |
-
2014
- 2014-04-22 CN CN201410161409.1A patent/CN103901419B/zh not_active Expired - Fee Related
Non-Patent Citations (3)
Title |
---|
The migration compensation methods for DTV based passive radar;Yuan Feng 等;《2013 IEEE Radar Conference》;20131231;1-4 * |
外辐射源雷达长时间积累目标探测性能分析;张南 等;《电子学报》;20080630;第36卷(第6期);1103-1107 * |
数字电视辐射源雷达的相参积累徙动补偿方法;杨金禄 等;《电子与信息学报》;20110228;第33卷(第2期);407-411 * |
Also Published As
Publication number | Publication date |
---|---|
CN103901419A (zh) | 2014-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103901419B (zh) | 一种基于频域相位校正的外辐射源雷达距离徙动补偿方法 | |
CN106970371B (zh) | 一种基于Keystone和匹配滤波的目标检测方法 | |
CN104062640B (zh) | 一种外辐射源雷达距离徙动补偿的快速实现方法 | |
CN103616679B (zh) | 基于差波束调制和波形分析的pd雷达测距测角方法 | |
Xing et al. | New parameter estimation and detection algorithm for high speed small target | |
CN103412301B (zh) | 米波雷达测角方法 | |
CN102323575B (zh) | Pd雷达在微弱信号检测过程中的距离走动校正方法 | |
CN103675759B (zh) | 一种改进的分数阶傅里叶变换机动弱目标检测方法 | |
CN102360070B (zh) | 超宽带脉冲信号接收装置及超宽带脉冲雷达*** | |
CN106896358A (zh) | 一种基于位置旋转变换的高速目标相参积累检测方法 | |
CN103823217A (zh) | 基于双频发射的双基地mimo雷达高速运动目标参数估计方法 | |
CN102156279A (zh) | 基于mimo的双基地雷达地面动目标检测方法 | |
CN110109091B (zh) | 一种针对高速目标的无源雷达参数估计方法及装置 | |
CN105158759A (zh) | 基于杂波相消的hrws sar通道相位偏差校正方法 | |
CN108931768A (zh) | 一种基于多信号融合的外辐射源微弱目标定位的方法 | |
CN104155647A (zh) | 基于机扫米波阵列雷达的估计目标方位角的方法 | |
CN104849708A (zh) | 基于频域多项式相位变换的高速机动目标参数估计方法 | |
CN103605112A (zh) | 多发多收干涉合成孔径雷达时频二维信号波形设计方法 | |
CN103116162A (zh) | 基于目标空间稀疏性的高分辨声呐定位方法 | |
CN104280720B (zh) | 前视机载雷达发射方向图的设计方法 | |
Zhang et al. | Joint high‐resolution range and DOA estimation via MUSIC method based on virtual two‐dimensional spatial smoothing for OFDM radar | |
CN113885002A (zh) | 基于尺度变换的外辐射源雷达捷变频信号相参积累方法 | |
Sun et al. | Underwater weak moving target detection method based on wideband multi-pulse coherent integration | |
Yang et al. | Implementation architecture of signal processing in pulse Doppler radar system based on FPGA | |
Wu et al. | Efficient parameter estimation method for maneuvering targets in discrete randomly-modulated radar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160330 Termination date: 20170422 |
|
CF01 | Termination of patent right due to non-payment of annual fee |