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

earthNutation

地球の章動運動を実装する

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

構文

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

説明

地球の章動角を実装する

angles = earthNutation(ephemerisTime) は、ユリウス日で表された ephemerisTime の国際天文学連合 (IAU) 1980 年章動シリーズを実装します。この関数は角度を返します。

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

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

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

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

例

地球の章動角度と速度を実装する

[angles,rates] = earthNutation(___) は、前の構文の入力引数の任意の組み合わせを使用して、国際天文学連合 (IAU) 1980 章動シリーズを実装します。この関数は角度と角速度を返します。

例

例

すべて折りたたむ

ライブ スクリプトを開く

1990 年 12 月 1 日の地球章動角度を実装します。エフェメリス モデルが指定されていないため、デフォルトの DE405 が使用されます。ユリウス日を指定するには、juliandate 関数を使用します。

angles = earthNutation(juliandate(1990,12,1))
angles = 1×2
10-4 ×

    0.6448    0.2083
ライブ スクリプトを開く

DE421 を使用して、2000 年 1 月 1 日の正午の地球の章動角度と角速度を実装します。

[angles,rates] = earthNutation([2451544.5 0.5],'421')
angles = 1×2
10-4 ×

   -0.6750   -0.2799


rates = 1×2
10-7 ×

    0.3687   -0.9937

入力引数

すべて折りたたむ

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

  • スカラー - 固定のユリウス日を 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 に関するこれらの天体暦を実装します。

データ型: double

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

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

データ型: char | string

出力引数

すべて折りたたむ

地球の章動角。M 行 2 列のベクトルとして返されます。ここで、M はユリウス日の数です。2 番目のベクトルには、ラジアン単位の d(psi) 角度と d(epsilon) 角度が含まれます。入力引数には複数のユリウス日またはエポックが含まれます。ベクトルの行数は、ephemerisTime 入力と同じです。

地球の章動角速度。M 行 2 列のベクトルとして返されます。ここで、M はユリウス日の数です。2 番目のベクトルには、ラジアン/日単位の d(psi) および d(epsilon) 角速度が含まれます。入力引数には複数のユリウス日またはエポックが含まれます。ベクトルの行数は、ephemerisTime 入力と同じです。

参照

[1] Folkner, W. M., J. G. Williams, and D. H. Boggs. "The Planetary and Lunar Ephemeris DE 421." JPL Interplanetary Network Progress Report 24-178, 2009.

[2] Vallado, David A. Fundamentals of Astrodynamics and Applications. McGraw-Hill, 1997.

バージョン履歴

R2013a で導入

参考

| |

外部の Web サイト