基于分数阶傅里叶变换的单站无源快速定位方法技术

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

【技术实现步骤摘要】

本专利技术属于信号处理
，特别涉及一种无源定位方法，可用于单站侦察。
技术介绍
无源定位是指单个侦察设备或者多个侦察设备在探测到散射源或者辐射源并获得有关定位参数的基础上，利用适当的数据处理手段，确定对散射源或辐射源在三维空间中的位置。无源定位可以简单分为多站定位和单站定位。多站定位通过部署在不同地理位置的站点接收到的目标辐射信号进行信号级别或者数据级别的融合，从而达到定位的目的。在通信设备定位中应用较广；而雷达发射天线往往波束较窄，很难保证多个接收站能接收到目标信号，同时多站的时间同步机制要求设备复杂，而单站定位由于几乎没有上述问题，因而得到了广泛的关注。单站无源定位通过单一观测平台被动接收辐射源信号完成目标定位，是获取目标位置信息的重要手段，已经应用到重要目标位置监视、紧急营救、智能交通以及反恐维稳等领域。单站定位与多站定位相比最大特点是只有一个观测平台，设备简单、部署灵活、机动性能强，因此在星载、机载、舰载等多种运动平台目标定位方面发挥了重要作用。目前单站定位方法的研究，包括测向定位、到达时间差定位、多普勒频率及其变化率定位、相位差及其变化率定位以及上述的组合定位法等。其中，多普勒频率变化率定位和相位差变化率定位较为常用，这两种方法为了提高目标定位跟踪能力，采用了如包括牛顿迭代方法、卡尔曼滤波及其扩展方法、贝叶斯估计理论的粒子滤波方法等。2004年开始，意大利在直升机上采用双天线实现多普勒变化率定位实验。该方法在飞机两侧放置2个干涉仪天线，利用其频率差测量结果定位。实验结果显示该方法的最佳位置在两侧斜视角±45度处，其他位置均将下降；在定位距离约5km，定位35秒时间条件下，目标定位误差将超过10％。上述的单站定位方法都是利用了目标与侦察设备位置变化的非线性变化特性，采用了非相干的积累方法实现目标位置的估计。但是随着现代雷达技术的发展，雷达的射频源稳定度得到了明显提高；雷达在成像、空时处理等模式下，都会有较长的信号稳定期，侦察获取信号的相干性明显增强。常规的单站定位方法不能有效的利用信号的相关性，需要很长的积累时间。
技术实现思路
本专利技术的目的在于针对上述已有技术的不足，提出一种基于分数阶傅里叶变换的单站无源快速定位方法，以减小运算量，改善定位精度，实现快速定位。本专利技术技术思路是，将合成孔径成像雷达的相关技术应用到侦察领域，利用辐射源的相干性，将长时间积累信号虚拟合成长的天线阵列，使其不满足远场条件，即天线阵列的二次相位无法忽略；利用分数阶傅里叶变换对该不可忽略的二次相位进行估计，进而利用二次相位与辐射源距离的反比关系实现对辐射源目标位置的估计。其实现步骤包括如下：(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)分数阶傅里叶变换的阶次α在[0,π/2]范围内等间隔的选择m个值：α1，α2…αk…αm，k＝1，2，3…m；(5b)对峰值采样复信号矩阵M分别进行α1，α2…αk…αm阶的分数阶傅里叶变换，记录分数阶傅里叶变换结果的峰值A1，A2…Ak…Am；(5c)比较分数阶傅里叶变换的峰值A1，A2…Ak…Am，利用最大的峰值Amax所对应的参数αmax计算最优调频斜率μ^2=tan(αmax);]]>(6)根据最优调频率获得补偿后的信号矩阵Mc：(6a)利用最优调频率构建补偿信号矩阵H＝[h(η1),h(η2)…h(ηi)…h(ηn)]，其中，h(ηi)为补偿信号，ηi＝iTpri，ηi表示辐射源脉冲到达时间，Tpri表示辐射源脉冲重复周期；(6b)将补偿信号矩阵H与峰值采样复信号矩阵M进行点乘，消除峰值采样复信号关于时间的二次相位，得到补偿后的信号矩阵：Mc＝[uu1(t)h(η1),…uui(t)h(ηi)…uun(t)h(ηn)]；(7)对补偿后的信号矩阵Mc进行傅里叶变换得到多普勒频率计算高精度的测向结果(8)利用高精度的测向结果计算辐射源和侦察设备的距离(9)结合高精度的测向结果和辐射源和侦察设备的距离得到辐射源位置。本专利技术与现有技术相比具有以下优点：第一，相对于传统的定位方法，本专利技术利用了辐射源的相干性，缩短定位时间，可以实现快速定位；第二，将合成孔径雷达的原理应用在雷达定位中，将长时间积累的信号虚拟合成长的天线阵列，实现对目标的精确角度测量；第三，利用虚拟长天线阵列不满足远场条件，即接收信号关于时间的二次相位无法忽略的特性，通过分数阶傅里叶变换进行时频域匹配搜索，得到距离参数，提高定位精度。附图说明图1是本专利技术的实现流程图；图2是本专利技术的定位误差图；图3是本专利技术方法定位时间与相位差分法定位时间的对比图。具体实施方式下面结合附图对本专利技术做进一步详细的描述。参照图1，本专利技术的具体实现步骤如下：步骤1、获取基带信号。本专利技术采用做匀速直线运动的侦察设备，对接收信号进行瞬时测频和测向，得到测频结果和测向结果并利用测频结果进行混频消除载频，得到基带信号ur(t)，测频误差要求较低，小于1MHz即可。步骤2、积累基带信号。对基带信号ur(t)进行积累，即每进行一次混频得到的基带信号进行储存，累计得到一组基带信号：ur1(t),ur2(t),…uri(t)…urn(t)，i＝1,2,3…n其中，n表示积累脉冲总个数。步骤3、对积累的基带信号进行匹配滤波。在积累的基带信号中选择其中的任意一个作为参考信号ur0(t)，将该参考信号ur0(t)与所有积累的信号ur1(t),ur2(t),…uri(t)…urn(t)分别做卷积，得到匹配滤本文档来自技高网...

