一种频率步进雷达引信速度补偿方法

申请号 CN201310567505.1 申请日 2013-11-14 公开(公告)号 CN103558596A 公开(公告)日 2014-02-05
申请人 上海电机学院; 发明人 胡秀娟;
摘要 本 发明 提供了一种 频率 步进雷达引信速度补偿方法,包括:将 波形 熵与免疫科隆选择 算法 进行结合,以波形熵为搜索的亲和度函数,并以速度补偿准则为约束条件,采用免疫克隆选择算法完成搜索,以实现 频率步进雷达 引信一维距离像速度补偿。所述速度补偿准则包括对一次 相位 项进行补偿的最大速度变化单元、以及二次相位项相位变化不超过时一维距离像不失真条件。
权利要求

1.一种频率步进雷达引信速度补偿方法,其特征在于,将波形熵与免疫科隆选择算法进行结合,以波形熵为搜索的亲和度函数,并以速度补偿准则为约束条件,采用免疫克隆选择算法完成搜索,以实现频率步进雷达引信一维距离像速度补偿。
2.根据权利要求1所述的频率步进雷达引信速度补偿方法,其特征在于,所述速度补偿准则包括对一次相位项进行补偿的最大速度变化单元、以及二次相位项相位变化不超过时一维距离像不失真条件。
3.根据权利要求1或2所述的频率步进雷达引信速度补偿方法,其特征在于,采用免疫克隆选择算法完成搜索包括:
第一步骤,用于初始化个体,产生一个初始化个体;
第二步骤,用于初始化种群,由初始化个体,通过变异产生第一代种群;
第三步骤,用于计算各种群个体的亲和度,将其值从大到小排列,且将多个亲和度值最大的个体放入记忆库中;
第四步骤,用于从记忆库中取出经第三步骤排序后的多点采样值,对其进行克隆复制,亲和度值最大的个体克隆的个数也相应最多,然后根据亲和度值从大到小排序;
第五步骤,用于将克隆的个体进行变异,从而获得新的抗体群落;
第六步骤,用于计算个体之间的亲和度,在个体之间相似度超过阈值的个体中只保留亲和度值最高的一个个体,其余从种群中排除,然后重新排序,并用亲和度最高的多个个体替代原记忆库中的个体;
第七步骤,在满足引信实时性要求和速度补偿准则的指导下,判断是否满足速度补偿精度要求,如果不满足精度要求,则重复执行第二步骤至第六步骤;如果满足速度补偿精度要求,则搜索结束。
4.根据权利要求1或2所述的频率步进雷达引信速度补偿方法,其特征在于,所述频率步进雷达引信速度补偿方法用于具有多模战斗部的空地巡飞弹。
5.根据权利要求1或2所述的频率步进雷达引信速度补偿方法,其特征在于,所述频率步进雷达引信速度补偿方法用于装载频率步进雷达引信的导弹。

说明书全文

一种频率步进雷达引信速度补偿方法

技术领域

[0001] 本发明涉及雷达信号处理技术领域,更具体地说,本发明涉及一种频率步进雷达引信速度补偿方法。

背景技术

[0002] 频率步进雷达信号的速度补偿方法很多,主要方法有时域互相关法、频域互相关法、循环补偿法、最小波形熵法等。现有技术也提出了其他很多方法,但大多由上述方法衍生而来。
[0003] 时域互相关法是利用两组回波数据互相关函数,求出目标在两之间的走动距离,进而利用走动距离和目标速度之间的关系来估计目标速度的,在速度估计范围内,该方法能够获得较好的速度估计值,并且具有较好的抗噪声和杂波干扰。
[0004] 频域互相关法是基于两组回波数据互相关函数在零时刻的值建立目标速度估计模型,速度补偿精度较高。循环补偿法则是假定速度在给定的区间变化的条件下,对一维距离像进行速度补偿的一种方法,该方法的目标函数是信噪比SNR,在循环结束后,比较目标函数的值,选出该函数值最大的距离像作为输出,速度精确估计值,依赖于速度步进量。
[0005] 最小波形熵法是一种闭环迭代算法,将衡量随机变量不确定性的熵引申用于衡量一个信号的能量沿其参数轴的发散程度,并定义一个波形熵,在速度轴上搜索该波形熵的全局最小值对应的运动参数值,即目标的运动参数,估计精度随信噪比SNR的降低而明显降低。
[0006] 在上述速度补偿方法中,每种方法都有其固有的缺点,如时域互相关法受逆快速傅里叶变化长度的影响较大,虽然能通过补零来减小估计误差,但付出的代价却是增大系统的计算量;频域互相关法虽精度高,但不模糊测速范围小,仅适用于弹目相对运动速度小的交会情况;循环速度补偿法除受速度步进量大小的限制外,还受所需补偿速度值是否为速度步进量的整数倍的影响;最小波形熵法计算量大,实时性较差,且存在局部最小值。

