查看“︁调日法”︁的源代码

'''调日法'''<ref>中國古时将天文数据的小数部分的[[分母]]称为「日」，「调日术」即是调节分母的意思。</ref>是[[南北朝]][[数学家]][[何承天 (南朝)|何承天]]发明的一种系统地寻找最佳逼近（[[帕德近似]]）以表示天文数据或[[数学常数]]的[[内插]]法。据宋史卷七十四：“宋世何承天，更以四十九分之二十六为强率，十七分之九为弱率；于强弱之际，以求日法……自后治历者，莫不因承天法，累强弱之数”调日法后來传入[[日本]]。
[[File:Diaorifa.GIF|thumb|right|300px|]]
中国有学者认为[[祖冲之]]可能利用何承天的调日法求得[[圆周率]]的约率和密率：
:圆周率的约率为 <math> \frac{22}{7}\ </math>
:圆周率的密率为 <math> \frac{355}{113}\ </math>

何承天的调日法是他对数学的一项重要贡献。一千年以后，15世纪法国数学家[[尼古拉·许凯]]（1455年 ━ 1488年），才使用相似的插入法。

==何承天调日法原理==
已知<math> \frac{a}{b}< \frac{c}{d} </math>

则<math> \frac{a}{b}< \frac{a+c}{b+d}< \frac{c}{d}</math> 

推而广之：
: <math> \frac{a}{b}< \frac{ma+kc}{mb+kd}< \frac{c}{d}</math>，其中 m,k 为正整数。 

欲求精确分数<math>f_n</math>使<math>|f-f_n|<\delta</math>，其中<math>\delta</math>为误差界限。

令<math>f_0= \frac{a}{b} </math>为弱率，<math>f_1=\frac{c}{d} </math>为强率。

第一步，根据下列方法求得一个近似分数

: <math>f_2= \frac{a+c}{b+d} </math> 

如果<math>f_2>f</math>，则将<math>f_2= \frac{a+c}{b+d}\ </math> 作为新的强分数，和旧弱分数 <math>\frac{a}{b}\ </math>'''调日'''得到近似分数：

: <math>f_3 = \frac{a+c+a}{b+d+b}\ </math>

如果<math>f_2<f</math>， 则将<math> f_2 = \frac{a+c}{b+d}\ </math> 作为新的弱分数，和旧强分数 <math> \frac{c}{d}\ </math> '''调日'''得到近似分数：

: <math>f_3 = \frac{a+c+c}{b+d+d}\ </math> 

反复操作，到<math>|f-f_n|<\delta</math> 为止。

另外，还可以直接求m,k的数值，加快逼近速度：
若<math>\frac{a}{b}<x<\frac{c}{d}</math>，且 <math>x-\frac{a}{b}=d_{1}</math>以及<math>\frac{c}{d}-x=d_{2}</math>

如果有正整数m,k满足：<math>\frac{kd}{mb}=\frac{d_{1}}{d_{2}}</math>

那么就有：<math>x=\frac{ma+kc}{mb+kd} </math>

证明如下：由条件可得

: <math>\begin {align}
bd_{1}&=bx-a\\
dd_{2}&=c-dx
\end{align}</math>

而根据<math>\frac{kd}{mb}=\frac{d_{1}}{d_{2}}</math>又有

:<math>mbd_{1}=kdd_{2}</math>

代入上面的两个关系式可得：

:<math>m(bx-a)=k(c-dx)</math>

解关于x的一元一次方程就有结果：
:<math>x=\frac{ma+kc}{mb+kd}</math>

==应用==
何承天调日法被同时代和后代数学家如[[赵爽]]，[[祖冲之]]，[[一行]]等运用。
===朔望月===
何承天将 <math> \frac{9}{17}=0.529412...\ </math>作为朔望月零数部分的弱率，以<math> 
\frac{26}{49}=0.530612...\ </math>作为朔望月零数部分的强率。运用调日法，最后得到 <math> \frac{399}{752}\ </math>，根据他的观测数值0.530585，首先计算d1,d2

:<math>\begin{align}
d_{1}&=0.530585-0.529412&=0.001173\\
d_{2}&=0.530612-0.530585&=0.000027
\end{align}
</math>

寻找满足以下关系的m,k值：
:<math>
\begin{align}
\frac{49k}{17m}&=\frac{1173}{27}\\
\frac{k}{m}&=\frac{1173\times17}{49\times27}&\approx{}15.07\ldots
\end{align}
</math>

可以令(m,k)=(1,15)

从而得到：

:<math>\frac{1\times9+15\times26}{1\times17+15\times49}=\frac{399}{752}</math>

727年[[唐朝]][[天文学家]][[一行]]在《[[大衍历]]》中用同样的弱率和强率求得  <math>\frac{1613}{3040}</math>

===闰周问题===
南北朝数学家祖冲之熟悉调日术，他以<math> \frac{4}{11}</math>为弱率， 
以<math> \frac{7}{19}</math>为强率，通过调日法得到<math> \frac{144}{391} </math>

===近点月===
何承天以<math>\frac{56}{101} </math>为弱率，以<math>\frac{5}{9}</math>为强率，用调日法求得近点月为<math>\frac{417}{752} </math>。祖冲之也得到高精度的数值<math>\frac{14631}{26377} </math>

===圆周率约率和密率===
祖冲之求圆周率约率和密率的方法已失传。有学者认为他用[[刘徽割圆术]]求得圆周率的约率和密率 ；也有学者认为祖冲之有可能用何承天的调日法求得圆周率的约率和密率的分数表示式<ref>[[吴文俊]] 主编 《[[中国数学史大系]]》第四卷 123页，ISBN7-300-0425-8/O</ref>。 祖冲之对调日法是熟悉的，他自己就用过调日法改进何承天近点月<math> \frac{417}{752} </math>为更加精确的<math>\frac{14631}{26377} </math>

取<math>\pi\approx3.1416</math>，先只考虑小数部分，根据<math>\frac{1}{8}<0.1416<\frac{1}{7}</math>，用调日法进行计算： 
:<math>
\begin{align}
d_{1}&=0.1416-0.125&=0.0166\\
d_{2}&=0.142857-0.1416&=0.001257
\end{align}
</math>

寻找满足以下关系的m,k值：

:<math>
\begin{align}
\frac{7k}{8m}=\frac{0.0166}{0.001257}\\
\frac{k}{m}=\frac{8\times0.0166}{7\times0.001257}&\approx{}15.09\ldots
\end{align}
</math>

所以可以令(m,k)=(1,15)，从而可以得到结果：

<math>3+\frac{1\times1+1\times15}{1\times8+15\times7}=3+\frac{16}{113}=\frac{355}{113}</math>

祖冲之密率 <math>\frac{355}{113}</math> 和π之误差为0.0000002668。下一个{{citation needed}}比之更為精确的分数为 <math>\frac{52163}{16604}= 3.1415923874</math> 误差为 -0.0000002662，分子、分母都比祖冲之密率的分子、分母複雜得多。

祖冲之很可能先用刘徽割圆术求出圆周率。刘徽割圆术计算需要多次开平方运算，例如用八次割圆术得到 <math>\pi \approx \frac{3927}{1250}=3.1416</math><ref>傅海伦编著 《中外数学史概论》 第四章 刘徽的割圆术  51页 科学出版社，ISBN978-7-03-018477-1</ref>， 无论分子分母都比祖冲之密率的分子分母复杂，但还不如密率的分数表示准确。用十一次割圆术可得到和密率相当精确但比较复杂的分数，再通过调日法求得准确而又简单的分数式。

调日法后传入[[日本]]。日本数学家[[关孝和]]（Seki, Takakazu, 1642-1708）在《括要算法》一书中称之为'''零约术'''，并用之得出圆周率的近似分数为 <math>  \frac{355}{113}</math><ref>
[[吴文俊]] 主编 《中国数学史大系》第四卷 125页，ISBN7-300-0425-8/O</ref>，正是祖冲之的密率。

==黄金分割与[[斐波那契数列]]==
[[黄金分割]]：
:<math>\varphi = \frac{\sqrt{5} + 1}{2} \approx 1.6180339887...</math>
用调日法求分数表示：
:<math>\frac{1}{1},\frac{2}{1},\frac{3}{2},\frac{5}{3},\frac{8}{5},\frac{13}{8},\frac{21}{13},\frac{34}{21},\frac{55}{34},\frac{89}{55},\frac{144}{89},\frac{233}{144},\frac{377}{233},\frac{610}{377},\frac{987}{610},\frac{1597}{987},\frac{2584}{1597},\frac{4181}{2584}</math>

分母1,2,3,5,8,13,21,....正是[[斐波那契数列]]。

==其他==
*√2=1.4142135623 ~=<math>{99 \over 70} </math>
*√3=1.7320508075 ~=<math>{71 \over 41} </math>
*√5=2.2360679775 ~=<math>{199 \over 89} </math>
*√10=3.162277660 ~=<math>{117 \over 37} </math>
*<math>\sqrt [12] {2}</math>=1.059463094~=<math>{107 \over 101}</math>
*e=2.718281828 ~=<math>{2721 \over 1001}</math>
*[[普朗克常数]] ~=<math>{53 \over 8}</math>x10<sup>-34</sup>
*[[万有引力常数]] G~=<math>{227 \over 34}</math>x10<sup>-11</sup>
*[[阿伏伽德罗常量]]~=<math>{241 \over 40}</math>x10<sup>23</sup>
*[[玻尔兹曼常数]]~=<math>{29 \over 21}</math>x10<sup>-23</sup>

==参考文献==
{{Reflist}}

{{中国数学史}}

[[Category:中國古代數學]]
[[Category:数学近似]]