CN114545324A - 一种适用于非均匀阵列的快速测向方法 - Google Patents

一种适用于非均匀阵列的快速测向方法 Download PDF

Info

Publication number
CN114545324A
CN114545324A CN202210433326.8A CN202210433326A CN114545324A CN 114545324 A CN114545324 A CN 114545324A CN 202210433326 A CN202210433326 A CN 202210433326A CN 114545324 A CN114545324 A CN 114545324A
Authority
CN
China
Prior art keywords
array
discrete
sequence
fourier transform
incoming wave
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
Application number
CN202210433326.8A
Other languages
English (en)
Other versions
CN114545324B (zh
Inventor
于明成
李继锋
李晃
朱文明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Yu'an Defense Technology Co ltd
Yangzhou Yuan Electronic Technology Co Ltd
Original Assignee
Nanjing Yu'an Defense Technology Co ltd
Yangzhou Yuan Electronic Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing Yu'an Defense Technology Co ltd, Yangzhou Yuan Electronic Technology Co Ltd filed Critical Nanjing Yu'an Defense Technology Co ltd
Priority to CN202210433326.8A priority Critical patent/CN114545324B/zh
Publication of CN114545324A publication Critical patent/CN114545324A/zh
Application granted granted Critical
Publication of CN114545324B publication Critical patent/CN114545324B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/14—Systems for determining direction or deviation from predetermined direction
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Abstract

本发明公开一种适用于非均匀阵列的快速测向方法，根据阵列参数确定离散空间序列的长度，构造并初始化离散空间序列，根据阵元数目及序列长度，来选择是使用离散傅里变换DFT还是选择快速傅里叶变换FFT计算，读取接收来波波长和接收通道并更新离散空间序列，对离散空间序列做DFT或FFT运算，获得与实部方向图最大值相对应的角度，然后迭代求解以及计算来波方向。本发明通过对阵元输出信号序列的适当处理，使得可以利用FFT快速获得方向图峰值所在的大概位置。基于此大概位置，利用迭代算法，可以获得方向图峰值的精确位置，进而实现精确测向。

Description

一种适用于非均匀阵列的快速测向方法
技术领域
本发明属于波达测向技术，具体涉及一种适用于非均匀阵列的快速测向方法。
背景技术
波达方向(DOA)估计是通过对分布在不同位置上的天线阵列接收的空间来波信号进行处理，获得来波方向的一种技术，是阵列信号处理研究的重要内容之一，在电子对抗、雷达、通信、声呐等领域得到了广泛应用。
根据阵元间距的特点，天线阵列可分为等间距的均匀阵和不等间距的非均匀阵。均匀阵列要求阵元间距不能大于来波波长的一半。非均匀阵列则突破了这个限制，能以较少的阵元获得较大的阵面有效口径，具有成本低、测向分辨力及测向精度高的优点。但同时，非均匀阵也存在测向算法的运算效率低和如果阵列流形设计不当容易引起测向模糊的问题。
非均匀阵测向模型如图1所示，
Figure 899810DEST_PATH_IMAGE001
为来波方向，
Figure 689911DEST_PATH_IMAGE002
为第
Figure 850765DEST_PATH_IMAGE003
个阵元距 最左边的参考阵元E0的以最短半波长为单位的距离；
Figure 587777DEST_PATH_IMAGE004
，
Figure 907900DEST_PATH_IMAGE005
；
式中，
Figure 806586DEST_PATH_IMAGE006
为正整数，
Figure 202539DEST_PATH_IMAGE007
为阵列须适应的最短波长。
当
Figure 743241DEST_PATH_IMAGE008
时，非均匀阵退化为均匀阵。
每个阵元的接收信号经正交下变频后可以表示为
Figure 917871DEST_PATH_IMAGE009
；
Figure 253037DEST_PATH_IMAGE010
；
Figure 122904DEST_PATH_IMAGE011
；
式中，
Figure 264035DEST_PATH_IMAGE012
为接收信号的相位函数，
Figure 496434DEST_PATH_IMAGE013
为第
Figure 923873DEST_PATH_IMAGE003
个阵元输出的接收信号分量，
Figure 874511DEST_PATH_IMAGE014
为噪声分量。对
Figure 553754DEST_PATH_IMAGE015
的一次采样称为一次快拍。
令：
Figure 640659DEST_PATH_IMAGE016
(3)
Figure 989732DEST_PATH_IMAGE017
则有
Figure 958825DEST_PATH_IMAGE018
对于上述非均匀阵测向模型，其无模糊测向方法的前提为：相邻阵元间距
Figure 379442DEST_PATH_IMAGE019
最大公约数为1。
到目前为止，阵列测向算法可以分为三大类：
1)基于波束形成的算法；
2)基于接收信号协方差矩阵的算法，这类算法具体还可以分为基于空间谱估计的算法和基于子空间分解的算法两种；
3)基于相位干涉仪解模糊的测向算法。
从工程实现以及测向性能两方面综合考虑，目前工程上实用的非均匀阵测向方法首选基于波束形成的算法。
所谓基于波束形成的测向算法的主要内容如下：
构造导向信号
Figure 9269DEST_PATH_IMAGE020
，
Figure 857139DEST_PATH_IMAGE021
和如下的波束形成器
Figure 579107DEST_PATH_IMAGE022
Figure 537836DEST_PATH_IMAGE023
(6)
式中，
Figure 537016DEST_PATH_IMAGE024
为
Figure 555788DEST_PATH_IMAGE025
的共轭，
Figure 499473DEST_PATH_IMAGE026
为方向图，
Figure 261893DEST_PATH_IMAGE027
为实部方向图，
Figure 364847DEST_PATH_IMAGE028
是第
Figure 351257DEST_PATH_IMAGE003
个阵元接收 的来波信号，
Figure 985501DEST_PATH_IMAGE029
是构造的与第
Figure 223715DEST_PATH_IMAGE003
个阵元相对应的导向信号。
考虑理想无噪声情况，有
Figure 994225DEST_PATH_IMAGE030
Figure DEST_PATH_IMAGE031
；
当
Figure 151537DEST_PATH_IMAGE032
=
Figure 698843DEST_PATH_IMAGE033
时，
Figure 599803DEST_PATH_IMAGE026
=
Figure 490398DEST_PATH_IMAGE027
达到其最大值
Figure 959557DEST_PATH_IMAGE034
，
Figure 302814DEST_PATH_IMAGE034
也称为信号处理增益。因此，如能获得
Figure 7464DEST_PATH_IMAGE027
的极大值所对应的
Figure 752567DEST_PATH_IMAGE035
，则由式(3)可得来波方向为
Figure 641894DEST_PATH_IMAGE036
，其中
Figure 472447DEST_PATH_IMAGE032
是导 向信号
Figure 715209DEST_PATH_IMAGE029
中的角度变量，
Figure 518080DEST_PATH_IMAGE033
为相关中间计算参数。
式(6)给出的波束形成器具有类傅里叶变换形式，对于非均匀阵列，只能采用类似 DFT的计算方式，在
Figure 391358DEST_PATH_IMAGE037
区间内，搜索使实部方向图达到极大值的
Figure 505945DEST_PATH_IMAGE035
。而搜索
Figure 755661DEST_PATH_IMAGE035
的计算 量取决于两个因素，一个是搜索步长，另一个是阵元数目。由式(6)，计算一个搜索点上的
Figure 898191DEST_PATH_IMAGE027
值需要
Figure 4687DEST_PATH_IMAGE038
次乘法、
Figure 544253DEST_PATH_IMAGE039
次加法共约
Figure 269764DEST_PATH_IMAGE040
次实数运算。设搜索步长为
Figure 640702DEST_PATH_IMAGE041
，则采用类似 DFT的方式，计算全部搜索点上的
Figure 855783DEST_PATH_IMAGE027
，需要的运算量为
Figure 335175DEST_PATH_IMAGE042
；
设最大来波方向为
Figure 661114DEST_PATH_IMAGE043
时，要求的测向精度为
Figure 417717DEST_PATH_IMAGE044
，则由式(3)，搜索步长必须
Figure 475803DEST_PATH_IMAGE045
才能满足测向精度要求。
例如，设
Figure 989961DEST_PATH_IMAGE043
=60°，
Figure 181908DEST_PATH_IMAGE044
=0.1°，有
Figure 730701DEST_PATH_IMAGE046
4.36×10-4，取
Figure 707491DEST_PATH_IMAGE041
=4×10-4，对于一个16阵元 的非均匀阵，共需要
Figure 771262DEST_PATH_IMAGE047
=1.6×105次运算，相当于2048个阵元构成的均匀阵利用FFT计算 时的运算量。
综上所述，现有技术方案对于非均匀阵来说，由于阵元的不等间距布置，使得传统的波束形成算法不能简单地利用FFT实现快速计算，当阵元数目较大时，导致计算量增大，不利于实时实现。
发明内容
发明目的：本发明的目的在于解决现有技术中存在的不足，提供一种适用于非均匀阵列的快速测向方法，通过对阵元输出信号序列的适当处理，使得可以利用FFT快速获得方向图峰值所在的大概位置。基于此大概位置，利用迭代算法，可以获得方向图峰值的精确位置，进而实现精确测向。
技术方案：本发明的一种适用于非均匀阵列的快速测向方法，所述非均匀阵元包 括N个阵元{
Figure 970162DEST_PATH_IMAGE048
}，来波方向为
Figure 373461DEST_PATH_IMAGE001
，
Figure 773350DEST_PATH_IMAGE049
为第
Figure 324417DEST_PATH_IMAGE003
个阵元距最左边的参考阵元E 0 的以最 短半波长为单位的距离，
Figure 61429DEST_PATH_IMAGE007
为该非均匀阵列的最短波长；
Figure 506185DEST_PATH_IMAGE011
； 该非均匀阵列中每个阵元的接收信号经正交下变频后表示为：
Figure 404871DEST_PATH_IMAGE050
；
Figure 646497DEST_PATH_IMAGE010
；
上式中，
Figure 983937DEST_PATH_IMAGE006
为正整数，
Figure 768354DEST_PATH_IMAGE012
为接收信号的相位函数，
Figure 103520DEST_PATH_IMAGE013
为第
Figure 363600DEST_PATH_IMAGE003
个阵元输出的接 收信号分量，
Figure 442415DEST_PATH_IMAGE014
为噪声分量；且
Figure 97649DEST_PATH_IMAGE004
，
Figure 603717DEST_PATH_IMAGE005
；
Figure 351093DEST_PATH_IMAGE051
为当前来波波长；
Figure 968019DEST_PATH_IMAGE052
和
Figure 992607DEST_PATH_IMAGE053
分别 是信号
Figure 669576DEST_PATH_IMAGE054
的实部和虚部；
Figure 904248DEST_PATH_IMAGE055
表示复数；
具体快速测向方法包括以下步骤：
步骤（1）、根据阵列参数确定离散空间序列的长度
Figure 324865DEST_PATH_IMAGE056
；
Figure 531856DEST_PATH_IMAGE057
；
Figure 566676DEST_PATH_IMAGE058
；
Figure 226328DEST_PATH_IMAGE059
为
Figure 247374DEST_PATH_IMAGE060
的最小整数，
Figure 308870DEST_PATH_IMAGE061
为高于副瓣电平的主瓣宽度；
步骤（2）、构造并初始化长度为
Figure 999746DEST_PATH_IMAGE056
的离散空间序列，
Figure 146694DEST_PATH_IMAGE062
；
步骤（3）、根据阵元数目
Figure 971430DEST_PATH_IMAGE034
及序列长度
Figure 887433DEST_PATH_IMAGE056
，来选择是使用离散傅里变换DFT还是选 择快速傅里叶变换FFT计算，若
Figure 780170DEST_PATH_IMAGE063
，选择FFT，否则，选择DFT；
步骤（4）、读取接收来波波长和接收通道输出
Figure 148835DEST_PATH_IMAGE054
，并更新离散空间序列Z m ，若m= p n , Z m =x n* conj(x 0)；
步骤（5）、根据步骤（3）的选择对离散空间序列Z m 做离散傅里变换DFT或快速傅里 叶变换FFT运算，获得实部方向图
Figure 511683DEST_PATH_IMAGE064
，并求得与
Figure 219876DEST_PATH_IMAGE064
最大值相对应的
Figure 846029DEST_PATH_IMAGE035
；
步骤（6）、基于
Figure 701990DEST_PATH_IMAGE035
进行迭代求解，获得使实部方向图
Figure 868529DEST_PATH_IMAGE027
最大的精确
Figure 759125DEST_PATH_IMAGE065
；
步骤（7）、计算来波方向
Figure 477551DEST_PATH_IMAGE066
。
进一步地，基于
Figure 820808DEST_PATH_IMAGE054
重新构造离散空间序列如下
Figure 259879DEST_PATH_IMAGE067
式中，
Figure 942664DEST_PATH_IMAGE024
为
Figure 645041DEST_PATH_IMAGE025
的共轭，则
Figure 537911DEST_PATH_IMAGE068
Figure DEST_PATH_IMAGE069
Figure 406772DEST_PATH_IMAGE070
上式中
Figure 803118DEST_PATH_IMAGE026
是离散空间序列
Figure DEST_PATH_IMAGE071
的傅里叶变换。
进一步地，通过DFT或FFT实现
Figure 348500DEST_PATH_IMAGE027
的快速计算，如下
Figure 666349DEST_PATH_IMAGE072
Figure DEST_PATH_IMAGE073
式中，
Figure 712802DEST_PATH_IMAGE074
为波束编号，
Figure DEST_PATH_IMAGE075
；设
Figure 885027DEST_PATH_IMAGE064
的极大值对应的波束编号为
Figure 929206DEST_PATH_IMAGE076
， 则有
Figure 672034DEST_PATH_IMAGE077
且
Figure 459862DEST_PATH_IMAGE035
位于
Figure 830800DEST_PATH_IMAGE027
主瓣。
进一步地，为提高精度，利用下式给出的迭代计算方法获得更加精确的
Figure 731367DEST_PATH_IMAGE065
。
初始化
Figure 23808DEST_PATH_IMAGE078
Figure 412064DEST_PATH_IMAGE079
；
Figure 840771DEST_PATH_IMAGE080
；
Figure 430015DEST_PATH_IMAGE081
；
Figure 944173DEST_PATH_IMAGE082
。
有益效果：当非均匀阵列的阵元数目较多时，本发明与类DFT的波束形成测向算法相比，极大的改进了计算效率。由于传统相位干涉仪也属于非均匀阵列，因此，本发明还适用于基线数目较多的相位干涉仪快速测向。
与现有技术相比，本发明具有以下优点：
1)、可以利用FFT实现快速计算；
2）、结合迭代求解，确保有高的测向精度；
3)、易于工程实现。
附图说明
图1为本发明的非均匀阵列示意图；
图2为本发明的工作流程示意图；
图3为实施例中测向误差随信噪比变化曲线示例图；
图4为实施例中来波方向为0°时，实部方向图示例图。
具体实施方式
下面对本发明技术方案进行详细说明，但是本发明的保护范围不局限于所述实施例。
如图1所示，本实施例的适用于非均匀阵列的快速测向方法，其非均匀阵元包括N 个阵元{
Figure 136120DEST_PATH_IMAGE048
}，来波方向为
Figure 684913DEST_PATH_IMAGE001
，
Figure 428747DEST_PATH_IMAGE049
为第
Figure 430201DEST_PATH_IMAGE003
个阵元距最左边的参考阵元E 0 的以最短 半波长为单位的距离，
Figure 425839DEST_PATH_IMAGE007
为该非均匀阵列的最短波长；
Figure 829138DEST_PATH_IMAGE011
；该 非均匀阵列中每个阵元的接收信号经正交下变频后表示为：
Figure 494606DEST_PATH_IMAGE083
；
Figure 983356DEST_PATH_IMAGE010
；
上式中，
Figure 720368DEST_PATH_IMAGE006
为正整数，
Figure 40491DEST_PATH_IMAGE012
为接收信号的相位函数，
Figure 204756DEST_PATH_IMAGE013
为第
Figure 869218DEST_PATH_IMAGE003
个阵元输出的接 收信号分量，
Figure 878762DEST_PATH_IMAGE014
为噪声分量；且
Figure 318971DEST_PATH_IMAGE004
，
Figure 388558DEST_PATH_IMAGE005
，
Figure 789584DEST_PATH_IMAGE051
为当前来波波长；
Figure 602819DEST_PATH_IMAGE052
和
Figure 897534DEST_PATH_IMAGE053
分别 是信号
Figure 403602DEST_PATH_IMAGE054
的实部和虚部；
Figure 806770DEST_PATH_IMAGE055
表示复数。
如图2所示，本实施例的快速测向方法包括以下步骤：
步骤（1）、根据阵列参数确定离散空间序列的长度
Figure 423696DEST_PATH_IMAGE056
；
Figure 448284DEST_PATH_IMAGE057
；
Figure 125253DEST_PATH_IMAGE058
；
Figure 359925DEST_PATH_IMAGE059
为
Figure 514963DEST_PATH_IMAGE060
的最小整数，
Figure 987533DEST_PATH_IMAGE061
为高于副瓣电平的主瓣宽度；
步骤（2）、构造并初始化长度为
Figure 261169DEST_PATH_IMAGE056
的离散空间序列
Figure 186400DEST_PATH_IMAGE084
；
基于
Figure 207445DEST_PATH_IMAGE085
重新构造离散空间序列如下
Figure 3363DEST_PATH_IMAGE086
；
式中，
Figure 959818DEST_PATH_IMAGE024
为
Figure 106765DEST_PATH_IMAGE025
的共轭，则
Figure 931502DEST_PATH_IMAGE087
Figure 847505DEST_PATH_IMAGE088
Figure 958549DEST_PATH_IMAGE089
上式中
Figure 592793DEST_PATH_IMAGE026
是离散空间序列
Figure 955641DEST_PATH_IMAGE071
的傅里叶变换；
步骤（3）、根据阵元数目
Figure 991730DEST_PATH_IMAGE034
及序列长度
Figure 24408DEST_PATH_IMAGE056
，选择DFT运算或FFT运算，若
Figure 145948DEST_PATH_IMAGE063
，选择FFT，否则，选择DFT；
实部方向图：
Figure 515750DEST_PATH_IMAGE090
Figure 203083DEST_PATH_IMAGE091
式中，
Figure 734558DEST_PATH_IMAGE074
为波束编号，
Figure 766231DEST_PATH_IMAGE075
；设
Figure 674144DEST_PATH_IMAGE064
的极大值对应的波束编号为
Figure 481563DEST_PATH_IMAGE076
， 则有
Figure 183940DEST_PATH_IMAGE092
；
由于
Figure 952175DEST_PATH_IMAGE027
关于
Figure 663779DEST_PATH_IMAGE032
的导数处处存在，因此，只要
Figure 60126DEST_PATH_IMAGE035
位于
Figure 933404DEST_PATH_IMAGE027
主瓣，就可以通过在小 范围内小步长细搜或解析求解的方法求得
Figure 516832DEST_PATH_IMAGE027
更精确的极大值位置
Figure 687919DEST_PATH_IMAGE065
；
步骤（4）、读取接收来波波长和接收通道输出
Figure 407613DEST_PATH_IMAGE054
，并更新离散空间序列Z m ，若m= p n , Z m =x n* conj(x 0)；
步骤（5）、基于步骤（3）的选择，对离散空间序列做DFT或FFT运算，获得实部方向图
Figure 514110DEST_PATH_IMAGE064
，并求得与
Figure 319255DEST_PATH_IMAGE064
最大值相对应的
Figure 44765DEST_PATH_IMAGE035
；
步骤（6）、基于
Figure 618966DEST_PATH_IMAGE035
进行迭代求解，获得使实部方向图
Figure 99626DEST_PATH_IMAGE027
最大的精确
Figure 188805DEST_PATH_IMAGE065
；
初始化
Figure 780323DEST_PATH_IMAGE078
；
Figure 894516DEST_PATH_IMAGE079
；
Figure 280498DEST_PATH_IMAGE080
；
Figure 122552DEST_PATH_IMAGE093
；
Figure 252182DEST_PATH_IMAGE082
；
步骤（7）、计算来波方向
Figure 473079DEST_PATH_IMAGE066
。
本实施例中，离散空间序列长度
Figure 295542DEST_PATH_IMAGE056
的确定原则为：确保序列
Figure 359313DEST_PATH_IMAGE094
的极大值
Figure 292634DEST_PATH_IMAGE095
位于实部方向图
Figure 695933DEST_PATH_IMAGE027
大于副瓣电平的主瓣区
Figure 876248DEST_PATH_IMAGE061
内。另一方面，阵列孔径
Figure 364998DEST_PATH_IMAGE096
相当于最快的空间频率，因此，根据采样理论
Figure 898747DEST_PATH_IMAGE097
应该小于
Figure 156553DEST_PATH_IMAGE098
，一般取
Figure 258501DEST_PATH_IMAGE056
=
Figure 968968DEST_PATH_IMAGE099
。
如图4所示，来波方向为0°时非均匀阵线阵实部方向图中，
Figure 306409DEST_PATH_IMAGE100
=186，
Figure 684301DEST_PATH_IMAGE099
=236，
Figure 442303DEST_PATH_IMAGE056
=256。此处
Figure 171225DEST_PATH_IMAGE056
取值
Figure 984460DEST_PATH_IMAGE099
即可。
另外，本实施例中若采用DFT计算
Figure 279175DEST_PATH_IMAGE064
需要约
Figure 519664DEST_PATH_IMAGE101
次实数运算；若采用FFT计算
Figure 407985DEST_PATH_IMAGE064
需要先求得
Figure 24911DEST_PATH_IMAGE102
，其计算量约
Figure 174133DEST_PATH_IMAGE103
次实数运算。因此，当阵元数目
Figure 851102DEST_PATH_IMAGE063
时，则适合利用FFT计算。如图4所示的12阵元非均匀阵列采用FFT计算
Figure 210408DEST_PATH_IMAGE064
其效率更高。
如图4所示，为
Figure 631025DEST_PATH_IMAGE001
=0°时的一个12阵元、阵元位于{0,3,7,13,19,23,29,37,43,47, 51,59}
Figure 838015DEST_PATH_IMAGE007
的非均匀线阵的实部方向图。
如图3所示，图3中给出来波方向为45°时，仿真得到的测向精度随信噪比的变化曲线。