发明内容

[0007] 从频率步进脉冲高分辨处理的度来看,要用到不同周期子脉冲的相位信息,而目标在各个周期之间的运动会使子脉冲回波的相位发生变化。相位的变化对子脉冲本身没有影响,但是在后面的合成处理中,它会破坏脉冲序列的相位关系,造成一维距离像输出产生误差。从频率步进波形模糊特性上来看,也需要解决这一问题。“距离——速度”耦合是频率步进脉冲体制不可回避的问题,因此弹目之间径向相对速度所产生的多普勒效应对目标一维距离像有非常大的影响,即波形发散,峰值降低,其程度与弹目相对速度大小有关,相对速度越大波形发散越严重,峰值降低越多,从而使距离分辨率下降越大;峰值位置产生时移,弹目相对速度越大,产生的时移越大;模糊函数时间轴切割图形与sinc函数相比进一步失真,而且主峰加宽,其程度均与弹目相对速度大小有关,速度越大,失真越严重;使信噪比降低,影响检测效果。
[0008] 因此,本发明旨在解决由于弹目径向相对速度引起的上述问题,而解决问题的关键就是速度补偿,即利用目标波形的最小波形熵进行速度补偿,速度补偿准则为约束条件,并将波形熵与免疫科隆选择算法进行结合来实现频率步进雷达引信的速度补偿。
[0009] 具体地说,对于频率步进脉冲信号,弹目相对速度越大,主瓣脉冲展宽、旁瓣电平抬高、波形趋于平坦,波形熵越大,弹目相对速度越小,波形锐化度越高,波形熵越小。当弹目径向相对速度补偿误差为零时,波形熵达到最小,此时得到的速度估计值为最佳速度补偿值。但最小波形熵法是一种迭代算法,其计算量大,实时性较差,且在局部区域进行。因此,在本发明中,将波形熵与免疫科隆选择算法有效结合,以波形熵为搜索的亲和度函数,并以速度补偿准则为约束条件,采用免疫克隆选择算法完成搜索,可减少计算量,加快搜索速度,从而满足引信实时性要求、引战良好配合的要求、及对目标实现精确打击的要求。
[0010] 本发明提供了一种频率步进雷达引信速度补偿方法,其中,将波形熵与免疫科隆选择算法进行结合,以波形熵为搜索的亲和度函数,并以速度补偿准则为约束条件,采用免疫克隆选择算法完成搜索,以实现频率步进雷达引信一维距离像速度补偿。
[0011] 优选地,所述速度补偿准则包括对一次相位项进行补偿的最大速度变化单元、以及二次相位项相位变化不超过时一维距离像不失真条件。
[0012] 优选地,采用免疫克隆选择算法完成搜索包括:
[0013] 第一步骤,用于初始化个体,产生一个初始化个体;
[0014] 第二步骤,用于初始化种群,由初始化个体,通过变异产生第一代种群;
[0015] 第三步骤,用于计算各种群个体的亲和度,将其值从大到小排列,且将多个亲和度值最大的个体放入记忆库中;
[0016] 第四步骤,用于从记忆库中取出经第三步骤排序后的多点采样值,对其进行克隆复制,亲和度值最大的个体克隆的个数也相应最多,然后根据亲和度值从大到小排序;
[0017] 第五步骤,用于将克隆的个体进行变异,从而获得新的抗体群落;
[0018] 第六步骤,用于计算个体之间的亲和度,在个体之间相似度超过阈值的个体中只保留亲和度值最高的一个个体,其余从种群中排除,然后重新排序,并用亲和度最高的多个个体替代原记忆库中的个体;
[0019] 第七步骤,在满足引信实时性要求和速度补偿准则的指导下,判断是否满足速度补偿精度要求,如果不满足精度要求,则重复执行第二步骤至第六步骤;如果满足速度补偿精度要求,则搜索结束。
[0020] 优选地,所述频率步进雷达引信速度补偿方法用于具有多模战斗部的空地巡飞弹。
[0021] 优选地,所述频率步进雷达引信速度补偿方法用于装载频率步进雷达引信的导弹。附图说明
[0022] 结合附图,并通过参考下面的详细描述,将会更容易地对本发明有更完整的理解并且更容易地理解其伴随的优点和特征,其中:
[0023] 图1示意性地示出了根据本发明优选实施例的采用免疫克隆选择算法完成搜索的步骤。
[0024] 需要说明的是,附图用于说明本发明,而非限制本发明。注意,表示结构的附图可能并非按比例绘制。并且,附图中,相同或者类似的元件标有相同或者类似的标号。

具体实施方式

[0025] 为了使本发明的内容更加清楚和易懂,下面结合具体实施例和附图对本发明的内容进行详细描述。
[0026] 为实现频率步进雷达引信一维距离像的速度补偿,并满足补偿精度和实时性要求,本发明优选实施例提出一种速度补偿方法,其中可采取下述手段之一或者下述手段的组合:
[0027] (1)结合克隆选择、免疫、最小波形熵思想,以最小波形熵为亲和函数,速度补偿准则为约束条件,免疫科隆选择进行局部、全局搜索,实现频率步进雷达引信一维距离像速度补偿;
[0028] (2)在结构上采用“初始种群产生——亲和度函数评价——种群个体信息交换——新种群产生”的循环过程,以寻求速度估计的最优解,且较好的群体多样性使群体中的某些个体能够跳出局部最优,更有效地达到全局最优解;
[0029] (3)设置了速度补偿准则,即设置了对一次相位项进行补偿的最大速度变化单元、以及二次相位项相位变化不超过π/2时一维距离像不失真条件;
[0030] (4)加入了克隆选择因子,对速度信息集合进行选取,可以得到充分体现速度信息的最小集合,减少了运算量,同时搜索解空间中的一系列的点,这使得它的搜索速度大大加快;
[0031] (5)该速度补偿方法是通过回波采样空间的扩张和压缩,将局部搜索和全局搜索结合起来实现问题的求解,即将一个低维的问题转化到更高维空间(扩张)去解决,然后将结果投影到低维空间(压缩)中,从而获得对问题更全面的认识,克服了传统方法搜索范围小的缺点;
[0032] (6)引入了克隆选择扩增、克隆删除、克隆变异等思想,这不仅使亲合度高的优秀抗体能够获得较大的克隆规模,提高了速度补偿精度,而且能加快方法的收敛速度,更好地满足引信实时性要求;
[0033] <具体实施方式>
[0034] 本发明具体实施方式以具有多模战斗部的空地巡飞弹为应用背景,该巡飞弹的战斗部是一种多用途/可选择EFP战斗部,这种战斗部可以根据攻击目标类型的不同,通过引信系统起爆方式的选择,分别形成攻击重型装甲目标的长延展杆、攻击轻型装甲目标的空气稳定的单个EFP和攻击中软目标(如雷达系统、器材、人员等)的多个EFP。因此,针对不同的攻击目标,要求引信对目标进行分类识别,故而选择不同的起爆方式。
[0035] 为实现引信对目标的分类识别,采用步进频雷达引信对目标的一维距离成像。但是,由于弹目之间存在径向相对速度,严重影响了成像质量,需要对其进行速度补偿,同时考虑引信对实时性的严格要求。为了精确补偿弹目径向相对速度,同时减少计算量,提高估计速度,本发明以波形熵为搜索的亲和度函数,采用免疫克隆选择算法来完成搜索。
[0036] 假 设 离 散 化 信 号 波 形 序 列 为X={x1, …,xN},令 ,借用统计学中概率的含义,把一维距离像中每个点组成的集合看作一个样本空间,pl则表示距离像中每一点出现的概率,则X的熵定义为定义搜索弹目径向运动参数的免疫克隆算法的亲和度函
数为f(dv,da)=-E(dv,da),其中,E(dv,da)是目标回波信号的波形熵,它是速度误差估计dv和加速度误差估计da的函数。
[0037] 图1示意性地示出了根据本发明优选实施例的采用免疫克隆选择算法完成搜索的步骤。
[0038] 如图1所示,采用免疫克隆选择算法完成搜索的步骤(免疫克隆选择算法的步骤)包括:
[0039] 第一步骤S1,用于初始化个体:假设有N个采样值,当k=0(k=0,1,2,…N-1),按照输入码书的索引值,产生一个初始化个体A(0);其中,N为适当正整数。
[0040] 第二步骤S2,用于初始化种群:由初始化个体A(0)通过变异产生第一代种群A1(0),A2(0),…,AM(0),其中M为种群规模;
[0041] 第三步骤S3,用于计算各种群个体的亲和度,将其值从大到小排列,且将N个亲和度值最大的个体放入记忆库中;
[0042] 第四步骤S4,用于从记忆库中取出经第三步骤S3排序后的N点采样值,对其进行克隆复制,亲和度值最大的个体克隆的个数也相应最多,然后根据亲和度值从大到小排序;
[0043] 第五步骤S5,用于将克隆的个体进行变异,从而获得新的抗体群落A(k);
[0044] 第六步骤S6,用于计算个体之间的亲和度,在个体之间相似度超过阈值Tr的个体中只保留亲和度值最高的一个个体(克隆选择),其余从种群中排除(克隆抑制),然后重新排序,并用亲和度最高的N个个体替代原记忆库中的个体(亚动过程);
[0045] 第七步骤S7,在满足引信实时性要求和速度补偿准则的指导下,判断是否满足速度补偿精度要求,如果不满足精度要求,则重复执行第二步骤S2至第六步骤S5;如果满足速度补偿精度要求,则搜索结束。
[0046] 例如,在具体实施例中,满足速度补偿精度要求的速度补偿准则可以是如下所述:
[0047] 假定弹目径向相对速度为vt,则回波走动的高分辨距离单元数为2vtNTr/c。一次相位变化主要产生耦合距离误差,其大小为 相应的回波延时误差为τe,即 当对采样序列进行补偿
时,若以距离走动不超过半个高分辨距离单元为条件,则 其中,
可见,为获得半个距离单
元的测距误差,测速误差Δvt1需满足 上式即为对一次相位项进行
补偿的最大速度变化单元。若以相参处理间隔NTr内二次相位项相位变化不超过π/2为合成像不失真条件,则测速误差Δvt2应满足关系 即
[0048] 本发明的频率步进雷达引信速度补偿方法可应用于具有多模战斗部的空地巡飞弹,该巡飞弹的战斗部是一种多用途/可选择爆炸成型弹丸(EFP)战斗部,这种战斗部可以根据攻击目标类型的不同,通过引信系统起爆方式的选择,分别形成攻击重型装甲目标的长延展杆、攻击轻型装甲目标的空气稳定的单个EFP和攻击中软目标的多个EFP。针对不同的攻击目标,频率步进雷达引信对目标进行分类识别,故而选择不同的起爆方式。除此之外,本发明还可应用于任何装载频率步进雷达引信的导弹,来处理由于弹目径向相对速度引起的雷达引信目标一维距离像的时移、散焦问题。
[0049] 上述应用可有效地提高频率步进雷达引信的定距精度和抗干扰性,从而使引信可选择最佳攻击时空来准确摧毁目标,实现引战配合,提高引信全天时、全天候作战能力。
[0050] 本发明具有以下积极效果:
[0051] (1)在背景噪声为零均值高斯白噪声、信噪比SNR=20dB、杂波为韦伯分布、信杂比SCR=20dB、目标具有10个强散射点、弹目径向相对速度为302m/s、逆快速傅里叶变换长度为512的情况下,利用本方法得到的速度估计值与真值十分接近,估计误差为0.1m/s;
[0052] (2)逆快速傅里叶变换长度为512的情况下,本方法在MATLAB下的运行所需时间与循环补偿法接近,在特定情况下,可满足引信的实时性要求;
[0053] (3)方法具有可移植性,通用性,该方法不仅可用于实时性要求较高的频率步进雷达引信,还可用于实时性要求并不是很苛刻的频率步进雷达制导和搜索跟踪中。
[0054] 此外,需要说明的是,除非特别说明或者指出,否则说明书中的术语“第一”、“第二”、“第三”等描述仅仅用于区分说明书中的各个组件、元素、步骤等,而不是用于表示各个组件、元素、步骤之间的逻辑关系或者顺序关系等。
[0055] 可以理解的是,虽然本发明已以较佳实施例披露如上,然而上述实施例并非用以限定本发明。对于任何熟悉本领域的技术人员而言,在不脱离本发明技术方案范围情况下,都可利用上述揭示的技术内容对本发明技术方案作出许多可能的变动和修饰,或修改为等同变化的等效实施例。因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化及修饰,均仍属于本发明技术方案保护的范围内。
QQ群二维码
意见反馈