プロが教える店舗＆オフィスのセキュリティ対策術

北斗七星の位置と角度

解決済

  • 質問者：2hosi2
  • 質問日時：
  • 回答数：4件

日時から北斗七星の位置と角度を計算で求める方法はありますか？
わかりやすくお願いしますｍ（＿＿）ｍ

通報する

この質問への回答は締め切られました。

質問の本文を隠す

A 回答 (4件)

No.4ベストアンサー

  • 回答者： equinox2
  • 回答日時：

JavaScriptで書くとこんな感じです。

　北極星と北斗七星の位置計算（東京）
------- ここから　------
<HTML>
<HEAD>
<TITLE>北斗七星の位置</TITLE>
<script language="JavaScript">
<!--

function Calcmjd(yy,mm,dd)
{ // 準ユリウス日(mjd)の計算
if( mm <= 2) {
y = yy - 1; m = mm + 12;
} else {
y = yy; m = mm;
}
ret = Math.floor(365.25 * y) + Math.floor(y / 400) - Math.floor(y / 100) ;
ret = ret + Math.floor(30.59 * ( m - 2)) + dd - 678912;
return ret;
}
// -->
</script>
</HEAD>

<BODY BGCOLOR="#E0E0EE">

<script language="JavaScript">
<!--
//
RA = new Array(8); // 赤経データ
DE = new Array(8); // 赤緯データ
NM = new Array(8); // 星名データ
PI = 3.14159265358979;
RAD = 180.0 / PI;
KEIDO = 139.75; // 経度　東経139度45分(東京)
IDO = 35.65; // 緯度　北緯35度39分(東京)

//
// 日時関連データの設定
//
CuDate = new Date(); // 現在の日時
yy = CuDate.getYear();
mm = CuDate.getMonth() + 1;
dd = CuDate.getDate();
hr = CuDate.getHours();
mn = CuDate.getMinutes();
sc = CuDate.getSeconds();

MJD = Calcmjd(yy,mm,dd) + hr / 24 + mn / 1440 + sc / 86400 - 0.375;
d = (0.671262 + 1.002737909 * (MJD - 40000) + KEIDO/360);
LST = 2*PI*(d - Math.floor(d));
document.write("<TT><CENTER><B>東京から見た北斗七星の位置<P>");
document.write(yy, "年", mm, "月", dd,"日 ", hr, "時", mn,"分",sc,"秒");
document.write("<P>　　　　　方位角　　 高度")
// 星のデータの設定
RA[0] = 2.5300; DE[0] = 89.267; NM[0] = "ＵＭｉ　α　";
RA[1] = 11.062; DE[1] = 61.750; NM[1] = "ＵＭａ　α　";
RA[2] = 11.030; DE[2] = 56.383; NM[2] = "ＵＭａ　β　";
RA[3] = 11.897; DE[3] = 53.700; NM[3] = "ＵＭａ　γ　";
RA[4] = 12.257; DE[4] = 57.033; NM[4] = "ＵＭａ　δ　";
RA[5] = 12.900; DE[5] = 55.967; NM[5] = "ＵＭａ　ε　";
RA[6] = 13.398; DE[6] = 54.933; NM[6] = "ＵＭａ　ζ　";
RA[7] = 13.792; DE[7] = 49.317; NM[7] = "ＵＭａ　η　";
// 星の表示位置の計算
document.write("<BR>");
srid = Math.sin(IDO / RAD);
crid = Math.cos(IDO / RAD);
for (i=0;i<8;i++) {
ra = 15 * RA[i] / RAD;
dc = DE[i] / RAD;
ha = LST - ra;
sdc = Math.sin(dc);
cdc = Math.cos(dc);
sha = Math.sin(ha);
cha = Math.cos(ha);
xs = sdc * srid + cdc * crid * cha;
h = Math.asin(xs);
s = cdc * sha;
c = cdc * srid * cha - sdc * crid;
if (c < 0) a = Math.atan(s/c) + PI;
else if (c > 0 && s <= 0) a = Math.atan(s/c) + 2*PI;
else a = Math.atan(s/c);
if (h == 0) h = 0.00001;
a = a * RAD;
h = h * RAD;
if (h < 0) continue;

if (a < 10) sa = "00" + a;
else if (a < 100) sa = "0" + a;
else sa = "" + a;

if (h < 10) sh = "0" + h;
else sh = "" + h;

if (i == 1) document.write("<BR>");
ssa = sa.substring(0,6);
ssh = sh.substring(0,5);
document.write(NM[i]," ",ssa,"　　",ssh,"<BR>");
}
// -->
</script>
</BODY>
</HTML>
    • good
    • 0
    • 件
通報する
この回答へのお礼

模範プログラムありがとうございました。
とても参考になりました＾＾

通報する
お礼日時：2009/02/24 20:38

No.3

  • 回答者： equinox2
  • 回答日時：

恒星の位置計算はこんな方法です。


http://star.gs/nyumon/keisan.htm

・北斗七星の各星の赤経・赤緯を調べる。
・観測地の緯度・経度を調べる。
・上記の計算方法で計算する。
    • good
    • 0
    • 件
通報する

No.2

  • 回答者： adasw-12
  • 回答日時：

私は８月の９時の位置を覚えているので暗算でしてしまうので


計算した事がありません。
http://www.nao.ac.jp/hoshizora/image/chart08_s.jpg

北斗七星は真西です。

９月には１５度斜めに（西に）傾き、
http://www.ksky.ne.jp/~tatsuo/gif/9gatu.gif

http://www.nao.ac.jp/hoshizora/chartlist.html

１１月に真下へ。
２月には真東へ。
５月には真上へ。
１ヶ月で１５度。３ヶ月で９０度を加えています。
（実際には３６０/恒星年）

これを計算すれば、
２月１５日の８時を基準点に取り、
観測月日（経過日数）を度に置き換えます。
（計算結果１）

時刻８時を基準点に取り、７時なら－１５度、９時なら＋１５度を加えます。

結果が９０度なら、真南。
１８０度なら真西。
２７０度なら真下。
となります。

真北を０度にしたい場合、
１１月１５日８時を０度基準にすれば良いです。

とりあえずは「雑な計算をし」合っているか確かめるべきです。
計算が合わないなら、その理由を考える事です。

太陽年と恒星年は別物です。
恒星は歳差の影響を受けます。

私も意地悪ですので、
たまに方程式にはに故意的な間違いを記載しておく事があります。
使えばすぐに気が付く物ですから。＾○＾

使わなければ永久に知る事はありません。
頑張って勉強してください。＾＾
    • good
    • 0
    • 件
通報する
この回答へのお礼

詳しい解説ありがとうござます＾o＾
気付いた点は一ヶ月30度ってところだけです・・・。

ありがとうございました。

通報する
お礼日時：2009/02/24 20:35

No.1

  • 回答者： houjutucho
  • 回答日時：

天文年間・天文ガイド等の雑誌を読むか


星座早見盤を、ご利用して下さい。
    • good
    • 0
    • 件
通報する

お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう！

Q質問する（無料）

関連するカテゴリからQ&Aを探す

ページトップページトップ

Q質問する（無料）

デイリーランキングこのカテゴリの人気デイリーQ&Aランキング

マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング

おすすめ情報