Claims (4)

1.一种适用于非均匀阵列的快速测向方法，其特征在于：所述非均匀阵元包括N个阵元 {
Figure 755007DEST_PATH_IMAGE001
}，来波方向为
Figure 527791DEST_PATH_IMAGE002
，
Figure 802915DEST_PATH_IMAGE003
为第
Figure 367757DEST_PATH_IMAGE004
个阵元距最左边的参考阵元E 0 的以最短半波长 为单位的距离，
Figure 437344DEST_PATH_IMAGE005
为该非均匀阵列适应的最短波长；
Figure 963003DEST_PATH_IMAGE006
；
该非均匀阵列中每个阵元的接收信号经正交下变频后表示为：
Figure 776239DEST_PATH_IMAGE007
；
Figure 946320DEST_PATH_IMAGE008
；
上式中，
Figure 186808DEST_PATH_IMAGE009
为正整数，
Figure 199764DEST_PATH_IMAGE010
为接收信号的相位函数，
Figure 816690DEST_PATH_IMAGE011
为第
Figure 592010DEST_PATH_IMAGE004
个阵元输出的接收信号 分量，
Figure 268979DEST_PATH_IMAGE012
为噪声分量；且
Figure 503651DEST_PATH_IMAGE013
，
Figure 658689DEST_PATH_IMAGE014
；
Figure 803363DEST_PATH_IMAGE015
为当前来波波长；
Figure 182392DEST_PATH_IMAGE016
和
Figure 107622DEST_PATH_IMAGE017
分别是信号
Figure 66351DEST_PATH_IMAGE018
的实部和虚部；
Figure 314799DEST_PATH_IMAGE019
表示复数；
具体快速测向方法包括以下步骤：
步骤（1）、根据阵列参数确定离散空间序列的长度
Figure 333570DEST_PATH_IMAGE020
；
Figure 542835DEST_PATH_IMAGE021
；
Figure 39675DEST_PATH_IMAGE022
；
Figure 893362DEST_PATH_IMAGE023
为
Figure 83034DEST_PATH_IMAGE024
的最小整数，
Figure 779595DEST_PATH_IMAGE025
为高于副瓣电平的主瓣宽度；
步骤（2）、构造并初始化长度为
Figure 80126DEST_PATH_IMAGE020
的离散空间序列，
Figure 530263DEST_PATH_IMAGE026
；
步骤（3）、根据阵元数目
Figure 890837DEST_PATH_IMAGE027
及序列长度
Figure 74694DEST_PATH_IMAGE020
，来选择是使用离散傅里变换DFT还是选择快 速傅里叶变换FFT计算，若
Figure 913337DEST_PATH_IMAGE028
，选择FFT，否则，选择DFT；
步骤（4）、读取接收来波波长和接收通道输出
Figure 741615DEST_PATH_IMAGE018
，并更新离散空间序列Z m ，若m=p n , Z m =x n* conj(x 0)；
步骤（5）、根据步骤（3）的选择对离散空间序列Z m 做离散傅里变换DFT或快速傅里叶变换 FFT运算，获得实部方向图
Figure 273091DEST_PATH_IMAGE029
，并求得与
Figure 678664DEST_PATH_IMAGE029
最大值相对应的
Figure 586577DEST_PATH_IMAGE030
；
步骤（6）、基于
Figure 518630DEST_PATH_IMAGE030
进行迭代求解，获得使实部方向图
Figure 221007DEST_PATH_IMAGE031
最大的精确
Figure 113877DEST_PATH_IMAGE032
；
步骤（7）、计算来波方向
Figure 294322DEST_PATH_IMAGE033
。
2.根据权利要求1所述的适用于非均匀阵列的快速测向方法，其特征在于：
基于
Figure 97193DEST_PATH_IMAGE034
重新构造离散空间序列如下
Figure DEST_PATH_IMAGE035
式中，
Figure 767209DEST_PATH_IMAGE036
为
Figure 85058DEST_PATH_IMAGE037
的共轭，则
Figure 757610DEST_PATH_IMAGE038
Figure DEST_PATH_IMAGE039
Figure 274042DEST_PATH_IMAGE040
上式中
Figure DEST_PATH_IMAGE041
是离散空间序列
Figure 990325DEST_PATH_IMAGE042
的傅里叶变换。
3.根据权利要求1所述的适用于非均匀阵列的快速测向方法，其特征在于：通过DFT或 FFT实现
Figure 795470DEST_PATH_IMAGE031
的快速计算，如下
Figure DEST_PATH_IMAGE043
Figure 380035DEST_PATH_IMAGE044
；
式中，
Figure 141187DEST_PATH_IMAGE045
为波束编号，
Figure 356267DEST_PATH_IMAGE046
；设
Figure 445446DEST_PATH_IMAGE029
的极大值对应的波束编号为
Figure DEST_PATH_IMAGE047
，则有
Figure 709068DEST_PATH_IMAGE048
且
Figure 137776DEST_PATH_IMAGE030
位于
Figure 586074DEST_PATH_IMAGE031
主瓣。
4.根据权利要求1所述的适用于非均匀阵列的快速测向方法，其特征在于：所述步骤 （6）求解
Figure 365812DEST_PATH_IMAGE032
的过程为：
初始化
Figure 180928DEST_PATH_IMAGE049
；
Figure 464141DEST_PATH_IMAGE050
；
Figure 348921DEST_PATH_IMAGE051
；
Figure 350375DEST_PATH_IMAGE052
；
Figure 955800DEST_PATH_IMAGE053
。
CN202210433326.8A 2022-04-24 2022-04-24 一种适用于非均匀阵列的快速测向方法 Active CN114545324B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210433326.8A CN114545324B (zh) 2022-04-24 2022-04-24 一种适用于非均匀阵列的快速测向方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210433326.8A CN114545324B (zh) 2022-04-24 2022-04-24 一种适用于非均匀阵列的快速测向方法

Publications (2)

Publication Number Publication Date
CN114545324A true CN114545324A (zh) 2022-05-27
CN114545324B CN114545324B (zh) 2022-07-12

Family

ID=81667213

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210433326.8A Active CN114545324B (zh) 2022-04-24 2022-04-24 一种适用于非均匀阵列的快速测向方法

Country Status (1)

Country Link
CN (1) CN114545324B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101349742A (zh) * 2008-08-29 2009-01-21 西安电子科技大学 用于宽带相控阵阵元间距优化和频域多目标测频测向方法
CN101470187A (zh) * 2007-12-26 2009-07-01 中国科学院声学研究所 一种用于线列阵的高精度测向方法
CN103091677A (zh) * 2012-07-08 2013-05-08 西北工业大学 一种基于时反的均匀线列阵波束形成方法
CN109597021A (zh) * 2018-11-05 2019-04-09 北京航天控制仪器研究所 一种波达方向估计方法及装置
US20190212411A1 (en) * 2018-01-08 2019-07-11 Hyundai Mobis Co., Ltd. Method and apparatus for estimating direction of arrival using generation of virtual received signals
CN110325880A (zh) * 2016-12-09 2019-10-11 俄罗斯联邦诺萨顿国家原子能公司 用于搜索和探测伽马辐射源头的方法
CN110361691A (zh) * 2019-07-24 2019-10-22 哈尔滨工程大学 基于非均匀阵列的相干信源doa估计fpga实现方法
CN111413675A (zh) * 2020-04-10 2020-07-14 扬州宇安电子科技有限公司 一种信号采集设备及其采集方法
US20210270925A1 (en) * 2020-02-27 2021-09-02 The Boeing Company Systems and methods for enhanced direction of arrival detection and calculation

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101470187A (zh) * 2007-12-26 2009-07-01 中国科学院声学研究所 一种用于线列阵的高精度测向方法
CN101349742A (zh) * 2008-08-29 2009-01-21 西安电子科技大学 用于宽带相控阵阵元间距优化和频域多目标测频测向方法
CN103091677A (zh) * 2012-07-08 2013-05-08 西北工业大学 一种基于时反的均匀线列阵波束形成方法
CN110325880A (zh) * 2016-12-09 2019-10-11 俄罗斯联邦诺萨顿国家原子能公司 用于搜索和探测伽马辐射源头的方法
US20190212411A1 (en) * 2018-01-08 2019-07-11 Hyundai Mobis Co., Ltd. Method and apparatus for estimating direction of arrival using generation of virtual received signals
CN109597021A (zh) * 2018-11-05 2019-04-09 北京航天控制仪器研究所 一种波达方向估计方法及装置
CN110361691A (zh) * 2019-07-24 2019-10-22 哈尔滨工程大学 基于非均匀阵列的相干信源doa估计fpga实现方法
US20210270925A1 (en) * 2020-02-27 2021-09-02 The Boeing Company Systems and methods for enhanced direction of arrival detection and calculation
CN111413675A (zh) * 2020-04-10 2020-07-14 扬州宇安电子科技有限公司 一种信号采集设备及其采集方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
TAKAHASHI, HIRATA.: "An antenna arrangement for phase comparison monopulse DOA estimation using nonuniform planar array", 《 RADAR CONFERENCE. IEEE》 *
张庭兰: "非均匀线阵测向算法研究", 《万方数据库》 *
谢远举: "非均匀阵列测向技术研究与实现", 《万方数据库》 *

