CN110879382B - 具有地形自适应的天线波束空间配置方法 - Google Patents

具有地形自适应的天线波束空间配置方法 Download PDF

Info

Publication number
CN110879382B
CN110879382B CN201911302423.8A CN201911302423A CN110879382B CN 110879382 B CN110879382 B CN 110879382B CN 201911302423 A CN201911302423 A CN 201911302423A CN 110879382 B CN110879382 B CN 110879382B
Authority
CN
China
Prior art keywords
antenna
gamma
angle
lambda
alpha
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
Application number
CN201911302423.8A
Other languages
English (en)
Other versions
CN110879382A (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.)
Shaanxi Changling Electronic Technology Co ltd
Original Assignee
Shaanxi Changling 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 Shaanxi Changling Electronic Technology Co ltd filed Critical Shaanxi Changling Electronic Technology Co ltd
Priority to CN201911302423.8A priority Critical patent/CN110879382B/zh
Publication of CN110879382A publication Critical patent/CN110879382A/zh
Application granted granted Critical
Publication of CN110879382B publication Critical patent/CN110879382B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/03Details of HF subsystems specially adapted therefor, e.g. common to transmitter and receiver

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种具有地形自适应的天线波束空间配置方法,主要解决现有多普勒天线在水面等地形下回波信号的频谱发生明显畸变的问题。其实现方案是:确定天线波束中心下压角γ0和天线入射角ψ0的取值范围;计算波束中心方向在X‑Y面的投影与Y轴的夹角α0;选择γ0的值,确定下压角γ方向的天线方向性系数函数f(γ);根据选择的α0值旋转坐标系,计算新直角坐标系下的波束中心横向倾斜角λ′0,并确定新直角坐标系下横向倾斜角λ′方向的天线方向性系数函数g(λ′);根据f(γ)和g(λ′)确定完整的天线方向性系数函数F(γ,λ′);根据F(γ,λ′)和α0值设计多普勒雷达天线。本发明设计简单、地形自适应能力好,用于多普勒导航系统。

Description

