CN111693979B - 一种基于泰勒展开的数字阵列单脉冲测角方法 - Google Patents
一种基于泰勒展开的数字阵列单脉冲测角方法 Download PDFInfo
- Publication number
- CN111693979B CN111693979B CN202010446467.4A CN202010446467A CN111693979B CN 111693979 B CN111693979 B CN 111693979B CN 202010446467 A CN202010446467 A CN 202010446467A CN 111693979 B CN111693979 B CN 111693979B
- Authority
- CN
- China
- Prior art keywords
- array
- vector
- angle
- theta
- signal
- 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
- 238000000691 measurement method Methods 0.000 title claims abstract description 11
- 239000013598 vector Substances 0.000 claims abstract description 46
- 238000012545 processing Methods 0.000 claims abstract description 10
- 238000003491 array Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 17
- 238000005457 optimization Methods 0.000 abstract 1
- 238000000034 method Methods 0.000 description 18
- 238000010586 diagram Methods 0.000 description 8
- 238000001514 detection method Methods 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
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
- 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
- G01S13/42—Simultaneous measurement of distance and other co-ordinates
- G01S13/44—Monopulse radar, i.e. simultaneous lobing
- G01S13/4463—Monopulse radar, i.e. simultaneous lobing using phased arrays
-
- 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
- G01S2013/0236—Special technical features
- G01S2013/0245—Radar with phased array antenna
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于泰勒展开的数字阵列单脉冲测角方法，该方法充分发挥数字阵列化处理方式灵活，能够对各个阵元的输出进行任意加权处理的特点，基于泰勒展开原理对和、差波束加权矢量进行优化设计，从而获得比传统单脉冲测角方法更加陡峭的鉴角曲线(单脉冲比)，有效提升了目标角度测量精度。
Description
技术领域
本发明属于角度测量技术领域，具体涉及一种基于泰勒展开的数字阵列单脉冲测角方法。
背景技术
单脉冲角度测量技术也被称为同时波束比较测角技术，是一种对辐射源或反射电磁能量的“目标”的角度定位技术，该技术利用同时产生的“和”、“差”波束来测量波束宽度范围内的目标角度，已经广泛应用于雷达系统中。单脉冲测角精度由单脉冲和差比的均方根误差推导得出，分别与回波信号的信噪比、3dB波束宽度以及归一化的单脉冲和差比斜率有关。在相控阵雷达设计过程中，综合考虑上述三个参数雷达将获得最佳的单脉冲测角性能，其中，回波信号的信噪比和3dB波束宽度与和波束方向图有关，单脉冲和差比斜率则分别与和波束方向图和差波束方向图有关。因此，为了获得最佳的单脉冲测角性能，需要对和、差波束方向图进行优化。
传统的相控阵雷达通过单脉冲和差网络来获得和差波束：各个子阵的输出相加形成和波束，左、右(或上、下)象限阵列输出相减形成差波束。显然，单脉冲和差网络限定了和、差方向图的优化范围。与此同时，和波束方向图与雷达探测性能有关，为了保证雷达探测性能，由天线口径决定的和波束方向图往往保持不变。因此，对于传统相控阵雷达而言，通过优化差波束方向图带来的单脉冲测角性能的提升十分有限。
现如今，先进的数字化阵列雷达通过数字波束形成技术替代单脉冲和差网络来获得和差波束，数字化处理方式对和、差波束方向图优化提供了极大的自由度，允许对各个子阵的输出进行数字加权处理，灵活地形成和、差波束，使得雷达系统能够在保证探测性能的同时，获得最佳的单脉冲测角精度。因此，如何优化单脉冲和、差方向图，使得数字阵列获得最佳的单脉冲测角性能需要进一步研究。
发明内容
有鉴于此，本发明的目的是提供一种于泰勒展开的数字阵列单脉冲测角方法，可以减小测角误差，提高测角精度。
一种泰勒展开的数字阵列单脉冲测角方法，包括如下步骤：
步骤1、假设相控阵雷达为一维线阵，天线单元非均匀排布，且所有天线单元均为各向同性阵元；
步骤2、假设远场处某个角度θt下有一期望信号st(t)以平面波入射，获得天线阵列的接收信号：
x＝a(θt)st+n (1)
其中，x表示阵列接收到的回波数据复向量，a(θt)表示阵列导向矢量，n表示高斯白噪声复
向量；
步骤3、导向矢量泰勒展开，具体为：
将角度θt与波束指向θ0差值近似为零，根据泰勒展开原理，在雷达波束指向角度θ＝θ0处完成导向矢量a(θt)的泰勒展开，忽略高阶项得：
x＝(a(θ0)+b(θ0)(θt-θ0))st+n (2)
其中，
选择某一加权矢量w对信号进行加权求和处理，可得：
yD表示阵列加权求和后的输出信号；为简化运算，要求该加权矢量w与波束中心的导向矢量
a(θ0)正交，即wHa(θ0)＝0，忽略噪声项的影响，式(3)化简为：
yD＝wHb(θ0)(θt-θ0)st (4)
步骤4、权矢量的选取，具体为：
从提升目标角度估计精度的角度出发，选取阵列加权矢量w的原则是使阵列加权求和后的输出信号yD的信噪比最大，求得权矢量为：
步骤5、估计期望信号的信号幅度参量，具体为：
未知的信号幅度参量st通过波束形成和通道信号数据来估计，和通路信号数据写作：
其中，aH(θ0)a(θt)是天线波束指向θ0角度时，天线阵列在θt角度下的响应，近似为θ0角度下的响应，即aH(θ0)a(θt)≈aH(θ0)a(θ0)；
忽略噪声项的影响，则式(6)化简为：
yS≈aH(θ0)a(θ0)st (7)
则幅度参量st的估计值写作：
步骤6、估计出目标角度，具体为：
确定阵列加权矢量w，并且估计出信号幅度参量st之后，将式(5)和式(8)代入式(4)，得：
令
由上式可估计出目标角度为：
本发明具有如下有益效果：
本发明提出的一种基于泰勒展开的数字阵列单脉冲测角方法，该方法充分发挥数字阵列化处理方式灵活，能够对各个阵元的输出进行任意加权处理的特点，基于泰勒展开原理对和、差波束加权矢量进行优化设计，从而获得比传统单脉冲测角方法更加陡峭的鉴角曲线(单脉冲比)，有效提升了目标角度测量精度。
附图说明
图1为本发明实施方式的信号处理流程图。
图2(a)和图2(b)分别为本发明方法与传统方法的和差方向图比较。
图3为本发明方法与传统方法的鉴角曲线比较。
具体实施方式
下面结合附图并举实施例，对本发明进行详细描述。
步骤1、天线阵列相对位置矢量计算。假设相控阵雷达为一维线阵，天线单元数为M，天线单元非均匀排布，且所有天线单元均为各向同性阵元。第m个天线单元的相位中心记为dm，称向量
d＝[d1,d2,…,dm,…,dM-1,dM]T (1)
为天线单元的绝对位置矢量。天线所有单元相位中心的平均值为
称为天线阵列的相位中心。天线的各单元相位中心相对于阵列相位中心的偏移量
组成的向量
称为天线单元的相对位置矢量。
步骤2、获得目标回波信号。假设远场处某个角度θt下有一期望信号st(t)以平面波入射，第m个天线单元接收信号为
xm(t)＝st(t-τm(θt))+nm(t)m＝1,2,…M (4)
其中，τm(θt)是信号到达第m个天线单元相对于天线阵列相位中心的时延，nm(t)是第m个天线单元接收的高斯白噪声。
假设期望信号st为窄带信号，则时间延迟可以近似为期望信号中心频点的相位偏移，式(14)改写为：
xm(t)＝am(θt)st(t)+nm(t)m＝1,2,…M (5)
其中，
c表示光速，f0表示期望信号的中心频率。
令
其中，x表示阵列接收到的回波数据复向量，a(θt)表示阵列导向矢量，n表示高斯白噪声复向量，阵列接收的窄带信号回波数据可表示为
x＝a(θt)st+n (7)
步骤3、导向矢量泰勒展开。根据阵列雷达工作原理可知，待测目标角度位于雷达探测波束主瓣范围内，因此，目标角度θt与波束指向θ0差值较小，可以近似为(θt-θ0)→0。根据泰勒展开原理，在雷达波束指向角度θ＝θ0处完成导向矢量a(θt)的泰勒展开，忽略高阶项可得：
a(θt)≈a(θ0)+b(θ0)(θt-θ0) (8)
其中，
代入阵列雷达接收的数据表达式(17)可得
x＝(a(θ0)+b(θ0)(θt-θ0))st+n (9)
选择某一加权矢量w对信号进行加权求和处理，可得
yD表示阵列加权求和后的输出信号。为简化运算，要求该加权矢量w与波束中心的导向矢量a(θ0)正交，即wHa(θ0)＝0，在信噪比足够的条件下可以忽略噪声项的影响，式(10)可以化简为
yD＝wHb(θ0)(θt-θ0)st (11)
因此，只需合理选取权矢量w，并估计出信号的复包络st，就可以从上式中估计出目标角度θt。
步骤4、权矢量的选取。从提升目标角度估计精度的角度出发，选取阵列加权矢量w的原则是使阵列加权求和后的输出信号yD的信噪比最大，求得权矢量为
步骤5、估计期望信号的复包络。未知的信号幅度参量st可以通过波束形成和通道信号数据来估计，和通路信号数据写作
其中，aH(θ0)a(θt)是天线波束指向θ0角度时，天线阵列在θt角度下的响应，可以近似为θ0角度下的响应，即aH(θ0)a(θt)≈aH(θ0)a(θ0)。
在信噪比足够的条件下，忽略噪声项的影响，则式(13)可以化简为yS≈aH(θ0)a(θ0)st (14)
则幅度参量st的估计值可以写作
步骤6、估计出目标角度。确定阵列加权矢量w，并且估计出信号幅度参量st之后，将式(12)和式(15)代入式(11)，可得：
令
由上式可估计出目标角度为：
下面给出应用本方法的仿真实例，具体仿真参数如表1所示。
表1均匀阵列仿真参数
传统差和比法和所提方法的和差天线方向图如图2所示。从图中可以看出，两种方法选用了相同的和通路加权矢量，因此，和方向图保持一致。两种方法选用了不同的差通路加权矢量，因此，差方向图存在差异。其中，所提方法的差方向图在和方向图主瓣范围内的取值比传统差和比法的差方向图要高，意味着所提方法拥有更陡的鉴角曲线。
两种方法在-1°～1°范围内获得的鉴角曲线(单脉冲比)如图3所示，从图中明显看出，所提方法获得的鉴角曲线比传统差和比法获得的鉴角曲线要陡。测角精度与单脉冲比K成反比，单脉冲比越大，测角精度越高。因此，所提方法拥有比传统方法更高的测量精度。
本发明的方法适用于数字阵列，是一种精确测量目标角度的方法，与传统方法相比，减小测角误差，提高测角精度。此外，测角方法对天线阵元排布没有约束，适用非均匀阵列。
综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (1)
1.一种泰勒展开的数字阵列单脉冲测角方法，其特征在于，包括如下步骤：
步骤1、假设相控阵雷达为一维线阵，天线单元非均匀排布，且所有天线单元均为各向同性阵元；
步骤2、假设远场处某个角度θt下有一期望信号st以平面波入射，获得天线阵列的接收信号：
x＝a(θt)st+n (1)
其中，x表示阵列接收到的回波数据复向量，a(θt)表示阵列导向矢量，n表示高斯白噪声复向量；
步骤3、导向矢量泰勒展开，具体为：
将角度θt与波束指向θ0差值近似为零，根据泰勒展开原理，在雷达波束指向角度θ＝θ0处完成导向矢量a(θt)的泰勒展开，忽略高阶项得：
x＝(a(θ0)+b(θ0)(θt-θ0))st+n (2)
其中，
选择某一加权矢量w对信号进行加权求和处理，可得：
yD表示阵列加权求和后的输出信号；为简化运算，要求该加权矢量w与波束中心的导向矢量a(θ0)正交，即wHa(θ0)＝0，忽略噪声项的影响，式(3)化简为：
yD＝wHb(θ0)(θt-θ0)st (4)
步骤4、权矢量的选取，具体为：
从提升目标角度估计精度的角度出发，选取阵列加权矢量w的原则是使阵列加权求和后的输出信号yD的信噪比最大，求得权矢量为：
步骤5、估计期望信号，具体为：
未知的期望信号st通过波束形成和通道信号数据来估计，和通路信号数据写作：
其中，aH(θ0)a(θt)是天线波束指向θ0角度时，天线阵列在θt角度下的响应，近似为θ0角度下的响应，即aH(θ0)a(θt)≈aH(θ0)a(θ0)；
忽略噪声项的影响，则式(6)化简为：
yS≈aH(θ0)a(θ0)st (7)
则期望信号st的估计值写作：
步骤6、估计出目标角度，具体为：
确定阵列加权矢量w，并且估计出期望信号st之后，将式(5)和式(8)代入式(4)，得：
令
由上式可估计出目标角度为：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010446467.4A CN111693979B (zh) | 2020-05-25 | 2020-05-25 | 一种基于泰勒展开的数字阵列单脉冲测角方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010446467.4A CN111693979B (zh) | 2020-05-25 | 2020-05-25 | 一种基于泰勒展开的数字阵列单脉冲测角方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111693979A CN111693979A (zh) | 2020-09-22 |
| CN111693979B true CN111693979B (zh) | 2022-10-11 |
Family
ID=72477537
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010446467.4A Active CN111693979B (zh) | 2020-05-25 | 2020-05-25 | 一种基于泰勒展开的数字阵列单脉冲测角方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111693979B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112782645B (zh) * | 2020-12-30 | 2021-11-19 | 无锡国芯微电子系统有限公司 | 一种四臂螺旋天线数据拟合测角方法 |
| CN113076670B (zh) * | 2021-03-26 | 2022-10-18 | 贵州航天电子科技有限公司 | 一种多物理场耦合的相控阵制导微系统协同优化方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5473331A (en) * | 1994-10-31 | 1995-12-05 | Hughes Aircraft Company | Combined SAR monopulse and inverse monopulse weapon guidance |
| CN101109810B (zh) * | 2007-08-24 | 2010-12-29 | 重庆大学 | 提高数字单脉冲系统跟踪精度的方法 |
| CN103728614B (zh) * | 2014-01-15 | 2015-10-28 | 西安电子科技大学 | 基于机扫米波雷达的改进单脉冲测角方法 |
| CN104166136A (zh) * | 2014-07-11 | 2014-11-26 | 河海大学 | 一种基于干扰子空间跟踪的高效自适应单脉冲测角方法 |
| CN108508423B (zh) * | 2018-01-25 | 2021-07-06 | 西安电子科技大学 | 基于异型阵的子阵数字和差单脉冲测角方法 |
| CN108398669B (zh) * | 2018-01-25 | 2021-08-06 | 电子科技大学 | 一种基于无需预延迟处理的空时宽带自适应单脉冲测角方法 |
-
2020
- 2020-05-25 CN CN202010446467.4A patent/CN111693979B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN111693979A (zh) | 2020-09-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108828522B (zh) | 一种利用垂直阵lcmv波束形成的水下目标辐射噪声测量方法 | |
| CN110488255B (zh) | 一种相控阵雷达单脉冲高分辨测角系统及方法 | |
| CN109946664B (zh) | 一种主瓣干扰下的阵列雷达导引头单脉冲测角方法 | |
| CN103235292B (zh) | 平面相控阵调零保形校正的全维和差测角方法 | |
| CN104237844B (zh) | 基于相位补偿的分布式米波平面阵列雷达方位角测量方法 | |
| RU2291464C2 (ru) | Способ измерения угла места целей при наличии отражений принимаемого эхосигнала от земной поверхности и импульсная наземная трехкоординатная радиолокационная станция для его реализации | |
| CN110673086A (zh) | 一种基于数字阵列雷达的二维角度超分辨方法 | |
| CN111239677B (zh) | 基于数字阵列的多波束被动单脉冲测角方法 | |
| CN110045334B (zh) | 低副瓣零陷波束形成方法 | |
| CN111693979B (zh) | 一种基于泰勒展开的数字阵列单脉冲测角方法 | |
| CN112147593B (zh) | 一种高速密集爆炸物破片目标的四维参数估计方法 | |
| CN110426670B (zh) | 基于tls-cs的外辐射源雷达超分辨doa估计方法 | |
| CN112881975B (zh) | 基于子阵特征矩阵联合对角化的单脉冲和差波束测角方法 | |
| CN110261826A (zh) | 一种零陷展宽的相干干扰抑制方法 | |
| CN109932679B (zh) | 一种传感器列系统最大似然角度分辨率估计方法 | |
| CN110196417B (zh) | 基于发射能量集中的双基地mimo雷达角度估计方法 | |
| US5528554A (en) | Linear array lateral motion compensation method | |
| CN109901131B (zh) | 基于斜投影的多径利用相干波束形成方法 | |
| CN109633563B (zh) | 基于多径信息的自适应相干波束形成方法 | |
| CN111812607A (zh) | 基于波束空间的米波mimo雷达低仰角估计方法 | |
| CN112965028B (zh) | 多波束相控阵差和比角度估计方法 | |
| Wu et al. | Array pattern synthesis and robust beamforming for a complex sonar system | |
| CN112711018B (zh) | 一种近场空间网状结构聚焦波束形成方法 | |
| Zhou et al. | Research on interferometer direction finding technology based on digital beam forming | |
| De et al. | Angle estimation using modified subarray level monopulse ratio algorithm and s-curve in digital phased array radar |
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 |