CN103576137B - 一种基于成像策略的多传感器多目标定位方法 - Google Patents
一种基于成像策略的多传感器多目标定位方法 Download PDFInfo
- Publication number
- CN103576137B CN103576137B CN201310452905.8A CN201310452905A CN103576137B CN 103576137 B CN103576137 B CN 103576137B CN 201310452905 A CN201310452905 A CN 201310452905A CN 103576137 B CN103576137 B CN 103576137B
- Authority
- CN
- China
- Prior art keywords
- echo
- receiver
- target
- max
- designated
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于成像策略的多传感器多目标的定位方法,它是通过将传感器网络建模为一个二维稀疏阵,把基于传感器网络的多目标定位问题转变为成像问题,采用像BP算法投影策略,用来解决繁杂的数据关联问题;通过贪婪算法搜寻三维图像的峰值,并利用剔除算法清除虚假目标定位多个目标。与传统多目标定位方法相比,本发明从成像的角度对目标的多个回波进行相参积累,从而能够更加提高目标信噪比,方便、准确的解算多目标;在成像空间方便的考虑多目标定位,从而提高多目标定位的概率。
Description
技术领域
本发明属于雷达组网***中的多目标定位的技术领域,它特别涉及到了传感器网络下的多目标定位的技术领域。
背景技术
相比于传统的双站雷达,传感器网络能够在一个较大的空间范围获得目标的多个双站延时,这在军事侦察领域具有非常大的优势。利用这些延时,通过解算有关目标位置的非线性距离方程(原理上类似于GPS***),我们能够在一个脉冲周期中定位目标,非常有利于跟踪任务。但是由于多目标,噪声以及杂波的存在,传感器网络必须在定位前分配包括噪声和杂波在内的不同目标的混合的延时,即完成目标和传感器之间的数据关联问题,从而创建正确的距离方程,这在实际处理中是一个极其复杂的任务。围绕上述多目标的数据关联的问题,一些大量精密的工作如联合概率数据关联(详见“XU Can,Li Zhi.Study on Joint Probability Denstiy Algorithm InMulti-Sensor Data Fusion,Computer Science and Automation Engineering(CSAE),2012IEEE International Conference onvol.3,May.2012,pp.:32-37”)等,不断的被探索。然而当目标的监视环境变得复杂,目标个数增多时,通过传统的数据关联算法来解决问题在计算量方面是非常棘手的。
发明内容:
为了克服传感器网络中多目标的数据关联问题,本发明提出了一种基于成像策略的多传感器多目标的定位方法,其特点是通过将传感器网络建模为一个二维稀疏阵,从而把基于传感器网络的多目标定位问题转变为成像问题,其中基于成像角度的多目标定位的示意图如图1所示;采用像BP算法(即后向投影算法)投影策略,用来解决繁杂的数据关联问题,其中BP算法的示意图如图2所示;通过贪婪算法搜寻三维图像的峰值,并利用剔除算法清除虚假目标定位多个目标;相对传统多目标定位方法,本发明能够在成像空间实现多目标定位,对目标的回波进行相参积累以提高信噪比,从而提高多目标定位的概率。
为了方便描述本发明的内容,首先作以下术语定义:
定义1、双站延时
由发射机发射信号开始计时,该信号经目标散射回接收机时终止计时,即发射机-目标-接收机,这段从发射到回波接收的时间即为该目标对此接收机的双站延时。双站延时可由安装在发射机和接收机上的GPS或者北斗定位***同步得到。
定义2、双站距离
双站延时乘以光速所得到的距离差即为由发射机-目标-接收机的双站距离。
定义3、双站距离方程
假设附图1的原理图中目标位置为x,发射机位置为T,接收机的数目为N,接收机的位置为ri(i=1,2,…,N),则信号从发射机经目标到达第i个接收机的距离为
||T-x||2+||x-ri||2=v*τi
其中,i=1,2,…,N,v表示光速,*表示乘号,τi表示信号从发射机经目标到达第i个接收机的双站延时,||·||2表示向量的2-范数,上式即为关于目标位置的双站距离方程。
定义4、距离单元编号
在雷达***中,将雷达测距的范围划分成若干小的区域并将其编号,雷达根据目标回波信号落入的区域编号计算目标与雷达之间的距离。在此定位***中,根据接收机所测得的各个目标的双站距离,结合距离向分辨率可得到各目标回波的量化编号。
定义5、合成孔径雷达反向投影算法
反向投影算法是基于匹配滤波原理的合成孔径雷达成像算法,其主要通过相干累加实现合成孔径雷达数据的聚焦成像。详细内容可参考文献:“Research on A novel fast back projection algorithmfor strip map bistatic SAR imaging”,Huang Yulin等。
定义6、接收波门
接收的数据从起始时刻到终止时刻的最短时间长度,即发射信号从发射到接收的时间大于这一最短时间长度时才对数据进行接收,否则不接收数据,这一最短时间长度即为接收时间波门。与接收的时间波门对应的则有接收距离波门,即从发射到接收的路程大于某一最小的距离波门时接收数据,否则不接收数据,这一最小的距离波门即为接收距离波门。
定义7、地理空间
地理空间指目标及杂波、噪声等所在的实际的三维场景空间。
定义8、成像空间
雷达成像空间是指将地理空间中的散射点投影到距离向-方位向的二维空间坐标系,该空间由合成孔径雷达成像空间中的两个相互正交的坐标基确定。目前典型合成孔径雷达的成像空间包括距离向-方位向投影空间。
定义9、笛卡尔坐标系
在数学里,笛卡儿坐标系(Cartesian坐标系),也称直角坐标系,是一种正交坐标系。二维的直角坐标系是由两条相互垂直、0点重合的数轴构成的。在平面内,任何一点的坐标是根据数轴上对应的点的坐标设定的。在平面内,任何一点与坐标的对应关系,类似于数轴上点与坐标的对应关系。三维笛卡儿坐标系是在二维笛卡儿坐标系的基础上根据右手定则增加第三维坐标(即Z轴)而形成的。
定义10、数据关联
数据关联指针对多目标多雷达的定位,需要分清各个被动站的测量数据中那些是来源于同一目标的,并把属于同一目标的测量数据组合起来,即完成测量数据关联。详见“王成,李少洪,黄槐.多站被动雷达***的多目标测量数据关联.电子学报.vol.30.No.12.Dec.2002”。
定义11、信噪比
信噪比是指信号的功率与环境噪声功率的比值。详见文献“丁鹭飞,耿富录.雷达原理(第三版).西安电子科技大学出版社.2009.8”。
定义12、单载频脉冲信号
单一载频脉冲信号可以表示为矩形信号与余弦信号的乘积,即其中fc是脉冲的载频,TP是脉冲宽度。详见“皮亦鸣,杨建宇,付毓生,杨晓波.合成孔径雷达成像原理.第一版.电子科技大学出版社.2007.3”。
定义13、相参积累
相参积累指在雷达***中将来自同一目标的回波进行相位的补偿并将回波数据进行累加以提高信噪比,增加信息量的过程。详见“王瑞军.基于非起伏目标的相参积累TBD方法研究.电子科技大学硕士学位论文.2011”。
定义14:动态存储矩阵
动态存储矩阵指能够根据实际的需要进行动态的更新(存储和清除)矩阵中各元素的值的矩阵。
本发明提出了一种基于成像策略的多传感器多目标定位方法,它包括如下步骤:
步骤1、相关参数的初始化
初始化的参数均已知,如下所示:雷达的空间距离分辨率为δr;雷达在距离向上的最大距离单元编号为L,L取值为整数;电磁波传播速度记为v;以发射机位置为原点,在目标,发射机,接收机所在的地理空间建立三维笛卡尔坐标系,则所有位置坐标都是以直角坐标表示;发射机位置记为Pt;空中目标个数记为Ntar,目标的散射截面积记为σ,第i个空中目标的位置记为Ptar(i)(i=1,2,…,Ntar),i为目标序号;地面接收机个数记为Nrec,第j个地面接收机位置记为Prec(j)(j=1,2,…,Nrec),j为接收机序号;每个接收机的接收距离波门均为Rrec0;每个接收机通道的信噪比记为SNR;雷达发射信号为单载频脉冲信号,其脉冲的载频为fc,脉冲宽度为Tp;单载频信号由发射机发射经目标i散射,被第j个接收机接收,经过去载频和离散采样后的回波信号记为Echo(i,j),此回波信号表示为Echo(i,j)=σi,j*exp(-pkRi,j),其中σi,j表示第i个目标对第j个接收机的雷达散射截面积(RCS),Ri,j=||Pt-Ptar(i)||2+||Ptar(i)-Prec(j)||2,表示发射信号经目标i散射到达接收机j的双站距离,exp(·)为自然指数为底的指数函数,||·||2表示向量的2-范数,*表示乘号,π为圆周率,v为光速,fc为载频,p为虚数单位,即p2=-1。
步骤2、创建回波存储矩阵
由安装在发射机和接收机上的全球定位***(GPS或者北斗定位***,收集到空中目标对接收机的双站延时τi,j及相应延时的回波信号Echo(i,j),τi,j表示第i个目标对第j个接收机的双站延时,Echo(i,j)表示第i个目标对第j个接收机相应延时的回波信号。由收集的延时τi,j和对应的回波Echo(i,j)建立一个存储目标回波的动态存储矩阵,回波的动态存储矩阵建立过程如下:
步骤2.1初始化矩阵
根据步骤1中的空中目标i、地面接收机j,定义一个L行Nrec列的动态存储矩阵,记为矩阵的每一单元记为 其中l表示矩阵的行号,c表示矩阵的列号;下标L为距离向上的最大距离单元编号,由所有接收机的最大探测距离而定;上标Nrec表示Nrec个接收机通道,即用于存储接收到的所有目标的回波,转到步骤2.2;
步骤2.2存储目标回波
针对步骤2.1中定义的动态存储矩阵进行动态的回波存储。由各个接收机所获得的所有目标的双站延时τi,j,计算所有目标的双站延时τi,j的距离单元编号,将所有目标的延时τi,j所对应的目标回波Echo(i,j)存储在动态存储矩阵的相应单元中;具体存储步骤见以下步骤2.2.1和步骤2.2.2;
步骤2.2.1根据延时信息计算相应回波的双站距离
根据接收机j所获得的目标i的延时τi,j及对应的目标回波信号Echo(i,j),由双站距离公式Ri,j=v*τi,j计算目标i对第j个接收机的双站距离,转到步骤2.2.2;
步骤2.2.2根据双站距离计算相应回波的距离单元编号
根据步骤2.2.1中所获得的双站距离Ri,j计算得到相应回波的距离单元编号,记为IDi,j,计算公式:其中IDi,j表示目标i的回波在接收机j中的距离单元编号,1≤IDi,j≤L,且IDi,j取值为正整数;其中L为最大的距离单元编号,转到步骤2.2.3;
步骤2.2.3根据距离单元编号将相应的回波进行存储
根据步骤2.2.2中得到的目标i对接收机j的距离单元编号IDi,j,将步骤1提供的去载频和离散采样后的回波信号Echo(i,j),存储于动态存储矩阵中的单元 即 采用传统的遍历方法遍历所有的目标序号i和接收机序号j,得到Ntar个目标对Nrec个接收机的回波信号Echo(i,j)的距离单元编号IDi,j,将去载频和离散采样后的回波信号Echo(i,j)根据编号IDi,j存储于动态存储矩阵的相应单元中,得到完整的Nrec个接收机存储Ntar个目标回波的动态存储矩阵回波的动态存储矩阵示意图如附图3所示;
步骤3、栅格划分
在步骤1中所建立笛卡尔坐标系中,设定一个包含所有目标的三维长方体搜索区域ΩSear,长方体的长、宽、高分别平行于坐标系的x轴、y轴、z轴。长方体的长为10×δr×Nx,长方体的宽为10×δr×Ny,长方体的高为10×δr×Nz,其中Nx,Ny,Nz分别为长方体沿x轴,y轴,z轴的采样点数,长方体搜索区域ΩSear的中心记为P0,根据步骤1中初始的空间距离分辨率δr将长方体搜索区域划分为一组微小的栅格,每一栅格大小为10δr×10δr×10δr,并为每一个栅格单元分配一个代表点,代表点坐标记为Pxyz=P0+[x-Nx/2,y-Ny/2,z-Nz/2]*10*δr,x=1,2…Nx,y=1,2…Ny,z=1,2…Nz,x,y,z分别是长方体沿x轴,y轴,z轴的采样点序号,即Pxyz表示采样点序号为x,y,z的代表点的三维坐标,相邻代表点之间相距10δr。
步骤4、信息投影
采用传统的遍历法在步骤3建立的整个搜索区域ΩSear中,寻找每一栅格代表点所对应各个接收机的回波,采用传统的相参积累方法对每一栅格代表点所对应各个接收机的回波进行相参积累;具体步骤如下:
步骤4.1
初始化定义一个动态存储目标成像信息的三维矩阵,记为ImaNx×Ny×Nz,矩阵的元素单元记为ImaNx×Ny×Nz(k1,k2,k3),k1、k2和k3分别表示矩阵ImaNx×Ny×Nz行号、列号和层号,ImaNx×Ny×Nz(k1,k2,k3)表示矩阵的任一元素,k1=1,2,…Nx,k2=1,2,…Ny,k3=1,2,…Nz,Nx表示矩阵的行数,Ny表示矩阵的列数,Nz表示矩阵的层数,Nx×Ny×Nz表示矩阵的大小。初始化接收机序号j=1,长方体沿x轴,y轴,z轴的采样点序号x=1,y=1,z=1,转到步骤4.2;
步骤4.2
针对步骤3中提供的栅格代表点Pxyz=P0+[x-Nx/2,y-Ny/2,z-Nz/2]*10*δr和接收机j,计算从发射机Pt经栅格代表点Pxyz到达接收机j的双站距离,记为Rxyz,j,则Rxyz,j=||Pt-Pxyz||2+||Pxyz-Prec(j)||2,即Rxyz,j表示发射信号经采样点序号为x,y,z的代表点Pxyz散射回接收机j的双站距离,并由Rxyz,j计算代表点Pxyz对接收机j的距离单元编号,记为IDxyz,j,则即IDxyz,j表示采样点序号为x,y,z的代表点Pxyz对接收机j的距离单元编号,Rrec0和δr分别为步骤1中初始化已知的接收机的接收距离波门和雷达的空间距离分辨率,转到步骤4.3;
步骤4.3
取出步骤2所建立的动态存储矩阵的第j列第IDxyz,j行中的数据作为代表点Pxyz所对应的回波,记为Echo(xyz,j),Echo(xyz,j)表示采样点序号为x,y,z的代表点对接收机j的回波,并将回波Echo(xyz,j)在成像空间相参积累,即ImaNx×Ny×Nz(k1,k2,k3)=ImaNx×Ny×Nz(k1,k2,k3)+Echo(xyz,j)*exp(pkRxyz,j),k1=x,k2=y,k3=z,接收机序号j增加1,并转到步骤4.4;
步骤4.4
如果j≤Nrec,则转到步骤4.2;
如果j>Nrec,则j置1,x增加1,转到步骤4.5;
步骤4.5
如果x≤Nx,则转到步骤4.2;
如果x>Nx,则x置1,y增加1,转到步骤4.6;
步骤4.6
如果y≤Ny,则转到步骤4.2;
如果y>Ny,则y置1,z增加1,转到步骤4.7;
步骤4.7
如果z≤Nz,则转到步骤4.2;
如果z>Nz,则遍历结束,得到成像空间中一个完整的三维图像矩阵ImaNx×Ny×Nz,转到步骤5;
步骤5、提取目标位置
首先,在步骤4得到的三维图像ImaNx×Ny×Nz中找出三维图像ImaNx×Ny×Nz的Ntar个最大值,记录找到的Ntar个最大值的序号为n(n=1,2…,Ntar);然后从三维图像ImaNx×Ny×Nz中剔除Ntar个最大值对应的全部信息;具体步骤如下:
步骤5.1
找出步骤4中得到的三维图像矩阵ImaNx×Ny×Nz的最大值,记为Max_n,转至步骤5.2;
步骤5.2
记录下最大值Max_n对应的三维矩阵中的位置,记为(k1Max_n,k2Max_n,k3Max_n),即ImaNx×Ny×Nz(k1,k2,k3)=Max_n,k1=k1Max_n,k2=k2Max_n,k3=k3Max_n,k1Max_n,k2Max_n,k3Max_n分别为三维图像矩阵中最大值Max_n所对应的行号,列号,层号,则该三维矩阵位置(k1Max_n,k2Max_n,k3Max_n)所对应的地理空间位置记为Pn=P0+[k1Max_n-Nx/2,k2Max_n-Ny/2,k3Max_n-Nz/2]*10*δr,即Pn表示获取的第n个目标的地理空间位置,初始化接收机序号j=1,转到步骤5.3;
步骤5.3
选择第n个位置Pn=P0+[k1Max_n-Nx/2,k2Max_n-Ny/2,k3Max_n-Nz/2]*10*δr和接收机j,由双站距离公式计算从发射机Pt经位置Pn到达接收机j的双站距离,记为Rn,j,则Rn,j=||Pt-Pn||2+||Pn-Prec(j)||2,即Rn,j表示第n个目标对接收机j的双站距离,并由Rn,j计算目标n的回波距离单元编号IDn,j表示目标n对接收机j的回波的距离单元编号,转到步骤5.4;
步骤5.4
在步骤2建立的动态存储矩阵的第j列第IDn,j个距离单元中找出目标n对应的回波,记为Echo(n,j),即 Echo(n,j)表示目标n对接收机j的回波,并将此回波从图像矩阵中剔除,即
ImaNx×Ny×Nz(k1,k2,k3)=ImaNx×Ny×Nz(k1,k2,k3)-Echo(n,j)*exp(pkRn,j)
其中k1=k1Max_n,k2=k2Max_n,k3=k3Max_n,j增加1,并转到步骤5.5;
步骤5.5
如果j≤Nrec,则转到步骤5.3;
如果j≥Nrec,则得到消除第n个最大值信息后的三维矩阵ImaNx×Ny×Nz,n增加1,并转到步骤5.6;
步骤5.6
如果n≤Ntar,则转到步骤5.1;
如果n>Ntar,则转到步骤6;
步骤6
多目标定位结束,得到Ntar个目标在三维地理空间的位置Pn(n=1,2,…,Ntar)。
本发明的创新点在于传统的基于传感器网络的多目标定位中,需要进行复杂的数据关联,而且又需要创建并解算复杂的多目标定位的非线性距离方程组,本发明提供一种基于成像策略的多传感器多目标的定位方法。该方法从成像的角度对目标的多个回波进行相参积累,从而能够更加提高目标信噪比,并方便、准确的解算多目标。
本发明的实质是通过将传感器网络建模为一个二维稀疏阵,从而把基于传感器网络的多目标定位问题看作成像问题,其中基于成像角度的多目标定位的示意图如图1所示;我们采用像BP算法,(详见“Jun,S.;Long,M.;Xiaoling,Z.,Streaming BP for Non-Linear Motion Compensation SAR ImagingBased on GPU,IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,vol.PP,Issue:99,2013,pp.1–16”),即后向投影算法这样的投影策略,能够很容易的被用来解决繁杂的数据关联问题,其中BP算法的示意图如图2所示;另外,投影方法的优势是三维图像的值包含了一些额外的目标信息,譬如目标的雷达散射截面积(RCS)或者存在概率等等,而这些信息在传统的基于数据关联的有关目标位置的双站距离方程组的创建中被舍弃了;通过贪婪算法搜寻三维图像的峰值,并利用剔除算法清除虚假目标,我们能够定位多个目标。
本发明的优点:本发明从成像的角度考虑多目标定位问题,并利用了目标多个方向的回波信号的相位信息,进行相参积累以提高信噪比,相对传统多目标定位方法,本发明能够在成像空间方便的考虑多目标定位,并能够利用成像的知识对目标的回波进行相参积累以提高信噪比,从而提高多目标定位的概率。
附图说明
图1为基于成像角度的多目标定位的模型
其中,多个传感器可看作接收目标散射信号的稀疏阵列。
图2为多传感器定位中的BP算法示意图
其中,A和B表示不同的目标的回波经过距离压缩后脉冲,F表示虚假目标的回波脉冲。N表示传感器接收机的个数,t表示时间轴。针对同一目标在不同传感器中的回波在成像空间进行相参累加,从而进行能量的积累。
图3为多目标的回波存储矩阵示意图
其中,L表示最大的距离单元编号,ID表示距离单元序号,回波存储矩阵为L行Nrec列,ID表示距离单元的序号,n表示接收机通道的序号。每一列对应一个接收机通道,用于存储各个目标的回波,。同一目标对不同传感器的散射回波放在不同的列,不同的目标对同一传感器的散射信号放在同一列的不同距离单元中。如目标1对接收机机1的回波放在回波矩阵的第一列对应的距离单元编号中。
图4为本发明方法的流程图
具体实施方式
本发明主要采用计算机仿真的方法进行验证,所有步骤、结论都在MATLAB-R2011a上验证正确。具体实施步骤如下:
步骤1、基于成像策略的传感器网络多目标定位方法相关参数的初始化
初始化目标定位的参数如下:雷达的空间距离分辨率为δr=10m;雷达在距离向上的距离单元个数为L=10000;电磁波传播速度v=3×108m/s;以发射机位置为原点,在目标,发射机,接收机所在的地理空间建立笛卡尔坐标系;发射机位置记为Pt=[0,0,0];空中目标个数Ntar=10,目标的散射截面积σ=(30×rand(1)+1)m2,目标随机分布在以[100,200,10]×1000为中心,边长两公里的立方体区域中,得第i个目标位置记为Ptar(i)=[100,200,10]×1000+[round((rand-0.5)/0.1),round((rand-0.5)/0.1),round((rand-0.5)/0.1)],i=1,2,…,10;地面接收机个数Nrec=40,接收机随机分布在以[100,50,0]×1000为中心的,边长20km的正方形区域中,得到第个接收机位置记为Prec(j)=[100,50,0]×1000+[rand(1)-0.5,rand(1)-0.5,0]×20×1000,j=1,2,…,40;每个接收机的接收距离波门为Rrec0=300km;每个接收机通道的信噪比为SNR=-20db;雷达发射信号为单载频信号,脉冲的载频fc=2GHz;脉冲宽度Tp=20μs;单载频信号由发射机发射经目标i散射,被第j个接收机接收,经过去载频和离散采样后的回波信号记为Echo(i,j), Ri,j=||Pt-Ptar(i)||2+||Ptar(i)-Prec(j)||2,表示发射信号经目标i散射到达接收机j的双站距离,i=1,2,…,10,j=1,2,…,40,exp(·)为自然指数为底的指数函数,||·||2表示向量的2-范数,π为圆周率,v为光速,p为虚数单位,即p2=-1。
步骤2、创建回波存储矩阵
由安装在发射机和接收机上的全球定位***(GPS)或者北斗定位***,收集到目标对接收机的双站延时τi,j及相应延时的回波信号Echo(i,j),i=1,2,…,10,j=1,2,…,40,i为目标序号,j为接收机序号,τi,j表示第i个目标对第j个接收机的双站延时,Echo(i,j)表示第i个目标对第j个接收机散射的回波信号。由收集到的延时τi,j和对应的回波Echo(i,j)建立一个存储目标回波的动态存储矩阵。回波的动态存储矩阵建立过程如下:
步骤2.1初始化矩阵
根据步骤1中的空中目标i(i=1,2,…10),地面接收机j(j=1,2,…40)以及距离向上的最大距离单元编号L=10000,定义一个10000行40列的动态存储矩阵,记为矩阵的每一单元记为 其中l表示矩阵的行号,c表示矩阵的列号,转到步骤2.2;
步骤2.2存储目标回波
针对步骤2.1中定义的矩阵进行动态的回波存储。数据处理中心根据各个接收机所获得的所有目标的双站延时τi,j,计算此延时的距离单元编号,并据此编号在相应接收机通道中将延时τi,j所对应的目标回波Echo(i,j)存储在矩阵的相应单元中,i=1,2,…10;j=1,2,…40,其中i为目标序号,j为接收机序号,具体步骤见以下步骤2.2.1和步骤2.2.2;
步骤2.2.1根据延时信息计算相应回波的双站距离
数据处理中心根据接收机j所获得的目标i的延时τi,j及对应的回波信号Echo(i,j),计算得到目标i对第j个接收机的双站距离Ri,j=3×108×τi,j,Ri,j则表示第i个目标对第j个接收机的双站距离,i=1,2,…10;j=1,2,…40,其中i为目标序号,j为接收机序号,10为目标数目,40为接收机数目,转到步骤2.2.2;
步骤2.2.2根据双站距离计算相应回波的距离单元编号
根据步骤2.2.1中所获得的双站距离Ri,j计算得到相应回波的距离单元编号,记为IDi,j,则IDi,j表示目标i的回波在接收机j中的距离单元编号,1≤IDi,j≤10000,且IDi,j取值为正整数;i=1,2,…,10,j=1,2,…,40,其中i为目标序号,j为接收机序号,10为目标数目40为接收机数目,转到步骤2.2.3;
步骤2.2.3根据距离单元编号将相应的回波进行存储
根据步骤2.2.2中得到的目标i对接收机j的距离单元编号IDi,j,将数据处理中心所收集得到的相应回波Echo(i,j),存储于矩阵中的单元即遍历所有的目标序号i(i=1,2,…,10)和接收机序号j(j=1,2,…,40),得到10个目标对40个接收机的回波Echo(i,j)的距离单元编号IDi,j(i=1,2,…,10,j=1,2,…,40),将回波Echo(i,j)根据编号IDi,j存储于动态存储矩阵的相应单元中,得到最终完整的10个接收机存储40个目标回波的动态存储矩阵其中1≤IDi,j≤10000且IDi,j取值为正整数,i=1,2,…,10,j=1,2,…,40,其中i为目标序号,j为接收机序号,目标数目为10个,接收机数目为40个,回波存储矩阵示意图如附图3所示;
步骤3、栅格划分
在步骤1中所建立笛卡尔坐标系中设定一个包含所有目标的三维长方体搜索区域ΩSear,长方体的长,宽,高分别平行于坐标系的x轴,y轴,z轴。长方体沿三个坐标轴的采样点数皆为40,长方体的长为10×10×40,长方体的宽为10×10×40,长方体的高为10×10×40,其长方体搜索区域ΩSear的中心记为P0=[100,200,10]×1000,根据步骤1中的空间分辨率δr=10m将搜索区域划分为一组微小的栅格,每一栅格大小为100×100×100,并为每一个栅格单元分配一个代表点,代表点坐标记为Pxyz=[100,200,10]×1000+[x-40/2,y-40/2,z-40/2]×10×10,x=1,2…40,y=1,2…40,z=1,2…40,x,y,z分别是立方体沿x轴,y轴,z轴的采样点序号,即Pxyz表示采样点为x,y,z的代表点的三维坐标,相邻代表点之间相距100m。
步骤4、信息投影
采用传统的遍历法寻找步骤3整个搜索区域ΩSear中的每一栅格代表点针对各个接收机的回波并相参积累;
步骤4.1
初始化定义一个动态存储目标成像信息的三维矩阵,记为ImaN40×40×40,矩阵的元素单元记为Ima40×40×40(k1,k2,k3),k1,k2和k3分别表示矩阵Ima40×40×40行号,列号和层号,Ima40×40×40(k1,k2,k3)表示矩阵的任一元素,k1=1,2,…40,k2=1,2,…40,k3=1,2,…40,40×40×40表示矩阵的大小。初始化接收机序号j=1,长方体沿x轴,y轴,z轴的采样序号为x=1,y=1,z=1,转到步骤4.2;
步骤4.2
选择点Pxyz=[100,200,10]×1000+[x-40/2,y-40/2,z-40/2]×10×10和接收机j,计算从发射机Pt=[0,0,0]经代表点Pxyz到达接收机j的双站距离Rxyz,j=||Pt-Pxyz||2+||Pxyz-Prec(j)||2,即Rxyz,j表示发射信号经采样序号为x,y,z的代表点Pxyz散射回接收机j的双站距离,||·||2表示向量的2-范数,并由此Rxyz,j计算代表点Pxyz对接收机j的距离单元编号,记为IDxyz,j,则即IDxyz,j表示样序号为x,y,z的代表点Pxyz对接收机j的距离单元编号,转到步骤4.3;
步骤4.3
如果0<IDxyz,j≤10000,则取出步骤2建立的回波矩阵的第j列第IDxyz,j行的数据作为代表点Pxyz所对应的回波,记为Echo(xyz,j),Echo(xyz,j)表示采样点序号为x,y,z的代表点对接收机j的回波,并将回波Echo(xyz,j)在成像空间相参积累,即 k1=x,k2=y,k3=z,Rxyz,j表示发射信号经采样序号为x,y,z的代表点Pxyz散射回接收机j的双站距离,*表示乘号,p为虚数单位,即p2=-1,exp(·)为自然指数为底的指数函数,j增加1并转到步骤4.4;
如果IDxyz,j≤0或者IDxyz,j>10000,则接收机序号j增加1,转到步骤4.4;
步骤4.4
如果j≤40,则转到步骤4.2;
如果j>40,则j置1,x增加1,转到步骤4.5;
步骤4.5
如果x≤40,则转到步骤4.2;
如果x>40,则x置1,y增加1,转到步骤4.6;
步骤4.6
如果y≤40,则转到步骤4.2;
如果y>40,则y置1,z增加1,转到步骤4.7;
步骤4.7
如果z≤40,则转到步骤4.2;
如果z>40,则遍历结束,得到成像空间中一个完整的三维图像Ima40×40×40,转到步骤5;
步骤5、提取目标位置
逐步找出步骤4中得到的三维图像的10个最大值,并从图像中剔除最大值对应的全部信息,记录找到的10个最大值的序号为n(i=1,2…,10),并初始化n=1,转到步骤5.1:
步骤5.1
找出步骤4中得到的最终完整的三维图像矩阵Ima40×40×40的最大值,记为Max_n,转至步骤5.2;
步骤5.2
记录下最大值Max_n对应的三维矩阵中的位置,记为(k1Max_n,k2Max_n,k3Max_n),即Ima40×40×40(k1,k2,k3)=Max_n,k1=k1Max_n,k2=k2Max_n,k3=k3Max_n,k1Max_n,k2Max_n,k3Max_n分别为成像矩阵中最大值Max_n所对应的行号,列号,层号,则该三维矩阵位置(k1Max_n,k2Max_n,k3Max_n)所对应的地理空间位置记为Pn=[100,200,10]×1000+[k1Max_n-40/2,k2Max_n-40/2,k3Max_n-40/2]×100,即Pn表示获取的第n个目标的地理空间位置,初始化接收机序号j=1,转到步骤5.3;
步骤5.3
选择第n个位置Pn=[100,200,10]×1000+[k1Max_n-Nx/2,k2Max_n-Ny/2,k3Max_n-Nz/2]×100和接收机j,计算从发射机Pt=[0,0,0]经位置Pn到达接收机j的双站距离,记为Rn,j,Rn,j=||Pt-PMax||2+||PMax-Prec(j)||2,即Rn,j表示第n个目标对接收机j的双站距离,||·||2表示向量的2-范数,并由Rn,j计算目标n的回波距离单元编号IDn,j表示目标n对接收机j的回波的距离单元编号,转到步骤5.4;
步骤5.4
如果0<IDn,j≤10000,则在步骤2建立的回波矩阵的第j列第IDn,j个距离单元中找出目标n对应的回波,记为Echo(n,j),即 Echo(n,j)表示目标n对接收机j的回波,并将此回波从图像矩阵中剔除,即
其中k1=k1Max_n,k2=k2Max_n,k3=k3Max_n,Rn,j表示第n个目标对接收机j的双站距离,*表示乘号,p为虚数单位,即p2=-1,exp(·)为自然指数为底的指数函数,j增加1并转到步骤5.5;
如果IDn,j≤0或者IDn,j>10000,则j增加1,转到步骤5.5;
步骤5.5
如果j≤40,则转到步骤5.3;
如果j>40,则得到消除第n个最大值信息后的三维矩阵Ima40×40×40,n增加1,并转到步骤5.6;
步骤5.6
如果n≤10,则j置1,求取三维矩阵Ima40×40×40的最大值赋给Max_n,转到步骤5.2;
如果n>10,则转到步骤6;
步骤6
多目标定位结束,并得到Ntar个目标在三维地理空间的位置Pn(n=1,2,…,10)。
通过本发明的具体实施可以看出,本发明通过成像的角度对基于传感器网络多目标回波进行了处理,较传统的基于数据关联建立非线性方程的方法快捷而且方便。该方法在基于上述初始化参数,并在接收机通道具有-20db的信噪比的情况下,经过多次蒙特卡洛实验,多目标定位的统计概率可以达到0.8265,可见该方法在低信噪比的情况下依然可以较为准确的对多目标进行定位。
Claims (1)
1.一种基于成像策略的多传感器多目标定位方法,其特征是该方法包括如下步骤:
步骤1、相关参数的初始化
初始化的参数均已知,如下所示:雷达的空间距离分辨率为δr;雷达在距离向上的最大距离单元编号为L,由所有接收机的最大探测距离而定,L取值为整数;电磁波传播速度记为v;以发射机位置为原点,在目标,发射机,接收机所在的地理空间建立三维笛卡尔坐标系,则所有位置坐标都是以直角坐标表示;发射机位置记为Pt;空中目标个数记为Ntar,目标的散射截面积记为σ,第i个空中目标的位置记为Ptar(i),i=1,2,…,Ntar,i为目标序号;地面接收机个数记为Nrec,第j个地面接收机位置记为Prec(j),j=1,2,…,Nrec,j为接收机序号;每个接收机的接收距离波门均为Rrec0;每个接收机通道的信噪比记为SNR;雷达发射信号为单载频脉冲信号,其脉冲的载频为fc,脉冲宽度为Tp;单载频信号由发射机发射经目标i散射,被第j个接收机接收,经过去载频和离散采样后的回波信号记为Echo(i,j),此回波信号表示为Echo(i,j)=σi,j*exp(-pkRi,j),其中σi,j表示第i个目标对第j个接收机的雷达散射截面积(RCS),Ri,j=||Pt-Ptar(i)||2+||Ptar(i)-Prec(j)||2,表示发射信号经目标i散射到达接收机j的双站距离,exp(·)为自然指数为底的指数函数,||·||2表示向量的2-范数,*表示乘号,π为圆周率,v为光速,fc为载频,p为虚数单位,即p2=-1;
步骤2、创建回波存储矩阵
由安装在发射机和接收机上的全球定位***(GPS)或者北斗定位***,收集到空中目标对接收机的双站延时τi,j及相应延时的回波信号Echo(i,j),τi,j表示第i个目标对第j个接收机的双站延时,Echo(i,j)表示第i个目标对第j个接收机相应延时的回波信号;由收集的延时τi,j和对应的回波Echo(i,j)建立一个存储目标回波的动态存储矩阵;回波的动态存储矩阵建立过程如下:
步骤2.1初始化矩阵
根据步骤1中的空中目标i,地面接收机j,以及距离向上的最大距离单元编号L,定义一个L行Nrec列的动态存储矩阵,记为矩阵的每一单元记为 l=1,2,…,L,c=1,2,…,Nrec,其中l表示矩阵的行号,c表示矩阵的列号;下标L为最大的距离单元编号,由所有接收机的最大探测距离而定;上标Nrec表示Nrec个接收机通道,即用于存储接收到的所有目标的回波,转到步骤2.2;
步骤2.2存储目标回波
针对步骤2.1中定义的动态存储矩阵进行动态的回波存储;由各个接收机所获得的所有目标的双站延时τi,j,计算所有目标的双站延时τi,j的距离单元编号,将所有目标的延时τi,j所对应的目标回波Echo(i,j)存储在动态存储矩阵的相应单元中;具体存储步骤见以下步骤2.2.1和步骤2.2.2;
步骤2.2.1根据延时信息计算相应回波的双站距离
根据接收机j所获得的目标i的延时τi,j及对应的目标回波信号Echo(i,j),由双站距离公式Ri,j=v*τi,j计算目标i对第j个接收机的双站距离,转到步骤2.2.2;
步骤2.2.2根据双站距离计算相应回波的距离单元编号
根据步骤2.2.1中所获得的双站距离Ri,j计算得到相应回波的距离单元编号,记为IDi,j,计算公式:其中IDi,j表示目标i的回波在接收机j中的距离单元编号,1≤IDi,j≤L,且IDi,j取值为正整数;其中L为步骤1中已知的最大的距离单元编号,转到步骤2.2.3;
步骤2.2.3根据距离单元编号将相应的回波进行存储
根据步骤2.2.2中得到的目标i对接收机j的距离单元编号IDi,j,将步骤1提供的去载频和离散采样后的回波信号Echo(i,j),存储于动态存储矩阵中的单元 即采用传统的遍历方法遍历所有的目标序号i和接收机序号j,得到Ntar个目标对Nrec个接收机的回波信号Echo(i,j)的距离单元编号IDi,j,将去载频和离散采样后的回波信号Echo(i,j)根据编号IDi,j存储于动态存储矩阵的相应单元中,得到完整的Nrec个接收机存储Ntar个目标回波的动态存储矩阵
步骤3、栅格划分
在步骤1中所建立笛卡尔坐标系中,设定一个包含所有目标的三维长方体搜索区域ΩSear,长方体的长、宽、高分别平行于坐标系的x轴、y轴、z轴;长方体的长为10×δr×Nx,长方体的宽为10×δr×Ny,长方体的高为10×δr×Nz,其中Nx,Ny,Nz分别为长方体沿x轴,y轴,z轴的采样点数,长方体搜索区域ΩSear的中心记为P0,根据步骤1中初始的空间距离分辨率δr将长方体搜索区域划分为一组微小的栅格,每一栅格大小为10δr×10δr×10δr,为每一个栅格单元分配一个代表点,代表点坐标记为Pxyz=P0+[x-Nx/2,y-Ny/2,z-Nz/2]*10*δr,x=1,2…Nx,y=1,2…Ny,z=1,2…Nz,x,y,z分别是长方体沿x轴,y轴,z轴的采样点序号,即Pxyz表示采样点序号为x,y,z的代表点的三维坐标,相邻代表点之间相距10δr;
步骤4、信息投影
采用传统的遍历法在步骤3建立的整个搜索区域ΩSear中,寻找每一栅格代表点所对应各个接收机的回波,采用传统的相参积累方法对每一栅格代表点所对应各个接收机的回波进行相参积累;具体步骤如下:
步骤4.1
初始化定义一个动态存储目标成像信息的三维矩阵,记为ImaNx×Ny×Nz,矩阵的元素单元记为ImaNx×Ny×Nz(k1,k2,k3),k1、k2和k3分别表示矩阵ImaNx×Ny×Nz行号、列号和层号,ImaNx×Ny×Nz(k1,k2,k3)表示矩阵的任一元素,k1=1,2,…Nx,k2=1,2,…Ny,k3=1,2,…Nz,Nx表示矩阵的行数,Ny表示矩阵的列数,Nz表示矩阵的层数,Nx×Ny×Nz表示矩阵的大小;初始化接收机序号j=1,长方体沿x轴,y轴,z轴的采样点序号x=1,y=1,z=1,转到步骤4.2;
步骤4.2
针对步骤3中提供的栅格代表点Pxyz=P0+[x-Nx/2,y-Ny/2,z-Nz/2]*10*δr和接收机j,计算从发射机Pt经栅格代表点Pxyz到达接收机j的双站距离,记为Rxyz,j,则Rxyz,j=||Pt-Pxyz||2+||Pxyz-Prec(j)||2,即Rxyz,j表示发射信号经采样点序号为x,y,z的代 表点Pxyz散射回接收机j的双站距离,并由Rxyz,j计算代表点Pxyz对接收机j的距离单元编号,记为IDxyz,j,则即IDxyz,j表示采样点序号为x,y,z的代表点pxyz对接收机j的距离单元编号,Rrec0和δr分别为步骤1中初始化已知的接收机的接收距离波门和雷达的空间距离分辨率,转到步骤4.3;
步骤4.3
取出步骤2所建立的动态存储矩阵的第j列第IDxyz,j行中的数据作为代表点Pxyz所对应的回波,记为Echo(xyz,j),Echo(xyz,j)表示采样点序号为x,y,z的代表点对接收机j的回波,并将回波Echo(xyz,j)在成像空间相参积累,即
k1=x,k2=y,k3=z;接收机序号j增加1,并转到步骤4.4;
步骤4.4
如果j≤Nrec,则转到步骤4.2;
如果j>Nrec,则j置1,x增加1,转到步骤4.5;
步骤4.5
如果x≤Nx,则转到步骤4.2;
如果x>Nx,则x置1,y增加1,转到步骤4.6;
步骤4.6
如果y≤Ny,则转到步骤4.2;
如果y>Ny,则y置1,z增加1,转到步骤4.7;
步骤4.7
如果z≤Nz,则转到步骤4.2;
如果z>Nz,则遍历结束,得到成像空间中一个完整的三维图像矩阵ImaNx×Ny×Nz,转到步骤5;
步骤5、提取目标位置
首先,在步骤4得到的三维图像ImaNx×Ny×Nz中找出三维图像ImaNx×Ny×Nz的Ntar个最大值,记录找到的Ntar个最大值的序号为n,n=1,2…,Ntar;然后从三维图像ImaNx×Ny×Nz中剔除Ntar个最大值对应的全部信息;具体步骤如下:
步骤5.1
找出步骤4中得到的三维图像矩阵ImaNx×Ny×Nz的最大值,记为Max_n,转至步骤5.2;
步骤5.2
记录下最大值Max_n对应的三维矩阵中的位置,记为(k1Max_n,k2Max_n,k3Max_n),即ImaNx×Ny×Nz(k1,k2,k3)=Max_n,k1=k1Max_n,k2Max_n,k3=k3Max_n,k1Max_n,k2Max_n,k3Max_n分别为三维图像矩阵中最大值Max_n所对应的行号,列号,层号,则该三维矩阵位置(k1Max_n,k2Max_n,k3Max_n)所对应的地理空间位置记为Pn=P0+[k1Max_n-Nx/2,k2Max_n-Ny/2,k3Max_n-Nz/2]*10*δr,即Pn表示获取的第n个目标的地理空间位置,初始化接收机序号j=1,转到步骤5.3;
步骤5.3
选择第n个位置Pn=P0+[k1Max_n-Nx/2,k2Max_n-Ny/2,k3Max_n-Nz/2]*10*δr和接收机j,由双站距离公式计算从发射机Pt经位置Pn到达接收机j的双站距离,记为Rn,j,则Rn,j=||Pt-Pn||2+||Pn-Prec(j)||2,即Rn,j表示第n个目标对接收机j的双站距离,并由Rn,j计算目标n的回波距离单元编号IDn,j表示目标n对接收机j的回波的距离单元编号,转到步骤5.4;
步骤5.4
在步骤2建立的动态存储矩阵的第j列第IDn,j个距离单元中找出目标n对应的回波,记为Echo(n,j),即Echo(n,j)表示目标n对接收机j的回波,并将此回波从图像矩阵中剔除,即
ImaNx×Ny×Nz(k1,k2,k3)=ImaNx×Ny×Nz(k1,k2,k3)-Echo(n,j)*exp(pkRn,j)
其中k1=k1Max_n,k2=k2Max_n,k3=k3Max_n,j增加1,并转到步骤5.5;
步骤5.5
如果j≤Nrec,则转到步骤5.3;
如果j>Nrec,则得到消除第n个最大值信息后的三维矩阵ImaNx×Ny×Nz,n增加1,并转到步骤5.6;
步骤5.6
如果n≤Ntar,则转到步骤5.1;
如果n>Ntar,则转到步骤6;
步骤6
多目标定位结束,得到Ntar个目标在三维地理空间的位置Pn,n=1,2,…,Ntar。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310452905.8A CN103576137B (zh) | 2013-09-27 | 2013-09-27 | 一种基于成像策略的多传感器多目标定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310452905.8A CN103576137B (zh) | 2013-09-27 | 2013-09-27 | 一种基于成像策略的多传感器多目标定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103576137A CN103576137A (zh) | 2014-02-12 |
CN103576137B true CN103576137B (zh) | 2015-05-27 |
Family
ID=50048317
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310452905.8A Active CN103576137B (zh) | 2013-09-27 | 2013-09-27 | 一种基于成像策略的多传感器多目标定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103576137B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104198987B (zh) * | 2014-09-09 | 2017-02-15 | 乐山师范学院 | 一种基于投影策略的分布式多传感器多目标无源定位方法 |
CN105425231B (zh) * | 2015-11-06 | 2017-10-17 | 电子科技大学 | 一种基于分层投影和泰勒展开的多传感器多目标定位方法 |
CN106371081B (zh) * | 2016-08-26 | 2018-11-20 | 电子科技大学 | 一种基于空间栅格数据对齐的多通道量测信息配置方法 |
US10139833B1 (en) * | 2017-05-17 | 2018-11-27 | GM Global Technology Operations LLC | Six-dimensional point cloud system for a vehicle |
CN110313016B (zh) * | 2017-06-15 | 2023-08-15 | 中山大学 | 一种基于稀疏正源分离模型的图像去模糊算法 |
CN107478981B (zh) * | 2017-08-22 | 2019-12-10 | 哈尔滨工业大学 | 一种基于贪婪算法的测试激励集合优选方法 |
CN109558568B (zh) * | 2018-11-21 | 2022-09-16 | 西安电子科技大学 | 一种基于cuda的目标rcs计算方法 |
CN109765555B (zh) * | 2018-12-28 | 2023-06-20 | 通鼎互联信息股份有限公司 | 通过构建虚拟sar图像实现wsn节点定位的方法及装置 |
CN110058194B (zh) * | 2019-04-01 | 2021-10-01 | 四川九洲防控科技有限责任公司 | 目标的方位确定方法及计算机可读存储介质 |
CN111273228B (zh) * | 2019-05-09 | 2023-02-07 | 哈尔滨工程大学 | 一种基于遍历搜索策略的多相干源定位方法 |
CN110426562B (zh) * | 2019-06-24 | 2021-06-04 | 乐山师范学院 | 基于分层搜索和距离空间投影的高精度闪电三维定位方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101049905B1 (ko) * | 2009-02-03 | 2011-07-15 | 국방과학연구소 | 다중분리 레이더 시스템에서 이동 표적의 위치 추적 방법 |
CN101799535A (zh) * | 2009-11-27 | 2010-08-11 | 西安电子科技大学 | Mimo雷达目标方向的估计方法 |
US8482452B2 (en) * | 2010-08-26 | 2013-07-09 | Lawrence Livermore National Security, Llc | Synthetic aperture integration (SAI) algorithm for SAR imaging |
CN102338874B (zh) * | 2011-06-24 | 2013-05-08 | 浙江大学 | 用于被动式多传感器目标跟踪的全局概率数据关联方法 |
CN102621527B (zh) * | 2012-03-20 | 2014-06-11 | 哈尔滨工程大学 | 基于数据重构的宽带相干源的方位估计方法 |
CN102914776B (zh) * | 2012-10-15 | 2014-07-23 | 西安电子科技大学 | 基于模糊c均值算法的多通道SAR动目标定位方法 |
CN103096465B (zh) * | 2013-01-11 | 2015-06-24 | 南京师范大学 | 一种环境自适应的多目标直接定位方法 |
-
2013
- 2013-09-27 CN CN201310452905.8A patent/CN103576137B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN103576137A (zh) | 2014-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103576137B (zh) | 一种基于成像策略的多传感器多目标定位方法 | |
CN103439693B (zh) | 一种线阵sar稀疏重构成像与相位误差校正方法 | |
CN105425231B (zh) | 一种基于分层投影和泰勒展开的多传感器多目标定位方法 | |
CN102749621B (zh) | 一种双基地合成孔径雷达频域成像方法 | |
CN101833095A (zh) | 一种基于空域展开的星机联合sar二维频域成像方法 | |
CN107167784A (zh) | 一种基于多通道比相定位的多人体目标定位跟踪方法 | |
CN104749570B (zh) | 一种移不变机载双基合成孔径雷达目标定位方法 | |
CN103698763A (zh) | 基于硬阈值omp的线阵sar稀疏成像方法 | |
CN102914773B (zh) | 一种多航过圆周sar三维成像方法 | |
CN102788979B (zh) | 一种基于后向投影InSAR成像配准的GPU实现方法 | |
CN102004250B (zh) | 基于频域展开的星机联合双基地合成孔径雷达成像方法 | |
CN103941243A (zh) | 一种基于sar三维成像的自旋式飞行器测高方法 | |
CN103954964B (zh) | 多角度合成孔径雷达数据获取的方法 | |
CN103995250A (zh) | 射频标签轨迹追踪方法 | |
CN103018740B (zh) | 一种基于曲面投影的InSAR成像方法 | |
Inggs et al. | Planning and design phases of a commensal radar system in the FM broadcast band | |
CN103336278A (zh) | 多视角观测下前视三维sar成像方法 | |
CN104833972A (zh) | 一种双基地调频连续波合成孔径雷达频率变标成像方法 | |
CN104698457A (zh) | 一种迭代曲面预测InSAR成像及高度估计方法 | |
CN103207387A (zh) | 一种机载相控阵pd雷达杂波的快速模拟方法 | |
CN102313887A (zh) | 一种星机联合双基地合成孔径雷达成像方法 | |
CN103376447A (zh) | 一种实现非合作双多基地雷达目标三维定位的方法 | |
CN104820221B (zh) | 多基合成孔径雷达的目标三维定位方法 | |
CN104535978A (zh) | 基于互信息的3DInISAR图像配准融合方法 | |
CN104166129A (zh) | 一种实波束雷达迭代最小均方误差角超分辨方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |