CN103235281A - 一种基于相位差增量的相关干涉仪测向方法 - Google Patents
一种基于相位差增量的相关干涉仪测向方法 Download PDFInfo
- Publication number
- CN103235281A CN103235281A CN2013101160501A CN201310116050A CN103235281A CN 103235281 A CN103235281 A CN 103235281A CN 2013101160501 A CN2013101160501 A CN 2013101160501A CN 201310116050 A CN201310116050 A CN 201310116050A CN 103235281 A CN103235281 A CN 103235281A
- Authority
- CN
- China
- Prior art keywords
- phase difference
- increment
- vector
- angle
- azimuth
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 76
- 239000013598 vector Substances 0.000 claims abstract description 82
- 238000005259 measurement Methods 0.000 claims description 24
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 1
- 230000006870 function Effects 0.000 description 10
- 238000004422 calculation algorithm Methods 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 206010028980 Neoplasm Diseases 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 201000011510 cancer Diseases 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 101100042630 Caenorhabditis elegans sin-3 gene Proteins 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 229910052741 iridium Inorganic materials 0.000 description 1
- GKOZUEZYRPOHIO-UHFFFAOYSA-N iridium atom Chemical compound [Ir] GKOZUEZYRPOHIO-UHFFFAOYSA-N 0.000 description 1
- 239000009438 liyan Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- SNICXCGAKADSCV-UHFFFAOYSA-N nicotine Chemical compound CN1CCCC1C1=CC=CN=C1 SNICXCGAKADSCV-UHFFFAOYSA-N 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Images
Landscapes
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
本发明涉及雷达中的相关干涉仪领域,本发明公开了一种基于相位差增量的相关干涉仪测向方法,步骤1.确定方位角增量;步骤2.得到实测相位差增量向量;步骤3.得到无模糊的实测相位差增量向量;步骤4.选出俯仰角为固定角度的所有相位差增量向量样本;步骤5.将实测相位差增量向量与步骤4选出的相位差增量向量样本求相似度;步骤6.得到方位角粗略估计值;步骤7.选取方位角为步骤6中的方位角粗略估计值的所有相位差向量样本;步骤8.将实测相位差向量与步骤7所有选出的相位差向量样本求相似度;步骤9.得到俯仰角、方位角粗略估计值;步骤10.构成新的样本库;步骤11.求相似度,得到俯仰角和方位角的估计值。通过本发明测向系统可在宽频率范围内实现快速、高精度测向。
Description
技术领域
本发明属于雷达通信中的相关干涉仪领域,尤其涉及一种基于相位差增量的相关干涉仪的测向方法。
背景技术
在现有的测向体制中,干涉仪测向具有测向精度高、方法简单、速度快等优点,因此在军事和民用领域都得到了广泛的应用。在军事领域,它可对雷达、通信、测控等辐射源进行定位,能在实施针对性电子干扰和军事打击中发挥重要作用;而在民用领域,它可用于交通管制、生命救援、无线电频谱管理等方面。
相位干涉仪通过测量接收阵元间的相位差来确定来波信号方向,在各种相位干涉仪测向技术中,相关干涉仪通过相关运算可以减小互耦、系统误差的影响,因而被实际测向系统广泛利用(见文献:K Struckman,CorrelationInterferometer Geolocation.IEEE Antennas and Propagation Society International Symposium,Albuquerque,NM,Unitedstates,9-14July2006;S Henault,YMM Antar,S Rajan,R Inkol,S Wang,Impact of Mutual Coupling on Wideband AdcockDirection Finders.Canadian Conference on Electrical and Computer Engineering,Niagara Falls,ON,Canada,4-7May2008)。由于均匀圆阵有着优越的结构特点,它可提供360度全方位测向,各方向都能获得良好的测向性能,使得均匀圆阵在实际阵列信号处理中被广泛采用(见文献:均匀圆阵测向技术研究,季宇;西安电子科技大学硕士学位论文,2009)。在实际许多应用中,如正交跳频网台分选(见文献:跳频信号的侦察技术研究,陈利虎;国防科技大学博士学位论文,2009),短波辐射源定位等,都需要同时获得方位角和俯仰角的信息。因此,必须研究相关干涉仪的二维测向方法。作为一维测向方法的一个自然推广,相关干涉仪的二维测向方法通过在方位角和俯仰角的变化范围内二维搜索相似性函数的极值,获得方位角和俯仰角估计值(见文献:改进的相关干涉仪测向处理方法,李淳,廖桂生,李艳斌;西安电子科技大学学报(自然科学版),2006,33(3):Page(s):400-403;相关处理在干涉仪测向中的应用,刘芬,明望,陶松;电子科学技术评论,2005,3,Page(s):31-37)。
在实际的宽带测向系统中,为测得感兴趣频段内可能的入射信号方向,常采用基于信道化的宽频段测向方法。该方法首先利用快速傅里叶变换(FFT)或多相滤波器组(见文献:Filter bank spectrum sensing for cognitive radios,F.B.Behrouz;IEEE Trans.on signal processing,Volume:56,Issue:5,2008,Page(s):1801-1811)将待测频段划分为许多子频带,然后在各子频带采用相关干涉仪测向方法进行测向(见文献:Direction of arrival estimation via extended phaseinterferometry,Y.W.Wu;S.Rhodes;E.H.Satorius;IEEE transactions on aerospace and electronic systems,Volume:31,Issue:1,1995,Page(s):375-381)。为了提高频率分辨率,测向系统往往要求划分许多子频带,此时如果继续采用同时对方位角和俯仰角进行搜索的二维测向方法,将极大地增加系统的运算量,如何减少相关干涉仪方法的运算量是保证系统实时测向的关键。已有不少学者对此问题进行了研究,提出了一些改进方法,如基于空间夹角的二维测向方法(见文献:相关处理在干涉仪测向中的应用,刘芬,明望,陶松;电子科学技术评论,2005,3,Page(s):31-37)、利用GPU实现的方法(见文献:基于GPU的宽带干涉仪测向算法实现,蒋林鸿,何子述,程婷,贾可新;现代雷达,2012,34(1):Page(s):35-39)、基于维度拆分的二维测向方法(见文献:一种基于维度拆分的相关干涉仪测向方法,程婷,贾可新,何子述,蒋林鸿;CN201110023635.X.电子科技大学.2011-09-07)等。这些方法都能通过减少系统运算量或提高处理速度而保证测向实时性,但是它们都存在一些不足之处。空间夹角测向方法虽然减小了系统运算量,但是牺牲了测向精度;GPU实现通过高效的并行运算可以提高处理速度,但是对硬件平台要求较高;基于维度拆分的测向方法在大幅减小运算量的同时能够保证较高的运算量,可是该方法因在存在相位模糊(见文献:基于相位干涉仪阵列多组解模糊的波达角估计方法研究,龚享铱,袁俊泉,苏令华;电子与信息学报,2006,28(1):Page(s):55-59)时失效而难以应用到宽频率范围测向中。
发明内容
针对现有技术中没有可在宽频率范围内实现快速、高精度测向的方法,提供一种基于相位差增量的相关干涉仪测向方法。
本发明的技术方案如下:一种基于相位差增量的相关干涉仪测向方法,其具体包含以下步骤:
步骤2.计算所选取的k个基线的实测相位差φ′1i,i=1,2,…K,K,得到实测相位差向量旋转转台后,再次得到实测相位差向量计算方位角增量所对应的实测相位差增量Δφ′i=φ′2i-φ′1i,i=1,2,…K,K,得到实测相位差增量向量
步骤4.在相位差增量样本库中选出俯仰角为固定角度、方位角不同时的所有相位差增量向量样本;
步骤6.选出步骤5中相似度最大的相位差增量向量样本所对应的方位角作为方位角粗略估计值,记为
步骤7.在相位差样本库中选取方位角为步骤6中的方位角粗略估计值、俯仰角不同时的所有相位差向量样本;
步骤10.根据步骤9所获得的俯仰角、方位角粗略估计值,选取由俯仰角为方位角为时,(2M+1)·(2N+1)种组合入射方向下的相位差向量样本构成新的样本库,其中2M+1<Q,2N+1<P,其中,Q为按照一定的角度间隔Δθ0将俯仰角测量范围进行划分,划分后所获得的离散俯仰角个数;P为按照一定的角度间隔将方位角测量范围进行划分,划分后所获得的离散方位角个数;
更进一步地,上述步骤5中求相似度的相似度函数为
其中,DΦ(i)为相位差增量样本库中的第i个相位差增量向量样本,P为按照一定角度间隔将方位角测量范围进行划分,划分后所获得的离散方位角的个数;[·]T代表转置运算。
更进一步地,上述步骤8中求相似度的相似度函数
其中,φk(i)代表第i个相位差向量样本Φ(i)的第K个元素,K为相位差向量的维数,Q为按照一定的角度间隔划分俯仰角测量范围所获得的离散俯仰角的个数。
通过采用以上的技术方案,本发明提出的技术方案有如下有益效果:,通过本发明测向系统可在宽频率范围内实现快速、高精度测向。基本不会影响最终的测向精度,相对于标准相关干涉仪运算量具有较大优势,本发明的基于相位差增量的相关干涉仪测向方法基本能达到标准相关干涉仪的测向精度。
附图说明
图1为M元均匀圆阵模型。
图2为某仿真场景下方位角-俯仰角-相关系数关系图。
图3为小半径波长比时,标准相关干涉仪方法、传统维度拆分方法与本发明的改进算法方法的方位角测角精度比较。
图4为小半径波长比时,标准相关干涉仪方法、传统维度拆分方法与本发明的改进算法方法的俯仰角测角精度比较。
图5为大半径波长比时,标准相关干涉仪方法、传统维度拆分方法与本发明的改进算法方法的方位角测角精度比较。
图6为大半径波长比时,标准相关干涉仪方法、传统维度拆分方法与本发明的改进算法方法的俯仰角测角精度比较。
具体实施方式
以下结合附图和实施例,对本发明进行进一步的详细说明。
为后文描述方便,首先进行如下定义:
相关干涉仪:相关干涉仪是一种通过将实测向量与已有的样本库中向量样本进行求相似度运算,选择出样本库中相似度最高的样本所对应的信号到达方向作为入射方向的测向设备。
相位差样本库:相位差样本库是在进行实际入射角度测量前,对各离散的入射频率、入射方位角、入射俯仰角下系统所测得的相位差向量进行有序排列,所获得的一个向量库。
相位差增量(Δφi):相位差增量是一定方位角改变量所对应的同一基线相位差的改变值。
相位差增量样本库:相位差增量样本库是在一定的方位角改变量下,对各离散的入射频率、入射方位角、入射俯仰角下理论计算得到的相位差增量向量进行有序排列,所获得的一个向量库。
相似度函数:相似度函数是用于定量地描述两个向量之间类似程度的映射关系。该映射的输入为待求相似度的两个向量,输出为相似程度。
本发明公开了一种基于相位差增量的相关干涉仪测向方法,其具体包括以下步骤。
步骤1选定一个转台的位置作为基准点,比如以转台的初始位置为基准点,就将转台复位到初始位置,根据阵元半径和入射信号所在频段确定可以确定方位角增量比如设定其中k为不大于1的正常数,d为所选基线的长度,λmin为入射信号所在频段的最小波长。所选基线长度m,n为阵元编号,如7阵元长基线基线长度公式为本领域的常用公式,在此不再赘述。
步骤2.先计算所选取的各基线实测相位差φ′1i,i,i=1,2,…K,K,其中k为所选取基线数目,得到实测相位差向量然后旋转转台度,再次计算所选取的各基线实测相位差φ′2i,i=1,2,…K,K,得到实测相位差向量接着,计算方位角度增量所对应的实测相位差Δφ′i=φ′2i-φ′1i,i=1,2,…K,K,得实测相位差增量向量最后可以把转台复位到初始位置以准备下次测向。
步骤4.在0至90度中任意选取固定的俯仰角值,比如35度,在相位差增量样本库中选出俯仰角为该角度、方位角不同时的所有相位差增量向量样本。
其中,ΔΦ(i)为相位差增量样本库中的第i个相位差增量向量样本,P为按照一定角度间隔将方位角测量范围进行划分,划分后所获得的离散方位角点数,比如为1度或者5度等值,就是将测量范围为360度的方位角分别划分为360个或者72个离散的方位角点数;[·]T代表转置运算。
步骤7.在相位差样本库中选取方位角为步骤6所测出角度、俯仰角不同时的所有相位差向量样本。
其中,φk(i)代表第i个相位差向量样本Φ(i)的第k个元素,Q为按照一定的角度间隔Δθ0划分俯仰角测量范围所获得的离散俯仰角点数。俯仰角的划分与方位角的划分方式类似。步骤2中的K表示相位差向量Φ(i)中含有的元素个数,如7元阵K=7。步骤8中的的k则相当于步骤2中的i,表示向量中的第k个元素,k取1,2,3,…,7,K为相位差向量的维数。
步骤10.根据步骤9所获得的俯仰角、方位角粗略估计值,选取由俯仰角为方位角为时,例如M、N为小于35的正整数,比如M=30,N=10,即存在(2M+1)·(2N+1)种组合入射方向下的相位差向量样本,将(2M+1)·(2N+1)种组合入射方向下的相位差向量样本构成新的样本库,其中2M+1<Q,2N+1<P。
步骤11.将实测相位差向量与步骤10新样本库中的向量样本分别求相似度,其中相似度函数可以为R2(i),i=1,2,…L,(2M+1)·(2N+1),选取最大相似度对应的入射角作为俯仰角和方位角的估计值
上述步骤中,步骤8和步骤11中选择利用相位差向量与其样本库进行相关而不是相位差增量和相位差增量样本库求解最终的方位角与俯仰角,是因为相位差φi比Δφi的等效半径波长比大,实测值比受噪声的影响要大,因而利用与Φ(i)得到的俯仰角、方位角的测向精度高;由于相位差样本库是以转台初始位置为参考坐标建立的,因而相关时用而不是
本发明的工作原理是:
考虑如图1所示的M元均匀圆阵,
圆阵半径为R,以圆心作为参考点。对于入射方向为的远场信号,载波频率为f0,波长为λ=c/f0,其中,c为光速,则第m个阵元相对于参考点的相位为
于是,第n个阵元与第m(≠n)个阵元之间的相位差可表示为
假设相位差向量由K个相位差构成,即Φ=[φ1,φ2,…,φi,…,φK],其中i代表第i种两阵元序号组合方式,记为(mi,ni),φi于是代表第ni个阵元与第mi(≠ni)个阵元之间的相位差;由式(4)知,全方位测向,相位差φm,n无模糊的条件为
以M=7为例进行分析,选取7根长基线(03,14,25,36,40,51,62)用来测向时,式(5)中 即基线相位差无模糊条件为
全方位测向时,相位差增量Δφm,n无模糊条件为
当式(5)相位差无模糊条件不成立时,只要选择满足式(8)的即可保证相位差增量无模糊。值得注意的是,式(8)只能保证以式(6)构建的相位差增量样本库即理论值Δφi在区间[-π,π]内,而相位差增量实测值即实测值可能存在模糊。如理论值φ1i=1.7π,φ2i=0.8π,Δφi=φ2i-φ1i=-0.9π,实测值φ′1i=-0.3π,φ′2i=0.8π,Δφ′i=φ′2i-φ′1i=1.1π,可见Δφi=Δφ′i-2π。易知,无噪声时
由式(9)知,[-π,π]是相位差增量实测值Δφ′i的无模糊区间,即当Δφ′i落在[-π,π]上时,不用解模糊处理,否则在步骤3中就要进行解模糊处理。
实际测向中,由于存在噪声影响,步骤3中解模糊可能会出现边界问题,即当Δφ′i在π或-π附近时,噪声可能使其取值从一侧跳到另一侧,如果再根据式(9)进行模糊处理就会出现错误。假设Δφ′i本在无模糊区间[-π,π]内且接近π或-π,如0.9π,而噪声使所测得的Δφ′i为1.05π,那么步骤3中进行解模糊处理后得到的最终值为-0.95π,显然这是错误的。为了解决此问题,我们可以把式(8)右侧的π减小到kπ(k≤1),如取0.8π,即无模糊区间为[-0.8π,0.8π],然后确定的取值。
比较式(7)与式(4),相位差增量Δφm,n可以视为半径为R'的阵列的基线相位差,由此可见,利用式(6)可避免式(4)高频段测向时的相位模糊问题,因而步骤5中可以利用式(1)进行相关求解初测方位角,此处体现了虚拟阵列变换(见文献:基于虚拟阵列变换的干涉仪测向算法,吴奉微,程婷,贾可新,何子述;2012,34(3):Page(s):42-46)的思想。由于虚拟阵列的等效半径波长比变小,因而利用虚拟阵列基线相位差即相位差增量获得的俯仰角、方位角初略估计值的偏差会增大,但是步骤7~11是利用原始阵列基线相位差获得最终的俯仰角、方位角估计值,因此基本不会影响最终的测向精度。
观察式(6),可见
将式(12)代入式(11),从而可得
即实测相位差增量向量和样本库中同一方位角、不同俯仰角对应的样本的相关系数是相同的,该相关系数的取值仅与方位角有关,因此通过固定俯仰角,对相位差增量样本库进行一维搜索就可以得到初测方位角。图2给出了半径为1米、入射信号频率为70MHz、入射角度为(103°,42°)时,7元均匀圆阵的实测增量相位差向量与相位差增量样本库中样本的相关系数随着方位角俯仰角变化的情况。
获得初测方位角后,就可以将实测相位差向量与相位差样本库中方位角等于估计值的所有样本求相似度,该相似度计算函数必须对俯仰角敏感,如式(2),找使得相似度最高的样本所对应的俯仰角作为俯仰角的粗略估计值。最后以该方位、俯仰角对为中心,选择其附近的几个方位角、俯仰角对,构成一个新的小样本库。将实测相位差与该新样本库中的样本求相似度,通过小范围的二维搜索寻找相似度最高的样本对应的方位角和俯仰角值。当来波信号频率较大时,由于样本库随方位角、俯仰角变化较快,初测方位角、俯仰角小的偏差可能会导致步骤7,10中提取的样本库中的各样本与真实样本之间存在较大的偏差,使最终的测向误差较大。为了减小初略估计值偏差带来的影响,步骤10中应该适当选取稍大点的M,N。
可以看到,若样本库中方位角的测角范围以Δθ0离散成P个值,俯仰角的测角范围以离散成Q个值,那么传统的相关干涉仪二维测向方法中的二维搜索量为P×Q次式(2)运算,传统维度拆分算法主要搜索量缩小为P次式(1)运算和Q+(2M+1)·(2N+1)次式(2)运算,其中式(1)运算时间要小于式(2)运算时间,算法运行时间大幅缩短。改进算法由于在步骤10中要适当增大M,N,其运算量会增加,但相对于标准相关干涉仪运算量还是具有较大优势。
下面具体举例说明本发明的的效果:
考虑一个7元均匀圆阵,阵列直径为1米。信号源为单频信号,入射频率分别为30MHz、750MHz,入射方位角和俯仰角分别为153°和45°。生成样本库时,方位角测量范围为[0°,359°],俯仰角测量范围为[0°,90°],方位角和俯仰角的离散间隔都为1度。的取值由式(8)确定,在信噪比从0dB变化到50dB的条件下,给出了相关干涉仪测向方法、传统维度拆分测向方法与本发明的基于相位差增量的方法获得的方位角和俯仰角测角标准差,其中每个信噪比下进行2000次蒙特卡洛实验。
图3,4所示为小半径波长比时(对应入射频率30MHz),即相位差未出现模糊时,相关干涉仪测向方法、传统维度拆分测向方法与本发明方法所获得的方位角、俯仰角测量标准差的比较。从图3,4中可知三者的测向精度基本完全一致,由于本发明方法求解初测方位角时所用的相位差增量的等效半径波长比大于原波长比,因而在低信噪比时测向精度略好。
图5,6所示为大半径波长比时(对应入射频率为750MHz),即相位差出现模糊时,相关干涉仪测向方法、传统维度拆分测向方法与本发明方法所获得的方位角、俯仰角测量标准差的比较。从图5,6中可见,此时传统维度拆分方法由于相位模糊而导致无法正确测向,本发明的方法基本能达到标准相关干涉仪的测向精度,低信噪比时俯仰角的测向精度比标准相关干涉仪稍差。
为了说明本发明所提出的改进方法的系统运算量相对于标准相关干涉仪能够减少,在CPU为Celeron G530,主频2.40GHz的计算机平台上,用C语言程序实现测向方法,统计10000次仿真下的平均测向时间。当方位角测量范围为[0°,359°],俯仰角测量范围为[0°,90°],方位角和俯仰角的离散间隔都为1度时,两方法的运行时间比较如表1所示。其中,步骤10中取M=30,N=10,此时测向精度基本能到达标准相关干涉仪的测向精度。
表1方法运行时间比较
方法 | 标准相关测向方法 | 本发明的方法 |
测向耗时(ms) | 22.3599 | 0.9664 |
加速比 | 23.14 | / |
上述表格中的测向时间加速比与测向范围和样本离散间隔有关,减小测向范围或增大样本离散间隔,加速比将减小。
这里已经通过具体的实施例子对本发明进行了详细描述,提供上述实施例的描述为了使本领域的技术人员制造或适用本发明,这些实施例的各种修改对于本领域的技术人员来说是容易理解的。本发明并不限于这些例子,或其中的某些方面。本发明的范围通过附加的权利要求进行详细说明。
上述说明示出并描述了本发明的一个优选实施例,但如前所述,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述发明构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。
Claims (7)
1.一种基于相位差增量的相关干涉仪测向方法,其具体包含以下步骤:
步骤4.在相位差增量样本库中选出俯仰角为固定角度、方位角不同时的所有相位差增量向量样本;
步骤7.在相位差样本库中选取方位角为步骤6中的方位角粗略估计值、俯仰角不同时的所有相位差向量样本;
步骤9.选出步骤8中相似度最大的相位差向量样本所对应的俯仰角、方位角作为俯仰角、方位角粗略估计值,记为;
步骤10.根据步骤9所获得的俯仰角、方位角粗略估计值,选取由俯仰角为,方位角为时,种组合入射方向下的相位差向量样本构成新的样本库,其中,,为按照一定的角度间隔将俯仰角测量范围进行划分,划分后所获得的离散俯仰角个数;为按照一定的角度间隔将方位角测量范围进行划分,划分后所获得的离散方位角个数;
6.如权利要求5所述的基于相位差增量的相关干涉仪测向方法,其特征在于所述步骤1中的基准点为转台的初始位置。
7.如权利要求6所述的基于相位差增量的相关干涉仪测向方法,其特征在于所述所述方法还包括:测向完成后将转台复位到初始位置以准备下次测向。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310116050.1A CN103235281B (zh) | 2013-04-03 | 2013-04-03 | 一种基于相位差增量的相关干涉仪测向方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310116050.1A CN103235281B (zh) | 2013-04-03 | 2013-04-03 | 一种基于相位差增量的相关干涉仪测向方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103235281A true CN103235281A (zh) | 2013-08-07 |
CN103235281B CN103235281B (zh) | 2015-01-21 |
Family
ID=48883332
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310116050.1A Expired - Fee Related CN103235281B (zh) | 2013-04-03 | 2013-04-03 | 一种基于相位差增量的相关干涉仪测向方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103235281B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020440A (zh) * | 2014-06-03 | 2014-09-03 | 西安电子科技大学 | 基于l型干涉式线性阵列的二维波达角估计方法 |
CN106199503A (zh) * | 2016-06-29 | 2016-12-07 | 成都中安频谱科技有限公司 | 相关干涉仪测向方法及装置 |
CN107219496A (zh) * | 2017-06-16 | 2017-09-29 | 中国电子科技集团公司第二十八研究所 | 一种改进的相关干涉仪鉴相方法 |
CN108196227A (zh) * | 2017-12-13 | 2018-06-22 | 北京华航无线电测量研究所 | 一种雷达信道化接收信号的前沿测频装置及方法 |
CN108872933A (zh) * | 2018-07-16 | 2018-11-23 | 电子科技大学 | 一种单站盲干涉仪定位方法 |
CN109633526A (zh) * | 2019-01-29 | 2019-04-16 | 杭州电子科技大学 | 基于方向函数的非均匀圆阵相位干涉仪测向解模糊方法 |
CN109959892A (zh) * | 2019-04-11 | 2019-07-02 | 中国电子科技集团公司第三十六研究所 | 一种均匀圆阵双通道干涉仪测向方法、装置和系统 |
CN110007267A (zh) * | 2019-01-29 | 2019-07-12 | 杭州电子科技大学 | 一种基于混合基线的均匀圆阵干涉仪测向解模糊方法 |
CN111427005A (zh) * | 2020-04-22 | 2020-07-17 | 中国人民解放军空军研究院战略预警研究所 | 一种旋转测向装置、系统及情报侦察设备 |
CN113376570A (zh) * | 2021-06-08 | 2021-09-10 | 嘉兴军创电子科技有限公司 | 一种基于信道化处理的相关干涉仪实时测向方法 |
CN113552534A (zh) * | 2021-08-05 | 2021-10-26 | 中国人民解放军战略支援部队航天工程大学 | 基于脉冲信号的旋转基线干涉仪相位标校方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100033377A1 (en) * | 2008-08-11 | 2010-02-11 | The Government of the United States as represented by the U. S. Navy | Sparse Interferometric Array for Radio Frequency Azimuth Direction Finding |
CN102175988A (zh) * | 2011-01-21 | 2011-09-07 | 电子科技大学 | 一种基于维度拆分的相关干涉仪测向方法 |
US8077091B1 (en) * | 2010-07-06 | 2011-12-13 | Intelligent Sciences, Ltd. | System and method for determining a position of a mobile device within a surveillance volume in the presence of multipath interference |
CN102411136A (zh) * | 2011-08-09 | 2012-04-11 | 电子科技大学 | 一种扩展基线解模糊的相位干涉仪测向方法 |
CN102419432A (zh) * | 2011-08-25 | 2012-04-18 | 电子科技大学 | 一种基于虚拟基线的圆阵相位干涉仪二维测向方法 |
CN102419430A (zh) * | 2011-08-17 | 2012-04-18 | 电子科技大学 | 一种基于平行基线的圆阵相位干涉仪二维测向方法 |
KR20120067711A (ko) * | 2010-12-16 | 2012-06-26 | 엘아이지넥스원 주식회사 | 위상차 페어 변경을 이용한 어레이 반경 확장 장치 및 방법 |
CN102540138A (zh) * | 2011-11-25 | 2012-07-04 | 华中科技大学 | 一种多基线相位搜索式二维空间谱测向方法 |
CN102749612A (zh) * | 2012-07-25 | 2012-10-24 | 中国人民解放军信息工程大学 | 一种相关干涉仪的快速测向方法 |
-
2013
- 2013-04-03 CN CN201310116050.1A patent/CN103235281B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100033377A1 (en) * | 2008-08-11 | 2010-02-11 | The Government of the United States as represented by the U. S. Navy | Sparse Interferometric Array for Radio Frequency Azimuth Direction Finding |
US8077091B1 (en) * | 2010-07-06 | 2011-12-13 | Intelligent Sciences, Ltd. | System and method for determining a position of a mobile device within a surveillance volume in the presence of multipath interference |
KR20120067711A (ko) * | 2010-12-16 | 2012-06-26 | 엘아이지넥스원 주식회사 | 위상차 페어 변경을 이용한 어레이 반경 확장 장치 및 방법 |
CN102175988A (zh) * | 2011-01-21 | 2011-09-07 | 电子科技大学 | 一种基于维度拆分的相关干涉仪测向方法 |
CN102411136A (zh) * | 2011-08-09 | 2012-04-11 | 电子科技大学 | 一种扩展基线解模糊的相位干涉仪测向方法 |
CN102419430A (zh) * | 2011-08-17 | 2012-04-18 | 电子科技大学 | 一种基于平行基线的圆阵相位干涉仪二维测向方法 |
CN102419432A (zh) * | 2011-08-25 | 2012-04-18 | 电子科技大学 | 一种基于虚拟基线的圆阵相位干涉仪二维测向方法 |
CN102540138A (zh) * | 2011-11-25 | 2012-07-04 | 华中科技大学 | 一种多基线相位搜索式二维空间谱测向方法 |
CN102749612A (zh) * | 2012-07-25 | 2012-10-24 | 中国人民解放军信息工程大学 | 一种相关干涉仪的快速测向方法 |
Non-Patent Citations (3)
Title |
---|
吴奉微等: "基于虚拟阵列变换的干涉仪测向算法", 《现代雷达》 * |
张春杰等: "非均匀圆阵天线模型解模糊误差研究", 《系统工程与电子技术》 * |
郑攀等: "二维干涉仪测向算法研究", 《现代电子技术》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020440A (zh) * | 2014-06-03 | 2014-09-03 | 西安电子科技大学 | 基于l型干涉式线性阵列的二维波达角估计方法 |
CN106199503A (zh) * | 2016-06-29 | 2016-12-07 | 成都中安频谱科技有限公司 | 相关干涉仪测向方法及装置 |
CN107219496A (zh) * | 2017-06-16 | 2017-09-29 | 中国电子科技集团公司第二十八研究所 | 一种改进的相关干涉仪鉴相方法 |
CN107219496B (zh) * | 2017-06-16 | 2019-06-11 | 中国电子科技集团公司第二十八研究所 | 一种改进的相关干涉仪鉴相方法 |
CN108196227A (zh) * | 2017-12-13 | 2018-06-22 | 北京华航无线电测量研究所 | 一种雷达信道化接收信号的前沿测频装置及方法 |
CN108872933A (zh) * | 2018-07-16 | 2018-11-23 | 电子科技大学 | 一种单站盲干涉仪定位方法 |
CN108872933B (zh) * | 2018-07-16 | 2022-02-08 | 电子科技大学 | 一种单站盲干涉仪定位方法 |
CN110007267A (zh) * | 2019-01-29 | 2019-07-12 | 杭州电子科技大学 | 一种基于混合基线的均匀圆阵干涉仪测向解模糊方法 |
CN110007267B (zh) * | 2019-01-29 | 2020-08-18 | 杭州电子科技大学 | 一种基于混合基线的均匀圆阵干涉仪测向解模糊方法 |
CN109633526B (zh) * | 2019-01-29 | 2020-09-01 | 杭州电子科技大学 | 基于方向函数的非均匀圆阵相位干涉仪测向解模糊方法 |
CN109633526A (zh) * | 2019-01-29 | 2019-04-16 | 杭州电子科技大学 | 基于方向函数的非均匀圆阵相位干涉仪测向解模糊方法 |
CN109959892A (zh) * | 2019-04-11 | 2019-07-02 | 中国电子科技集团公司第三十六研究所 | 一种均匀圆阵双通道干涉仪测向方法、装置和系统 |
CN109959892B (zh) * | 2019-04-11 | 2023-03-10 | 中国电子科技集团公司第三十六研究所 | 一种均匀圆阵双通道干涉仪测向方法、装置和系统 |
CN111427005A (zh) * | 2020-04-22 | 2020-07-17 | 中国人民解放军空军研究院战略预警研究所 | 一种旋转测向装置、系统及情报侦察设备 |
CN113376570A (zh) * | 2021-06-08 | 2021-09-10 | 嘉兴军创电子科技有限公司 | 一种基于信道化处理的相关干涉仪实时测向方法 |
CN113552534A (zh) * | 2021-08-05 | 2021-10-26 | 中国人民解放军战略支援部队航天工程大学 | 基于脉冲信号的旋转基线干涉仪相位标校方法 |
CN113552534B (zh) * | 2021-08-05 | 2022-02-01 | 中国人民解放军战略支援部队航天工程大学 | 基于脉冲信号的旋转基线干涉仪相位标校方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103235281B (zh) | 2015-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103235281B (zh) | 一种基于相位差增量的相关干涉仪测向方法 | |
CN102175988B (zh) | 一种基于维度拆分的相关干涉仪测向方法 | |
CN102411136B (zh) | 一种扩展基线解模糊的相位干涉仪测向方法 | |
CN103323811B (zh) | 基于虚拟同心圆环阵列的参数估计方法 | |
CN102419432A (zh) | 一种基于虚拟基线的圆阵相位干涉仪二维测向方法 | |
CN106501770B (zh) | 基于幅相误差阵列的远近场宽带混合源中近场源定位方法 | |
CN104914408B (zh) | 基于中国余数定理的频率、doa联合测量方法以及装置 | |
CN101963659A (zh) | 利用三维正交干涉仪实现宽频段测向的方法 | |
CN102419430A (zh) | 一种基于平行基线的圆阵相位干涉仪二维测向方法 | |
CN106980104B (zh) | 用于传感器阵列的信号波达方向自校正方法 | |
Chen et al. | Two-dimensional direction of arrival estimation for improved archimedean spiral array with MUSIC algorithm | |
CN106950529A (zh) | 声矢量近场源esprit和music参数估计方法 | |
CN101252382B (zh) | 一种宽频段信号极化与doa估计方法及装置 | |
CN106997037A (zh) | 声矢量传感器阵列空间旋转解相干到达角估计方法 | |
CN103869298A (zh) | 一种分布式mimo天波超视距雷达海杂波仿真方法 | |
CN102445680A (zh) | 短波宽带相关干涉仪投影技术 | |
Sabibolda et al. | Improving the Accuracy And Performance Speed Of The Digital Spectral-Correlation Method For Measuring Delay In Radio Signals And Direction Finding | |
CN108983145A (zh) | 电磁矢量传感器阵列宽带相干源定位方法 | |
CN107271959B (zh) | 基于多天线正交频分复用调制信号的定位方法 | |
CN102200573A (zh) | 一种对近场目标信号来波方向进行测定的方法 | |
CN111190135A (zh) | 一种适用于任意阵列的校准方法 | |
CN112968303B (zh) | 阵列天线、定位方法、定位系统、ble定位装置及ble设备 | |
CN113270729B (zh) | 立体阵列天线、定位方法和系统、ble芯片及无线mcu芯片 | |
Shi et al. | Algorithm of signal processing with five-channel interferometer | |
Mollai et al. | Wideband two dimensional interferometric direction finding algorithm using base-triangles and a proposed minimum planar array |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150121 Termination date: 20150403 |
|
EXPY | Termination of patent right or utility model |