具有地形自适应的天线波束空间配置方法
技术领域
本发明属于雷达技术领域,特别是一种天线波束空间配置方法,可用于多普勒雷达的天线设计。
背景技术
多普勒雷达是一种利用多普勒效应测量载体坐标系下三轴向速度的速度传感器。多普勒雷达通过天线发射和接收电磁波,根据回波信号的频谱计算接收和发射的电磁波频率之差,即多普勒频率,再根据多普勒频率计算得到载体和地面之间的相对运动速度值。
回波信号的频谱形状直接影响多普勒雷达测速精度,而多普勒雷达回波信号的频谱形状由天线方向性系数和地面散射系数共同决定。
天线波束空间配置会直接影响不同地形条件下回波信号的频谱形状,从而影响多普勒雷达的测速精度。现有的多普勒雷达天线波束空间配置方法为:在X-Y坐标系下,将波束方向与X轴的空间夹角设为下压角γ,将波束最大值方向与X轴的空间夹角设为波束中心下压角γ0,将波束方向与Y轴夹角的余角设为横向倾斜角λ,将波束最大值方向与Y轴夹角的余角设为波束中心横向倾斜角λ0,波束方向与Z轴的夹角设为入射角ψ。得到多普勒雷达天线方向性系数是下压角γ和横向倾斜角λ的函数:
Figure BDA0002322181280000021
而散射系数是入射角ψ的函数,所以天线方向性系数变化与散射系数变化的方向不一致。这种方法具有以下缺点:在散射系数随入射角的改变变化较大的地形下,尤其是水面,回波信号的频谱发生明显畸变,导致多普勒雷达的测速误差增大。过大的测速误差将会导致计算得出的多普勒雷达载体即时速度与实际速度存在偏差,计算得出的多普勒雷达载体位置与实际位置存在偏差,测速误差越大,偏差越大,这将导致多普勒雷达的载体失效、失控、偏离航线、损毁甚至失踪。因此,必须对多普勒雷达的测速误差进行修正。
现有的多普勒雷达采用软件修正的方法对回波信号的频谱畸变导致的测速误差进行修正。这种修正方法存在以下不足:第一是对海面畸变测速误差进行统一的固定值修正,而不同海情的海面其回波信号的频谱畸变并不相同,该修正方法的修正不准确;第二是仅针对海面畸变测速误差进行了修正,而在陆地的不同地形下如沙漠等,其回波信号的频谱仍可能存在较大的畸变,该修正方法无法对陆地地形导致的畸变测速误差进行修正;第三是需要进行陆海模式的手动切换,增加了操作人员的工作量,容易出现遗漏和误操作,造成修正不正确。
发明内容
本发明的目的在于针对上述现有技术的不足,提出一种具有地形自适应的天线波束空间配置方法,以减小多普勒雷达对所有地形下,尤其是海面的测速误差,保证载体的正常飞行。
本发明的技术思路是:通过旋转天线坐标系,使得旋转后坐标系上的波束中心横向倾斜角λ′0与入射角ψ0相等,设计天线方向图为下压角γ和旋转后坐标系上横向倾斜角λ′的函数,并使得天线方向图变化方向与散射系数变化方向接近一致,从而减少散射系数变化导致的回波信号的频谱畸变,达到地形自适应的目的。具体实现步骤包括如下:
1)设定天线波束中心下压角γ0的取值范围为Δγ±3°,其中Δγ的取值范围为68°~82°;
2)设定天线入射角ψ0的取值范围为Δψ±3°,其中Δψ的取值范围为18°~32°;
3)根据天线波束中心下压角γ0和天线入射角ψ0计算波束中心方向在X-Y面的投影与Y轴的夹角α0
Figure BDA0002322181280000031
4)选择满足上述条件的γ0和α0,使α0与γ0的差值最小;
5)根据γ0的值确定下压角γ方向的天线方向性系数函数为:
Figure BDA0002322181280000032
6)将直角坐标系X-Y绕Z轴沿顺时针方向旋转α0角度,得到新直角坐标系X'-Y';
7)根据α0的值计算新直角坐标系X'-Y'下的波束中心横向倾斜角λ′0
8)根据λ′0的值确定新直角坐标系X'-Y'下横向倾斜角λ′方向的天线方向性系数函数为:
Figure BDA0002322181280000041
9)根据下压角γ方向的天线方向性系数函数f(γ)和横向倾斜角λ′方向的天线方向性系数函数g(λ′),确定完整的天线方向性系数函数为:
Figure BDA0002322181280000042
10)设计多普勒雷达天线:
10a)根据α0的值确定天线的外形为锐角为α0的平行四边形;
10b)根据完整的天线方向性系数函数F(γ,λ′)和天线具体形式,确定天线辐射单元X方向的间距dx与Y方向的间距dy
10c)根据天线具体形式和天线辐射单元X方向的间距dx与Y方向的间距dy以及天线的外形,得到天线波束地面投影形状,完成天线波束的空间配置。
本发明具有如下优点:
本发明通过旋转坐标系,得到的天线方向性系数函数,该函数可以使λ′和ψ的值在天线波束中心范围内近似相等,天线的方向性系数近似为下压角γ和入射角ψ的函数,当散射系数随入射角ψ变化时,回波信号的频谱畸变很小,从而使得多普勒雷达的地形畸变测速误差很小。与现有技术相比本发明具有如下优点:
1.能实现多普勒雷达速度测量的各类地形自适应能力,无需软件进行速度修正,无需人员手动操作;
2.仅通过改变天线空间配置方法即可实现,无需增加额外硬件或软件,制作的天线外形及安装可与原有天线兼容。
附图说明
图1为本发明的实现流程图;
图2为现有的多普勒雷达天线四波束空间配置图;
图3为本发明的多普勒雷达天线四波束空间配置图;
图4为现有的多普勒雷达天线B波束地面投影图;
图5为本发明的多普勒雷达天线B波束地面投影图。
具体实施方式
以下参照附图对本发明做进一步说明。
参照图1,本发明的实现步骤如下:
步骤1.确定天线波束中心下压角γ0的取值范围。
由于多普勒雷达天线是利用多普勒效应来测量多普勒雷达载体的速度的,其波束方向不能垂直向下,也不能水平向前,而是在载体坐标系下与XYZ三个坐标轴均有一定的夹角才能测量载体的三轴向速度。一般要求载体的前向速度,即X方向速度测量更精确,而天线波束中心下压角γ0越小,对X方向速度测量更精确;另一方面,天线波束中心下压角γ0越小,回波信号越弱。因此,确定天线波束中心下压角γ0的取值范围通常为68°~82°,本实例在该范围内确定一个±3°的取值范围,即γ0选取71°±3°。
步骤2.确定天线入射角ψ0的取值范围。
与步骤1的原因相同,天线入射角ψ0的取值范围通常为18°~32°。本实例选择ψ0的取值范围为26°±3°。
步骤3.计算波束中心方向在X-Y面的投影与Y轴的夹角α0
参照图2的几何关系,根据
Figure BDA0002322181280000061
的计算公式,确定α0的变化范围为44.87°~50.60°。
步骤4.确定γ0和α0的值。
利用计算机按0.01°的精度计算取值γ0和ψ0范围内所有组合对应的α0,并计算所有组合下γ0和α0的差值;在所有的差值中,按γ0和α0的差值最小的原则,选择γ0的值为72°,ψ0的值为25.91°,计算出α0的值为45°,此时,γ0和α0的差值为27°,在所有组合中最小。
步骤5.确定下压角γ方向的天线方向性系数函数。
为了达到多普勒雷达天线的波束中心下压角γ0的值为72°的目的,需要构造天线下压角γ方向的天线方向性系数函数。如果天线不要求副瓣,则天线可用等幅分布实现,此时下压角γ方向的天线方向性系数函数为
Figure BDA0002322181280000071
该函数在γ=γ0时取得最大值。
如果天线对副瓣有要求,则需要在该函数上加窗函数,如泰勒窗函数或切比雪夫窗函数,可以有效降低天线的副瓣。
步骤6.构造新直角坐标系X'-Y'。
为了达到天线的方向性系数可近似为下压角γ和入射角ψ的函数的目的,需要将坐标系进行旋转,具体的做法是将直角坐标系X-Y绕Z轴沿顺时针方向旋转α0角度,得到新直角坐标系X'-Y'。
步骤7.计算波束中心横向倾斜角λ′0
参照图3的几何关系,按照波束中心横向倾斜角公式:λ′=sin-1(sinλ·sinα0+cosγ·cosα0),计算得到λ′0的值为25.91°,与ψ0的值相同,与设计思路吻合。
步骤8.确定横向倾斜角λ′方向的天线方向性系数函数。
为了达到多普勒雷达天线的横向倾斜角λ′0的值为25.91°的目的,需要构造天线下压角λ′0方向的天线方向性系数函数。
如果天线不要求副瓣,则天线可用等幅分布实现,此时下压角λ′方向的天线方向性系数函数为
Figure BDA0002322181280000081
该函数在λ′=λ′0时取得最大值;
如果天线对副瓣有要求,则需要在该函数上加窗函数,如泰勒窗函数或切比雪夫窗函数,可以有效降低天线的副瓣。
步骤9.确定完整的天线方向性系数函数。
由于下压角γ和横向倾斜角λ′是可分离变量,所以平面天线的方向性系数函数F(γ,λ′)为两个方向的方向性系数函数相乘,即:
Figure BDA0002322181280000082
步骤10.设计多普勒雷达天线。
10a)根据α0的值确定天线的外形为锐角为α0的平行四边形;
由于下压角γ和横向倾斜角λ′不正交,故天线外形需要布局为平行四边形,而平行四边形中的锐角与步骤6中坐标系旋转的角度α0相等,则平行四边形的锐角为α0
10b)确定天线辐射单元X方向的间距dx与Y方向的间距dy
具体的天线形式可以是多样的,常见的天线形式有波导裂缝阵列天线和微带阵列天线,而载体要求设备重量轻,本实例中采用微带阵列天线来实现。
根据完整的天线方向性系数函数F(γ,λ′),通过如下公式确定微带阵列天线辐射单元X方向的间距dx与Y方向的间距dy
Figure BDA0002322181280000091
Figure BDA0002322181280000092
其中,fr为天线的工作频率,c为电磁波在自由空间中的传播速度,εr为微带阵列天线材料的相对介电常数,m和n为空间波数,取值为整数。
10c)确定天线波束地面投影形状
微带阵列天线辐射单元X方向的间距dx与Y方向的间距dy保证了天线的波束最大值方向;而微带阵列天线的平行四边形外形,保证了波束中心横向倾斜角λ′0与入射角ψ0在同一平面且大小相等,根据确定的天线的波束最大值方向和波束中心横向倾斜角λ′0与入射角ψ0的几何关系,即可确定天线的空间波束形状,从而确定天线波束地面投影形状,如图5所示,从图5可见λ′和ψ0的曲线夹角很小,从而完成天线波束的空间配置。
将图5与图4所示的现有雷达天线波束的地面投影形状进行比较可见,由于现有雷达天线满足天线方向性系数函数
Figure BDA0002322181280000093
其λ和ψ0的曲线夹角很大,当散射系数随入射角ψ变化时,回波信号的频谱畸变很大,增加了多普勒雷达的地形畸变测速误差。本发明得到的天线波束地面投影形状,其λ′和ψ0的曲线夹角很小,当散射系数随入射角ψ变化时,回波信号的频谱畸变很小,减小了多普勒雷达因地形因素影响的测速误差。

Claims (3)

1.一种具有地形自适应的天线波束空间配置方法,包括:
1)设定天线波束中心下压角γ0的取值范围为Δγ±3°,其中Δγ的取值范围为68°~82°;
2)设定天线入射角ψ0的取值范围为Δψ±3°,其中Δψ的取值范围为18°~32°;
3)根据天线波束中心下压角γ0和天线入射角ψ0计算波束中心方向在X—Y面的投影与Y轴的夹角α0
Figure FDA0002322181270000011
4)选择满足上述条件的γ0和α0,使α0与γ0的差值最小;
5)根据γ0的值确定下压角γ方向的天线方向性系数函数为:
Figure FDA0002322181270000012
6)将直角坐标系X-Y绕Z轴沿顺时针方向旋转α0角度,得到新直角坐标系X’-Y’;
7)根据α0的值计算新直角坐标系X’-Y’下的波束中心横向倾斜角λ′0
8)根据λ′0的值确定新直角坐标系X’-Y’下横向倾斜角λ′方向的天线方向性系数函数为:
Figure FDA0002322181270000021
9)根据下压角γ方向的天线方向性系数函数f(γ)和横向倾斜角λ′方向的天线方向性系数函数g(λ′),确定完整的天线方向性系数函数为:
Figure FDA0002322181270000022
10)设计多普勒天线:
10a)根据α0的值确定天线的外形为锐角为α0的平行四边形;
10b)根据完整的天线方向性系数函数F(γ,λ′)和天线具体形式,确定天线辐射单元X方向的间距dx与Y方向的间距dy
10c)根据天线具体形式和天线辐射单元X方向的间距dx与Y方向的间距dy以及天线的外形,得到天线波束地面投影形状,完成天线波束的空间配置。
2.根据权利要求1所述的方法,其中步骤7)中的波束中心横向倾斜角λ′0,通过如下公式计算:
Figure FDA0002322181270000023
3.根据权利要求1所述的方法,其中步骤10b)中确定天线辐射单元X方向的间距dx与Y方向的间距dy,分别由如下公式确定:
Figure FDA0002322181270000031
Figure FDA0002322181270000032
其中,fr为天线的工作频率,c为电磁波在自由空间中的传播速度,εr为微带天线的相对介电常数,m和n为空间波数,取值为整数。
CN201911302423.8A 2019-12-17 2019-12-17 具有地形自适应的天线波束空间配置方法 Active CN110879382B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911302423.8A CN110879382B (zh) 2019-12-17 2019-12-17 具有地形自适应的天线波束空间配置方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911302423.8A CN110879382B (zh) 2019-12-17 2019-12-17 具有地形自适应的天线波束空间配置方法

