CN102542159A - 一种在轨航天器状态预测方法 - Google Patents
一种在轨航天器状态预测方法 Download PDFInfo
- Publication number
- CN102542159A CN102542159A CN2011104124090A CN201110412409A CN102542159A CN 102542159 A CN102542159 A CN 102542159A CN 2011104124090 A CN2011104124090 A CN 2011104124090A CN 201110412409 A CN201110412409 A CN 201110412409A CN 102542159 A CN102542159 A CN 102542159A
- Authority
- CN
- China
- Prior art keywords
- spacecraft
- telemetry
- state
- parameter
- model
- 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
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种在轨航天器状态预测方法,首先接收航天器的遥测数据,剔除遥测数据中的野值并对剔除野值后的遥测数据进行采样,生成具有时间间隔的数据。然后判断预测参数是否与航天器状态或者测控事件有关,如果有关则采用与航天器状态或者遥控事件有关的参数预测方法进行预测;如果无关则开始累积与预测参数对应的遥测数据,然后判断遥测数据是否为周期性变化,如果遥测数据为周期变化的,则采用周期变化规律的预测模型进行参数预测;如果遥测数据为非周期变化的,则采用与航天器状态或者遥控事件无关的参数预测方法进行预测。
Description
技术领域
本发明涉及一种在轨航天器状态预测方法。
背景技术
随着在轨航天器数量、设计寿命和型号种类的不断增加,维持其安全稳定运行变得越来越重要。在轨航天器长期运行在空间环境中,受到多种不确定性因素的作用,其性能与功能可能会出现变化。因此,分析在轨航天器的遥测参数变化趋势,对遥测数据进行预测并在此基础上实现在轨航天器状态预测,可以在早期及时发现遥测数据的异常变化,有效避免可能发生的重大故障,降低航天器在轨运行的风险,对于提高航天器在轨运行的安全性和可靠性具有重要的意义。
现有虽然已经有技术成熟的状态预测方法,但在轨航天器状态变化情况复杂,现有方法不能完全适应在轨航天器状态的多种变化规律,在用于在轨航天器状态进行预测时其预测的性能无法满足要求,对于某些特殊情况预测性能甚至可能出现严重下降。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供了一种采用动态预测策略的在轨航天器状态预测方法,根据数据的变化规律建立匹配机制,实现数据变化规律的识别,选择相适应的预测方法并给出高精度的预测结果,满足在轨航天器状态预测的需求。
本发明的技术解决方案是:一种在轨航天器状态预测方法,步骤如下:
(1)接收航天器的遥测数据,剔除遥测数据中的野值并对剔除野值后的遥测数据进行采样,生成具有时间间隔的数据,所述的遥测数据包括历史遥测数据和实时遥测数据;
(2)判断预测参数是否与航天器状态或者测控事件有关,如果有关则转步骤(3),采用与航天器状态或者遥控事件有关的参数预测方法进行预测;如果无关则开始累积与预测参数对应的遥测数据,然后判断遥测数据是否为周期性变化,如果遥测数据为周期变化的,则转步骤(4)采用周期变化规律的预测模型进行参数预测;如果遥测数据为非周期变化的,则转步骤(5)采用与航天器状态或者遥控事件无关的参数预测方法进行预测;
(3)判断预测参数为AR模型或是ARMA模型,通过模型检测建立时间序列模型,选取能够表征系统状态的参数来组成向量并利用时间序列模型进行状态预测;当出现连续N个奇异点后,获取与奇异点参数对应的航天器实时状态,航天器状态预测结束;如果未能获取实时状态,则确定所述奇异点是否由于航天器状态正常切换而导致,如果奇异点不是航天器状态切换而导致的,则转步骤(5)采用与航天器状态或者遥控事件无关的参数预测方法进行预测;如果奇异点是航天器状态切换而导致的,则累积数据,重新建立时间序列模型,直至预测出航天器的状态后结束;
(4)通过对各遥测参数的历史数据进行分析并建立待匹配的模型,然后用待预测参数的抽样点对待匹配的模型进行匹配和调整,直至找到最终匹配的模型并用该匹配的模型进行航天器状态的预测,航天器状态预测结束;
(5)根据遥测数据的类型选择时间序列法或者曲线拟合法进行处理;如果是采用时间序列法,则转步骤(3)进行航天器的状态预测;如果采用曲线拟合法,则利用遥测数据拟合出曲线拟合法中的相关参数后进行航天器的状态预测。
本发明与现有技术相比的优点在于:本发明方法在对遥测数据变化规律深入分析的基础上,将航天器的状态预测参数划分为与航天器运行状态和测控事件有关的参数、周期变化规律的参数、与航天器运行状态和测控事件无关的参数三种情况,并根据三种情况分别选择相适应的数据预测方法,并采取动态的预测方式,可在预测过程中根据预测的效果修改预测方法的配置,提高航天器状态预测的精度。
附图说明
图1为本发明方法的流程图;
图2为本发明与状态或事件有关的状态预测方法流程;
图3为本发明依据周期变化规律参数进行航天器状态预测的流程;
图4为本发明与状态或事件无关的状态预测方法流程。
具体实施方式
开始接收航天器状态数据以后,需要对数据进行预处理,剔除野值,并提供特定时间间隔的数据。然后需要进行状态解读,判断预测参数是否与状态或事件有关,如果有关,则调用与状态有关参数预测方法(图1中的算法1)进行预测;如果无关,累积数据,对当前的数据是否为周期变化规律进行判断。如果经过判断,判决当前数据为周期变化,则建立周期变化规律的预测模型,并根据该模型进行预测(图中1的算法2)。如果在给定的数据量内,未能判断出该数据为周期变化规律,则按与状态和事件无关的实时预测方法(图1中的算法3)的流程进行处理。
与状态有关参数预测方法流程如图2所示。首先通过数据预处理、平稳性检测,计算出AR模型或是ARMA模型的参数,并通过模型检测以及优化得到建立好的时间序列模型。之后,对遥测数据进行实时预测,当出现连续N个奇异点(拐点)后,获取实时状态,并在没有成功获取实时状态时利用模式特征识别以及模式匹配分析该拐点是否由于状态正常切换而导致。如果不是状态切换而导致的拐点,则采用与状态无关的实时预测方法进行预测,如果是,则累计样本,切换模型,如此循环下去,直到结束。
与状态有关参数预测方法步骤如下:
1)获取与状态和测控事件有关的参数设置信息,包括卫星型号、分系统、参数、预测的置信区间、距离判别函数、实时预测模型采样点数等相关的信息,以及参数状态模式配置信息等。
2)建立参数变化的时间序列模型。包括平稳性检验、白噪声检验、时间序列建模等步骤。
3)形成模式向量。选择时序方法中能够表征系统状态的参数来组成向量,用以反映系统的状态。不同方法模型下(如AR模型、ARMA模型等)选取的模式向量不同。
4)提取特征值。由于模式向量中各参数的重要性不同,同时,它们也不一定相互独立的,因此,需通过恰当的变换,选择较少的参数来形成低维向量,即特征值。
5)构造距离函数。利用时序模型的参数构造距离判别函数,作为模式特征判据。距离判别函数包括Euclide距离判别函数、残差偏移距离判别函数、Mahalaobis距离判别函数、Mann距离判别函数等。
6)确定门限值。根据系统状态的变动情况,合理的确定能将状态进行分类的门限值。例如门限值可选取为拟合值的95%置信区间。
7)状态监测和参数预测。根据前面定义好的距离判别函数和门限值,识别状态是否出现异常。分为两种情况:当状态未发现变化时,可利用前面建立的时间序列预测模型来预测配置文件中规定的时间段内的参数变化情况;当状态发现变化时,根据配置文件信息记录变化记录(如100秒或100次记录),累积数据,并结合状态解读得到的实时状态,判断是否发生正常的状态切换(如从地影转入光照状态),如果是正常的状态切换,则直接用状态对应的预测模型进行预测,如果不是正常的状态切换(即不是通过状态解读得到状态变化的信息),说明发生了变化,则调用状态无关实时预测方法进行预测(奇异点检测、实时预测模型建模预测等)。
如图3所示,建立周期变化规律的预测模型进行航天器状态预测的如下:
1)模式匹配。(1)通过对各参数历史数据的分析建立若干待匹配的模型。如多项式、正弦、余弦、指数等变化规律。(2)对于模型未知的待预测参数,提取一定数量的抽样点并在待匹配模型库中寻找并建立匹配模型。
2)模式验证。用待预测参数的抽样点检验上步计算出的模型,如果在置信区间内则进行下一步,否则返回上步继续寻找匹配模型。
3)分支匹配。确定上步验证后的模型参数。
4)再次验证。用待测参数非抽样点检验上述模型,如在置信区间内则进行下一步,否则返回到“模式匹配”寻找相匹配的模型。
5)预测。用待预测参数的已知数据和上步确定的模型进行预测。
如图4所示,与状态无关的航天器状态预测主要过程为:读入遥测数据后,通过判断是否选择的是时间序列法,如果不是,则进行数据的预处理,计算出曲线拟合法中的相关参数;如果是采用时间序列法,则进行数据预处理,并且判断处理后的数据是否平稳,平稳后,则利用该数据计算出选择的AR模型法或是ARMA模型法中的参数,并且检测该模型并进行优化。在得到模型后,根据设定的门限值进行预测。在预测过程中,根据新的实时数据,进行模型参数的优化,以提高预测的精度,如此循环,直到结束。
与状态有关参数预测方法步骤如下:
1)获取与状态和测控事件无关的参数设置信息,包括卫星型号、分系统、参数、预测的置信区间、距离判别函数、实时预测模型采样点数等相关的信息设置。
2)建立参数变化的自回归(AR)模型。包括平稳性检验、白噪声检验、时间序列建模等步骤。
3)奇异点检测。定义距离判别函数,识别是否出现异常。这里直接采用欧几里得距离函数。若根据定义好的距离函数判断出序列的相邻若干点的值连续超越置信区间,找出发生异常的时间和该时间点上的数值提前报告。
4)利用故障点附近处的观测值,建立新的预测模型,同时不断加入新的观测值,对模型进行滚动优化。
5)参数值预测。采用前述的多项式拟合外推预测或其他方法进行预测,给出一定时间后的参数预测值。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (1)
1.一种在轨航天器状态预测方法,其特征在于步骤如下:
(1)接收航天器的遥测数据,剔除遥测数据中的野值并对剔除野值后的遥测数据进行采样,生成具有时间间隔的数据,所述的遥测数据包括历史遥测数据和实时遥测数据;
(2)判断预测参数是否与航天器状态或者测控事件有关,如果有关则转步骤(3),采用与航天器状态或者遥控事件有关的参数预测方法进行预测;如果无关则开始累积与预测参数对应的遥测数据,然后判断遥测数据是否为周期性变化,如果遥测数据为周期变化的,则转步骤(4)采用周期变化规律的预测模型进行参数预测;如果遥测数据为非周期变化的,则转步骤(5)采用与航天器状态或者遥控事件无关的参数预测方法进行预测;
(3)判断预测参数为AR模型或是ARMA模型,通过模型检测建立时间序列模型,选取能够表征系统状态的参数来组成向量并利用时间序列模型进行状态预测;当出现连续N个奇异点后,获取与奇异点参数对应的航天器实时状态,航天器状态预测结束;如果未能获取实时状态,则确定所述奇异点是否由于航天器状态正常切换而导致,如果奇异点不是航天器状态切换而导致的,则转步骤(5)采用与航天器状态或者遥控事件无关的参数预测方法进行预测;如果奇异点是航天器状态切换而导致的,则累积数据,重新建立时间序列模型,直至预测出航天器的状态后结束;
(4)通过对各遥测参数的历史数据进行分析并建立待匹配的模型,然后用待预测参数的抽样点对待匹配的模型进行匹配和调整,直至找到最终匹配的模型并用该匹配的模型进行航天器状态的预测,航天器状态预测结束;
(5)根据遥测数据的类型选择时间序列法或者曲线拟合法进行处理;如果是采用时间序列法,则转步骤(3)进行航天器的状态预测;如果采用曲线拟合法,则利用遥测数据拟合出曲线拟合法中的相关参数后进行航天器的状态预测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110412409.0A CN102542159B (zh) | 2011-12-08 | 2011-12-08 | 一种在轨航天器状态预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110412409.0A CN102542159B (zh) | 2011-12-08 | 2011-12-08 | 一种在轨航天器状态预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102542159A true CN102542159A (zh) | 2012-07-04 |
CN102542159B CN102542159B (zh) | 2014-10-08 |
Family
ID=46349031
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110412409.0A Active CN102542159B (zh) | 2011-12-08 | 2011-12-08 | 一种在轨航天器状态预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102542159B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105608060A (zh) * | 2015-12-21 | 2016-05-25 | 中国银联股份有限公司 | 基于离群点检测的行业趋势预测方法及其系统 |
CN108563142A (zh) * | 2017-12-11 | 2018-09-21 | 南京航空航天大学 | 一种基于sha的航天器acs任务寿命预测方法 |
CN110348132A (zh) * | 2019-07-15 | 2019-10-18 | 长光卫星技术有限公司 | 基于Bayes估计的轨控效果快速评估方法 |
CN111612048A (zh) * | 2020-04-30 | 2020-09-01 | 中国西安卫星测控中心 | 无监督聚类异常检测的方法 |
CN111680748A (zh) * | 2020-06-08 | 2020-09-18 | 中国人民解放军63920部队 | 航天器状态模式的识别方法及识别装置 |
CN112257754A (zh) * | 2020-09-24 | 2021-01-22 | 北京航天测控技术有限公司 | 航天器运行状态的分析方法和装置 |
CN113110548A (zh) * | 2021-04-21 | 2021-07-13 | 北京控制工程研究所 | 一种航天器椭球集合演化的设计方法 |
CN115563888A (zh) * | 2022-12-06 | 2023-01-03 | 南京航空航天大学 | 一种航天器剩余寿命预测方法、系统、电子设备及介质 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106815402A (zh) * | 2016-12-13 | 2017-06-09 | 上海卫星工程研究所 | 一种全尺寸航天器结构在轨热变形预示方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6126117A (en) * | 1997-12-25 | 2000-10-03 | Nec Corporation | State control device of moving body and its state control method |
CN101859146A (zh) * | 2010-07-16 | 2010-10-13 | 哈尔滨工业大学 | 一种基于预测滤波和经验模态分解的卫星故障预测方法 |
KR101059043B1 (ko) * | 2008-12-04 | 2011-08-24 | 한국전자통신연구원 | 정지궤도 인공위성 관제를 위한 비행역학 분석계획 장치 |
CN102208028A (zh) * | 2011-05-31 | 2011-10-05 | 北京航空航天大学 | 一种适用于动态复杂系统的故障预测和诊断方法 |
-
2011
- 2011-12-08 CN CN201110412409.0A patent/CN102542159B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6126117A (en) * | 1997-12-25 | 2000-10-03 | Nec Corporation | State control device of moving body and its state control method |
KR101059043B1 (ko) * | 2008-12-04 | 2011-08-24 | 한국전자통신연구원 | 정지궤도 인공위성 관제를 위한 비행역학 분석계획 장치 |
CN101859146A (zh) * | 2010-07-16 | 2010-10-13 | 哈尔滨工业大学 | 一种基于预测滤波和经验模态分解的卫星故障预测方法 |
CN102208028A (zh) * | 2011-05-31 | 2011-10-05 | 北京航空航天大学 | 一种适用于动态复杂系统的故障预测和诊断方法 |
Non-Patent Citations (3)
Title |
---|
TAKEHISA YAIRI ET AL: "Fault Detection by Mining Association Rules from House-keeping Data", 《PROCEEDING OF THE 6TH INTERNATIONAL SYMPOSIUM ON ARTIFICIAL INTELLIGENCE AND ROBOTICS & AUTOMATION IN SPACE:I-SAIRAS 2001》 * |
房红征等: "一种与状态和事件相关的卫星实时故障预警方法", 《电子测量与仪器学报》 * |
秦巍: "一种基于历史遥测数据的在轨卫星故障预警系统", 《航天器工程》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105608060A (zh) * | 2015-12-21 | 2016-05-25 | 中国银联股份有限公司 | 基于离群点检测的行业趋势预测方法及其系统 |
CN108563142A (zh) * | 2017-12-11 | 2018-09-21 | 南京航空航天大学 | 一种基于sha的航天器acs任务寿命预测方法 |
CN108563142B (zh) * | 2017-12-11 | 2020-07-24 | 南京航空航天大学 | 一种基于sha的航天器acs任务寿命预测方法 |
CN110348132A (zh) * | 2019-07-15 | 2019-10-18 | 长光卫星技术有限公司 | 基于Bayes估计的轨控效果快速评估方法 |
CN111612048B (zh) * | 2020-04-30 | 2023-03-10 | 中国西安卫星测控中心 | 无监督聚类异常检测的方法 |
CN111612048A (zh) * | 2020-04-30 | 2020-09-01 | 中国西安卫星测控中心 | 无监督聚类异常检测的方法 |
CN111680748A (zh) * | 2020-06-08 | 2020-09-18 | 中国人民解放军63920部队 | 航天器状态模式的识别方法及识别装置 |
CN111680748B (zh) * | 2020-06-08 | 2024-02-02 | 中国人民解放军63920部队 | 航天器状态模式的识别方法及识别装置 |
CN112257754A (zh) * | 2020-09-24 | 2021-01-22 | 北京航天测控技术有限公司 | 航天器运行状态的分析方法和装置 |
CN112257754B (zh) * | 2020-09-24 | 2023-07-28 | 北京航天测控技术有限公司 | 航天器运行状态的分析方法和装置 |
CN113110548A (zh) * | 2021-04-21 | 2021-07-13 | 北京控制工程研究所 | 一种航天器椭球集合演化的设计方法 |
CN115563888A (zh) * | 2022-12-06 | 2023-01-03 | 南京航空航天大学 | 一种航天器剩余寿命预测方法、系统、电子设备及介质 |
CN115563888B (zh) * | 2022-12-06 | 2023-04-07 | 南京航空航天大学 | 一种航天器剩余寿命预测方法、系统、电子设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN102542159B (zh) | 2014-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102542159A (zh) | 一种在轨航天器状态预测方法 | |
Wen et al. | Degradation modeling and RUL prediction using Wiener process subject to multiple change points and unit heterogeneity | |
CN110678816B (zh) | 用于控制技术系统的方法和控制机构 | |
CN103646167A (zh) | 一种基于遥测数据的卫星异常状态检测系统 | |
CN107895014B (zh) | 一种基于MapReduce框架的时间序列桥梁监测数据分析方法 | |
CN109767054A (zh) | 基于深度神经网络算法的能效云评估方法及边缘能效网关 | |
CN104156615A (zh) | 基于ls-svm的传感器检测数据点异常检测方法 | |
KR101223898B1 (ko) | 태양광 발전모듈의 고장 진단 및 예지 방법 | |
CN102222151A (zh) | 一种基于自回归滑动平均的模拟电路故障预测方法 | |
Cui et al. | An optimized swinging door algorithm for wind power ramp event detection | |
Ak et al. | Data analytics and uncertainty quantification for energy prediction in manufacturing | |
Zha et al. | Selection of time window for wind power ramp prediction based on risk model | |
CN113075512A (zh) | 基于声检测的变压器放电故障诊断方法和系统 | |
Armon et al. | Algorithmic network monitoring for a modern water utility: a case study in Jerusalem | |
Hosseini et al. | Traffic flow prediction using MI algorithm and considering noisy and data loss conditions: An application to Minnesota traffic flow prediction | |
CN111080484A (zh) | 一种配电网异常数据监测方法及装置 | |
CN103914630A (zh) | 建筑能耗监管系统 | |
CN102799108A (zh) | 基于低频能耗信息的能耗解聚合方法及家庭能源管理系统 | |
CN116720983A (zh) | 一种基于大数据分析的供电设备异常检测方法及系统 | |
CN104680010A (zh) | 一种汽轮机组稳态运行数据筛选方法 | |
CN117172601A (zh) | 一种基于残差全卷积神经网络的非侵入式负荷监测方法 | |
CN115469585B (zh) | 一种基于大数据的水电机组状态监测方法及系统 | |
Yang et al. | Two-stage prediction technique for rolling bearings based on adaptive prediction model | |
CN101923605B (zh) | 铁路防灾风预警方法 | |
Gu et al. | Process analysis for performance evaluation of prognostics methods orienting to engineering application |
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 |