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

moonLibration

月の献杯

ページ内をすべて折りたたむ

構文

angles = moonLibration(ephemerisTime)
angles = moonLibration(ephemerisTime,ephemerisModel)
angles = moonLibration(ephemerisTime,ephemerisModel,action)
[angles,rates] = moonLibration(___)

説明

月の秤動角を実装する

angles = moonLibration(ephemerisTime) は、ユリウス日で表された ephemerisTime の月の秤動角を実装します。

この関数は、NASA ジェット推進研究所が提供するチェビシェフ係数を使用します。

この機能を使用するには、アドオン エクスプローラーを使用してエフェメリス データをダウンロードする必要があります。詳細については、aeroDataPackageを参照してください。

angles = moonLibration(ephemerisTime,ephemerisModel) は ephemerisModel 係数を使用してこれらの値を実装します。

angles = moonLibration(ephemerisTime,ephemerisModel,action) は action を使用してエラー報告を決定します。

例

月の秤動角度と速度を実装する

[angles,rates] = moonLibration(___) は、前の構文の入力引数の任意の組み合わせを使用して、月の秤動角度と秤動率を実装します。

例

例

すべて折りたたむ

DE405 を使用して、1990 年 12 月 1 日の月の秤動角を実装します。juliandate 関数を使用して、入力されたユリウス日付値を計算します。

angles = moonLibration(juliandate(1990,12,1))
angles =
   1.0e+03 *
    0.0001    0.0004    1.8010
ライブ スクリプトを開く

天体暦 (DE421) を指定し、日付 (2000 年 1 月 1 日) にユリウス暦関数を使用して、月の秤動角度と秤動率の両方を計算します。

[angles,rates] = moonLibration([2451544.5 0.5],'421')
angles = 1×3
103 ×

   -0.0001    0.0004    2.5643


rates = 1×3

   -0.0001    0.0000    0.2301

入力引数

すべて折りたたむ

位置を計算するユリウス日。次のいずれかとして指定します。

  • スカラー - 固定のユリウス日を 1 つ指定します。

  • 2 要素ベクトル - ユリウス日を複数の部分で指定します。最初の要素は、補間エポックまたはその前の最新の真夜中である特定のエポックのユリウス日です。2 番目の要素は、最初の要素とエポックの間で経過した日の小数部分です。2 番目の要素は正の値である必要があります。最初の要素と 2 番目の要素の合計値は、最大ユリウス日を超えることはできません。

  • 列ベクトル - M 要素を持つ列ベクトルを指定します。ここで、M はユリウス日の数です。

  • M 行 2 列の行列 — 行列を指定します。ここで、M はユリウス日の数で、2 番目の列には経過日数 (ユリウス暦の日付/経過日のペア) が含まれます。

データ型: double

ジェット推進研究所によって定義された次のエフェメリスのいずれかとして指定されるエフェメリス係数:

  • '405' — 1998 年にリリースされました。この暦は、ユリウス暦の日付範囲 2305424.50 (1599 年 12 月 9 日) から 2525008.50 (2201 年 2 月 20 日) を考慮に入れています。

    この関数は、1998 年に採用された国際天体基準フレーム バージョン 1.0 に基づいてこれらの天体暦を計算します。

  • '421' — 2008 年にリリースされました。この暦は、ユリウス暦の日付範囲 2414992.5 (1899 年 12 月 4 日) から 2469808.5 (2050 年 1 月 2 日) を考慮に入れています。

    この関数は、1998 年に採用された国際天体基準フレーム バージョン 1.0 に基づいてこれらの天体暦を計算します。

  • '423' — 2010 年にリリースされました。この暦は、ユリウス暦の日付範囲 2378480.5 (1799 年 12 月 16 日) から 2524624.5 (2200 年 2 月 1 日) を考慮に入れています。

    この関数は、2010 年に採用された国際天体基準フレーム バージョン 2.0 に基づいてこれらの天体暦を計算します。

  • '430' — 2013 年にリリースされました。この暦は、ユリウス暦の日付範囲 2287184.5 (1549 年 12 月 21 日) から 2688976.5 (2650 年 1 月 25 日) を考慮に入れています。

    この関数は、2010 年に採用された国際天体基準フレーム バージョン 2.0 に関するこれらの天体暦を実装します。

  • '432t' — 2014 年 4 月にリリースされました。この暦は、ユリウス暦の日付範囲 2287184.5 (1549 年 12 月 21 日) から 2688976.5 (2650 年 1 月 25 日) を考慮に入れています。

    この関数は、2010 年に採用された国際天体基準フレーム バージョン 2.0 に関するこれらの天体暦を実装します。

データ型: double

入力が範囲外の場合の関数の動作。次のいずれかの値として指定されます。

値説明
'None'何も起こりません。
'Warning'MATLAB® コマンド ウィンドウに警告が表示され、モデル シミュレーションが続行されます。
'Error'MATLAB は例外を返し、モデルのシミュレーションが停止します。

データ型: char | string

出力引数

すべて折りたたむ

月の秤動角は、M 行 3 列の数値配列として返されます。M は、行数で表したユリウス日の数です。列には、月の姿勢のオイラー角 (φ θ ψ) がラジアンで含まれています。

入力引数に複数のユリウス日またはエポックが含まれている場合、この配列の行数は ephemerisTime 入力と同じになります。

月の秤動角速度は、M 行 3 列の数値配列として返されます。M は、行数で表したユリウス日の数です。列には、月の秤動オイラー角速度 (ω) がラジアン/日で含まれています。

入力引数に複数のユリウス日またはエポックが含まれている場合、この配列の行数は ephemerisTime 入力と同じになります。

バージョン履歴

R2013a で導入

参考

| |

外部の Web サイト