CN112986939B - 一种多径环境下机载相控阵雷达信号检测方法 - Google Patents
一种多径环境下机载相控阵雷达信号检测方法 Download PDFInfo
- Publication number
- CN112986939B CN112986939B CN202110144082.7A CN202110144082A CN112986939B CN 112986939 B CN112986939 B CN 112986939B CN 202110144082 A CN202110144082 A CN 202110144082A CN 112986939 B CN112986939 B CN 112986939B
- Authority
- CN
- China
- Prior art keywords
- multipath
- echo
- direct
- target
- signal
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000001514 detection method Methods 0.000 claims abstract description 66
- 238000009825 accumulation Methods 0.000 claims abstract description 26
- 238000005070 sampling Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims description 13
- 230000005540 biological transmission Effects 0.000 claims description 6
- 230000006835 compression Effects 0.000 claims description 6
- 238000007906 compression Methods 0.000 claims description 6
- RVRCFVVLDHTFFA-UHFFFAOYSA-N heptasodium;tungsten;nonatriacontahydrate Chemical compound O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.[Na+].[Na+].[Na+].[Na+].[Na+].[Na+].[Na+].[W].[W].[W].[W].[W].[W].[W].[W].[W].[W].[W] RVRCFVVLDHTFFA-UHFFFAOYSA-N 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000001427 coherent effect Effects 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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/411—Identification of targets based on measurements of radar reflectivity
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
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
本发明公开了一种多径环境下机载相控阵雷达信号检测方法,包括:得到直达回波功率和多径回波功率;根据直达回波功率和多径回波功率得到综合多径反射系数;得到第m个阵元的快时间采样的目标回波数据;得到同一目标的直达回波距离和多径回波距离对应的目标直达数据和目标多径数据;得到快时间数据;对每个脉冲的所述快时间数据进行恒虚警检测以得到检测结果。本发明方法分析了在多径环境下机载相控阵雷达的综合多径反射系数的特点;分析了在多径环境下机载相控阵雷达的直达回波和多径回波的复包络特点;利用非相干积累方法提高了机载相控阵雷达在多径环境下在目标距离单元处的输出增益,提高了***的检测概率。
Description
技术领域
本发明属于机载雷达信号处理技术领域,具体涉及一种多径环境下机载相控阵雷达信号检测方法。
背景技术
雷达的工作机制是先通过发射电磁波,然后接收从目标表面后向散射回的电磁回波,再对电磁回波数据进行分析、检测的一种工作过程。其中对电磁回波数据进行分析、检测过程的目的是检测目标,然后对目标进行参数估计,得到目标的方位、高度以及速度等信息。因而信号检测是雷达进行信号处理的重要内容。
机载雷达是将天线阵列放置在空中飞行的载体上的一种雷达体制。当天线阵放置在高空中飞行的载机上时,此时机载雷达的视角大,可以探测到地基雷达无法探测到的目标。高精度的机载雷达可以快速精确锁定地面目标,配合地面力量摧毁敌方目标。因此机载雷达性能好坏是未来战场上夺取制空权的关键指标。
多径现象是指电磁波从雷达某个固定方向发射出去后,经过目标表面的后向散射后,有两个或多个目标方向的回波信号反射回来。如果雷达信号处理时并不知道此时的回波中有多径回波信号时,将导致在信号处理时无法和接收到的回波信号模型完全匹配,直接影响目标的检测性能,一般地基米波雷达在低仰角测高时容易发生多径现象。实际上,在机载雷达中也经常有多径现象发生,特别是机载雷达在平静的海面上飞行时,平静海面相当于一个强的反射体,目标的直达信号除了通过直达路径返回到雷达外,还通过海面的强镜面反射回雷达。如果机载雷达仍然认为此时只有直达信号,那么在多径信号的影响下,最后将导致雷达可能无法正确的检测到目标信号,直接影响到后续雷达***对目标的定位、跟踪。
机载雷达在多径环境下,还存在着直达信号和多径回波在目标处的夹角大的特点。直达回波和多径回波之间的几何关系决定了直达回波信号和多径回波信号在目标处的夹角大小,过大的直达回波和多径回波在目标处的夹角,就有可能导致在多径方向上等效的雷达散射截面积大于直达方向上等效的散射截面积,进而导致多径回波方向的接收功率大于直达回波方向上的接收功率,因此在信号检测时若仍然只用直达信号作为检测数据,那么在恒虚警检测(CFAR,Constant False-Alarm Rate)时很有可能导致漏检。
因此,机载相控阵雷达在多径环境下的信号检测面临的情况更复杂,是一个需要重点关注的问题。
发明内容
为了解决现有技术中存在的上述问题,本发明提供了一种多径环境下机载相控阵雷达信号检测方法。本发明要解决的技术问题通过以下技术方案实现:
一种多径环境下机载相控阵雷达信号检测方法,包括如下步骤:
步骤1、根据机载雷达的参数得到直达回波功率Pd和多径回波功率Pi;
步骤2、基于综合多径反射系数计算模型,根据所述直达回波功率Pd和所述多径回波功率Pi得到综合多径反射系数;
步骤3、对每个阵元的回波数据进行脉冲压缩,根据第m个阵元的直达回波信号、第m个阵元的多径回波信号和所述综合多径反射系数得到第m个阵元的快时间采样的目标回波数据;
步骤4、根据所述目标回波数据和仰角波位得到目标的多径回波数据所在的距离单元,以得到同一目标的直达回波距离和多径回波距离对应的目标直达数据和目标多径数据;
步骤5、基于平方律检波,将所述目标多径数据加到直达距离单元中,实现目标直达数据和目标多径数据的非相干积累,得到快时间数据;
步骤6、对每个脉冲的所述快时间数据进行恒虚警检测以得到检测结果。
在本发明的一个实施例中,所述直达回波功率Pd和所述多径回波功率Pi的计算公式分别为:
其中,Pt,d、Pt,i分别是雷达直达方向的发射功率和多径方向的发射功率,且Pt,d=Pt,i,Gt,d、Gt,i都是雷达直达方向的发射天线增益和多径方向的发射天线增益,且Gt,d=Gt,i,Gr,d、Gr,i分别是直达方向的接收增益和多径方向的接收增益,且Gr,i=εGr,d,ε是接收增益的损耗,|ε|<1,λ是波长,σd和σi分别是直达方向和多径方向的散射截面积,Rd为目标的直达回波距离,Ri为目标的多径回波距离。
在本发明的一个实施例中,所述综合多径反射系数计算模型为:
其中,ρ为综合多径反射系数。
在本发明的一个实施例中,所述目标回波数据的计算公式为:
Xm(t)=sm(t-τrd(θd,fd))+ρsm(t-τri(θi,fi))
其中,Xm(t)是第m个阵元的快时间采样的目标回波数据,sm(t-τrd(θd,fd))是第m个阵元的直达回波信号,ρsm(t-τri(θi,fi))是第m个阵元的多径回波信号,t是时间序列,τrd是直达回波信号的时延,θd是直达回波信号的仰角,fd是直达回波信号的多普勒频率,ρ是综合多径反射系数,τri是多径回波信号的时延,θi是多径回波信号的仰角,fi是多径回波信号的多普勒频率。
在本发明的一个实施例中,所述多径回波信号与所述直达回波信号的时延差Δτ为:
其中,τi是多径回波信号与参考阵元的时延,τd是直达回波信号回波到阵列的参考阵元的时延,Ri为目标的多径回波距离,Rd为目标的直达回波距离,c是光速;
并且,
其中,gm(t-τrd(θd,fd))是直达回波信号的复包络,gm(t-τri(θi,fi))是多径回波信号的复包络。
在本发明的一个实施例中,所述快时间数据的计算公式为:
|Xd_i|2=|Xd|2+|Xi|2
其中,|Xd_i|2每个脉冲的快时间数据,|Xd|2和|Xi|2分别是每个脉冲在直达距离单元的输出和多径距离单元的输出。
在本发明的一个实施例中,所述恒虚警为平均恒虚警。
在本发明的一个实施例中,所述恒虚警的标称因子T为:
其中,Pfa是虚警概率,N是参考单元数。
在本发明的一个实施例中,在所述步骤6之后,还包括:
步骤7、计算单脉冲非相干积累后的检测概率。
在本发明的一个实施例中,所述检测概率的计算公式为:
其中,Pd是检测概率,S′是阈值,Xi是第i个快时间数据,T是标称因子,N是参考单元数,μ为背景噪声功率,n′=2N,z是计算检测概率的中间变量,Ps是直达和多径回波功率大小,A是信号幅度大小,ρ是综合多径反射系数,QM(·)是Marcum的Q函数,I1(·)是修正的一阶贝塞尔函数。
本发明的有益效果:
本发明方法给出了多径环境下机载相控阵雷达的信号模型;分析了在多径环境下机载相控阵雷达的综合多径反射系数的特点;分析了在多径环境下机载相控阵雷达的直达回波和多径回波的复包络特点;利用非相干积累方法提高了机载相控阵雷达在多径环境下在目标距离单元处的输出增益,进而提高了***的检测概率,同时给出了检测概率的表达式。
以下将结合附图及实施例对本发明做进一步详细说明。
附图说明
图1是本发明实施例提供的一种多径环境下机载相控阵雷达信号检测方法的流程示意图;
图2是本发明实施例提供的一种多径环境下机载相控阵雷达的信号模型图;
图3是本发明实施例提供的一种多径环境下机载相控阵雷达非相干积累的信号检测流程图;
图4是本发明实施例提供的一种多径环境下直达回波信号和多径回波信号在目标处的夹角图;
图5是本发明实施例提供的一种单脉冲下原始方法和非相干积累方法在目标处的输出增益图;
图6是本发明实施例提供的一种随着综合多径反射系数变化的非相干积累信号检测概率图;
图7是本发明实施例提供的一种随着虚警概率变化的非相干积累信号检测概率图。
具体实施方式
下面结合具体实施例对本发明做进一步详细的描述,但本发明的实施方式不限于此。
实施例一
请参见图1,图1是本发明实施例提供的一种多径环境下机载相控阵雷达信号检测方法的流程示意图。本实施例提供一种多径环境下机载相控阵雷达信号检测方法,该多径环境下机载相控阵雷达信号检测方法包括步骤1-步骤6,其中:
步骤1、根据机载雷达的参数得到直达回波功率Pd和多径回波功率Pi。
具体地,机载雷达的多径问题相比于地基雷达而言,最大的特点就是直达回波信号和多径回波信号之间的距离差大,导致直达回波信号和多径回波信号的数据不在同一距离单元中,并且导致了直达回波信号和多径回波信号在目标处的夹角变大,过大的夹角就会导致直达回波信号等效的散射截面积(RCS,Radar Cross Section)比多径回波信号等效的RCS小,进而导致多径回波信号功率大于直达回波信号功率的情况。
当电磁波辐射到目标后,面对的是一个复杂的散射体,散射截面积不仅和电磁波照射到目标的角度有关,还和飞行姿态、目标的表面材料等有关。因此很难通过一个固定的模型来表征直达回波信号和多径回波信号的回波功率大小,但是它们之间的关系可以通过经典的雷达方程推导出来。接下来通过雷达方程分析在多径环境下机载相控阵雷达的综合多径反射系数特点。
其中,直达回波功率Pd和多径回波功率Pi的计算公式分别为:
其中,Pt,d、Pt,i分别是雷达直达方向的发射功率和多径方向的发射功率,且Pt,d=Pt,i,Gt,d、Gt,i都是雷达直达方向的发射天线增益和多径方向的发射天线增益,且Gt,d=Gt,i,Gr,d、Gr,i分别是直达方向的接收增益和多径方向的接收增益,且Gr,i=εGr,d,ε是接收增益的损耗,|ε|<1,λ是波长,σd和σi分别是直达方向和多径方向的散射截面积,Rd为目标的直达回波距离,Ri为目标的多径回波距离。
步骤2、基于综合多径反射系数计算模型,根据直达回波功率Pd和多径回波功率Pi得到综合多径反射系数。
具体地,综合多径反射系数计算模型为:
其中,ρ为综合多径反射系数。
也就是说,|ρ|是否大于1由上式决定。ε主要由地面反射区的性质决定,在地基雷达中,一般米波雷达低仰角测高时会出现多径现象,此时的σi/σd和都约为1,因此此时多径的反射系数的幅度主要由反射面的损耗决定。在机载相控阵雷达的环境下,此时σi/σd和/>不再等效为1了,需要具体考虑。
步骤3、对每个阵元的回波数据进行脉冲压缩,根据第m个阵元的直达回波信号、第m个阵元的多径回波信号和综合多径反射系数得到第m个阵元的快时间采样的目标回波数据。
具体地,图2是机载相控阵雷达存在多径时的直达回波信号和多径回波信号的信号模型图。通过图2可以看到,直达回波信号和多径回波信号在不同的距离单元,因此在快时间采样后,对采样的回波数据进行脉冲压缩,可以检测到直达信号距离Rd和多径信号距离Ri。进而确定直达回波信号与接收天线中参考阵元的时延τd=2Rd/c以及多径回波信号与接收天线中参考阵元的时延τi=(Rd+R)i/c,其中,c是光速,Ri=R1+R2,R1、R2请参见图2。
接下来分析机载相控阵雷达在多径环境下直达回波信号和多径回波信号的信号模型特点。在不考虑噪声的情况下,第m个阵元的快时间采样的目标回波数据为:
Xm(t)=sm(t-τrd(θd,fd))+ρsm(t-τri(θi,fi))
其中,Xm(t)是第m个阵元的快时间采样的目标回波数据,受直达回波RCS、多径回波RCS、信号的损耗等影响,sm(t-τrd(θd,fd))是第m个阵元的直达回波信号,ρsm(t-τri(θi,fi))是第m个阵元的多径回波信号,t是时间序列,τrd是直达回波信号的时延,θd是直达回波信号的仰角,fd是直达回波信号的多普勒频率,ρ是综合多径反射系数,τri是多径回波信号的时延,θi是多径回波信号的仰角,fi是多径回波信号的多普勒频率。
其中,第m个阵元的直达回波信号sm(t-τrd(θd,fd))可以进一步表示为:
其中,gm(t-τrd(θd,fd))是发射LFM(线性调频)信号的第m个阵元的直达回波信号的复包络,fc是发射载频,τsm(θd)是直达回波信号在第m个阵元与参考阵元间的波程差产生的时延,τd,k(fd)是第k个脉冲与参考脉冲间在第τd时刻时直达回波信号与雷达之间的径向速度Vd引起的路程差导致的时延,c是光速,是相邻阵元间的距离,τd是直达回波信号到阵列的参考阵元的时延,gm(t)是发射LFM信号的复包络,θd是直达回波信号的仰角,fd是直达回波信号的多普勒频率。
其中,第m个阵元的多径回波信号ρsm(t-τri(θi,fi))可以进一步表示为:
其中,gm(t-τri(θi,fi))是发射LFM信号的第m个阵元的多径回波信号的复包络,θi是多径回波信号的仰角,fi是多径回波信号的多普勒频率,τsm(θi)是多径回波信号在第m个阵元与参考阵元间的波程差产生的时延,τi,k(fi)是第k个脉冲与参考脉冲间在第τi时刻多径回波信号与雷达之间的径向速度Vi引起的路程差导致的时延,τi是多径回波信号与参考阵元的时延。
根据直达回波信号和多径回波信号的距离,可知多径回波信号与直达回波信号的时延差Δτ为:
其中,τi是多径回波信号与参考阵元的时延,τd是直达回波信号回波到阵列的参考阵元的时延,Ri为目标的多径回波距离,Rd为目标的直达回波距离,c是光速。
由于直达回波信号和多径回波信号在不同的距离单元中,因此多普勒频率fd、fi可以通过慢时间采样数据求得。因此通过直达信号模型、多径信号模型、估计的多普勒频率和时延差Δτ,可以确定直达回波信号和多径回波信号之间的信号相位关系。然而,在机载雷达中,由于直达回波信号和多径回波信号的距离差很大,因此Δτ>1/fs,其中fs是采样频率。这样就导致在离散采样时,直达回波信号与多径回波信号的复包络延迟导致复包络中的随机初相是未知的,进而导致它们的复包络不再近似相等,即因此即使能确定直达回波信号和多径回波信号的相位关系,但是由于复包络受时延的影响,将不能直接通过联合导向矢量类的方法检测目标信号了。
步骤4、根据目标回波数据和仰角波位得到目标的多径回波数据所在的距离单元,以得到同一目标的直达回波距离和多径回波距离对应的目标直达数据和目标多径数据。
具体地,在实际信号检测时,是无法直接知道直达回波信号的多径回波所在的距离单元,但是可以通过目标扫描的仰角波位和脉冲压缩后的距离信息对目标多径距离单元进行预判。具体过程是,通过对目标回波数据进行匹配滤波后,可以将直达距离和多径距离均检测出来,根据目标回波数据、***工作时的仰角波位为机载雷达高度为H、脉冲压缩后的距离为Rd,结合图2中的几何关系以及《雷达手册(第二版)》第二版第49面到第51面的内容,可以计算出目标对应的距离Ri,然后将Ri对应的距离单元中的数据经过平方律检波后加到直达距离单元中,具体的工作流程见图3。
进一步地,利用直达回波数据和多径回波数据的非相干积累。检测到直达回波信号和多径回波信号对应的距离单元Rd和Ri,然后在快时间域对一个时间段的所有数据进行平方律检波,得到每个脉冲在直达距离单元的输出|Xd|2和多径距离单元的输出|Xi|2,然后对直达回波数据和多径回波数据进行非相干积累得到快时间数据,即:
|Xd_i|2=|Xd|2+|Xi|2
在单脉冲检测中将每个脉冲对应的|Xd_i|2数据整理到直达距离单元中,因此直达距离单元中的目标信号功率为Ps=(1+|ρ|2)A2,A2是目标信号包络的平方。
步骤6、对每个脉冲的快时间数据进行恒虚警检测以得到检测结果。
一般雷达在工作时,一个驻留内会发射多个脉冲,但是为了叙述方便,本实施例只讨论一个驻留内只发射一个脉冲的情况。
在图3中,将非相干积累后的快时间数据进行恒虚警检测(CFAR,Constant False-Alarm Rate)。CFAR的自适应判决准则为:
其中,H1表示目标存在,H0表示目标不存在,D表示检测统计量,T表示标称化因子,Z表示估计出的杂波、噪声功率,也就是对杂波、噪声包络的估计。为了简化过程,假设机载相控阵雷达在多径环境下的回波数据中除了直达回波信号和多径波信号外,只有高斯白噪声。
图3中采用的CFAR方法是单元平均CFAR(CA-CFAR,CA-CFAR单元平均恒虚警检测)。假设CA-CFAR检测器的背景噪声都服从复高斯分布,因而背景噪声的数据包络x都是服从瑞利分布的变量。经过平方律检波器后,背景噪声的变量Xi=xi 2,i=1,2,…,L服从指数分布,即Xi~G(1,μ),其中G(α,β)是伽马函数,且α=1时为指数分布,μ是噪声功率,L是快时间采样数。
在CA-CFAR检测器中,假设参考单元个数为N=2n,那么经过平方律检波器后的背景噪声的平方和为:
因此Z的概率密度函数(PDF)为:
其中,*是卷积符号,f(·)表示概率密度函数。结合伽马函数性质,可以进一步得到Z~G(N,μ)。按照虚警概率的定义,对Z的PDF从阈值S到+∞进行积分,可以得到:
其中,Pfa是虚警概率,,是不完全伽马函数,τ是计算虚警概率的中间变量。对于N=1的独立数据样本采样,可得:
Pfa(S)=e-S/μ
因此在N=1时的采样,可以得到阈值S和虚警概率Pfa的关系为:
S=-μlnPfa
在CA-CFAR检测器中,阈值S的表达式为:
由于标称因子T为常数,因此阈值S的概率密度函数和相同,N个指数分布变量的和为伽马分布G(N,μ)。因此可以得到虚警概率的数学期望P′fa为:
在雷达***设计时固定了虚警概率的数学期望P′fa后,就可以求出标称因子T为:
在一个具体实施例中,在所述步骤6之后,还包括:
步骤7、计算单脉冲非相干积累后的检测概率。
从步骤6的标称因子计算公式可以看到虚警概率不依赖背景噪声功率μ,这也是CFAR的重要特点。
进一步可以计算出单个目标的检测概率为:
其中,Pd是检测概率,λ′是功率信噪比。
机载相控阵雷达存在多径时,直达回波信号和多径回波信号的角度和多普勒频率不同,并且单个脉冲快时间采样的回波信号复包络不对齐,因此直达回波信号和多径回波信号无法进行相干检测,但是机载雷达在地面反射区较强的环境下,多径现象是存在的,特别是平静的海面上。机载相控阵雷达在多径环境下探测目标时,由于多径信号实际上是直达回波信号的“镜像”,因此考虑通过非相干积累的方式,将多径回波信号对应的距离单元中的数据在平方律检波后积累到直达距离单元中,然后再对非相干积累后的数据进行CA-CFAR。
具体实施步骤是先进行脉冲压缩,检测到直达回波信号和多径回波信号对应的距离单元Rd和Ri,然后在快时间域对数据进行平方律检波,得到每个脉冲在直达距离单元的输出|Xd|2和多径距离单元的输出|Xi|2,然后对直达数据和多径数据进行非相干积累,得到|Xd_i|2=|Xd|2+|Xi|2。
在非相干积累完成后,对非相干积累后的数据进行CA-CFAR,最后对CA-CFAR处理后的数据按照自适应判决准则进行判断,检测出目标。
通过上面分析可以发现在非相干检测时,虚警概率Pfa和标称因子T之间的关系没有变,只是单脉冲的目标检测概率Pd发生了变化。具体来说,雷达***设定了虚警概率Pfa后,CA-CFAR的标称因子T为
此时直达回波信号和多径回波信号的单脉冲非相干积累相当于两个目标的包络平方和,因此检测概率为:
其中,Pd是检测概率,S′是阈值,Xi是第i个快时间数据,T是标称因子,N是参考单元数,μ为背景噪声功率,n′=2N,z是计算检测概率的中间变量,Ps是直达和多径回波功率大小,A是信号幅度大小,ρ是综合多径反射系数,QM(·)是Marcum的Q函数,I1(·)是修正的一阶贝塞尔函数。
1.仿真条件
1.1机载相控阵雷达有M=20个收发共置天线,机载相控阵雷达的高度为8公里,雷达工作的波长λ为1米,目标的直达回波信号仰角为θd=7°,对应的多径回波信号仰角为θi=-25.87°。直达回波信号的距离Rd∈[25,50]km。
1.2在1.1的条件下,直达回波信号的距离为Rd=45km,ε=0.8,假定σi/σd∈[0.45,4.1],进一步通过计算可以得到ρ=[-1.49,-0.49]。因此在不失一般性的情况下,假设回波信号的信噪比为12dB,多径反射系数ρ=-1.2,相控阵发射线性调频信号,时宽为Tp=40us,带宽为Bs=4MHz,脉冲重复频率为fr=2000Hz,载频为fc=300MHz,光速为c=3×108m/s,采样频率为fs=Bs,仿真了在单脉冲检测中,保护单元为1个,参考单元N=16个时,快时间采样上常规CA-CFAR和改进的CA-CFAR的输出增益对比。
1.3在1.2的条件下,仿真了Pfa=10-6,ρ=[-1.5,-1,-0.5]取不同值时,随着输入信噪比变化的检测概率。
1.4在1.2的条件下,仿真了ρ=-1.2,Pfa=[10-5,10-6,10-7,10-8]时,随着输入信噪比变化的检测概率。
2.仿真内容
2.1图4是直达回波信号与多径回波信号在目标处的夹角,横坐标是直达回波信号的距离,纵坐标就是它们在目标处的夹角。通过图4可以看到,夹角都超过了17度,根据《雷达手册(第二版)》P427中提供的数据,可以知道这么大的角度是很有可能出现反射系数|ρ|大于1的情况,进而说明前面分析反射系数的情况是符合实际情况的。
2.2图5中的左图中的原始方法就是指单脉冲检测中常规的CA-CFAR方法,右图中的非相干积累方法指的是单脉冲检测中改进的CA-CFAR方法。通过1.2的条件可以看到,直达距离单元在第1200距离单元中。对比图5的左图和右图可以看到,在单脉冲下非相干积累把多径距离单元的信号能量也集中到直达距离单元中,因此右图在直达距离单元处超过了CA-CFAR的门限值,而左图的包络平方没有超过CA-CFAR的门限值,导致无法检测到目标。
2.3图6中仿真了随着ρ分别为[-1.5,-1,-0.5]时的检测概率,通过与原始方法的对比可以发现,没有采用非相干积累时,单脉冲快采样数据在检测时在直达距离单元只有直达回波信号,而采用了非相干积累时在直达距离单元中不仅有直达回波信号的能量,而且还有多径回波信号的能量,因此用非相干积累的检测概率都高于没有用非相干积累的原始方法,并且随着|ρ|的增加,非相干积累的信号能量越多,因此对应的检测概率也越高。
2.4图7中仿真了当ρ=-1.2时随着虚警概率Pfa的增加,虚警概率的变化情况。从图中可以看到采用非相干积累后的检测概率比原始方法的检测概率要高,并且随着Pfa的增加,对应的检测概率也会更大,这点也和CFAR思想相吻合。结合图6的结果,可以发现在多径环境下,机载相控阵雷达在信号检测时用非相干积累的方法可以有效的提高目标的检测概率。
本发明方法给出了多径环境下机载相控阵雷达的信号模型;分析了在多径环境下机载相控阵雷达的综合多径反射系数的特点;分析了在多径环境下机载相控阵雷达的直达回波和多径回波的复包络特点;利用非相干积累方法提高了机载相控阵雷达在多径环境下在目标距离单元处的输出增益,进而提高了***的检测概率,同时给出了检测概率的表达式。
在本发明的描述中,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中,“多个”的含义是两个或两个以上,除非另有明确具体的限定。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特数据点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特数据点可以在任何的一个或多个实施例或示例中以合适的方式结合。此外,本领域的技术人员可以将本说明书中描述的不同实施例或示例进行接合和组合。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (10)
1.一种多径环境下机载相控阵雷达信号检测方法,其特征在于,包括如下步骤:
步骤1、根据机载雷达的参数得到直达回波功率Pd和多径回波功率Pi;
步骤2、基于综合多径反射系数计算模型,根据所述直达回波功率Pd和所述多径回波功率Pi得到综合多径反射系数;
步骤3、对每个阵元的回波数据进行脉冲压缩,根据第m个阵元的直达回波信号、第m个阵元的多径回波信号和所述综合多径反射系数得到第m个阵元的快时间采样的目标回波数据;
步骤4、根据所述目标回波数据和仰角波位得到目标的多径回波数据所在的距离单元,以得到同一目标的直达回波距离和多径回波距离对应的目标直达数据和目标多径数据;
步骤5、基于平方律检波,将所述目标多径数据加到直达距离单元中,实现目标直达数据和目标多径数据的非相干积累,得到快时间数据;
步骤6、对每个脉冲的所述快时间数据进行恒虚警检测以得到检测结果。
2.根据权利要求1所述的多径环境下机载相控阵雷达信号检测方法,其特征在于,所述直达回波功率Pd和所述多径回波功率Pi的计算公式分别为:
其中,Pt,d、Pt,i分别是雷达直达方向的发射功率和多径方向的发射功率,且Pt,d=Pt,i,Gt,d、Gt,i都是雷达直达方向的发射天线增益和多径方向的发射天线增益,且Gt,d=Gt,i,Gr,d、Gr,i分别是直达方向的接收增益和多径方向的接收增益,且Gr,i=εGr,d,ε是接收增益的损耗,|ε|<1,λ是波长,σd和σi分别是直达方向和多径方向的散射截面积,Rd为目标的直达回波距离,Ri为目标的多径回波距离。
3.根据权利要求2所述的多径环境下机载相控阵雷达信号检测方法,其特征在于,所述综合多径反射系数计算模型为:
其中,ρ为综合多径反射系数。
4.根据权利要求1所述的多径环境下机载相控阵雷达信号检测方法,其特征在于,所述目标回波数据的计算公式为:
Xm(t)=sm(t-τrd(θd,fd))+ρsm(t-τri(θi,fi))
其中,Xm(t)是第m个阵元的快时间采样的目标回波数据,sm(t-τrd(θd,fd))是第m个阵元的直达回波信号,ρsm(t-τri(θi,fi))是第m个阵元的多径回波信号,t是时间序列,τrd是直达回波信号的时延,θd是直达回波信号的仰角,fd是直达回波信号的多普勒频率,ρ是综合多径反射系数,τri是多径回波信号的时延,θi是多径回波信号的仰角,fi是多径回波信号的多普勒频率。
5.根据权利要求1所述的多径环境下机载相控阵雷达信号检测方法,其特征在于,所述多径回波信号与所述直达回波信号的时延差Δτ为:
其中,τi是多径回波信号与参考阵元的时延,τd是直达回波信号回波到阵列的参考阵元的时延,Ri为目标的多径回波距离,Rd为目标的直达回波距离,c是光速;
并且,
其中,gm(t-τrd(θd,fd))是直达回波信号的复包络,gm(t-τri(θi,fi))是多径回波信号的复包络。
6.根据权利要求1所述的多径环境下机载相控阵雷达信号检测方法,其特征在于,所述快时间数据的计算公式为:
|Xd_i|2=|Xd|2+|Xi|2
其中,|Xd_i|2每个脉冲的快时间数据,|Xd|2和|Xi|2分别是每个脉冲在直达距离单元的输出和多径距离单元的输出。
7.根据权利要求1所述的多径环境下机载相控阵雷达信号检测方法,其特征在于,所述恒虚警为平均恒虚警。
8.根据权利要求6所述的多径环境下机载相控阵雷达信号检测方法,其特征在于,所述恒虚警的标称因子T为:
其中,Pfa是虚警概率,N是参考单元数。
9.根据权利要求1所述的多径环境下机载相控阵雷达信号检测方法,其特征在于,在所述步骤6之后,还包括:
步骤7、计算单脉冲非相干积累后的检测概率。
10.根据权利要求9所述的多径环境下机载相控阵雷达信号检测方法,其特征在于,所述检测概率的计算公式为:
其中,Pd是检测概率,S′是阈值,Xi是第i个快时间数据,T是标称因子,N是参考单元数,μ为背景噪声功率,n′=2N,z是计算检测概率的中间变量,Ps是直达和多径回波功率大小,A是信号幅度大小,ρ是综合多径反射系数,QM(·)是Marcum的Q函数,I1(·)是修正的一阶贝塞尔函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110144082.7A CN112986939B (zh) | 2021-02-02 | 2021-02-02 | 一种多径环境下机载相控阵雷达信号检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110144082.7A CN112986939B (zh) | 2021-02-02 | 2021-02-02 | 一种多径环境下机载相控阵雷达信号检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112986939A CN112986939A (zh) | 2021-06-18 |
CN112986939B true CN112986939B (zh) | 2023-12-22 |
Family
ID=76346207
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110144082.7A Active CN112986939B (zh) | 2021-02-02 | 2021-02-02 | 一种多径环境下机载相控阵雷达信号检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112986939B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2452415A (en) * | 2007-08-31 | 2009-03-04 | Raymarine Uk Ltd | Digital radar or sonar apparatus |
CN107607938A (zh) * | 2017-08-28 | 2018-01-19 | 电子科技大学 | 一种频率分集阵列雷达低空目标检测方法 |
CN108445461A (zh) * | 2018-01-29 | 2018-08-24 | 中国人民解放军国防科技大学 | 一种多径条件下雷达目标检测方法 |
CN110865362A (zh) * | 2019-11-29 | 2020-03-06 | 桂林电子科技大学 | 一种基于fda-mimo雷达的低慢小目标检测方法 |
CN111751800A (zh) * | 2020-06-04 | 2020-10-09 | 南昌大学 | 一种频控阵雷达角度-距离参数解耦合方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10802130B2 (en) * | 2017-11-28 | 2020-10-13 | Viettel Group | Marine target detection in cluttered environments |
-
2021
- 2021-02-02 CN CN202110144082.7A patent/CN112986939B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2452415A (en) * | 2007-08-31 | 2009-03-04 | Raymarine Uk Ltd | Digital radar or sonar apparatus |
CN107607938A (zh) * | 2017-08-28 | 2018-01-19 | 电子科技大学 | 一种频率分集阵列雷达低空目标检测方法 |
CN108445461A (zh) * | 2018-01-29 | 2018-08-24 | 中国人民解放军国防科技大学 | 一种多径条件下雷达目标检测方法 |
CN110865362A (zh) * | 2019-11-29 | 2020-03-06 | 桂林电子科技大学 | 一种基于fda-mimo雷达的低慢小目标检测方法 |
CN111751800A (zh) * | 2020-06-04 | 2020-10-09 | 南昌大学 | 一种频控阵雷达角度-距离参数解耦合方法 |
Non-Patent Citations (2)
Title |
---|
低空多径环境下雷达目标检测性能研究;周豪;胡国平;匡旭斌;师俊朋;;现代雷达(第02期);全文 * |
脉冲雷达回波信号检测性能比较分析;洪兴勇;洪一;李文谨;胡国华;;电路与***学报(第02期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112986939A (zh) | 2021-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2932296B1 (en) | Subsurface imaging radar | |
Sjögren et al. | Suppression of clutter in multichannel SAR GMTI | |
US20030142011A1 (en) | Surface wave radar | |
US6580392B2 (en) | Digital beamforming for passive detection of target using reflected jamming echoes | |
Vespe et al. | Radar target classification using multiple perspectives | |
Hu et al. | Optimal signal processing in ground-based forward scatter micro radars | |
CN110658514B (zh) | 一种水下静态目标的分类识别方法 | |
Li et al. | Radar and inverse scattering | |
CN103197295B (zh) | 利用先验信息的宽带频率捷变角度超分辨方法 | |
CN109581366B (zh) | 一种基于目标导向矢量失配的离散旁瓣杂波识别方法 | |
CN112986939B (zh) | 一种多径环境下机载相控阵雷达信号检测方法 | |
CN112684425B (zh) | 一种恒虚警检测后的目标二次甄别方法 | |
CN115616502A (zh) | 无人飞行器机载雷达目标检测的杂波抑制方法 | |
Cao et al. | Target detection for low angle radar based on multi-frequency order-statistics | |
CN113093174A (zh) | 基于phd滤波雷达起伏微弱多目标的检测前跟踪方法 | |
Vu et al. | Experimental results on moving target detection by focusing in UWB low frequency SAR | |
CN112986899B (zh) | 一种机载mimo雷达在多径环境下的doa估计方法 | |
RU2773818C1 (ru) | Способ оценки эффективной площади рассеяния аэродинамической цели | |
Choi et al. | Use of Stepped Carrier Frequency in Fast Scan Mode to Detect Small Targets on the Sea Surface | |
Onunka et al. | Modelling Radar Signal Error Performance under Atmospheric Refraction and Clutter Attenuation. | |
Belhi et al. | Nouvelles méthodes pour l’estimation des paramètres des distributions gaussiennes composées. | |
Antipov et al. | Estimation of a constant false alarm rate processing loss for a high-resolution maritime radar system | |
Marcus | Balancing Antenna Performance vs. Radar Cross Section for a Passive Radar-Detecting Sensor on an Aircraft | |
Hu et al. | Incoherent Integration Detection Method of Airborne Phased Array Radar in a Multipath Environment | |
Elhoshy | Signal processing techniques for modern radar systems |
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 |