一种基于地月锥时空剖分的频谱信息表征方法及系统

1.本发明属于地月空间科学探索及资源开发利用领域，尤其涉及一种基于地月锥时空剖分的频谱信息表征方法及系统。


背景技术：

2.地月空间已成为国际上极为活跃的科技探索前沿热点，许多国家陆续开展探月任务，在国际频率划分规定和“先登先用”的使用规则约束下，全球地月空间科学探索和开发利用的可用频率资源愈发紧张，供需不平衡加剧。由于地月空间飞行器的轨道特性复杂，搭载的无线电类型载荷的电磁波辐射范围和辐射路径存在大时空变化，为确保地月空间科学任务顺利实施，迫切需要地球至月球大时空场景下的表征架构和更为精细的计算方法，支撑技术和产业应用进一步发展。
3.目前，针对地球高低轨空间的表征与计算架构主要采用基于经纬度的栅格化的方法，然而，由于地球、月球的自转及公转，基于地球经纬度的剖分架构在地月锥形空间中具有时变性，难以延展应用。面向地月空间大时空场景下的时、空、频、能等多维频谱信息，缺少统一尺度下的计算与表征方法，制约了其对地月空间科学卫星通信能量、干扰等频谱信息的定量化计算与时空立体表达。


技术实现要素：

