CN102745343B - 同步卫星早期测控段干扰弧段快速预报方法 - Google Patents
同步卫星早期测控段干扰弧段快速预报方法 Download PDFInfo
- Publication number
- CN102745343B CN102745343B CN201210239598.0A CN201210239598A CN102745343B CN 102745343 B CN102745343 B CN 102745343B CN 201210239598 A CN201210239598 A CN 201210239598A CN 102745343 B CN102745343 B CN 102745343B
- Authority
- CN
- China
- Prior art keywords
- satellite
- point
- synchronous
- synchronous satellite
- fixed
- 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.)
- Expired - Fee Related
Links
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种同步卫星早期测控段干扰弧段快速预报方法,属于航天测量与控制领域,涉及一种基于测站波束范围和同步星点频的卫星空间几何干扰计算方法。该方法根据待定点同步卫星的精密轨道计算预报弧段内卫星的位置
Description
技术领域
本发明专利属于航天测量与控制领域,涉及一种基于测站波束范围和同步星点频的卫星空间几何干扰计算方法。
背景技术
同步卫星在定点之前的准同步段,由于其轨道高度和周期都很接近同步卫星,加之所有的同步卫星都是定点在赤道上空,这样一来可能导致卫星相对某一测控站的向量和已定点的某卫星相对这个测控站的向量之间的夹角在某一时间段内非常小,这时测控站发令就可能会对已定点的同步卫星的信号造成不同程度的干扰,给该星的用户带来许多不便。为了成功地对卫星进行控制并不对其所经区域的同步卫星造成干扰,就应避免在可能产生干扰的时间段内进行控制。如何准确的计算这些可能产生干扰的弧段,以前的方法是根据卫星当前的位置和附近同步卫星的定点经度来计算,而同步卫星在定点之后,受一些摄动因素的影响会偏离原来的位置,根据国际卫星组织的规定同步卫星的定点范围为±0.05度,所以同步卫星在定点后必须定期对其进行控制,以保证卫星的正常运行,这样一来用定点经度计算的结果和实际的偏差较大,可能会导致一些和实际正好相反的结果,造成不必要的麻烦,鉴于上述原因本项目给出一种根据TLE来进行预报干扰弧段的计算方法。
发明内容
本发明的目的是:为了解决现有技术中同步卫星早期测控段干扰弧段快速预报方法存在的偏差大的问题,本发明给出一种根据TLE来进行预报干扰弧段的方法。
本发明的技术方案是:同步卫星早期测控段干扰弧段快速预报方法,包括如下步骤:
步骤一:根据待定点同步卫星的精密轨道计算预报弧段内卫星的位置
步骤二:i=1;
步骤三:根据第i个已定点卫星的TLE中卫星一天运行的圈数筛选出同步卫星,如果第i个卫星是同步卫星,则进入步骤四,否则,认为待定点同步卫星对第i个卫星不会造成干扰,进入步骤七;
步骤四:计算第i个已定点同步卫星的定点经度ai,如果|ai-bi|≤5°,进入步骤五;否则,认为待定点同步卫星对第i个已定点同步卫星不会造成干扰,进入步骤七;其中,bi为待定点同步卫星的当前经度;
步骤五:根据第i个已定点同步卫星的TLE计算该卫星在预报弧段内的位置根据待定点同步卫星位置计算t时刻对第i个已定点同步卫星在某测站的仰角,如果仰角大于零度,即跟踪可见,进入步骤六,否则,认为待定点同步卫星对第i个已定点同步卫星不会造成干扰,进入步骤七;
步骤六:根据t时刻第i个已定点同步卫星和待定点同步卫星的位置和计算两个卫星相对测站的矢量夹角:
上式中为测站的位置;
如果θi(t)≤θ0,给定值θ0满足0.01°≤θ0≤0.1°,则认为该时刻待定点同步卫星会对第i个已定点同步卫星造成干扰;
步骤七:如果i<N,则:i=i+1,进入步骤三;否则:同步卫星早期测控段干扰弧段快速预报完成。
本发明的有益效果是:快速准确地进行干扰弧段预报,确定可能受干扰的卫星和时间段,对测控事件的安排进行适当调整,在顺利完成实时测控任务的同时,也避免此类纠纷的发生,确保了同步卫星定点工作的顺利完成。
具体实施方式
本实施例中已定点卫星个数为60个,需要计算干扰的测站为Sta_A,本实施例中同步卫星早期测控段干扰弧段快速预报方法包括如下步骤:
步骤一:根据待定点同步卫星的精密轨道计算预报弧段内卫星的位置
步骤二:i=1;
步骤三:根据第i个已定点卫星的TLE中卫星一天运行的圈数筛选出同步卫星,如果第i个卫星是同步卫星,则进入步骤四,否则,认为待定点同步卫星对第i个卫星不会造成干扰,进入步骤七;
步骤四:计算第i个已定点同步卫星的定点经度ai,如果|ai-bi|≤5°,进入步骤五,否则,认为待定点同步卫星对第i个已定点同步卫星不会造成干扰,进入步骤七;其中,bi为待定点同步卫星的当前经度;
步骤五:根据第i个已定点同步卫星的TLE计算该卫星在预报弧段内的位置根据待定点同步卫星位置计算t时刻对第i个已定点同步卫星在某测站的仰角,如果仰角大于零度,即跟踪可见,进入步骤六,否则,认为待定点同步卫星对第i个已定点同步卫星不会造成干扰,进入步骤七;
步骤六:根据t时刻第i个已定点同步卫星和待定点同步卫星的位置和计算两个卫星相对测站的矢量夹角:
上式中为Sta_A测站的位置;
如果θi(t)≤θ0,本实施例中给定值θ0=0.05°,则认为该时刻待定点同步卫星会对第i个已定点同步卫星造成干扰;
步骤七:如果i<60,则:i=i+1,进入步骤三;否则:同步卫星早期测控段干扰弧段快速预报完成。
经对下载的所有TLE进行计算,待定点卫星会对第12、35个已定点同步卫星造成干扰,干扰弧段见下表:
本实施例中干扰弧段预报结果
Claims (1)
1.同步卫星早期测控段干扰弧段快速预报方法,包括如下步骤:
步骤一:根据待定点同步卫星的精密轨道计算预报弧段内卫星的位置
步骤二:i=1;
步骤三:根据第i个已定点卫星的TLE中卫星一天运行的圈数筛选出同步卫星,如果第i个卫星是同步卫星,则进入步骤四,否则,认为待定点同步卫星对第i个卫星不会造成干扰,进入步骤七;
步骤四:计算第i个已定点同步卫星的定点经度ai,如果|ai-bi|≤5°,进入步骤五;否则,认为待定点同步卫星对第i个已定点同步卫星不会造成干扰,进入步骤七;其中,bi为待定点同步卫星的当前经度;
步骤五:根据第i个已定点同步卫星的TLE计算该卫星在预报弧段内的位置根据待定点同步卫星位置计算t时刻对第i个已定点同步卫星在某测站的仰角,如果仰角大于零度,即跟踪可见,进入步骤六,否则,认为待定点同步卫星对第i个已定点同步卫星不会造成干扰,进入步骤七;
步骤六:根据t时刻第i个已定点同步卫星和待定点同步卫星的位置和计算两个卫星相对测站的矢量夹角:
上式中为测站的位置;
如果θi(t)≤θ0,给定值θ0满足0.01°≤θ0≤0.1°,则认为该时刻待定点同步卫星会对第i个已定点同步卫星造成干扰;
步骤七:如果i<N,则:i=i+1,进入步骤三;否则:同步卫星早期测控段干扰弧段快速预报完成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210239598.0A CN102745343B (zh) | 2012-07-12 | 2012-07-12 | 同步卫星早期测控段干扰弧段快速预报方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210239598.0A CN102745343B (zh) | 2012-07-12 | 2012-07-12 | 同步卫星早期测控段干扰弧段快速预报方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102745343A CN102745343A (zh) | 2012-10-24 |
CN102745343B true CN102745343B (zh) | 2014-08-13 |
Family
ID=47025942
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210239598.0A Expired - Fee Related CN102745343B (zh) | 2012-07-12 | 2012-07-12 | 同步卫星早期测控段干扰弧段快速预报方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102745343B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6078286A (en) * | 1998-06-03 | 2000-06-20 | Motorola, Inc. | Method and apparatus for efficient acquisition and tracking of satellites |
FR2809831A1 (fr) * | 2000-06-06 | 2001-12-07 | Loral Space Systems Inc | Correction de l'argument du perigee avec une commande de la longitude pour des satellites geosynchrones en orbite excentrique inclinee |
CN1486545A (zh) * | 2000-11-13 | 2004-03-31 | Լ����E����ķ | 一种用于实施一个不干扰地球同步卫星群的非地球同步卫星星座的改进的系统和方法 |
DE102007041994A1 (de) * | 2007-09-04 | 2009-04-02 | Schepelmann, Jürgen | Mechanismus zum Ausschleudern einer Stoßmasse zur Ablenkung eines störenden Satelliten auf eine tolerierbare Bahn durch Kollision |
CN102126564A (zh) * | 2010-01-19 | 2011-07-20 | 泰勒斯公司 | 用于最优化人造卫星的质量的方法和装置 |
CN102494676A (zh) * | 2011-12-12 | 2012-06-13 | 中国科学院长春光学精密机械与物理研究所 | 复杂背景下的卫星自动识别装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8041118B2 (en) * | 2007-02-16 | 2011-10-18 | The Boeing Company | Pattern recognition filters for digital images |
-
2012
- 2012-07-12 CN CN201210239598.0A patent/CN102745343B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6078286A (en) * | 1998-06-03 | 2000-06-20 | Motorola, Inc. | Method and apparatus for efficient acquisition and tracking of satellites |
FR2809831A1 (fr) * | 2000-06-06 | 2001-12-07 | Loral Space Systems Inc | Correction de l'argument du perigee avec une commande de la longitude pour des satellites geosynchrones en orbite excentrique inclinee |
CN1486545A (zh) * | 2000-11-13 | 2004-03-31 | Լ����E����ķ | 一种用于实施一个不干扰地球同步卫星群的非地球同步卫星星座的改进的系统和方法 |
DE102007041994A1 (de) * | 2007-09-04 | 2009-04-02 | Schepelmann, Jürgen | Mechanismus zum Ausschleudern einer Stoßmasse zur Ablenkung eines störenden Satelliten auf eine tolerierbare Bahn durch Kollision |
CN102126564A (zh) * | 2010-01-19 | 2011-07-20 | 泰勒斯公司 | 用于最优化人造卫星的质量的方法和装置 |
CN102494676A (zh) * | 2011-12-12 | 2012-06-13 | 中国科学院长春光学精密机械与物理研究所 | 复杂背景下的卫星自动识别装置 |
Also Published As
Publication number | Publication date |
---|---|
CN102745343A (zh) | 2012-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20100138063A1 (en) | Systems and methods for interfacing renewable power sources to a power grid | |
CN106841778B (zh) | 基于pmu实现的次同步和超同步谐波参数的处理方法 | |
CN102269809B (zh) | 一种基于地形高程数据的机载气象雷达地杂波剔除方法 | |
CN110838864A (zh) | 无人值守的卫星地面站跟踪控制系统 | |
CN108802788A (zh) | 一种航向偏差的确定方法、装置、设备以及存储介质 | |
CN102901975A (zh) | 一种移动终端和获取移动终端位置信息的方法 | |
CN105626389B (zh) | 用于监测风力发电机组的塔筒状态的系统和方法 | |
CN105738918A (zh) | 基于卫星导航信息融合的输电线路位移监测系统及方法 | |
CN102565451A (zh) | 通用飞机航行风向、风速的测算方法 | |
CN112540390A (zh) | 一种航天器的跟踪预报计算方法及装置 | |
CN114397009A (zh) | 一种基于北斗差分定位的输电导线舞动监测系统及方法 | |
CN105575154A (zh) | 一种车辆gps定位丢失数据补偿方法 | |
CN102745343B (zh) | 同步卫星早期测控段干扰弧段快速预报方法 | |
CN104535974A (zh) | 一种飞机雷达系统校靶装置及其使用方法 | |
CN109856616A (zh) | 一种雷达定位相对系统误差修正方法 | |
CN106647704B (zh) | 一种geo轨道数传天线跟踪精度确定方法 | |
CN103248063B (zh) | 一种基于pmu的多直流协调广域阻尼控制方法 | |
ZHANG et al. | Integrated correction method of multi-point numerical weather prediction wind speed and irradiance | |
CN103129701A (zh) | 导流罩在分段上的定位方法 | |
JP6429523B2 (ja) | 太陽光発電量予測システム | |
CN103414184A (zh) | 一种频率变化情况下序分量的计算方法 | |
CN202216688U (zh) | 一种天线方位角检测装置 | |
CN113487100A (zh) | 光伏发电出力全局精准预测方法及系统 | |
CN105510950A (zh) | 一种基于共时钟接收机的实时精密定姿方法 | |
CN114675092B (zh) | 一种输电线路绝缘子空间电场测量方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140813 |