疯了，竟然算不出昼长。。

qinkong

在气象学书上找到的昼长计算公式：
cost=-tgδ·tgφ(其中δ为太阳赤纬，φ为纬度)

昼长T=2t=2arccos(-tgδ*tgφ)

太阳赤纬δ=[0.006918-0.399912cos(b)+0.070257sin(b)-0.006758cos(2b)+0.000907sin(2b)-0.002697cos(3b)+0.00148sin(3b)](180/pi)；其中b=2pi(N-1)/365；N为天数，1月1日为0.

相应的我写了个matlab函数：

function T=daylength(lat,day)   %其中lat为站点纬度(单位为弧度)，day为天数（1月1日为0）

x=tan(lat);b=2*pi*(day-1)/365; %b单位为弧度

delta=(0.006918-0.399912*cos(b)+0.070257*sin(b)-0.006758*cos(2*b)+0.000907*sin(2*b)-0.002697*cos(3*b)+0.00148*sin(3*b))*(180/pi);   %单位为度

y=tan(delta/180*pi);

T=2*acos(-x*y);
但是结果一直不对啊，比如赤道处昼长为12小时。可计算出来的不对啊，如下。

>> daylength(0,100)
ans =
    3.1416
其他的结果也不对啊，acos函数的值域为0-pi，那算出来的昼长最大也才2pi啊，这个不可能啊。迷惑了。。。求解！

qinkong

自顶！！！！！

topmad

我在计算的时候  使用的日序是1月1日为1
虽然没有仔细看你的程序 但是里面有个 day-1  如果是1月1日为0  不知道会不会有大问题

我用fortran算过 没有问题

qinkong

topmad 发表于 2014-7-19 09:51
我在计算的时候  使用的日序是1月1日为1
虽然没有仔细看你的程序 但是里面有个 day-1  如果是1月1日为0   ...

。。这个1月1日为0没啥大影响，比如我算赤道第100天昼长的话，无非就是前后差一天的关系。但是我直观的从这个教科书中给出的公式就觉得不对啊，昼长=2arccos(**)。arcos函数的值域是0-pi，那么昼长最多也就2pi，6个多小时。。。这根本不对啊。。。。你怎么算对了。。。

topmad

qinkong 发表于 2014-7-19 14:52
。。这个1月1日为0没啥大影响，比如我算赤道第100天昼长的话，无非就是前后差一天的关系。但是我直观的从 ...

等会我把我的方法贴上来

topmad

Ws=arcos（-tanφtanδ）

.δ=0.409sin（0.0172J-1.39）

N=(24/π)Ws.

φ是纬度   
N是日照时数  与昼长应该是一样的   

看公式的时候 突然想起来一个东西  计算的时候 要注意  有一个地方是要取整的  具体忘了是哪里  

这个公式貌似比较简单..。比你那个简化了一些东西

qinkong

topmad 发表于 2014-7-20 09:59
Ws=arcos（-tanφtanδ）.δ=0.409sin（0.0172J-1.39）N=(24/π)Ws.

φ是纬度   

咦，你的这一步N=(24/π)Ws好像和我在书上找到的公式不一样，能解释下这个公式么，括号里面的小n是什么含义。为什么我查到的就是N=2Ws呢？。。。你查到的公式是在网上的么？方便的话能贴个链接么？

谢谢！

topmad

qinkong 发表于 2014-7-20 19:02
咦，你的这一步N=(24/π)Ws好像和我在书上找到的公式不一样，能解释下这个公式么，括号里面的小n是什么含 ...

呀 貌似我给你弄错了  

一个是可照时数 一个是日照时数.....

qinkong

topmad 发表于 2014-7-21 08:11
呀 貌似我给你弄错了  

一个是可照时数 一个是日照时数.....

啥意思，可照时数应该就是昼长吧。您能不能帮我看看我的问题到底出在哪里。。我觉得我现在是当局者迷，云山雾罩啊。

qinkong

顶一下，这个问题还没解决。。。