Текущая версия |
Ваш текст |
Строка 4: |
Строка 4: |
| {{Определение | | {{Определение |
| |definition= | | |definition= |
− | Действительное число <tex> \alpha </tex> называется '''вычислимым''' (англ. ''computable number''), если существует [[Вычислимые функции|вычислимая функция]] <tex> a(\varepsilon): \mathbb Q \rightarrow \mathbb Q </tex>, такая, что <tex>|\alpha - a(\varepsilon)| < \varepsilon </tex> для любого рационального <tex> \varepsilon > 0 </tex>. | + | Действительное число <tex> \alpha </tex> называется '''вычислимым''' (computable number), если существует вычислимая функция <tex> a: \mathbb Q \rightarrow \mathbb Q </tex>, такая, что <tex>|\alpha - a(\varepsilon)| < \varepsilon </tex> для любого рационального <tex> \varepsilon > 0 </tex>. |
| }} | | }} |
| | | |
Строка 11: |
Строка 11: |
| {{Теорема | | {{Теорема |
| |statement= | | |statement= |
− | Число <tex> \alpha </tex> вычислимо <tex>\iff</tex> множество <tex>A = \{x \in \mathbb Q \mid x < \alpha \} </tex> [[Разрешимые (рекурсивные) языки|разрешимо]]. | + | Число <tex> \alpha </tex> вычислимо тогда и только тогда, когда множество <tex> \{a \in Q \mid a < \alpha \} </tex>, разрешимо. |
| |proof= | | |proof= |
− | <tex> \Longrightarrow </tex>
| |
− |
| |
− | : Если <tex> \alpha </tex> {{---}} рациональное, то существует тривиальный разрешитель для <tex> A </tex>, который просто сравнивает полученный элемент с <tex> \alpha </tex>.
| |
− |
| |
− | : В противном случае, построим разрешитель для <tex> A </tex>:
| |
− |
| |
− | '''function''' <tex> p(x)</tex>:
| |
− | '''for''' <tex> n = 1</tex> '''to''' <tex>\infty </tex>
| |
− | '''if''' <tex> x < a \left(\dfrac{1}{n} \right) - \dfrac{1}{n} </tex>
| |
− | '''return''' <tex>1</tex>
| |
− | '''if''' <tex> x > a \left(\dfrac{1}{n} \right) + \dfrac{1}{n} </tex>
| |
− | '''return''' <tex>0</tex>
| |
− |
| |
− | <tex> \Longleftarrow </tex>
| |
− |
| |
− | : Построим функцию <tex> a(\varepsilon) </tex>
| |
− |
| |
− | '''function''' <tex> a(\varepsilon) </tex>:
| |
− | '''for''' <tex> x \in A </tex>
| |
− | '''if''' <tex> x + \varepsilon \notin A </tex>
| |
− | '''return''' <tex>x</tex>
| |
− |
| |
− | : Так как <tex> A </tex> разрешимо, <tex> \alpha = \sup A </tex> и для любого <tex> x \in A </tex> проверка в условном операторе завершается за конечное время, то функция <tex> a(\varepsilon) </tex> вычислима для любого рационального <tex> \varepsilon </tex>.
| |
− | }}
| |
− |
| |
− | '''Важное замечание''': построенное нами доказательство неконструктивно, так как мы не знаем наперед, рационально ли число <tex> \alpha </tex>, и уж тем более не пытаемся понять в случае его рациональности, чему именно оно равно. Но, так как мы ставим целью исследование свойств вычислимых чисел, а не явное построение соответствующих этим свойствам программ, то нам это доказательство полностью подходит.
| |
− |
| |
− | С учетом только что доказанной теоремы, далее при проверке на принадлежность числа <tex> x </tex> множеству <tex> A </tex> будем писать просто <tex> x < \alpha </tex>.
| |
− |
| |
− | {{Теорема
| |
− | |statement=
| |
− | Число <tex> \alpha </tex> вычислимо <tex>\iff</tex> последовательность знаков представляющей его двоичной записи вычислима.
| |
− | |proof=
| |
− | <tex> \Longrightarrow </tex>
| |
− |
| |
− | : Если число <tex> \alpha </tex> {{---}} рациональное, то необходимую последовательность можно получить, воспользовавшись стандартным алгоритмом перевода числа в двоичную систему счисления. Рассмотрим случай, когда <tex> \alpha \in \mathbb R \setminus \mathbb Q </tex>.
| |
− |
| |
− | : Очевидно, двоичная запись целой части <tex> \alpha </tex> всегда вычислима (так как множество чисел, меньших <tex> \alpha </tex>, разрешимо, то можно перебрать все целые числа в порядке возрастания их абсолютных величин и найти наибольшее число, меньшее <tex> \alpha </tex>), поэтому будем считать, что <tex> \alpha \in (0; 1) </tex>.
| |
− |
| |
− | : Напишем программу, которая по числу <tex> n </tex> вычисляет <tex> n </tex>-ный знак после запятой в двоичном представлении числа <tex> \alpha </tex> :
| |
− |
| |
− | '''function''' <tex>p(n)</tex>:
| |
− | <tex>l = 0, r = 1 </tex>
| |
− | '''for''' <tex> k = 1</tex> '''to''' <tex>n </tex>
| |
− | <tex> m = \dfrac{l+r}{2} </tex>
| |
− | '''if''' <tex> m < \alpha</tex>
| |
− | <tex>l = m, t = 1</tex>
| |
− | '''else'''
| |
− | <tex>r = m, t = 0</tex>
| |
− | '''return''' <tex>t</tex>
| |
− |
| |
− | <tex> \Longleftarrow </tex>
| |
− | : Для любого рационального <tex> \varepsilon > 0 </tex>, найдем <tex> n: 2^{-n} < \varepsilon </tex>, тогда в качестве значения функции <tex> a(\varepsilon) </tex> можно взять часть последовательности знаков двоичной записи <tex> \alpha </tex> с <tex> n </tex> знаками после запятой.
| |
− | }}
| |
− |
| |
− | {{Определение
| |
− | |definition=
| |
− | Последовательность рациональных чисел <tex> \{ r_n \} </tex> '''вычислимо сходится''' к <tex> \alpha </tex>, если существует вычислимая функция <tex> N: \mathbb Q \rightarrow \mathbb N </tex>, такая, что для любого рационального <tex> \varepsilon > 0 </tex> выполняется <tex>\forall n > N(\varepsilon): |r_n - \alpha| < \varepsilon </tex>.
| |
| }} | | }} |
| | | |
| {{Теорема | | {{Теорема |
| |statement= | | |statement= |
− | Число <tex> \alpha </tex> вычислимо <tex>\iff</tex> существует вычислимая последовательность рациональных чисел, вычислимо сходящаяся к <tex> \alpha </tex>. | + | Число <tex> \alpha </tex> вычислимо тогда и только тогда, когда существует вычислимая последовательность рациональных чисел, вычислимо сходящаяся к <tex> \alpha </tex>. {{TODO|t=вычислимо сходящаяся}} |
| |proof= | | |proof= |
− | <tex> \Longrightarrow </tex>
| |
− |
| |
− | : Так как <tex> \alpha </tex> вычислимо, то, согласно предыдущей теореме, вычислима и его двоичная запись. Пусть <tex> r_n </tex> {{---}} часть последовательности знаков двоичной записи <tex> \alpha </tex> с <tex> n </tex> знаками после запятой. Данная последовательность вычислима, а также вычислимо сходится к <tex> \alpha </tex>, так как <tex> N(\varepsilon) = \lceil -\log_2 \varepsilon \rceil</tex>.
| |
− |
| |
− | <tex> \Longleftarrow </tex>
| |
− |
| |
− | : Пусть <tex> a(\varepsilon) = r_{N(\varepsilon)} </tex>, тогда <tex> \alpha </tex> вычислимо по определению.
| |
| }} | | }} |
| | | |
| {{Теорема | | {{Теорема |
| |statement= | | |statement= |
− | | + | Пусть числа <tex> \alpha, \beta </tex> вычислимы. Тогда также вычислимы числа <tex> \alpha + \beta </tex>, <tex> \alpha - \beta </tex>, <tex> \alpha \beta </tex> и <tex> \frac \alpha \beta </tex>. |
− | Пусть числа <tex> \alpha, \beta </tex> вычислимы. Тогда также вычислимы числа <tex> \alpha + \beta </tex>, <tex> \alpha - \beta </tex>, <tex> \alpha \beta </tex> и <tex> \dfrac {\alpha} {\beta} </tex>. | |
| |proof= | | |proof= |
− | В пределах этого доказательства будем обозначать функцию-приближение для произвольного вычислимого числа <tex> x </tex> как <tex> f_x </tex>.
| |
− |
| |
− | Для того, чтобы получить приближение для результата операции, нам нужно выразить функцию-результат через приближения для операндов.
| |
− |
| |
− | Заметим, что <tex> |(\alpha \pm \beta) - (a \pm b)| \leqslant |\alpha - a| \pm |\beta - b| </tex>, для произвольных рациональных <tex> a, b </tex>, значит, в качестве необходимых функций для <tex> \alpha + \beta </tex> и <tex> \alpha - \beta </tex> можно взять
| |
− |
| |
− | <tex > f_{\alpha + \beta}(\varepsilon) = f_{\alpha}\left(\dfrac {\varepsilon} {2}\right) + f_{\beta}\left(\dfrac {\varepsilon} {2}\right) </tex>
| |
− |
| |
− | и
| |
− |
| |
− | <tex> f_{\alpha - \beta}(\varepsilon) = f_{\alpha}\left(\dfrac {\varepsilon} {2}\right) - f_{\beta}\left(\dfrac {\varepsilon} {2}\right) </tex>
| |
− |
| |
− | соответственно (при подстановке в неравенство <tex> f_{\alpha} </tex> и <tex> f_{\beta} </tex> вместо <tex> a </tex> и <tex> b </tex> каждый модуль в правой части не превосходит <tex> \dfrac {\varepsilon} {2} </tex>, поэтому <tex> f_{\alpha \pm \beta} </tex> не превосходит <tex> \varepsilon </tex>).
| |
− |
| |
− | Далее, так как
| |
− |
| |
− | <tex> |\alpha \beta - ab| = |(\alpha \beta - a \beta) + (a \beta - ab)| \leqslant |\beta||\alpha - a| + |a||\beta - b| \leqslant |b_\beta||\alpha - a| + |a||\beta - b|</tex>,
| |
− |
| |
− | где <tex>b_\beta</tex> {{---}} наименьшее рациональное число, большее <tex>\beta</tex> по модулю (т.е. <tex> b_\beta \in \mathbb Q, |b_\beta| > |\beta| </tex>), то
| |
− |
| |
− | <tex> f_{\alpha \beta}(\varepsilon) = f_{\alpha}\left(\dfrac {\varepsilon} {b_\beta}\right) f_{\beta}\left(\dfrac {\varepsilon} {a}\right) </tex>.
| |
− |
| |
− | Убедимся в вычислимости числа <tex> \dfrac {1} {\alpha} </tex>:
| |
− |
| |
− | <tex> \left|\dfrac {1} {\alpha} - \dfrac{1}{a}\right| \leqslant \dfrac {|a - \alpha|}{a a_\alpha} </tex>, где <tex> a_\alpha \in \mathbb Q, |a_\alpha| < |\alpha| </tex>.
| |
− |
| |
− | <tex> f_{\frac {1} {\alpha}}(\varepsilon) = f_{\alpha}(\varepsilon a a_\alpha) </tex>.
| |
− |
| |
− | Отсюда, <tex> \dfrac {\alpha} {\beta} = \dfrac{1} {\alpha} \beta </tex> также вычислимо.
| |
− |
| |
| }} | | }} |
| | | |
Строка 127: |
Строка 31: |
| Корень многочлена с вычислимыми коэффициентами вычислим. | | Корень многочлена с вычислимыми коэффициентами вычислим. |
| |proof= | | |proof= |
− | Пусть <tex> x </tex> {{---}} корень многочлена <tex> P(x) </tex> с вычислимыми коэффициентами.
| |
− |
| |
− | Если <tex> x \in \mathbb Q </tex>, то его можно найти точно, перебрав все рациональные числа.
| |
− |
| |
− | Иначе, выберем некоторый интервал <tex> [a; b]: x \in [a; b] </tex> (<tex> a, b </tex> {{---}} вычислимы), достаточно малый, чтобы полином <tex> P(x) </tex> был монотонным на отрезках <tex> [a; x] </tex> и <tex> [x; b] </tex>.
| |
− |
| |
− | Заметим, что для вычислимого <tex> t </tex> значение <tex> P(t) </tex> также вычислимо, так как в процессе его вычисления используются только операции, вычислимость значений которых уже была ранее доказана.
| |
− |
| |
− | Теперь, если полином имеет разные знаки на отрезках <tex> [a; x] </tex> и <tex> [x; b] </tex>, то для поиска <tex> x </tex> можно воспользоваться двоичным поиском нуля на <tex> [a; b] </tex>, иначе {{---}} троичным поиском экстремума на <tex> [a; b] </tex>.
| |
− |
| |
− | Останавливая тот или иной алгоритм, когда текущая длина интервала становится меньше <tex> \varepsilon </tex> и возвращая левую границу в качестве ответа, получаем функцию <tex> a_x(\varepsilon) </tex>.
| |
| }} | | }} |
| | | |
| {{Теорема | | {{Теорема |
| |statement= | | |statement= |
− | Предел вычислимо сходящейся вычислимой последовательности <tex> r_n </tex> вычислимых чисел вычислим. | + | Предел вычислимо сходящейся вычислимой последовательности вычислимых чисел вычислим. |
| |proof= | | |proof= |
− | Пусть <tex> \alpha = \lim\limits_{n \to \infty} r_n </tex>. Запишем формально данные нам условия:
| |
− |
| |
− | <tex> \forall \varepsilon_1 > 0\ \forall n > (\varepsilon_1): |r(n) - \alpha| < \varepsilon_1 </tex>
| |
− |
| |
− | <tex> \forall n\ \forall \varepsilon_2 > 0\ |r(n) - p_n(\varepsilon_2)| < \varepsilon_2 </tex>
| |
− |
| |
− | Здесь функции <tex> N(\varepsilon_1) </tex>, <tex> r(n) </tex> и все <tex> p_n(\varepsilon_2) </tex> вычислимы.
| |
− |
| |
− | Построим функцию <tex> a(\varepsilon) </tex>, которая дает приближение к <tex> \alpha </tex> с точностью до <tex> \varepsilon </tex>:
| |
− |
| |
− | '''function''' <tex> a(\varepsilon) </tex>:
| |
− | n = <tex> N\left(\dfrac {\varepsilon} {2}\right) + 1 </tex>
| |
− | '''return''' <tex> p_n\left(\dfrac {\varepsilon} {2}\right) </tex>
| |
− |
| |
− | Так как <tex> \left|\alpha - p_{n}\left(\dfrac {\varepsilon} {2}\right)\right| < |\alpha - r(n)| + \left|r(n) - p_n\left(\dfrac {\varepsilon} {2}\right)\right| </tex>, оба слагаемых меньше <tex> \dfrac {\varepsilon} {2} </tex> (первое {{---}} по выбору <tex> n </tex>, второе {{---}} в силу вычислимости <tex> p_n </tex>), то <tex> \left|\alpha - p_n\left(\dfrac \varepsilon 2\right)\right| < \varepsilon </tex>, и <tex> a(\varepsilon) </tex> действительно вычисляет требуемое приближение.
| |
| }} | | }} |
| | | |
Строка 165: |
Строка 43: |
| {{Определение | | {{Определение |
| |definition= | | |definition= |
− | Действительное число <tex> \alpha </tex> называется '''перечислимым снизу''' (англ. ''recursively enumerable number''), если множество <tex> \{ a \in \mathbb Q \mid a < \alpha \} </tex> [[Перечислимые языки|перечислимо]]. | + | Действительное число <tex> \alpha </tex> называется '''перечислимым снизу''' (recursively enumerable number), если множество <tex> \{ a \in Q \mid a < \alpha \} </tex> перечислимо. |
| }} | | }} |
| | | |
− | {{Определение
| + | Аналогично определяются перечислимые сверху числа. |
− | |definition=
| |
− | Действительное число <tex> \alpha </tex> называется '''перечислимым сверху''', если множество <tex> \{ a \in \mathbb Q \mid a > \alpha \} </tex> перечислимо.
| |
− | }}
| |
| | | |
| === Свойства === | | === Свойства === |
Строка 177: |
Строка 52: |
| {{Теорема | | {{Теорема |
| |statement= | | |statement= |
− | Число <tex> \alpha </tex> перечислимо снизу <tex>\iff</tex> существует вычислимая возрастающая последовательность рациональных чисел, пределом которой является <tex> \alpha </tex>. | + | Число <tex> \alpha </tex> перечислимо снизу тогда и только тогда, когда существует вычислимая возрастающая последовательность рациональных чисел, пределом которой является <tex> \alpha </tex>. |
| |proof= | | |proof= |
− | <tex>\Longrightarrow</tex>
| |
− |
| |
− | : По определению <tex> \alpha </tex>, множество <tex> A = \{ a \in \mathbb Q \mid a < \alpha \} </tex> перечислимо. Кроме того, <tex> \sup A = \alpha </tex>.
| |
− |
| |
− | : По определению нижней грани, <tex> \forall \varepsilon > 0\ \exists x_\varepsilon \in A: \varepsilon > \alpha - x_\varepsilon </tex>. Тогда можно взять, например, последовательность <tex> a_n = x_{\frac{1}{n}} </tex>.
| |
− |
| |
− | <tex>\Longleftarrow</tex>
| |
− |
| |
− | : Построим полуразрешитель для множества <tex> A </tex>:
| |
− |
| |
− | '''function''' <tex>p(x)</tex>:
| |
− | '''for''' n = <tex>1</tex> '''to''' <tex>\infty</tex>
| |
− | '''if''' <tex> x < a_n </tex>
| |
− | '''return''' <tex>1</tex>
| |
− |
| |
− | : Если <tex> x \in A</tex>, то <tex> \alpha - x = t > 0 </tex>, и так как <tex> \exists N:\ \forall n > N |a_n - \alpha| < t </tex>, то программа вернет ответ при <tex> n > N </tex>.
| |
| }} | | }} |
| | | |
| {{Теорема | | {{Теорема |
| |statement= | | |statement= |
− | Число <tex> \alpha </tex> вычислимо <tex>\iff</tex> оно перечислимо сверху и снизу. | + | Число <tex> \alpha </tex> вычислимо тогда и только тогда, когда оно перечислимо сверху и снизу. |
− | |proof=
| |
− | Обозначим множества <tex> \{a \in \mathbb Q \mid a < \alpha \} </tex> и <tex> \{a \in \mathbb Q \mid a > \alpha \} </tex> за <tex> A </tex> и <tex> B </tex> соответственно.
| |
− | | |
− | Если <tex> \alpha </tex> рационально, то необходимые (полу)разрешители строятся тривиально.
| |
− | | |
− | В противном случае, так как <tex> B = \mathbb Q \setminus A</tex>, то перечислимость множеств <tex> A </tex> и <tex> B </tex> равносильна разрешимости множества <tex> A </tex>, которая, в свою очередь, равносильна вычислимости <tex> \alpha </tex>.
| |
| }} | | }} |
| | | |
| == Последовательность Шпеккера == | | == Последовательность Шпеккера == |
− |
| |
− | Множество всех программ счётно, поэтому множество вычислимых чисел также счётно. Однако, множество вещественных чисел несчётно, значит, существуют невычислимые вещественные числа. Построим явно пример такого числа.
| |
| | | |
| {{Определение | | {{Определение |
| |definition= | | |definition= |
− | Пусть <tex> A </tex> {{---}} некоторое перечислимое, но неразрешимое множество натуральных чисел. Пронумеруем его элементы. '''Последовательностью Шпеккера''' <tex> \{q_n\} </tex> называется последовательность рациональных чисел, <tex>n</tex>-ный член которой определяется как <tex> q_n = \sum\limits_{k=1}^{n} 2^{-A(n)-1} </tex>. | + | Пусть <tex> A </tex> {{---}} некоторое перечислимое, но неразрешимое множество натуральных чисел. Занумеруем его элементы. '''Последовательностью Шпеккера''' <tex> \{q_n\} </tex> называется последовательность рациональных чисел, <tex>n</tex>-ный член которой определяется как <tex> q_n = \sum\limits_{k=1}^{n} 2^{-A(n)-1} </tex>. |
| }} | | }} |
| | | |
− | Данная последовательность строго возрастает и ограничена числом <tex> 1 </tex>, следовательно, она сходится по признаку Вейерштрасса. | + | Данная последовательность строго возрастает и ограничена числом <tex> 1 </tex>, следовательно, по признаку Вейерштрасса, она сходится. |
| | | |
| {{Теорема | | {{Теорема |
| |statement= | | |statement= |
| Число <tex> q = \lim\limits_{n \to \infty} q_n </tex> перечислимо снизу, но невычислимо. | | Число <tex> q = \lim\limits_{n \to \infty} q_n </tex> перечислимо снизу, но невычислимо. |
− | |proof=
| |
− | <tex> q </tex> перечислимо снизу, как предел возрастающей вычислимой последовательности рациональных чисел.
| |
− |
| |
− | Допустим теперь, что <tex> q </tex> {{---}} вычислимо.
| |
− |
| |
− | Пусть <tex> A = \{p \mid p(p) = 1\}</tex>. Рассмотрим двоичную запись числа <tex> q </tex>, если ее <tex> n </tex>-ный знак после запятой равен 1, то <tex> n \in A </tex>, иначе {{---}} <tex> n \notin A </tex>. Мы построили разрешитель для множества <tex> A </tex>. Тем не менее, известно, что <tex> A </tex> {{---}} неразрешимое множество, а это невозможно, значит, <tex> q </tex> {{---}} невычислимо.
| |
| }} | | }} |
| | | |
− | == См. также == | + | == Ссылки == |
− | * [[Рекурсивные функции]]
| + | * ''Верещагин Н. К., Шень А.'' Лекции по математической логике и теории алгоритмов. Часть 3. Вычислимые функции — М.: МЦНМО, 1999 - С. 14 |
− | * [[Вычислимые функции]]
| + | * http://en.wikipedia.org/wiki/Computable_number |
− | | + | * http://en.wikipedia.org/wiki/Specker_sequence |
− | == Источники информации ==
| |
− | * ''Верещагин Н. К., Шень А.'' {{---}} Лекции по математической логике и теории алгоритмов. Часть 3. Вычислимые функции — М.: МЦНМО, 1999 {{---}} стр. 14 | |
− | * [http://en.wikipedia.org/wiki/Computable_number Computable number] | |
− | * [http://en.wikipedia.org/wiki/Specker_sequence Specker sequence] | |
| | | |
| + | [[Категория: Теория вычислимости]] |
| [[Категория: Теория формальных языков]] | | [[Категория: Теория формальных языков]] |
− | [[Категория: Теория вычислимости]]
| |
− | [[Категория: Разрешимые и перечислимые языки]]
| |