CN105137406B - 天基雷达效能优化方法 - Google Patents

天基雷达效能优化方法 Download PDF

Info

Publication number
CN105137406B
CN105137406B CN201510443121.8A CN201510443121A CN105137406B CN 105137406 B CN105137406 B CN 105137406B CN 201510443121 A CN201510443121 A CN 201510443121A CN 105137406 B CN105137406 B CN 105137406B
Authority
CN
China
Prior art keywords
mrow
msub
msup
mfrac
inclination angle
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
CN201510443121.8A
Other languages
English (en)
Other versions
CN105137406A (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.)
Shanghai Institute of Satellite Engineering
Original Assignee
Shanghai Institute of Satellite Engineering
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 Shanghai Institute of Satellite Engineering filed Critical Shanghai Institute of Satellite Engineering
Priority to CN201510443121.8A priority Critical patent/CN105137406B/zh
Publication of CN105137406A publication Critical patent/CN105137406A/zh
Application granted granted Critical
Publication of CN105137406B publication Critical patent/CN105137406B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/40Means for monitoring or calibrating

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)
  • Radio Relay Systems (AREA)

Abstract

本发明提供了一种天基雷达效能优化方法,包括以下:根据目标区域的纬度带和天基雷达的载荷作用距离确定轨道倾角;根据天基雷达的载荷作用距离和轨道倾角计算轨道高度;根据轨道倾角将多颗天基雷达卫星同轨道面内相位角均布或N个轨道面均布。本发明根据目标区域的纬度带,结合天基雷达的载荷作用距离,可选择0°轨道倾角或非0°轨道倾角,即低倾角轨道,可实现利用相对较少的卫星颗数,实现对目标区域的多次探测和高系统效能,以满足飞机目标的高时间分辨率探测需求。

Description

天基雷达效能优化方法
技术领域
本发明涉及卫星总体设计,具体地,涉及一种天基雷达效能优化方法。
背景技术
天基雷达是主动探测手段,具有雷达站位高,覆盖范围大,定位精度高,下视探测有利于反隐身等优势,填补天基对飞机探测手段的空白。
考虑飞机等目标为运动目标,具有运动速度快、机动能力强的特点,因此天基雷达需要具备较高的系统效能,即利用尽可能少的卫星颗数实现对关注区域的高时间分辨率重访,从而保证对飞机目标的连续探测。
经过对现有技术的检索,发现申请号为201210000464.3,发明名称为天基相控阵雷达空间多目标定轨方法的发明主要解决现有技术对空间弱目标不能有效估计和目标定轨精度低的问题。其过程为:对目标回波数据采用调零保形算法进行处理,运用距离脉冲压缩原理得到距离先验信息,并据此对回波信号进行分段处理;根据目标数将回波信号分段,每段数据为目标信号及其临近单元信号;对各段数据运用和差多波束测角原理实现多目标二维角度估计;基于空间目标跟踪结果和探测星轨道信息,进行坐标转换;采用Laplace型迭代算法实现对空间远近不同目标定轨,并运用最小二乘算法提高定轨精度。本发明能够降低强信号对弱目标的影响,对弱目标参数进行准确估计,可用于空间态势感知、轨道资源管理等实际应用领域。但是此发明不适用速度快机动能力强的目标。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种天基雷达效能优化方法。
根据本发明提供的天基雷达效能优化方法,包括以下步骤:
步骤1:根据目标区域的纬度带和天基雷达的载荷作用距离确定轨道倾角;
步骤2:根据天基雷达的载荷作用距离和轨道倾角计算轨道高度;
步骤3:根据轨道倾角将多颗天基雷达卫星同轨道面内相位角均布或多个轨道面均布。
优选地,在步骤1中,当目标区域在天基雷达的载荷作用距离内时,则轨道倾角为0°;当目标区域在天基雷达的载荷作用距离时,则轨道倾角为目标区域的纬度加上10°~15°。
优选地,在步骤2中,当轨道倾角为0°时,则根据载荷作用距离计算轨道高度h:
其中,l为载荷作用距离,Re为地球半径;
当轨道倾角不为0°时,则根据载荷作用距离和轨道倾角,计算得到轨道高度:
其中,Re为地球半径、μ为地球引力常数、ωe为地球自转角速度、J2为地球非球形带谐项、D为轨道回归周期、M为回归圈数、e为偏心率、h为轨道高度、i为轨道倾角。
优选地,在步骤3中,当轨道倾角为0°时,则将多颗天基雷达卫星同轨道面内相位角均布;当轨道倾角不为0°时,将多颗天基雷达卫星多个轨道面均布,即相邻轨道面天基雷达卫星的相位角差值相等。
优选地,当轨道倾角为0°时,则目标区域的最大重访时间Trev max为:
其中,N为天基雷达的卫星数量,Te为轨道交点周期。
优选地,当轨道倾角不为0°时,将N颗天基雷达卫星采用Walker星座N/N/F构型,其中N为2n,n为常数,n=1、2、3……,即N个轨道面均布,每个轨道面1颗卫星,相邻轨道面卫星相位角相差2πF/N,其中,F为常数,0≤F≤N-1,当回归圈数为奇数时,F为1;当回归圈数为偶数时,在n=1时,F为0、在n≠1时,F为2,则目标区域最大重访时间Trev max为:
与现有技术相比,本发明具有如下的有益效果:
本发明根据目标区域的纬度带,结合天基雷达的载荷作用距离,可选择0°轨道倾角或非0°轨道倾角,即低倾角轨道,可实现利用相对较少的卫星颗数,实现对目标区域的多次探测和高系统效能,以满足飞机目标的高时间分辨率探测需求。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明中目标区域和低倾角轨道星下点位置关系;
图2为本发明中低倾角轨道天基雷达卫星组网示意图;
图3为本发明中0°轨道倾角时轨道高度和载荷作用距离的变化曲线;
图4为本发明中0°轨道倾角时可探测最大纬度带和载荷作用距离的变化曲线;
图5为本发明中0°轨道倾角时1颗天基雷达卫星最大重访时间和载荷作用距离的变化曲线;
图6为本发明的步骤流程图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
本发明所要解决的技术问题是提供一种针对关注区域的天基雷达系统效能优化方法包括以下步骤:
步骤1:根据目标区域的纬度带和天基雷达的载荷作用距离确定轨道倾角;
步骤2:根据天基雷达的载荷作用距离和轨道倾角计算轨道高度;
步骤3:根据轨道倾角将多颗天基雷达卫星同轨道面内相位角均布或多个轨道面均布。
其中,步骤1中,轨道倾角应结合雷达载荷能力适当选择,当雷达作用距离可覆盖关注区域纬度时,选择0°倾角,除此之外,轨道倾角通常需要比关注区域纬度大10°~15°,具体数值由雷达作用距离决定。
其中,步骤2中,如果轨道倾角为0°,可根据载荷作用距离选择轨道高度h,由下式计算:
式中:l为雷达作用距离,Re为地球平均半径。
如果轨道倾角不为0°,优先选择1天回归轨道,结合雷达作用距离和轨道倾角,计算得到轨道高度,由下式计算:
式中Re为地球半径、μ为地球引力常数、ωe为地球自转角速度、J2为地球非球形带谐项、D为轨道回归周期、M为回归圈数、e为偏心率、h为轨道高度、i为轨道倾角,由步骤1决定,结合雷达作用距离,选择卫星运行15圈、14圈或13圈,近似1天地面轨迹重复,即D=1,M=15、14、13……。
其中,步骤中,如果轨道倾角为0°,通过N颗天基雷达同轨道面内相位角均布,可对纬度Lat范围内关注区域最大重访时间为:
式中Te为轨道交点周期。
如果轨道倾角不为0°,通过N颗天基雷达采用Walker星座N/N/F构型,其中N为2n,n=1、2、3……,即N个轨道面均布,每个轨道面1颗卫星,相邻轨道面卫星相位角相差2πF/N,其中F为常数,0≤F≤N-1,当回归圈数为奇数时,F为1;当回归圈数为偶数时,在n=1时,F为0、在n≠1时,F为2,使得卫星地面轨迹一致,可实现关注区域最大重访时间为:
如图1所示,以中国沿海某目标区域为例,假设雷达作用距离为2500km时,可选择轨道倾角为38°,轨道高度494km,15圈地面轨迹重复轨道,目标区域应在星下点菱形区域的中间位置,2颗卫星采用Walker2/2/1组网构型,可对该目标区域实现每天探测14次以上,最大重访时间小于2小时。
如图2所示,以中国沿海某目标区域为例,假设雷达作用距离为5000km时,可选择倾角为0°,轨道高度1728km,6颗星相位角均布,可对该目标区域实现每天探测66次,最大重访时间为22分钟(不考虑载荷探测时间)。
如图3所示,选择倾角为0°,雷达作用距离和轨道高度的关系,如雷达作用距离为5000km时,对应的轨道高度为1728km。
如图4所示,选择倾角为0°,雷达作用距离和可探测最大纬度的关系,如雷达作用距离为5000km时,对应的可探测最大纬度为38°。
如图5所示,选择倾角为0°,雷达作用距离和1颗星重访时间的关系,如雷达作用距离为5000km时,对应的1颗星重访时间为2.2小时,按照步骤(3)计算公式,6颗星组网时最大重访时间即为22分钟。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。

