日出日落時(shí)間計(jì)算程序(c語(yǔ)言)-文庫(kù)吧
2025-06-22 14:18 本頁(yè)面
【正文】 else { for(i=0。imonth1。i++){ a=a+days_of_month_1[i]。 } } a=a+date。 return a。 } //求從格林威治時(shí)間公元2000年1月1日到計(jì)算日天數(shù)days long double t_century(int days, long double UTo){ return ((long double)days+UTo/360)/36525。 } //求格林威治時(shí)間公元2000年1月1日到計(jì)算日的世紀(jì)數(shù)t long double L_sun(long double t_century){ return (+*t_century)。 } //求太陽(yáng)的平黃徑 long double G_sun(long double t_century){ return (+*t_century)。 } //求太陽(yáng)的平近點(diǎn)角 long double ecliptic_longitude(long double L_sun,long double G_sun){ return (L_sun+*sin(G_sun*PI/180)+*sin(2*G_sun*PI/180))。 } //求黃道經(jīng)度 long double earth_tilt(long double t_century){ return (*t_century)。 } //求地球傾角 long double sun_deviation(long double earth_tilt, long double ecliptic_longitude){ return (180/PI*asin(sin(PI/180*earth_tilt)*sin(PI/180*ecliptic_longitude)))。 } //求太陽(yáng)偏差 long double GHA(long double UTo, long double G_sun, long double ecliptic_longitude){ return (*sin(G_sun*PI/180)*sin(2*G_sun*PI/180)+*sin(2*ecliptic_longitude*PI/180)*sin(4*ecliptic_longitude*PI/180))。 } //求格林威治時(shí)間的太陽(yáng)時(shí)間角GHA long double e(long double h, long double glat, long double sun_deviation){ return 180/PI*acos((sin(h*PI/180)sin(gl
數(shù)學(xué)相關(guān)推薦
鄂ICP備17016276號(hào)-1