CN106053978A - Satellite non-periodic remote measurement analog quantity interpretation method based on windows - Google Patents
Satellite non-periodic remote measurement analog quantity interpretation method based on windows Download PDFInfo
- Publication number
- CN106053978A CN106053978A CN201610346317.XA CN201610346317A CN106053978A CN 106053978 A CN106053978 A CN 106053978A CN 201610346317 A CN201610346317 A CN 201610346317A CN 106053978 A CN106053978 A CN 106053978A
- Authority
- CN
- China
- Prior art keywords
- window
- data
- analog quantity
- stable state
- satellite
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
Abstract
Provided is a satellite non-periodic remote measurement analog quantity interpretation method based on windows, comprising: setting a window width; creating continuous windows, searching for a first stable state and calculating a mean value; detecting subsequent windows, and searching for a first change state window; continuing detecting subsequent windows, finding a second stable state and calculating a mean value; calculating the mean value between the first and second stable states to be used as the variation of analog quantity; retrieving external stimulation logs, determining whether external stimulation corresponding to the variation of analog quantity exists, and determining the variation of analog quantity to be anomalous change if external stimulation corresponding to the variation of analog quantity does not exist; and using the second stable state as a new stable state, and repeating the steps to realize automatic interpretation of satellite non-periodic remote measurement analog quantity. The satellite non-periodic remote measurement analog quantity interpretation method can automatically identify the change (increase or decrease) and variation of non-periodic remote measurement analog quantity through software, and discover abnormal change caused by equipment faults.
Description
Technical field
The present invention relates to a kind of satellite aperiodicity remote measurement analog quantity interpretation method based on window, be used for realizing satellite distant
Survey the intelligent interpretation of analog quantity, belong to satellite test field.
Background technology
Satellite, during electrical property integration test, passes down the telemetry of whole star by observing and controlling passage.Defend
The telemetry of star includes digital quantity and the big class of analog quantity two, and wherein analog quantity is divided into again static simulation amount, periodicity modulus analog quantity
With aperiodicity analog quantity.Tester needs these data to carry out real-time supervision and interpretation to determine whether satellite works
Normally.Digital quantity remote measurement represents the duty of on-board equipment, and meaning is obvious, direct, and numerical value is generally solid at a limited number of
Definite value changes.The span of static simulation amount is constant, then shows occur extremely beyond threshold value.Periodicity modulus analog quantity then has
Strict period of change and rule.
At present, satellite telemetry analog quantity is mainly based on artificial interpretation, higher to knowledge, the skill requirement of tester,
And it is easily generated during interpretation and fails to judge, judge by accident.Simultaneously along with the development of data analysis technique, satellite telemetering data interpretation
Gradually change to software automation, intelligent interpretation.Along with the development of data analysis technique, satellite telemetering data interpretation by
Gradually change to software automation, intelligent interpretation.Digital quantity, static simulation amount and periodicity modulus analog quantity remote measurement have significantly
Data characteristics or Changing Pattern, be easier to realize intelligent interpretation by software.
And for aperiodicity analog quantity, owing to changing random feature, it is achieved software intelligent interpretation is the most difficult.With whole
As a example by star load current remote measurement, as an important analog quantity remote measurement of the whole star duty of sign, its change in value has two
Kind of reason, one is outside normal excitation, and including control instruction, environmental change etc., two is equipment fault, due to external drive or set
Standby fault is random, irregular, utilizes software automatically to identify aperiodicity analog quantity so there is presently no
Method.
Summary of the invention
The technology of the present invention solves problem: overcome the deficiencies in the prior art, it is provided that a kind of satellite of based on window aperiodic
Property remote measurement analog quantity interpretation method, can automatically identify the change of aperiodicity remote measurement analog quantity by software and (increase or subtract
Little) and variable quantity, and it can be found that the ANOMALOUS VARIATIONS that caused by equipment fault.
The technical solution of the present invention: a kind of satellite aperiodicity remote measurement analog quantity interpretation method based on window, step
Rapid as follows:
(1) base area arranges window width in the face of the sampling period of satellite telemetry analog quantity, and described window width is n and adopts
The sample cycle, n >=3；
(2) in time series data, continuous print window, the satellite telemetry number gathered according to ground are created at each window
According to finding first stable state, described first stable state refers to judge first the window of the data stabilization obtained, and calculates first surely
The data mean value of state；
(3) continue detection subsequent window, when certain window data instability, i.e. find a change state window, according to
The equal value difference of this window and previous window obtains the change information of data；
(4) repeating step (3), until certain window data is stable, this window is second stable state, calculates second surely
The data mean value of state；
(5) calculate first stable state and the equal value difference of second stable state, obtain change and the variable quantity of this analog quantity；
(6) retrieval external drive daily record, if without the external drive of corresponding this analog quantity change, then this is abnormal change
Changing, if there being the external drive of corresponding this analog quantity change, then this is normal variation；
(7) using second stable state as first new stable state, step (3) to (6) is repeated, it is achieved to satellite aperiodicity
The automatic interpretation of remote measurement analog quantity.
Window data is the most stable to utilize adjacent window apertures equal value difference method to judge, specific implementation is as follows:
Making Lim is steady-state deviation threshold value set in advance, when Lim is absolute value, judges according to following method:
If | Ew2-Ew1| < Lim, then second window data is stable, otherwise second window data instability；
Wherein Ew1It is the data mean value of first window, Ew2It is the data mean value of second window, first window and
Two windows are adjacent window apertures；
When Lim is percentage ratio, judge according to following method:
If | Ew2-Ew1|/Ew1* 100% < Lim, then second window data is stable, otherwise second window data shakiness
Fixed.
Window data is the most stable to utilize reverse test method to judge, specific implementation is as follows:,
(3.1) data sequence in current window is divided into M section, then obtains the average of every segment data, be designated as y1,
y2..., yM；
(3.2) y is calculatediPermutation number Ai, described AiEqual to yi+1To yMIntermediate value is more than yiData amount check；
(3.3) calculating backward sum A, described backward sum A is A1, A2..., AM-1Sum；
(3.4) calculating the expectation of A, variance and statistic, described A is desired for E (A)=M (M-1)/4, and variance is D (A)
=M (2M2+ 3M-5)/72, statistic Z=[A+0.5-E (A)]/(D (A))1/2, Z submits to N (0,1) distribution；
(3.5) in the case of level of significance α=0.05, if | Z | < 1.96, then it is assumed that this sequence is stationary sequence, i.e. when
Front window data stabilization.
Present invention advantage compared with prior art is:
(1) present invention can recognize that the change (being increased or decreased) of analog data, and can be given from first steady
State is to the variable quantity of second stable state, it is possible to rejects the change of normal analog data, identifies abnormal change, it is simple to prompting and
Report to the police, it is achieved that the intelligent interpretation of aperiodicity remote measurement analog quantity, it is possible to substitute artificial interpretation, it is to avoid fail to judge.
(2) present invention uses adjacent window apertures equal value difference method to judge stablizing of window data for window width in the case of less
Property, the method is applicable to change satellite telemetry analog data faster, and during such as satellite speed becomes remote measurement, the sampling period is 1 second
Load current, appearance control electric current etc., such data complete a next state change and only need several or ten several sampling periods, by arranging
Window width is less value, uses adjacent window apertures equal value difference method can identify stable state and change state fast and accurately, in time
Make prompting or report to the police.
(3) present invention uses reverse test method to judge the stability of window data for window width in the case of bigger.Should
Method is applicable to change the remote temperature sensing in the gradual data of satellite telemetry analog data slowly, such as satellite, sampling period
Generally 16 seconds or more, (sample 112 times) in 30 minutes and may only generate less than the change of 2 degree.This situation i.e. can set
Put big width window (n=100, M=10).Change can be identified accurately slow by big width window and reverse test method
The variation tendency of slow telemetry, can be in the way of substituting existing artificial interpretation.
Accompanying drawing explanation
Fig. 1 is the method flow diagram of the present invention；
Fig. 2 is the data analysis schematic diagram of embodiment of the present invention Satellite load current；
Detailed description of the invention
The invention provides a kind of Satellite Simulation amount interpretation method based on window, by setting up time series data window
Aperiodicity analog quantity is analyzed, identifies change (increase or reduce) and the variable quantity of analog quantity, and can send out
The ANOMALOUS VARIATIONS now caused by equipment fault.
As it is shown in figure 1, specifically comprise the following steps that
(1) window width and steady-state deviation threshold value are set.Base area was arranged in the face of the sampling period of satellite telemetry analog quantity
Window width, described window width is n sampling period, n >=3.When analog quantity pace of change is very fast, window width should take
One less value, such interpretation reaction can be very fast.Such as whole star load current, its sampling period is 1s, can arrange
Window width is 5s, and after change occurs in electric current, 5 sampling intervals of delay i.e. can complete to identify and interpretation.When analog quantity changes
When speed is slower, window width should take a bigger value, and such interpretation just can be more accurate.During the i.e. stable state of steady-state deviation threshold value
Adjacent two admissible mean bias of window, can be a percentage ratio, it is also possible to be an absolute value, when adjacent window apertures
When all value difference is in this threshold range, it is determined that in window, data are in stable state.
(2) in time series data, create continuous window, search first stable state, the average of first stable state of record.
In judging window, data the most stable (being in stable state) have two kinds of methods, i.e. adjacent window apertures equal value difference method and reverse test method.
Adjacent window apertures equal value difference method is applicable to the situation that window width is less.Steady-state deviation threshold value is formula during absolute value
As follows:
If | Ew2-Ew1| < Lim, then second window data is stable, otherwise second window data instability；
Wherein Ew1It is the data mean value of first window, Ew2It is the data mean value of second window, first window and
Two windows are adjacent window apertures, and Lim is steady-state deviation threshold value absolute value set in advance.
Sometimes steady-state deviation threshold value is given with percents, and steady-state deviation threshold value is that formula during percentage ratio is as follows:
|Ew2-Ew1|/Ew1* 100% < Lim, then second window data is stable, otherwise second window data instability.
Reverse test method is applicable to the situation that window width is bigger.First data sequence in window is divided into M section, then asks
Go out the average of every segment data, be designated as y1, y2..., yM。yiPermutation number AiEqual to yi+1To yMMiddle numerical value is more than yiNumber, i.e. Ai
Equal to yjNumber, yjMeet yj>yi(j > i)) (i.e. with from big to small as standard order).Calculate A1, A2..., AM-1Sum is
Backward sum, is designated as A.A is desired for
E (A)=M (M-1)/4,
Variance is D (A)=M (2M2+ 3M-5)/72,
Statistic Z=[A+0.5-E (A)]/(D (A))1/2,
Z approximation submits to N (0,1) distribution, in the case of level of significance α=0.05, if | Z | < 1.96, then it is assumed that this sequence
Be classified as stationary sequence, i.e. current window is stable state.
(3) detection subsequent window is continued, when adjacent window apertures average not grade or reverse test method result video data are unstable
Time, it is i.e. to find a change state window, is appreciated that according to the equal value difference of this window Yu previous window the change of data (increases
Or reduce)；
(4) step (3) is repeated, until adjacent window apertures average is equal or reverse test method result display window data steady,
It is i.e. to find second stable state, records second stable state average；
(5) change and the variable quantity of this analog quantity is drawn according to the equal value difference of first stable state, second stable state；
(6) retrieval external drive daily record, if without the external drive of corresponding this analog quantity change, then this is abnormal change
Changing, otherwise this is normal variation；
(7) using second stable state as first new stable state, (3) to (6) step is repeated, it is achieved to satellite aperiodicity
The automatic interpretation of remote measurement analog quantity.
Embodiment:
Fig. 2 is a change curve of certain satellite load electric current, 1 second sampling period, wherein the data of 1～50s such as table 1 institute
Showing, data unit is ampere (A).Begin with adjacent window apertures equal value difference method from data 11 and carry out data interpretation, window width is set
Degree is 5 sampling periods, and threshold value is absolute value 0.2, and the result of window calculation is as shown in table 2.
The sampled data of certain satellite load electric current of table 1
| 1～10s | 1.6 | 1.7 | 1.75 | 1.6 | 1.7 | 1.75 | 1.6 | 1.7 | 1.75 | 1.6 |
| 11～20s | 1.7 | 1.75 | 1.6 | 1.7 | 1.75 | 1.6 | 1.7 | 1.75 | 1.6 | 1.7 |
| 21～30s | 1.75 | 2 | 2.1 | 2.2 | 3.5 | 3.6 | 3.8 | 3.7 | 3.8 | 3.85 |
| 31～40s | 3.9 | 4.1 | 4.5 | 4.9 | 5.2 | 5.3 | 5.2 | 5.3 | 5.1 | 5.2 |
| 41～50s | 5.3 | 5.1 | 5.2 | 5.3 | 5.1 | 5.2 | 5.3 | 5.1 | 5.2 | 5.3 |
The average of table 2 window and equal value difference
Equal value difference according to adjacent window apertures, it is possible to determine that the 2nd window is first stable state, the 3rd～6 window is change
State, the 7th window is second stable state.The equal value difference recording first stable state and second stable state is+3.53, then this satellite is born
Carrying the telemetering of current interpretation conclusion in this time period is " load current increases 3.53A ".Can be by retrieval external drive day
Will, checks whether this is ANOMALOUS VARIATIONS.Using second stable state as first new stable state.Repeat said method, to 45s with
After data continue calculate stable state and change state, it is achieved the interpretation of satellite load electric current.
The inventive method can carry out intelligentized interpretation to satellite telemetry analog data aperiodic, detects analog quantity
Variation tendency and variable quantity, in combination with external drive daily record, reject normal data variation, ANOMALOUS VARIATIONS carried
Show warning, it is possible to well meet the data monitoring requirement during satellite test, thus reach satellite telemetry simulation aperiodic
Automatization's interpretation of amount.
The non-detailed description of the present invention is known to the skilled person technology.
Claims (3)
1. a satellite aperiodicity remote measurement analog quantity interpretation method based on window, it is characterised in that step is as follows:
(1) base area arranges window width in the face of the sampling period of satellite telemetry analog quantity, and described window width is n sampling week
Phase, n >=3；
(2) in time series data, create continuous print window, seek according to the satellite telemetering data that ground gathers at each window
Look for the window of the data stabilization that first stable state, described first stable state refer to judge first obtain, calculate first stable state
Data mean value；
(3) continue detection subsequent window, when certain window data instability, i.e. find a change state window, according to this window
Mouthful with the change information of the equal value difference of previous window acquisition data；
(4) repeating step (3), until certain window data is stable, this window is second stable state, calculates second stable state
Data mean value；
(5) calculate first stable state and the equal value difference of second stable state, obtain change and the variable quantity of this analog quantity；
(6) retrieval external drive daily record, if without the external drive of corresponding this analog quantity change, then this is ANOMALOUS VARIATIONS, if
Have the external drive of corresponding this analog quantity change, then this is normal variation；
(7) using second stable state as first new stable state, step (3) to (6) is repeated, it is achieved to satellite aperiodicity remote measurement
The automatic interpretation of analog quantity.
A kind of satellite aperiodicity remote measurement analog quantity interpretation method based on window the most according to claim 1, its feature
It is: window data is the most stable to utilize adjacent window apertures equal value difference method to judge, specific implementation is as follows:
Making Lim is steady-state deviation threshold value set in advance, when Lim is absolute value, judges according to following method:
If | Ew2-Ew1| < Lim, then second window data is stable, otherwise second window data instability；
Wherein Ew1It is the data mean value of first window, Ew2It is the data mean value of second window, first window and second
Window is adjacent window apertures；
When Lim is percentage ratio, judge according to following method:
If | Ew2-Ew1|/Ew1* 100% < Lim, then second window data is stable, otherwise second window data instability.
A kind of satellite aperiodicity remote measurement analog quantity interpretation method based on window the most according to claim 1, its feature
It is: window data is the most stable to utilize reverse test method to judge, specific implementation is as follows:,
(3.1) data sequence in current window is divided into M section, then obtains the average of every segment data, be designated as y1, y2..., yM；
(3.2) y is calculatediPermutation number Ai, described AiEqual to yi+1To yMIntermediate value is more than yiData amount check；
(3.3) calculating backward sum A, described backward sum A is A1, A2..., AM-1Sum；
(3.4) calculating the expectation of A, variance and statistic, described A is desired for E (A)=M (M-1)/4, and variance is D (A)=M
(2M2+ 3M-5)/72, statistic Z=[A+0.5-E (A)]/(D (A))1/2, Z submits to N (0,1) distribution；
(3.5) in the case of level of significance α=0.05, if | Z | < 1.96, then it is assumed that this sequence is stationary sequence, i.e. current window
Mouth data stabilization.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610346317.XA CN106053978B (en) | 2016-05-23 | 2016-05-23 | A kind of satellite aperiodicity telemetering analog quantity interpretation method based on window |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610346317.XA CN106053978B (en) | 2016-05-23 | 2016-05-23 | A kind of satellite aperiodicity telemetering analog quantity interpretation method based on window |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106053978A true CN106053978A (en) | 2016-10-26 |
| CN106053978B CN106053978B (en) | 2018-12-21 |
Family
ID=57174320
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610346317.XA Active CN106053978B (en) | 2016-05-23 | 2016-05-23 | A kind of satellite aperiodicity telemetering analog quantity interpretation method based on window |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106053978B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106778839A (en) * | 2016-11-30 | 2017-05-31 | 哈尔滨工业大学 | The satellite telemetering data time series extracted based on local extremum optimization represents method |
| CN107525686A (en) * | 2017-08-21 | 2017-12-29 | 中国人民解放军63789部队 | A kind of more telemetry parameter analysis window division methods of spacecraft |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4980844A (en) * | 1988-05-27 | 1990-12-25 | Victor Demjanenko | Method and apparatus for diagnosing the state of a machine |
| CN101718864A (en) * | 2009-11-13 | 2010-06-02 | 航天东方红卫星有限公司 | System for real-time monitoring abnormal change of satellite telemetry parameters |
| CN103646167A (en) * | 2013-11-22 | 2014-03-19 | 北京空间飞行器总体设计部 | Satellite abnormal condition detection system based on telemeasuring data |
| CN104155130A (en) * | 2014-07-21 | 2014-11-19 | 航天东方红卫星有限公司 | Comprehensive test intelligent interpretation system for small satellite |
| CN104483836A (en) * | 2014-11-17 | 2015-04-01 | 国家电网公司 | Online steady state processing method of telemetry data |
| CN105242534A (en) * | 2015-09-11 | 2016-01-13 | 中国人民解放军国防科学技术大学 | Telemetry parameter and correlation with satellite control behavior-based satellite condition monitoring method |
| CN105637502A (en) * | 2013-08-14 | 2016-06-01 | 视觉股份公司 | Method and apparatus for determining a smoothed data point within a stream of data points |
-
2016
- 2016-05-23 CN CN201610346317.XA patent/CN106053978B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4980844A (en) * | 1988-05-27 | 1990-12-25 | Victor Demjanenko | Method and apparatus for diagnosing the state of a machine |
| CN101718864A (en) * | 2009-11-13 | 2010-06-02 | 航天东方红卫星有限公司 | System for real-time monitoring abnormal change of satellite telemetry parameters |
| CN105637502A (en) * | 2013-08-14 | 2016-06-01 | 视觉股份公司 | Method and apparatus for determining a smoothed data point within a stream of data points |
| CN103646167A (en) * | 2013-11-22 | 2014-03-19 | 北京空间飞行器总体设计部 | Satellite abnormal condition detection system based on telemeasuring data |
| CN104155130A (en) * | 2014-07-21 | 2014-11-19 | 航天东方红卫星有限公司 | Comprehensive test intelligent interpretation system for small satellite |
| CN104483836A (en) * | 2014-11-17 | 2015-04-01 | 国家电网公司 | Online steady state processing method of telemetry data |
| CN105242534A (en) * | 2015-09-11 | 2016-01-13 | 中国人民解放军国防科学技术大学 | Telemetry parameter and correlation with satellite control behavior-based satellite condition monitoring method |
Non-Patent Citations (1)
| Title |
|---|
| 韩冬等: "基于实时数据库的小卫星数据监视工具设计", 《航天器工程》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106778839A (en) * | 2016-11-30 | 2017-05-31 | 哈尔滨工业大学 | The satellite telemetering data time series extracted based on local extremum optimization represents method |
| CN106778839B (en) * | 2016-11-30 | 2021-02-26 | 哈尔滨工业大学 | Satellite telemetry data time sequence representation method based on local extremum optimization extraction |
| CN107525686A (en) * | 2017-08-21 | 2017-12-29 | 中国人民解放军63789部队 | A kind of more telemetry parameter analysis window division methods of spacecraft |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106053978B (en) | 2018-12-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109657937A (en) | A kind of Reliability Assessment and life-span prediction method based on degraded data | |
| US10482204B2 (en) | System for processing data and modelling for analysis of the energy consumption of a site | |
| CN111669123B (en) | Method and device for fault diagnosis of photovoltaic string | |
| CN110209560A (en) | Data exception detection method and detection device | |
| CN109345408A (en) | Abnormal data discriminating method based on shot and long term memory network | |
| CN104880605A (en) | Battery management system and battery monomer DC internal resistance measuring method | |
| CN109063885A (en) | A kind of substation's exception metric data prediction technique | |
| CN110580492A (en) | Track circuit fault precursor discovery method based on small fluctuation detection | |
| CN106053978A (en) | Satellite non-periodic remote measurement analog quantity interpretation method based on windows | |
| CN106935038B (en) | Parking detection system and detection method | |
| CN107135230A (en) | A kind of network flow abnormal detecting method and equipment based on INFIGARCH models | |
| CN107222169A (en) | A kind of photovoltaic array based on outlier detection is monitored and fault detection method in real time | |
| CN115469585B (en) | Hydropower unit state monitoring method and system based on big data | |
| CN110552682A (en) | Descending analysis method for vertical well yield of carbonate rock triple medium gas reservoir | |
| CN106646106A (en) | Power grid fault detection method based on change point detection technology | |
| CN105158607B (en) | A kind of dependent multichannel analog amount method of real-time | |
| CN117411436B (en) | Photovoltaic module state detection method, system and storage medium | |
| CN113158474B (en) | Uncertainty circuit reliability analysis method based on Caputo type fractional order differentiation | |
| Eisenbarth et al. | A data-driven method for estimating parameter uncertainty in PMU-based Power Plant Model Validation | |
| CN113221492B (en) | Method for analyzing internal transient state of digital circuit under electromagnetic interference | |
| CN113395031B (en) | Complex system fault diagnosis method and system based on multistage model | |
| CN116245357B (en) | Fault diagnosis method and system for intelligent energy-saving cabinet | |
| CN112464848B (en) | Information flow abnormal data monitoring method and device based on density space clustering | |
| CN109556649A (en) | A kind of signal acquisition monitoring method of intelligence sensor | |
| CN115982665B (en) | Quality anomaly auditing method and system for water turbine measurement data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |