一种基于投影策略的分布式多传感器多目标无源定位方法 |
|||||||
申请号 | CN201410454491.7 | 申请日 | 2014-09-09 | 公开(公告)号 | CN104198987A | 公开(公告)日 | 2014-12-10 |
申请人 | 乐山师范学院; | 发明人 | 樊玲; 周昌海; | ||||
摘要 | 本 发明 公开了一种基于投影策略的分布式多 传感器 多目标无源 定位 方法,涉及无源定位技术领域。本方法将传感器获得的观测量通过距离差方程投影到投影空间,在投影空间中积累不同传感器获得的目标回波的幅度值。通过对投影空间积累值采用 门 限检测判断目标的有无及目标的 位置 ,实现对多个目标的检测和定位。本方法解决了量测数据关联的难题,计算复杂度低、实时性强。 | ||||||
权利要求 | 1.一种基于投影策略的分布式多传感器多目标无源定位方法,其特征在于:其步骤如下: |
||||||
说明书全文 | 一种基于投影策略的分布式多传感器多目标无源定位方法技术领域[0001] 本发明属于无源定位技术领域,更具体地说是涉及一种基于投影策略的分布式多传感器多目标无源定位方法。 背景技术[0002] 分布式多传感器无源定位利用多个传感器接收来自目标反射的信号进行定位。由于传感器自身不向外辐射电磁波,因此具有良好的隐蔽性、抗电磁干扰和电磁隐身、低空探测能力强等优点,在军事侦察、监视等领域具有广泛的应用前景。 [0003] 现有的分布式无源定位方法主要有:测向交叉定位法、到达时间差定位法、多普勒频差定位法等。 [0004] 到达时间差定位法具有较强的工程可实现性和较高的定位精度,因此被广泛应用于雷达、声纳及辐射源定位系统中。 [0005] 在单目标情况下,到达时间差定位法的原理如下。 [0006] 如图1所示,设定三维空间中任意分布 个传感器,各传感器的位置分别为,目标位置为 ,这里 表示矩阵转置。不失一般性,假定传感器1位于坐标原点,即 。目标到各传感器的距离为 ,目标到 原点的距离为 , 表示Euclidean范数。因此,目标到传感器 和传感器1的距离差为: (1) 距离差 除以信号传播速度就是目标信号到达两个传感器的时间差,因此,通常将距离差作为观测量建立观测方程,称为距离差方程。 [0007] 公式(l)是一个非线性方程组,采用解非线性方程组的方法(例如:最小二乘法、泰勒级数展开法等)可以解算出单目标情况下的目标位置,实现对单个目标的定位。 [0008] 在多目标情况下,采用距离差方程解算目标位置时,首先需要确定不同传感器中获得的哪些观测量来源于同一目标,即量测数据关联。多传感器多目标量测数据关联问题是 维分配问题,当 时是NP难题,计算复杂度随传感器个数的增加呈指数增长。因此,量测数据关联问题使得分布式多传感器多目标无源定位技术的实用性受到极大限制。 发明内容[0009] 本发明的目的在于提供一种基于投影策略的分布式多传感器多目标无源定位方法。本发明的无源定位方法不需要量测数据关联,本发明的无源定位方法首先将探测区域等间隔量化为投影空间,将传感器获得的观测量通过距离差方程投影到投影空间中,在投影空间中积累不同传感器获得的目标回波的幅度值。通过对投影空间积累值的门限检测判决目标的有无及目标的位置,实现对多个目标的检测和定位。本发明的无源定位方法解决了分布式多传感器多目标无源定位中面临的量测数据关联问题,具有实时性强、定位精度高的优点。 [0010] 为了解决上述现有技术中的不足,本发明是通过下述技术方案实现的:一种基于投影策略的分布式多传感器多目标无源定位方法,其特征在于:其步骤如下: A.准备传感器:在特定探测区域范围内分散布置 个传感器,传感器个数满足 ; B.设定坐标系:设定一个 笛卡尔直角坐标系,坐标系原点位于第1个传感器处,传感器在 笛卡尔直角坐标系的位置坐标分别为 , ,其中, ; C.获取数据:以第1个传感器为主站,其它传感器为辅站,从传感器 , 中读取数据,获得观测量集合;所述观测量集合是传感器接收到的距离差 以及所述距离差对应的幅度值 的集合; D.划分投影空间,令投影值为零:以距离分辨率 为量化间隔划分探测区域,得到三维投影空间;所述三维投影空间由若干个量化后的投影单元和代表每个投影单元的投影值组成,可以表示为: ; 初始化投影空间投影单元值为零: ; E.坐标转换:将步骤D得到的三维投影空间中每一个投影单元坐标转换为在 笛卡尔直角坐标系下的坐标 ; F.投影准备:将E步骤中转换得到的 笛卡尔直角坐标系下的坐标 代 入距离差方程中,计算每个投影单元转换坐标 到传感器 , 与传感器1 之间的距离差 ; G.投影:计算F步骤中得到的距离差 与传感器观测量集合中距离差 的 距离 ,并将该距离 与距离分辨率 的一半作数值对比;若该距离 小于或等于距离分辨率的一半 ,则对应的投影单元投影值加上该距离差幅度值;若该距离大于距离分辨率的一半,则对应的投影单元值不变; H:提取目标:将G步骤中得到的投影值与门限值 作对比,若投影值小于门限值则判没有目标;若投影值大于或等于门限值则判有目标,投影值大于门限值的总个数即为目标的总个数,投影值大于门限值的投影单元的 笛卡尔直角坐标值,即为目标所在的位置。 [0011] 所 述 的 距 离 差 是 指 目 标 到 传 感 器 、 的 距 离 之 差,其 中,;距离差除以信号传播速度就是目标反射的电磁波到达两个传感器的时间差;所述距离差计算的方程为: ; 其中, 表示目标位置, 表示传感器 的位置, 表示矩阵转置, 表示Euclidean范数, 表示目标到传感器 与传感器 的距离差;传感器1位于坐标原点,即 ,因此,距离差方程重写为: 其中, 表示目标到传感器 与传感器1的距离差。 [0012] 所述观测量集合具体是指传感器接收到的距离差及该距离差对应的幅度值,第个传感器的观测量集合可以表示为:其中,为传感器序号,为距离差序号, 为传感器 接收到的距离差的总个数, 为传感器 接收到的第 个距离差, 为距离差 对应的幅度值。 [0013] 所述距离分辨率 是指第三方辐射源能区分的最小的距离;所述第三方辐射源具体是指民用的电视、广播、通信、卫星以及手机基站等民用辐射源等。 [0014] 所述投影空间具体是指用于传感器观测量投影的三维空间,它是将 笛卡尔直角坐标等间隔量化后形成的,可以表示为:其中, 和 分别为 轴、 轴和 轴量化后的值; 、 和 分别是 轴、 轴和 轴量化单元的总个数,分别等于 、 和 , 分别表示传感器在 方向、 方向和 方向的最远探测距离, 表示向上取 整, 表示量化间隔; 表示投影空间中一个投影单元的坐标, 表示投影单元 的投影值。 [0015] 所述坐标转换具体是指将步骤D得到的三维投影空间中每一个投影单元坐标转换为在 笛卡尔直角坐标系下的坐标 ;所述投影具体是指将投影单元的 笛卡尔直角坐标代入距离差方程中,计算每一个投影单元的直角坐标到传感器 , 与传感器1之间的距离差 ,将该距离差与观测量集合 中的距离差进行比较,找出满足条件 的距离差 ,其中 为量化间隔。如果存在满足条件的距离差 ,说明对传感器 而言,该投影单元中存在目标,找到距离差 的幅度值 ,令该投影单元处的投影值等于原来的投影值加上幅度值 。 [0016] 所述投影值具体是指投影单元处通过投影积累的传感器观测量集合中距离差的幅度值。 [0017] 与现有技术相比,本发明所带来的有益的技术效果表现在:1、与现有技术相比,本发明的创新点在于:针对分布式多传感器多目标无源定位技术中复杂的量测数据关联问题,本发明提供了一种无需量测数据关联即可实现分布式多传感器多目标无源定位的新方法。本方法将传感器获得的观测量通过距离差方程投影到投影空间,在投影空间中积累不同传感器获得的目标回波的幅度值。通过对投影空间积累值采用门限检测判断目标的有无及目标的位置,实现对多个目标的检测和定位。本方法解决了量测数据关联的难题,计算复杂度低、实时性强。 [0018] 2、本发明的优点在于基于投影策略,将多传感器观测量投影到投影空间,在投影空间中检测定位目标,避免了分布式多传感器多目标无源定位技术中的量测数据关联问题,具有计算复杂度低、实时性强、定位精度高的优点。 [0019] 3、本发明先设定一个 笛卡尔直角坐标系,然后将该 笛卡尔直角坐标系等间隔量化后得到三维投影空间,对 笛卡尔直角坐标系等间隔量化后使得B步骤中 笛卡尔直角坐标系中连续值,转换成一系列离散的整数值,方便对目标的提取和定位;而本发明中的量化间隔取值越大,投影单元越大,定位精度越低;量化间隔取值越小,投影单元越小,定位精度越高。因此,量化间隔取值不应太大,但量化间隔取值不应过小,量化间隔取值过小,会导致投影单元数过多,计算量过大,难以实现实时处理,因此,量化间隔的取值通常设定为第三方辐射源所能区分的最小距离,即第三方辐射源的距离分辨率 。附图说明 [0020] 图1 为定位系统目标和传感器位置关系示意图。 [0021] 图2 为本发明的流程图。 [0022] 图3 为本发明实施例1的目标定位结果示意图。其中,*表示目标的真实位置,O表示目标的定位结果。 [0023] 图4 为本发明实施例2的目标定位结果示意图。其中,*表示目标的真实位置,O表示目标的定位结果。 [0024] 图5 为本发明实施例3的目标定位结果示意图。其中,*表示目标的真实位置,O表示目标的定位结果。 [0025] 图6 为本发明实施例4的目标定位结果示意图。其中,*表示目标的真实位置,O表示目标的定位结果。 具体实施方式[0026] 实施例1作为本发明一较佳实施例,本实施例主要采用仿真实验的方法进行验证,所有步骤、结论都在MATLABR2012a上验证正确。具体实施例步骤如下: A.准备传感器:在特定探测区域范围内分散布置 个传感器,传感器个数满足 ; B.设定坐标系:设定一个 笛卡尔直角坐标系,坐标系原点位于第1个传感器处,传感器在 笛卡尔直角坐标系的位置坐标分别为 , ,其中, ;传感器2、3、4分别位于 km, km和 km 处,传感器5~11随机散布于由传感器2、3和4决定的三角形区域内;传感器的探测区域为[XX,YY,ZZ]=[50,50,10]km; C.获取数据:以第1个传感器为主站,其它传感器为辅站,从传感器 , 中读取数据,获得观测量集合;所述观测量集合是传感器接收到的距离差 以及所述距离差对应的幅度值 的集合; D.划分投影空间,令投影值为零:以距离分辨率 为量化间隔划分探测区域,得到三维投影空间;所述三维投影空间由若干个量化后的投影单元和代表每个投影单元的投影值组成,可以表示为: ; 初始化投影空间投影单元值为零: ; E.坐标转换:将步骤D得到的三维投影空间中每一个投影单元坐标 转换为在笛卡尔直角坐标系下的坐标 ,其中 ; F.投影准备:将E步骤中转换得到的 笛卡尔直角坐标系下的坐标 代 入距离差方程中,计算每个投影单元转换坐标 到传感器 , 与传感器1 之间的距离差 : , 表示Euclidean范数; G.投影:计算F步骤中得到的距离差 与传感器观测量集合中距离差 的距离 ,并将该距离 与距离分辨率的一半 作数值对比,其中 ; 若 ,则在投影单元 处进行幅度积累 , 其中 为距离差 的幅度值; 若 ,则投影单元 的投影值不变 ; H:提取目标:将G步骤中得到的投影值 与门限值 作对比, 如果 ,判没有目标; 如果 ,判有目标,投影值大于门限的总个数即为目标的总个数,投影值大于门限的投影单元坐标 对应的X-Y-Z笛卡尔直角坐标系中的坐标 ,即为目标所在的位置。 [0027] 实施例2作为本发明一较佳实施例,本实施例主要采用仿真实验的方法进行验证,所有步骤、结论都在MATLABR2012a上验证正确。具体实施例步骤如下: A.准备传感器:在特定探测区域范围内分散布置 个传感器,传感器个数满足 ; B.设定坐标系:设定一个 笛卡尔直角坐标系,坐标系原点位于第1个传感器处,传感器在 笛卡尔直角坐标系的位置坐标分别为 , ,其中, ;传感器2、3、4分别位于 km, km和 km 处,传感器5~21随机散布于由传感器2、3和4决定的三角形区域内;传感器的探测区域为[XX,YY,ZZ]=[50,50,10]km; C.获取数据:以第1个传感器为主站,其它传感器为辅站,从传感器 , 中读取数据,获得观测量集合;所述观测量集合是传感器接收到的距离差 以及所述距离差对应的幅度值 的集合; D.划分投影空间,令投影值为零:以距离分辨率 为量化间隔划分探测区域,得到三维投影空间;所述三维投影空间由若干个量化后的投影单元和代表每个投影单元的投影值组成,可以表示为: ; 初始化投影空间投影单元值为零: ; E.坐标转换:将步骤D得到的三维投影空间中每一个投影单元坐标 转换为在笛卡尔直角坐标系下的坐标 ,其中 ; F.投影准备:将E步骤中转换得到的 笛卡尔直角坐标系下的坐标 代 入距离差方程中,计算每个投影单元转换坐标 到传感器 , 与传感器1 之间的距离差 : , 表示Euclidean范数; G.投影:计算F步骤中得到的距离差 与传感器观测量集合中距离差 的距离 ,并将该距离 与距离分辨率的一半 作数值对比,其中 ; 若 ,则在投影单元 处进行幅度积累 , 其中 为距离差 的幅度值; 若 ,则投影单元 的投影值不变 ; H:提取目标:将G步骤中得到的投影值 与门限值 作对比, 如果 ,判没有目标; 如果 ,判有目标,投影值大于门限的总个数即为目标的总个数,投影值大于门限的投影单元坐标 对应的X-Y-Z笛卡尔直角坐标系中的坐标 ,即为目标所在的位置。 [0028] 实施例3作为本发明一较佳实施例,本实施例主要采用仿真实验的方法进行验证,所有步骤、结论都在MATLABR2012a上验证正确。具体实施例步骤如下: A.准备传感器:在特定探测区域范围内分散布置 个传感器,传感器个数满足 ; B.设定坐标系:设定一个 笛卡尔直角坐标系,坐标系原点位于第1个传感器处,传感器在 笛卡尔直角坐标系的位置坐标分别为 , ,其中, ;传感器2、3、4分别位于 km, km和 km 处,传感器5~31随机散布于由传感器2、3和4决定的三角形区域内;传感器的探测区域为[XX,YY,ZZ]=[50,50,10]km; C.获取数据:以第1个传感器为主站,其它传感器为辅站,从传感器 , 中读取数据,获得观测量集合;所述观测量集合是传感器接收到的距离差 以及所述距离差对应的幅度值 的集合; D.划分投影空间,令投影值为零:以距离分辨率 为量化间隔划分探测区域,得到三维投影空间;所述三维投影空间由若干个量化后的投影单元和代表每个投影单元的投影值组成,可以表示为: ; 初始化投影空间投影单元值为零: ; E.坐标转换:将步骤D得到的三维投影空间中每一个投影单元坐标 转换为在笛卡尔直角坐标系下的坐标 ,其中 ; F.投影准备:将E步骤中转换得到的 笛卡尔直角坐标系下的坐标 代 入距离差方程中,计算每个投影单元转换坐标 到传感器 , 与传感器1 之间的距离差 : , 表示Euclidean范数; G.投影:计算F步骤中得到的距离差 与传感器观测量集合中距离差 的距离 ,并将该距离 与距离分辨率的一半 作数值对比,其中 ; 若 ,则在投影单元 处进行幅度积累 , 其中 为距离差 的幅度值; 若 ,则投影单元 的投影值不变 ; H:提取目标:将G步骤中得到的投影值 与门限值 作对比, 如果 ,判没有目标; 如果 ,判有目标,投影值大于门限的总个数即为目标的总个数,投影值大于门限的投影单元坐标 对应的X-Y-Z笛卡尔直角坐标系中的坐标 ,即为目标所在的位置。 [0029] 实施例4作为本发明一较佳实施例,本实施例主要采用仿真实验的方法进行验证,所有步骤、结论都在MATLABR2012a上验证正确。具体实施例步骤如下: A.准备传感器:在特定探测区域范围内分散布置 个传感器,传感器个数满足 ; B.设定坐标系:设定一个 笛卡尔直角坐标系,坐标系原点位于第1个传感器处,传感器在 笛卡尔直角坐标系的位置坐标分别为 , ,其中, ;传感器2、3、4分别位于 km, km和 km 处,传感器5~21随机散布于由传感器2、3和4决定的三角形区域内;传感器的探测区域为[XX,YY,ZZ]=[20,20,10]km; C.获取数据:以第1个传感器为主站,其它传感器为辅站,从传感器 , 中读取数据,获得观测量集合;所述观测量集合是传感器接收到的距离差 以及所述距离差对应的幅度值 的集合; D.划分投影空间,令投影值为零:以距离分辨率 为量化间隔划分探测区域,得到三维投影空间;所述三维投影空间由若干个量化后的投影单元和代表每个投影单元的投影值组成,可以表示为: ; 初始化投影空间投影单元值为零: ; E.坐标转换:将步骤D得到的三维投影空间中每一个投影单元坐标 转换为在笛卡尔直角坐标系下的坐标 ,其中 ; F.投影准备:将E步骤中转换得到的 笛卡尔直角坐标系下的坐标 代 入距离差方程中,计算每个投影单元转换坐标 到传感器 , 与传感器1 之间的距离差 : , 表示Euclidean范数; G.投影:计算F步骤中得到的距离差 与传感器观测量集合中距离差 的距离 ,并将该距离 与距离分辨率的一半 作数值对比,其中 ; 若 ,则在投影单元 处进行幅度积累 , 其中 为距离差 的幅度值; 若 ,则投影单元 的投影值不变 ; H:提取目标:将G步骤中得到的投影值 与门限值 作对比, 如果 ,判没有目标; 如果 ,判有目标,投影值大于门限的总个数即为目标的总个数,投影值大于门限的投影单元坐标 对应的X-Y-Z笛卡尔直角坐标系中的坐标 ,即为目标所在的位置。 |