基于最小二乘法的单站无源快速定位方法技术

本发明专利技术公开了一种基于最小二乘法的单站无源快速定位方法，主要解决现有技术对辐射源相关性利用不足、定位精度较低、定位时间长的问题。其实现步骤是：1获取基带信号并对其进行积累，得到基带信号串；2、对基带信号串进行匹配滤波；3、对匹配滤波结果进行多倍插值，提取峰值采样复信号；4、提取峰值采样复信号的相位信息进行最小二乘估计，获得最优调频率；5、利用最优调频率峰值采样复信号进行相位补偿；6、对补偿后的信号进行傅里叶变换，得到高精度测向结果，根据高精度测向结果，计算辐射源距离，进而得到辐射源位置。本发明专利技术将合成孔径雷达的概念应用到侦察定位方面，具有定位精度高，定位时间短的优点，可用于目标侦察和干扰源定位。

【技术实现步骤摘要】

本专利技术属于信号处理
，特别涉及一种无源定位方法，可用于单站侦察。
技术介绍
无源定位是指由单个侦察设备或者多个侦察设备在探测到散射源或者辐射源并 获得有关定位参数的基础上，利用适当的数据处理手段，确定出散射源或辐射源在三维空 间中的位置。 无源定位可以简单分为多站定位和单站定位。多站定位通过部署在不同地理位置 的站点接收到的目标辐射信号进行信号级别或者数据级别的融合，从而达到定位的目的。 在通信设备定位中应用较广；而雷达发射天线往往波束较窄，很难保证多个接收站能接收 到目标信号，同时多站的时间同步机制要求设备复杂，而单站定位由于几乎没有上述问题， 因而得到了广泛的关注。 单站无源定位通过单一观测平台被动接收辐射源信号完成目标定位，是获取目标 位置信息的重要手段，已经应用到重要目标位置监视、紧急营救、智能交通以及反恐维稳等 领域。单站定位与多站定位相比最大特点是只有一个观测平台，设备简单、部署灵活、机动 性能强，因此在星载、机载、舰载等多种运动平台目标定位方面发挥了重要作用。 目前单站定位方法的研究，包括测向定位、到达时间差定位、多普勒频率及其变化 率定位、相位差及其变化率定位以及上述的组合定位法等。其中，多普勒频率变化率定位和 相位差变化率定位较为常用，这两种方法为了提高目标定位跟踪能力，采用了如包括牛顿 迭代方法、卡尔曼滤波及其扩展方法、贝叶斯估计理论的粒子滤波方法等。2004年开始，意 大利在直升机上采用双天线实现多普勒变化率定位实验。该方法在飞机两侧放置2个干涉 仪天线，利用其频率差测量结果定位。实验结果显示该方法的最佳位置在两侧斜视角±45 度处，其他位置均将下降；在定位距离约5km，定位35秒时间条件下，目标定位误差将超过 10%〇上述的单站定位方法都是利用了目标与侦察设备位置变化的非线性变化特性，采 用了非相干的积累方法实现目标位置的估计。但是随着现代雷达技术的发展，雷达的射频 源稳定度得到了明显提高;雷达在成像、空时处理等模式下，都会有较长的信号稳定期，侦 察获取信号的相干性明显增强。常规的单站定位方法不能有效的利用信号的相关性，需要 很长的积累时间。
技术实现思路
本专利技术的目的在于针对上述已有技术的不足，提出一种基于最小二乘法的单站无 源快速定位方法，以减小运算量，改善定位精度，实现快速定位。 本专利技术技术思路是，将合成孔径成像雷达的相关技术应用到侦察领域，利用辐射 源的相干性，将长时间积累信号虚拟合成长的天线阵列，使其不满足远场条件，即天线阵列 的二次相位无法忽略；利用最小二乘法对该不可忽略的二次相位进行估计，进而利用二次 相位与辐射源距离的反比关系实现对辐射源目标位置的估计。其实现步骤包括如下： (1)侦察设备做匀速直线运动，对接收信号进行瞬时测频和测向，得到测频结果/ 和测向结果#，.并利用测频结果/进行混频消除载频，得到基带信号Ur(t); (2)对基带信号Ur(t)进行积累，得到一组积累的基带信号： Url(t) ,Ur2(t) , 1''Uri ( t) ''Urn (t ) , ? = 1,2,3 ·''Π ,其中，n表示积累脉冲总个数； (3)对积累的基带信号进行匹配滤波： (3a)在积累的基带信号中选择其中一个作为参考信号urQ(t); (3b)将参考信号urQ(t)与所有积累脉冲信号进行卷积，得到一组卷积结果m(t)， U2(t) , ???Ui(t)-'-Un(t)； (4)对每个卷积结果m(t)进行多倍插值，得到采样信号^(〇,提取插值结果Vl(t) 峰值处的采样复信号um(t)，得到峰值采样复信号矩阵： M= ； (5)利用峰值采样复信号矩阵M的相位信息进行最小二乘估计，得到最优调频率 (5a)利用测向结果&构建多普勒补偿信号； 其中，ru = iTpri，ru表示辐射源脉冲到达时间，Tpri表示辐射源脉冲重复周期； (5b)用η个多普勒补偿信号g(m)，g(n2)"_g(ni)…g(nn)，构建多普勒补偿信号矩 阵 G=; (5c)将多普勒补偿信号矩阵G与峰值采样复信号矩阵Μ进行点乘，获得乘积矩阵 Mg，提取乘积矩阵Mg的相位信息，得到相位矩阵： Φ = ? (5d)利用相位矩阵Φ的元素和辐射源脉冲到达时间ru，按照下式估计出最优调频 斜率Λ: (6)根据最优调频率Λ获得补偿后的信号矩阵M。： (6a)利用最优调频率爲构建补偿信号矩阵H=，其 中，h(ru)为补偿信号； (6b)将补偿信号矩阵Η与峰值采样复信号矩阵Μ进行点乘，消除峰值采样复信号关 于时间的二次相位，得到补偿后的信号矩阵： Mc= ; (7)对补偿后的信号矩阵Μ。进行傅里叶变换得到多普勒频率$，计算高精度的测 向结果邊_; (8)利用高精度的测向结果4^，计算辐射源和侦察设备的距离h (9)结合高精度的测向结果^_和辐射源和侦察设备的距离卩，得到辐射源位置。 本专利技术与现有技术相比具有以下优点： 第一，相对于传统的定位方法，本专利技术利用了辐射源的相干性，缩短定位时间，可 以实现快速定位； 第二，将合成孔径雷达的原理应用在雷达定位中，将长时间积累的信号虚拟合成 长的天线阵列，实现对目标的精确角度测量； 第三，利用虚拟长天线阵列不满足远场条件，即接收信号关于时间的二次相位无 法忽略的特性，利用最小二乘法得到距离参数，提高定位精度。【附图说明】图1是本专利技术的实现流程图； 图2是本专利技术的定位误差图； 图3是本专利技术方法的定位时间与相位差分法定位时间的对比图。【具体实施方式】 下面结合附图对本专利技术做进一步详细的描述。 参照图1，本专利技术的具体实现步骤如下：步骤1、获取基带信号。 本专利技术采用做匀速直线运动的侦察设备，对接收信号进行瞬时测频和测向，得到 测频结果/和测向结果A并利用测频结果/进行混频消除载频，得到基带信号u r(t)，测频 误差要求较低，小于1MHz即可。 步骤2、积累基带信号。 对基带信号ur(t)进行积累，即每进行一次混频得到的基带信号进行储存，累计得 到一组基带信号： Url(t) ,Ur2(t) , 1''Uri ( t) ''Urn (t ) , ? = 1,2,3 ·''Π 其中，n表示积累脉冲总个数。 步骤3、对积累的基带信号进行匹配滤波。 在积累的基带信号中选择其中的任意一个作为参考信号urQ(t)，将该参考信号ur0 ⑴与所有积累的信号Url⑴，1!:2(〇，~11:办卜_11" 1(〇分别做卷积，得到匹配滤波结果： Ul(t) ,U2(t), ???Ui(t)---Un(t)步骤4、多倍插值，提取峰值采样复信号。 (4a)对每个卷积结果m(t)分别进行傅里叶变换，将卷积结果m(t)转换到频域，并 在频域中心插入(K-l)Nf个零，得到插值后的频域信号，其中，Nf为傅里叶变换后的数据个 数，K为插值倍数，取值为2的整数次幂，本实例K取8; (4b)对插值后的频域信号进行逆傅里叶变换，将信号恢复到时域，得到采样信号 Vl(t);记录采样结果Vl(t)达到峰值时的时间tmax，提取时间t max对应的峰值采样复信号um (t) - V i ( tmax )； (4c)对所有的卷积结果ui(t)，U2(t)，…Ui(t)'"un(t)进行(4本文档来自技高网...