4.本发明的目的在于克服现有技术缺陷，提出了一种基于地月锥时空剖分的频谱信息表征方法及系统。为地月空间电磁频谱态势的分析、计算、应用提供统一的表征架构，解决地月空间电磁频谱大数据同化和表征的问题。
5.为了实现上述目的，本发明提出了一种基于地月锥时空剖分的频谱信息表征方法，所述方法包括：
6.步骤1)采用地月会合坐标系为时空基准，划分地月空间锥形；
7.步骤2)计算地月空间锥形剖分架构；
8.步骤3)基于每个剖分点进行频谱信息的计算与表征。
9.作为上述方法的一种改进，所述步骤1)具体包括：
10.采用地月会合坐标系为时空基准，标定地球高轨空间至月球l2点轨道簇最远轨迹，且覆盖l1点轨道簇最远轨迹范围内为地月空间锥形；其中，地月空间锥形的左侧边界点l(x
l
,0,0)，右侧边界点r(xr,0,0)，l1点轨道簇上的点为(x
l1
,0,z
l1
)。
11.作为上述方法的一种改进，所述步骤2)具体包括：
12.步骤2-1)令沿x轴的剖分截面间的步长为δx，令x轴坐标为x的截面内的剖分角度步长为δθ
x
、剖分半径步长为δρ
x
，l点坐标为(x
l
,0,0)，r点坐标为(xr,0,0)，地心坐标为(xe,0,0)，月心坐标为(xm,0,0)，令x
pre
＝x
l-δx,re为地球半径，rm为月球半径；
13.步骤2-2)计算当前剖面的横坐标x＝x
pre
+δx，判断x＜＝xr是否成立，判断为是，则执行步骤2-3)，判断为否，执行步骤2-10)；
14.步骤2-3)计算横坐标为x时的剖面半径r
x
，令θ
pre
＝2π+δθ
x
+δbias，ρ
pre
＝r
x
+δρ
x
；δbias为相邻截面间剖分角度的偏移值；
15.步骤2-4)计算ρ＝ρ
pre-δρ
x
，判断ρ≥0是否成立，判断为是，则执行步骤2-5)，判断为否，则令x
pre
＝x，执行步骤2-2)；
16.步骤2-5)计算θ＝θ
pre-δθ
x
，判断θ≥0是否成立，判断为是，则执行步骤2-6)，判断为否，令ρ
pre
＝ρ，执行步骤2-4)；
17.步骤2-6)计算剖分点(ρ,θ)的坐标值y＝ρcosθ，z＝ρsinθ，若需要剔除地月星体内部点，则执行步骤2-7)，否则，执行步骤2-8)；
18.步骤2-7)判断|x-xe|≤re是否成立，判断成立，则计算点(x,y,z)与地心的距离d
pe
，若满足d
pe
＜re，则执行步骤2-9)；判断|x-xm|≤rm是否成立，若成立，则计算点(x,y,z)与月心的距离d
pm
，若满足d
pm
＜rm，则执行步骤2-9)；
19.步骤2-8)记录当前剖分点的三维坐标(x,y,z)；
20.步骤2-9)令θ
pre
＝θ，执行步骤2-5)；
21.步骤2-10)记录所有剖分点坐标，结束。
22.作为上述方法的一种改进，所述步骤2-3)中计算横坐标为x时的剖面半径r
x
，具体包括：
23.根据下式计算xs和zs，其中，xs和zs分别表示椎体与地球高轨轨道相切的切点在x轴的坐标和在z轴的坐标：
24.xs＝-rgcos(γ)+xe，zs＝rgsin(γ)；
25.其中，γ＝π-α-β；
26.当x≤xs时，
27.当x＞xs时，r
x
＝(x-xs)cot(γ)+zs；
28.其中，rg为地球静止轨道半径。
29.作为上述方法的一种改进，所述步骤3)包括：
30.根据计算得到的剖分点坐标，计算在某时刻，某卫星的通信链路在划分的地月空间锥形中的功率通量密度分布。
31.作为上述方法的一种改进，所述步骤3)包括：
32.根据计算得到的剖分点坐标，计算在某时刻，某卫星的通信链路在划分的地月空间锥形中的自由空间链路损耗分布。
33.另一方面，本发明提出了一种基于地月锥时空剖分的频谱信息表征系统，所述系统包括：
34.地月空间锥形划分模块，用于采用地月会合坐标系为时空基准，划分地月空间锥形；
35.架构计算模块，用于计算地月空间锥形剖分架构；和
36.频谱计算表征模块，用于基于每个剖分点进行频谱信息的计算与表征。
37.与现有技术相比，本发明的优势在于：
38.本发明提出一种面向地月空间电磁频谱态势表征需求的锥形即地月锥时空剖分
方法，支撑对任意剖分点进行频谱相关信息的计算与表征，解决地月空间多尺度空间电磁频谱数据缺乏统一组织架构的问题。
附图说明
39.图1是地月会合坐标系下地月锥剖分示意图；
40.图2是沿-x方向的地月锥界面剖分示意图；
41.图3是本发明的基于地月锥时空剖分的频谱信息表征方法流程图；
42.图4是计算剖面半径的几何方法示意图；
43.图5是基于orion卫星网络资料的b2aa波束重构的dro卫星参考天线方向图；
44.图6是dro轨道卫星基于地月锥时空剖分方法的功率通量密度表征示意，其中图6(a)是基于地月锥时空剖分的pfd值表征，图6(b)是补全空值的地月锥剖分点；
45.图7是基于capstone卫星网络资料的b4波束重构的nrho卫星参考天线方向图；
46.图8是nrho轨道卫星基于本发明剖分方法的功率通量密度表征示意，其中图8(a)是基于地月锥时空剖分的pfd值表征，图8(b)是补全空值的地月锥剖分点。
具体实施方式
47.一、地月空间锥形范围定义
48.采用地月会合坐标系为时空基准，标定地球高轨空间至月球l2点轨道簇最远轨迹，且覆盖l1点轨道簇最远轨迹范围内为地月空间锥形，简称地月锥，如图1所示。
49.其中，边界节点及地心、月心坐标变量及含义说明如表1所示：
50.表1地月空间锥形相关点含义
[0051][0052]
二、地月空间锥形时空剖分方法
[0053]
地月锥在时间维度的剖分是指以指定时间步长进行切分，在空间维度的剖分，首先是以地月会合坐标系x轴上的固定分割步长进行剖分，如图1所示，其次对剖面上以极坐标下的(ρ,θ)进行剖分，如图2所示。
[0054]
基于地月锥的时空剖分，可以对任意剖分点进行频谱相关信息等的计算与表征。地月空间锥形时空剖分架构计算方法，具体步骤如下：如图3所示，
[0055]
1)令沿x轴的剖分截面间的步长为δx，令x轴坐标为x的截面内的剖分角度步长为
δθ
x
、剖分半径步长为δρ
x
，l点坐标为(x
l
,0,0)，r点坐标为(xr,0,0)，地心坐标为(xe,0,0)，月心坐标为(xm,0,0)，令x
pre
＝x
l-δx,re为地球半径，rm为月球半径；
[0056]
2)计算当前剖面的横坐标x＝x
pre
+δx，判断x＜＝xr是否成立，若是，则执行步骤3)，否则，执行步骤10)；
[0057]
3)计算横坐标为x时的剖面半径r
x
，令θ
pre
＝2π+δθ
x
+δbias，ρ
pre
＝r
x
+δρ
x
；
[0058]
4)计算ρ＝ρ
pre-δρ
x
，判断ρ≥0是否成立，若是，则执行步骤5)，否则，则令x
pre
＝x，执行步骤2)；
[0059]
5)计算θ＝θ
pre-δθ
x
，判断θ≥0是否成立，若是，则执行步骤6)，否则，令ρ
pre
＝ρ，执行步骤4)；
[0060]
6)计算剖分点(ρ,θ)的坐标值y＝ρcosθ，z＝ρsinθ，若需要剔除地月星体内部点，则执行步骤7)，否则，执行步骤8)；
[0061]
7)判断|x-xe|≤re是否成立，若成立，则计算点(x,y,z)与地心的距离d
pe
，若满足d
pe
＜re，则执行步骤9)；判断|x-xm|≤rm是否成立，若成立，则计算点(x,y,z)与月心的距离d
pm
，若满足d
pm
＜rm，则执行步骤9)；
[0062]
8)记录当前剖分点的三维坐标(x,y,z)；
[0063]
9)令θ
pre
＝θ，执行步骤5)；
[0064]
10)记录所有剖分点坐标，结束。
[0065]
其中，步骤3)中计算剖面横坐标为x的剖面半径r
x
方法步骤如下：
[0066]
1)当x≤xs时，
[0067]
2)当x＞xs时，r
x
＝(x-xs)cot(γ)+zs；
[0068]
其中，rg为地球静止轨道半径；
[0069]
xs＝-rgcos(γ)+xe，zs＝rgsin(γ)；
[0070]
γ＝π-α-β；
[0071][0072]
如图4所示是计算剖面半径的几何方法示意图。
[0073]
三、基于每个剖分点进行频谱信息的计算与表征
[0074]
下面结合附图和实施例对本发明的技术方案进行详细的说明。
[0075]
实施例1
[0076]
本发明的实施例1提出了一种基于地月锥时空剖分的频谱信息表征方法，包括：
[0077]
步骤1)采用地月会合坐标系为时空基准，划分地月空间锥形；
[0078]
步骤2)计算地月空间锥形剖分架构；
[0079]
步骤3)基于每个剖分点进行频谱信息的计算与表征。
[0080]
基于所提方法，以表征某dro轨道卫星在地月锥中的s频段的功率分布为例，步骤3)即计算其某一时刻下，通信链路在地月锥中的功率通量密度(pfd)分布。其中，地月锥边界参数选取如表2所示，坐标值基于地月会合坐标系计算。
[0081]
表2地月锥边界参数取值
[0082]
点坐标变量坐标值
l(x
l
,0,0)(-46915.7,0,0)r(xr,0,0)(461586.4,0,0)l
1n
(x
l1
,0,z
l1
)(363204.1,0,113899.2)月心坐标(xm,0,0)(386127.9,0,0)地心坐标(xe,0,0)(-4749.4,0,0)
[0083]
剖分架构计算采用参数如表3所示：
[0084]
表3实施例1中剖分点计算采用的参数
[0085][0086][0087]
仿真中采用的卫星系统为dro轨道，平均轨道高度65000公里，卫星天线指向堪培拉站，链路空口参数如表4所示，卫星天线方向图基于itu登记的orion卫星网络资料的b2aa波束，并针对波束旁瓣进行拟合重构，其参考天线方向图及重构模型如图5所示。
[0088]
表4dro卫星通信链路空口参数
[0089][0090]
基于所提方法，表征剖分点pfd值如图6(a)所示，由于计算中考虑了地球、月球遮挡因素，且卫星天线指向范围为[-90
°
,90
°
]，因此，存在部分剖分点值为空，将空值的剖分点补全如图6(b)所示。
[0091]
此外，步骤3)也可以是根据计算得到的剖分点坐标，计算在某时刻，dro卫星的通信链路在划分的地月空间锥形中的自由空间链路损耗分布。具体自由空间链路损耗分布属于本领域技术人员已知的方法，不再赘述。
[0092]
需要说明的是，步骤3)还可以是其他根据计算得到的剖分点坐标，进行相关能量、干扰等频谱信息的定量化计算与时空立体表达，均不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。
[0093]
实施例2
[0094]
本发明的实施例1提出了一种基于地月锥时空剖分的频谱信息表征方法，包括：
[0095]
步骤1)采用地月会合坐标系为时空基准，划分地月空间锥形；
[0096]
步骤2)计算地月空间锥形剖分架构；
[0097]
步骤3)基于每个剖分点进行频谱信息的计算与表征。
[0098]
基于所提方法，以表征l2点nrho轨道卫星在地月锥中的功率分布为例，步骤3)即计算其某一时刻下，通信链路在地月锥中的功率通量密度(pfd)分布。其中，地月锥边界参数选取参考表2，剖分架构计算采用参数如表5所示：
[0099]
表5实施例2中剖分点计算采用的参数
[0100][0101]
仿真中采用的卫星系统为l2点的nrho轨道，远月点约73000km，近月点约1900km，卫星天线指向堪培拉站，链路空口参数参考itu登记的capstone卫星网络资料参数，如表6所示，卫星天线方向图参考capstone卫星网络资料的b4波束，并针对波束旁瓣进行拟合重构，其参考天线方向图及重构模型如图7所示
[0102]
表6l2点nrho卫星通信链路空口参数
[0103][0104]
基于所提方法，仿真计算剖分点pfd值如图8(a)所示。由于计算中考虑了地球、月球遮挡因素，存在部分剖分点值为空，将空值的剖分点补全如图8(b)所示。
[0105]
此外，步骤3)也可以是根据计算得到的剖分点坐标，计算在某时刻，nrho卫星的通信链路在划分的地月空间锥形中的自由空间链路损耗分布。具体自由空间链路损耗分布属于本领域技术人员已知的方法，不再赘述。
[0106]
需要说明的是，步骤3)还可以是其他根据计算得到的剖分点坐标，进行相关能量、干扰等频谱信息的定量化计算与时空立体表达，均不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。
[0107]
实施例3
[0108]
本发明的实施例3提出了一种基于地月锥时空剖分的频谱信息表征系统，基于实施例1或实施例2的方法实现，该系统包括：
[0109]
地月空间锥形划分模块，用于采用地月会合坐标系为时空基准，划分地月空间锥形；
[0110]
架构计算模块，用于计算地月空间锥形剖分架构；和
[0111]
频谱计算表征模块，用于基于每个剖分点进行频谱信息的计算与表征。
[0112]
最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。

技术特征：
1.一种基于地月锥时空剖分的频谱信息表征方法，所述方法包括：步骤1)采用地月会合坐标系为时空基准，划分地月空间锥形；步骤2)计算地月空间锥形剖分架构；步骤3)基于每个剖分点进行频谱信息的计算与表征。2.根据权利要求1所述的基于地月锥时空剖分的频谱信息表征方法，其特征在于，所述步骤1)具体包括：采用地月会合坐标系为时空基准，标定地球高轨空间至月球l2点轨道簇最远轨迹，且覆盖l1点轨道簇最远轨迹范围内为地月空间锥形；其中，地月空间锥形的左侧边界点l(x
l
,0,0)，右侧边界点r(x
r
,0,0)，l1点轨道簇上的点为(x
l1
,0,z
l1
)。3.根据权利要求2所述的基于地月锥时空剖分的频谱信息表征方法，其特征在于，所述步骤2)具体包括：步骤2-1)令沿x轴的剖分截面间的步长为δx，令x轴坐标为x的截面内的剖分角度步长为δθ
x
、剖分半径步长为δρ
x
，l点坐标为(x
l
,0,0)，r点坐标为(x
r
,0,0)，地心坐标为(x
e
,0,0)，月心坐标为(x
m
,0,0)，令x
pre
＝x
l-δx,r
e
为地球半径，r
m
为月球半径；步骤2-2)计算当前剖面的横坐标x＝x
pre
+δx，判断x＜＝x
r
是否成立，判断为是，则执行步骤2-3)，判断为否，执行步骤2-10)；步骤2-3)计算横坐标为x时的剖面半径r
x
，令θ
pre
＝2π+δθ
x
+δbias，ρ
pre
＝r
x
+δρ
x
；δbias为相邻截面间剖分角度的偏移值；步骤2-4)计算ρ＝ρ
pre-δρ
x
，判断ρ≥0是否成立，判断为是，则执行步骤2-5)，判断为否，则令x
pre
＝x，执行步骤2-2)；步骤2-5)计算θ＝θ
pre-δθ
x
，判断θ≥0是否成立，判断为是，则执行步骤2-6)，判断为否，令ρ
pre
＝ρ，执行步骤2-4)；步骤2-6)计算剖分点(ρ,θ)的坐标值y＝ρcosθ，z＝ρsinθ，若需要剔除地月星体内部点，则执行步骤2-7)，否则，执行步骤2-8)；步骤2-7)判断|x-x
e
|≤r
e
是否成立，判断成立，则计算点(x,y,z)与地心的距离d
pe
，若满足d
pe
＜r
e
，则执行步骤2-9)；判断|x-x
m
|≤r
m
是否成立，若成立，则计算点(x,y,z)与月心的距离d
pm
，若满足d
pm
＜r
m
，则执行步骤2-9)；步骤2-8)记录当前剖分点的三维坐标(x,y,z)；步骤2-9)令θ
pre
＝θ，执行步骤2-5)；步骤2-10)记录所有剖分点坐标，结束。4.根据权利要求3所述的基于地月锥时空剖分的频谱信息表征方法，其特征在于，所述步骤2-3)中计算横坐标为x时的剖面半径r
x
，具体包括：根据下式计算x
s
和z
s
，其中，x
s
和z
s
分别表示锥体与地球高轨轨道相切的切点在x轴的坐标和在z轴的坐标：x
s
＝-r
g
cos(γ)+x
e
，z
s
＝r
g
sin(γ)；其中，γ＝π-α-β；当x≤x
s
时，
当x＞x
s
时，r
x
＝(x-x
s
)cot(γ)+z
s
；其中，r
g
为地球静止轨道半径。5.根据权利要求4所述的基于地月锥时空剖分的频谱信息表征方法，其特征在于，所述步骤3)包括：根据计算得到的剖分点坐标，计算在某时刻，某卫星的通信链路在划分的地月空间锥形中的功率通量密度分布。6.根据权利要求4所述的基于地月锥时空剖分的频谱信息表征方法，其特征在于，所述步骤3)包括：根据计算得到的剖分点坐标，计算在某时刻，某卫星的通信链路在划分的地月空间锥形中的自由空间链路损耗分布。7.一种基于地月锥时空剖分的频谱信息表征系统，其特征在于，所述系统包括：地月空间锥形划分模块，用于采用地月会合坐标系为时空基准，划分地月空间锥形；架构计算模块，用于计算地月空间锥形剖分架构；和频谱计算表征模块，用于基于每个剖分点进行频谱信息的计算与表征。

技术总结
本发明属于地月空间科学探索及资源开发利用领域，尤其涉及一种基于地月锥时空剖分的频谱信息表征方法及系统，所述方法包括：采用地月会合坐标系为时空基准，划分地月空间锥形；计算地月锥剖分架构；基于每个剖分点进行频谱信息的计算与表征。本发明提出一种面向地月空间电磁频谱态势表征需求的锥形即地月锥时空剖分方法，支撑对任意剖分点进行频谱相关信息的计算与表征，解决地月空间多尺度空间电磁频谱数据缺乏统一组织架构的问题。磁频谱数据缺乏统一组织架构的问题。磁频谱数据缺乏统一组织架构的问题。


技术研发人员：王静 姚秀娟 李震 高翔 智佳 李雪
受保护的技术使用者：中国科学院国家空间科学中心
技术研发日：2023.05.11
技术公布日：2023/8/23