久久久久久久999_99精品久久精品一区二区爱城_成人欧美一区二区三区在线播放_国产精品日本一区二区不卡视频_国产午夜视频_欧美精品在线观看免费

 找回密碼
 立即注冊

QQ登錄

只需一步,快速開始

搜索
查看: 7430|回復(fù): 9
打印 上一主題 下一主題
收起左側(cè)

51單片機可否實現(xiàn)上下五千年農(nóng)歷計算(非查表法)

[復(fù)制鏈接]
跳轉(zhuǎn)到指定樓層
樓主
ID:67796 發(fā)表于 2014-12-6 11:39 | 只看該作者 回帖獎勵 |倒序瀏覽 |閱讀模式
本帖最后由 exv 于 2014-12-6 12:58 編輯

看了論壇中幾個有農(nóng)歷的時鐘程序,無一例外使用查表法實現(xiàn)的,且只能查上下100年(因為DS1302只有上下一百年)。

看過莆田第十中學(xué)許劍偉老師的農(nóng)歷算法,可實現(xiàn)上下五千年農(nóng)歷的計算(含節(jié)氣),
與主函數(shù)直接相關(guān)的C代碼就有6000余行;
關(guān)于日月星辰軌道計算的代碼有數(shù)十組,
每一組換算往往涉及幾十個64位double。

不知道51單片機在 合理 配置的情況下,能否實現(xiàn)這樣的計算?

如果不行,實現(xiàn)這一功能的最低配置是?

許劍偉老師的代碼這兒就不貼了(估計也貼不下),百度下 莆田第十中學(xué)許劍偉老師 即可。

51hei人才濟濟,下面這段計算圓周率的代碼大家應(yīng)該非常熟悉,古人幾十本書還沒說清的一兩行就搞定了:
  1. long a=10000,b,c=2800,d,e,f[2801],g;
  2. void main(){for(;b-c;)f[b++]=a/5;for(;d=0,g=c*2;c-=14,cout<<e+d/a,e=d%a)for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);}
復(fù)制代碼


古人制農(nóng)歷最多也就是算盤,難道 51單片機 真敵不過古人算盤?

另外,許老師的日歷程序除了可定農(nóng)歷日期/節(jié)氣,還可依此結(jié)合設(shè)定的經(jīng)緯度換算出每天的日落日出,非常實用。

分享到:  QQ好友和群QQ好友和群 QQ空間QQ空間 騰訊微博騰訊微博 騰訊朋友騰訊朋友
收藏收藏2 分享淘帖 頂 踩

相關(guān)帖子

回復(fù)

使用道具 舉報

沙發(fā)
ID:67992 發(fā)表于 2014-12-6 12:25 來自觸屏版 | 只看該作者
不會,太高深了。算盤結(jié)果可暫記紙上,記多少都可以,容量無限大,今天記不完明天接著記,速度可以無限慢。如果比這些,51單片機還真拍馬莫及。

評分

參與人數(shù) 1黑幣 +50 收起 理由
admin + 50 回帖助人的獎勵!

查看全部評分

回復(fù)

使用道具 舉報

板凳
ID:67796 發(fā)表于 2014-12-6 12:45 | 只看該作者
dgahz 發(fā)表于 2014-12-6 12:25
不會,太高深了。算盤結(jié)果可暫記紙上,記多少都可以,容量無限大,今天記不完明天接著記,速度可以無限慢。 ...

程序本身并不高深,有高中水平即能看懂(不然 莆田第十中學(xué)的學(xué)生們可就悲催了),也就是些牛頓方程什么的。

只是這么大的計算量,放在RAM有限的51單片機確實是一個麻煩,也難怪目前能見到的都是查表法實現(xiàn)的。

不知如果借助一個2M的EEPROM能否解決這個問題?

評分

參與人數(shù) 1威望 +50 黑幣 +50 收起 理由
admin + 50 + 50 回帖助人的獎勵!

查看全部評分

回復(fù)

使用道具 舉報

地板
ID:67992 發(fā)表于 2014-12-6 13:40 來自觸屏版 | 只看該作者
農(nóng)歷是怎么算出來的我真搞不懂,和牛頓有關(guān)?印象中牛頓是力學(xué)噸
回復(fù)

使用道具 舉報

5#
ID:67992 發(fā)表于 2014-12-6 13:41 來自觸屏版 | 只看該作者
出社會太久,知識都丟了。
回復(fù)

使用道具 舉報

6#
ID:69706 發(fā)表于 2014-12-6 16:31 | 只看該作者
應(yīng)該是不管用吧
回復(fù)

使用道具 舉報

7#
ID:67796 發(fā)表于 2014-12-6 16:40 | 只看該作者
本帖最后由 exv 于 2014-12-6 16:46 編輯

還是貼一些莆田十中許老師的代碼給大家溫習(xí)下相關(guān)知識:

1.算出太陽/月亮/地球的位置(牛頓力學(xué))
2.反算太陽月亮視黃經(jīng),相等時刻定月初一
3.根據(jù)各個月初一導(dǎo)出整個農(nóng)歷,節(jié)氣可從太陽黃經(jīng)直接導(dǎo)出。

