黄道傾斜角

黄道傾斜角

 黄道傾斜角とは簡単に言えば地軸の傾きのことだ。
 どの天体に対しても黄道傾斜角はあるが、単に「黄道傾斜角」と言う場合には"地球の"黄道傾斜角である ことが多い。


黄道傾斜角の計算方法

 ここでは黄道傾斜角の計算方法として知られている公式のうち2つを挙げる。

黄道傾斜角 = 23.45 + (8.26 - 46.845*T - 0.0059*T*T + 0.00181*T*T*T) / 3600 ただしTは1900年1月1日12:00からの経過世紀(小数点を含む数値)

黄道傾斜角 = 23.433333333 + (21.448 - 46.8150*T - 0.00059*T*T + 0.001813*T*T*T)/3600 ただしTは2000年1月1日12:00からの経過世紀(小数点を含む数値)


ソースコード

//
//黄道傾斜角の取得
//
//1900年1月1日12:00からの経過世紀Tを利用した公式
//23.45 + (8.26 - 46.845*T - 0.0059*T*T + 0.00181*T*T*T) / 3600
//を利用
//
double	GetObliquity1900(double dJD)
{
	double	ret;
	double	c;

	c = (dJD - 2415021) / (100 * 365.25);
	ret = 23.45 + (8.26 - 46.845*c - 0.0059*c*c + 0.00181*c*c*c) / 3600;

	return	ret;
}


//
//黄道傾斜角の取得
//
//2000年1月1日12:00からの経過世紀Tを利用した公式
//23.433333333 + (21.448 - 46.8150*T - 0.00059*T*T + 0.001813*T*T*T)/3600
//を利用
//
double	GetObliquity2000(double dJD)
{
	double	ret;
	double	c;

	c = (dJD - 2451545) / (100 * 365.25);
	ret = 23.433333333 + (21.448 - 46.8150*c - 0.00059*c*c + 0.001813*c*c*c)/3600;

	return	ret;
}


カテゴリー「VC++ TIPS」 のエントリー