【技术保护点】
一种基于最小二乘法的单站无源快速定位方法，包括如下步骤：(1)侦察设备做匀速直线运动，对接收信号进行瞬时测频和测向，得到测频结果和测向结果并利用测频结果进行混频消除载频，得到基带信号ur(t)；(2)对基带信号ur(t)进行积累，得到一组积累的基带信号：ur1(t),ur2(t),…uri(t)…urn(t)，i＝1,2,3…n，其中，n表示积累脉冲总个数；(3)对积累的基带信号进行匹配滤波：(3a)在积累的基带信号中选择其中一个作为参考信号ur0(t)；(3b)将参考信号ur0(t)与所有积累脉冲信号进行卷积，得到一组卷积结果u1(t),u2(t),…ui(t)…un(t)；(4)对每个卷积结果ui(t)进行多倍插值，得到采样信号vi(t)，提取插值结果vi(t)峰值处的采样复信号uui(t)，得到峰值采样复信号矩阵：M＝[uu1(t),uu2(t),…uui(t)…uun(t)]；(5)利用峰值采样复信号矩阵M的相位信息进行最小二乘估计，得到最优调频率(5a)利用测向结果构建多普勒补偿信号；g(ηi)=exp(-j2πvλcosθ^ηi),]]>其中，ηi＝iTpri，ηi表示辐射源脉冲到达时间，Tpri表示辐射源脉冲重复周期；(5b)用n个多普勒补偿信号g(η1),g(η2)…g(ηi)…g(ηn)，构建多普勒补偿信号矩阵G＝[g(η1),g(η2)…g(ηi)…g(ηn)]；(5c)将多普勒补偿信号矩阵G与峰值采样复信号矩阵M进行点乘，获得乘积矩阵MG，提取乘积矩阵MG的相位信息，得到相位矩阵：(5d)利用相位矩阵Φ的元素和辐射源脉冲到达时间ηi，按照下式估计出最优调频斜率(6)根据最优调频率获得补偿后的信号矩阵Mc：(6a)利用最优调频率构建补偿信号矩阵H＝[h(η1),h(η2)…h(ηi)…h(ηn)]，其中，h(ηi)为补偿信号；(6b)将补偿信号矩阵H与峰值采样复信号矩阵M进行点乘，消除峰值采样复信号关于时间的二次相位，得到补偿后的信号矩阵：Mc＝[uu1(t)h(η1),…uui(t)h(ηi)…uun(t)h(ηn)]；(7)对补偿后的信号矩阵Mc进行傅里叶变换得到多普勒频率计算高精度的测向结果(8)利用高精度的测向结果计算辐射源和侦察设备的距离(9)结合高精度的测向结果和辐射源和侦察设备的距离得到辐射源位置。...

【技术特征摘要】

【专利技术属性】
技术研发人员：刘高高，鲍丹，朱光前，武斌，秦国栋，蔡晶晶，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61