このページは機械翻訳を使用して翻訳されました。最新版の英語を参照するには、ここをクリックします。

ecef2eci

地球中心の慣性平均赤道平均春分点における位置と速度ベクトル

構文

[r_eci] = ecef2eci(utc,r_ecef)

[r_eci,v_eci] = ecef2eci(___,v_ecef)

[r_eci,v_eci,a_eci] = ecef2eci(___,a_ecef)

[r_eci,v_eci,a_eci] = ecef2eci(___,Name,Value)

説明

[r_eci] = ecef2eci(utc,r_ecef) は、特定の協定世界時(UTC) における地球中心地球中心の地球固定(ECI) 座標系の位置ベクトルを計算します。地球中心の地球固定座標系の詳細については、「アルゴリズム」を参照してください。

[r_eci,v_eci] = ecef2eci(___,v_ecef) は、指定された位置ベクトルと速度ベクトルを計算します。

[r_eci,v_eci,a_eci] = ecef2eci(___,a_ecef) は、指定された位置、速度、加速度ベクトルの位置、速度、加速度ベクトルを計算します。

[r_eci,v_eci,a_eci] = ecef2eci(___,Name,Value) は、地球姿勢パラメーターを使用して、位置、速度、加速度のベクトルをより高い精度で計算します。地球姿勢パラメーターが指定されていない場合、関数はそれを 0 に設定します。

例

ECEFの位置と速度をECIに変換する

この例では、2019 年 1 月 4 日の 12:00 における ECEF の位置と速度を ECI に変換する方法を示します。

r_ecef = [-5762640 -1682738 3156028];
v_ecef = [3832 -4024 4837];
utc = [2019 1 4 12 0 0];
[r_eci, v_eci] = ecef2eci(utc, r_ecef, v_ecef);

極運動の影響を考慮したECEF位置をECIに変換する

この例では、極運動の影響を含めて、2019 年 1 月 4 日 12:00 の ECEF の位置を ECI に変換する方法を示します。

r_ecef = [-5762640 -1682738 3156028];
utc = [2019 1 4 12 0 0];
mjd = mjuliandate(utc);
pm = polarMotion(mjd, 'action', 'none')*180/pi;
r_eci = ecef2eci(utc, r_ecef, 'pm', pm);

datetime配列を使用してECEF位置をECIに変換する

この例では、datetime 配列 utcDT を使用して、2019 年 1 月 4 日 12:00 の ECEF の位置と速度を ECI に変換する方法を示します。

r_ecef = [-5762640 -1682738 3156028];
v_ecef = [3832 -4024 4837];
utcDT = datetime(2019, 1, 4, 12, 0, 0)

utcDT = datetime
   04-Jan-2019 12:00:00

[r_eci, v_eci] = ecef2eci(utcDT, r_ecef, v_ecef)

r_eci = 3×1
106 ×

   -2.9818
    5.2070
    3.1616

v_eci = 3×1
103 ×

   -3.3837
   -4.8870
    4.8430

入力引数

utc — 協定世界時

1行6列の配列 | 1行6列の行列 | スカラー datetime 配列

協定世界時(UTC) は次のいずれかとして指定されます。

• 年、月、日、時、分、秒の順序で並べられた UTC 値の 1 行 6 列の配列:

  時間値	入力
  年	2013 など、1 より大きい整数である double 値。
  月	1 から 12 の範囲内の、0 より大きい整数である double 値。
  日	1 から 31 の範囲内の、0 より大きい整数である double 値。
  時	1 から 24 の範囲内の、0 より大きい整数である double 値。
  分と秒	1 から 60 の範囲内の、0 より大きい整数である double 値。

• スカラー datetime 配列。配列を作成するには、datetime 関数を使用します。

例: [2000 1 12 4 52 12.4]

データ型: double

r_ecef — 位置コンポーネント

3行1列の配列

ECEF 位置コンポーネントの配列。3 行 1 列の配列として指定されます。

データ型: double

v_ecef — 速度成分

3行1列の配列

ECEF 速度成分。3 行 1 列の配列として指定されます。

データ型: double

a_ecef — 加速コンポーネント

3行1列の配列

ECEF 加速コンポーネント。3 行 1 列の配列として指定されます。

データ型: double

名前と値の引数

オプションの引数のペアを Name1=Value1,...,NameN=ValueN として指定します。ここで、Name は引数名で、Value は対応する値です。名前と値の引数は他の引数の後に指定しなければなりませんが、ペアの順序は重要ではありません。

