技术领域
[0001] 本
发明涉及面目标回波处理技术领域,特别是一种地外天体着陆测量雷达面目标回波信号测距处理方法。
背景技术
[0002] 着陆器的轨迹特点要求着陆雷达能够测量不同波束指向的距离。传统高度计根据波束垂直入射的回波特性,提取回波前沿计算垂向距离;与传统高度计不同,着陆雷达测距
算法必须兼顾垂直入射与斜入射的回波特性。天线发射的波束具有一定的宽度,在着陆器下降过程中天线波束照射一片区域,与照射面的各散射点之间的径向距离不同,具有不同延时的分量相互
叠加构成雷达接收到的回波信号。面目标回波处理算法要准确提取回波的斜距值,其重点是最大程度减小主瓣
波形由于面目标回波矢量叠加造成的剧烈
波动的影响,并且精确提取回波主瓣
重心。
发明内容
[0003] 对此,本发明提供一种地外天体着陆测量雷达面目标回波
信号处理方法,其能够解决传统海面回波处理算法无法有效应用于着陆面目标处理的问题,提取主瓣的中心
位置和重心位置,从而计算得到捕获值和
跟踪值,并且减少由于出现回波毛刺而失去跟踪的情况。
[0004] 根据本发明提供的一种地外天体着陆测量雷达面目标回波信号测距处理方法,其用于提取面目标回波主瓣的中心位置与重心位置。
[0005] 该方法包括:对中频回波进行预处理,其包括A/D
采样、解调滤波
抽取和FFT;所述测距捕获通道包括捕获CFAR和捕获回波中心计算,从而得到回波主瓣中心位置;所述跟踪通道包括跟踪波
门截取、跟踪CFAR、跟踪回波中心计算、主瓣截取和OCOG,从而得到重心位置以提取测距信息。
[0006] 其中,对中频回波进行预处理,经过A/D采样之后的雷达接收信号公式为:
[0007]
[0008] 其中:t为时间,Td为接收信号相对发射信号的延迟,Pr为接收信号功率,ω0为线性调频的中心
角频率,fd为多普勒频率, 为发射信号的初相, 为接收信号的初相,NTzq为第N个调频周期;k——调频斜率, 为调频带宽,TM=5ms为调频周期;
[0009] 对经过A/D采样之后的雷达接收信号进行抽取滤波,将抽取滤波之后的数据进行解调,解调公式如下:
[0010]
[0011] 其中 为正调频解调信号, 为正调频接收信号;
[0013]
[0015]
[0016] 对解调之后的信号进行滤波,滤波后的信号:
[0017]
[0018]
[0019] 负调频抽取滤波和解调滤波处理与正调频处理方法一样,得到负调频解调滤波后的信号:
[0020]
[0021] 分别对正调频和负调频回波信号进行NFFT点FFT处理
[0022]
[0023]
[0024] 其中,测距捕获通道包括捕获CFAR和捕获回波中心计算;
[0025] 对于FFT处理后的信号采用捕获CFAR检测,检测到峰信号后,如果峰值信号超过检测门限,则判定存在回波信号;其中捕获CFAR检测门限定为NCFAR_n×Pn,4≤NCFAR_n≤10,Pn为噪声幅值,在FFT处理后的信号幅值上分段统计噪声,分别求均值后,将均值最小的噪声统计结果作为噪声幅值;
[0026] 捕获回波中心计算以FFT处理后的信号峰值为中心开一个窗口,检测超过门限的幅度值,从而计算回波的中心,其操作如下:
[0027] 以检测窗口峰值为中心,左右宽度为w1和w2进行检测,检测是否存在幅度超过门限hs的点,其中w1=wid×γ1,w2=wid×γ2,wid=fp×α,其中0.15≤α≤0.3,fp为超过门限的幅度值所对应的频率绝对值,如果为正斜率回波,则如果为负斜率回波,则
门限hs=Ap/Np, 其中Ap为峰值的幅度;
[0028] 将检测到的满足门限的信号中,提取幅度最大的10个信号,将其对应的频率与幅度做加权处理得到的
频谱位置,将其作为截取回波的中心fr,其中An和fn分别是各点对应的幅度和频谱位置;
[0029]
[0030] 将fr+和fr-分别作为正调频和负调频的波门中心;
[0031] 其中,捕获距离Rc计算公式如下:
[0032]
[0033] 其中fs为抽取后的信号采样率,c为光速。
[0034] 所述跟踪通道包括跟踪波门截取、跟踪CFAR、跟踪回波中心计算、主瓣截取和OCOG,从而得重心位置以提取测距信息:
[0035] 测距情况下跟踪通道与跟踪通道使能相连,如果跟踪尚未启动,但捕获通道已经得到有效的捕获距离,则在下一周期启动跟踪,以本周期的有效捕获距离Rc作为下一周期的跟踪递推数据Rtrack;
[0036] 按照上一周期的距离跟踪递推数据Rtrack和速度跟踪递推数据vtrack设置跟踪波门中心和宽度;
[0037] 负调频的跟踪波门中心
[0038]
[0040] 正调频的跟踪波门中心
[0041]
[0042]
[0043] 跟踪波门中心Wgz,跟踪波门宽度为上一周期的距离跟踪递推数据的Nbm倍,0.15≤Nbm≤0.3:
[0044]
[0045] 负调频截取时按照截取中心左侧 右侧 的比例截取;
[0046] 正调频截取时按照截取中心左侧 右侧 的比例截取;
[0047] 对于FFT处理后的信号采用跟踪CFAR检测,检测到峰信号后,如果峰值信号超过检测门限,则判定存在回波信号;其中跟踪CFAR检测门限定为NCFAR_n×Pn;
[0048] 跟踪回波中心计算以FFT处理后的信号峰值为中心开一个窗口,检测超过门限的幅度值,从而计算回波的中心,其操作如下:
[0049] 以检测窗口峰值为中心,左右宽度为w1和w2进行检测,检测是否存在幅度超过门限hs的点,其中w1=wid×γ1,w2=wid×γ2,wid=fp×α,如果为正斜率回波,则如果为负斜率回波,则
[0050] 将检测到的满足门限的信号中,提取幅度最大的10个信号,将其对应的频率与幅度做加权处理得到的频谱位置,将其作为截取回波的中心fr,其中An和fn分别是各点对应的幅度和频谱位置;
[0051]
[0052] 将 和 分别作为正调频和负调频的跟踪波门中心位置;
[0053] 主瓣截取是以跟踪CFAR输出的回波中心位置为
基础,在其上截取
指定宽度的回波,宽度Wzx:
[0054]
[0055] 其中,0.5≤Nzp≤0.9;
[0056] 负调频截取时按照截取中心左侧 右侧 的比例截取;
[0057] 正调频截取时按照截取中心左侧 右侧 的比例截取;
[0058] 在截取的主瓣宽度内,再进行二次判决截取主瓣,判决条件如下:
[0059] 主瓣内的信号幅值超过门限,所述门限为Nocog_n×Pn,2≤Nocog_n≤4;从波门中心往左右两边各自搜索,如果连续Nlx点信号幅值都没有超过门限,则主瓣下限/上限截止,其中6≤Nlx≤20;
[0060] 重心估计值S:
[0061]
[0062] 其中:
[0063] s[tr(m)]为二次截取后的FFT处理后的主瓣信号;
[0064] m为二次截取后的主瓣信号的FFT位置点。
[0065] 将OCOG中正调频重心作为 负调频回波中心作为
[0066] 跟踪距离Rt计算公式如下:
[0067]
[0068] 跟踪启动后,以本周期的有效跟踪距离Rt作为下一周期的跟踪递推数据Rtrack。
[0069] 本发明中,捕获通道与跟踪通道协同处理可以减少测量过程中野值的影响,因为着陆器下降过程,距离在一定的动态范围内,通过跟踪通道可以避免回波的毛刺的影响,使跟踪稳定在跟踪截取门限内。即使存在大动态的情况,通过捕获通道可以及时解算出测距值,与跟踪通道进行比较,可以及时退出跟踪,重新确定跟踪波门中心位置。这样,能够实现对地外天体面目标回波的有效测量,并且合理分配了
硬件计算资源,提高了处理速度和测距的
稳定性。
附图说明
[0070] 在以下参照附图对本发明的非限制性
实施例所做的详细描述中,本发明的其他特性和优点将会变得更明显,其中:
[0071] 图1示出根据本发明实施例的地外天体着陆测量雷达面目标回波信号测距处理方法的示意图。
[0072] 附图标记说明:
[0073] 1中频回波 2 AD采样
[0074] 3解调滤波抽取 4 FFT
[0075] 5预处理 6捕获CFAR
[0076] 7捕获回波中心计算 8捕获通道
[0077] 9跟踪通道使能 10跟踪波门截取
[0078] 11跟踪CFAR 12跟踪回波中心计算
[0079] 13主瓣截取 14 OCOG
[0080] 15波门中心及波门宽度 16跟踪通道
[0081] 17中心位置 18重心位置
[0082] 具体实例方式
[0083] 下面结合附图和实施例对本
申请做进一步的详细说明。可以理解的是,此处所描述的具体实施例仅用于解释本发明实施例,而非对本发明的限制。另外还需说明的是,为了便于描述,附图中仅示出与实施例密切相关的部分。实施例能够以多种形式实施,而不应该理解为限于在此说明的形式。提供这些实施例是为了使得本申请公开的内容更全面和完整,并将本发明的构思全面地传达给本领域技术人员。在图中相同的附图标记表示相同或相似的部分。
[0084] 此外,所描述的特征、结构或特性可以通过任何合适的方式结合在一个或多个实施例中。在下面的描述中,提供许多具体细节从而给出对实施例的充分理解。然而,本领域技术人员将意识到,可以实践本发明的技术方案而没有特定细节中的一个或更多,或者可以采用其他替代方式。可以理解,附图中所示的
框图不一定必须与物理上独立的实体相对应。
[0085] 针对地外天体着陆面目标回波展宽且存在毛刺的问题,本
发明人提出一种技术构思,对地外天体着陆面目标回波进行处理提取主瓣回波中心位置和重心位置。
[0086] 如图1所示,根据本发明的地外天体着陆测量雷达面目标回波信号测距处理方法包括以下步骤:
[0087] 对中频回波1进行预处理5,其包括A/D采样2、解调滤波抽取3和FFT4;
[0088] 所述测距捕获通道8包括捕获CFAR6和捕获回波中心计算7,从而得到回波主瓣中心位置17;
[0089] 所述测距跟踪通道16包括跟踪波门截取10、跟踪CFAR11、跟踪回波中心计算12、主瓣截取13和OCOG14,从而得重心位置18以提取测距信息。
[0090] 下面对本发明优选实施例的各步骤的执行方法做出详细说明:
[0091] 一、对中频回波1进行预处理5,其包括A/D采样2、解调滤波抽取3和FFT4。
[0092] 经过A/D采样2之后的雷达接收信号公式为:
[0093]
[0094] 其中:t为时间,Td为接收信号相对发射信号的延迟,Pr为接收信号功率,ω0为线性调频的中心角频率,fd为多普勒频率, 为发射信号的初相, 为接收信号的初相,NTzq为第N个调频周期;
[0095] k——调频斜率, 为调频带宽,TM=5ms为调频周期;
[0096] 对经过A/D采样2之后的雷达接收信号进行抽取滤波,将抽取滤波之后的数据进行解调,解调公式如下:
[0097]
[0098] 其中 为正调频解调信号, 为正调频接收信号;
[0099] 解调实部信号:
[0100]
[0101] 解调虚部信号:
[0102]
[0103] 对解调之后的信号进行滤波,滤波后的信号:
[0104]
[0105] 负调频抽取滤波和解调滤波处理与正调频处理方法一样,得到负调频解调滤波后的信号:
[0106]
[0107] 分别对正调频和负调频回波信号进行NFFT点FFT4处理
[0108]
[0109]
[0110] 第二步中,测距捕获通道8包括捕获CFAR6和捕获回波中心计算7的具体处理如下:
[0111] 对于FFT4处理后的信号采用捕获CFAR6检测,检测到峰信号后,如果峰值信号超过检测门限,则判定存在回波信号;其中捕获CFAR6检测门限定为NCFAR_n×Pn,4≤NCFAR_n≤10,Pn为噪声幅值,在FFT4处理后的信号幅值上分段统计噪声,分别求均值后,将均值最小的噪声统计结果作为噪声幅值;
[0112] 捕获回波中心计算7以FFT4处理后的信号峰值为中心开一个窗口,检测超过门限的幅度值,从而计算回波的中心,其操作如下:
[0113] 以检测窗口峰值为中心,左右宽度为w1和w2进行检测,检测是否存在幅度超过门限hs的点,其中w1=wid×γ1,w2=wid×γ2,wid=fp×α,其中0.15≤α≤0.3,fp为超过门限的幅度值所对应的频率绝对值,如果为正斜率回波,则如果为负斜率回波,则
门限hs=Ap/Np, 其中Ap为峰值的幅度;
[0114] 将检测到的满足门限的信号中,提取幅度最大的10个信号,将其对应的频率与幅度做加权处理得到的频谱位置,将其作为截取回波的中心fr,其中An和fn分别是各点对应的幅度和频谱位置;
[0115]
[0116] 将fr+和fr-分别作为正调频和负调频的波门中心15;
[0117] 其中,捕获距离Rc计算公式如下:
[0118]
[0119]
[0120] 其中fs为抽取后的信号采样率,c为光速。
[0121] 第三步中,跟踪通道16包括跟踪波门截取10、跟踪CFAR11、跟踪回波中心计算12、主瓣截取13和OCOG14,从而得重心位置18以提取测距信息具体包括:
[0122] 测距情况下跟踪通道16与跟踪通道使能9相连,如果跟踪尚未启动,但捕获通道8已经得到有效的捕获距离,则在下一周期启动跟踪,以本周期的有效捕获距离Rc作为下一周期的跟踪递推数据Rtrack;
[0123] 按照上一周期的距离跟踪递推数据Rtrack和速度跟踪递推数据vtrack设置跟踪波门中心和宽度15;
[0124] 负调频的跟踪波门中心
[0125]
[0126] 其中,λ为测速载波波长。
[0127] 正调频的跟踪波门中心
[0128]
[0129] 跟踪波门中心Wgz,跟踪波门宽度为上一周期的距离跟踪递推数据的Nbm倍,0.15≤Nbm≤0.3:
[0130]
[0131] 负调频截取时按照截取中心左侧 右侧 的比例截取;
[0132] 正调频截取时按照截取中心左侧 右侧 的比例截取;
[0133] 对于FFT4处理后的信号采用跟踪CFAR11检测,检测到峰信号后,如果峰值信号超过检测门限,则判定存在回波信号;其中跟踪CFAR11检测门限定为NCFAR_n×Pn;
[0134] 跟踪回波中心计算12以FFT4处理后的信号峰值为中心开一个窗口,检测超过门限的幅度值,从而计算回波的中心,其操作如下:
[0135] 以检测窗口峰值为中心,左右宽度为w1和w2进行检测,检测是否存在幅度超过门限hs的点,其中w1=wid×γ1,w2=wid×γ2,wid=fp×α,如果为正斜率回波,则如果为负斜率回波,则
[0136] 将检测到的满足门限的信号中,提取幅度最大的10个信号,将其对应的频率与幅度做加权处理得到的频谱位置,将其作为截取回波的中心fr,其中An和fn分别是各点对应的幅度和频谱位置;
[0137]
[0138] 将 和 分别作为正调频和负调频的跟踪波门中心位置;
[0139] 主瓣截取13是以跟踪CFAR11输出的回波中心位置为基础,在其上截取指定宽度的回波,宽度Wzx:
[0140]
[0141] 其中,0.5≤Nzp≤0.9;
[0142] 负调频截取时按照截取中心左侧 右侧 的比例截取;
[0143] 正调频截取时按照截取中心左侧 右侧 的比例截取;
[0144] 在截取的主瓣宽度内,再进行二次判决截取主瓣,判决条件如下:
[0145] 主瓣内的信号幅值超过门限,所述门限为Nocog_n×Pn,2≤Nocog_n≤4;从波门中心往左右两边各自搜索,如果连续Nlx点信号幅值都没有超过门限,则主瓣下限/上限截止,其中6≤Nlx≤20;
[0146] 重心估计值S:
[0147]
[0148] 其中:
[0149] s[tr(m)]为二次截取后的FFT4处理后的主瓣信号;
[0150] m为二次截取后的主瓣信号的FFT位置点。
[0151] 将OCOG中正调频重心作为 负调频回波中心作为
[0152] 跟踪距离Rt计算公式如下:
[0153]
[0154] 跟踪启动后,以本周期的有效跟踪距离Rt作为下一周期的跟踪递推数据Rtrack。
[0155] 综上,根据本发明实施例的技术方案采用捕获通道和跟踪通道同时处理方式,首先对中频信号进行预处理,包括AD采样、解调滤波抽取和FFT;然后捕获通道与跟踪通道同时进行处理。捕获通道包括捕获CFAR和捕获回波中心计算,从而得到中心位置;跟踪通道包括跟踪波门截取、跟踪CFAR、跟踪回波中心计算、主瓣截取和OCOG,从而得到回波主瓣的重心位置提取测距信息。捕获通道与跟踪通道协同处理可以减少测量过程中野值的影响,因为着陆器下降过程,距离在一定的动态范围内,通过跟踪通道可以避免回波的毛刺的影响,使跟踪稳定在跟踪截取门限内。即使存在大动态的情况,通过捕获通道可以及时解算出测距测速值,与跟踪通道进行比较,可以及时退出跟踪,重新确定跟踪波门中心位置。这样,能够实现对地外天体面目标回波的有效测量,并且合理分配了硬件计算资源,提高了处理速度和测距的稳定性。
[0156] 以上描述了本发明的技术构思及根据本发明技术构思的实施例。本领域技术人员在阅读
说明书及实践这里给出的实施例后,将容易想到本发明的其他实施方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化。这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未说明的本领域中的公知常识或者惯用技术手段。说明书和实施例仅为示例性的,本发明的保护范围由
权利要求限定。应当理解,本发明并不局限于上面已经描述及在附图中示出的内容,本领域技术人员可以在不脱离本申请公开的范围内进行各种
修改和变型。