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

dcmeci2ecef

地球中心慣性座標（ECI）を地球中心の地球固定座標（ECEF）に変換する

構文

dcm=dcmeci2ecef(reduction,utc)

dcm = dcmeci2ecef(reduction,utc,deltaAT)

dcm = dcmeci2ecef(reduction,utc,deltaAT,deltaUT1)

dcm = dcmeci2ecef(reduction,utc,deltaAT,deltaUT1,polarmotion)

dcm = dcmeci2ecef(reduction,utc,deltaAT,deltaUT1,polarmotion,Name,Value)

説明

dcm=dcmeci2ecef(reduction,utc) は、位置方向余弦行列(ECI から ECEF) を 3 x 3 x M 配列として計算します。計算は、指定された削減方法と協定世界時(UTC) に基づいて行われます。

dcm = dcmeci2ecef(reduction,utc,deltaAT) は、国際原子時と UTC の差を使用して、位置方向余弦行列を計算します。

dcm = dcmeci2ecef(reduction,utc,deltaAT,deltaUT1) は、UTC と世界時(UT1) の差を使用します。

dcm = dcmeci2ecef(reduction,utc,deltaAT,deltaUT1,polarmotion) は極移動を使用します。

dcm = dcmeci2ecef(reduction,utc,deltaAT,deltaUT1,polarmotion,Name,Value) は、1 つ以上の名前と値のペアの引数によって指定された追加オプションを使用します。

例

IAU-2000/2006縮小法を使用して変換する

2000 年 1 月 12 日、4 時間 52 分 12.4 秒と 2000 年 1 月 12 日、4 時間 52 分 13 秒の地球中心慣性 (ECI) 座標を地球中心の地球固定(ECEF) 座標に変換します。IAU-2000/2006 削減を使用します。削減方法とUTCのみを指定します。

dcm = dcmeci2ecef('IAU-2000/2006',[2000 1 12 4 52 12.4;2000 1 12 4 52 13])

dcm = 
dcm(:,:,1) =

   -0.9975   -0.0708   -0.0000
    0.0708   -0.9975   -0.0000
   -0.0000   -0.0000    1.0000


dcm(:,:,2) =

   -0.9975   -0.0709   -0.0000
    0.0709   -0.9975   -0.0000
   -0.0000   -0.0000    1.0000

IAU-76/FK5縮約法を使用して変換する

この例では、2000 年 1 月 12 日 4 時 52 分 12.4 秒の地球中心慣性 (ECI) 座標を地球中心の地球固定(ECEF) 座標に変換する方法を示します。IAU-76/FK5 削減を使用します。極運動などのオプションの引数も含め、すべての引数を指定します。

dcm = dcmeci2ecef('IAU-76/FK5',[2000 1 12 4 52 12.4],32,0.234,[-0.0682e-5 ...
0.1616e-5],'dNutation',[-0.2530e-6 -0.0188e-6])

dcm = 3×3

   -0.9975   -0.0708   -0.0000
    0.0708   -0.9975   -0.0000
   -0.0000   -0.0000    1.0000

IAU-76/FK5縮約とdatetime配列を使用した変換

この例では、日時配列 utcDT として表される 2 つの UTC 日付に対して、IAU-2000/2006 縮小を使用して、地球中心慣性 (ECI) 座標を地球中心の地球固定(ECEF) 座標に変換する方法を示します。その他のすべてのパラメーターはデフォルトで null 配列になります。

utcDT = datetime([2000 1 12 4 52 12.4;2000 1 12 4 52 13])

utcDT = 2×1 datetime
   12-Jan-2000 04:52:12
   12-Jan-2000 04:52:13

DCM = dcmeci2ecef('IAU-2000/2006', utcDT)

DCM = 
DCM(:,:,1) =

   -0.9975   -0.0708   -0.0000
    0.0708   -0.9975   -0.0000
   -0.0000   -0.0000    1.0000


DCM(:,:,2) =

   -0.9975   -0.0709   -0.0000
    0.0709   -0.9975   -0.0000
   -0.0000   -0.0000    1.0000

入力引数

reduction — 削減方法

'IAU-76/FK5' | 'IAU-2000/2006'

方向余弦行列を計算するための削減方法。次のいずれかとして指定します。

• IAU-76/FK5

  国際天文学連合 (IAU)-76/第 5 基本カタログ (FK5) (IAU-76/FK5) 参照システムを使用して計算を削減します。変換の参照座標系が FK5 の場合は、この削減方法を選択します。この削減では、'dNutation' の名前と値のペアを使用できます。

  メモ

  この方法では、IAU 1976 の歳差運動モデルと IAU 1980 の章動理論を使用して計算を削減します。このモデルと理論はもはや最新ではありませんが、ソフトウェアは既存の実装に対してこの削減方法を提供します。この削減法が使用する極運動近似のため、dcmeci2ecef は方向余弦行列ではなく変換行列を計算します。

• IAU-2000/2006

  国際天文学連合 (IAU)-2000/2005 参照システムを使用して計算を削減します。変換の基準座標系が IAU-2000 の場合は、この削減方法を選択します。この削減方法では、P03 歳差運動モデルを使用して計算を削減します。この削減では、'dCIP' の名前と値のペアを使用できます。