【技术保护点】
一种基于分数阶傅里叶变换的单站无源快速定位方法，包括如下步骤：(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)分数阶傅里叶变换的阶次α在[0,π/2]范围内等间隔的选择m个值：α1,α2…αk…αm,k＝1,2,3…m；(5b)对峰值采样复信号矩阵M分别进行α1,α2…αk…αm阶的分数阶傅里叶变换，记录分数阶傅里叶变换结果的峰值A1,A2…Ak…Am；(5c)比较分数阶傅里叶变换的峰值A1,A2…Ak…Am，利用最大的峰值Amax所对应的参数αmax计算最优调频斜率μ^2=tan(αmax);]]>(6)根据最优调频率获得补偿后的信号矩阵Mc：(6a)利用最优调频率构建补偿信号矩阵H＝[h(η1),h(η2)…h(ηi)…h(ηn)]，其中，h(ηi)为补偿信号，ηi＝iTpri，ηi表示辐射源脉冲到达时间，Tpri表示辐射源脉冲重复周期；(6b)将补偿信号矩阵H与峰值采样复信号矩阵M进行点乘，消除峰值采样复信号关于时间的二次相位，得到补偿后的信号矩阵：Mc＝[uu1(t)h(η1),…uui(t)h(ηi)…uun(t)h(ηn)]；(7)对补偿后的信号矩阵Mc进行傅里叶变换得到多普勒频率计算高精度的测向结果(8)利用高精度的测向结果计算辐射源和侦察设备的距离(9)结合高精度的测向结果和辐射源和侦察设备的距离得到辐射源位置。...

【技术特征摘要】
1.一种基于分数阶傅里叶变换的单站无源快速定位方法，包括如下步骤：
(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)分数阶傅里叶变换的阶次α在[0,π/2]范围内等间隔的选择m个值：
α1,α2…αk…αm,k＝1,2,3…m；
(5b)对峰值采样复信号矩阵M分别进行α1,α2…αk…αm阶的分数阶傅里叶变换，
记录分数阶傅里叶变换结果的峰值A1,A2…Ak…Am；
(5c)比较分数阶傅里叶变换的峰值A1,A2…Ak…Am，利用最大的峰值Amax所对应
的参数αmax计算最优调频斜率μ^2=tan(αmax);]]>(6)根据最优调频率获得补偿后的信号矩阵Mc：
(6a)利用最优调频率构建补偿信号矩阵H＝[h(η1),h(η2)…h(ηi)…h(ηn)]，其中，
h(ηi)为补偿信号，ηi＝iTpri，ηi表示辐射源脉冲到达时间，Tpri表示辐射源脉冲重复周
期；
(6b)将补偿信号矩阵H与峰值采样复信号矩阵M进行点乘，消除峰值采样复信
号关于时间的二次相位，得到补偿后的信号矩阵：
Mc＝[uu1(t)h(η...

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