CN109840360A - 一种卫星临边探测方式下星座最小规模设计方法 - Google Patents
一种卫星临边探测方式下星座最小规模设计方法 Download PDFInfo
- Publication number
- CN109840360A CN109840360A CN201910037300.XA CN201910037300A CN109840360A CN 109840360 A CN109840360 A CN 109840360A CN 201910037300 A CN201910037300 A CN 201910037300A CN 109840360 A CN109840360 A CN 109840360A
- Authority
- CN
- China
- Prior art keywords
- satellite
- airspace
- constellation
- geocentric angle
- latitude
- 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
Classifications
-
- 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
Landscapes
- Instructional Devices (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种卫星临边探测方式下星座最小规模设计方法,包括如下步骤:步骤A:计算临边探测方式下单星覆盖空域对应的地心角和单星不能覆盖的空白空域对应的地心角,步骤B:计算填补空白空域需要的卫星数量,构成一组无覆盖空白空域的卫星,步骤C:代入地理纬度要求计算这一组卫星覆盖空域对应的地心角,步骤D:对照全球360°全覆盖的要求计算需要的总的卫星数量。采用本发明确定的卫星星座数量,能够适应全球不同地理纬度带不同海拔高度以上空域全覆盖观测需求,确保了星座满足应用要求的同时使用的卫星数量最少,能够显著优化卫星星座方案。
Description
技术领域
本发明涉及卫星星座设计领域,还涉及卫星轨道设计领域。
背景技术
随着卫星技术的不断发展,卫星服务于国民经济和国防建设的能力越来越强,各行各业对卫星的任务功能要求也日趋复杂,往往需要多颗卫星组网协同配合实现既定功能。当前大多数卫星均为直接对地探测方式,在装载成像相机或合成孔径雷达等有效载荷时,卫星对地面的可探测区域直观,对于覆盖探测需求设计方法已较为成熟。但是临边探测方式与直接对地探测方式不同,是对关注区域上空空域的探测,使得卫星探测时覆盖区域的计算方法与直接对地探测有所不同,因此,需要改变以往的设计方法,实现卫星星座规模优化设计。
发明内容
本发明针对卫星采用临边探测方式,对全球指定纬度区域的一定高度以上空域进行信息获取的需求,提出了一种确定卫星最小数量规模的设计方法,该方法能够满足某一地理纬度带全覆盖的要求,而且能够应用于任意海拔高度空域的探测任务,通过不同卫星轨道高度的迭代计算,能够得到随卫星轨道高度不同卫星星座最小规模的变化曲线,便于开展卫星星座的优化设计。目前没有发现同本发明类似的技术说明或报道。
本发明具体通过以下方案实现:
一种卫星临边探测方式下星座最小规模设计方法,适用于在任意地理纬度和海拔高度空域要求下的卫星星座最小数量设计,采用多颗卫星视场拼接的方式实现全球一定地理纬度空域的全覆盖,并且对拼接时是否满足地理纬度要求进行了有效性判断,星座内所有卫星的轨道高度一致,且均采用0°倾角轨道,通过代入不同轨道高度值迭代计算得到不同轨道高度下卫星星座最小规模;具体包括如下步骤:
步骤A:计算临边探测方式下单星覆盖空域对应的地心角和单星不能覆盖的空白空域对应的地心角;
步骤B:计算填补空白空域需要的卫星数量,构成一组无覆盖空白空域的卫星;
步骤C:代入地理纬度要求计算这一组卫星覆盖空域对应的地心角;
步骤D:对照全球360°全覆盖的要求计算需要的总的卫星数量。
所述步骤A中通过以下公式计算单星覆盖空域对应的地心角α:
式中,ψ2为观测空域高度对应的地心角的余角,R为地球半径,h为需要观测的空域高度;
单星未覆盖空域对应的地心角β:
β=2(α1-α)。
式中,α1为单星临边探测可覆盖最高纬度值。
所述步骤B分为两种类型:
类型一:当最高纬度值α1小于90°,即覆盖不超过极点的情况
采用卫星视场紧密衔接部署方式后,所能够覆盖的地理纬度为gama:
式中,R为地球半径,α1为单星临边探测可覆盖最高纬度值,α为单星覆盖空域对应的地心角;
若多星视场拼接后能够覆盖的地理纬度gama能够满足需要的覆盖的地理纬度X,即gama>X时,可证明可通过多星拼接满足要求;为了填补单星的未覆盖空域,需要填补的卫星数量L为:L=[β/α]+1,其中β为单星未覆盖空域对应的地心角;
因此,L颗卫星加上原本的1颗卫星构成无空白空域的卫星组,卫星组的卫星数量M为:M=L+1=[β/α]+2;
类型二:当最高纬度值α1达到90°及以上,即覆盖到极点的情况
卫星能够覆盖的纬度必然满足X的要求,仅需要通过多颗卫星视场拼接将单颗卫星的探测空白区域填补完整,即可满足全球一定地理纬度X的全覆盖要求;因此这种情况下需要的卫星总数量K为:K=[β/α]+2,其中β为单星未覆盖空域对应的地心角;
至此类型二所需要的所有步骤完成,后续步骤C、步骤D仅针对类型一。
所述步骤C中M颗卫星总共能够覆盖的空域对应的地心角δ为:δ=2α+β+Lα,其中,α为单星覆盖空域对应的地心角,β为单星未覆盖空域对应的地心角,L为需要填补的卫星数量;
但是,上述覆盖范围的两侧存在部分空域不满足地理纬度X的要求,其对应的地心角ψ需要去除,ψ的计算方法为:
式中,α1为单星临边探测可覆盖最高纬度值,R为地球半径,X为需要的覆盖的地理纬度。
因此,一个卫星组M颗卫星在一定地理纬度要求下有效覆盖范围对应的地心角为θ,其计算方法为:
θ=δ-2ψ。
式中,δ为一组M颗卫星总共能够覆盖的空域对应的地心角,ψ为不满足地理纬度要求的空域对应的地心角。
所述步骤D中对于整个地球360°地心角,需要的卫星组数量N为N=[360/θ];
由于N为向下取整,因此部署N组卫星后仍然留有一定的空白区域待填补;设这部分空白区域对应的地心角为Δ,那么:
若Δ≤α-ψ,则需要再增加1颗卫星,即卫星总数K为:K=N×M+1;
若Δ≥α+β-ψ,则需要再增加1组卫星,即卫星总数K为:K=(N+1)×M;
若α-ψ<Δ<α+β-ψ,则需要增加Q颗卫星和卫星总数K为:K=N×M+Q+1。
本发明具有以下有益效果:
利用本发明所述方法,能够实现满足特定地理纬度带上方特定高度空域全覆盖要求的临边探测卫星星座最小数量规模设计,解决星座设计中数量优化的问题。本发明涉及的方法通过简单数值计算即可实施,不需要借助复杂的仿真工具。
附图说明
图1为临边探测方式覆盖示意图。
图2为本发明中的卫星观测覆盖原理图(南北极剖面视角)。
图3为本发明中的卫星观测覆盖原理图(赤道剖面视角)。
图4为本发明中的两星视场紧密衔接示意图。
图5为本发明中的多颗卫星拼接覆盖空域示意图。
图6为本发明中的一定地理纬度覆盖对应的地心角。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
本发明中涉及的卫星探测方式是临边探测方式,卫星运行轨道倾角为0°。由于卫星对地覆盖能力与探测有效载荷的工作模式有关,本发明涉及的卫星探测载荷是指具备360°环扫探测工作模式的有效载荷,其对地探测示意图如图1所示。
具体设计方法如下:
设卫星运行高度为H,地球半径为R,需要观测的空域高度为h,要求覆盖的地理纬度带为X。各个相关参数如图2~图5所示。
步骤A:计算临边探测方式下单星覆盖空域对应的地心角和单星不能覆盖的空白空域对应的地心角。
单星临边探测可覆盖最高纬度值α1:
单星覆盖空域对应的地心角α:
单星未覆盖空域对应的地心角β:
β=2(α1-α)
步骤B:计算填补空白空域需要的卫星数量,构成一组无覆盖空白空域的卫星。
类型一:当最高纬度值α1小于90°,即覆盖不超过极点的情况
为填补单星未覆盖空域,需要多颗卫星视场紧密衔接填补未覆盖空域。由于视场是弧形,最高可覆盖地理纬度为视场交点。采用卫星视场紧密衔接部署方式后,所能够覆盖的地理纬度为gama:
若多星视场拼接后能够覆盖的地理纬度gama能够满足需要的覆盖的地理纬度X,即gama>X时,可证明可通过多星拼接满足要求。为了填补单星的未覆盖空域,需要填补的卫星数量L为:
L=[β/α]+1
因此,L颗卫星加上原本的1颗卫星构成无空白空域的卫星组,卫星组的卫星数量M为:
M=L+1=[β/α]+2
类型二:当最高纬度值α1达到90°及以上,即覆盖到极点的情况
卫星能够覆盖的纬度必然满足X的要求,仅需要通过多颗卫星视场拼接将单颗卫星的探测空白区域填补完整,即可满足全球一定地理纬度X的全覆盖要求。因此这种情况下需要的卫星总数量K为:
K=[β/α]+2
至此类型二所需要的所有步骤完成,后续步骤C、步骤D仅针对类型一。
步骤C:代入地理纬度要求计算这一组卫星覆盖空域对应的地心角。
那么,M颗卫星总共能够覆盖的空域对应的地心角δ为:
δ=2α+β+Lα
但是,上述覆盖范围的两侧存在部分空域不满足地理纬度X的要求,其对应的地心角ψ需要去除,ψ的计算方法为:
因此,一个卫星组M颗卫星在一定地理纬度要求下有效覆盖范围对应的地心角为θ
θ=δ-2ψ
步骤D:对照全球360°全覆盖的要求计算需要的总的卫星数量。
对于整个地球360°地心角,需要的卫星组数量N为N=[360/θ]
由于N为向下取整,因此部署N组卫星后仍然留有一定的空白区域待填补。设这部分空白区域对应的地心角为Δ,那么:
若Δ≤α-ψ,则需要再增加1颗卫星,即卫星总数K为:
K=N×M+1
若Δ≥α+β-ψ,则需要再增加1组卫星,即卫星总数K为:
K=(N+1)×M
若α-ψ<Δ<α+β-ψ,则需要增加Q颗卫星和卫星总数K为:
K=N×M+Q+1。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。
Claims (6)
1.一种卫星临边探测方式下星座最小规模设计方法,其特征在于:采用多颗卫星视场拼接的方式实现全球一定地理纬度空域的全覆盖,并且对拼接时是否满足地理纬度要求进行了有效性判断;星座内所有卫星的轨道高度一致,且均采用0°倾角轨道,通过代入不同轨道高度值迭代计算得到不同轨道高度下卫星星座最小规模。
2.如权利要求1所述的一种卫星临边探测方式下星座最小规模设计方法,其特征在于:包括如下步骤:
步骤A:计算临边探测方式下单星覆盖空域对应的地心角和单星不能覆盖的空白空域对应的地心角;
步骤B:计算填补空白空域需要的卫星数量,构成一组无覆盖空白空域的卫星;
步骤C:代入地理纬度要求计算这一组卫星覆盖空域对应的地心角;
步骤D:对照全球360°全覆盖的要求计算需要的总的卫星数量。
3.如权利要求1所述的一种卫星临边探测方式下星座最小规模设计方法,其特征在于,所述步骤A中通过以下公式计算单星覆盖空域对应的地心角α:
式中,ψ2为观测空域高度对应的地心角的余角,R为地球半径,h为需要观测的空域高度;
单星未覆盖空域对应的地心角β:
β=2(α1-α);
式中,α1为单星临边探测可覆盖最高纬度值。
4.如权利要求1所述的一种卫星临边探测方式下星座最小规模设计方法,其特征在于,所述步骤B分为两种类型:
类型一:当最高纬度值α1小于90°,即覆盖不超过极点的情况
采用卫星视场紧密衔接部署方式后,所能够覆盖的地理纬度为gama:
式中,R为地球半径,α1为单星临边探测可覆盖最高纬度值,α为单星覆盖空域对应的地心角;
若多星视场拼接后能够覆盖的地理纬度gama能够满足需要的覆盖的地理纬度X,即gama>X时,可证明可通过多星拼接满足要求;为了填补单星的未覆盖空域,需要填补的卫星数量L为:L=[β/α]+1,其中β为单星未覆盖空域对应的地心角;
因此,L颗卫星加上原本的1颗卫星构成无空白空域的卫星组,卫星组的卫星数量M为:M=L+1=[β/α]+2;
类型二:当最高纬度值α1达到90°及以上,即覆盖到极点的情况
卫星能够覆盖的纬度必然满足X的要求,仅需要通过多颗卫星视场拼接将单颗卫星的探测空白区域填补完整,即可满足全球一定地理纬度X的全覆盖要求;因此这种情况下需要的卫星总数量K为:K=[β/α]+2,其中β为单星未覆盖空域对应的地心角;
至此类型二所需要的所有步骤完成,后续步骤C、步骤D仅针对类型一。
5.如权利要求1所述的一种卫星临边探测方式下星座最小规模设计方法,其特征在于,所述步骤C中M颗卫星总共能够覆盖的空域对应的地心角δ为:δ=2α+β+Lα,其中,α为单星覆盖空域对应的地心角,β为单星未覆盖空域对应的地心角,L为需要填补的卫星数量;
但是,上述覆盖范围的两侧存在部分空域不满足地理纬度X的要求,其对应的地心角ψ需要去除,ψ的计算方法为:
式中,α1为单星临边探测可覆盖最高纬度值,R为地球半径,X为需要的覆盖的地理纬度;
因此,一个卫星组M颗卫星在一定地理纬度要求下有效覆盖范围对应的地心角为θ,其计算方法为:
θ=δ-2ψ;
式中,δ为一组M颗卫星总共能够覆盖的空域对应的地心角,ψ为不满足地理纬度要求的空域对应的地心角。
6.如权利要求1所述的一种卫星临边探测方式下星座最小规模设计方法,其特征在于,所述步骤D中对于整个地球360°地心角,需要的卫星组数量N为N=[360/θ];
由于N为向下取整,因此部署N组卫星后仍然留有一定的空白区域待填补;设这部分空白区域对应的地心角为Δ,那么:
若Δ≤α-ψ,则需要再增加1颗卫星,即卫星总数K为:K=N×M+1;
若Δ≥α+β-ψ,则需要再增加1组卫星,即卫星总数K为:K=(N+1)×M;
若α-ψ<Δ<α+β-ψ,则需要增加Q颗卫星和卫星总数K为:K=N×M+Q+1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910037300.XA CN109840360B (zh) | 2019-01-15 | 2019-01-15 | 一种卫星临边探测方式下星座最小规模设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910037300.XA CN109840360B (zh) | 2019-01-15 | 2019-01-15 | 一种卫星临边探测方式下星座最小规模设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109840360A true CN109840360A (zh) | 2019-06-04 |
CN109840360B CN109840360B (zh) | 2023-02-17 |
Family
ID=66883844
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910037300.XA Active CN109840360B (zh) | 2019-01-15 | 2019-01-15 | 一种卫星临边探测方式下星座最小规模设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109840360B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110954088A (zh) * | 2019-12-16 | 2020-04-03 | 中国科学院微小卫星创新研究院 | 一种具有高覆盖率的用于观测空间目标的方法 |
CN110954087A (zh) * | 2019-12-16 | 2020-04-03 | 中国科学院微小卫星创新研究院 | 一种具有高覆盖率的用于观测空间目标的系统 |
CN111800182A (zh) * | 2020-06-17 | 2020-10-20 | 北京理工大学 | 一种实现柔性覆盖全球通信星座的设计方法 |
CN113176072A (zh) * | 2021-03-30 | 2021-07-27 | 上海卫星工程研究所 | 适用于对地光学探测星座阳光干扰的通用分析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1633087A (zh) * | 2004-12-09 | 2005-06-29 | 中国科学院软件研究所 | 一种低轨道极星座卫星通信系统星座参数测量方法 |
WO2010096592A2 (en) * | 2009-02-19 | 2010-08-26 | Korb C Laurence | Methods for optimizing the performance, cost and constellation design of satellites for full and partial earth coverage |
CN105721040A (zh) * | 2016-01-25 | 2016-06-29 | 张雅声 | 一种区域覆盖推扫星座及其构建方法 |
WO2017023621A1 (en) * | 2015-07-31 | 2017-02-09 | Viasat, Inc. | Flexible capacity satellite constellation |
-
2019
- 2019-01-15 CN CN201910037300.XA patent/CN109840360B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1633087A (zh) * | 2004-12-09 | 2005-06-29 | 中国科学院软件研究所 | 一种低轨道极星座卫星通信系统星座参数测量方法 |
WO2010096592A2 (en) * | 2009-02-19 | 2010-08-26 | Korb C Laurence | Methods for optimizing the performance, cost and constellation design of satellites for full and partial earth coverage |
WO2017023621A1 (en) * | 2015-07-31 | 2017-02-09 | Viasat, Inc. | Flexible capacity satellite constellation |
CN105721040A (zh) * | 2016-01-25 | 2016-06-29 | 张雅声 | 一种区域覆盖推扫星座及其构建方法 |
Non-Patent Citations (2)
Title |
---|
WEILIAN SU ET AL.: ""Global communication coverage using Cubesats"", 《2017 IEEE 7TH ANNUAL COMPUTING AND COMMUNICATION》 * |
张景成等: ""基于IBEA的低轨区域覆盖星座优化设计"", 《计算机应用与软件》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110954088A (zh) * | 2019-12-16 | 2020-04-03 | 中国科学院微小卫星创新研究院 | 一种具有高覆盖率的用于观测空间目标的方法 |
CN110954087A (zh) * | 2019-12-16 | 2020-04-03 | 中国科学院微小卫星创新研究院 | 一种具有高覆盖率的用于观测空间目标的系统 |
CN110954087B (zh) * | 2019-12-16 | 2020-09-01 | 中国科学院微小卫星创新研究院 | 一种具有高覆盖率的用于观测空间目标的系统 |
CN112083445A (zh) * | 2019-12-16 | 2020-12-15 | 中国科学院微小卫星创新研究院 | 一种减少用于观测空间目标的观测卫星的数目的方法 |
CN112083445B (zh) * | 2019-12-16 | 2024-01-26 | 中国科学院微小卫星创新研究院 | 一种减少用于观测空间目标的观测卫星的数目的方法 |
CN111800182A (zh) * | 2020-06-17 | 2020-10-20 | 北京理工大学 | 一种实现柔性覆盖全球通信星座的设计方法 |
CN113176072A (zh) * | 2021-03-30 | 2021-07-27 | 上海卫星工程研究所 | 适用于对地光学探测星座阳光干扰的通用分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109840360B (zh) | 2023-02-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109840360A (zh) | 一种卫星临边探测方式下星座最小规模设计方法 | |
CN106249253A (zh) | 低轨通信和导航增强混合星座的优化设计方法 | |
US10732292B2 (en) | Method and system for dealing with antenna blockage in a low earth orbit constellation | |
CN106767730A (zh) | 用静态网格描述的卫星动态条带区域拆分方法 | |
CN106788671B (zh) | 一种可变构形的混轨卫星星座系统 | |
CN110727903B (zh) | 一种有限覆盖资源实现最大观测面积的卫星任务规划方法 | |
CN109639338B (zh) | 一种适用于通导遥一体化应用的全球覆盖星座的设计方法 | |
Cook et al. | Relative accuracy of rectifications using coordinates determined from maps and the global positioning system | |
CN112083445B (zh) | 一种减少用于观测空间目标的观测卫星的数目的方法 | |
CN109413662A (zh) | 一种低轨通信卫星星座与用户站连通规划方法 | |
CN113108776B (zh) | 一种具有高覆盖率的空间目标观测方法 | |
Sai et al. | SuperView-1-China’s first commercial remote sensing satellite constellation with a high resolution of 0.5 m | |
CN110717673B (zh) | 一种资源充足情形下最小观测成本的卫星任务规划方法 | |
CN112504278B (zh) | 一种机动测控站导航规划方法 | |
CN113758467A (zh) | 基于区域划分和局部网格嵌套的遥感卫星区域观测方法 | |
CN112130147A (zh) | 一种基于海陆目标位置信息的成像波位确定方法 | |
CN111800182B (zh) | 一种实现柔性覆盖全球通信星座的设计方法 | |
CN116208236B (zh) | 一种sar卫星星座任务规划方法 | |
CN111637881B (zh) | 一种基于地表模型图的星下点地表属性预报方法 | |
Petrie | An Alternative Source of Very High-resolution Imagery: The Resurs-DK1 Satellite | |
CN105629232A (zh) | 用于geosar相位定标的多地面定标站的布局方法 | |
Xue et al. | Research on a Decomposition Method for Remote Sensing Satellite Monitoring Tasks in Striped Regions | |
CN118153915A (zh) | 基于椭球模型的卫星最大覆盖范围计算方法、设备及介质 | |
Codou et al. | Operational scheduling of direct tasking innovative concept to improve reactivity on earth observation system | |
da Silva Curiel et al. | Real-time mosaic-rapid response, high resolution imaging from space |
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 |