Поставьте и решите задачу: для данного выпуклого многоугольника найти все точки P такие, что сумма расстояний от P до всех вершин минимальна; сравните решения для различных чисел вершин (треугольник, квадрат, правильный n-угольник) и обсудите методы доказательства оптимальности (геометрический, аналитический, вариационный)
Постановка. Пусть даны вершины выпуклого многоугольника A1,…,An⊂R2A_1,\dots,A_n\subset\mathbb R^2A1,…,An⊂R2. Найти все точки P∈R2P\in\mathbb R^2P∈R2, минимизирующие сумму расстояний до вершин f(P)=∑i=1n∣P−Ai∣.
f(P)=\sum_{i=1}^n |P-A_i|. f(P)=i=1∑n∣P−Ai∣. Общие свойства и условие оптимальности. - Функция fff является суммой выпуклых функций, поэтому выпукла; следовательно множество минимизаторов выпукло и (при неналичии коллинеарности всех точек) минимум единственен. - Необходимое и достаточное условие для внутренней точки PPP (звено типа «градиент = 0»): ∑i=1nP−Ai∣P−Ai∣=0.
\sum_{i=1}^n \frac{P-A_i}{|P-A_i|}=0. i=1∑n∣P−Ai∣P−Ai=0.
Это означает, что в точке минимума векторы единичных направлений к вершинам векторно компенсируются. - Если оптимум достигается в одной из вершин, скажем P=AkP=A_kP=Ak, то субградиентный критерий даёт условие ∣∑i≠kAk−Ai∣Ak−Ai∣∣≤1,
\left|\sum_{i\ne k}\frac{A_k-A_i}{|A_k-A_i|}\right|\le 1, i=k∑∣Ak−Ai∣Ak−Ai≤1,
что эквивалентно тому, что нулевая векторная сумма может быть получена добавлением некоторого вектора длины не более 1 для слагаемого с модулем неопределённым в точке совпадения. Частные случаи и сравнение. 1) Два пункта (n=2n=2n=2). Минимум на отрезке A1A2A_1A_2A1A2 достигается в середине. Формула: P=A1+A22P=\frac{A_1+A_2}{2}P=2A1+A2. 2) Треугольник (n=3n=3n=3) — задача Ферма–Торричелли. - Если в треугольнике есть угол ≥120∘\ge 120^\circ≥120∘ (то есть ≥2π/3\ge 2\pi/3≥2π/3), то точка минимума — соответствующая вершина. - Если все углы <120∘<120^\circ<120∘, то существует единственная внутренняя точка PPP (Ферма), такая что три отрезка PAiPA_iPAi попарно составляют углы по 120∘120^\circ120∘; аналитически это удовлетворяет ∑i=13P−Ai∣P−Ai∣=0.
\sum_{i=1}^3 \frac{P-A_i}{|P-A_i|}=0. i=1∑3∣P−Ai∣P−Ai=0.
Геометрическая конструкция: на сторонах строят правильные треугольники и проводят соединяющие от вершины противоположные вершины — пересечение этих прямых есть точка Ферма (доказательство через поворот на 60∘60^\circ60∘). 3) Квадрат и правильный многоугольник. - Для квадрата (и для любого правильного nnn-угольника) по симметрии центр OOO правильного многоугольника даёт стационарную точку и, ввиду выпуклости fff, является единственным минимумом: P=O.
P=O. P=O.
- Для большого nnn (правильный nnn-угольник) центр минимизирует сумму из-за симметрии; для равномерно расположенных точек центровая точка удовлетворяет векторному условию (векторы взаимно компенсируются). Замечание: в общем случае выпуклого (но не регулярного) многоугольника геометрическая медиана обычно не совпадает с центроидом (ср. центр масс) и нет простых закрытых формул. Методы доказательства оптимальности. 1) Геометрический (конструктивный): для треугольника используют поворот на 60∘60^\circ60∘ или построение наружных равносторонних треугольников — показывает, что точка с углами 120∘120^\circ120∘ минимизирует сумму. Для симметричных конфигураций (правильные многоугольники) симметрия даёт очевидную минимальную точку. 2) Аналитический: рассматривают функцию f(P)f(P)f(P), её (суб)градиент. Для внутренней оптимальной точки необходимо и достаточно условие ∑i=1nP−Ai∣P−Ai∣=0.
\sum_{i=1}^n \frac{P-A_i}{|P-A_i|}=0. i=1∑n∣P−Ai∣P−Ai=0.
Второй шаг — использовать выпуклость fff (глобальный минимум при исчезновении градиента). Можно также оценивать изменение при малом сдвиге hhh: f(P+h)−f(P)=∑i=1n(∣P+h−Ai∣−∣P−Ai∣)=∑i=1nP−Ai∣P−Ai∣⋅h+o(∣h∣).
f(P+h)-f(P)=\sum_{i=1}^n \left(|P+h-A_i|-|P-A_i|\right)=\sum_{i=1}^n \frac{P-A_i}{|P-A_i|}\cdot h+o(|h|). f(P+h)−f(P)=i=1∑n(∣P+h−Ai∣−∣P−Ai∣)=i=1∑n∣P−Ai∣P−Ai⋅h+o(∣h∣). 3) Вариационный и численный: применяют субградиентную теорию и численные алгоритмы (Weiszfeld): P(t+1)=∑i=1nAi∣P(t)−Ai∣∑i=1n1∣P(t)−Ai∣.
P^{(t+1)}=\frac{\sum_{i=1}^n \frac{A_i}{|P^{(t)}-A_i|}}{\sum_{i=1}^n \frac{1}{|P^{(t)}-A_i|}}. P(t+1)=∑i=1n∣P(t)−Ai∣1∑i=1n∣P(t)−Ai∣Ai.
Weiszfeld сходится к геометрической медиане при корректной обработке случая, когда итерация попадает точно в одну из точек данных. Краткие выводы. - Минимizer суммы расстояний — геометрическая медиана; она единственна для ненаправленных конфигураций точек и удовлетворяет суммарному векторному условию. - Для треугольника: либо внутренняя точка с углами 120∘120^\circ120∘, либо вершина при угле ≥120∘\ge 120^\circ≥120∘. - Для квадрата и любых правильных nnn-угольников — центр по симметрии. - Доказательства используют: чисто геометрические конструкции (треугольник), аналитические методы (субградиент/выпуклость) и вариационные/численные методы (Weiszfeld) для общего случая.
f(P)=∑i=1n∣P−Ai∣. f(P)=\sum_{i=1}^n |P-A_i|.
f(P)=i=1∑n ∣P−Ai ∣.
Общие свойства и условие оптимальности.
- Функция fff является суммой выпуклых функций, поэтому выпукла; следовательно множество минимизаторов выпукло и (при неналичии коллинеарности всех точек) минимум единственен.
- Необходимое и достаточное условие для внутренней точки PPP (звено типа «градиент = 0»):
∑i=1nP−Ai∣P−Ai∣=0. \sum_{i=1}^n \frac{P-A_i}{|P-A_i|}=0.
i=1∑n ∣P−Ai ∣P−Ai =0. Это означает, что в точке минимума векторы единичных направлений к вершинам векторно компенсируются.
- Если оптимум достигается в одной из вершин, скажем P=AkP=A_kP=Ak , то субградиентный критерий даёт условие
∣∑i≠kAk−Ai∣Ak−Ai∣∣≤1, \left|\sum_{i\ne k}\frac{A_k-A_i}{|A_k-A_i|}\right|\le 1,
i=k∑ ∣Ak −Ai ∣Ak −Ai ≤1, что эквивалентно тому, что нулевая векторная сумма может быть получена добавлением некоторого вектора длины не более 1 для слагаемого с модулем неопределённым в точке совпадения.
Частные случаи и сравнение.
1) Два пункта (n=2n=2n=2). Минимум на отрезке A1A2A_1A_2A1 A2 достигается в середине. Формула: P=A1+A22P=\frac{A_1+A_2}{2}P=2A1 +A2 .
2) Треугольник (n=3n=3n=3) — задача Ферма–Торричелли.
- Если в треугольнике есть угол ≥120∘\ge 120^\circ≥120∘ (то есть ≥2π/3\ge 2\pi/3≥2π/3), то точка минимума — соответствующая вершина.
- Если все углы <120∘<120^\circ<120∘, то существует единственная внутренняя точка PPP (Ферма), такая что три отрезка PAiPA_iPAi попарно составляют углы по 120∘120^\circ120∘; аналитически это удовлетворяет
∑i=13P−Ai∣P−Ai∣=0. \sum_{i=1}^3 \frac{P-A_i}{|P-A_i|}=0.
i=1∑3 ∣P−Ai ∣P−Ai =0. Геометрическая конструкция: на сторонах строят правильные треугольники и проводят соединяющие от вершины противоположные вершины — пересечение этих прямых есть точка Ферма (доказательство через поворот на 60∘60^\circ60∘).
3) Квадрат и правильный многоугольник.
- Для квадрата (и для любого правильного nnn-угольника) по симметрии центр OOO правильного многоугольника даёт стационарную точку и, ввиду выпуклости fff, является единственным минимумом:
P=O. P=O.
P=O. - Для большого nnn (правильный nnn-угольник) центр минимизирует сумму из-за симметрии; для равномерно расположенных точек центровая точка удовлетворяет векторному условию (векторы взаимно компенсируются).
Замечание: в общем случае выпуклого (но не регулярного) многоугольника геометрическая медиана обычно не совпадает с центроидом (ср. центр масс) и нет простых закрытых формул.
Методы доказательства оптимальности.
1) Геометрический (конструктивный): для треугольника используют поворот на 60∘60^\circ60∘ или построение наружных равносторонних треугольников — показывает, что точка с углами 120∘120^\circ120∘ минимизирует сумму. Для симметричных конфигураций (правильные многоугольники) симметрия даёт очевидную минимальную точку.
2) Аналитический: рассматривают функцию f(P)f(P)f(P), её (суб)градиент. Для внутренней оптимальной точки необходимо и достаточно условие
∑i=1nP−Ai∣P−Ai∣=0. \sum_{i=1}^n \frac{P-A_i}{|P-A_i|}=0.
i=1∑n ∣P−Ai ∣P−Ai =0. Второй шаг — использовать выпуклость fff (глобальный минимум при исчезновении градиента). Можно также оценивать изменение при малом сдвиге hhh:
f(P+h)−f(P)=∑i=1n(∣P+h−Ai∣−∣P−Ai∣)=∑i=1nP−Ai∣P−Ai∣⋅h+o(∣h∣). f(P+h)-f(P)=\sum_{i=1}^n \left(|P+h-A_i|-|P-A_i|\right)=\sum_{i=1}^n \frac{P-A_i}{|P-A_i|}\cdot h+o(|h|).
f(P+h)−f(P)=i=1∑n (∣P+h−Ai ∣−∣P−Ai ∣)=i=1∑n ∣P−Ai ∣P−Ai ⋅h+o(∣h∣).
3) Вариационный и численный: применяют субградиентную теорию и численные алгоритмы (Weiszfeld):
P(t+1)=∑i=1nAi∣P(t)−Ai∣∑i=1n1∣P(t)−Ai∣. P^{(t+1)}=\frac{\sum_{i=1}^n \frac{A_i}{|P^{(t)}-A_i|}}{\sum_{i=1}^n \frac{1}{|P^{(t)}-A_i|}}.
P(t+1)=∑i=1n ∣P(t)−Ai ∣1 ∑i=1n ∣P(t)−Ai ∣Ai . Weiszfeld сходится к геометрической медиане при корректной обработке случая, когда итерация попадает точно в одну из точек данных.
Краткие выводы.
- Минимizer суммы расстояний — геометрическая медиана; она единственна для ненаправленных конфигураций точек и удовлетворяет суммарному векторному условию.
- Для треугольника: либо внутренняя точка с углами 120∘120^\circ120∘, либо вершина при угле ≥120∘\ge 120^\circ≥120∘.
- Для квадрата и любых правильных nnn-угольников — центр по симметрии.
- Доказательства используют: чисто геометрические конструкции (треугольник), аналитические методы (субградиент/выпуклость) и вариационные/численные методы (Weiszfeld) для общего случая.