1, 0, 2,-2, 1, 5800, 10, 2, -3045, 0, -1,
2, 0, 2,-2, 2, 6443, 0, -7, -2768, 0, -4,
-2, 0, 0, 2, 1, -5774, -11, -15, 3041, 0, -5,
2, 0, 2, 0, 1, -5350, 0, 21, 2695, 0, 12,
0,-1, 2,-2, 1, -4752, -11, -3, 2719, 0, -3,
0, 0, 0,-2, 1, -4940, -11, -21, 2720, 0, -9,
-1,-1, 0, 2, 0, 7350, 0, -8, -51, 0, 4,
2, 0, 0,-2, 1, 4065, 0, 6, -2206, 0, 1,
1, 0, 0, 2, 0, 6579, 0, -24, -199, 0, 2,
0, 1, 2,-2, 1, 3579, 0, 5, -1900, 0, 1,
1,-1, 0, 0, 0, 4725, 0, -6, -41, 0, 3,
-2, 0, 2, 0, 2, -3075, 0, -2, 1313, 0, -1,
3, 0, 2, 0, 2, -2904, 0, 15, 1233, 0, 7,
0,-1, 0, 2, 0, 4348, 0, -10, -81, 0, 2,
1,-1, 2, 0, 2, -2878, 0, 8, 1232, 0, 4,
0, 0, 0, 1, 0, -4230, 0, 5, -20, 0, -2,
-1,-1, 2, 2, 2, -2819, 0, 7, 1207, 0, 3,
-1, 0, 2, 0, 0, -4056, 0, 5, 40, 0, -2,
0,-1, 2, 2, 2, -2647, 0, 11, 1129, 0, 5,
-2, 0, 0, 0, 1, -2294, 0, -10, 1266, 0, -4,
1, 1, 2, 0, 2, 2481, 0, -7, -1062, 0, -3,
2, 0, 0, 0, 1, 2179, 0, -2, -1129, 0, -2,
-1, 1, 0, 1, 0, 3276, 0, 1, -9, 0, 0,
1, 1, 0, 0, 0, -3389, 0, 5, 35, 0, -2,
1, 0, 2, 0, 0, 3339, 0, -13, -107, 0, 1,
-1, 0, 2,-2, 1, -1987, 0, -6, 1073, 0, -2,
1, 0, 0, 0, 2, -1981, 0, 0, 854, 0, 0,
-1, 0, 0, 1, 0, 4026, 0, -353, -553, 0, -139,
0, 0, 2, 1, 2, 1660, 0, -5, -710, 0, -2,
-1, 0, 2, 4, 2, -1521, 0, 9, 647, 0, 4,
-1, 1, 0, 1, 1, 1314, 0, 0, -700, 0, 0,
0,-2, 2,-2, 1, -1283, 0, 0, 672, 0, 0,
1, 0, 2, 2, 1, -1331, 0, 8, 663, 0, 4,
-2, 0, 2, 2, 2, 1383, 0, -2, -594, 0, -2,
-1, 0, 0, 0, 2, 1405, 0, 4, -610, 0, 2,
1, 1, 2,-2, 2, 1290, 0, 0, -556, 0, 0
);
var dL=0,dE=0,c;
for(i=0;i<77*11;i+=11){ //周期项取和计算
c = (nu[i]*l + nu[i+1]*l1 + nu[i+2]*F +nu[i+3]*D +nu[i+4]*Om)/rad;
dL += (nu[i+5]+nu[i+6]*t)*Math.sin(c) + nu[i+7] *Math.cos(c);
dE += (nu[i+8]+nu[i+9]*t)*Math.cos(c) + nu[i+10]*Math.sin(c);
}
dL/=10000000, dE/=10000000; //IAU2000B章动值, dL是黄经章动,dE是交角章动
var NU=new Array( //IAU1980章动序列
//l l F D Om psin t_sin, ecos t_cos
0, 0, 0, 0, 1,-171996,-1742,92025, 89,//1
0, 0, 2,-2, 2, -13187, -16, 5736,-31,//2
0, 0, 2, 0, 2, -2274, -02, 977, -5,//3
0, 0, 0, 0, 2, 2062, 02, -895, 5,
0,-1, 0, 0, 0, -1426, 34, 54, -1,
1, 0, 0, 0, 0, 712, 01, -7, 0,
0, 1, 2,-2, 2, -517, 12, 224, -6,
0, 0, 2, 0, 1, -386, -4, 200, 0,
1, 0, 2, 0, 2, -301, 0, 129, -1,
0,-1, 2,-2, 2, 217, -5, -95, 3,//10
-1, 0, 0, 2, 0, 158, 0, -1, 0,
0, 0, 2,-2, 1, 129, 1, -70, 0,
-1, 0, 2, 0, 2, 123, 0, -53, 0,//13
1, 0, 0, 0, 1, 63, 1,-33,0,//14
0, 0, 0, 2, 0, 63, 0, -2,0,
-1, 0, 2, 2, 2,-59, 0, 26,0,
-1, 0, 0, 0, 1,-58,-1, 32,0,
1, 0, 2, 0, 1,-51, 0, 27,0,
-2, 0, 0, 2, 0,-48, 0, 1,0,
-2, 0, 2, 0, 1, 46, 0,-24,0,
0, 0, 2, 2, 2,-38, 0, 16,0,
2, 0, 2, 0, 2,-31, 0, 13,0,
2, 0, 0, 0, 0, 29, 0, -1,0,
1, 0, 2,-2, 2, 29, 0,-12,0,
0, 0, 2, 0, 0, 26, 0, -1,0,
0, 0, 2,-2, 0,-22, 0, 0,0,
-1, 0, 2, 0, 1, 21, 0,-10,0,
0, 2, 0, 0, 0, 17,-1, 0,0,
0, 2, 2,-2, 2,-16, 1, 7,0,
-1, 0, 0, 2, 1, 16, 0, -8,0,
0, 1, 0, 0, 1,-15, 0, 9,0,
1, 0, 0,-2, 1,-13, 0, 7,0,
0,-1, 0, 0, 1,-12, 0, 6,0,
2, 0,-2, 0, 0, 11, 0, 0,0,
-1, 0, 2, 2, 1,-10, 0, 5,0,
1, 0, 2, 2, 2, -8, 0, 3,0,
0,-1, 2, 0, 2, -7, 0, 3,0,
0, 0, 2, 2, 1, -7, 0, 3,0,
1, 1, 0,-2, 0, -7, 0, 0,0,
0, 1, 2, 0, 2, 7, 0, -3,0,//40
-2, 0, 0, 2, 1, -6, 0, 3,0,
0, 0, 0, 2, 1, -6, 0, 3,0,
2, 0, 2,-2, 2, 6, 0, -3,0,
1, 0, 0, 2, 0, 6, 0, 0,0,
1, 0, 2,-2, 1, 6, 0, -3,0,
0, 0, 0,-2, 1, -5, 0, 3,0,
0,-1, 2,-2, 1, -5, 0, 3,0,
2, 0, 2, 0, 1, -5, 0, 3,0,
1,-1, 0, 0, 0, 5, 0, 0,0,
1, 0, 0,-1, 0, -4, 0, 0,0,
0, 0, 0, 1, 0, -4, 0, 0,0,
0, 1, 0,-2, 0, -4, 0, 0,0,
1, 0,-2, 0, 0, 4, 0, 0,0,
2, 0, 0,-2, 1, 4, 0, -2,0,
0, 1, 2,-2, 1, 4, 0, -2,0,
1, 1, 0, 0, 0, -3, 0, 0,0,
1,-1, 0,-1, 0, -3, 0, 0,0,
-1,-1, 2, 2, 2, -3, 0, 1,0,
0,-1, 2, 2, 2, -3, 0, 1,0,
1,-1, 2, 0, 2, -3, 0, 1,0,//60
3, 0, 2, 0, 2, -3, 0, 1,0,
-2, 0, 2, 0, 2, -3, 0, 1,0,
1, 0, 2, 0, 0, 3, 0, 0,0,
-1, 0, 2, 4, 2, -2, 0, 1,0,
1, 0, 0, 0, 2, -2, 0, 1,0,
-1, 0, 2,-2, 1, -2, 0, 1,0,
0,-2, 2,-2, 1, -2, 0, 1,0,
-2, 0, 0, 0, 1, -2, 0, 1,0,
2, 0, 0, 0, 1, 2, 0, -1,0,
3, 0, 0, 0, 0, 2, 0, 0,0,
1, 1, 2, 0, 2, 2, 0, -1,0,
0, 0, 2, 1, 2, 2, 0, -1,0,
1, 0, 0, 2, 1, -1, 0, 0,0,
1, 0, 2, 2, 1, -1, 0, 1,0,
1, 1, 0,-2, 1, -1, 0, 0,0,
0, 1, 0, 2, 0, -1, 0, 0,0,
0, 1, 2,-2, 0, -1, 0, 0,0,
0, 1,-2, 2, 0, -1, 0, 0,0,
1, 0,-2, 2, 0, -1, 0, 0,0,
1, 0,-2,-2, 0, -1, 0, 0,0,
1, 0, 2,-2, 0, -1, 0, 0,0,
1, 0, 0,-4, 0, -1, 0, 0,0,
2, 0, 0,-4, 0, -1, 0, 0,0,
0, 0, 2, 4, 2, -1, 0, 0,0,
0, 0, 2,-1, 2, -1, 0, 0,0,
-2, 0, 2, 4, 2, -1, 0, 1,0,
2, 0, 2, 2, 2, -1, 0, 0,0,
0,-1, 2, 0, 1, -1, 0, 0,0,
0, 0,-2, 0, 1, -1, 0, 0,0,
0, 0, 4,-2, 2, 1, 0, 0,0,
0, 1, 0, 0, 2, 1, 0, 0,0,
1, 1, 2,-2, 2, 1, 0, -1,0,
3, 0, 2,-2, 2, 1, 0, 0,0,
-2, 0, 2, 2, 2, 1, 0, -1,0,
-1, 0, 0, 0, 2, 1, 0, -1,0,
0, 0,-2, 2, 1, 1, 0, 0,0,
0, 1, 2, 0, 1, 1, 0, 0,0,
-1, 0, 4, 0, 2, 1, 0, 0,0,
2, 1, 0,-2, 0, 1, 0, 0,0,
2, 0, 0, 2, 0, 1, 0, 0,0,
2, 0, 2,-2, 1, 1, 0, -1,0,
2, 0,-2, 0, 1, 1, 0, 0,0,
1,-1, 0,-2, 0, 1, 0, 0,0,
-1, 0, 0, 1, 1, 1, 0, 0,0,
-1,-1, 0, 2, 1, 1, 0, 0,0, //105
0, 1, 0, 1, 0, 1, 0, 0,0 //106
);
var dL2=0,dE2=0;
for(i=0;i<106*9;i+=9){
c = ( NU[i]*l + NU[i+1]*l1 + NU[i+2]*F +NU[i+3]*D +NU[i+4]*Om )/rad;
dL2 += ( NU[i+5] + NU[i+6]*t/10 )*Math.sin(c);
dE2 += ( NU[i+7] + NU[i+8]*t/10 )*Math.cos(c);
}
dL2/=10000, dE2/=10000; //IAU1980的章动值
//输出,显示
document.write("IAU2000B与IAU1980章动序列的比较,单位是角秒:<br><br>");
document.write("J2000起算的儒略世纪数:t=" + t + "<br><br>");
ddL=dL-dL2; ddE=dE-dE2; //二者之差
document.write("IAU2000B黄经章动:" + dL + " 交角章动:" + dE + "<br>");
document.write("IAU1980 黄经章动:" + dL2 + " 交角章动:" + dE2 + "<br><br>");
document.write("二者之差:" + ddL + " 二者之差:" + ddE + "<br>");
</script>