utc — 協定世界時

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

関数が方向余弦行列を計算するときの基準となる、年、月、日、時、分、秒の順序の協定世界時(UTC)。次のいずれかとして指定します。

• 年の値には、2013 など、1 より大きい整数の double 値を入力します。

• 月の値には、1 から 12 の範囲内の 0 より大きい整数の double 値を入力します。

• 日の値には、1 から 31 の範囲内の 0 より大きい整数の double 値を入力します。

• 時間の値には、1 から 24 の範囲内の 0 より大きい整数の double 値を入力します。

• 分と秒の値には、1 から 60 の範囲内の 0 より大きい整数の double 値を入力します。

これらの値は、次のいずれかの形式で指定します。

• 1行1列の配列

  1 行 6 列の datetime 配列を指定します。

• M 行 6 列の行列

  UTC 値の M 行 6 列の配列を指定します。ここで、M は計算する方向余弦または変換行列の数です。各行は 1 セットの UTC 値に対応します。

• 1行1列の配列

  datatime 配列の 1 行 1 列の配列を指定します。配列を作成するには、datetime 関数を使用します。

• M行1列の配列

  M 変換行列の M 行 1 列の datatime 配列を、UTC 日付ごとに 1 つ指定します。配列を作成するには、datetime 関数を使用します。

例: [2000 1 12 4 52 12.4] は、UTC 値の 1 行 6 列の配列です。

例: [2000 1 12 4 52 12.4;2010 6 5 7 22 0] は UTC 値の M 行 6 列の配列で、M は 2 です。

データ型: double

deltaAT — 国際原子時とUTCの違い

スカラー | 1 次元配列

関数が方向余弦または変換行列を計算する、国際原子時(TAI) と UTC の差 (秒単位)。デフォルトでは、関数は M 行 1 列のゼロの配列を想定します。

• スカラー

  1 つの方向余弦または変換行列を計算するには、1 つの差分時間値を指定します。

• 1 次元配列

  M 要素を持つ 1 次元配列を指定します。ここで、M は計算する方向余弦または変換行列の数です。各行は 1 セットの UTC 値に対応します。

例: 32

TAI と UTC の差として 32 秒を指定します。

データ型: double

deltaUT1 — UTCと世界時（UT1）の違い

スカラー | 1 次元配列

関数が方向余弦または変換行列を計算する秒単位の UTC と世界時(UT1) の差。デフォルトでは、この関数は M 行 1 列のゼロの配列を想定します。

• スカラー

  1 つの方向余弦または変換行列を計算するには、1 つの差分時間値を指定します。

• 1 次元配列

  M 要素の差分時間値を持つ 1 次元配列を指定します。ここで、M は計算する方向余弦または変換行列の数です。各行は 1 セットの UTC 値に対応します。

例: 0.234

UTC と UT1 の差として 0.234 秒を指定します。

データ型: double

polarmotion — 極移動

1行2列の配列 | M 行 2 列の配列

x 軸と y 軸に沿った地殻の動きによる地球の極移動(ラジアン)。デフォルトでは、この関数は M 行 2 列のゼロ配列を想定します。

• 1行2列の配列

  1 行 2 列の極移動値の配列を指定して、1 つの方向余弦または変換行列を変換します。

• M 行 2 列の配列

  極移動値の M 行 2 列の配列を指定します。ここで、M は変換する方向余弦または変換行列の数です。各行は 1 セットの UTC 値に対応します。

例: [-0.0682e-5 0.1616e-5]

データ型: double

名前と値の引数

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

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

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

dNutation — 経度（dDeltaPsi）と黄道傾斜角（dDeltaEpsilon）の調整

M 行 2 列の配列

経度 (dDeltaPsi) と黄道傾斜角(dDeltaEpsilon) の調整 (ラジアン単位) は、dNutation と M 行 2 列の配列で構成されるコンマ区切りのペアです。この名前と値のペアを IAU-76/FK5 削減で使用します。デフォルトでは、この関数は M 行 2 列のゼロ配列を想定します。

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

• M 行 2 列の配列

  調整値の M 行 2 列の配列を指定します。ここで、M は変換する方向余弦または変換行列の数です。各行は、経度と黄道傾斜角の値の 1 セットに対応します。

データ型: double

dCIP — 天体中間極（CIP）の位置の調整

M 行 2 列の配列

天体中間極 (CIP) の位置の調整 (ラジアン単位)。dCIP と M 行 2 列の配列で構成されるコンマ区切りのペアとして指定します。この位置 (dDeltaX、dDeltaY) は、x 軸と y 軸に沿っています。この議論を IAU-200/2006 の削減と組み合わせて使用します。デフォルトでは、この関数は M 行 2 列のゼロの配列を想定します。

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

• M 行 2 列の配列

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

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

データ型: double

出力引数

dcm — 方向余弦または変換行列

3行3列のM配列

方向余弦または変換行列。3 行 3 列の M 配列として返されます。

バージョン履歴

R2013b で導入

参考

ecef2lla | ecef2eci | eci2ecef | geoc2geod | geod2geoc | lla2ecef | datetime