CN103529437A - 系留气球载相控阵雷达在多目标下分辨空地目标的方法 - Google Patents
系留气球载相控阵雷达在多目标下分辨空地目标的方法 Download PDFInfo
- Publication number
- CN103529437A CN103529437A CN201310482277.8A CN201310482277A CN103529437A CN 103529437 A CN103529437 A CN 103529437A CN 201310482277 A CN201310482277 A CN 201310482277A CN 103529437 A CN103529437 A CN 103529437A
- Authority
- CN
- China
- Prior art keywords
- angle
- target
- clutter
- depression
- array
- 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
- 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
- G01S7/414—Discriminating targets with respect to background clutter
-
- 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
- G01S13/00—Systems 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/52—Discriminating between fixed and moving objects or between objects moving at different speeds
-
- 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/021—Auxiliary means for detecting or identifying radar signals or the like, e.g. radar jamming signals
-
- 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
- G01S7/415—Identification of targets based on measurements of movement associated with the target
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、将各目标的俯仰角与地杂波下俯角范围进行比较,得到各目标的分辨结果。本发明能精确估计各目标的俯仰角,并利用目标俯仰信息实现空地多目标分辨,可用于提高系留气球载相控阵雷达对低空目标的侦察能力。
Description
技术领域
本发明涉及雷达技术领域,具体说是一种多目标情况下系留气球载相控阵雷达分辨空地目标的方法,可用于预警探测、通信中继、电子干扰、防灾减灾和公共安全。
背景技术
系留气球工作时悬浮在数千米的高空,可作为浮空平台搭载预警雷达系统,组成系留气球载雷达,进行下视探测。其克服了地基雷达易受地球曲率和地物遮蔽等影响的固有缺点,大大提高了侦察探测距离,同时,相比于预警飞机,可获得更高的效费比,是现今侦查巡航导弹、低空突防飞行器等低空低速目标的重要手段之一。
由于系留气球载相控阵雷达通常需要对检测到的低空目标和地面目标进行不同的后续处理,因此,系留气球载相控阵雷达首先要分辨出多个目标,再判断出多个目标中,哪些属于空中目标,哪些属于地面目标。常规雷达系统中,对多个目标进行分辨主要依靠的处理方式有:距离维处理、多普勒维处理、角度维处理等。利用上述处理方式实现多目标分辨,首先要求多个目标的信息能在某一维或某几维被分离,再根据分离后的各个目标的信息,利用单目标判别方法分别对每个目标进行判别,以实现多目标分辨。
但是,由于系留气球载相控阵雷达工作在下视方式,车辆等地面运动目标同样会被雷达探测到。对于距离维分辨,由于系留气球载相控阵雷达的俯仰维距离环经过地面,若地面目标与低空目标处于同一个俯仰维距离环内,即与雷达的径向距离相同,将无法通过距离维的信号处理分离各个目标的信息,因此,无法依靠距离维处理实现空地目标分辨。
同时,低空目标距离地面高度一般在几百米左右,在非向站飞行时,低空目标相对系留气球载相控阵雷达的径向速度与地面运动目标相对系留气球载相控阵雷达的径向速度通常很接近,空地目标的多普勒频率可能会落入同一个多普勒通道中,即无法通过多普勒维的信号处理分离各个目标的回波信息。此时,在上述距离维处理无法分辨空地目标的情况下,同样也无法依靠多普勒维处理实现空地目标的分辨。
在上述距离维处理和多普勒维处理均不能实现空地目标分辨的情况下,考虑到同一个俯仰维距离环中,低空目标与地面运动目标的俯仰角不同,而地面运动目标的俯仰角与同一距离单元中的地面下俯角相同,故可尝试在角度维分离多个目标的信息,将问题转化为单目标判别问题。但实际情况中,各目标俯仰角差异较小,现有系留气球载相控阵雷达所使用的大多是常规的相控阵角度估计方法,此方法不具备足够的分辨率,无法精确估计各目标的俯仰角,即不能很好的在角度维分离各目标的信息,因此,现有系留气球载相控阵雷达无法利用各目标的俯仰角信息实现角度维的多目标分辨。
发明内容
本发明的目的在于针对上述现有技术的不足,提出一种系留气球载相控阵雷达在多目标下分辨空地目标的方法,以精确估计各目标的俯仰角,实现对空地多目标的分辨。
本发明的技术方案是这样实现的:
一、技术原理
相控阵雷达在俯仰维拥有较多的自由度,在对各个接收通道积累多次相关的脉冲重复周期数据后,首先进行距离维和多普勒维的处理,再通过距离和多普勒二维检测,分离目标与地杂波信息,对地杂波下俯角的范围和各目标进行角度估计,最后将每个目标的俯仰角与地杂波下俯角的范围进行比较,得出判别结果,实现空地多目标的分辨。本发明中,由于对地面下俯角的范围进行估计时,不存在分辨的要求,故应用常规的相控阵角度估计方法对地面下俯角进行估计,而在估计各目标的俯仰角时,则使用超分辨相控阵角度估计方法,以精确估计各目标的俯仰角。
二、技术方案
根据上述原理,本发明的技术方案包括如下步骤:
(1)对相控阵雷达俯仰维各个阵元接收到的回波信号进行脉冲压缩,得到每一个脉冲重复周期的脉冲压缩结果Xm(n)=[x1m(n),x2m(n),...,xim(n),...,xLm(n)]T,其中,m表示脉冲重复周期序号,m=1,2,...,M,M为一个相干积累周期中所包含的脉冲个数,n为距离单元的序号,n=1,2,...,N,N为一个脉冲重复周期中的距离单元数,xim(n)表示第m次脉冲重复周期中第i个阵元的脉冲压缩结果,i为阵元编号,i=1,2,...,L,L为等距均匀线阵的阵元个数,[·]T表示转置;
(2)对每一个脉冲重复周期中的脉冲压缩结果Xm(n)在俯仰维做多约束方向的数字波束形成,得到每个约束方向的M路合成信号其中,ψe表示一个俯仰维约束方向,e=1,2,...,E,E表示俯仰维约束方向的个数;
(4)对每个约束方向的R-D图,沿距离维和多普勒维进行二维恒虚警率检测,记录出现目标的R-D图所对应的波束指向ψt和目标所在距离单元nt;
(5)在每一个脉冲重复周期的脉冲压缩结果Xm(n)中,提取目标所在距离单元nt的脉冲压缩结果Xm(nt),记为快拍zm=Xm(nt),将所有快拍组成快拍矩阵Z=[z1z2...zm...zM],其中,m=1,2,...,M;
(6)将快拍矩阵Z中的M列数据直接相加,得到L个阵元在多普勒维的零通道数据,用此L个数据组成用于估计地杂波下俯角的一组数据向量V=[v1,v2,...,vi,...,vL]T,其中vi表示第i个阵元的多普勒零通道数据,i=1,2,...,L;
(7)根据距离单元nt和雷达距离分辨率,得到以距离单元nt为中心的分辨单元,记为nr,由目标所在的距离单元nt与雷达距离分辨率,计算分辨单元nr中的远端下俯角θf与近端下俯角θc;
(8)根据远端下俯角θf与近端下俯角θc,对数据向量V使用阵列角度估计方法,得到距离单元nt中地杂波下俯角的范围:θmin~θmax,其中,θmin=θg-ε,θmax=θg+ε,θg是由阵列角度估计方法直接估计得到的地杂波下俯角,ε为角度偏差量,其取值与该距离单元的杂噪比以及实际中关于地杂波空间分布的先验知识有关;
(9)对步骤5中的快拍矩阵Z做三脉冲对消杂波抑制,得到不含杂波的快拍矩阵U,其中,不含杂波的快拍矩阵U是L行T列的矩阵,L是阵元数,T表示三脉冲对消杂波抑制后得到的脉冲个数,T=M-2;
(10)根据不含杂波的快拍矩阵U,得到目标的协方差矩阵W,应用信源数估计方法,计算距离单元nt中所包含的空地目标个数Q;
(11)根据空地目标个数Q,估计每个目标的俯仰角θq,其中,q=1,2,...,Q;
(12)将每个目标的俯仰角θq分别与步骤8得到的地杂波下俯角范围θmin~θmax进行比较,若θmin≤θq≤θmax,则判定该目标为地面目标,否则,判定该目标为空中目标。
本发明与现有技术相比具有以下优点:
1、本发明由于在系留气球载相控阵雷达上使用阵列超分辨角度估计算法,使系留气球载相控阵雷达具备了对多个目标的俯仰角进行精确估计的能力。
2、本发明由于利用俯仰信息分辨空地目标,有效解决了现有系留气球载相控阵雷达无法在多目标情况下实现空地目标分辨的问题。
3、本发明由于使用零多谱勒通道的数据估计地杂波的下俯角,充分利用了零多普勒通道中地杂波的高杂噪比,因此可以较准确的估计出地杂波的下俯角。
4、本发明由于将实际地杂波下俯角作为对应距离单元的地面下俯角,有效地克服了采用三角函数关系计算地面下俯角时,结果受地形起伏影响较大的缺点。
附图说明
图1是本发明使用的场景布局图;
图2是本发明的实现流程图;
图3是本发明约束多个波束的空域模拟示意图;
图4是本发明中分辨单元与距离单元关系示意图;
图5是用本发明方法和常规阵列角度估计方法对角度的估计结果比较图;
图6是本发明方法实现多目标分辨的统计性能仿真图;
图7是本发明中空中目标被正确分辨的概率随空中目标下俯角变化的仿真结果图。
具体实施方式
参照图1,本发明使用的场景包括:一个系留气球载相控阵雷达、一个低空目标、一个地面目标和目标所在距离单元中的地杂波,主要场景参数有:系留气球载相控阵雷达的高度H,雷达阵元数L,目标距离雷达的径向距离R;雷达工作参数有:雷达发射信号的带宽B,雷达采样频率fs。
参照图2,本发明的实施步骤如下:
步骤1、对相控阵雷达俯仰维各阵元接收信号进行脉冲压缩。
对相控阵雷达俯仰维各个阵元接收到的回波信号进行脉冲压缩,得到每一个脉冲重复周期的脉冲压缩结果Xm(n)=[x1m(n),x2m(n),...,xim(n),...,xLm(n)]T,其中,m表示脉冲重复周期序号,m=1,2,...,M,M为一个相干积累周期中所包含的脉冲个数,n为距离单元的序号,n=1,2,...,N,N为一个脉冲重复周期中的距离单元数,xim(n)表示第m次脉冲重复周期中第i个阵元的脉冲压缩结果,i为阵元编号,i=1,2,...,L,L为等距均匀线阵的阵元个数,[·]T表示转置。
步骤2、对每一个脉冲重复周期中的脉冲压缩结果Xm(n)在俯仰维做多约束方向的数字波束形成。
2a)设系留气球载雷达在俯仰维的探测范围为ψmin~ψmax,在该探测范围内以半功率波束宽度ψ0.5为间隔,均匀选择多个俯仰方向ψe,其对应的空域导向矢量a(ψe):
其中,m=1,2,...,M,M为一个相干积累周期中所包含的脉冲个数,[·]H表示共轭转置,根据阵列信号处理可知,每一个空域导向矢量a(ψe)均对应一个天线方向图,该方向图的主瓣指向由ψe决定,经过空域导向矢量a(ψe)加权求和的合成信号代表了约束方向ψe的回波信息,图3给出了指向各个约束方向ψe的方向图的主瓣。
实际应用中有多种杂波抑制方法可供选择,如脉冲对消杂波抑制方法,自适应杂波抑制方法等。本发明使用脉冲对消杂波抑制方法中常见的三脉冲对消杂波抑制方法为例说明处理过程,在条件允许情况下,也可以使用其他的杂波抑制方法。
本步骤的具体实现如下:
3a)记三脉冲对消杂波抑制所用的滤波器F=[h0,h1,h2]T,滤波器F的系数分别为:h0=1,h1=-2,h2=1,其中,[·]T表示转置;
其中,p=1,2,...,T,T表示三脉冲对消后得到的脉冲个数,T=M-2,M为一个相干积累周期中所包含的脉冲个数;
步骤4、对每个约束方向的R-D图进行二维恒虚警率检测,记录出现目标的R-D图所对应的波束指向ψt和目标所在距离单元nt。
4a)对每个约束方向的R-D图,按照奈曼-皮尔逊准则设定检测门限ξ;
4b)将各R-D图中每个距离单元的幅度值与门限ξ进行比较,结合雷达距离分辨率,在幅度值超过门限的距离单元中确定目标,记录出现目标的波束指向ψt和目标所在距离单元nt。
步骤5、根据目标所在距离单元nt构造快拍矩阵Z。
在每一个脉冲重复周期的脉冲压缩结果Xm(n)中,提取目标所在距离单元nt的脉冲压缩结果Xm(nt),记为快拍zm=Xm(nt),将所有快拍组成快拍矩阵Z=[z1z2...zm...zM],其中,m=1,2,...,M。
步骤6、将快拍矩阵Z中的M列数据直接相加,得到L个阵元在多普勒维的零通道数据,用此L个数据组成用于估计地杂波下俯角的一组数据向量V=[v1,v2,...,vi,...,vL]T,其中vi表示第i个阵元的多普勒零通道数据,i=1,2,...,L。
步骤7、计算分辨单元nr中的远端下俯角θf与近端下俯角θc。
根据距离单元nt和雷达距离分辨率,得到以距离单元nt为中心的分辨单元,记为nr,距离单元nt与分辨单元nr的关系如图4所示。将该分辨单元nr中远离雷达的端点与阵列天线相位中心的连线相对雷达法线方向的夹角,记为远端下俯角θf,该分辨单元nr中远离雷达的端点与阵列天线相位中心的连线相对雷达法线方向的夹角,记为远端下俯角θc,其计算步骤如下:
7a)根据目标所在距离单元nt,计算该距离单元的径向距离R:
R=ntc/2fs,
其中,c为光速,fs为雷达采样频率;
7b)根据径向距离R,计算远端下俯角θf:
其中,H为雷达距地面高度,ΔR表示雷达的距离分辨率,ΔR=c/2B,B为雷达发射信号的带宽;
7c)根据径向距离R,计算近端下俯角θc:
步骤8、根据远端下俯角θf、近端下俯角θc和步骤6中的数据向量V,使用阵列角度估计方法计算地杂波下俯角的范围θmin~θmax。
场景参数中系留气球高度和目标距离雷达的径向距离为三角关系,通过该三角关系可计算出目标所在距离单元的地面下俯角,用于对后续目标的判定,但由于该计算方法未考虑实际中地形的起伏特性与杂波的空间分布特性,会使得计算得到的地面下俯角存在较大误差,故本发明选择从杂波的回波信息中估计地杂波的俯仰角,并结合实际中能够得到的杂波空间分布特性,得到地杂波下俯角的范围,作为地面下俯角的范围,避免了上述误差的产生。其具体实现步骤如下:
8a)由远端下俯角θf和近端下俯角θc,构造地杂波下俯角的搜索范围θf~θc;
其中,d为系留气球载相控阵雷达的阵元间距,λ为雷达发射信号的波长,j为虚数符号, [·]T表示转置;
其中,[·]H表示共轭转置;
8e)根据地杂波下俯角θg,得到地杂波下俯角的范围:θmin~θmax,其中,θmin=θg-ε,θmax=θg+ε,ε为角度偏差量,它的取值与该距离单元的杂噪比以及实际中关于地杂波空间分布的先验知识有关。
步骤9、对步骤5中的快拍矩阵Z做三脉冲对消杂波抑制,得到不含杂波的快拍矩阵U。
本步骤中所述的三脉冲对消杂波抑制,其处理方式与步骤3中所述的三脉冲对消杂波抑制相同。
步骤10、根据不含杂波的快拍矩阵U,得到目标的协方差矩阵W,应用信源数估计方法,计算距离单元nt中所包含的空地目标个数Q。
常见信源数估计方法有很多,如最小描述长度方法,信息论方法,盖尔圆法等,本发明使用信息论方法为例说明处理过程,在条件允许情况下,也可以使用其他的信源数估计方法。使用信源数估计方法有两个前提条件:一是要求各个目标是不相干的;二是要求目标个数小于阵元数。这两个条件在实际中一般容易满足,故本发明假设多个目标之间均不相干,同时,目标个数最多不超过阵元数。本步骤的具体实现方式如下:
10a)根据不含杂波的快拍矩阵U,得到目标的协方差矩阵W:
其中,协方差矩阵W为L维方阵,L是阵元数,T表示三脉冲对消杂波抑制后得到的脉冲个数,T=M-2,[·]H表示共轭转置;
10b)对协方差矩阵W进行特征分解,得到L个特征值,记为ατ,其中,τ=1,2,...,L,L是阵元数;
10c)从1~L的所有整数中搜索目标的个数βσ,其中,σ=1,2,...,L;
其中,似然函数
步骤11、根据空地目标个数Q,估计每个目标的俯仰角θq,其中,q=1,2,...,Q。
现有的阵列超分辨角度估计方法有:线性约束最小方差方法、多重信号分类方法、旋转因子不变法、最大似然方法等,本发明以最大似然方法为例说明处理过程。当然,在条件允许的情况下,也可以采用其它的超分辨测角方法。
最大似然方法是一种基于多维搜索的角度估计方法,每个目标均具有一个角度搜索范围,当没有目标角度的先验知识时,一般令所有目标具有相同的角度搜索范围。
本步骤的具体实现如下:
11a)根据出现目标的波束指向ψt和雷达半功率波束宽度ψ0.5构造最大似然方法的目标角度搜索范围:φmin~φmax,并从此范围中均匀选出D个角度,其中,φmin=ψt-ψ0.5/2,φmax=ψt+ψ0.5/2;
11b)记各目标的角度搜索变量为每个目标均从上述的D个角度中任意选取一个角度,共得到DQ个不同的角度组合其中,kq=1,2,...,D,q=1,2,...,Q,Q为目标个数,Υ=1,2,...,DQ,DQ表示D的Q次方;
其中,d为系留气球载相控阵雷达的阵元间距,λ为雷达发射信号的波长,j为虚数符号, 表示转置;
11d)根据目标的协方差矩阵W,对每一个俯仰角组合ΓΥ,计算其代价值J(ΓΥ):
J(ΓΥ)=tr(PAW),
11e)记所有俯仰角组合的代价值J(ΓΥ)中的最大值为Jmax,Jmax对应的俯仰角组合中的Q个角度,即为各目标的俯仰角θq,其中,q=1,2,...,Q。
步骤12、将每个目标的俯仰角θq分别与步骤8得到的地杂波下俯角范围θmin~θmax进行比较,若θmin≤θq≤θmax,则判定该目标为地面目标,否则,判定该目标为空中目标。
本发明的可行性和优点可通过以下仿真进一步说明。
一、实验场景:
1.采用图1所示的场景布局图,整个场景由一个系留气球载相控阵雷达、一个低空目标、一个地面目标和目标所在距离单元中的地杂波组成,主要仿真参数有:系留气球载相控阵雷达的高度H=3km,雷达阵元数L=30,两目标距离雷达的水平距离S=20km,低空目标距地面高度H1=500m,雷达发射信号带宽B=1MHz,雷达采样频率fs=2MHz,雷达距离分辨率ΔR=300m,地杂波下俯角的角度偏差量ε=0.1°,低空目标的俯仰角θ1=7.7°,地面目标的俯仰角θ2=8.53°,地杂波下俯角θg=8.53°,不含杂波的脉冲数T=60,两个目标的信噪比SNR均为15dB,杂噪比CNR=40dB,噪声为加性高斯白噪声。
二、实验内容:
实验一:为证明本发明能够精确估计各个目标的俯仰角,并实现多目标分辨,用本发明方法和常规阵列角度估计方法对上述参数进行仿真,结果如5所示。
从图5中还可以看到,由于空地目标的俯仰角差异很小,常规的阵列角度估计方法无法对各个角度进行精确估计,而本发明方法对各目标的俯仰角实现了较准确的估计。同时,由于判定此目标为空中目标,判定此目标为地面目标,此判别结果与场景设置一致,证明了本发明能够精确估计各个目标的俯仰角信息,并利用此信息实现空地多目标分辨。
实验二:为说明本发明方法在进行空地多目标分辨时的统计性能,令空中目标下俯角θ1=8.51°,其他参数不变,用本发明方法对场景中的空中目标和地杂波进行100次仿真,结合最大似然角度估计方法的估计结果服从高斯分布这一原理,得到空中目标俯仰角和地杂波下俯角的估计结果分布曲线,其仿真结果如图6所示,图6中虚线表示的曲线是空中目标俯仰角的的估计结果分布曲线,实线表示的曲线是地杂波下俯角的估计结果分布曲线,两根竖直实线之间的角度分别表示和同时,在和区间内,对空中目标的高斯分布概率密度函数进行积分,即计算图中阴影部分的面积,得到当前参数下,空中目标被正确分辨的概率为0.91。
本实验中,由于空中目标和地面目标的信噪比相同,空地目标俯仰角的估计结果具有相似的分布曲线,差别仅是曲线中心对应的俯仰角不同,考虑到系留气球载相控阵雷达通常对空中目标更感兴趣,故仅对空中目标的估计结果分布曲线及正确分辨的概率进行了仿真,对于地面目标,可以由类似的仿真过程得到对应的结果。
由图6可以看出,由于高斯白噪声的存在,本发明方法的角度估计结果呈高斯分布,各次的估计结果与真实值存在一定的偏差,但本发明方法仍然能在较高的概率下实现空中目标的分辨。
实验三:为进一步说明本发明方法进行空地多目标分辨的性能,固定地面目标的俯仰角不变,改变场景中空中目标的俯仰角,设其俯仰角变化范围为:8.05°~8.52°,变化间隔为0.01°,每变化一次,利用实验二的方法计算空中目标被正确分辨的概率,得到空中目标被正确分辨的概率随其俯仰角变化的曲线,其仿真结果如图7所示,图7中的虚线表示地杂波下俯角的估计值。由于系留气球载相控阵雷达通常对空中目标更感兴趣,故本实验仅对空中目标的变化进行了仿真,对于地面目标,可以由类似的仿真过程得到对应的结果。
从图7中可以看到,在俯仰角的变化范围8.05°~8.50°内,空中目标被正确分辨的概率保持为1,随着空中目标的俯仰角进一步增大,空中目标的俯仰角与地杂波下俯角的差逐渐减小,空中目标被正确分辨的概率开始下降。根据实验二可知,这是因为随着空中目标俯仰角的与地杂波下俯角逐渐接近,对空中目标俯仰角的估计值有更大的概率落入地杂波下俯角的范围中,导致正确分辨的概率下降。记正确分辨概率为0.9时对应的空中目标俯仰角为临界角θr,从图7中可知,在本实验参数下,临界角θr约为8.51°。
在实际应用中,可以根据雷达工作参数和工作环境确定某以距离单元的临界角θr,本发明方法在空中目标的俯仰角小于临界角θr时能较好的实现空地多目标分辨。
Claims (5)
1.一种系留气球载相控阵雷达在多目标下分辨空地目标的方法,包括如下步骤:
(1)对相控阵雷达俯仰维各个阵元接收到的回波信号进行脉冲压缩,得到每一个脉冲重复周期的脉冲压缩结果Xm(n)=[x1m(n),x2m(n),...,xim(n),...,xLm(n)]T,其中,m表示脉冲重复周期序号,m=1,2,...,M,M为一个相干积累周期中所包含的脉冲个数,n为距离单元的序号,n=1,2,...,N,N为一个脉冲重复周期中的距离单元数,xim(n)表示第m次脉冲重复周期中第i个阵元的脉冲压缩结果,i为阵元编号,i=1,2,...,L,L为等距均匀线阵的阵元个数,[·]T表示转置;
(2)对每一个脉冲重复周期中的脉冲压缩结果Xm(n)在俯仰维做多约束方向的数字波束形成,得到每个约束方向的M路合成信号其中,ψe表示一个俯仰维约束方向,e=1,2,...,E,E表示俯仰维约束方向的个数;
(3)对每个约束方向的M路合成信号进行三脉冲对消杂波抑制及相干积累,得到每个约束方向的距离-多普勒二维图,即R-D图;
(4)对每个约束方向的R-D图,沿距离维和多普勒维进行二维恒虚警率检测,记录出现目标的R-D图所对应的波束指向ψt和目标所在距离单元nt;
(5)在每一个脉冲重复周期的脉冲压缩结果Xm(n)中,提取目标所在距离单元nt的脉冲压缩结果Xm(nt),记为快拍zm=Xm(nt),将所有快拍组成快拍矩阵Z=[z1z2...zm...zM],其中,m=1,2,...,M;
(6)将快拍矩阵Z中的M列数据直接相加,得到L个阵元在多普勒维的零通道数据,用此L个数据组成用于估计地杂波下俯角的一组数据向量V=[v1,v2,...,vi,...,vL]T,其中vi表示第i个阵元的多普勒零通道数据,i=1,2,...,L;
(7)根据距离单元nt和雷达距离分辨率,得到以距离单元nt为中心的分辨单元,记为nr,由目标所在的距离单元nt与雷达距离分辨率,计算分辨单元nr中的远端下俯角θf与近端下俯角θc;
(8)根据远端下俯角θf与近端下俯角θc,对数据向量V使用阵列角度估计方法,得到距离单元nt中地杂波下俯角的范围:θmin~θmax,其中,θmin=θg-ε,θmax=θg+ε,θg是由阵列角度估计方法直接估计得到的地杂波下俯角,ε为角度偏差量,其取值与该距离单元的杂噪比以及实际中关于地杂波空间分布的先验知识有关;
(9)对步骤5中的快拍矩阵Z做三脉冲对消杂波抑制,得到不含杂波的快拍矩阵U,其中,不含杂波的快拍矩阵U是L行T列的矩阵,L是阵元数,T表示三脉冲对消杂波抑制后得到的脉冲个数,T=M-2;
(10)根据不含杂波的快拍矩阵U,得到目标的协方差矩阵W,应用信源数估计方法,计算距离单元nt中所包含的空地目标个数Q;
(11)根据空地目标个数Q,估计每个目标的俯仰角θq,其中,q=1,2,...,Q;
(12)将每个目标的俯仰角θq分别与步骤8得到的地杂波下俯角范围θmin~θmax进行比较,若θmin≤θq≤θmax,则判定该目标为地面目标,否则,判定该目标为空中目标。
3.根据权利要求1所述的系留气球载相控阵雷达在多目标下分辨空地目标的方法,其中步骤(8)所述的根据远端下俯角θf与近端下俯角θc,对数据向量V使用阵列角度估计方法,得到距离单元nt中地杂波下俯角的范围:θmin~θmax,按如下步骤进行:
8a)由远端下俯角θf和近端下俯角θc,构造地杂波下俯角的搜索范围θf~θc;
其中,d为系留气球载相控阵雷达的阵元间距,λ为雷达发射信号的波长,j为虚数符号, [·]T表示转置;
其中,[·]H表示共轭转置;
8e)根据地杂波下俯角θg,得到地杂波下俯角的范围:θmin~θmax,其中,θmin=θg-ε,θmax=θg+ε,ε为角度偏差量,它的取值与该距离单元的杂噪比以及实际中关于地杂波空间分布的先验知识有关。
4.根据权利要求1所述的系留气球载相控阵雷达在多目标下分辨空地目标的方法,其中步骤(10)所述的根据不含杂波的快拍矩阵U,得到目标的协方差矩阵W,按如下公式计算:
其中,协方差矩阵W为L维方阵,L是阵元数,T表示三脉冲对消杂波抑制后得到的脉冲个数,T=M-2,[·]H表示共轭转置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310482277.8A CN103529437B (zh) | 2013-10-15 | 2013-10-15 | 系留气球载相控阵雷达在多目标下分辨空地目标的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310482277.8A CN103529437B (zh) | 2013-10-15 | 2013-10-15 | 系留气球载相控阵雷达在多目标下分辨空地目标的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103529437A true CN103529437A (zh) | 2014-01-22 |
CN103529437B CN103529437B (zh) | 2015-07-15 |
Family
ID=49931592
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310482277.8A Expired - Fee Related CN103529437B (zh) | 2013-10-15 | 2013-10-15 | 系留气球载相控阵雷达在多目标下分辨空地目标的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103529437B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103905131A (zh) * | 2014-03-11 | 2014-07-02 | 上海永畅信息科技有限公司 | 一种基于空地信道探测的高速公路车辆感知系统及感知方法 |
CN104515981A (zh) * | 2014-12-08 | 2015-04-15 | 广西大学 | 一种基于压缩感知的雷达信号处理方法和装置 |
CN105785354A (zh) * | 2016-03-11 | 2016-07-20 | 中国人民解放军国防科学技术大学 | 一种相控阵雷达探测高速运动目标的搜索参数优化方法 |
CN106199551A (zh) * | 2016-07-26 | 2016-12-07 | 深圳中钧空间科学技术有限公司 | 面阵相扫低空雷达信号处理系统 |
CN108120976A (zh) * | 2017-12-08 | 2018-06-05 | 中国船舶重工集团公司第七二四研究所 | 一种基于多普勒通道特性的地杂波谱泄露抑制方法 |
CN106093935B (zh) * | 2016-08-26 | 2018-08-14 | 西安电子科技大学 | 俯仰宽波束雷达地面和低空目标分类方法 |
CN108983214A (zh) * | 2018-05-03 | 2018-12-11 | 西安电子工程研究所 | 一种雷达导引头目标选择方法 |
CN111220955A (zh) * | 2019-11-04 | 2020-06-02 | 南京航空航天大学 | 基于垂直阵列协方差矩阵特征分解的机载气象雷达地杂波抑制方法 |
CN111600641A (zh) * | 2020-05-22 | 2020-08-28 | 电子科技大学 | 一种中继探测中的波束宽度优化方法 |
CN113219430A (zh) * | 2021-04-14 | 2021-08-06 | 加特兰微电子科技(上海)有限公司 | 地杂波处理方法、装置、计算机设备和存储介质 |
CN113391274A (zh) * | 2021-06-15 | 2021-09-14 | 电子科技大学 | 一种机载相控阵雷达低空目标最优搜索方法 |
CN113447896A (zh) * | 2021-06-07 | 2021-09-28 | 重庆大学 | 一种基于动态遮挡判断的起伏地形sar回波仿真方法 |
CN114217278A (zh) * | 2021-12-14 | 2022-03-22 | 齐鲁空天信息研究院 | 一种基于动态目标的相控阵雷达幅相校准方法及系统 |
CN114879156A (zh) * | 2022-03-18 | 2022-08-09 | 北京遥感设备研究所 | 射频探测主波束内多目标分辨方法、装置及电子设备 |
CN115308730A (zh) * | 2022-08-15 | 2022-11-08 | 中国船舶重工集团公司第七二四研究所 | 基于阵元回波能量匹配的导向矢量重构超分辨测角方法 |
CN116047462A (zh) * | 2023-03-31 | 2023-05-02 | 中国人民解放军空军预警学院 | 端射阵机载雷达最优阵元数目和阵元间距选取方法与装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6420995B1 (en) * | 1965-04-05 | 2002-07-16 | Systems Information And Electronic Systems Integration, Inc. | Radar and IFF system |
CN101465071A (zh) * | 2009-01-08 | 2009-06-24 | 上海交通大学 | 多平台目标跟踪与分布交互仿真系统 |
US7714782B2 (en) * | 2004-01-13 | 2010-05-11 | Dennis Willard Davis | Phase arrays exploiting geometry phase and methods of creating such arrays |
CN103018727A (zh) * | 2011-09-27 | 2013-04-03 | 中国科学院电子学研究所 | 一种基于样本训练的机载雷达非平稳杂波抑制方法 |
CN103197301A (zh) * | 2013-03-19 | 2013-07-10 | 中国人民解放军海军航空工程学院 | 海面微动目标Radon-线性正则变换长时间相参积累检测方法 |
-
2013
- 2013-10-15 CN CN201310482277.8A patent/CN103529437B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6420995B1 (en) * | 1965-04-05 | 2002-07-16 | Systems Information And Electronic Systems Integration, Inc. | Radar and IFF system |
US7714782B2 (en) * | 2004-01-13 | 2010-05-11 | Dennis Willard Davis | Phase arrays exploiting geometry phase and methods of creating such arrays |
CN101465071A (zh) * | 2009-01-08 | 2009-06-24 | 上海交通大学 | 多平台目标跟踪与分布交互仿真系统 |
CN103018727A (zh) * | 2011-09-27 | 2013-04-03 | 中国科学院电子学研究所 | 一种基于样本训练的机载雷达非平稳杂波抑制方法 |
CN103197301A (zh) * | 2013-03-19 | 2013-07-10 | 中国人民解放军海军航空工程学院 | 海面微动目标Radon-线性正则变换长时间相参积累检测方法 |
Non-Patent Citations (1)
Title |
---|
虞庆庆等: "系留气球载雷达系统雷达结构总体分析", 《现代雷达》 * |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103905131B (zh) * | 2014-03-11 | 2016-11-02 | 上海永畅信息科技有限公司 | 一种基于空地信道探测的高速公路车辆感知系统及感知方法 |
CN103905131A (zh) * | 2014-03-11 | 2014-07-02 | 上海永畅信息科技有限公司 | 一种基于空地信道探测的高速公路车辆感知系统及感知方法 |
CN104515981A (zh) * | 2014-12-08 | 2015-04-15 | 广西大学 | 一种基于压缩感知的雷达信号处理方法和装置 |
CN105785354A (zh) * | 2016-03-11 | 2016-07-20 | 中国人民解放军国防科学技术大学 | 一种相控阵雷达探测高速运动目标的搜索参数优化方法 |
CN105785354B (zh) * | 2016-03-11 | 2018-02-02 | 中国人民解放军国防科学技术大学 | 一种相控阵雷达探测高速运动目标的搜索参数优化方法 |
CN106199551A (zh) * | 2016-07-26 | 2016-12-07 | 深圳中钧空间科学技术有限公司 | 面阵相扫低空雷达信号处理系统 |
CN106093935B (zh) * | 2016-08-26 | 2018-08-14 | 西安电子科技大学 | 俯仰宽波束雷达地面和低空目标分类方法 |
CN108120976A (zh) * | 2017-12-08 | 2018-06-05 | 中国船舶重工集团公司第七二四研究所 | 一种基于多普勒通道特性的地杂波谱泄露抑制方法 |
CN108120976B (zh) * | 2017-12-08 | 2021-02-23 | 中国船舶重工集团公司第七二四研究所 | 一种基于多普勒通道特性的地杂波谱泄露抑制方法 |
CN108983214A (zh) * | 2018-05-03 | 2018-12-11 | 西安电子工程研究所 | 一种雷达导引头目标选择方法 |
CN108983214B (zh) * | 2018-05-03 | 2022-04-08 | 西安电子工程研究所 | 一种雷达导引头目标选择方法 |
CN111220955A (zh) * | 2019-11-04 | 2020-06-02 | 南京航空航天大学 | 基于垂直阵列协方差矩阵特征分解的机载气象雷达地杂波抑制方法 |
CN111220955B (zh) * | 2019-11-04 | 2022-10-11 | 南京航空航天大学 | 基于垂直阵列协方差矩阵特征分解的机载气象雷达地杂波抑制方法 |
CN111600641B (zh) * | 2020-05-22 | 2021-07-02 | 电子科技大学 | 一种中继探测中的波束宽度优化方法 |
CN111600641A (zh) * | 2020-05-22 | 2020-08-28 | 电子科技大学 | 一种中继探测中的波束宽度优化方法 |
CN113219430A (zh) * | 2021-04-14 | 2021-08-06 | 加特兰微电子科技(上海)有限公司 | 地杂波处理方法、装置、计算机设备和存储介质 |
CN113447896A (zh) * | 2021-06-07 | 2021-09-28 | 重庆大学 | 一种基于动态遮挡判断的起伏地形sar回波仿真方法 |
CN113447896B (zh) * | 2021-06-07 | 2023-03-14 | 重庆大学 | 一种基于动态遮挡判断的起伏地形sar回波仿真方法 |
CN113391274A (zh) * | 2021-06-15 | 2021-09-14 | 电子科技大学 | 一种机载相控阵雷达低空目标最优搜索方法 |
CN113391274B (zh) * | 2021-06-15 | 2022-08-05 | 电子科技大学 | 一种机载相控阵雷达低空目标最优搜索方法 |
CN114217278A (zh) * | 2021-12-14 | 2022-03-22 | 齐鲁空天信息研究院 | 一种基于动态目标的相控阵雷达幅相校准方法及系统 |
CN114879156A (zh) * | 2022-03-18 | 2022-08-09 | 北京遥感设备研究所 | 射频探测主波束内多目标分辨方法、装置及电子设备 |
CN114879156B (zh) * | 2022-03-18 | 2024-09-24 | 北京遥感设备研究所 | 射频探测主波束内多目标分辨方法、装置及电子设备 |
CN115308730A (zh) * | 2022-08-15 | 2022-11-08 | 中国船舶重工集团公司第七二四研究所 | 基于阵元回波能量匹配的导向矢量重构超分辨测角方法 |
CN115308730B (zh) * | 2022-08-15 | 2024-10-25 | 中国船舶集团有限公司第七二四研究所 | 基于阵元回波能量匹配的导向矢量重构超分辨测角方法 |
CN116047462A (zh) * | 2023-03-31 | 2023-05-02 | 中国人民解放军空军预警学院 | 端射阵机载雷达最优阵元数目和阵元间距选取方法与装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103529437B (zh) | 2015-07-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103529437B (zh) | 系留气球载相控阵雷达在多目标下分辨空地目标的方法 | |
CN103869311B (zh) | 实波束扫描雷达超分辨成像方法 | |
Capraro et al. | Implementing digital terrain data in knowledge-aided space-time adaptive processing | |
CN101915911B (zh) | 基于相消积累空时谱的空间任意构型分布式sar动目标参数估计方法 | |
CN102156279B (zh) | 基于mimo的双基地雷达地面动目标检测方法 | |
CN105223560B (zh) | 基于杂波俯仰方位谱稀疏恢复的机载雷达目标检测方法 | |
CN110412559A (zh) | 分布式无人机mimo雷达的非相参融合目标检测方法 | |
CN103197294B (zh) | 多频融合最大似然低空目标仰角估计方法 | |
CN109358325B (zh) | 起伏地形背景下雷达高度表的地形反演方法 | |
KR102013205B1 (ko) | 레이더 신호 처리 모의 장치 및 방법 | |
CN106483516A (zh) | 基于先验知识的雷达杂波空时自适应处理方法 | |
CN113253230B (zh) | 基于子孔径处理的天基预警雷达空中动目标检测方法及系统 | |
CN102998672A (zh) | 基于相干化处理的步进频率isar成像方法 | |
CN112612006B (zh) | 基于深度学习的机载雷达非均匀杂波抑制方法 | |
CN104166129A (zh) | 一种实波束雷达迭代最小均方误差角超分辨方法 | |
CN104166134A (zh) | 一种实波束前视扫描雷达目标二维定位方法 | |
Oveis et al. | Coherent method for ground‐moving target indication and velocity estimation using Hough transform | |
CN109061640B (zh) | 一种用于顺轨干涉sar海流反演的方位模糊抑制方法 | |
Watson et al. | Non-line-of-sight radar | |
CN103630899A (zh) | 地面运动目标高分辨雷达压缩感知成像的方法 | |
Page et al. | Detection and tracking of moving vehicles with Gotcha radar systems | |
Gao et al. | Static background removal in vehicular radar: Filtering in azimuth-elevation-doppler domain | |
CN100585429C (zh) | 一种基于非直线天线阵列的无源通道校正方法 | |
CN114355306B (zh) | 一种利用小孔径高频雷达海洋回波反演浪高的方法 | |
CN112835006B (zh) | 一种基于帧间积累的跟踪雷达海上小目标检测方法及系统 |
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: 20150715 Termination date: 20211015 |
|
CF01 | Termination of patent right due to non-payment of annual fee |