Опять МФТИ

zykov
Сообщений: 1393
Зарегистрирован: 06 янв 2016, 17:41

Опять МФТИ

Сообщение zykov » 23 дек 2019, 09:56

Что-то много понаписали. Путаница вышла. Попробую кратко подвести резюме.
(У меня через спектр, но конечно есть и другие пути.)

$$|Ax| \to \min, \; |x|=1$$
У нас оператор $A$ легко диагонализируется (оператор дискретной или обычной производной имеет диагональный вид в Фурье базисе).
По теореме Min-max, минимум нормы будет для собственного вектора с минимальным по абсолютной величине собственным значением.
(Без других ограничений это была бы константа для которой нулевое собственное значение.)

Единственное осложнение, это дополнительное ограничение. Для дискретного $x_0=x_n=0$, для непрерывного $x(0)=x(1)=0$.
Ограничение линейное - $x$ должен лежать в подпространстве. Можно рассмотреть это подпространство уже в новом базисе (где оператор диагонален), что уже проще, но всё равно несколько сложно.

В данном случае простой фокус позволяет значительно упростить решение.
Сначала дополним вектор/функцию нечётным образом. Для дискретного $x_{-k}=-x_k$ при $k=0..n$, для непрерывного $x(-t)=-x(t)$ при $t \in [-1,1]$. Потом дополним везде периодическим образом. Для дискретного $x_{k+2n}=x_k$, для непрерывного $x(t+2)=x(t)$.
Ограничения $x_0=0$ и $x(0)=0$ следуют из нечётности. Ограничения $x_{-n}=x_n=0$ и $x(-1)=x(1)=0$ следуют из нечётности и периодичности (например, $x(-1)=x(1)$ из периодичности, $x(-1)=-x(1)$ из нечётности, значит оба нули).
Т.е. дополнительных ограничений нет, кроме периодичности и нечётности.
Периодичность даёт дискретный спектр, т.е. ряд Фурье (плюс, в дискретном случае ещё сам ряд периодичен, т.е. имеет всего $2n$ измерений).
Нечётность при переходе в Фурье базис тоже имеет простую интерпретацию - амплитуды косинусов нулевые, остаются только синусы.
Теперь уже константа с её нулевым собственным значением запрещена нечётностю.
Следующие собственные вектора с минимальными (по абсолютной величине) собственными значениями будут $e^{\pm i\pi k/n}$ (для непрерывного $e^{\pm i\pi t}$). Они комбинируются в $\sin(\pi k/n)$ (для непрерывного в $\sin(\pi t)$).

peregoudov
Сообщений: 620
Зарегистрирован: 29 дек 2015, 13:17

Опять МФТИ

Сообщение peregoudov » 23 дек 2019, 11:48

zykov писал(а):Source of the post И как этот минимум ищется?
Напрямую, поиском минимального собственного значения. Спектр в яме известен.

zykov писал(а):Source of the post Замечание: там $A$, это наше $A^TA$.
Ну то есть ровно та матрица, что я и рассматривал?

Тут ведь, как в непрерывном случае, принципиальны граничные условия. Какой спектр имеет оператор $d^2\!/dx^2$ на [0,1]? Да хз, от граничных условий зависит. Так и в дискретной цепочке: собственные частоты зависят от граничных условий. Выписать решение для бесконечной цепочки --- не штука, и циркулянты тут мало что упрощают. И ниоткуда вы не получите этот несчастный $\sin(\pi/2n)$, кроме как честно решая уравнения для конечной цепочки с учетом граничных условий.

zykov
Сообщений: 1393
Зарегистрирован: 06 янв 2016, 17:41

Опять МФТИ

Сообщение zykov » 23 дек 2019, 14:07

peregoudov писал(а):Source of the post И ниоткуда вы не получите этот несчастный $\sin(\pi/2n)$,

Эээ... Нет слов.
В предыдущем посте всё чётко расписано, что и откуда получается.

peregoudov
Сообщений: 620
Зарегистрирован: 29 дек 2015, 13:17

Опять МФТИ

Сообщение peregoudov » 23 дек 2019, 20:19

Да я понимаю, как вы рассуждаете. Вы фактически определяете минимальное собственное число трехдиагональной теплицевой матрицы таким вот вычурным методом. Просто мне кажется, что стандартный дуболомный подход ничуть не сложнее, а в непрерывном случае --- так даже явно проще.

zykov
Сообщений: 1393
Зарегистрирован: 06 янв 2016, 17:41

Опять МФТИ

Сообщение zykov » 23 дек 2019, 21:13

peregoudov писал(а):Source of the post Просто мне кажется, что стандартный дуболомный подход ничуть не сложнее, а в непрерывном случае --- так даже явно проще.

Стандартный подход "в лоб" потребует бумажки или компьютера.
Здесь же "устный счёт", всё в уме решается.

Аватар пользователя
Ian
Сообщений: 960
Зарегистрирован: 18 янв 2016, 19:42

Опять МФТИ

Сообщение Ian » 25 дек 2019, 10:33

Ну у меня тоже близко вокруг этой матрицы.
Зная, что вариационную задачу иначе как уравнением Лагранжа мы решать не умеем, в дискретном случае вместо вариаций будут малые отклонения этих [math] по отдельности.Придется найти все локальные экстремумы, интересно же.
Достаточно доказывать неравенство при ограничении[math]
[math]

Единичная сфера в [math]компактна, существует точка, где достигается минимум F>0. В этой точке существует градиент F (потому что он везде существует) и должен быть направлен наружу сферы перпендикулярно ей, в крайнем случае равен 0. Значит (уточнение правила множителей Лагранжа), существует [math], что [math]
[math]

