CN107229036B - 基于信号处理的多通道阵列雷达幅相误差在线检测方法 - Google Patents

基于信号处理的多通道阵列雷达幅相误差在线检测方法 Download PDF

Info

Publication number
CN107229036B
CN107229036B CN201710390917.0A CN201710390917A CN107229036B CN 107229036 B CN107229036 B CN 107229036B CN 201710390917 A CN201710390917 A CN 201710390917A CN 107229036 B CN107229036 B CN 107229036B
Authority
CN
China
Prior art keywords
pulse
reference target
pulses
kth
amplitude
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
Application number
CN201710390917.0A
Other languages
English (en)
Other versions
CN107229036A (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN201710390917.0A priority Critical patent/CN107229036B/zh
Publication of CN107229036A publication Critical patent/CN107229036A/zh
Application granted granted Critical
Publication of CN107229036B publication Critical patent/CN107229036B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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/40Means for monitoring or calibrating

Abstract

本发明公开了一种基于信号处理的多通道阵列雷达幅相误差在线检测方法,其主要思路为:确定多通道阵列雷达,所述多通道阵列雷达的检测范围内存在目标,将该目标记为参照目标,多通道阵列雷达向参照目标发射P个脉冲并接收参照目标反射回来的回波信号,记为参照目标反射的回波信号数据;分别得到参照目标反射P个脉冲的回波信号数据集合s和P个脉冲的波束扫描图;分别获知第1个脉冲是否正常至第P个脉冲是否正常;分别记P'为正常的脉冲个数,P”为不正常的脉冲个数,0≤P'≤P,0≤P”≤P,P'+P”=P,P'、P”分别为大于等于0的正整数;根据参照目标反射P个脉冲的回波信号数据集合s对P”个脉冲与P'个脉冲分别进行幅相误差分析,进而得到基于信号处理的多通道阵列雷达幅相误差在线检测结果。

Description

基于信号处理的多通道阵列雷达幅相误差在线检测方法
技术领域
本发明属于雷达信号处理技术领域,特别涉及一种基于信号处理的多通道阵列雷达幅相误差在线检测方法,适用于对多通道阵列雷达的阵列通道幅度失配和相位失配进行检测和校正。
背景技术
雷达是利用电磁波发现目标并测定目标具体位置的设备,在军事和生活上都有重要的应用,雷达的性能状态决定了其工作的可信度,雷达的自检系统能够对雷达的性能好坏进行自动检测;当雷达的状态出现问题导致性能下降,该雷达的自检系统就会反馈给雷达的使用者,雷达的使用者会根据自检结果对雷达进行维修,进而保证雷达能够正常工作。
雷达的自检是对整个雷达工作状态的一种测评,且雷达的自检必须可靠、稳定,而且具有全面性;特别是现在的雷达系统都很复杂,自检和校正就都显得尤为重要;而目前雷达的自检方法都很简单,归纳如下:
(1)状态采样法,也称为旁检模式,该状态采样法通过对雷达里面的各个元器件和节点分别进行温度、电压、电流等采样,进而判断雷达是否能够工作正常。
但该状态采样法方法只能检测出雷达能否接收到回波信号,浮于表面,不能检测出回波信号中的相位误差、幅度误差等更深层次的问题;随着雷达技术的不断发展,对雷达探测性能的要求越来越高,使得采用数字波束形成技术的多通道雷达,成为雷达发展的一个重要方向;多通道雷达在接收到回波信号后,形成多通道雷达回波数据,如果出现多通道的相位、幅度不一致,相干处理多通道回波数据就会收到影响,严重影响多通道雷达的性能;为了检测出相位和幅度的误差,需要对回波信号进行信号处理,并对信号处理后的回波信号进行分析,找到误差并校正。
(2)模拟信号法,也称为自检模式,该模拟信号法通过雷达接收机接收模拟雷达回波信号,并对接收到的模拟雷达回波信号进行分析,将分析的结果与接收的模拟雷达回波信号的参数进行比对,获得模拟雷达回波信号中的相位和幅度误差,从而判断雷达是否故障。
该模拟信号法虽然能检测出模拟雷达回波信号中的相位和幅度误差,但不能实时监测雷达的工作状态,不能及时诊断出雷达的故障,也会影响雷达系统的工作。
发明内容
针对上述现有技术存在的不足,本发明的目的在于提出一种基于信号处理的多通道阵列雷达幅相误差在线检测方法,该种基于信号处理的多通道阵列雷达幅相误差在线检测方法能够提高雷达的测向精度。
本发明的主要思路:本发明利用广播式自动相关监视(Automatic DependentSurveillance-Broadcas t,ADS-B)向雷达提供已知目标的参数,雷达向该目标发射脉冲并接收回波信号,通过对回波信号进行信号处理得到参数,将得到的参数与已知目标的参数进行比对,进而检测出雷达在幅度和相位上的误差。
为达到上述技术目的,本发明采用如下技术方案予以实现。
一种基于信号处理的多通道阵列雷达幅相误差在线检测方法,包括以下步骤:
步骤1,确定多通道阵列雷达,所述多通道阵列雷达的检测范围内存在目标,将该目标记为参照目标;多通道阵列雷达向参照目标发射P个脉冲并接收参照目标反射回来的回波信号,记为参照目标反射的回波信号数据;
初始化:令k表示第k个脉冲,k∈{1,2,…,P},P表示多通道阵列雷达向参照目标发射的脉冲总个数,P为大于0的正整数;
步骤2,将多通道阵列雷达对参照目标发射第k个脉冲并接收参照目标反射回来的回波信号,记为参照目标反射第k个脉冲的回波信号数据sk
步骤3,令k分别取1至P,重复执行步骤2,进而分别得到参照目标反射第1个脉冲的回波信号数据s1至参照目标反射第P个脉冲的回波信号数据sP,记为参照目标反射P个脉冲的回波信号数据集合s,s={s1,s2,…,sP};然后将k置为1;
步骤4,对参照目标反射第k个脉冲的回波信号数据sk进行数字信号处理,得到数字信号处理后参照目标反射第k个脉冲的回波信号数据Xk(θ),记为第k个脉冲的波束扫描图;
步骤5,令k加1,重复执行步骤4,进而分别得到第1个脉冲的波束扫描图至第P个脉冲的波束扫描图,记为P个脉冲的波束扫描图;然后将k置为1;
步骤6,根据第k个脉冲的波束扫描图,获知第k个脉冲是否正常;其中第k个脉冲正常为第k个脉冲的相位和第k个脉冲的幅值都不存在误差,第k个脉冲不正常为第k个脉冲的相位和第k个脉冲的幅值中至少一个存在误差;
步骤7,令k加1,重复执行步骤6,进而分别获知第1个脉冲是否正常至第P个脉冲是否正常;
分别记P'为正常的脉冲个数,P”为不正常的脉冲个数,0≤P'≤P,0≤P”≤P,P'+P”=P,P'、P”分别为大于等于0的正整数;
步骤8,根据参照目标反射P个脉冲的回波信号数据集合s对P”个脉冲分别进行幅相误差分析,进而得到基于信号处理的多通道阵列雷达幅相误差在线检测结果。
本发明与现有的技术相比有以下优点:
第一,本发明方法根据广播式自动相关监视(Automatic DependentSurveillance-Broadcas t,ADS-B)信息系统共享的特点,合理利用ADS-B的数据资源,能够实时监测雷达的幅相误差,实现空中监管的协同能力。
第二,针对多通道阵列雷达中幅相产生误差的问题,利用波束扫描和比对校正,适用阵型广泛,简单可行,能找到其他自检方法检测不到的更深层次的问题。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细说明。
图1是本发明的一种基于信号处理的多通道阵列雷达幅相误差在线检测方法实现流程图;
图2是利用仿真条件1对第1个脉冲的波束扫描图;
图3是利用仿真条件1对第3个脉冲的波束扫描图;
图4是利用仿真条件1对第1个脉冲第10个阵元的接收回波示意图;
图5是利用仿真条件1对第3个脉冲第10个阵元的接收回波示意图;
图6是利用仿真条件2对第1个脉冲的波束扫描图;
图7是利用仿真条件2对第3个脉冲的波束扫描图;
图8是利用仿真条件2对第1个脉冲第20个阵元的接收回波示意图;
图9是利用仿真条件2对第3个脉冲第20个阵元的接收回波示意图;
图10是利用实验条件对第49个脉冲的波束扫描图;
图11是利用实验条件对第50个脉冲的波束扫描图;
图12(a)是利用实验条件对第35个阵元分别接收第49个脉冲的幅值对比图;
图12(b)是利用实验条件对第35个阵元分别接收第50个脉冲的幅值对比图;
图12(c)是利用实验条件对第35个阵元分别接收第51个脉冲的幅值对比图;
图13(a)利用实验条件对第36个阵元分别接收第49个脉冲的幅值对比图;
图13(b)利用实验条件对第36个阵元分别接收第50个脉冲的幅值对比图;
图13(c)利用实验条件对第36个阵元分别接收第51个脉冲的幅值对比图;
图14(a)是利用实验条件对第37个阵元分别接收第49个脉冲的幅值对比图;
图14(b)是利用实验条件对第37个阵元分别接收第50个脉冲的幅值对比图;
图14(c)是利用实验条件对第37个阵元分别接收第51个脉冲的幅值对比图。
具体实施方式
参照图1,为本发明的一种基于信号处理的多通道阵列雷达幅相误差在线检测方法实现流程图;其中所述基于信号处理的多通道阵列雷达幅相误差在线检测方法,包括以下步骤:
步骤1,在无杂波区选取合适的大目标,所述合适的大目标为雷达散射截面(RadarCross Section,RCS)在设定范围以上的目标,其中设定范围以上为100m2以上,一般为民航飞机;目前民航飞机基本都装载了广播式自动相关监视(Automatic DependentSurveillance-Broadcas t,ADS-B)设备,地面上的雷达也都配备有广播式自动相关监视(Automatic Dependent Surveillance-Broadcas t,ADS-B)设备。
步骤2,多通道阵列雷达是信号的产生和接收全部采用多通道并行处理的数字化阵列雷达,多通道阵列雷达是由N个阵元组成的等距线阵,任意两个相邻阵元之间的阵元间距d'均相等;多通道阵列雷达能接收到两方面的数据,一是多通道阵列雷达发射脉冲串后接收到的回波数据;二是参照目标通过ADS-B数据链向多通道阵列雷达告知其自身的精确位置,多通道阵列雷达即可获得该参照目标的方位、距离、速度、航向信息。
多通道阵列雷达向参照目标发射P个脉冲并接收参照目标反射回来的回波信号,记为参照目标反射的回波信号数据。
初始化:令k表示第k个脉冲,k∈{1,2,…,P},P表示多通道阵列雷达向参照目标发射的脉冲总个数,P为大于0的正整数。
步骤3,将多通道阵列雷达对参照目标发射第k个脉冲并接收参照目标反射回来的回波信号,记为参照目标反射第k个脉冲的回波信号数据sk,其表达式为:
sk=akexp(j2πfdk△t)a(θ0)+nk
其中,k∈{1,2,…,P},a(θ0)表示参照目标反射的回波信号数据的导向矢量,
Figure BDA0001307492740000051
exp表示指数函数操作,j表示虚数单位,θ0表示参照目标所在方向与多通道阵列雷达法线方向之间的夹角,λ表示参照目标反射的回波信号数据的波长;akexp(j2πfdk△t)表示sk的复包络,sk表示参照目标反射第k个脉冲的回波信号数据,ak表示sk的幅值,fd为sk的多普勒频率,△t为多通道阵列发射脉冲的脉冲间隔,nk表示多通道阵列雷达接收sk时伴有的噪声,nk=[nk1,nk2,…,nkl,,…,nkN],nkl,表示多通道阵列雷达接收第k个脉冲时第l个阵元伴有的噪声,l∈{1,2,…,N},N表示多通道阵列雷达包含的阵元个数,N、P、d'分别为大于0的正整数,exp表示指数函数操作,j表示虚数单位,∈表示属于。
步骤4,令k分别取1至P,重复执行步骤3,进而分别得到参照目标反射第1个脉冲的回波信号数据s1至参照目标反射第P个脉冲的回波信号数据sP,记为参照目标反射P个脉冲的回波信号数据集合s,s={s1,s2,…,sP};然后将k置为1。
步骤5,对参照目标反射第k个脉冲的回波信号数据sk进行数字信号处理,具体为:对参照目标反射第k个脉冲的回波信号数据sk进行数字波束形成(DBF)处理,即对参照目标反射第k个脉冲的回波信号数据sk进行加权求和后形成θ方向的接收波束,θ表示对参照目标反射第k个脉冲的回波信号数据sk进行波束扫描的扫描角度,θ的取值范围为[α12],α1表示参照目标所在角度范围最小值,α2表示参照目标所在角度范围最大值;进而得到数字信号处理后参照目标反射第k个脉冲的回波信号数据Xk(θ),其表达式为:
Xk(θ)=|wH(θ)sk|
其中,上标H表示共轭转置操作,| |表示求绝对值操作,w(θ)表示对参照目标反射第k个脉冲的回波信号数据sk进行波束扫描的扫描权矢量,
Figure BDA0001307492740000052
λ表示参照目标反射的回波信号数据的波长,d'表示多通道阵列雷达中任意两个相邻阵元之间的阵元间距,sin表示求正弦操作,exp表示指数函数操作,j表示虚数单位,N表示多通道阵列雷达包含的阵元个数。
所述数字信号处理后参照目标反射第k个脉冲的回波信号数据Xk(θ)为数字波束形成(DBF)处理后参照目标反射第k个脉冲的回波信号数据,记为第k个脉冲的波束扫描图。
步骤6,令k加1,重复执行步骤5,分别得到第1个脉冲的波束扫描图至第P个脉冲的波束扫描图,记为P个脉冲的波束扫描图,P表示多通道阵列雷达向参照目标发射的脉冲总个数;然后将k置为1。
步骤7,根据第k个脉冲的波束扫描图,获知第k个脉冲是否正常;其中第k个脉冲正常为第k个脉冲的相位和第k个脉冲的幅值都不存在误差,第k个脉冲不正常为第k个脉冲的相位和第k个脉冲的幅值中至少一个存在误差。
具体地,根据第k个脉冲的波束扫描图,分别得到第k个脉冲最大增益时对应的角度
Figure BDA0001307492740000061
第k个脉冲的波束扫描图的主瓣宽度
Figure BDA0001307492740000062
以及第k个脉冲主瓣与第k个脉冲副瓣之比
Figure BDA0001307492740000063
所述第k个脉冲主瓣为第k个脉冲的波束扫描图的第一峰值,主瓣宽度为-3dB宽度;所述第k个脉冲副瓣为第k个脉冲的波束扫描图的第二峰值。
(1)将第k个脉冲最大增益时对应的角度
Figure BDA0001307492740000064
与θ0做比较,θ0表示参照目标真实角度,所述参照目标真实角度为参照目标所在方向与多通道阵列雷达法线方向之间的夹角;进而得到第k个脉冲最大增益时对应的角度与参照目标真实角度之间的误差εk
Figure BDA0001307492740000065
若εk0,ε0为设置的实际情况下最大角度误差值,则第k个脉冲最大增益时对应的角度
Figure BDA0001307492740000066
在正常范围内,所述第k个脉冲最大增益时对应的角度
Figure BDA0001307492740000067
在正常范围内为第k个脉冲的波束扫描图最大增益时对应的角度与参照目标方位一致。
(2)将第k个脉冲的波束扫描图的主瓣宽度
Figure BDA0001307492740000068
与B0作比较,B0为设置的理想情况下自适应波束形成的主瓣宽度,得到第k个脉冲的波束扫描图的主瓣宽度与理想情况下自适应波束形成的主瓣宽度的误差δk
Figure BDA0001307492740000069
若δk0,δ0表示设置的实际情况下最大主瓣宽度误差,则第k个脉冲的波束扫描图的主瓣宽度
Figure BDA00013074927400000610
在正常范围内,其中正常范围为δk0
(3)将第k个脉冲主瓣与第k个脉冲副瓣之比
Figure BDA00013074927400000611
与A0作比较,A0为设置的理想情况下主副瓣比,
Figure BDA0001307492740000071
得到第k个脉冲的主副瓣比与理想情况下主副瓣比的误差σk
Figure BDA0001307492740000072
其中第k个脉冲的主副瓣比为第k个脉冲主瓣与第k个脉冲副瓣之比
Figure BDA0001307492740000073
若σk0,σ0表示设置的实际情况下最大主副瓣之比误差,则第k个脉冲主瓣与第k个脉冲副瓣之比
Figure BDA0001307492740000074
在正常范围内。
若第k个脉冲最大增益时对应的角度
Figure BDA0001307492740000075
与参照目标方位一致、第k个脉冲的波束扫描图的主瓣宽度
Figure BDA0001307492740000076
在正常范围内,以及第k个脉冲主瓣与第k个脉冲副瓣之比
Figure BDA0001307492740000077
也在正常范围内,则第k个脉冲正常,即第k个脉冲没有幅值和相位的误差;若第k个脉冲最大增益时对应的角度
Figure BDA0001307492740000078
第k个脉冲的波束扫描图的主瓣宽度
Figure BDA0001307492740000079
第k个脉冲主瓣与第k个脉冲副瓣之比
Figure BDA00013074927400000710
至少一个不在正常范围内,则第k个脉冲不正常,所述第k个脉冲不正常为第k个脉冲的相位和第k个脉冲的幅值中至少一个存在误差。
步骤8,令k分别取1至P,重复步骤7,进而分别获知第1个脉冲是否正常至第P个脉冲是否正常。
分别记P'为正常的脉冲个数,P”为不正常的脉冲个数,0≤P'≤P,0≤P”≤P,P'+P”=P,P'、P”分别为大于等于0的正整数。
步骤9,对方向不一致的脉冲做进一步分析:根据参照目标反射P个脉冲的回波信号数据集合s对P”个脉冲分别进行幅相误差分析,如果P”个脉冲与P'个脉冲的幅值分别相同,则说明P”个脉冲各自的相位分别产生了误差;如果P”个脉冲与P'个脉冲的幅值分别不相同,则说明P”个脉冲各自的相位或P”个脉冲各自幅值或P”个脉冲各自的相位和幅值都产生了误差,记为基于信号处理的多通道阵列雷达幅相误差在线检测结果。
9.1令syl表示P”个不正常的脉冲中第y个脉冲、第l个阵元接收的回波数据,y∈{1,2,…,P”},l∈{1,2,…,N};令sy'l表示P'个正常的脉冲中第y'个脉冲、第l个阵元接收的回波数据,y'∈{1,2,…,P'},l∈{1,2,…,N},N表示多通道阵列雷达包含的阵元个数,y和l的初始值分别为1。
9.2获取P”个不正常的脉冲中第y个脉冲、第l个阵元接收的回波数据的误差分析结果,其过程为:
如果P”个不正常的脉冲中第y个脉冲、第l个阵元接收的回波数据幅值与P'个正常的脉冲中第y'个脉冲、第l个阵元接收的回波数据幅值相等,则说明P”个不正常的脉冲中第y个脉冲、第l个阵元接收的回波数据相位产生了误差:如果P”个不正常的脉冲中第y个脉冲、第l个阵元接收的回波数据幅值与P'个正常的脉冲中第y'个脉冲、第l个阵元接收的回波数据幅值不相等,则说明P”个不正常的脉冲中第y个脉冲、第l个阵元接收的回波数据相位和幅值中至少一个产生了误差。
9.3令l分别取1至N,重复9.2,分别得到P”个不正常的脉冲中第y个脉冲、第1个阵元接收的回波数据的误差分析结果至P”个不正常的脉冲中第y个脉冲、第N个阵元接收的回波数据的误差分析结果,记为P”个不正常的脉冲中第y个脉冲的误差分析结果,所述P”个不正常的脉冲中第y个脉冲的误差分析结果包含N个误差,如果N个误差都是相位误差,则P”个不正常的脉冲中第y个脉冲是相位误差;如果N个误差不都是相位误差,则P”个不正常的脉冲中第y个脉冲是相位误差和幅值误差。
9.4令y分别取1至P”,重复9.3,分别得到P”个不正常的脉冲中第1个脉冲的误差分析结果至P”个不正常的脉冲中第P”个脉冲的误差分析结果,记为基于信号处理的多通道阵列雷达幅相误差在线检测结果。
通过以下计算机仿真对本发明效果作进一步验证说明。
(一)仿真条件:
仿真条件1,假设多通道阵列雷达为由50个阵元组成的等距线阵,阵元间距d'等于半波长,参照目标与多通道阵列雷达阵列法线方向的夹角为20°,信噪比为0dB,波束扫描的扫描范围为[10°,30°];多通道阵列雷达向参照目标发射的脉冲总个数为5,每个脉冲有1000个距离单元,其中第3个脉冲的幅值和相位都存在误差。
仿真条件2,假设多通道阵列雷达为由50个阵元组成的等距线阵,阵元间距d'等于半波长,参照目标与多通道阵列雷达阵列法线方向的夹角为20°,信噪比为0dB,波束扫描的扫描范围为[10°,30°];多通道阵列雷达向参照目标发射的脉冲总个数为5,每个脉冲有1000个距离单元,其中第3个脉冲仅相位存在误差。
(二)仿真内容
仿真1,在仿真条件1下,对每个脉冲进行波束扫描,可以得到5个脉冲的波束扫描图,选取第1个脉冲的波束扫描图,如图2所示,图2中横坐标为角度,单位为度,纵坐标为幅度,单位为dB。第3个脉冲没有检测到目标,如图3所示,图3中横坐标为角度,单位为度,纵坐标为幅度,单位为dB。
将每个阵元接收的第3个回波脉冲幅值与其他脉冲进行对比,第3个脉冲与其他脉冲的幅值不同;选取第1个脉冲第10个阵元的接收回波,如图4所示,图4中横坐标为距离单元,单位为1,纵坐标为幅度,单位为1;选取第3个脉冲第10个阵元的接收回波,如图5所示,图5中横坐标为距离单元,单位为1,纵坐标为幅度,单位为1。
由图2、图3、图4、图5可见,第3个脉冲的相位和幅度产生了误差。
仿真2,在仿真条件2下,对每个脉冲进行波束扫描,可以得到5个脉冲的波束扫描图,选取第1个脉冲的波束扫描图,如图6所示,图6中横坐标为角度,单位为度,纵坐标为幅度,单位为dB;第3个脉冲没有检测到目标,如图7所示,图7中横坐标为角度,单位为度,纵坐标为幅度,单位为dB。
将每个阵元接收的第3个回波脉冲幅值与其他脉冲进行对比,第3个脉冲与其他脉冲的幅值相同;选取第1个脉冲第20个阵元的接收回波,如图8所示,图8中横坐标为距离单元,单位为1,纵坐标为幅度,单位为1。选取第3个脉冲第20个阵元的接收回波,如图9所示,图9中横坐标为距离单元,单位为1,纵坐标为幅度,单位为1。
由图6、图7、图8、图9可见,第3个脉冲的相位产生了误差。
通过实测数据对本发明效果作进一下验证说明。
(三)实验条件
多通道阵列雷达为圆周阵列,N个阵元均匀分布在半径为a的圆周上,多通道阵列雷达共接收到P个脉冲,每个脉冲分别包含L个距离单元。
(四)实验内容
在实验条件(三)下,对每个脉冲进行波束扫描,可以得到P个脉冲的波束扫描图,选取第49个脉冲的波束扫描图,如图10所示,图10中横坐标为角度,单位为度,纵坐标为幅度,单位为dB;第23、第50、第88、第132、第201、第245个脉冲都没有检测到目标。
选取第50个脉冲的波束扫描图,如图11所示,图11中横坐标为角度,单位为度,纵坐标为幅度,单位为dB。
将每个阵元接收的第23、第50、第88、第132、第201、第245个脉冲的幅值与其他脉冲进行比较,这6个脉冲在幅值上与其他脉冲相比有较大的误差,选取第49、第50、第51个连续脉冲进行对比,阵元35、第49、第50、第51个脉冲幅值相同,如图12(a)、图12(b)和图12(c)所示,图12(a)、图12(b)和图12(c)中的横坐标分别为距离单元,单位为1,纵坐标分别为幅度,单位为1;阵元36、第49个脉冲幅值大于第50个脉冲,如图13(a)、图13(b)和图13(c)所示,图13(a)、图13(b)和图13(c)中的横坐标分别为距离单元,单位为1,纵坐标分别为幅度,单位为1;阵元37、第49个脉冲匹配滤波后幅值小于第50个脉冲,如图14(a)、图14(b)和图14(c)所示,图14(a)、图14(b)和图14(c)中的横坐标分别为距离单元,单位为1,纵坐标分别为幅度,单位为1。
由图10、图11、图12(a)、图12(b)、图12(c)、图13(a)、图13(b)、图13(c)、图14(a)、图14(b)和图14(c)可见,第50、第88、第132、第201、第245个脉冲的相位和幅值都产生了误差。
综上所述,仿真实验验证了本发明的正确性,有效性和可靠性。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围;这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (4)

1.一种基于信号处理的多通道阵列雷达幅相误差在线检测方法,其特征在于,包括以下步骤:
步骤1,确定多通道阵列雷达,所述多通道阵列雷达的检测范围内存在目标,将该目标记为参照目标;多通道阵列雷达向参照目标发射P个脉冲并接收参照目标反射回来的回波信号,记为参照目标反射的回波信号数据;
初始化:令k表示第k个脉冲,k∈{1,2,…,P},P表示多通道阵列雷达向参照目标发射的脉冲总个数,P为大于0的正整数;
步骤2,将多通道阵列雷达对参照目标发射第k个脉冲并接收参照目标反射回来的回波信号,记为参照目标反射第k个脉冲的回波信号数据sk
步骤3,令k分别取1至P,重复执行步骤2,进而分别得到参照目标反射第1个脉冲的回波信号数据s1至参照目标反射第P个脉冲的回波信号数据sP,记为参照目标反射P个脉冲的回波信号数据集合s,s={s1,s2,…,sP};然后将k置为1;
步骤4,对参照目标反射第k个脉冲的回波信号数据sk进行数字信号处理,得到数字信号处理后参照目标反射第k个脉冲的回波信号数据Xk(θ),记为第k个脉冲的波束扫描图;
步骤5,令k加1,重复执行步骤4,进而分别得到第1个脉冲的波束扫描图至第P个脉冲的波束扫描图,记为P个脉冲的波束扫描图;然后将k置为1;
步骤6,根据第k个脉冲的波束扫描图,获知第k个脉冲是否正常;其中第k个脉冲正常为第k个脉冲的相位和第k个脉冲的幅值都不存在误差,第k个脉冲不正常为第k个脉冲的相位和第k个脉冲的幅值中至少一个存在误差;
过程为:
根据第k个脉冲的波束扫描图,分别得到第k个脉冲最大增益时对应的角度
Figure FDA0002379962070000011
第k个脉冲的波束扫描图的主瓣宽度
Figure FDA0002379962070000012
以及第k个脉冲主瓣与第k个脉冲副瓣之比
Figure FDA0002379962070000013
所述第k个脉冲主瓣为第k个脉冲的波束扫描图的第一峰值,主瓣宽度为-3dB宽度;所述第k个脉冲副瓣为第k个脉冲的波束扫描图的第二峰值;
(1)将第k个脉冲最大增益时对应的角度
Figure FDA0002379962070000014
与θ0做比较,θ0表示参照目标真实角度,所述参照目标真实角度为参照目标所在方向与多通道阵列雷达法线方向之间的夹角;进而得到第k个脉冲最大增益时对应的角度与参照目标真实角度之间的误差εk
Figure FDA0002379962070000021
若εk0,ε0为设置的实际情况下最大角度误差值,则第k个脉冲最大增益时对应的角度
Figure FDA0002379962070000022
在正常范围内;
(2)将第k个脉冲的波束扫描图的主瓣宽度
Figure FDA0002379962070000023
与B0作比较,B0为设置的理想情况下自适应波束形成的主瓣宽度,得到第k个脉冲的波束扫描图的主瓣宽度与理想情况下自适应波束形成的主瓣宽度的误差δk
Figure FDA0002379962070000024
若δk0,δ0表示设置的实际情况下最大主瓣宽度误差,则第k个脉冲的波束扫描图的主瓣宽度
Figure FDA0002379962070000025
在正常范围内;
(3)将第k个脉冲主瓣与第k个脉冲副瓣之比
Figure FDA0002379962070000026
与A0作比较,A0为设置的理想情况下主副瓣比,
Figure FDA0002379962070000027
得到第k个脉冲的主副瓣比与理想情况下主副瓣比的误差σk
Figure FDA0002379962070000028
其中第k个脉冲的主副瓣比为第k个脉冲主瓣与第k个脉冲副瓣之比
Figure FDA0002379962070000029
若σk0,σ0表示设置的实际情况下最大主副瓣之比误差,则第k个脉冲主瓣与第k个脉冲副瓣之比
Figure FDA00023799620700000210
在正常范围内;
若第k个脉冲最大增益时对应的角度
Figure FDA00023799620700000211
与参照目标方位一致、第k个脉冲的波束扫描图的主瓣宽度
Figure FDA00023799620700000212
在正常范围内,以及第k个脉冲主瓣与第k个脉冲副瓣之比
Figure FDA00023799620700000213
也在正常范围内,则第k个脉冲正常,即第k个脉冲没有幅值和相位的误差;若第k个脉冲最大增益时对应的角度第k个脉冲的波束扫描图的主瓣宽度
Figure FDA00023799620700000215
第k个脉冲主瓣与第k个脉冲副瓣之比
Figure FDA00023799620700000216
至少一个不在正常范围内,则第k个脉冲不正常,所述第k个脉冲不正常为第k个脉冲的相位和第k个脉冲的幅值中至少一个存在误差;
步骤7,令k加1,重复执行步骤6,进而分别获知第1个脉冲是否正常至第P个脉冲是否正常;
分别记P'为正常的脉冲个数,P”为不正常的脉冲个数,0≤P'≤P,0≤P”≤P,P'+P”=P,P'、P”分别为大于等于0的正整数;
步骤8,根据参照目标反射P个脉冲的回波信号数据集合s对P”个脉冲分别进行幅相误差分析,进而得到基于信号处理的多通道阵列雷达幅相误差在线检测结果。
2.如权利要求1所述的一种基于信号处理的多通道阵列雷达幅相误差在线检测方法,其特征在于,在步骤2中,所述参照目标反射第k个脉冲的回波信号数据sk,其表达式为:
sk=akexp(j2πfdk△t)a(θ0)+nk
其中,k∈{1,2,…,P},a(θ0)表示参照目标反射的回波信号数据的导向矢量,
Figure FDA0002379962070000031
exp表示指数函数操作,j表示虚数单位,θ0表示参照目标所在方向与多通道阵列雷达法线方向之间的夹角,λ表示参照目标反射的回波信号数据的波长;ak exp(j2πfdk△t)表示sk的复包络,sk表示参照目标反射第k个脉冲的回波信号数据,ak表示sk的幅值,fd为sk的多普勒频率,△t为多通道阵列发射脉冲的脉冲间隔,nk表示多通道阵列雷达接收sk时伴有的噪声,nk=[nk1,nk2,…,nkl,…,nkN],nkl,表示多通道阵列雷达接收第k个脉冲时第l个阵元伴有的噪声,l∈{1,2,…,N},N表示多通道阵列雷达包含的阵元个数,N、P、d'分别为大于0的正整数,exp表示指数函数操作,j表示虚数单位,∈表示属于。
3.如权利要求1所述的一种基于信号处理的多通道阵列雷达幅相误差在线检测方法,其特征在于,在步骤4中,所述对参照目标反射第k个脉冲的回波信号数据sk进行数字信号处理,具体为对参照目标反射第k个脉冲的回波信号数据sk进行数字波束形成处理;
所述数字信号处理后参照目标反射第k个脉冲的回波信号数据Xk(θ),其表达式为:
Xk(θ)=|wH(θ)sk|
其中,上标H表示共轭转置操作,||表示求绝对值操作,w(θ)表示对参照目标反射第k个脉冲的回波信号数据sk进行波束扫描的扫描权矢量,
Figure FDA0002379962070000032
λ表示参照目标反射的回波信号数据的波长,d'表示多通道阵列雷达中任意两个相邻阵元之间的阵元间距,sin表示求正弦操作,exp表示指数函数操作,j表示虚数单位,N表示多通道阵列雷达包含的阵元个数。
4.如权利要求1所述的一种基于信号处理的多通道阵列雷达幅相误差在线检测方法,其特征在于,步骤8的子步骤为:
8.1令syl表示P”个不正常的脉冲中第y个脉冲、第l个阵元接收的回波数据,y∈{1,2,…,P”},l∈{1,2,…,N};令sy'l表示P'个正常的脉冲中第y'个脉冲、第l个阵元接收的回波数据,y'∈{1,2,…,P'},l∈{1,2,…,N},N表示多通道阵列雷达包含的阵元个数,y和l的初始值分别为1;
8.2获取P”个不正常的脉冲中第y个脉冲、第l个阵元接收的回波数据的误差分析结果,其过程为:
如果P”个不正常的脉冲中第y个脉冲、第l个阵元接收的回波数据幅值与P'个正常的脉冲中第y'个脉冲、第l个阵元接收的回波数据幅值相等,则说明P”个不正常的脉冲中第y个脉冲、第l个阵元接收的回波数据相位产生了误差;如果P”个不正常的脉冲中第y个脉冲、第l个阵元接收的回波数据幅值与P'个正常的脉冲中第y'个脉冲、第l个阵元接收的回波数据幅值不相等,则说明P”个不正常的脉冲中第y个脉冲、第l个阵元接收的回波数据相位和幅值中至少一个产生了误差;
8.3令l分别取1至N,重复8.2,分别得到P”个不正常的脉冲中第y个脉冲、第1个阵元接收的回波数据的误差分析结果至P”个不正常的脉冲中第y个脉冲、第N个阵元接收的回波数据的误差分析结果,记为P”个不正常的脉冲中第y个脉冲的误差分析结果,所述P”个不正常的脉冲中第y个脉冲的误差分析结果包含N个误差,如果N个误差都是相位误差,则P”个不正常的脉冲中第y个脉冲是相位误差;如果N个误差中至少一个不是相位误差,则P”个不正常的脉冲中第y个脉冲是相位误差和幅值误差;
8.4令y分别取1至P”,重复8.3,分别得到P”个不正常的脉冲中第1个脉冲的误差分析结果至P”个不正常的脉冲中第P”个脉冲的误差分析结果,记为基于信号处理的多通道阵列雷达幅相误差在线检测结果。
CN201710390917.0A 2017-05-27 2017-05-27 基于信号处理的多通道阵列雷达幅相误差在线检测方法 Active CN107229036B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710390917.0A CN107229036B (zh) 2017-05-27 2017-05-27 基于信号处理的多通道阵列雷达幅相误差在线检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710390917.0A CN107229036B (zh) 2017-05-27 2017-05-27 基于信号处理的多通道阵列雷达幅相误差在线检测方法

Publications (2)

Publication Number Publication Date
CN107229036A CN107229036A (zh) 2017-10-03
CN107229036B true CN107229036B (zh) 2020-07-03

Family

ID=59933822

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710390917.0A Active CN107229036B (zh) 2017-05-27 2017-05-27 基于信号处理的多通道阵列雷达幅相误差在线检测方法

Country Status (1)

Country Link
CN (1) CN107229036B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108732553B (zh) * 2018-06-01 2022-02-01 北京航空航天大学 一种激光雷达波形时刻鉴别方法与在线测距系统
CN109633598B (zh) * 2019-01-25 2022-12-02 西安电子科技大学 基于特征分析的阵列雷达目标检测方法
CN109946665B (zh) * 2019-03-07 2023-04-21 西安电子科技大学 基于阵列雷达的获取真实目标的方法
CN112578353A (zh) * 2020-02-28 2021-03-30 加特兰微电子科技(上海)有限公司 测量目标角度的装置及方法、传感器和设备
CN112578345A (zh) * 2020-11-20 2021-03-30 福瑞泰克智能系统有限公司 一种雷达遮挡检测方法、装置、设备及存储介质
CN113759329B (zh) * 2021-07-23 2023-06-27 西安电子科技大学 基于内外场联合的频率分集阵列雷达幅相误差校正方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102053247A (zh) * 2009-10-28 2011-05-11 中国科学院电子学研究所 一种多基线合成孔径雷达三维成像的相位校正方法
CN102129068A (zh) * 2009-12-08 2011-07-20 中国科学院声学研究所 一种对合成孔径声纳系统相位误差估计的测试系统及方法
CN104076337A (zh) * 2014-07-04 2014-10-01 西安电子科技大学 基于阵元幅相误差校正的机载雷达杂波抑制方法
CN104346956A (zh) * 2013-08-06 2015-02-11 四川信能科技发展有限公司 广播式自动相关监视ads-b自主防虚假目标欺骗办法
CN104865571A (zh) * 2015-04-02 2015-08-26 西安电子科技大学 多通道多子带滑动聚束模式sar成像方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102053247A (zh) * 2009-10-28 2011-05-11 中国科学院电子学研究所 一种多基线合成孔径雷达三维成像的相位校正方法
CN102129068A (zh) * 2009-12-08 2011-07-20 中国科学院声学研究所 一种对合成孔径声纳系统相位误差估计的测试系统及方法
CN104346956A (zh) * 2013-08-06 2015-02-11 四川信能科技发展有限公司 广播式自动相关监视ads-b自主防虚假目标欺骗办法
CN104076337A (zh) * 2014-07-04 2014-10-01 西安电子科技大学 基于阵元幅相误差校正的机载雷达杂波抑制方法
CN104865571A (zh) * 2015-04-02 2015-08-26 西安电子科技大学 多通道多子带滑动聚束模式sar成像方法

Also Published As

Publication number Publication date
CN107229036A (zh) 2017-10-03

Similar Documents

Publication Publication Date Title
CN107229036B (zh) 基于信号处理的多通道阵列雷达幅相误差在线检测方法
CN107121670B (zh) 一种基于合成孔径雷达的反无人机防御方法
Melvin Space-time adaptive radar performance in heterogeneous clutter
EP0301089B1 (en) Method and apparatus for detecting an out of beam condition in a monopulse radar receiver
CN105445701B (zh) Ddma‑mimo雷达目标的单脉冲角度估计方法
CN110488255A (zh) 一种相控阵雷达单脉冲高分辨测角系统及方法
KR20190074224A (ko) 세컨더리 레이더의 안테나 패턴들을 측정하기 위한 방법 및 이러한 방법을 구현하는 세컨더리 레이더
CN104007413A (zh) 考虑信源方位误差的阵列位置误差校正方法
CN112666543B (zh) 一种稀疏阵列tdm-mimo雷达及其校正方法
Gerlach et al. Robust adaptive matched filtering using the FRACTA algorithm
CN111337919A (zh) 用于测量二次雷达的主天线波瓣的方位角准确度和图案的方法以及实现这种方法的雷达
Bosse et al. Model-based multifrequency array signal processing for low-angle tracking
Dumas et al. Self-calibration and antenna grouping for bistatic oceanographic High-Frequency Radars
RU2746799C1 (ru) Статистический способ имитационной защиты загоризонтной радиолокационной станции поверхностной волны
CN109581366B (zh) 一种基于目标导向矢量失配的离散旁瓣杂波识别方法
Zhou et al. Angle measurement accuracy analysis of sum-difference amplitude-comparison monopulse in onshore or shipborne ISAR
Hersey et al. Adaptive ground clutter suppression for conformal array radar systems
Xinghua et al. Performance gain bounds of coherently combining multiple radars in a target-based calibration manner
CN113625265A (zh) 基于波束空间的方位超分辨方法
Howard et al. Experimental results of the complex indicated angle techique for multipath correction
RU2327181C1 (ru) Способ определения местоположения источника радиоизлучения
CN111812608A (zh) 基于mtd脉冲积累及模态分解的雷达目标方位角估计方法
Klintberg et al. Mitigation of ground clutter in airborne bistatic radar systems
Cann Range gate straddling loss and joint probability with partial correlation
CN111142097A (zh) 一种调幅广播外辐射源雷达中的目标方向估计方法

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