Continued Fraction
连分数¶
连分数是一种记号。例如,长为
只是为了形式上简洁,才记成等号左边的样子。这里的四个变元可以任意取值。
连分数中,从前面一直到第
渐进分数有递推关系:
这个式子和 Farey 数列的递推关系很像。
上面的例子是有限连分数。如果分式无限地写下去,有无限个变元,就得到无限连分数。无限连分数收敛等价于渐进分数收敛。
有定理:
无限连分数,如果各变元均大于等于
因为只要各变元为正,无限连分数的偶渐进分数单调递增(都比它小),奇渐进分数单调递减(都比它大)。而在均大于等于
显然可以看到,连分数关于下标为偶数的变元单调递增,关于下标为奇数的变元单调递减。这无论它有限或无限都成立。
还有一个事实,对于有限连分数,全体结尾为
简单连分数:连分数从第
简单连分数的值,一定大于偶数的渐进分数,一定小于奇数的渐进分数。无限简单连分数一定收敛。
仿照一般分数的概念,第
如果“规定第
如果两个无限简单连分数的值相等,必然逐项相等。
如果两个有限简单连分数的值相等,不仅要逐项相等,而且必然项数也相同。
无限简单连分数不能与有限简单连分数值相等。有理数与有限简单连分数具有一一对应关系,因此无限简单连分数全都是无理数。
如果要求
- 取倒数,得到的数大于
1 - 取整,得到的数在
0 1 - 对余数重复上述操作。
这样就得到了相应的表示。
循环连分数¶
仿照循环小数的概念,如果在连分数后面形成了循环,则形成 循环连分数。如果循环节从第 0 项开始,称为 纯循环连分数,否则称为 混循环连分数。例如纯循环连分数:
混循环连分数:
混循环连分数后面循环的部分,一般要找最早循环的部分,称为它的“纯循环部分”。循环节一般取最小正整数。
二次有理数¶
形如
对于二次有理数
以下的二次有理数,特指实二次有理数,即
对二次有理数执行“无限简单连分数”计算,即取倒数、取整交替,得到的余数还是二次有理数。
这里考虑一类特殊形式的二次有理数。设
称
只要让分母整除分子的范数,取倒数、取整的交替过程会变得简单很多。事实上如果这条件不成立,只要分子分母同时乘以分母的绝对值,并强行压入根号,条件就成立了。因此,任何二次有理数都能写成这形式。
如果这条件成立,那么每个余数都满足这条件。
根据上文的简单连分数算法:对余数取整可以得到
取整后得到的新的
由于范数为负,取倒数之后
余数循环显然是循环连分数的循环条件。由于计算方法固定,只要余数重复必循环。
余数
余数有限必重复。因此 所有二次有理数都可以写成循环连分数。反过来容易证明 循环连分数一定是二次有理数。
由于无限简单连分数写法唯一,因此就找到了二次有理数的一种唯一表示。
最佳有理逼近¶
讨论如何用有理数“最佳地”逼近无理数,不妨假设无理数落入
“最佳”一词的概念:选定的有理数必须保证,比它与无理数的距离更近的有理数,分母都比它大。不存在分母比它小的有理数,到给定无理数的距离更近。
最佳有理数:在 Farey 数列的某一行中,与给定无理数距离最近的那个有理数。
这个有理数可能在下面几行依旧与无理数“距离最近”,但一定有某一行,会找到一个新的有理数,与无理数距离更近。因此去重复后可以得到 最佳逼近有理数列,分母严格递增,距离递减。
比无理数大的称为上逼近,否则为下逼近。由于无理数和有理数之间一定有有理数,最佳逼近有理数列必然为若干个上逼近,之后若干个下逼近,交替进行的形式。
有结论:
渐进分数必然为最佳逼近。偶项渐进分数全都是下逼近,奇项渐进分数全都是上逼近。渐进分数列是下上交错的逼近。
在最佳逼近有理数列中,渐进分数是下上关系改变之前的倒数第一个数。如果将最佳逼近有理数列都写成有限简单连分数,那么在渐进分数之后(下上关系改变之后),连分数长度加一。
例如,
那么它的渐进分数列是:
它的最佳逼近有理数列是:
从每个渐进分数(不包含)开始,到下一个渐进分数(包含)为止,同为上逼近,或同为下逼近。
在最佳逼近列中,每一个最佳分数是上一个最佳分数与再往前一个渐进分数的分子分母对应求和。
Dirichlet 逼近定理¶
Dirichlet 逼近定理是说,对于任意的一个无理数
更强的结论是,在右边的分母上能放一个
更进一步的定理是 Kronecker 的逼近定理。
如果
这两个定理都可以用抽屉原理来解决。事实上,历史上第一次正式提出抽屉原理,就是 Dirichlet 为了解决 Pell 方程 而研究这个有理数逼近条件才正式提出来抽屉原理的。
当然,采用抽屉原理的证明可以发现,上文中提到的最佳逼近有理数列,每项满足定理中右边改成分母为一次式的不等式。
进一步有结论,渐进有理数列中,每一项均满足 Dirichlet 逼近定理的不等式。
buildLast update and/or translate time of this article,Check the history
editFound smelly bugs? Translation outdated? Wanna contribute with us? Edit this Page on Github
peopleContributor of this article OI-wiki
translateTranslator of this article Visit the original article!
copyrightThe article is available under CC BY-SA 4.0 & SATA ; additional terms may apply.