Also Published As

Publication number Publication date
CN114545324B (zh) 2022-07-12

Similar Documents

Publication Publication Date Title
CN104730491B (zh) 一种基于l型阵的虚拟阵列doa估计方法
CN106546948B (zh) 基于幅相误差阵列的远近场宽带混合源中远场源测向方法
CN109061554B (zh) 一种基于空间离散网格动态更新的目标到达角度估计方法
CN106501770B (zh) 基于幅相误差阵列的远近场宽带混合源中近场源定位方法
CN110007267B (zh) 一种基于混合基线的均匀圆阵干涉仪测向解模糊方法
CN112130111B (zh) 一种大规模均匀十字阵列中单快拍二维doa估计方法
CN108398669B (zh) 一种基于无需预延迟处理的空时宽带自适应单脉冲测角方法
CN111965591B (zh) 一种基于四阶累积量矢量化dft的测向估计方法
CN106707250B (zh) 基于互耦校正的雷达阵列自适应波束形成方法
CN112051540B (zh) 一种快速高精度测向方法
CN113189592A (zh) 考虑幅相互耦误差的车载毫米波mimo雷达测角方法
CN103323832A (zh) 一种相控阵三维摄像声纳系统换能器阵列的幅相误差校正方法
CN113835063B (zh) 一种无人机阵列幅相误差与信号doa联合估计方法
CN113671439B (zh) 基于非均匀智能超表面阵列的无人机集群测向系统及方法
CN108761380A (zh) 一种用于提高精度的目标波达方向估计方法
CN111368256B (zh) 一种基于均匀圆阵的单快拍测向方法
CN109541526A (zh) 一种利用矩阵变换的圆环阵方位估计方法
CN114545324B (zh) 一种适用于非均匀阵列的快速测向方法
CN114487992B (zh) 无须解模糊的多基线相位干涉仪测向方法
CN114563760B (zh) 一种基于sca阵型的二阶超波束形成方法、设备及介质
CN115826004A (zh) 一种基于二维角度及时差联合的三星协同直接定位方法
CN110531310A (zh) 基于子空间和内插变换的远场相干信号波达方向估计方法
CN109870670A (zh) 一种基于阵列重构的混合信号参数估计方法
CN114609580A (zh) 一种基于非圆信号的无孔互质阵列设计方法
CN111366891B (zh) 一种基于伪协方差矩阵的均匀圆阵单快拍测向方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 210000 building A1, Huangpu science and technology building, No. 2 Huangpu Road, Xuanwu District, Nanjing, Jiangsu Province

Patentee after: Nanjing Yu'an Defense Technology Co.,Ltd.

Country or region after: Zhong Guo

Patentee after: Yangzhou Yuan Electronic Technology Co.,Ltd.

Address before: 210000 building A1, Huangpu science and technology building, No. 2 Huangpu Road, Xuanwu District, Nanjing, Jiangsu Province

Patentee before: Nanjing Yu'an Defense Technology Co.,Ltd.

Country or region before: Zhong Guo

Patentee before: YANGZHOU YU'AN ELECTRONIC TECHNOLOGY CO.,LTD.