Рекуррентность [math]
Так как [math] , характеристические корни могут быть либо 2 действительных положительных [math], либо комплексные [math]В первом случае [math](чтобы [math]) единственный вариант [math],t=0;2. Мы увидим что в этих случаях локальные максимумы на сфере.

Пусть [math]
[math]
[math]
[math]-это точки локальных минимумов при разных m. Найдем [math] из ограничения[math][math]

Сравним в них значения [math]
[math]

При [math] сумма меньше, так как [math] не меняют знаки

Значение суммы можно найти тем же приемом понижения степени, должно получиться [math]

guryev
Сообщений: 3
Зарегистрирован: 25 дек 2019, 15:04

Опять МФТИ

Сообщение guryev » 25 дек 2019, 16:28

Возвращаясь к первой задаче - можно применить индукцию.

Очевидно, что если точек всего две, в единичном шаре они находятся на расстоянии $$d_1 \leqslant \frac{2}{\sqrt{1}}$$, и их центры масс, естественно, тоже.

Теперь, пусть есть два непересекающиеся множества, $P$ и $Q$, состоящие каждое из $n$ точек единичного шара и расстояние между их центрами масс

$$d_n \leqslant \frac{2}{\sqrt{n}}$$

Обозначим вектор центра масс множества $P$ как $\mathbf{P}$, а множества $Q$ - через $\mathbf{Q}$. Очевидно,

$$d_n^2=(\mathbf{P}-\mathbf{Q})^2$$

Берём две новые точки из шара и добавляем одну из них в множество $P$, а другую - в $Q$. Больше ничего не меняем, вдруг получится.

Итак, для определённости, добавляем точку $\mathbf{a}$ в $P$, а точку $\mathbf{b}$ в $Q$. Теперь квадрат расстояния между центрами масс новых множеств будет

$$d_{n+1}^2=\left(\frac{n\mathbf{P}+\mathbf{a}}{n+1}-\frac{n\mathbf{Q}+\mathbf{b}}{n+1}\right)^2$$

То есть,

$$d_{n+1}^2=\frac{1}{(n+1)^2}(n^2(\mathbf{P}-\mathbf{Q})^2+2n(\mathbf{P}-\mathbf{Q})\cdot(\mathbf{a}-\mathbf{b})+(\mathbf{a}-\mathbf{b})^2)$$

Найдём максимальные значения для трёх слагаемых второго множителя. Для первого слагаемого это, как следует из гипотезы индукции, $4n$. Для третьего это $4$, поскольку $\mathbf{a}$ и $\mathbf{b}$ лежат в одном единичном шаре.

А для второго слагаемого можно считать, что максимум - ноль. Почему? Потому, что если оно окажется положительным, можно повторить вышеупомянутые рассуждения, только на этот раз поместить $\mathbf{b}$ в $P$, а $\mathbf{a}$ в $Q$. Всё будет точно так же, только $\mathbf{a}$ и $\mathbf{b}$ поменяются местами, и соответственно, у второго слагаемого поменяется знак. Таким образом, не теряя общности, можно считать, что второе слагаемое меньше или равно нулю. Надо только правильно выбрать, куда помещать $\mathbf{a}$, а куда $\mathbf{b}$.

Теперь, после несложных выкладок можно получить $$d_{n+1}^2 \leqslant \frac{4}{n+1}$$, т.е., $$d_{n+1} \leqslant \frac{2}{\sqrt{n+1}}$$.

zykov
Сообщений: 1393
Зарегистрирован: 06 янв 2016, 17:41

Опять МФТИ

Сообщение zykov » 25 дек 2019, 18:04

guryev писал(а):Source of the post Возвращаясь к первой задаче - можно применить индукцию.

Наверно авторы задачи это и имели ввиду.
У peregoudov получилась полее строгая оценка, но и метод сложнее.

Аватар пользователя
Ian
Сообщений: 960
Зарегистрирован: 18 янв 2016, 19:42

Опять МФТИ

Сообщение Ian » 25 дек 2019, 23:54

guryev писал(а):Возвращаясь к первой задаче
, Вы ее решили просто.А Вы случайно не участвовали в самой олимпиаде?
А для второго слагаемого можно считать, что максимум - ноль.
Как я и предсказывал в посте 1, худший для такой оценки случай, когда a и b лежат в подпространстве, перпендикулярном гиперплоскости, содержащей остальные 2m-2 точки, тогда ровно ноль.

zykov
Сообщений: 1393
Зарегистрирован: 06 янв 2016, 17:41

Опять МФТИ

Сообщение zykov » 26 дек 2019, 11:53

zykov писал(а):Source of the post У peregoudov получилась полее строгая оценка, но и метод сложнее.

Что интересно, при размерности от $2m-1$ и выше это самая строгая оценка.
Экстремальный пример - правильный $2m-1$ мерный симплекс (имеет $2m$ вершин) вписанный в $2m-1$ мерный шар единичного радиуса. Тут, в силу симметрии, как ни разбивай - будет одинаковое расстояние. Значит минимум равен максимуму и равен среднему.

Вот интересно, какая наиболее строгая оценка при размерности менее $2m-1$? (Она должна быть меньше чем $2/\sqrt{2m-1}$.)
Например в одномерном случае экстремальный пример - одна точка на одном конце, все другие точки на другом конце. Тут наиболее строгая оценка: $2/m$.


Вернуться в «Математика»

Кто сейчас на форуме

Количество пользователей, которые сейчас просматривают этот форум: нет зарегистрированных пользователей и 6 гостей