Publications (2)

Publication Number Publication Date
CN110879382A CN110879382A (zh) 2020-03-13
CN110879382B true CN110879382B (zh) 2023-03-28

Family

ID=69730971

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911302423.8A Active CN110879382B (zh) 2019-12-17 2019-12-17 具有地形自适应的天线波束空间配置方法

Country Status (1)

Country Link
CN (1) CN110879382B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001228243A (ja) * 2000-02-18 2001-08-24 Mitsubishi Electric Corp ドップラーレーダ装置及びそのビーム走査方法
FR2809186A1 (fr) * 2000-05-22 2001-11-23 Celine Corbrion Procede et dispositif pour mesurer la vitesse d'un mobile
CN102854496A (zh) * 2012-09-19 2013-01-02 中国民航大学 基于双门限控制的机载气象雷达地杂波抑制方法
CN108549079A (zh) * 2018-04-13 2018-09-18 西安电子科技大学 三天线干涉延迟多普勒雷达高度表地形特显点定位方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001228243A (ja) * 2000-02-18 2001-08-24 Mitsubishi Electric Corp ドップラーレーダ装置及びそのビーム走査方法
FR2809186A1 (fr) * 2000-05-22 2001-11-23 Celine Corbrion Procede et dispositif pour mesurer la vitesse d'un mobile
CN102854496A (zh) * 2012-09-19 2013-01-02 中国民航大学 基于双门限控制的机载气象雷达地杂波抑制方法
CN108549079A (zh) * 2018-04-13 2018-09-18 西安电子科技大学 三天线干涉延迟多普勒雷达高度表地形特显点定位方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
机车多普勒测速雷达波束方向校正方法;孙多等;《火控雷达技术》(第01期);全文 *
机载雷达多普勒波束锐化算法改进;杨波;《现代雷达》(第11期);全文 *

