CN106597411B - 雷达信号处理方法 - Google Patents
雷达信号处理方法 Download PDFInfo
- Publication number
- CN106597411B CN106597411B CN201611261993.3A CN201611261993A CN106597411B CN 106597411 B CN106597411 B CN 106597411B CN 201611261993 A CN201611261993 A CN 201611261993A CN 106597411 B CN106597411 B CN 106597411B
- Authority
- CN
- China
- Prior art keywords
- signal
- echo
- target
- airbound target
- echo 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
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
本发明公开了一种雷达信号处理方法,属于信号处理领域。该方法包括接收回波信号,所述回波信号是发射信号受到飞行目标反射后形成的信号;将所述回波信号转换为基带信号;对所述基带信号进行预定信号处理操作,得到与所述飞行目标对应的目标信号;所述预定信号处理操作至少包括快速傅里叶变换FFT处理和二维恒虚警率CFAR检测;根据所述目标信号确定所述飞行目标的运动轨迹和/或运动趋势;解决了雷达在监测飞行高度低、体积小、速度慢的飞行目标时漏检率和虚警概率高的问题,达到了降低漏检率和虚警概率,提高对飞行目标的定位及跟踪精度的效果。
Description
技术领域
本发明实施例涉及信号处理领域,特别涉及一种雷达信号处理方法。
背景技术
随着诸如无人机、三角翼、滑翔伞、热气球、航模等飞行目标的普及和应用,对该类飞行目标的探测和跟踪具有迫切需求,
相关技术中,雷达向外发射发射信号,再接收由飞行目标反射回来的回波信号,利用脉冲压缩技术和相参积累等技术对回波信号进行处理,得到飞行目标的信息,再根据飞行目标的信息进行预警。
然而,该类飞行目标具有飞行高度低、体积小、速度慢的特点,容易隐匿于地杂波中,且RCS(Radar Cross-Section,雷达散射截面积)小也令该类飞行目标不易于被雷达发现,导致监测时漏检率和虚警概率高。
发明内容
为了解决现有技术的问题,本发明实施例提供了一种雷达信号处理方法及装置。该技术方案如下:
第一方面,提供了一种雷达信号处理方法,该方法包括:
接收回波信号,回波信号是发射信号受到飞行目标反射后形成的信号;
将回波信号转换为基带信号;
对基带信号进行预定信号处理操作,得到与飞行目标对应的目标信号;预定信号处理操作至少包括快速傅里叶变换FFT处理和二维恒虚警率CFAR检测;
根据目标信号确定飞行目标的运动轨迹和/或运动趋势。
可选的,将回波信号转换为基带信号,包括:
对回波信号进行模拟数字AD采样,得到两个正交的数字信号;
对两个正交的数字信号进行数字下变频处理,得到两个正交的基带信号。
可选的,对基带信号进行预定信号处理操作,得到与飞行目标对应的目标信号,包括:
对基带信号进行正交解调处理,得到线性调频信号;
对线性调频信号进行脉冲压缩处理,得到第一信号;
对第一信号进行FFT处理以及相参积累,得到第二信号;
在时域和频域上同时对第二信号进行二维CFAR检测,得到与飞行目标对应的目标信号。
可选的,在时域和频域上同时对第二信号进行二维CFAR检测,得到与飞行目标对应的目标信号,包括:
检测第二信号的检测单元M(k,l)是否满足M(k,l)>Ka×Y(k,l)且M(k,l)>Kb;
若检测单元M(k,l)满足M(k,l)>Ka×Y(k,l)且M(k,l)>Kb,则确定检测单元M(k,l)对应的第二信号的信号幅度过门限;
对信号幅度超过门限的第二信号,依据预设的虚警率进行门限调整,得到目标信号;
其中,M(k,l)和Y(k,l)按照如下公式计算:
Y(k,l)为二维CFAR滑动参考窗内的信号平均模,k表示距离门号,l表示滤波器号,Ka表示固定门限,Kb表示浮动门限,Rwidth表示二维CFAR滑窗的距离向长度;Fwidth表示二维CFAR滑窗的频率向长度;Rwb表示二维CFAR滑窗保护区域的距离向长度;Fwb表示二维CFAR滑窗保护区域的频率向长度。
可选的,当目标信号的数量为N个时,N≥2,方法还包括:
对N个目标信号进行目标跨接处理,得到实际目标信号。
可选的,对N个目标信号进行目标跨接处理,得到实际目标信号,包括:
针对第i个目标信号,检测Nr(i)-Nr(i+1)是否小于rvar且Nf(i)-Nf(i+1)是否小于fvar;
若Nr(i)-Nr(i+1)小于rvar且Nf(i)-Nf(i+1)小于fvar,则检测第i个目标信号对应的飞行目标的功率是否小于第i+1个目标信号对应的飞行目标的功率;
若第i个目标信号对应的飞行目标的功率小于第i+1个目标信号对应的飞行目标的功率,则删除第i个目标信号;
在i<N-1时,令i=i+1,重复执行针对第i个目标信号,检测Nr(i)-Nr(i+1)是否小于rvar且Nf(i)-Nf(i+1)是否小于fvar的步骤;
其中,Nr表示目标的距离门号,Nf表示目标的频率门号,rvar和fvar表示可调系统参数。
可选的,该方法还包括:
若第i个目标信号对应的飞行目标的功率不小于第i+1个目标信号对应的飞行目标的功率,则删除第i+1个目标信号。
可选的,当飞行目标的数量至少为两个时,根据目标信号确定飞行目标的运动轨迹和/或运动趋势,包括:
检测目标信号对应的点迹数据是否与预定波门关联,预定波门是根据飞行目标已经出现过的位置预测得到的;
若目标信号对应的点迹数据与预定波门关联,则对目标信号进行航迹维持得到第一航迹,对第一航迹进行滤波更新;
若目标信号对应的点迹数据未与预定波门关联,则根据目标信号生成第二航迹;
将经过滤波更新后的第一航迹和第二航迹按时间顺序排序;
将在预定时间内未出现更新的航迹删除,将属于同一飞行目标的航迹合并,得到第三航迹;
根据第三航迹得到飞行目标的运动轨迹和/或运动趋势。
可选的,该方法还包括:
将飞行目标的运动轨迹和/或运动趋势发送至监控设备,监控设备用于显示飞行目标的运动轨迹和/或运动趋势。
第二方面,提供了一种雷达信号处理装置,该装置包括:
接收模块,用于接收回波信号,所述回波信号是发射信号受到飞行目标反射后形成的信号;
转换模块,用于将所述回波信号转换为基带信号;
信号处理模块,用于对所述基带信号进行预定信号处理操作,得到与所述飞行目标对应的目标信号;所述预定信号处理操作至少包括快速傅里叶变换FFT处理和二维恒虚警率CFAR检测;
确定模块,用于根据所述目标信号确定所述飞行目标的运动轨迹和/或运动趋势。
可选的,转换模块,包括:
采样单元,用于对回波信号进行模拟数字AD采样,得到两个正交的数字信号;
数字下变频处理单元,用于对两个正交的数字信号进行数字下变频处理,得到两个正交的基带信号。
可选的,处理模块,包括:
第一处理单元,用于对基带信号进行正交解调处理,得到线性调频信号;
第二处理单元,用于对线性调频信号进行脉冲压缩处理,得到第一信号;
第三处理单元,用于对第一信号进行FFT处理以及相参积累,得到第二信号;
第四处理单元,用于在时域和频域上同时对第二信号进行二维CFAR检测,得到与飞行目标对应的目标信号。
可选的,第四处理单元,具体用于:
检测第二信号的检测单元M(k,l)是否满足M(k,l)>Ka×Y(k,l)且M(k,l)>Kb;
若检测单元M(k,l)满足M(k,l)>Ka×Y(k,l)且M(k,l)>Kb,则确定检测单元M(k,l)对应的第二信号的信号幅度过门限;
对信号幅度超过门限的第二信号,依据预设的虚警率进行门限调整,得到目标信号;
其中,M(k,l)和Y(k,l)按照如下公式计算:
Y(k,l)为二维CFAR滑动参考窗内的信号平均模,k表示距离门号,l表示滤波器号,Ka表示固定门限,Kb表示浮动门限,Rwidth表示二维CFAR滑窗的距离向长度;Fwidth表示二维CFAR滑窗的频率向长度;Rwb表示二维CFAR滑窗保护区域的距离向长度;Fwb表示二维CFAR滑窗保护区域的频率向长度。
可选的,当目标信号的数量为N个时,N≥2,装置还包括:
跨接处理模块,用于对N个目标信号进行目标跨接处理,得到实际目标信号。
可选的,跨接处理模块,具体用于:
针对第i个目标信号,检测Nr(i)-Nr(i+1)是否小于rvar且Nf(i)-Nf(i+1)是否小于fvar;
若Nr(i)-Nr(i+1)小于rvar且Nf(i)-Nf(i+1)小于fvar,则检测第i个目标信号对应的飞行目标的功率是否小于第i+1个目标信号对应的飞行目标的功率;
若第i个目标信号对应的飞行目标的功率小于第i+1个目标信号对应的飞行目标的功率,则删除第i个目标信号;
在i<N-1时,令i=i+1,重复执行针对第i个目标信号,检测Nr(i)-Nr(i+1)是否小于rvar且Nf(i)-Nf(i+1)是否小于fvar的步骤;
其中,Nr表示目标的距离门号,Nf表示目标的频率门号,rvar和fvar表示可调系统参数。
可选的,跨接处理模块,还用于:
若第i个目标信号对应的飞行目标的功率不小于第i+1个目标信号对应的飞行目标的功率,则删除第i+1个目标信号。
可选的,当飞行目标的数量至少为两个时,确定模块,具体用于:
检测目标信号对应的点迹数据是否与预定波门关联,预定波门是根据飞行目标已经出现过的位置预测得到的;
若目标信号对应的点迹数据与预定波门关联,则对目标信号进行航迹维持得到第一航迹,对第一航迹进行滤波更新;
若目标信号对应的点迹数据未与预定波门关联,则根据目标信号生成第二航迹;
将经过滤波更新后的第一航迹和第二航迹按时间顺序排序;
将在预定时间内未出现更新的航迹删除,将属于同一飞行目标的航迹合并,得到第三航迹;
根据第三航迹得到飞行目标的运动轨迹和/或运动趋势。
可选的,装置还包括:
发送模块,用于飞行目标的运动轨迹和/或运动趋势发送至监控设备,监控设备用于显示飞行目标的运动轨迹和/或运动趋势。
本发明实施例提供的技术方案带来的有益效果是:
本发明实施例提供的雷达信号处理装置,通过接收回波信号,将回波信号转换为基带信号,对基带信号至少进行FFT处理和二维CFAR处理,解决了雷达在监测飞行高度低、体积小、速度慢的飞行目标时漏检率和虚警概率高的问题,达到了降低漏检率和虚警概率,提高对飞行目标的定位及跟踪精度的效果。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是根据一示例性实施例示出的一种雷达信号处理方法的流程图;
图2是根据另一示例性实施例示出的一种雷达信号处理方法的流程图;
图3是根据另一示例性实施例示出的一种二维CFAR滑窗的示意图;
图4是根据另一示例性实施例示出的一种雷达信号处理方法的实施示意图;
图5是根据另一示例性实施例示出的一种雷达信号处理装置的框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
请参考图1,其示出了本发明一个实施例提供的雷达信号处理方法的流程图。该在雷达信号处理方法适用于对地监视雷达和/或低空监视雷达的处理板中。如图1所示,该雷达信号处理方法可以包括以下步骤:
步骤101,接收回波信号,回波信号是发射信号受到飞行目标反射后形成的信号。
可选的,飞行目标的数量为一个,或者,飞行目标的数量为至少两个。
雷达通过天线向外发射发射信号,发射信号受到飞行目标反射后形成回波信号,通过天线接收回波信号。
步骤102,将回波信号转换为基带信号。
步骤103,对基带信号进行预定信号处理操作,得到与飞行目标对应的目标信号。
预定信号处理操作至少包括FFT(Fast Fourier Transformation,快速傅里叶变换)处理和二维CFAR(Constant false alarm rate,恒虚警率)检测。
由于监测环境复杂,在接收回波信号时还会接收到杂波,需要将目标信号从杂波中分离出来,因此需要对基带信号进行预订信号处理操作。
步骤104,根据目标信号确定飞行目标的运动轨迹和/或运动趋势。
可选的,根据目标信号确定出飞行目标的方位角、距离信息、高度信息。
可选的,根据目标信号确定出飞行目标的点迹数据,根据点迹数据确定出飞行目标的运动轨迹和/或运动趋势。
综上所述,本发明实施例提供的雷达信号处理方法,通过接收回波信号,将回波信号转换为基带信号,对基带信号至少进行FFT处理和二维CFAR处理,解决了雷达在监测飞行高度低、体积小、速度慢的飞行目标时漏检率和虚警概率高的问题,达到了降低漏检率和虚警概率,提高对飞行目标的定位及跟踪精度的效果。
请参考图2,其示出了本发明另一个实施例提供的雷达信号处理方法的流程图。该在雷达信号处理方法适用于对地监视雷达和/或低空监视雷达的处理板中。如图2所示,该雷达信号处理方法可以包括以下步骤:
步骤201,接收回波信号。
回波信号是发射信号受到飞行目标反射后形成的信号。
可选的,回波信号为中频信号。
可选的,回波信号通过公式表示,式中为矩阵函数,为一个时间宽度为T的周期性方波信号,A表示信号幅度,f0表示发射信号的中心频率,fd表示目标多普勒频率。
步骤202,对回波信号进行AD(analog to digital,模拟数字)采样,得到两个正交的数字信号。
通过AD采样将回波信号变换到数字域。
步骤203,对两个正交的数字信号进行数字下变频处理,得到两个正交的基带信号。
可选的,通过FPGA(Field Programmable Gate Array,现场可编程门阵列)对两个正交的数字信号进行数字下变频处理。
可选的,按照30倍对两个正交的数字信号进行抽取,将240M采样信号降到8M,使得基带信号的速率降到DSP芯片(Digital Signal Processing,数字信号处理)可以处理的范围内,大约10MHz以内。
步骤204,对基带信号进行正交解调处理,得到线性调频信号。
线性调频信号相比回波信号去除了发射信号的中心频率f0,保留了目标多普勒频率fd,线性调频信号的表达式如下:
对线性调频信号S(t)进行“0”中频解调后,其数学表达式如下:
其中,Km表示线性调频斜率。
步骤205,对线性调频信号进行脉冲压缩处理,得到第一信号。
由于解调后的线性调频信号具有较高的脉冲压缩信噪比增益,因此,系统对线性调频信号进行脉冲压缩处理,得到第一信号。
由于线性调频信号的宽带宽积为49(脉宽×带宽=7MHz×7us),远远大于1,因此带来了约16dB(16≈10×log10(49))的脉冲压缩信噪比增益,有效提高目标散射体回波的信噪比,方便后续检测。另外,由于脉冲压缩处理将回波信号在距离维度上大幅收窄,从而将相距较近的两个飞行目标产生的回波信号分离开来,提高后的距离分辨率为21.43m(光速/(2×带宽)=3×108/(2×7×10(-6)))。
步骤206,对第一信号进行FFT处理以及相参积累,得到第二信号。
具体地,对第一信号进行FFT处理,将第一信号离散化,得到信号s(n),s(n)的信号表达式如下所示:
其中,M=T/△,△表示ADC(Analog-to-Digital Converter,模数转换器)的采样间隔。s(n)的表达式中加0.5是因为线性调频信号是一个对称信号,其对称中心子在两个采样之间。
再将信号s(n)与匹配滤波器在频域进行相乘,完成相参积累。
匹配滤波器的数学表达式为线性调频信号的复共轭函数:
h(n)=exp[j2πKm(-M/2+n+0.5)2△2],n=0,1,2,...M-1;
信号s(n)与匹配滤波器在频域进行相乘的表达式如下:
S(f)=FFT[s(n)×h(n)],n=0,1,2,...N-1。
通过相参积累可以提高回波信号的信噪比,依据天线扫描速度、天线波束宽度、PRF(pulse recurrence frequency,脉冲重复频率)、速度分辨率等参数指标设计出积累脉冲个数为128个,因此会带来约21dB的积累信噪比改善,从而进一步提高了回波信号的信噪比。另外,通过FFT处理可将不同速度的飞行目标在多普勒维度上区分开来,特别是将低空飞行目标与固定的地面、山脉、建筑物、树木等干扰物体引起的地杂波区分开来。
步骤207,在时域和频域上同时对第二信号进行CFAR检测,得到与飞行目标对应的目标信号。
具体地,该步骤由如下几个步骤实现:
步骤1,检测第二信号的检测单元M(k,l)是否满足M(k,l)>Ka×Y(k,l)且M(k,l)>Kb。
步骤2,若所述检测单元M(k,l)满足M(k,l)>Ka×Y(k,l)且M(k,l)>Kb,则确定检测单元M(k,l)对应的第二信号的信号幅度过门限。
其中,M(k,l)和Y(k,l)按照如下公式计算:
Y(k,l)为二维CFAR滑动参考窗内的信号平均模,k表示距离门号,l表示滤波器号,Ka表示固定门限,Kb表示浮动门限,Rwidth表示二维CFAR滑窗的距离向长度;Fwidth表示二维CFAR滑窗的频率向长度;Rwb表示二维CFAR滑窗保护区域的距离向长度;Fwb表示二维CFAR滑窗保护区域的频率向长度。图3示例性地示出了二维CFAR滑窗的示意图。
步骤3,对信号幅度超过门限的第二信号,依据预设的虚警率进行门限调整,得到目标信号。
可选的,预设的虚警率为85%。
由于经过CFAR处理得到的第二信号中包括真实飞行目标的第二信号和虚假飞行目标的第二信号,调整门限后,将虚假飞行目标的第二信号过滤,只留下了真实飞行目标的第二信号,也即目标信号。
虚假飞行目标是由杂波、干扰引起的,真实飞行目标是飞行高度低、体积小、速度慢的飞行目标。
当目标信号的数量为N个,N≥2时,执行步骤208。
步骤208,对N个目标信号进行目标跨接处理,得到实际目标信号。
由于当回波信号中存在强目标时,在二维CFAR处理时,可能出现目标分类,将一个飞行目标作为多个飞行目标输出,因此需要对目标信号进行目标跨接处理
该步骤具体实现方式如下:
针对第i个目标信号,检测Nr(i)-Nr(i+1)是否小于rvar且Nf(i)-Nf(i+1)是否小于fvar。
其中,Nr表示目标的距离门号,Nf表示目标的频率门号,rvar和fvar表示可调系统参数。
若Nr(i)-Nr(i+1)小于rvar且Nf(i)-Nf(i+1)小于fvar,则检测第i个目标信号对应的飞行目标的功率是否小于第i+1个目标信号对应的飞行目标的功率。
若第i个目标信号对应的飞行目标的功率小于第i+1个目标信号对应的飞行目标的功率,则删除第i个目标信号。
若第i个目标信号对应的飞行目标的功率不小于第i+1个目标信号对应的飞行目标的功率,则删除第i+1个目标信号。
在i<N-1时,令i=i+1,重复执行针对第i个目标信号,检测Nr(i)-Nr(i+1)是否小于rvar且Nf(i)-Nf(i+1)是否小于fvar的步骤。
步骤209,根据目标信号确定飞行目标的运动轨迹和/或运动趋势。
可选的,根据目标信号确定出飞行目标的方位角、距离信息、高度信息。
为了提高飞行目标的测距和测速精度,还需要对飞行目标的距离和速度坐标进行求质心处理,计算公式如下:
其中,N表示离散化的坐标;Ncent为求完质心的浮点坐标;P(·)表示对应点功率。
通过求质心处理检测出飞行目标的实际的点迹数据,并提取飞行目标的点迹的距离门信息、方位信息、俯仰角度及频率门信息。
飞行目标的距离信息根据距离门信息计算得到,计算公式为距离=CFAR距离门/基带采样率×150。
方位信息也即飞行目标的目标质心所处的雷达实时扫描中心角度;俯仰角度参考俯仰向电扫描交底以及和差通道计算的目标俯仰角。
飞行目标的速度根据动目标的多普勒频移计算出目标的径向速度,再根据雷达扫描俯仰角即可解算出飞行目标的真实运动速度,计算公式为:速度=CFAR频域速度门×(PRF/128)*波长/2。
可选的,根据目标信号确定出飞行目标的点迹数据,根据点迹数据确定出飞行目标的运动轨迹和/或运动趋势。
采用TWS(track-while-scan,边扫描边跟踪)的方法实现对多个飞行目标的跟踪。
该步骤具体由如下方式实现:
如图4所示,通过测量数据输入模块31对目标信号进行预处理,得到TES所需的信息。在同一天线帧内,波束扫过时,同一个飞行目标会在多个雷达帧内有回波信号返回,点迹预处理模块32将这些雷达帧内的同一个飞行目标的书籍进行相关处理,形成综合点迹发送给航迹关联模块33。
检测目标信号是否与预定波门关联,预定波门是根据飞行目标已经出现过的位置预测得到的。
具体地,根据飞行目标已经出现过得位置进行目标位置预测34,生成预定波门35,航迹关联模块33检测目标信号是否与预定波门35关联。
若目标信号与预定波门关联,则对目标信号进行航迹维持得到第一航迹,对第一航迹进行滤波更新。
具体地,将与预定波门关联的目标信号放入航迹维持模块41,进行机动检测36得到第一航迹,并对第一航迹进行kalman(卡尔曼)自适应滤波42。
若目标信号未与预定波门关联,则根据目标信号生成第二航迹。
若目标信号未与预定波门关联,则将目标信号放入航迹起始模块37,由航迹起始模块37生产临时的第二航迹。
将经过滤波更新后的第一航迹和第二航迹按时间顺序排序。
第一航迹和第二航迹一起进入航迹管理模块38,由航迹管理模块38对第一航迹和第二航迹按时间顺序进行排序。
将在预定时间内未出现更新的航迹删除,将属于同一飞行目标的航迹合并,得到第三航迹。
可选的,如果在一定距离窗、角度窗、速度窗内则认为是同一飞行目标的产生航迹,并进行航迹合并。
根据所述第三航迹得到所述飞行目标的运动轨迹和/或运动趋势。
航迹管理模块38将第三航迹发送至目标状态输出模块39和/或目标位置预测模块40,由目标状态输出模块39根据第三航迹输出飞行目标的运动轨迹,和/或,由目标位置预测模块40根据第三航迹预测飞行目标的运动趋势并输出。
步骤210,将飞行目标的运动轨迹和/或运动趋势发送至监控设备,监控设备用于显示飞行目标的运动轨迹和/或运动趋势。
可选的,还通过处理板将雷达探测图像、目标点迹数据、航迹数据发送至监控设备,监控设备显示雷达探测图像、目标点迹数据、航迹数据、飞行目标的运动轨迹和/或运动趋势。
综上所述,本发明实施例提供的雷达信号处理方法,通过接收回波信号,将回波信号转换为基带信号,对基带信号至少进行FFT处理和二维CFAR处理,解决了雷达在监测飞行高度低、体积小、速度慢的飞行目标时漏检率和虚警概率高的问题,达到了降低漏检率和虚警概率,提高对飞行目标的定位及跟踪精度的效果。
下述为本发明装置实施例,可以用于执行本发明方法实施例。对于本发明装置实施例中未披露的细节,请参照本发明方法实施例。
请参照图5,其示出了本发明一个实施例提供的雷达信号处理装置的结构方框图。该雷达信号处理装置可以通过软件、硬件或者两者的结合实现成为上述可提供雷达信号处理方法的低空监视雷达和/或对地监视雷达的全部或者一部分。该装置包括:
接收模块510,用于接收回波信号,所述回波信号是发射信号受到飞行目标反射后形成的信号。
转换模块520,用于将所述回波信号转换为基带信号。
信号处理模块530,用于对所述基带信号进行预定信号处理操作,得到与所述飞行目标对应的目标信号;所述预定信号处理操作至少包括快速傅里叶变换FFT处理和二维恒虚警率CFAR检测。
确定模块540,用于根据所述目标信号确定所述飞行目标的运动轨迹和/或运动趋势。
综上所述,本发明实施例提供的雷达信号处理装置,通过接收回波信号,将回波信号转换为基带信号,对基带信号至少进行FFT处理和二维CFAR处理,解决了雷达在监测飞行高度低、体积小、速度慢的飞行目标时漏检率和虚警概率高的问题,达到了降低漏检率和虚警概率,提高对飞行目标的定位及跟踪精度的效果。
可选的,转换模块,包括:
采样单元,用于对回波信号进行模拟数字AD采样,得到两个正交的数字信号;
数字下变频处理单元,用于对两个正交的数字信号进行数字下变频处理,得到两个正交的基带信号。
可选的,处理模块,包括:
第一处理单元,用于对基带信号进行正交解调处理,得到线性调频信号。
第二处理单元,用于对线性调频信号进行脉冲压缩处理,得到第一信号。
第三处理单元,用于对第一信号进行FFT处理以及相参积累,得到第二信号。
第四处理单元,用于在时域和频域上同时对第二信号进行二维CFAR检测,得到与飞行目标对应的目标信号。
可选的,第四处理单元,具体用于:
检测第二信号的检测单元M(k,l)是否满足M(k,l)>Ka×Y(k,l)且M(k,l)>Kb;
若检测单元M(k,l)满足M(k,l)>Ka×Y(k,l)且M(k,l)>Kb,则确定检测单元M(k,l)对应的第二信号的信号幅度过门限;
对信号幅度超过门限的第二信号,依据预设的虚警率进行门限调整,得到目标信号;
其中,M(k,l)和Y(k,l)按照如下公式计算:
Y(k,l)为二维CFAR滑动参考窗内的信号平均模,k表示距离门号,l表示滤波器号,Ka表示固定门限,Kb表示浮动门限,Rwidth表示二维CFAR滑窗的距离向长度;Fwidth表示二维CFAR滑窗的频率向长度;Rwb表示二维CFAR滑窗保护区域的距离向长度;Fwb表示二维CFAR滑窗保护区域的频率向长度。
可选的,当目标信号的数量为N个时,N≥2,装置还包括:
跨接处理模块,用于对N个目标信号进行目标跨接处理,得到实际目标信号。
可选的,跨接处理模块,具体用于:
针对第i个目标信号,检测Nr(i)-Nr(i+1)是否小于rvar且Nf(i)-Nf(i+1)是否小于fvar;
若Nr(i)-Nr(i+1)小于rvar且Nf(i)-Nf(i+1)小于fvar,则检测第i个目标信号对应的飞行目标的功率是否小于第i+1个目标信号对应的飞行目标的功率;
若第i个目标信号对应的飞行目标的功率小于第i+1个目标信号对应的飞行目标的功率,则删除第i个目标信号;
在i<N-1时,令i=i+1,重复执行针对第i个目标信号,检测Nr(i)-Nr(i+1)是否小于rvar且Nf(i)-Nf(i+1)是否小于fvar的步骤;
其中,Nr表示目标的距离门号,Nf表示目标的频率门号,rvar和fvar表示可调系统参数。
可选的,跨接处理模块,还用于:
若第i个目标信号对应的飞行目标的功率不小于第i+1个目标信号对应的飞行目标的功率,则删除第i+1个目标信号。
可选的,当飞行目标的数量至少为两个时,确定模块,具体用于:
检测目标信号对应的点迹数据是否与预定波门关联,预定波门是根据飞行目标已经出现过的位置预测得到的;
若目标信号对应的点迹数据与预定波门关联,则对目标信号进行航迹维持得到第一航迹,对第一航迹进行滤波更新;
若目标信号对应的点迹数据未与预定波门关联,则根据目标信号生成第二航迹;
将经过滤波更新后的第一航迹和第二航迹按时间顺序排序;
将在预定时间内未出现更新的航迹删除,将属于同一飞行目标的航迹合并,得到第三航迹;
根据第三航迹得到飞行目标的运动轨迹和/或运动趋势。
可选的,装置还包括:
发送模块,用于飞行目标的运动轨迹和/或运动趋势发送至监控设备,监控设备用于显示飞行目标的运动轨迹和/或运动趋势。
需要说明的是,采用本发明实施例提供的雷达信号处理方法的低空监视雷达和/或对地监视雷达,可以达到如下性能指标:
1、距离量程:最高32km;
2、作用距离(Pc=0.8,平均虚警率不多于1个航迹/min):
对RCS≥5m2的空中目标(如:固定翼无人机、直升机):29km;
对RCS≥0.1m2的空中小目标(如:小型无人机):11km;
3、最小可检测目标径向速度(绝对值):5m/s;
4、不模糊测速范围:±200m/s;
5、雷达探测盲区:≤600m;
6、距离分辨率:≤30m;
7、测距精度:≤50m;
8、测角精度:≤0.5°;
9、测高精度:在10km的距离上,优于300m;
10、目标处理容量:多目标跟踪数不小于64个。
需要说明的是:上述实施例提供的雷达信号处理装置在执行雷达信号处理方法时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模块完成,即将设备的内部结构划分成不同的功能模块,以完成以上描述的全部或者部分功能。另外,上述实施例提供的雷达信号处理装置与雷达信号处理方法实施例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种雷达信号处理方法,其特征在于,所述方法包括:
接收回波信号,所述回波信号是发射信号受到飞行目标反射后形成的信号;
将所述回波信号转换为基带信号;
对所述基带信号进行预定信号处理操作,得到与所述飞行目标对应的目标信号;所述预定信号处理操作至少包括快速傅里叶变换FFT处理和二维恒虚警率CFAR检测;
当所述飞行目标的数量至少为两个时,检测所述目标信号对应的点迹数据是否与预定波门关联,所述预定波门是根据所述飞行目标已经出现过的位置预测得到的;
若所述目标信号对应的点迹数据与所述预定波门关联,则对所述目标信号进行航迹维持得到第一航迹,对所述第一航迹进行滤波更新;
若所述目标信号对应的点迹数据未与所述预定波门关联,则根据所述目标信号生成第二航迹;
将经过滤波更新后的所述第一航迹和所述第二航迹按时间顺序排序;
将在预定时间内未出现更新的航迹删除,将属于同一飞行目标的航迹合并,得到第三航迹;
根据所述第三航迹得到所述飞行目标的运动轨迹和/或运动趋势。
2.根据权利要求1所述的方法,其特征在于,所述将所述回波信号转换为基带信号,包括:
对所述回波信号进行模拟数字AD采样,得到两个正交的数字信号;
对两个正交的所述数字信号进行数字下变频处理,得到两个正交的所述基带信号。
3.根据权利要求1所述的方法,其特征在于,所述对所述基带信号进行预定信号处理操作,得到与所述飞行目标对应的目标信号,包括:
对所述基带信号进行正交解调处理,得到线性调频信号;
对所述线性调频信号进行脉冲压缩处理,得到第一信号;
对所述第一信号进行所述FFT处理以及相参积累,得到第二信号;
在时域和频域上同时对所述第二信号进行所述二维CFAR检测,得到与所述飞行目标对应的目标信号。
4.根据权利要求3所述的方法,其特征在于,所述在时域和频域上同时对所述第二信号进行所述二维CFAR检测,得到与所述飞行目标对应的目标信号,包括:
检测所述第二信号的检测单元M(k,l)是否满足M(k,l)>Ka×Y(k,l)且M(k,l)>Kb;
若所述检测单元M(k,l)满足M(k,l)>Ka×Y(k,l)且M(k,l)>Kb,则确定所述检测单元M(k,l)对应的所述第二信号的信号幅度过门限;
对信号幅度超过所述门限的第二信号,依据预设的虚警率进行门限调整,得到所述目标信号;
其中,所述M(k,l)和所述Y(k,l)按照如下公式计算:
所述Y(k,l)为二维CFAR滑动参考窗内的信号平均模,k表示距离门号,l表示滤波器号,Ka表示固定门限,Kb表示浮动门限,Rwidth表示二维CFAR滑窗的距离向长度;Fwidth表示二维CFAR滑窗的频率向长度;Rwb表示二维CFAR滑窗保护区域的距离向长度;Fwb表示二维CFAR滑窗保护区域的频率向长度。
5.根据权利要求4所述的方法,其特征在于,当所述目标信号的数量为N个时,N≥2,所述方法还包括:
对N个所述目标信号进行目标跨接处理,得到实际目标信号。
6.根据权利要求5所述的方法,其特征在于,所述对N个所述目标信号进行目标跨接处理,得到实际目标信号,包括:
针对第i个目标信号,检测Nr(i)-Nr(i+1)是否小于rvar且Nf(i)-Nf(i+1)是否小于fvar;
若所述Nr(i)-Nr(i+1)小于所述rvar且所述Nf(i)-Nf(i+1)小于所述fvar,则检测所述第i个目标信号对应的飞行目标的功率是否小于所述第i+1个目标信号对应的飞行目标的功率;
若所述第i个目标信号对应的飞行目标的功率小于所述第i+1个目标信号对应的飞行目标的功率,则删除所述第i个目标信号;
在i<N-1时,令i=i+1,重复执行所述针对第i个目标信号,检测Nr(i)-Nr(i+1)是否小于rvar且Nf(i)-Nf(i+1)是否小于fvar的步骤;
其中,Nr表示目标的距离门号,Nf表示目标的频率门号,rvar和fvar表示可调系统参数。
7.根据权利要求6所述的方法,其特征在于,所述方法还包括:
若所述第i个目标信号对应的飞行目标的功率不小于所述第i+1个目标信号对应的飞行目标的功率,则删除所述第i+1个目标信号。
8.根据权利要求1至7任一所述的方法,其特征在于,所述方法还包括:
将所述飞行目标的运动轨迹和/或运动趋势发送至监控设备,所述监控设备用于显示所述飞行目标的运动轨迹和/或运动趋势。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611261993.3A CN106597411B (zh) | 2016-12-30 | 2016-12-30 | 雷达信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611261993.3A CN106597411B (zh) | 2016-12-30 | 2016-12-30 | 雷达信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106597411A CN106597411A (zh) | 2017-04-26 |
CN106597411B true CN106597411B (zh) | 2019-06-04 |
Family
ID=58582110
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611261993.3A Active CN106597411B (zh) | 2016-12-30 | 2016-12-30 | 雷达信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106597411B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3486680B1 (en) * | 2017-11-21 | 2024-01-03 | Rohde & Schwarz GmbH & Co. KG | Testing system and method for testing the performance of a detector |
CN109154655A (zh) * | 2017-12-18 | 2019-01-04 | 深圳市大疆创新科技有限公司 | 目标信号检测方法、设备、无人机及农业无人机 |
CN108519511B (zh) * | 2018-03-28 | 2019-12-27 | 电子科技大学 | 一种线性调频信号频率特征参数的时域测量方法 |
CN108919213B (zh) * | 2018-08-06 | 2022-07-29 | 中国航空工业集团公司雷华电子技术研究所 | 一种机载雷达同步在线分析系统 |
CN110873877B (zh) | 2019-04-25 | 2021-04-23 | 北京航空航天大学 | 目标运动轨迹的确定方法及装置 |
CN110531332B (zh) * | 2019-07-02 | 2023-02-21 | 中国航空工业集团公司雷华电子技术研究所 | 一种基于分段门限的低空慢速小目标检测方法 |
CN110501684B (zh) * | 2019-08-23 | 2022-12-23 | 北京航天朗智科技有限公司 | 雷达数据处理装置及雷达数据处理方法 |
CN110764081A (zh) * | 2019-11-05 | 2020-02-07 | 北京理工大学 | 一种精密跟踪测量雷达信号的处理系统 |
CN112014820A (zh) * | 2020-08-27 | 2020-12-01 | 南京矽典微系统有限公司 | 基于运动趋势意图判断的信号处理方法及系统 |
CN113406618B (zh) * | 2021-06-22 | 2022-10-25 | 哈尔滨工业大学 | 一种tws雷达多目标连续跟踪方法 |
CN113341377B (zh) * | 2021-08-02 | 2021-11-05 | 网络通信与安全紫金山实验室 | 雷达基带模块及雷达系统 |
CN113985393B (zh) * | 2021-10-25 | 2024-04-16 | 南京慧尔视智能科技有限公司 | 一种目标检测方法、装置及系统 |
CN116383717B (zh) * | 2023-03-30 | 2024-04-30 | 中国人民解放军93209部队 | 一种智能化综合无人机识别系统及方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07159515A (ja) * | 1993-12-08 | 1995-06-23 | Mitsubishi Electric Corp | レーダ信号処理方式 |
CN102288941A (zh) * | 2011-05-19 | 2011-12-21 | 北京航空航天大学 | 一种基于fpga和dsp的中频lfm-pd雷达信号实时处理系统及处理方法 |
CN102798863A (zh) * | 2012-07-04 | 2012-11-28 | 西安电子科技大学 | 基于汽车防撞雷达的道路中央隔离带检测方法 |
CN104076352A (zh) * | 2014-06-27 | 2014-10-01 | 电子科技大学 | 低截获测速方法及雷达装置 |
CN104330791A (zh) * | 2014-10-24 | 2015-02-04 | 上海无线电设备研究所 | 一种基于频域切变的相参积累方法 |
JP2015230284A (ja) * | 2014-06-06 | 2015-12-21 | 株式会社東芝 | レーダ装置及びそのレーダ信号処理方法 |
CN105445716A (zh) * | 2015-11-25 | 2016-03-30 | 上海无线电设备研究所 | 一种基于旋翼回波信号时域自相关的悬停直升机检测方法 |
CN105548970A (zh) * | 2015-12-11 | 2016-05-04 | 无锡市雷华科技有限公司 | 一种飞鸟探测雷达处理机 |
CN105572670A (zh) * | 2015-12-11 | 2016-05-11 | 无锡市雷华科技有限公司 | 一种飞鸟探测雷达系统 |
-
2016
- 2016-12-30 CN CN201611261993.3A patent/CN106597411B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07159515A (ja) * | 1993-12-08 | 1995-06-23 | Mitsubishi Electric Corp | レーダ信号処理方式 |
CN102288941A (zh) * | 2011-05-19 | 2011-12-21 | 北京航空航天大学 | 一种基于fpga和dsp的中频lfm-pd雷达信号实时处理系统及处理方法 |
CN102798863A (zh) * | 2012-07-04 | 2012-11-28 | 西安电子科技大学 | 基于汽车防撞雷达的道路中央隔离带检测方法 |
JP2015230284A (ja) * | 2014-06-06 | 2015-12-21 | 株式会社東芝 | レーダ装置及びそのレーダ信号処理方法 |
CN104076352A (zh) * | 2014-06-27 | 2014-10-01 | 电子科技大学 | 低截获测速方法及雷达装置 |
CN104330791A (zh) * | 2014-10-24 | 2015-02-04 | 上海无线电设备研究所 | 一种基于频域切变的相参积累方法 |
CN105445716A (zh) * | 2015-11-25 | 2016-03-30 | 上海无线电设备研究所 | 一种基于旋翼回波信号时域自相关的悬停直升机检测方法 |
CN105548970A (zh) * | 2015-12-11 | 2016-05-04 | 无锡市雷华科技有限公司 | 一种飞鸟探测雷达处理机 |
CN105572670A (zh) * | 2015-12-11 | 2016-05-11 | 无锡市雷华科技有限公司 | 一种飞鸟探测雷达系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106597411A (zh) | 2017-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106597411B (zh) | 雷达信号处理方法 | |
CN110161474B (zh) | 一种基于数字波束形成技术的低慢小目标探测方法 | |
US7154434B1 (en) | Anti-personnel airborne radar application | |
US7898457B2 (en) | System and method for processing imagery from synthetic aperture systems | |
CN108196240B (zh) | 一种适用于csar成像的地面动目标轨迹重构方法 | |
US20120093359A1 (en) | Batch Detection Association for Enhanced Target Descrimination in Dense Detection Environments | |
CN101000376A (zh) | 双基地合成孔径雷达的双门限恒虚警运动目标检测方法 | |
CN106468772B (zh) | 一种基于距离-多普勒量测的多站雷达人体跟踪方法 | |
Bluestein et al. | Doppler radar analysis of the Northfield, Texas, tornado of 25 May 1994 | |
CN115508862A (zh) | 一种基于激光雷达的机场风切变预警方法、电子设备及可读介质 | |
Radhakrishna et al. | Fuzzy logic algorithms to identify birds, precipitation, and ground clutter in S-band radar data using polarimetric and nonpolarimetric variables | |
CN112198490B (zh) | 一种机载毫米波测云雷达地杂波抑制方法、装置及雷达 | |
CN113406639A (zh) | 基于车载移动式雷达的fod检测方法、系统及介质 | |
McCann et al. | Visualising the aspect-dependent radar cross section of seabirds over a tidal energy test site using a commercial marine radar system | |
CN112684425B (zh) | 一种恒虚警检测后的目标二次甄别方法 | |
CN108508413A (zh) | 一种基于概率统计的低信噪比条件下目标检测方法 | |
CN111796270A (zh) | 周界安防雷达横向穿越目标检测方法、系统、介质及设备 | |
Kelly et al. | Challenges and solutions of remote sensing at offshore wind energy developments | |
Hong et al. | Single-channel UWB SAR ground moving targets detection method using change detection based on single-pass sub-aperture images | |
Tian et al. | A novel DP-TBD algorithm for tracking slowly maneuvering targets using ViSAR image sequences | |
CN114063057A (zh) | 一种对海多功能雷达信号处理方法 | |
CN112147609B (zh) | 多假设模糊匹配Radon变换的高重频雷达高速机动目标检测方法 | |
JP7389629B2 (ja) | 長時間積分による検出方法および対応するレーダシステム | |
Putri et al. | Development of FMCW Radar Signal Processing for High-Speed Railway Collision Avoidance | |
Yeşil et al. | A radar target classification approach for bird elimination |
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 |