許老師的代碼雖然多,有相當部分是固定不變的天文/地理參數(shù),包括一大堆校正因子集,故而前面的帖子提出使用EEPROM作存儲輔助。

既然圓周率求算都能從幾十本書變?yōu)閮尚写a,農(nóng)歷又是極為常用的東西,必定有非常多人去研究,一個簡潔高效能在單片機上運行的代碼應(yīng)該是已經(jīng)有人寫出過的。

下面代碼是通過地球軌道方程反算太陽黃經(jīng),從而計算出各個節(jié)氣。
(51hei只見到一個歷法程序有節(jié)氣,只是人家沒給代碼,不知道是不是用下面方法算出的)


  1. double qi_low(double W)
  2. {
  3.         //最大誤差小于30分鐘,平均5分
  4.         double t,L,v= 628.3319653318;
  5.         t =  ( W - 4.895062166 )/v; //第一次估算,誤差2天以內(nèi)
  6.         t -= ( 53*t*t + 334116*cos( 4.67+628.307585*t) + 2061*cos( 2.678+628.3076*t)*t )/v/10000000; //第二次估算,誤差2小時以內(nèi)
  7.         
  8.         L = 48950621.66 + 6283319653.318*t + 53*t*t //平黃經(jīng)
  9.                 +334166 * cos( 4.669257+  628.307585*t) //地球橢圓軌道級數(shù)展開
  10.                 +3489 * cos( 4.6261  + 1256.61517*t ) //地球橢圓軌道級數(shù)展開
  11.                 +2060.6 * cos( 2.67823 +  628.307585*t ) * t  //一次泊松項
  12.                 - 994 - 834*sin(2.1824-33.75705*t); //光行差與章動修正
  13.         
  14.         t -= (L/10000000 -W )/628.332 + (32*(t+1.8)*(t+1.8)-20)/86400/36525;
  15.         return t*36525 + (double)8/24;
  16. }
復(fù)制代碼


評分

參與人數(shù) 1威望 +50 黑幣 +50 收起 理由
admin + 50 + 50 回帖助人的獎勵!

查看全部評分

回復(fù)

使用道具 舉報

8#
ID:74921 發(fā)表于 2015-3-20 11:26 | 只看該作者
好高深的有沒有?反正我是看迷糊了
回復(fù)

使用道具 舉報

9#
ID:67796 發(fā)表于 2015-3-20 17:48 | 只看該作者
a719563181 發(fā)表于 2015-3-20 11:26
好高深的有沒有?反正我是看迷糊了

其實,是有的問題被復(fù)雜化了。

農(nóng)歷本來是用來指導(dǎo)一年的農(nóng)時,大多數(shù)時候也不需要精確到分秒。
如同買鞋子的時候最多用皮尺之類的工具量量尺寸,沒人會去用游標卡尺。

這一步
1.算出太陽/月亮/地球的位置(牛頓力學(xué))

為了得到精確的值,使用近代天文學(xué)算法(多體),耗費了大量計算資源。
現(xiàn)代算法是從《史記歷書》演化而來,只是司馬遷直接把太陽月亮假設(shè)為勻速圓周運動,通過這樣一個近似,計算量極大減小,所得的歷法也滿足指導(dǎo)農(nóng)時之用。

當然,現(xiàn)代算法也是一種學(xué)問上的嚴謹。日月食的時間能精確到秒(參看許老師代碼)。

評分

參與人數(shù) 1黑幣 +50 收起 理由
admin + 50 回帖助人的獎勵!

查看全部評分

回復(fù)

使用道具 舉報

10#
ID:199295 發(fā)表于 2017-5-11 11:01 | 只看該作者
很好的資料
回復(fù)

使用道具 舉報

您需要登錄后才可以回帖 登錄 | 立即注冊

本版積分規(guī)則

小黑屋|51黑電子論壇 |51黑電子論壇6群 QQ 管理員QQ:125739409;技術(shù)交流QQ群281945664

Powered by 單片機教程網(wǎng)

快速回復(fù) 返回頂部 返回列表
主站蜘蛛池模板: 久久国产免费 | 成年人视频免费看 | 国产福利在线视频 | 艳妇臀荡乳欲伦交换h漫 | 天天色播 | 国产日韩免费 | 日韩欧美一区二区三区久久婷婷 | 在线观看国产小视频 | av在线播放不卡 | 精品黄色片 | 97精品在线 | 黄色网在线 | 国产精品伦一区二区三级视频 | 日日干狠狠干 | 美女视频一区二区 | 午夜专区| 男女啪啪网站 | 日韩在线成人 | av一二三 | 亚洲精品久久久久avwww潮水 | 精品日韩一区二区 | 精品国产一区二 | 欧美成人一区二区 | 精品视频免费在线观看 | 日韩欧美一区二区在线 | 80日本xxxxxxxxx96| 免费一级片 | 亚洲欧美国产毛片在线 | 午夜成人免费视频 | 少妇高潮久久久久久潘金莲 | 蜜臀av性久久久久av蜜臀妖精 | 干少妇视频 | www国产| 国产精品羞羞答答 | 亚洲免费久久 | 成人午夜激情视频 | 四虎永久在线 | 欧美日韩在线看 | www免费视频 | 亚洲精品不卡 | 在线免费黄色 |