Продолжим?
Вот мы, применяя алгоритм Евклида, получили последовательность P_j, начинающуюся с наших P_1=P и P_2=P' (ну или абы каких P_1 и P_2, если мы ищем топологическую степень произвольной рациональной функции R=P_1/P_2). И её можно определить вот такой цепочкой:
P_1+P_3 = P_2 * Q_1,
P_2+P_4 = P_3 * Q_2,
P_3+P_5 = P_4 * Q_3,
...
и заканчивается она на каком-то P_k=const, и P_{k+1}=0.
Давайте для каждой точки x рассмотрим последовательность значений P_1(x), P_2(x), ..., P_k(x),
и обозначим через N(x) число перемен знака в этой последовательности: сколько раз + меняется на - и наоборот.
Классический способ формулировать теорему Штурма — как раз в терминах числа перемен знака:
Теорема Штурма. Пусть P — многочлен без кратных корней, не обращающийся в ноль в точках a и b. Тогда число корней P на отрезке [a,b] равно N(a)-N(b).
(А полное число вещественных корней это N(-\infty)-N(\infty), благо что знак полинома на плюс-минус бесконечности определить не проблема.)
Пример: возьмём P(x)=x^2-1. Тогда P'(x)=2x, и мы получаем последовательность
P_1=x^2-1,
P_2=2x,
P_3=1 (помните про смену знака у остатка!).
Последовательность знаков в разных точках и число перемен знака:
-\infty: +,-,+; 2 перемены;
-2: +,-,+; 2 перемены;
-0.5: -,-,+; 1 перемена;
-0.1: -,-,+; 1 перемена;
0.1: -,+,+; 1 перемена;
1.5: +,+,+; 0 перемен;
\infty: +,+,+; 0 перемен.
Уже из бесконечностей мы видим, что число перемен знака меняется между +\infty и -\infty на 2, так что у многочлена x^2-1 два вещественных корня. Более того, такое же число перемен знака между -2 и 1.5, так что они живут на этом отрезке. Причём один из них живёт на отрезке [-2,-0.5], а второй на отрезке [0.1,1.5]. Большой прогресс!
Но если серьёзно — мы только что научились алгоритмически (и довольно просто!) искать (с любой нужной точностью) все корни вещественного многочлена. Начинаем с какого-нибудь отрезка, на котором они гарантированно все лежат (за пределами которого старший член a_n*x^n превышает по модулю сумму модулей всех остальных), а потом запускаем деление пополам, последовательно деля каждый отрезок, на котором число перемен знака меняется.
Осталось
а) доказать;
б) понять, как этот подход "склеивается" с предыдущим.
Вот мы, применяя алгоритм Евклида, получили последовательность P_j, начинающуюся с наших P_1=P и P_2=P' (ну или абы каких P_1 и P_2, если мы ищем топологическую степень произвольной рациональной функции R=P_1/P_2). И её можно определить вот такой цепочкой:
P_1+P_3 = P_2 * Q_1,
P_2+P_4 = P_3 * Q_2,
P_3+P_5 = P_4 * Q_3,
...
и заканчивается она на каком-то P_k=const, и P_{k+1}=0.
Давайте для каждой точки x рассмотрим последовательность значений P_1(x), P_2(x), ..., P_k(x),
и обозначим через N(x) число перемен знака в этой последовательности: сколько раз + меняется на - и наоборот.
Классический способ формулировать теорему Штурма — как раз в терминах числа перемен знака:
Теорема Штурма. Пусть P — многочлен без кратных корней, не обращающийся в ноль в точках a и b. Тогда число корней P на отрезке [a,b] равно N(a)-N(b).
(А полное число вещественных корней это N(-\infty)-N(\infty), благо что знак полинома на плюс-минус бесконечности определить не проблема.)
Пример: возьмём P(x)=x^2-1. Тогда P'(x)=2x, и мы получаем последовательность
P_1=x^2-1,
P_2=2x,
P_3=1 (помните про смену знака у остатка!).
Последовательность знаков в разных точках и число перемен знака:
-\infty: +,-,+; 2 перемены;
-2: +,-,+; 2 перемены;
-0.5: -,-,+; 1 перемена;
-0.1: -,-,+; 1 перемена;
0.1: -,+,+; 1 перемена;
1.5: +,+,+; 0 перемен;
\infty: +,+,+; 0 перемен.
Уже из бесконечностей мы видим, что число перемен знака меняется между +\infty и -\infty на 2, так что у многочлена x^2-1 два вещественных корня. Более того, такое же число перемен знака между -2 и 1.5, так что они живут на этом отрезке. Причём один из них живёт на отрезке [-2,-0.5], а второй на отрезке [0.1,1.5]. Большой прогресс!
Но если серьёзно — мы только что научились алгоритмически (и довольно просто!) искать (с любой нужной точностью) все корни вещественного многочлена. Начинаем с какого-нибудь отрезка, на котором они гарантированно все лежат (за пределами которого старший член a_n*x^n превышает по модулю сумму модулей всех остальных), а потом запускаем деление пополам, последовательно деля каждый отрезок, на котором число перемен знака меняется.
Осталось
а) доказать;
б) понять, как этот подход "склеивается" с предыдущим.
К прошедшему вчера Pi Day хочется вспомнить несколько видео: три видео 3Blue1Brown,
*) https://www.youtube.com/watch?v=d-o3eB9sfls — очень классное доказательство формулы для суммы обратных квадратов, где окружность с расставленными на ней маяками появляется явно;
*) https://www.youtube.com/watch?v=8GPy_UMV-08 — а это формула Валлиса (произведение 2/1 * 2/3 * 4/3 * 4/5 * ... = π/2); опять окружность, выход в комплексные числа, многочлены, произведение расстояний — и всё получается;
*) https://www.youtube.com/watch?v=NaL_Cb42WyY — как связать сумму
1 - 1/3 + 1/5 - 1/7 + ... = π/4
с окружностью (спойлер: никакого Тейлора и арктангенса! Разложение чисел в сумму квадратов, характеры, подсчёт числа целых точек).
И рождественскую лекцию Дональда Кнута в Стенфорде, "Why Pi":
https://www.youtube.com/watch?v=cI6tt9QfRdo (которая, на самом деле, заслуживает отдельного пересказа).
*) https://www.youtube.com/watch?v=d-o3eB9sfls — очень классное доказательство формулы для суммы обратных квадратов, где окружность с расставленными на ней маяками появляется явно;
*) https://www.youtube.com/watch?v=8GPy_UMV-08 — а это формула Валлиса (произведение 2/1 * 2/3 * 4/3 * 4/5 * ... = π/2); опять окружность, выход в комплексные числа, многочлены, произведение расстояний — и всё получается;
*) https://www.youtube.com/watch?v=NaL_Cb42WyY — как связать сумму
1 - 1/3 + 1/5 - 1/7 + ... = π/4
с окружностью (спойлер: никакого Тейлора и арктангенса! Разложение чисел в сумму квадратов, характеры, подсчёт числа целых точек).
И рождественскую лекцию Дональда Кнута в Стенфорде, "Why Pi":
https://www.youtube.com/watch?v=cI6tt9QfRdo (которая, на самом деле, заслуживает отдельного пересказа).
YouTube
Why is pi here? And why is it squared? A geometric answer to the Basel problem
A most beautiful proof of the Basel problem, using light.
Help fund future projects: https://www.patreon.com/3blue1brown
An equally valuable form of support is to simply share some of the videos.
Special thanks to these supporters: http://3b1b.co/basel-thanks…
Help fund future projects: https://www.patreon.com/3blue1brown
An equally valuable form of support is to simply share some of the videos.
Special thanks to these supporters: http://3b1b.co/basel-thanks…
Forwarded from Непрерывное математическое образование
а про формулу Валлиса для π и про появление π в комбинаторике лет 10 назад интересно рассказывал Д.Кнут — можно посмотреть видео или почитать конспект
YouTube
Stanford Lecture: Donald Knuth—"Why Pi?"(2010)
Don Knuth's 16th Annual Christmas Tree Lecture
December 6th, 2010
Professor Donald Knuth discusses recent discoveries that have uncovered a fascinating relationship between circles and the theory of trees.
Learn more: http://scpd.stanford.edu/knuth/index.jsp
December 6th, 2010
Professor Donald Knuth discusses recent discoveries that have uncovered a fascinating relationship between circles and the theory of trees.
Learn more: http://scpd.stanford.edu/knuth/index.jsp
Forwarded from Непрерывное математическое образование
если хочется объяснения и/или еще формул — можно заглянуть в калейдоскоп
Математические байки
Продолжим? Вот мы, применяя алгоритм Евклида, получили последовательность P_j, начинающуюся с наших P_1=P и P_2=P' (ну или абы каких P_1 и P_2, если мы ищем топологическую степень произвольной рациональной функции R=P_1/P_2). И её можно определить вот такой…
Посмотрим сначала на классическое доказательство теоремы Штурма. Если мы говорим о том, что число корней P на отрезке [a,b] выражается через изменение на этом отрезке числа перемен знака — так очень естественно спросить себя, а как вообще число перемен знака N(x) меняется, когда меняется x?
Естественно, что меняться оно может только если меняет знак какой-то из многочленов P_j. Но — и это замечательно — если меняет знак (переходя через ноль) любой P_j, кроме P_1, то число перемен знака не меняется.
Действительно, пусть в точке x_0 обращается в ноль какой-то P_j, и j>1. Его соседи тогда в этой точке не обращаются в ноль: иначе там обратились бы в ноль и все полиномы вообще, а у P_1 и P_2 нет общих корней.
Но заметим, что его соседи противоположного знака: ведь
P_{j-1} + P_{j+1} = P_j * Q_{j-1},
раз правая часть обратилась в ноль, то и сумма соседей P_{j-1}+P_{j+1} в точке x_0 обращается в ноль.
А тогда неважно, какой у P_j рядом с x_0 знак — всё равно на P_{j-1},P_j,P_{j+1} приходится ровно одна перемена знака!
И мы, собственно, такое уже видели, когда смотрели на пример P(x)=x^2: рядом с x_0=0, где обращается в ноль производная P_2=2x, число перемен знака не менялось:
===
-0.1: -,-,+; 1 перемена;
0.1: -,+,+; 1 перемена;
===
Естественно, что меняться оно может только если меняет знак какой-то из многочленов P_j. Но — и это замечательно — если меняет знак (переходя через ноль) любой P_j, кроме P_1, то число перемен знака не меняется.
Действительно, пусть в точке x_0 обращается в ноль какой-то P_j, и j>1. Его соседи тогда в этой точке не обращаются в ноль: иначе там обратились бы в ноль и все полиномы вообще, а у P_1 и P_2 нет общих корней.
Но заметим, что его соседи противоположного знака: ведь
P_{j-1} + P_{j+1} = P_j * Q_{j-1},
раз правая часть обратилась в ноль, то и сумма соседей P_{j-1}+P_{j+1} в точке x_0 обращается в ноль.
А тогда неважно, какой у P_j рядом с x_0 знак — всё равно на P_{j-1},P_j,P_{j+1} приходится ровно одна перемена знака!
И мы, собственно, такое уже видели, когда смотрели на пример P(x)=x^2: рядом с x_0=0, где обращается в ноль производная P_2=2x, число перемен знака не менялось:
===
-0.1: -,-,+; 1 перемена;
0.1: -,+,+; 1 перемена;
===
Так что переходы значений "внутренних" полиномов через ноль на число перемен знака не влияют. P_n=1 через ноль перейти вообще не может, это константа. Остаётся только P_1.
А значение многочлена P переходит через ноль — правее корня знак у P такой же, как у его производной P' в этом корне, а левее корня противоположный. Поэтому как раз при переходе через ноль P (и только там) число перемен знака уменьшается на единицу. Теорема Штурма доказана.
Осталось это склеить с тем взглядом, который у нас был раньше.
А значение многочлена P переходит через ноль — правее корня знак у P такой же, как у его производной P' в этом корне, а левее корня противоположный. Поэтому как раз при переходе через ноль P (и только там) число перемен знака уменьшается на единицу. Теорема Штурма доказана.
Осталось это склеить с тем взглядом, который у нас был раньше.
Да, если P_1 и P_2 произвольные без общих корней, то большая часть рассуждения выше тоже работает. Только разница в числе перемен знака на отрезке [a,b] считает не просто корни P_1 — а корни P_1 со знаком: если в корне P_1' и P_2 одного знака, то мы его считаем с плюсом, а если разного, то с минусом. (Потому что именно так пересечение корня P_1 меняет число перемен знака — а пересечение корней всех остальных P_j, как и раньше, ничего не меняет.)
Математические байки
Да, если P_1 и P_2 произвольные без общих корней, то большая часть рассуждения выше тоже работает. Только разница в числе перемен знака на отрезке [a,b] считает не просто корни P_1 — а корни P_1 со знаком: если в корне P_1' и P_2 одного знака, то мы его считаем…
Собственно, это буквально соответствует вычислению части суммы для топологической степени deg_{top} P_1/P_2 по прообразу нуля (то есть по пересечению с вертикальной прямой на плоскости (P_1,P_2)), когда мы ограничиваемся прообразами из отрезка параметров [a,b].
Математические байки
Photo
Понять, что выкладки через цепную дробь приводят к тому же результату, не очень сложно. Действительно, когда мы делим с остатком, то знаки P_j и P_{j+1} отличаются на знак неполного частного Q_j (а P_{j+1} слишком мал, и ни на что не влияет). Соответственно, если многочлен Q_j в данной (плюс или минус) бесконечности отрицателен, то перемена знака между P_j и P_{j+1} там есть, а если положителен, то нет. Вот и получается, что когда мы считаем разницу между числом перемен знака в минус и в плюс бесконечности, то Q_j чётной степени вклада не дают (ибо они в бесконечностях одного знака), а Q_j нечётной степени дают вклад, равный знаку старшего коэффициента.
Так что ответы "через цепную дробь" и "через перемены знаков" для числа всех вещественных нулей (и для топологической степени рациональной функции R=P_1/P_2) совпадают. Уже хорошо.
Но хотелось бы как-то понять, как о переменах знака думать в топологических терминах — при том, что для отображения отрезка о степени уже говорить нельзя, у него концы есть.
Так что ответы "через цепную дробь" и "через перемены знаков" для числа всех вещественных нулей (и для топологической степени рациональной функции R=P_1/P_2) совпадают. Уже хорошо.
Но хотелось бы как-то понять, как о переменах знака думать в топологических терминах — при том, что для отображения отрезка о степени уже говорить нельзя, у него концы есть.
На самом деле, история довольно простая — и параллельная тому, что мы уже делали. А именно — разложение рациональной функции
R=P_1/P_2
в цепную дробь Хирцебруха ("с минусами") разбивается в череду операций:
*) повернуть картинку на 90 градусов (P_1,P_2) -> (P_2,-P_1)
*) добавить к второй координате первую, умноженную на Q_1
(P_2,-P_1) -> (P_2, Q_1*P_2-P_1) = (P_2,P_3)
*) повернуть картинку на 90 градусов (P_2,P_3) -> (P_3,-P_2)
*) добавить к второй координате первую, умноженную на Q_2
(P_3,-P_2) -> (P_3, Q_2*P_3-P_2) = (P_3,P_4)
*) и так далее.
R=P_1/P_2
в цепную дробь Хирцебруха ("с минусами") разбивается в череду операций:
*) повернуть картинку на 90 градусов (P_1,P_2) -> (P_2,-P_1)
*) добавить к второй координате первую, умноженную на Q_1
(P_2,-P_1) -> (P_2, Q_1*P_2-P_1) = (P_2,P_3)
*) повернуть картинку на 90 градусов (P_2,P_3) -> (P_3,-P_2)
*) добавить к второй координате первую, умноженную на Q_2
(P_3,-P_2) -> (P_3, Q_2*P_3-P_2) = (P_3,P_4)
*) и так далее.
Математические байки
И добавили ко второй координате первую, умноженную на Q_1 — из (P_2,-P_1) получив (P_2,P_3).
После чего продолжили в том же духе. Так вот — на операциях "добавить" мы пересечения с осью ординат не меняем (потому что добавляем первую координату, на что-то там умноженную). И их тип (по или против часовой стрелки) тоже.
Математические байки
Photo
А вот на операциях "повернуть" — можем поменять. Потому что сначала мы считали (с учётом знака) пересечения с одной прямой, а потом с другой. Если бы мы работали с замкнутой кривой, то мы бы считали буквально степень отображения окружности, и ничего поменяться бы не могло, потому что при повороте прямой прообразы бы только сливались "плюс с минусом", и ничего бы не менялось. Но если мы работаем с отрезком [a,b] параметров — то у отрезка кривой есть концы. И вот тут-то мы с числом перемен знака и столкнёмся!
Уже звучит очень естественно, что важно, пересекут ли в процессе поворота (кривой или системы координат — кому как удобнее) концы кривой ось ординат, пересечения с которой мы считаем. Но проще всего в этом убедиться, уронив всю ситуацию обратно на проективную прямую (= окружность, образуемую всеми прямыми через ноль).
Математические байки
Представьте себе теперь, что вдоль круглого стадиона бежит атлет. А мы должны сказать, сколько кругов он сделал. Понятно, что бежать за ним весьма трудозатратно; гораздо проще встать в одной точке и считать, сколько раз он мимо нас пробежал. Но считать надо…
Если по окружности-стадиону бежит атлет и возвращается в итоге в начальную точку — число кругов, которые он пробежал, можно посчитать, встав в любую точку стадиона и считая пересечения (если он пробежал куда надо, то с плюсом, а если обратно, то с минусом).
А если траектория атлета незамкнута — то арбитры, стоящие в разных точках круга, могут получить разные суммы. Но не сильно: если атлет начал и закончил бежать на одной и той же дуге, то их результаты совпадут — потому что он после этого может пройтись от конца своего пути к началу, не проходя мимо арбитров, и получить уже замкнутый путь.
А если траектория атлета незамкнута — то арбитры, стоящие в разных точках круга, могут получить разные суммы. Но не сильно: если атлет начал и закончил бежать на одной и той же дуге, то их результаты совпадут — потому что он после этого может пройтись от конца своего пути к началу, не проходя мимо арбитров, и получить уже замкнутый путь.
А если атлет начинал и закончил путь на разных дугах — то если добавить к пути замыкание "в положительном направлении", после которого суммы у арбитров совпадут, то один из двух арбитров его увидит, а другой нет. И уже совсем легко увидеть, что общий ответ такой:
число пересечений (с учётом знака) для арбитра A1 минус
число пересечений (с учётом знака) для арбитра А2 =
F(начала пути)-F(конца пути),
где F это 1, если мы начали на дуге от A2 к A1 (обозначенной тут +-), и 0 иначе.
число пересечений (с учётом знака) для арбитра A1 минус
число пересечений (с учётом знака) для арбитра А2 =
F(начала пути)-F(конца пути),
где F это 1, если мы начали на дуге от A2 к A1 (обозначенной тут +-), и 0 иначе.