这个公式里所有的都是弧度值,h,φ,δ,都是弧度值.
t表示小时,也就是几点,24小时制.
φ直接可知,
δ要用公式δ=0.3723+23.2567sinθ+0.1149sin2θ-0.1712sin3θ-0.758cosθ+0.3656cos2θ+0.0201cos3θ计算,
其中θ=2πt/365.2422,
这里t又由两部分组成,即 t=N-N0,N0=79.6764+0.2422×(年份-1985)-INT〔(年份-1985)/4〕(式中INT表示取整数部分,例如INT(3.25)=3)