Also Published As

Publication number Publication date
CN110879382A (zh) 2020-03-13

Similar Documents

Publication Publication Date Title
CN109669178A (zh) 一种星载三阵元单脉冲两维测向方法
US7136012B2 (en) Approach radar with array antenna having rows and columns skewed relative to the horizontal
CN106654507B (zh) 基于宽带双极化天线阵列的相位干涉仪的测向方法
CN111830332B (zh) 全空域曲面阵多标校区域联合标校方法
CN110518366B (zh) 基于3d打印技术的圆极化介质反射单元及反射阵天线
CN113239524B (zh) 基于移动旋转参考系仿真计算电扫阵列天线相位中心方法
CN115061156A (zh) 一种基于组合导航的阵列天线抗卫星导航欺骗方法及系统
CN111624565A (zh) 大型共形相控阵系统多区域联合波束跟踪方法
CN113884977B (zh) 一种对一维干涉仪测向圆锥效应的校正方法
CN114996965B (zh) 一种小口径低频段阵列搭建及测向校准方法
CN115360518B (zh) 一种相控阵天线的通道幅相标校方法
CN114578349A (zh) 用于堤坝隐患检测的无人机载地质雷达系统与巡检方法
CN108427104A (zh) 一种全极化多输入多输出合成孔径雷达(fp-mimo-sar)的极化校准方法
Ryzhikov et al. Measurement of angular coordinates of point targets in the onboard weather navigation radar based on a multi-channel phased antenna array with an assimetic pattern
CN111736120B (zh) 一种基于天波传播校正源信号的阵列误差校正方法
CN110879382B (zh) 具有地形自适应的天线波束空间配置方法
CN115113156B (zh) 用于双极化相控阵气象雷达的标校方法和系统
CN107015065B (zh) 窄波束天线电轴、相位中心和时延的远场联合标定方法
RU136588U1 (ru) Устройство для измерения малых углов места целей в радиолокационных станциях
CN112305512B (zh) 一种双频干涉成像高度计高精度基线估计方法
CN115618650B (zh) 一种修正球面相控阵天线的虚拟球心位置坐标的方法
CN118534447B (zh) 基于超短基线十阵元声基阵及其定位方法
CN115549821B (zh) 一种阵列天线外场差分相位差标定方法
Matvienko Estimation of the Practically Attainable Accuracy of Modern Ultrashort Baseline Hydroacoustic Navigation Systems for Underwater Robots
CN116930862B (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