R2021a より前では、コンマを使用して名前と値をそれぞれ区切り、Name を引用符で囲みます。

例: 'dUT1',0.234

dAT — TAIとUTCの違い

0 (既定値) | スカラー

国際原子時(TAI) と UTC の差をスカラー値 (秒) で指定します。

例: 32

データ型: double

dUT1 — UTCと世界時の違い

0 (既定値) | スカラー

UTC と世界時(UT1) の差をスカラー値 (秒単位) で指定します。

例: 0.234

データ型: double

pm — 極移動

ゼロの配列 (既定値) | 1行2列の配列

x 軸と y 軸に沿った地殻の動きによる極移動(度)。

ヒント

変位を計算するには、polarMotion 関数を使用します。

例: pm = polarMotion(mjd, 'action', 'none')*180/pi;

データ型: double

dCIP — CIP位置の調整

1行2列の配列

天の中間極 (CIP) の位置の調整 (度単位)。dCIP と M 行 2 列の配列で構成されるコンマ区切りのペアとして指定します。この位置 (dDeltaX、dDeltaY) は、x 軸と y 軸に沿っています。デフォルトでは、この関数は 1 行 2 列のゼロ配列を想定します。

履歴値については、国際地球回転・基準系サービス Web サイト (https://www.iers.org) を参照し、地球姿勢データのデータ/製品ページに移動してください。

• M 行 2 列の配列

  位置調整値の M 行 2 列の配列を指定します。ここで、M は変換する方向余弦または変換行列の数です。各行は、dDeltaX 値と dDeltaY 値の 1 セットに対応します。

例: [-0.2530e-6 -0.0188e-6]

データ型: double

lod — 暦日の長さが長すぎる

0 (既定値) | スカラー

暦日の長さ（天文学的に決定された一日の長さと 86400 SI 秒の差）。スカラー値として秒単位で指定します。

例: 32

データ型: double

出力引数

r_eci — 位置コンポーネント

3行1列の配列

ECI 位置コンポーネント。3 行 1 列の配列として指定されます。

v_eci — 速度成分

3行1列の配列

ECI 速度成分。3 行 1 列の配列として指定されます。

a_eci — 加速コンポーネント

3行1列の配列

ECI 加速コンポーネント。3 行 1 列の配列として指定されます。

アルゴリズム

ecef2eci 関数は、次の地球中心の座標系を使用します。

• 地球中心慣性系 (ECI) — 使用される慣性系は国際天体基準座標系 (ICRF) です。このフレームは、J2000 (2000 年 1 月 1 日 12:00:00 TT) で実現された ECI 座標系と同等に扱うことができます。詳細については、ECI and ECEF Coordinatesを参照してください。

• 地球中心地球固定座標系 (ECEF) — 使用される固定座標系は、国際地球基準座標系 (ITRF) です。この基準座標系は、ICRF 座標系からの IAU2000/2006 縮小によって実現されます。詳細については、ECEF Coordinatesを参照してください。

参照

[1] Vallado, D. A. Fundamentals of Astrodynamics and Applications. alg. 4. New York: McGraw-Hill, 1997.

[2] Gottlieb, R. G., "Fast Gravity, Gravity Partials, Normalized Gravity, Gravity Gradient Torque and Magnetic Field: Derivation, Code and Data," Technical Report NASA Contractor Report 188243, NASA Lyndon B. Johnson Space Center, Houston, Texas, February 1993.

[3] Konopliv, A. S., S. W. Asmar, E. Carranza, W. L. Sjogen, D. N. Yuan., "Recent Gravity Models as a Result of the Lunar Prospector Mission, Icarus", Vol. 150, no. 1, pp 1–18, 2001.

[4] Lemoine, F. G., D. E. Smith, D.D. Rowlands, M.T. Zuber, G. A. Neumann, and D. S. Chinn, "An improved solution of the gravity field of Mars (GMM-2B) from Mars Global Surveyor", Journal Of Geophysical Research, Vol. 106, No. E10, pp 23359-23376, October 25, 2001.

[5] Seidelmann, P.K., Archinal, B.A., A’hearn, M.F. et al. "Report of the IAU/IAG Working Group on cartographic coordinates and rotational elements: 2006." Celestial Mech Dyn Astr 98, 155–180 (2007).

バージョン履歴

R2019a で導入

参考

eci2ecef | dcmeci2ecef | aeroReadIERSData | deltaCIP | polarMotion | deltaUT1 | siderealTime | datetime | CubeSat Vehicle (Aerospace Blockset)