CN104062640B - 一种外辐射源雷达距离徙动补偿的快速实现方法 - Google Patents

一种外辐射源雷达距离徙动补偿的快速实现方法 Download PDF

Info

Publication number
CN104062640B
CN104062640B CN201410306628.4A CN201410306628A CN104062640B CN 104062640 B CN104062640 B CN 104062640B CN 201410306628 A CN201410306628 A CN 201410306628A CN 104062640 B CN104062640 B CN 104062640B
Authority
CN
China
Prior art keywords
time dimension
signal
czt
dimensional matrix
fast
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
Application number
CN201410306628.4A
Other languages
English (en)
Other versions
CN104062640A (zh
Inventor
单涛
陶然
冯远
张志鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201410306628.4A priority Critical patent/CN104062640B/zh
Publication of CN104062640A publication Critical patent/CN104062640A/zh
Application granted granted Critical
Publication of CN104062640B publication Critical patent/CN104062640B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S13/581Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of interrupted pulse modulated waves and based upon the Doppler effect resulting from movement of targets
    • G01S13/582Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of interrupted pulse modulated waves and based upon the Doppler effect resulting from movement of targets adapted for simultaneous range and velocity measurements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/285Receivers
    • G01S7/292Extracting wanted echo-signals
    • G01S7/2923Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/285Receivers
    • G01S7/288Coherent receivers
    • G01S7/2883Coherent receivers using FFT processing

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及一种外辐射源雷达距离徙动补偿的快速实现方法,属于雷达目标探测技术领域。该方法包括如下步骤:(1)将外辐射源雷达天线接收的直达波信号与回波信号分段处理,重排成二维矩阵;(2)将直达波和回波的每一行分别变换到频域;(3)将变换后的直达波矩阵与回波矩阵进行共轭相乘;(4)每一列经过低通滤波器之后下抽取;(5)对每一列分别进行CZT变换;(6)将每一行变换到时域,得到经过徙动补偿的互模糊函数结果。与传统的基于CZT+IFFT法实现keystone变换补偿距离徙动的方法相比,本发明在积累增益相当的情况下,减少了运算步骤,大大节省了运算量。

Description

一种外辐射源雷达距离徙动补偿的快速实现方法
技术领域
本发明涉及一种外辐射源雷达距离徙动补偿的快速实现方法,属于雷达目标探测技术领域。
背景技术
外辐射源雷达是一种特殊的双(多)基地无源雷达系统,其本身不发射信号,而是利用非合作辐射源(如广播、电视、通信信号等)对目标进行定位和跟踪。由于远距离目标的回波十分微弱,通常被噪声和对消后残余所淹没,因而需要采用长时间相参积累来提高积累增益。通常当积累时间较小时,随着积累时间的延长,积累增益不断提高,雷达的探测威力也随之提高。而当积累时间达到一定值后,目标在相参积累时间内位置和运动状态的变化使得互模糊函数中目标峰能量在时延维发生分散,产生了距离徙动现象,该现象如附图1所示。距离徙动的存在影响积累增益的提高,使信噪比下降,影响目标检测效果。所以在外辐射源雷达系统中需要研究如何对距离徙动进行补偿。
keystone变换是应用于合成孔径雷达领域中常用的一种距离徙动校正技术。合成孔径雷达的回波支撑域是个二维平面,平面坐标分别为快时间和慢时间。对快时间做傅里叶变换后,回波就变换到了快时间维频域-慢时间维时域,然后进行keystone变换,即对此二维数据平面进行慢时间尺度变换。下面以合成孔径雷达为例,分析keystone变换的原理。设点目标的基带回波信号表示为:
s r ( t ) = s r ( t ^ , t n ) = Ap ( t ^ , t n ) exp ( - j 4 π f c c R ( t n ) ) - - - ( 1 )
式中和tn=nTr分别表示快时间和慢时间,Tr为脉冲重复周期,A为目标回波幅度,R(tn)为tn时刻目标的距离,为归一化回波包络,fc为载波频率,c为光速。假定目标在tn时刻里近似匀速运动,则R(tn)=R0+vtn,v为目标径向速度。频域方法实现匹配滤波得到:
S r ( f , t n ) = A | P ( f ) | 2 exp ( - j 4 π c f c R 0 ) exp ( - j 4 π c fR 0 ) exp ( - j 4 π c fvt n ) exp ( - j 4 π c f c vt n ) - - - ( 2 )
式中P(f)为回波包络信号的傅里叶变换,-B/2≤f≤B2,B表示基带信号带宽。前两个指数项表示目标在t0时刻的距离R0对应的载波和各频率分量的相移;第三个指数项表示包络移位造成各频率分量的相移;第四个指数项为多普勒效应引起的载波相位变化,即多普勒频移fdc=-2v/λc,λc为载波波长。其中,第三个指数项是包络移位造成的,表现为不同频率分量的相位随慢时间tn的变化快慢不同,f越高,变化越快,即不同频率分量的多普勒频率不同,这是目标运动造成距离徙动在频域的体现。采用keystone变换对慢时间维进行尺度变换的方法,定义新的虚拟慢时间tn′:
t n ′ = f c + f f c t n - - - ( 3 )
进行变量替换可以得到:
S r ( f , t n ′ ) = A | P ( f ) | 2 exp [ - j 4 π c ( f c + f ) R 0 ] exp ( - j 4 π c f c vt n ′ ) - - - ( 4 )
以tn′为虚拟慢时间来度量,相位的变化率是恒定的,即不同脉冲信号之间的包络不再有平移,从而解决了距离徙动问题。keystone变换的原理如附图2所示。keystone变换前的数据格式如附图2(a)白色圆圈所示,变换后的数据格式如附图2(b)白色圆圈所示。由于变换后的数据需要进行快速傅里叶逆变换(简称IFFT)运算变换回时域,所以需要调整变换后的数据的采样间隔,使得不同频率对应的数据的采样间隔相等。调整后的数据如附图2(b)中黑色圆点所示。
对采样后的离散信号进行keystone变换,可以采用SINC插值法、DFT(离散傅里叶变换)+IFFT法、CZT(线性调频z变换)+IFFT法等方法实现。其中SINC插值法和DFT+IFFT法运算量很大,不利于实时处理,一般采用CZT+IFFT法实现keystone变换,其应用于外辐射源雷达信号处理的流程如附图3所示。keystone变换应用于外辐射源雷达时能很好的解决距离徙动的问题,但是运算量非常大,尽管CZT+IFFT法比SINC插值法和DFT+IFFT法运算量有一定程度的减小,但是算法运算量仍然很大,因此需要寻找快速算法,进一步减少运算量。
发明内容
本发明的目的针对外辐射源雷达系统检测微弱目标时,需要进行长时间相参积累而出现的距离徙动现象造成的积累增益降低的问题,提出一种外辐射源雷达距离徙动补偿的快速实现方法。
本发明的思想是首先通过采用低通滤波下抽的方法,节省了CZT变换的运算量,然后将CZT+IFFT法与多普勒滤波结合,通过CZT和快时间维IFFT得到互模糊结果,在不影响积累增益的前提下减少运算量,实现快速的距离徙动补偿。
本发明的目的是通过以下技术方案实现的:
一种外辐射源雷达距离徙动补偿的快速实现方法,包括以下步骤:
首先对算法中出现的符号进行说明:
B:基带信号带宽;
fs:外辐射源雷达系统中,接收到的信号经过A/D采样及数字下变频后的基带采样率;
fc:回波信号的载波频率;
Ts:相参积累时间;
P:低通滤波器的阶数;
N:分段重排后的二维矩阵慢时间维的采样点数,即分段数;
Nr:分段时每段重叠点数,即计算时延点数;
M:分段重排后的二维矩阵快时间维采样点数,即fsTs/N+Nr
D:下抽取倍数;
L:CZT变换的中间参数,L取满足L≥(2N/D)-1的2的整次幂的最小整数;
数据长度为N的快速傅里叶(简称FFT)变换所需复乘次数为:N(log2N)/2;
数据长度为N的CZT变换所需复乘次数为:3L(log2L)/2+6N+L;
步骤一,信号分段
根据系统需求得到雷达设计参数,确定信号分段数N和计算时延点数Nr,对回波信号和直达波信号进行分段处理。
对总长度为fsTs的一维直达波信号采用等长分段方式分段,把包含fsTs/N个采样点的每一段依次存为矩阵的一行构成N行fsTs/N列的二维矩阵,之后再在二维矩阵的每一行的末尾补上Nr个零值,构成N行fsTs/N+Nr列的二维矩阵,即N行M列的矩阵,记为r(n,m)(n=1,2...,N,m=1,2...,M)。对总长度为fsTs的一维回波信号采用重叠方式分段,构成N行M列的二维矩阵,记为e(n,m)(n=1,2...,N,m=1,2...,M),每一行数据包括一段之中等长部分fsTs/N个采样点和重叠部分Nr个采样点,每一段依次存为矩阵的一行,共N行。直达波和回波二维矩阵的行为快时间维,列为慢时间维。分段方式如附图4所示。
步骤二,快时间维FFT变换
将重构后的直达波和回波二维信号进行快时间维FFT,分别变换至快时间维频域,可以用下式表示:
R(n,k)=FFT[r(n,m)](5)
E(n,k)=FFT[e(n,m)](6)
其中r(n,m)表示直达波信号二维矩阵,e(n,m)表示回波信号二维矩阵,m表示快时间维时域每一行的采样点顺序,n表示慢时间维时域每一列的采样点顺序;R(n,k)表示直达波二维信号变换到快时间维频域的信号,E(n,k)表示回波二维信号变换到快时间维频域的信号,k表示快时间维频域每一行的采样点顺序。
步骤三,共轭相乘
将上一步得到的回波矩阵和直达波矩阵在快时间维频域共轭相乘,即每一行共轭相乘,实现频域的匹配滤波,公式如下式所示:
C(n,k)=R(n,k)·E*(n,k)(7)
C(n,k)表示直达波与回波矩阵每一行共轭相乘的结果二维矩阵。
步骤四,慢时间维低通滤波下抽取
为了降低运算量,在满足实际需求的参数设置情况下,通过对共轭相乘的结果低通滤波然后进行下抽取来降低采样率,从而降低了后续运算量。设低通滤波器系数的长度为P,系数为h(i),抽取倍数为D,则滤波后信号Cfir(n,k)可以表示为:
C fir ( n , k ) = Σ i = 0 P - 1 C ( n - i , k ) h ( i ) - - - ( 8 )
其中当n-i<0或n-i+1-P>0时,C(n-i,k)=0。抽取之后的信号CD(nD,k)可以表示为:
CD(nD,k)=Cfir(nDD,k)(9)
其中nD表示抽取之后慢时间维的采样点顺序,即列的采样点顺序,其满足nDD≤N。低通滤波下抽取并不会改变雷达的频率分辨率,并且选择合适的滤波器阶数、系数和抽取倍数,可以在不影响积累增益的前提下,大幅降低运算量。
步骤五,慢时间维CZT变换
常规通过CZT+IFFT法实现keystone变换需要两步,即先对慢时间维进行CZT变换,然后对慢时间维进行IFFT变换。常规方法中为了得到距离信息和多普勒信息,在完成keystone变换之后还需要对快时间维频域进行IFFT完成脉冲压缩,获得距离信息,接着对慢时间维进行FFT变换获取多普勒信息。本方法先进行慢时间维CZT变换,接着直接进行快时间维频域的IFFT变换,省略中间慢时间维的IFFT变换和完成脉冲压缩之后的慢时间维FFT。由于FFT和IFFT都是线性的,两步相“抵消”之后,就直接得到了多普勒信息,后续只需要对快时间维频域进行IFFT变换到时域就可以获得距离信息,从而得到经过徙动补偿的互模糊函数结果。这样减少了运算步骤和运算量。慢时间维CZT变换如下式所示:
Cczt(fd,k)=CZT[CD(nD,k)](10)
其中Cczt(fd,k)表示CZT变换之后的二维矩阵的结果,其中包含了目标的多普勒信息,fd表示慢时间维频域的采样点顺序。
步骤六,快时间维IFFT
对快时间维进行IFFT变换到时域,获取目标的距离信息,得到互模糊函数结果的二维数据矩阵C(fd,τ)。
C(fd,τ)=IFFT[Cczt(fd,k)](11)
其中τ表示Cczt(fd,k)经过IFFT变换之后快时间维时域的采样点顺序。
步骤七,目标检测
根据系统需求,选择合适的后续处理算法,对上一步得到的包含运动目标距离信息和多普勒信息的二维数据进行处理,得到目标的距离和速度信息。
有益效果
与传统使用CZT+IFFT法实现keystone变换方法进行距离徙动补偿方法相比,本发明提出的一种外辐射源雷达距离徙动补偿的快速实现方法计算量小,运算步骤简单,易于实现。算法流程步骤如附图5所示。
具体运算量为:CZT+IFFT法实现keystone变换方法总的复乘次数为:
3 2 NM log 2 M + 3 2 LM log 2 L + N 2 ( Nr + M ) log 2 N + ( 7 N + L ) M
本发明提出的方法复乘次数为:
5 4 NM log 2 M + 3 4 LM log 2 L 2 + 4 NM + NP + 1 2 LM
运算量对比分析如附图6所示。综上可得,本发明的运算复杂度大大降低,既节省了运算量,又易于实现。
附图说明
图1是背景技术中外辐射源雷达长时间相参积累的距离徙动现象示意图。
图2是背景技术中keystone变换原理示意图,其中(a)为变换前的数据示意图,(b)为变换后的数据示意图。
图3是背景技术中CZT+IFFT法实现keystone变换算法步骤示意图。
图4是本发明中的信号分段方式示意图,其中(a)为回波信号分段方法示意图,(b)为直达波信号分段方式示意图。
图5是本发明中的算法步骤流程示意图。
图6是本发明中的运算量对比示意图,其中(a)为当分段数固定,本发明快速计算法与CZT+IFFT法的运算量随着计算时延点数变化情况示意图;(b)为本发明快速计算法与CZT+IFFT法运算量的比值示意图。
图7是具体实施方式中互模糊函数侧视图,其中(a)为采用CZT+IFFT方法实现keystone变换的互模糊函数侧视图,(b)为本发明采用低通滤波下抽取结合CZT快速方法实现keystone变换的互模糊侧视图。
图8是本发明算法流程结构图。
具体实施方式
下面将结合具体实例及附图对本发明作进一步详细说明。
在外辐射雷达系统中,直达波天线接收直达波信号,回波天线接收目标反射的电视回波信号,使用检测民航飞机时采集的数据验证算法的有效性。
实验参数:中央电视广播塔发射信号载波频率为fc=674MHz,带宽B=7.56MHz,基带采样率fs=10MHz,相参积累时间Ts=1s,低通滤波器系数的长度P=21,其系数为:
h(l)=[0.0000,0.0036,-0.0000,-0.0122,0.0000,0.0343,-0.0000,
-0.0858,0.0000,0.3106,0.4991,0.3106,0.0000,-0.0858,-0.0000,
0.0343,0.0000,-0.0122,-0.000,0.0036,0.0000]
慢时间维采样点数即分段数N=2500,每段不重叠的采样点数fsTs/N=4000,分段时重叠的点数Nr=4000,下抽取倍数D=2,L=4096,检测目标为民航飞机,其多普勒频率fd=499Hz。
如图8所示,本发明算法具体流程如下:
1)信号分段
直达波与回波的分段数均为N=2500,即行数为2500;回波重叠分段每个脉冲的不重叠的数据点数4000,重叠点数为Nr=4000,故回波分段二维重排之后每一行的数据点数为M=8000。直达波每段补零个数为4000,故直达波分段二维重排之后每一行点数为M=8000,行为快时间维,列为慢时间维。
2)快时间维FFT变换
分别依次对直达波和回波的每一行做8000点的FFT变换。
3)共轭相乘
把做完FFT变换的直达波2500行8000列二维矩阵与做完FFT变换的回波2500行8000列二维矩阵的对应值共轭相乘,得到2500行8000列的二维矩阵C(n,k)。
4)慢时间维低通滤波下抽取
把得到的C(n,k)二维矩阵的慢时间维即每一列与阶数P=21的低通滤波器系数卷积,之后对每一列进行D=2的2倍下抽取,得到1250行8000列的二维矩阵CD(nD,k)。
5)慢时间维CZT变换
对CD(nD,k)的每一列进行N/D即1250点的CZT变换,其中CZT变换中的参数L=4096,共8000列,得到1250行8000列的矩阵Cczt(fd,k)。
6)快时间维IFFT变换
对矩阵Cczt(fd,k)的每一行做8000点的IFFT,共1250行,得到包含距离和多普勒信息的1250行8000列的互模糊函数数据矩阵C(fd,τ),其侧视图如附图7(b)所示。
7)目标检测
对上一步得到的包含运动目标距离信息和多普勒信息的结果矩阵C(fd,τ),根据需要进行后续处理,如采用恒虚警算法,判断目标有无,最后得到目标的距离和速度。
采用CZT+IFFT法实现keystone变换总的复乘次数为2.0418×109,多普勒单元1250对应多普勒频率为0Hz,目标处于1749多普勒单元,因此目标的多普勒频率fd=499Hz,其信噪比为26.59dB。采用低通滤波下抽取结合CZT快速算法实现keystone变换总的复乘次数为1.0268×109,多普勒单元625对应多普勒频率为0Hz,目标处于1124多普勒单元,则目标的多普勒频率fd=499Hz,其信噪比为26.59dB。由此可见本发明采用的方法在不影响积累增益的前提下节省了运算量,减少了运算步骤,降低了算法复杂度。
以上所述仅为本发明的具体实施例,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (1)

