急求日落时角的计算方法！

对于与地球纬度相关的运算，凡是想用平均值来计算的都是错的，无一例外，昼长增加量不是每天等量，而是服从三角函数规律。你应该计算那个纬度的“白昼经角”。

这一步我7年前就已经推导出来了：白昼经角 = 360 - 2arccos[tg太阳赤纬/tg(90-当地纬度)]，记住这个“白昼经角”。公式中有arccos ，当arccos中括号中的数值大于1或小于-1都是无意义，大于1的是极昼，小于-1的是极夜。

计算白昼经角就是计算某个纬度看见太阳所在的那个赤纬圈有多少度在地平线上，度数越多日照越长。

就拿你题目的20°N和40°N作例子，在冬至日。代公式360 - 2arccos[tg太阳赤纬/tg(90-当地纬度)]，太阳位于-23.5度(南纬)，即360 - 2arccos[tg(-23.5)/tg(90-20)]，得出20°N的白昼经度是161.788度；而40°N同样代公式，得出白昼经度是137.2度。

地球每小时自转15度，所以两地的白昼经度分别除以15，得出20°N的白昼时间是161.788/15 = 10.78小时，即10小时47分；40°N的白昼时间是137.2/15 = 9.15小时，即9小时09分，完全与题目吻合。

回到题目，问题就是1月15日太阳所在的赤纬坐标是多少。这一天距离春分有16+28+21=65天，所以太阳距离春分点的角距是65天/365天*360度 = 64度。

接着是另一组公式，知道太阳与春分点的角距求太阳所在的赤纬：sin赤纬 = sin黄赤交角·sin太阳黄经。这里的“太阳黄经”就是上面的太阳与春分点的距离，不过黄经的范围是0-360度，所以上面得出的64度应该是取补角，即296度。

然后代公式，sin赤纬 = sin23.5·sin296 = -0.3584 ，从而得出太阳的赤纬坐标是-21度，即直射21°S 。

之后就将太阳赤纬的-21度和当地纬度35度代入之前那个公式：360 - 2arccos[tg太阳赤纬/tg(90-当地纬度)]，即360 - 2arccos[tg(-21°)/tg(90-35)]，得出白昼经角是148.816度，再除以每小时15度