WO2019231039A1 - 통제 불가능한 인공우주물체 재진입 예측 방법 - Google Patents

통제 불가능한 인공우주물체 재진입 예측 방법 Download PDF

Info

Publication number
WO2019231039A1
WO2019231039A1 PCT/KR2018/007053 KR2018007053W WO2019231039A1 WO 2019231039 A1 WO2019231039 A1 WO 2019231039A1 KR 2018007053 W KR2018007053 W KR 2018007053W WO 2019231039 A1 WO2019231039 A1 WO 2019231039A1
Authority
WO
WIPO (PCT)
Prior art keywords
time point
orbital
space object
artificial space
value
Prior art date
Application number
PCT/KR2018/007053
Other languages
English (en)
French (fr)
Inventor
최은정
Original Assignee
한국 천문 연구원
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 한국 천문 연구원 filed Critical 한국 천문 연구원
Publication of WO2019231039A1 publication Critical patent/WO2019231039A1/ko

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G3/00Observing or tracking cosmonautic vehicles
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G7/00Simulating cosmonautic conditions, e.g. for conditioning crews
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13Differential equations
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/10Artificial satellites; Systems of such satellites; Interplanetary vehicles
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G99/00Subject matter not provided for in other groups of this subclass
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/24Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/62Systems for re-entry into the earth's atmosphere; Retarding or landing devices
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Definitions

  • the present invention relates to an artificial space object reentry prediction method, and more particularly, to an uncontrollable artificial space object reentry prediction method by applying a Draw Scale Factor Estimation (DSFE) method.
  • DSFE Draw Scale Factor Estimation
  • the atmosphere reentry process starts within about one month, and when the artificial space object weighing more than 1 ton falls, about 10-40% of the debris reaches the surface.
  • the fall of uncontrollable artificial space objects is difficult to predict, which can lead to damage to life and property on the ground. Therefore, in order to prepare for space risks, the technology of predicting the fall risk of space objects is essential to minimize such risks.
  • An object of the present invention to solve the above problems is to provide an uncontrollable artificial space object re-entry prediction method that accurately predicts the predicted time and point of the fall by applying the atmospheric correction ratio estimation method (DSFE).
  • DSFE atmospheric correction ratio estimation method
  • the uncontrollable artificial space object reentry prediction method of the present invention calculates an average long radius and a latitude factor by inputting a contact orbital element or a mean orbital element (TLE) of two different viewpoints of an artificial space object. Doing; Average long radius at the second time point predicted by performing orbital propagation through the Cowell method's high-precision orbital propagator using numerical integration to the second time point by applying the initial atmospheric correction ratio value to the second time point, Calculating latitude factor and atmospheric drag; The optimum atmospheric correction ratio value is changed while changing the atmospheric correction ratio value until the error is smaller than an arbitrary convergence value by comparing the predicted average long radius or latitude factor value with the average long radius or latitude factor value of the second input time.
  • TLE mean orbital element
  • the average orbital element TLE may be converted into a contact orbital element to calculate an average orbit in a true of date coordinate system.
  • the convergence value may be a position error arbitrarily determined by the user.
  • the present invention by applying the atmospheric correction ratio estimation, it is possible to accurately predict the time and point of the atmosphere reentry of the uncontrollable artificial space object.
  • FIG. 1 is a flowchart illustrating a method for predicting the re-entry of an uncontrolled artificial space object according to a first embodiment of the present invention.
  • FIG. 2 is a flowchart illustrating a method for predicting reentry of an uncontrolled artificial space object according to a second embodiment of the present invention.
  • FIG. 3 is a flowchart illustrating a method for predicting reentry of an uncontrolled artificial space object according to a third embodiment of the present invention.
  • 1 to 3 are flowcharts illustrating an uncontrollable artificial space object re-entry prediction method according to the first to third embodiments of the present invention.
  • an uncontrollable artificial space object reentry prediction method first includes two different viewpoints of the artificial space object ( , ) Orbital element , ) To enter your average long radius ( , ) And latitude arguments ( , ) Is calculated (S100, S200, S300).
  • the track element is a contact track element or an average track element (TLE), in the case of an average track element can be converted to a contact track element to calculate the average trajectory in the TOD (True of Date) coordinate system.
  • the first time point ( Orbital information of the )
  • a second time point To perform orbit propagation (S110, S210, S310).
  • the orbital propagation at this time is the second time point predicted through the high precision orbital propagator of the Cowell method using numerical integration ( Orbital element in ) Mean longitude () ) And latitude arguments ( ) And atmospheric drag ( Calculate here, Is the drag coefficient, A is the cross-sectional area, m is the mass, ⁇ is the airtightness, Vector speed, Is the velocity vector magnitude.
  • the high precision orbital propagator of the Cowell method considers the position of the artificial space object at any time in consideration of all the perturbation forces affecting the artificial space object such as the earth's gravitational field, atmospheric influence, solar and moon attraction, and solar radiation pressure. Since the algorithm for obtaining the speed is well known in the art, a detailed description thereof will be omitted.
  • the optimum atmospheric correction ratio value S220, S320
  • S140, S240, S340 otherwise change the atmospheric correction ratio value (S130, S230, S330), and repeats the process.
  • the first point in time At the second point in time ( Atmospheric correction percentage value () ) Reflects the predicted average long radius ( ) Or latitude argument ( ) And the second time point initially entered ( ) 'S average long radius ( ) Or latitude argument ( ) And change the atmospheric correction ratio value until the error is smaller than the convergence value, Find).
  • the convergence value is a position error arbitrarily determined by the user, for example, 10 ⁇ 4 km.
  • the optimized atmospheric correction rate value ( ) is predicted (S150, S250, S350).

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Mathematical Physics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Astronomy & Astrophysics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Automation & Control Theory (AREA)
  • Algebra (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

본 발명은 통제 불가능한 인공우주물체 재진입 예측 방법에 관한 것으로, 인공우주물체의 다른 두 시점의 접촉궤도요소 또는 평균궤도요소(TLE)를 입력하여 평균 장반경과 위도 인수를 계산하는 단계; 상기 첫 번째 시점의 궤도정보를 임의 상수인 초기 대기보정 비율 값을 적용하여 두 번째 시점까지 수치적분을 이용한 코웰 방법의 고정밀 궤도 전파기를 통해 궤도 전파를 수행하여 예측된 두 번째 시점에서의 평균 장반경, 위도인수 및 대기항력을 계산하는 단계; 상기 예측된 평균 장반경 또는 위도인수 값과 초기 입력된 두 번째 시점의 평균 장반경 또는 위도 인수 값을 비교하여 오차가 임의의 수렴 값보다 작아질 때까지 대기 보정 비율 값을 변경하면서 최적의 대기 보정 비율 값을 결정하는 단계; 및 상기 최적화된 대기 보정 비율 값을 적용하여 두 번째 시점부터 재진입 시점까지 수치적분을 이용한 코웰 방법의 고정밀 궤도 전파기를 통해 궤도 예측을 수행하여 인공우주물체의 대기권 재진입 시점 및 지점을 예측하는 단계;를 포함하는 것을 특징으로 한다.

Description

통제 불가능한 인공우주물체 재진입 예측 방법
본 발명은 인공우주물체 재진입 예측 방법에 관한 것으로, 더욱 상세하게는 대기 보정 비율 추정 방법(Drag Scale Factor Estimation; DSFE)을 적용하여 통제 불가능한 인공우주물체 재진입 예측 방법에 관한 것이다.
1톤 이상의 인공우주물체 중 통제가 불가능한 인공우주물체의 추락은 지상 피해가 발생할 가능성이 높기 때문에 인공우주물체의 고도가 250km이하로 되는 경우 국가적으로도 인공우주물체 추락 및 충돌 대응 매뉴얼에 따라 우주물체 추락 상황에 대한 위기경보를 발령한다. 따라서 인공우주물체 추락 상황 및 위험도를 예측하기 위해서는 신속하고 정확한 재진입 예측 정보를 제공하는 것이 매우 중요하다.
특히, 고도 250km에 이르면 약 한 달 이내에 대기권 재진입 과정을 시작하게되며, 1톤 이상의 무게를 갖는 인공우주물체의 추락 시에는 약 10~40%의 파편들이 지표면에 도달하게 된다. 그 중에서도 통제가 불가능한 인공우주물체의 추락은 예측이 어려워, 지상에서의 인명 및 자산 피해로 이어질 수 있다. 따라서 우주위험 대비를 위해서는 이러한 위험을 최소화하기 위한 우주물체의 추락 위험 예측 기술이 반드시 필요하다.
그러나, 이러한 통제 불가능한 인공우주물체 재진입 예측의 기존 방법은 TLE 평균궤도를 이용한 SGP4(Simplified General Perturbations 4) 궤도 전파기를 사용하여 재진입 시점을 예측하였으나 실제 추락 예측 시각 및 지점과 비교하면 정밀도가 매우 떨어져 실제 우주물체 추락 상황에서는 사용할 수 없는 문제점이 있다.
상기와 같은 문제점을 해결하기 위한 본 발명의 목적은 대기 보정 비율 추정 방법(DSFE)을 적용하여 추락 예상 시각 및 지점을 정밀하게 예측하는 통제 불가능한 인공우주물체 재진입 예측 방법을 제공하는데 있다.
상기와 같은 목적을 달성하기 위한 본 발명의 통제가 불가능한 인공우주물체 재진입 예측 방법은, 인공우주물체의 다른 두 시점의 접촉궤도요소 또는 평균궤도요소(TLE)를 입력하여 평균 장반경과 위도 인수를 계산하는 단계; 상기 첫 번째 시점의 궤도정보를 임의 상수인 초기 대기보정 비율 값을 적용하여 두 번째 시점까지 수치적분을 이용한 코웰 방법의 고정밀 궤도 전파기를 통해 궤도 전파를 수행하여 예측된 두 번째 시점에서의 평균 장반경, 위도인수 및 대기항력을 계산하는 단계; 상기 예측된 평균 장반경 또는 위도인수 값과 초기 입력된 두 번째 시점의 평균 장반경 또는 위도 인수 값을 비교하여 오차가 임의의 수렴 값보다 작아질 때까지 대기 보정 비율 값을 변경하면서 최적의 대기 보정 비율 값을 결정하는 단계; 및 상기 최적화된 대기 보정 비율 값을 적용하여 두 번째 시점부터 재진입 시점까지 수치적분을 이용한 코웰 방법의 고정밀 궤도 전파기를 통해 궤도 예측을 수행하여 인공우주물체의 대기권 재진입 시점 및 지점을 예측하는 단계;;를 포함하는 것을 특징으로 한다.
상기 평균궤도요소(TLE)는 접촉궤도요소로 변환하여 TOD(True of Date)좌표계에서의 평균궤도를 계산하여 사용할 수 있다.
상기 수렴 값은 사용자가 임의로 정하는 위치 오차일 수 있다.
상기와 같이, 본 발명에 따르면 대기 보정 비율 추정을 적용하여 통제가 불가능한 인공우주물체의 대기권 재진입 시점 및 지점을 정밀하게 예측할 수 있다.
도 1은 본 발명의 제1 실시예에 따른 통제 불가능한 인공우주물체 재진입 예측 방법을 나타낸 흐름도이다.
도 2는 본 발명의 제2 실시예에 따른 통제 불가능한 인공우주물체 재진입 예측 방법을 나타낸 흐름도이다.
도 3은 본 발명의 제3 실시예에 따른 통제 불가능한 인공우주물체 재진입 예측 방법을 나타낸 흐름도이다.
아래에서는 첨부한 도면을 참고로 하여 본 발명의 실시 예에 대하여 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자가 용이하게 실시할 수 있도록 상세히 설명한다. 그러나 본 발명은 여러 가지 상이한 형태로 구현될 수 있으며 여기에서 설명하는 실시 예에 한정되지 않는다.
이하에서는 본 발명의 실시예에 따른 통제 불가능한 인공우주물체 재진입 예측 방법에 대하여 설명한다.
도 1 내지 도 3은 본 발명의 제1 내지 제3 실시예에 따른 통제 불가능한 인공우주물체 재진입 예측 방법을 나타낸 흐름도이다.
도 1 내지 도 3을 참조하면, 본 발명에 따른 통제 불가능한 인공우주물체 재진입 예측 방법은, 먼저 인공우주물체의 다른 두 시점(
Figure PCTKR2018007053-appb-I000001
,
Figure PCTKR2018007053-appb-I000002
)에서 초기 궤도요소(
Figure PCTKR2018007053-appb-I000003
,
Figure PCTKR2018007053-appb-I000004
)를 입력하여 평균 장반경(
Figure PCTKR2018007053-appb-I000005
,
Figure PCTKR2018007053-appb-I000006
)과 위도 인수(
Figure PCTKR2018007053-appb-I000007
,
Figure PCTKR2018007053-appb-I000008
)를 계산한다(S100,S200,S300). 여기서, 상기 궤도요소는 접촉궤도요소 또는 평균궤도요소(TLE)이며, 평균궤도요소일 경우 접촉궤도요소로 변환하여 TOD(True of Date)좌표계에서의 평균궤도를 계산하여 사용할 수 있다.
다음으로, 첫 번째 시점(
Figure PCTKR2018007053-appb-I000009
)의 궤도정보를 임의 상수인 초기 대기 보정 비율 값(
Figure PCTKR2018007053-appb-I000010
)을 적용하여 두 번째 시점(
Figure PCTKR2018007053-appb-I000011
)까지 궤도 전파를 수행한다(S110,S210,S310). 이때의 궤도 전파는 수치적분을 이용한 코웰(Cowell) 방법의 고정밀 궤도 전파기를 통해 예측된 두 번째 시점(
Figure PCTKR2018007053-appb-I000012
)에서의 궤도요소(
Figure PCTKR2018007053-appb-I000013
)에 따른 평균 장반경(
Figure PCTKR2018007053-appb-I000014
) 및 위도인수(
Figure PCTKR2018007053-appb-I000015
)와 대기항력(
Figure PCTKR2018007053-appb-I000016
)을 계산한다. 여기서,
Figure PCTKR2018007053-appb-I000017
는 항력계수, A는 단면적, m은 질량, ρ는 기밀도,
Figure PCTKR2018007053-appb-I000018
는 속도벡터,
Figure PCTKR2018007053-appb-I000019
는 속도벡터 크기이다.
상기 코웰(Cowell) 방법의 고정밀 궤도 전파기는 지구 중력장, 대기 영향, 태양과 달의 인력, 태양 복사압 등 인공우주물체에 영향을 주는 모든 섭동력을 고려하여 임의의 시각에서의 인공우주물체의 위치와 속도를 얻는 알고리즘으로 이 분야에서 널리 알려진 기술이므로 자세한 설명은 생략하기로 한다.
다음에, 도 1과 같이 평균 장반경 비교 값, 도 2와 같이 위도인수 비교 값, 도 3과 같이 평균 장반경 비교 값 및 위도인수 비교 값 중 어느 하나의 비교 값의 오차가 수렴 값과 비교하여(S120,S220,S320) 최소가 되면 최적의 대기보정 비율 값(
Figure PCTKR2018007053-appb-I000020
)으로 결정하고(S140,S240,S340), 그렇지 않을 경우 대기 보정 비율 값을 변화시키며(S130,S230,S330), 해당 과정을 반복한다. 즉, 첫 번째 시점(
Figure PCTKR2018007053-appb-I000021
)에서 두 번째 시점(
Figure PCTKR2018007053-appb-I000022
)까지 대기보정 비율 값(
Figure PCTKR2018007053-appb-I000023
)을 반영하여 예측된 평균 장반경(
Figure PCTKR2018007053-appb-I000024
) 또는 위도인수(
Figure PCTKR2018007053-appb-I000025
) 값과 초기 입력된 두 번째 시점(
Figure PCTKR2018007053-appb-I000026
)의 평균 장반경(
Figure PCTKR2018007053-appb-I000027
) 또는 위도인수(
Figure PCTKR2018007053-appb-I000028
) 값을 비교하여, 그 오차가 수렴 값보다 작아질 때까지 대기 보정 비율 값을 변경하면서 최적의 대기 보정 비율 값(
Figure PCTKR2018007053-appb-I000029
)을 찾는다. 여기서, 상기 수렴 값은 사용자가 임의로 정하는 위치 오차로 예를 들면 10-4km 등으로 정하는 값이다.
다음에, 최적화된 대기보정 비율 값(
Figure PCTKR2018007053-appb-I000030
)을 적용하여, 두 번째 시점(
Figure PCTKR2018007053-appb-I000031
)부터 재진입 시점까지 수치적분을 이용한 코웰 방법의 고정밀 궤도 전파기를 통해 궤도 예측을 수행하여 고도 100km 이내의 재진입 시점 및 지점 예측 정밀도를 향상시켜 통제 불가능한 인공우주물체의 대기권 재진입 시점 및 지점(위도, 경도, 고도)을 예측한다(S150,S250,S350).
이상에서 본 발명의 실시예에 대하여 상세하게 설명하였지만 본 발명의 권리범위는 이에 한정되는 것은 아니고 다음의 청구범위에서 정의하고 있는 본 발명의 기본 개념을 이용한 당업자의 여러 변형 및 개량 형태 또한 본 발명의 권리범위에 속하는 것이다.

Claims (3)

  1. 인공우주물체의 다른 두 시점의 접촉궤도요소 또는 평균궤도요소(TLE)를 입력하여 평균 장반경과 위도 인수를 계산하는 단계;
    상기 첫 번째 시점의 궤도정보를 임의 상수인 초기 대기 보정 비율 값을 적용하여 두 번째 시점까지 수치적분을 이용한 코웰 방법의 고정밀 궤도 전파기를 통해 궤도 전파를 수행하여 예측된 두 번째 시점에서의 평균 장반경, 위도인수 및 대기항력을 계산하는 단계;
    상기 예측된 평균 장반경 또는 위도인수 값과 초기 입력된 두 번째 시점의 평균 장반경 또는 위도 인수 값을 비교하여 오차가 임의의 수렴 값보다 작아질 때까지 대기 보정 비율 값을 변경하면서 최적의 대기 보정 비율 값을 결정하는 단계; 및
    상기 최적화된 대기 보정 비율 값을 적용하여 두 번째 시점부터 재진입 시점까지 수치적분을 이용한 코웰 방법의 고정밀 궤도 전파기를 통해 궤도 예측을 수행하여 인공우주물체의 대기권 재진입 시점 및 지점을 예측하는 단계;를 포함하는 것을 특징으로 하는 통제 불가능한 인공우주물체 재진입 예측 방법.
  2. 제1항에 있어서,
    상기 평균궤도요소(TLE)는 접촉궤도요소로 변환하여 TOD(True of Date)좌표계에서의 평균궤도를 계산하여 사용하는 것을 특징으로 하는 통제 불가능한 인공우주물체 재진입 예측 방법.
  3. 제1항에 있어서,
    상기 수렴 값은 사용자가 임의로 정하는 위치 오차인 것을 특징으로 하는 통제 불가능한 인공우주물체 재진입 예측 방법.
PCT/KR2018/007053 2018-06-01 2018-06-22 통제 불가능한 인공우주물체 재진입 예측 방법 WO2019231039A1 (ko)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
KR1020180063148A KR101888339B1 (ko) 2018-06-01 2018-06-01 통제 불가능한 인공우주물체 재진입 예측 방법
KR10-2018-0063148 2018-06-10

Publications (1)

Publication Number Publication Date
WO2019231039A1 true WO2019231039A1 (ko) 2019-12-05

Family

ID=63719851

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/KR2018/007053 WO2019231039A1 (ko) 2018-06-01 2018-06-22 통제 불가능한 인공우주물체 재진입 예측 방법

Country Status (4)

Country Link
US (1) US20190377770A1 (ko)
EP (1) EP3578467A1 (ko)
KR (1) KR101888339B1 (ko)
WO (1) WO2019231039A1 (ko)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111695081A (zh) * 2020-06-12 2020-09-22 北京卫星环境工程研究所 空间目标获取方法、装置、设备和存储介质
CN114460952A (zh) * 2022-01-17 2022-05-10 上海卫星工程研究所 椭圆轨道伴飞构型初始化的双星协同变轨方法及系统

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110595485B (zh) * 2019-08-08 2023-08-11 上海航天电子有限公司 基于两行根数的低轨卫星长期轨道预报方法
CN114114359B (zh) * 2022-01-27 2022-04-26 中国人民解放军32035部队 单星与地基设备联合的再入预报方法、装置和电子设备
CN116384599B (zh) * 2023-06-06 2023-08-22 中国空气动力研究与发展中心超高速空气动力研究所 基于能量分析的航天器leo圆轨道衰降过程参数预报方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005096702A (ja) * 2003-09-26 2005-04-14 Fujitsu Ltd 人工衛星の軌道制御計画方法、プログラム及び装置
JP2009220622A (ja) * 2008-03-13 2009-10-01 Fujitsu Ltd 中高度人工衛星捕捉方法及び装置
JP2009227019A (ja) * 2008-03-21 2009-10-08 Fujitsu Ltd 準天頂衛星の制御計画生成プログラム、方法及び装置
JP2013167606A (ja) * 2012-02-17 2013-08-29 Nec Toshiba Space Systems Ltd 軌道位置推定方法、軌道位置推定装置及び軌道位置推定プログラム
JP2016043891A (ja) * 2014-08-26 2016-04-04 三菱電機株式会社 宇宙飛翔物体監視装置、宇宙飛翔物体監視方法、及び、プログラム

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3978952B2 (ja) * 1999-10-15 2007-09-19 三菱電機株式会社 人工衛星搭載用時刻補正装置及びコマンド装置
US7443340B2 (en) * 2001-06-06 2008-10-28 Global Locate, Inc. Method and apparatus for generating and distributing satellite tracking information
JP6273164B2 (ja) * 2014-05-09 2018-01-31 三菱重工業株式会社 観測支援装置、観測支援方法、及び観測支援プログラム
KR101877382B1 (ko) 2017-08-16 2018-07-12 한국항공우주연구원 Tle 보상기법을 적용한 tle 생성방법 및 이를 이용한 생성장치

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005096702A (ja) * 2003-09-26 2005-04-14 Fujitsu Ltd 人工衛星の軌道制御計画方法、プログラム及び装置
JP2009220622A (ja) * 2008-03-13 2009-10-01 Fujitsu Ltd 中高度人工衛星捕捉方法及び装置
JP2009227019A (ja) * 2008-03-21 2009-10-08 Fujitsu Ltd 準天頂衛星の制御計画生成プログラム、方法及び装置
JP2013167606A (ja) * 2012-02-17 2013-08-29 Nec Toshiba Space Systems Ltd 軌道位置推定方法、軌道位置推定装置及び軌道位置推定プログラム
JP2016043891A (ja) * 2014-08-26 2016-04-04 三菱電機株式会社 宇宙飛翔物体監視装置、宇宙飛翔物体監視方法、及び、プログラム

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111695081A (zh) * 2020-06-12 2020-09-22 北京卫星环境工程研究所 空间目标获取方法、装置、设备和存储介质
CN111695081B (zh) * 2020-06-12 2021-09-07 北京卫星环境工程研究所 空间目标获取方法、装置、设备和存储介质
CN114460952A (zh) * 2022-01-17 2022-05-10 上海卫星工程研究所 椭圆轨道伴飞构型初始化的双星协同变轨方法及系统
CN114460952B (zh) * 2022-01-17 2023-03-24 上海卫星工程研究所 椭圆轨道伴飞构型初始化的双星协同变轨方法及系统

Also Published As

Publication number Publication date
KR101888339B1 (ko) 2018-09-20
US20190377770A1 (en) 2019-12-12
EP3578467A1 (en) 2019-12-11

Similar Documents

Publication Publication Date Title
WO2019231039A1 (ko) 통제 불가능한 인공우주물체 재진입 예측 방법
CN111721289B (zh) 自动驾驶中车辆定位方法、装置、设备、存储介质及车辆
CN109556632B (zh) 一种基于卡尔曼滤波的ins/gnss/偏振/地磁组合导航对准方法
CN111102981B (zh) 一种基于ukf的高精度卫星相对导航方法
Li et al. A novel algorithm for marine strapdown gyrocompass based on digital filter
US7347090B1 (en) Methods and systems for calculating atmospheric vehicle air data
CN111680462B (zh) 基于空间目标在光学相平面位置变化的制导方法和系统
RU2454631C1 (ru) Способ автономной навигации и ориентации космических аппаратов на основе виртуальных измерений зенитных расстояний звезд
JP2016114603A (ja) Gpsドリフトの処理方法、およびそのシステム
US10809718B2 (en) Systems and methods for controlling aerial vehicles
JP7195588B2 (ja) 予測装置、予測方法、及び予測プログラム
CN112161632A (zh) 一种基于相对位置矢量测量的卫星编队初始定位算法
CN104864875B (zh) 一种基于非线性h∞滤波的航天器自主定位方法
CN113916217B (zh) 一种基于分区平流层大气折射模型的星光定位方法
JP2009220622A (ja) 中高度人工衛星捕捉方法及び装置
CN110209185B (zh) 一种利用星敏信息姿态重置的航天器姿态机动稳定方法
Guo et al. Analysis and design of an attitude calculation algorithm based on elman neural network for SINS
US20190033884A1 (en) Systems and methods for controlling aerial vehicles
Gao et al. SIMU/Triple star sensors integrated navigation method of HALE UAV based on atmospheric refraction correction
CN115420280A (zh) 一种航天器脉冲星/星光折射组合导航方法
CN114001730A (zh) 融合定位方法、装置、计算机设备和存储介质
Chen et al. On SINS/GPS integrated navigation filtering method aided by radial basis function neural network
Lu et al. Optimal scheme of star observation of missile-borne inertial navigation system/stellar refraction integrated navigation
CN116358606B (zh) 惯导系统的初始粗对准方法、装置、设备及介质
CN114675276B (zh) 一种空间目标雷达轨道改进方法及装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18920955

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 18920955

Country of ref document: EP

Kind code of ref document: A1