一种冻结轨道的偏心率的卫星轨道确定方法技术

本发明专利技术公开了一种冻结轨道的偏心率的卫星轨道确定方法，该方法包括有初始轨道参数设置，轨道要素的平根数到瞬根数的转换，利用STK软件辅助计算出寿命期内的轨道要素，提取偏心率的极值处理和冻结轨道更新的步骤。本发明专利技术提出的方法解决了原有冻结轨道在确定偏心率过程中，在采用高精度计算方法计算时出现计算复杂且时间长，或者在采用低精度计算方法时出现计算过程简单、且精度低的缺陷。

【技术实现步骤摘要】
【专利摘要】本专利技术公开了，该方法包括有初始轨道参数设置，轨道要素的平根数到瞬根数的转换，利用STK软件辅助计算出寿命期内的轨道要素，提取偏心率的极值处理和冻结轨道更新的步骤。本专利技术提出的方法解决了原有冻结轨道在确定偏心率过程中，在采用高精度计算方法计算时出现计算复杂且时间长，或者在采用低精度计算方法时出现计算过程简单、且精度低的缺陷。【专利说明】
本专利技术涉及一种卫星轨道的确定方法，更特别地说，是指一种冻结轨道的偏心率的确定方法。
技术介绍
卫星在轨飞行期间，地球非球型摄动会带来轨道偏心率和近地点幅角的短周期项和长期项变化:短周期项变化由卫星所处纬度决定，无法通过适当的轨道设计予以消除；近地点幅角的长期项将造成轨道拱线的在赤道面内不断进动(以645km高度太阳同步轨道为例，进动周期约为100天)，将对卫星应用产生不利影响，希望通过卫星轨道设计予以消除。轨道偏心率的变化极其重要，它是制约各类(不同高度)空间飞行体轨道寿命的关键因素之一。由于地球重力场同时存在奇数阶和偶数阶带谐函数系数，即南北半球分布不对称，若将卫星的近地点幅角平根设置为90°以及偏心率平根设置为适当数值，可以消除偏心率和近地点幅角的长期项；该偏心率数值即为冻结轨道的偏心率。对于对地观测遥感和测绘卫星来说，冻结轨道的偏心率和近地点幅角可以保证卫星经过任一纬度时具有相同的对地高度，这对于相机成像是极其有利的。中国专利申请号CN200910089663.4，申请日2009年7月23日，专利技术名称“全天候覆盖卫星的重访轨道确定方法”。该专利申请通过考察升/降交点的分布顺序随地理纬度变化的规律，可针对双轨覆盖卫星进行重访轨道设计，弥补现有方法仅可针对单轨覆盖进行设计的不足。冻结轨道的偏心率的数值与轨道半长轴、倾角以及重力场模型带谐函数系数等有关，求解方法大都采用解析算法:只采用重力场模型带谐函数前两阶系数JjPJ3求解，算法简单但精度较差(杨维廉，基于Brouwer平根数的冻结轨道，中国空间科学技术，1998年10月，第5期，16-21);采用高阶谐函数求解冻结偏心率，可将计算精度提高，但需要展开勒让德多项式导致算法复杂(杨维廉，火星卫星的冻结轨道研究，航天器工程，2011年5月，第3期，20-24)。此外，上述算法均没有考虑重力场模型的田谐函数以及日月引力等因素。
技术实现思路
为了解决原有冻结轨道在确定偏心率过程中，在采用高精度计算方法计算时出现计算复杂且时间长，或者在采用低精度计算方法时出现计算过程简单、且精度低的缺陷，本专利技术提出了。本专利技术提出的，该方法包括有初始轨道要素参数设置步骤；轨道要素平根数到瞬根数的转换步骤；采用STK软件计算寿命期内的轨道要素的步骤；偏心率极值处理步骤和卫星冻结轨道更新截止步骤。本专利技术提出的卫星冻结轨道的偏心率的确定方法优点在于:①本专利技术的卫星冻结轨道确定方法通过考虑了重力场模型的田谐函数以及日月引力等因素影响，可解决原有确定冻结偏心率方法精度低的不足。②本专利技术的卫星冻结轨道确定方法辅助利用STK软件对不同时刻下的偏心率进行获取，缩短了轨道计算时间，并使得寿命期内的偏心率差值小于阈值，提高了偏心率的稳定性。③本专利技术的卫星冻结轨道确定方法在获得冻结偏心率的过程中，避免了展开勒让德多项式等复杂操作。【专利附图】【附图说明】图1改进前寿命期内偏心率的平根数曲线图。图2改进后寿命期内偏心率的平根数曲线图。【具体实施方式】下面将结合附图和实施例对本专利技术做进一步的详细说明。在1995年12月第I版《航天器飞行动力学原理》，肖业伦编著，第44页的图2_13轨道要素定义中，一般卫星轨道包括有六个要素，其中:轨道偏心率记为e，轨道偏心率e的平根数记为Ae，轨道偏心率e的瞬根数记为Be ；轨道半长轴记为a，轨道`半长轴a的平根数记为Aa，轨道半长轴a的瞬根数记为Ba ；近地点幅角记为ω，近地点幅角ω的平根数记为A ω，近地点幅角ω的瞬根数记为Βω ;纬度幅角记为U，纬度幅角u的平根数记为Au，纬度幅角记为u的瞬根数记为Bu ；轨道倾角记为i，轨道倾角i的平根数记为Ai，轨道倾角i的瞬根数记为Bi ；轨道升交点赤经记为Ω，轨道升交点赤经Ω的平根数记为A Ω，轨道升交点赤经Ω的瞬根数记为ΒΩ。为了方便表述说明，卫星轨道要素采用集合形式表达为D= {e,a, ω,π,?, Ω}，卫星轨道要素的平根数采用集合形式表达为AD = {Ae，Aa，A?，AU，Ai，AQ}，卫星轨道要素的瞬根数采用集合形式表达为BD = {Be, Ba,Βω, Bu, Bi，B Ω }。在本专利技术中，将平均轨道根数(或者平均轨道要素)简称为平根数，密切轨道根数(或者密切轨道要素)简称为瞬根数。平根数转换成瞬根数采用S G P 4 (简化普适摄动)算法，请参照:刘光明，文援兰，廖瑛，“基于无奇异变换的双行轨道根数生成算法”，系统工程与电子技术，2011, 33 (5):1104-1123, Liu Guangming, Wen Yuanlan, Liao Ying.Fitting Algorithm of TLE Parameters Based on Non-singular Transformation .Systems Engineering and Electronics。在本专利技术中，卫星寿命期记为T，卫星寿命期T内的初始时刻记为h，第一时刻记为ti，第二时刻记为t2，卫星寿命期T内的任意一时刻记为ti;ti时刻的前一时刻记为ti_1;ti时刻的后一时刻记为ti+1，卫星寿命期T内的最后时刻记为tn，η表示偏心率采集次数。在本专利技术中，冻结轨道(Frozen Orbit)的轨道要素采用矩阵表达为:【权利要求】1.，其特征在于该方法包括有下列步骤: 步骤1:初始参数设置 在初始时刻h下，设置: 步骤1 一 1:轨道半长轴的平根数Aattl、轨道倾角的平根数Aittl、纬度幅角的平根数AuV轨道升交点赤经的平根数A Ω t0 ; 步骤1 一 2:近地点幅角的平根数A ω \，且A ω \等于90度； 步骤1 — 3:轨道偏心率的平根数Aettl,且 2.根据权利要求1所述的冻结轨道的偏心率的卫星轨道确定方法，其特征在于:所述步骤3的STK的计算处理包括有下列步骤: 步骤3 — 1:将步骤二得到的t0时刻的瞬根数BDt0 = (Bet0, Bat0, B ω ; But0, Bit0, B Ω tJ加载至STK软件中进行轨道计算，获得STK软件计算输出的h时刻的轨道偏心率Cettl、轨道半长轴Cattl、近地点幅角Co V纬度幅角Cuttl、轨道倾角Cittl、轨道升交点赤经CQ V为了方便表述说明，h时刻的STK软件计算输出的卫星轨道要素采用集合形式表达为⑶\ = ICetο, Cat0lCwt0, Cut0, Cit0, CQ V}。 步骤3 — 2:将步骤二得到的t0时刻的瞬根数BDt0 = (Bet0, Bat0, B ω V But0, Bit0, B Ω tJ加载至STK软件中进行轨道计算，获得STK软件计算输出的^时刻的轨道偏心率Cet1、轨道半长轴Cat1、近地点幅角Ccot1、纬度幅角Cut1、轨道倾角Cit1、本文档来自技高网...

【技术保护点】
一种冻结轨道的偏心率的卫星轨道确定方法，其特征在于该方法包括有下列步骤：步骤1：初始参数设置在初始时刻t0下，设置：步骤1－1：轨道半长轴的平根数Aat0、轨道倾角的平根数Ait0、纬度幅角的平根数Aut0、轨道升交点赤经的平根数AΩt0；步骤1－2：近地点幅角的平根数Aωt0，且Aωt0等于90度；步骤1－3：轨道偏心率的平根数Aet0，且Aet0=J3RsinAi2J2Aa;]]>其中，R为地球赤道半径，J2为重力场模型带谐函数的第二阶系数，J3为重力场模型带谐函数的第三阶系数；步骤1－4：应用初始时刻t0的设置更新冻结轨道FT=efafωfufifΩf,]]>得到第一更新冻结轨道Aet0Aat0Aωt0Aut0Ait0AΩt0,]]>且Aet0=J3RsinAi2J2Aa;]]>Aωt0＝90°。执行步骤2；步骤2：平根数到瞬根数的转换将轨道偏心率的平根数Aet0转换为轨道偏心率的瞬根数Bet0；将轨道半长轴的平根数Aat0转换为轨道半长轴的瞬根数Bat0；将近地点幅角的平根数Aωt0转换为近地点幅角的瞬根数Bωt0；将纬度幅角的平根数Aut0转换为纬度幅角瞬根数But0；将轨道倾角的平根数Ait0转换为轨道倾角的瞬根数Bit0；将轨道升交点赤经的平根数AΩt0转换为轨道升交点赤经的瞬根数BΩt0；t0时刻的卫星轨道要素的平根数ADt0＝{Aet0,Aat0,Aωt0,Aut0,Ait0,AΩt0}转换成瞬根数记为BDt0＝{Bet0,Bat0,Bωt0,But0,Bit0,BΩt0}；执行步骤3；步骤3：采用STK计算寿命期内的轨道要素对BDt0＝{Bet0,Bat0,Bωt0,But0,Bit0,BΩt0}进行STK的处理，输出轨道偏心率CeT＝{Cet0,Cet1,Cet2,…,Ceti,…,Cetn}；步骤4：偏心率的极值处理从STK输出轨道偏心率CeT＝{Cet0,Cet1,Cet2,…,Ceti,…,Cetn}中选出轨道偏心率的最大值记为最小值记为并计算STK输出偏心率的差值记为及STK输出偏心率的均值记为执行步骤5－1；步骤5：冻结轨道更新截止条件步骤5－1：设置偏心率阈值e阈值，且e阈值＝6×10‑6；步骤5－2：若Δe≤e阈值，不对第一更新冻结轨道Aet0Aat0Aωt0Aut0Ait0AΩt0]]>进行更新，从而结束卫星冻结轨道的偏心率的确定；步骤5－3：若Δe＞e阈值，则同时将步骤1－2所述的近地点幅角平根数Aωt0赋值给冻结轨道的近地点幅角ωf，所述的STK输出偏心率均值赋值给冻结轨道的轨道偏心率ef，即对第一更新冻结轨道Aet0Aat0Aωt0Aut0Ait0AΩt0]]>进行更新，得到第二更新冻结轨道且Aet0=e‾=Σi=0nCetin+1,]]>Aωt0＝90°；执行步骤5－4；步骤5－4：重复步骤2、步骤3和步骤4直至STK输出偏心率的差值小于等于偏心率阈值e阈值时，结束卫星冻结轨道的偏心率的确定。...

【技术特征摘要】

【专利技术属性】
技术研发人员：徐明，魏延，汪作鹏，徐世杰，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11