1.一种外辐射源雷达距离徙动补偿的快速实现方法,包括以下步骤:
步骤一,信号分段
根据系统需求得到雷达设计参数,确定信号分段数N和计算时延点数Nr,对回波信号和直达波信号进行分段处理;
对总长度为fsTs的一维直达波信号采用等长分段方式分段,把包含fsTs/N个采样点的每一段依次存为矩阵的一行构成N行fsTs/N列的二维矩阵,之后再在二维矩阵的每一行的末尾补上Nr个零值,构成N行fsTs/N+Nr列的二维矩阵,即N行M列的矩阵,记为r(n,m)(n=1,2...,N,m=1,2...,M);其中fs表示外辐射源雷达系统中,接收到的信号经过A/D采样及数字下变频后的基带采样率;Ts表示相参积累时间;
对总长度为fsTs的一维回波信号采用重叠方式分段,构成N行M列的二维矩阵,记为e(n,m),每一行数据包括一段之中等长部分fsTs/N个采样点和重叠部分Nr个采样点,每一段依次存为矩阵的一行,共N行;
记直达波和回波二维矩阵的行为快时间维,列为慢时间维;
步骤二,快时间维FFT
将重构后的直达波和回波二维信号进行快时间维FFT,分别变换至快时间维频域,可以用下式表示:
R(n,k)=FFT[r(n,m)];
E(n,k)=FFT[e(n,m)];
其中r(n,m)表示直达波信号二维矩阵,e(n,m)表示回波信号二维矩阵,m表示快时间维时域每一行的采样点顺序,n表示慢时间维时域每一列的采样点顺序;
R(n,k)表示直达波二维信号变换到快时间维频域的信号,E(n,k)表示回波二维信号变换到快时间维频域的信号,k表示快时间维频域每一行的采样点顺序;
步骤三,共轭相乘
将上一步得到的R(n,k)和E(n,k)在快时间维频域共轭相乘,即每一行共轭相乘,实现频域的匹配滤波,公式如下式所示:
C(n,k)=R(n,k)·E*(n,k);
C(n,k)表示R(n,k)与E(n,k)每一行共轭相乘的结果二维矩阵;
步骤四,慢时间维低通滤波下抽取
设低通滤波器系数的长度为P,系数为h(i),抽取倍数为D,则滤波后信号Cfir(n,k)可以表示为:
C f i r ( n , k ) = Σ i = 0 P - 1 C ( n - i , k ) h ( i ) ;
其中当n-i<0或n-i+1-P>0时,C(n-i,k)=0;抽取之后的信号CD(nD,k)可以表示为:
CD(nD,k)=Cfir(nDD,k);
其中nD表示抽取之后慢时间维的采样点顺序,即列的采样点顺序,其满足
nDD≤N;
步骤五,慢时间维CZT变换
慢时间维CZT变换如下式所示:
Cczt(fd,k)=CZT[CD(nD,k)];
其中Cczt(fd,k)表示CD(nD,k)经过CZT变换之后的二维矩阵,其中包含了目标的多普勒信息,fd表示慢时间维频域的采样点顺序;
步骤六,快时间维IFFT
快时间维IFFT变换如下式所示:获取目标的距离信息,得到C(fd,τ)。
C(fd,τ)=IFFT[Cczt(fd,k)]
其中C(fd,τ)表示Cczt(fd,k)经过IFFT变换之后的互模糊函数结果二维矩阵,τ表示IFFT变换之后快时间维时域的采样点顺序;
步骤七,目标检测
根据系统需求,选择合适的后续处理算法,对上一步得到的包含运动目标距离信息和多普勒信息的二维数据进行处理,得到目标的距离和速度信息。
CN201410306628.4A 2014-06-30 2014-06-30 一种外辐射源雷达距离徙动补偿的快速实现方法 Expired - Fee Related CN104062640B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410306628.4A CN104062640B (zh) 2014-06-30 2014-06-30 一种外辐射源雷达距离徙动补偿的快速实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410306628.4A CN104062640B (zh) 2014-06-30 2014-06-30 一种外辐射源雷达距离徙动补偿的快速实现方法

Publications (2)

Publication Number Publication Date
CN104062640A CN104062640A (zh) 2014-09-24
CN104062640B true CN104062640B (zh) 2016-03-30

Family

ID=51550441

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410306628.4A Expired - Fee Related CN104062640B (zh) 2014-06-30 2014-06-30 一种外辐射源雷达距离徙动补偿的快速实现方法

Country Status (1)

Country Link
CN (1) CN104062640B (zh)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104793186A (zh) * 2015-01-21 2015-07-22 佛山市智海星空科技有限公司 基于多抽样率和速度通道划分的徒动补偿方法及系统
CN104714213A (zh) * 2015-03-30 2015-06-17 宁波高新区宁源科技服务有限公司 一种基于多抽样率和速度通道划分的徒动补偿方法
CN105388467B (zh) * 2015-11-03 2017-08-08 南京信息工程大学 一种修正多普勒天气雷达回波衰减的方法
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 电子科技大学 一种基于信号融合的外辐射源微弱目标定位的方法
CN108549064A (zh) * 2018-07-24 2018-09-18 电子科技大学 基于脉内多普勒频率模糊补偿的外辐射源动目标探测方法
CN111580057A (zh) * 2019-02-19 2020-08-25 北京华航无线电测量研究所 一种运动目标回波距离徙动校正方法
CN111650563B (zh) * 2020-06-15 2022-05-31 桂林电子科技大学 一种外辐射源雷达同频干扰时延和能量快速估计系统及方法
CN112130140B (zh) * 2020-08-24 2023-12-19 南京理工大学 一种伪码调相引信系统提高测速测距精度的方法
US11668790B2 (en) 2021-05-25 2023-06-06 Nxp B.V. Radar communications with oversampling
CN113884999A (zh) * 2021-09-27 2022-01-04 中国人民解放军32802部队 一种外辐射源雷达的跨距离单元补偿方法
US20230194705A1 (en) * 2021-12-20 2023-06-22 Aptiv Technologies Limited Movement Compensation for Multi-Dimensional Radar Signal Processing
CN114384484B (zh) * 2022-01-24 2023-01-24 电子科技大学 一种基于分段处理的匀加速运动目标快速相参积累方法
CN117111109A (zh) * 2023-08-28 2023-11-24 南京威翔科技有限公司 一种用于低空监视的时序控制方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2677343A1 (en) * 2011-02-18 2013-12-25 Mitsubishi Electric Corporation Passive radar device

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5865794B2 (ja) * 2012-07-13 2016-02-17 三菱電機株式会社 レーダ装置
WO2014092052A1 (ja) * 2012-12-14 2014-06-19 三菱電機株式会社 レーダ装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2677343A1 (en) * 2011-02-18 2013-12-25 Mitsubishi Electric Corporation Passive radar device

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种高效的外辐射源雷达高径向速度目标实时检测方法;关欣等;《电子与信息学报》;20130331;第35卷(第3期);581-588 *
数字电视辐射源雷达的相参积累徙动补偿方法;杨金禄等;《电子与信息学报》;20110228;第33卷(第2期);407-411 *

Also Published As

Publication number Publication date
CN104062640A (zh) 2014-09-24

Similar Documents

Publication Publication Date Title
CN104062640B (zh) 一种外辐射源雷达距离徙动补偿的快速实现方法
CN103616679B (zh) 基于差波束调制和波形分析的pd雷达测距测角方法
CN103675759B (zh) 一种改进的分数阶傅里叶变换机动弱目标检测方法
CN102288941B (zh) 一种基于fpga和dsp的中频lfm-pd雷达信号实时处理系统及处理方法
CN104698459B (zh) 一种应用于缺失数据的条带sar压缩感知成像方法
CN103760558B (zh) 一种太赫兹雷达isar成像方法
CN102879783B (zh) 基于稀疏探测频率信号的isar成像方法
CN102608602B (zh) 一种基于完全互补序列的超低旁瓣合成孔径雷达成像方法
CN103454632B (zh) 一站固定式调频连续波双基地sar成像方法
CN103901419B (zh) 一种基于频域相位校正的外辐射源雷达距离徙动补偿方法
Li et al. ISAR imaging of nonuniformly rotating target based on the multicomponent CPS model under low SNR environment
CN104502898B (zh) 将修正rft和修正mdcft相结合的机动目标参数估计方法
CN103777178B (zh) 一种同步误差补偿方法、设备及系统
CN104215948A (zh) 一种基于参考信号的星载sar回波校正和脉冲压缩方法
CN106507958B (zh) 外辐射源雷达信号实时相干积累的方法
Tian et al. High-speed maneuvering target detection approach based on joint RFT and keystone transform
CN104076343A (zh) 星载三通道sar-gmti自适应杂波抑制方法
CN102608603B (zh) 一种基于完全互补序列的多通道合成孔径雷达成像方法
CN102445691A (zh) 一种多通道星载合成孔径雷达方位频谱稀疏重建方法
CN103293528B (zh) 扫描雷达超分辨成像方法
CN103323818A (zh) 多通道合成孔径雷达系统非均匀采样奇异点的方法和装置
CN103630878A (zh) 一种基于迭代自适应方法的频带融合方法及装置
CN104849708A (zh) 基于频域多项式相位变换的高速机动目标参数估计方法
CN104345300A (zh) 杂波空时谱线性补偿的机载非正侧视阵雷达stap方法
CN109031299B (zh) 低信噪比条件下基于相位差分的isar平动补偿方法

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: 20170630

CF01 Termination of patent right due to non-payment of annual fee