画太阳算法

㈠ 太阳高度角的赤纬算法

上述式子中都涉及太阳赤纬，太阳赤纬的算法如下：
由于太阳赤纬角在周年运动中任何时刻的具体值都是严格已知的，所以它（ED）也可
以用与式（1）相类似的表达式表述，即：
ED=0.3723+23.2567sinθ+0.1149sin2θ－0.1712sin3θ－0.758cosθ+0.3656cos
2θ+0.0201cos3θ（5）
式中θ称日角，即 θ=2πt/365.2422（2）
这里t又由两部分组成，即 t=N－N0 （3）
式中N为积日，所谓积日，就是日期在年内的顺序号，例如，1月1日其积日为1，平年12月31日的积日为365，闰年则为366，等等。
N0=79.6764+0.2422×（年份－1985）－INT〔（年份－1985）/4〕
（式中INT表示取整数部分，例如INT(3.25)=3）