Claims (4)

1.一种天基雷达效能优化方法,其特征在于,包括以下步骤:
步骤1:根据目标区域的纬度带和天基雷达的载荷作用距离确定轨道倾角;
步骤2:根据天基雷达的载荷作用距离和轨道倾角计算轨道高度;
步骤3:根据轨道倾角将多颗天基雷达卫星同轨道面内相位角均布或多个轨道面均布;
在步骤1中,当目标区域在天基雷达的载荷作用距离内时,则轨道倾角为0°;当目标区域在天基雷达的载荷作用距离时,则轨道倾角为目标区域的纬度加上10°~15°;
在步骤2中,当轨道倾角为0°时,则根据载荷作用距离计算轨道高度h:
<mrow> <mi>h</mi> <mo>=</mo> <msqrt> <mrow> <msup> <mi>l</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>R</mi> <mi>e</mi> </msub> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <msub> <mi>R</mi> <mi>e</mi> </msub> </mrow>
其中,l为载荷作用距离,Re为地球半径;
当轨道倾角不为0°时,则根据载荷作用距离和轨道倾角,计算得到轨道高度:
<mrow> <mn>2</mn> <mi>&amp;pi;</mi> <mi>D</mi> <mo>=</mo> <mn>2</mn> <mi>&amp;pi;</mi> <msqrt> <mfrac> <msup> <mrow> <mo>(</mo> <mi>Re</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mn>3</mn> </msup> <mi>&amp;mu;</mi> </mfrac> </msqrt> <mo>/</mo> <mo>&amp;lsqb;</mo> <mn>1</mn> <mo>+</mo> <mfrac> <mn>3</mn> <mn>2</mn> </mfrac> <msub> <mi>J</mi> <mn>2</mn> </msub> <msup> <mrow> <mo>(</mo> <mfrac> <mi>Re</mi> <mrow> <mo>(</mo> <mi>Re</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mn>3</mn> </msup> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mn>3</mn> <mo>-</mo> <mn>4</mn> <msup> <mi>sin</mi> <mn>2</mn> </msup> <mi>i</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mrow> <mo>(</mo> <msub> <mi>&amp;omega;</mi> <mi>e</mi> </msub> <mo>+</mo> <mfrac> <mrow> <mn>3</mn> <msub> <mi>J</mi> <mn>2</mn> </msub> </mrow> <mn>2</mn> </mfrac> <msup> <mrow> <mo>(</mo> <mfrac> <mi>Re</mi> <mrow> <mo>(</mo> <mi>Re</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msqrt> <mfrac> <mi>&amp;mu;</mi> <msup> <mrow> <mo>(</mo> <mi>Re</mi> <mo>+</mo> <mi>h</mi> <mo>)</mo> </mrow> <mn>3</mn> </msup> </mfrac> </msqrt> <mi>cos</mi> <mi> </mi> <mi>i</mi> <mo>)</mo> </mrow> <mi>M</mi> </mrow>
其中,Re为地球半径、μ为地球引力常数、ωe为地球自转角速度、J2为地球非球形带谐项、D为轨道回归周期、M为回归圈数、e为偏心率、h为轨道高度、i为轨道倾角。
2.根据权利要求1所述的天基雷达效能优化方法,其特征在于,在步骤3中,当轨道倾角为0°时,则将多颗天基雷达卫星同轨道面内相位角均布;当轨道倾角不为0°时,将多颗天基雷达卫星多个轨道面均布,即相邻轨道面天基雷达卫星的相位角差值相等。
3.根据权利要求2所述的天基雷达效能优化方法,其特征在于,当轨道倾角为0°时,则目标区域的最大重访时间Trev max为:
<mrow> <msub> <mi>T</mi> <mrow> <mi>r</mi> <mi>e</mi> <mi>v</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <mi>&amp;pi;</mi> </mrow> <mrow> <mi>N</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>2</mn> <mi>&amp;pi;</mi> </mrow> <msub> <mi>T</mi> <mi>e</mi> </msub> </mfrac> <mo>-</mo> <msub> <mi>&amp;omega;</mi> <mi>e</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>
其中,N为天基雷达的卫星数量,Te为轨道交点周期。
4.根据权利要求3所述的基雷达效能优化方法,其特征在于,当轨道倾角不为0°时,将N颗天基雷达卫星采用Walker星座N/N/F构型,其中N为2n,n为常数,n=1、2、3……,即N个轨道面均布,每个轨道面1颗卫星,相邻轨道面卫星相位角相差2πF/N,其中,F为常数,0≤F≤N-1,当回归圈数为奇数时,F为1;当回归圈数为偶数时,在n=1时,F为0、在n≠1时,F为2,则目标区域最大重访时间Trev max为:
<mrow> <msub> <mi>T</mi> <mrow> <mi>r</mi> <mi>e</mi> <mi>v</mi> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <mi>&amp;pi;</mi> </mrow> <mrow> <mi>n</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>2</mn> <mi>&amp;pi;</mi> </mrow> <msub> <mi>T</mi> <mi>e</mi> </msub> </mfrac> <mo>-</mo> <msub> <mi>&amp;omega;</mi> <mi>e</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>.</mo> </mrow> 1
CN201510443121.8A 2015-07-24 2015-07-24 天基雷达效能优化方法 Active CN105137406B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510443121.8A CN105137406B (zh) 2015-07-24 2015-07-24 天基雷达效能优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510443121.8A CN105137406B (zh) 2015-07-24 2015-07-24 天基雷达效能优化方法

Publications (2)

Publication Number Publication Date
CN105137406A CN105137406A (zh) 2015-12-09
CN105137406B true CN105137406B (zh) 2017-12-01

Family

ID=54722812

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510443121.8A Active CN105137406B (zh) 2015-07-24 2015-07-24 天基雷达效能优化方法

Country Status (1)

Country Link
CN (1) CN105137406B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108572362B (zh) * 2018-01-30 2020-07-31 哈尔滨工业大学 一种tws雷达空时联合关联跟踪方法及装置
CN108469604B (zh) * 2018-02-07 2019-04-09 哈尔滨工业大学 Tws雷达空时联合数据模拟方法及系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7744037B2 (en) * 2007-04-11 2010-06-29 National Defense University Calculating method for deducing possibilities of all possible scenarios of satellite members in low earth satellite constellation
CN101614813B (zh) * 2009-07-23 2011-08-17 航天东方红卫星有限公司 全天候覆盖卫星的重访轨道确定方法
US20130062471A1 (en) * 2010-01-14 2013-03-14 Wah L. Lim Inclined orbit satellite communication system
CN103675832B (zh) * 2013-11-29 2015-11-25 中国空间技术研究院 一种快速重访离散目标的低轨卫星轨道设计方法
CN104038272B (zh) * 2014-06-10 2017-04-26 哈尔滨工业大学 一种光照约束下的中轨全球覆盖星座

Also Published As

Publication number Publication date
CN105137406A (zh) 2015-12-09

Similar Documents

Publication Publication Date Title
CN101887122B (zh) 可消除地面高程误差的星载sar图像目标定位方法
CN101403790B (zh) 单频gps接收机的精密单点定位方法
CN101968542B (zh) 一种利用地球站对月球探测器进行跟踪的方法
CN108181630B (zh) 一种北斗双天线旋转快速定向方法
CN102221688B (zh) 一种雷达系统误差估计方法
CN106405589A (zh) 一种全球电离层格网模型的确定方法及装置
CN102565834A (zh) 一种单频gps测向系统及其测向定位方法
CN103744052B (zh) 一种应用于空中目标定位的双星测时差测向方法及装置
CN107576939A (zh) 一种基于虚拟测距信标的单信标测距定位方法
CN105203023A (zh) 一种车载三维激光扫描系统安置参数的一站式标定方法
CN103033188A (zh) 基于综合孔径观测的导航卫星自主时间同步方法
CN103728647A (zh) 一种基于卫星载波信号调制的弹体滚转角测量方法
CN110672871B (zh) 一种基于组合惯导信息和光电转塔视频跟踪的运动目标测速测向方法
CN102749621A (zh) 一种双基地合成孔径雷达频域成像方法
CN102169173A (zh) 一种倾斜地球同步轨道合成孔径雷达模糊度分析方法
CN102736073B (zh) 一种通用模式下星载sar距离向模糊度的计算方法
CN106772326A (zh) 一种多子阵合成孔径声纳相位误差分析方法
CN105137406B (zh) 天基雷达效能优化方法
CN103713300B (zh) 一种准静态双星定位的方法及其应用
CN102508279A (zh) 卫星导航系统gnss定姿测量值的处理方法及gnss定姿测量仪
CN104792321A (zh) 一种基于辅助定位的土地信息采集系统及方法
Guo et al. Research on spatial-variant property of bistatic ISAR imaging plane of space target
CN116008945B (zh) 一种基于车载四面二维相控阵雷达航迹关联方法
CN1303431C (zh) 机载合成孔径雷达测量区域定位系统
CN110440791A (zh) 一种机载雷达rcs测量航线设计方法和装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant