畫太陽演算法

㈠ 太陽高度角的赤緯演算法

上述式子中都涉及太陽赤緯，太陽赤緯的演算法如下：
由於太陽赤緯角在周年運動中任何時刻的具體值都是嚴格已知的，所以它（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）