Рассмотрите задачу нахождения наибольшего произведения чисел при заданной сумме (задача оптимизации с ограничением); сравните методы Лагранжа, неявных производных и неравенств вида AM-GM
Задача (в стандартной форме): максимизировать произведение P=x1x2⋯xnP=x_1x_2\cdots x_nP=x1x2⋯xn при фиксированной сумме ∑i=1nxi=S\sum_{i=1}^n x_i=S∑i=1nxi=S. Предполагаем xi>0x_i>0xi>0 (иначе логарифм и AM–GM требуют оговорок). Решение и сравнение методов. 1) Метод Лагранжа (наиболее формальный) - Берём логарифм: lnP=∑i=1nlnxi\ln P=\sum_{i=1}^n\ln x_ilnP=∑i=1nlnxi. Максимизация PPP эквивалентна максимизации ∑lnxi\sum\ln x_i∑lnxi при ∑xi=S\sum x_i=S∑xi=S. - Лагранжиан: L=∑i=1nlnxi−λ(∑i=1nxi−S)\mathcal L=\sum_{i=1}^n\ln x_i-\lambda\Big(\sum_{i=1}^n x_i-S\Big)L=∑i=1nlnxi−λ(∑i=1nxi−S). - Условия стационарности: ∂L∂xi=1xi−λ=0⇒xi=1λ\frac{\partial\mathcal L}{\partial x_i}=\frac{1}{x_i}-\lambda=0\Rightarrow x_i=\frac{1}{\lambda}∂xi∂L=xi1−λ=0⇒xi=λ1 для всех iii. - Значит x1=⋯=xn=Snx_1=\dots=x_n=\frac{S}{n}x1=⋯=xn=nS. Максимум: Pmax=(Sn)nP_{\max}=\Big(\frac{S}{n}\Big)^nPmax=(nS)n. - Достаточность: ∑lnxi\sum\ln x_i∑lnxi — вогнутая функция, поэтому стационарная точка на допустимом множестве даёт глобальный максимум; границы (xi→0x_i\to0xi→0) дают P→0P\to0P→0, хуже. 2) Метод неявных производных (исключение переменной) - Исключаем xn=S−∑i=1n−1xix_n=S-\sum_{i=1}^{n-1}x_ixn=S−∑i=1n−1xi. Рассматриваем f(x1,…,xn−1)=∑i=1n−1lnxi+ln (S−∑i=1n−1xi)f(x_1,\dots,x_{n-1})=\sum_{i=1}^{n-1}\ln x_i+\ln\!\Big(S-\sum_{i=1}^{n-1}x_i\Big)f(x1,…,xn−1)=∑i=1n−1lnxi+ln(S−∑i=1n−1xi). - Ставим частные производные в ноль: ∂f∂xi=1xi−1S−∑j=1n−1xj=0\frac{\partial f}{\partial x_i}=\frac{1}{x_i}-\frac{1}{S-\sum_{j=1}^{n-1}x_j}=0∂xi∂f=xi1−S−∑j=1n−1xj1=0. - Отсюда все xix_ixi равны между собой и равны S/nS/nS/n. Результат тот же: Pmax=(S/n)nP_{\max}=(S/n)^nPmax=(S/n)n. - Комментарий: этот метод эквивалентен Лагранжу для одного линейного ограничения, иногда проще технически. 3) Метод неравенства AM–GM (наиболее простой и элементарный) - AM–GM: x1+⋯+xnn≥(x1⋯xn)1/n\frac{x_1+\dots+x_n}{n}\ge (x_1\cdots x_n)^{1/n}nx1+⋯+xn≥(x1⋯xn)1/n. - Подставляя сумму SSS: Sn≥(P)1/n⇒P≤(Sn)n\frac{S}{n}\ge (P)^{1/n}\Rightarrow P\le\Big(\frac{S}{n}\Big)^nnS≥(P)1/n⇒P≤(nS)n. - Равенство достигается при x1=⋯=xn=S/nx_1=\dots=x_n=S/nx1=⋯=xn=S/n. Получаем тот же максимум. - Преимущество: прямой, не требует дифференцирования; даёт глобальный результат сразу. Ограничение: применим при xi≥0x_i\ge0xi≥0 и когда цель имеет мультипликативную структуру. Короткое сравнение методов - Лагранж: системный, универсальный для дифференцируемых целей и нескольких ограничений; требует проверки сущности стационарной точки (но в нашем случае вогнутость решает). - Исключение переменной/неявные производные: проще при одном ограничении, даёт те же уравнения; по сути частный случай Лагранжа. - AM–GM: наиболее короткий и интуитивный для этой задачи; даёт явную глобальную оценку и условие равенства без анализа границ и вторых производных; но применим лишь для неотрицательных переменных и специальной формы задачи. Замечание по условиям: если допускаются отрицательные xix_ixi, поведение продукта и оптимум меняются (может не существовать конечного максимума), поэтому обычно ставят xi≥0x_i\ge0xi≥0 или xi>0x_i>0xi>0.
Решение и сравнение методов.
1) Метод Лагранжа (наиболее формальный)
- Берём логарифм: lnP=∑i=1nlnxi\ln P=\sum_{i=1}^n\ln x_ilnP=∑i=1n lnxi . Максимизация PPP эквивалентна максимизации ∑lnxi\sum\ln x_i∑lnxi при ∑xi=S\sum x_i=S∑xi =S.
- Лагранжиан: L=∑i=1nlnxi−λ(∑i=1nxi−S)\mathcal L=\sum_{i=1}^n\ln x_i-\lambda\Big(\sum_{i=1}^n x_i-S\Big)L=∑i=1n lnxi −λ(∑i=1n xi −S).
- Условия стационарности: ∂L∂xi=1xi−λ=0⇒xi=1λ\frac{\partial\mathcal L}{\partial x_i}=\frac{1}{x_i}-\lambda=0\Rightarrow x_i=\frac{1}{\lambda}∂xi ∂L =xi 1 −λ=0⇒xi =λ1 для всех iii.
- Значит x1=⋯=xn=Snx_1=\dots=x_n=\frac{S}{n}x1 =⋯=xn =nS . Максимум: Pmax=(Sn)nP_{\max}=\Big(\frac{S}{n}\Big)^nPmax =(nS )n.
- Достаточность: ∑lnxi\sum\ln x_i∑lnxi — вогнутая функция, поэтому стационарная точка на допустимом множестве даёт глобальный максимум; границы (xi→0x_i\to0xi →0) дают P→0P\to0P→0, хуже.
2) Метод неявных производных (исключение переменной)
- Исключаем xn=S−∑i=1n−1xix_n=S-\sum_{i=1}^{n-1}x_ixn =S−∑i=1n−1 xi . Рассматриваем f(x1,…,xn−1)=∑i=1n−1lnxi+ln (S−∑i=1n−1xi)f(x_1,\dots,x_{n-1})=\sum_{i=1}^{n-1}\ln x_i+\ln\!\Big(S-\sum_{i=1}^{n-1}x_i\Big)f(x1 ,…,xn−1 )=∑i=1n−1 lnxi +ln(S−∑i=1n−1 xi ).
- Ставим частные производные в ноль:
∂f∂xi=1xi−1S−∑j=1n−1xj=0\frac{\partial f}{\partial x_i}=\frac{1}{x_i}-\frac{1}{S-\sum_{j=1}^{n-1}x_j}=0∂xi ∂f =xi 1 −S−∑j=1n−1 xj 1 =0.
- Отсюда все xix_ixi равны между собой и равны S/nS/nS/n. Результат тот же: Pmax=(S/n)nP_{\max}=(S/n)^nPmax =(S/n)n.
- Комментарий: этот метод эквивалентен Лагранжу для одного линейного ограничения, иногда проще технически.
3) Метод неравенства AM–GM (наиболее простой и элементарный)
- AM–GM: x1+⋯+xnn≥(x1⋯xn)1/n\frac{x_1+\dots+x_n}{n}\ge (x_1\cdots x_n)^{1/n}nx1 +⋯+xn ≥(x1 ⋯xn )1/n.
- Подставляя сумму SSS: Sn≥(P)1/n⇒P≤(Sn)n\frac{S}{n}\ge (P)^{1/n}\Rightarrow P\le\Big(\frac{S}{n}\Big)^nnS ≥(P)1/n⇒P≤(nS )n.
- Равенство достигается при x1=⋯=xn=S/nx_1=\dots=x_n=S/nx1 =⋯=xn =S/n. Получаем тот же максимум.
- Преимущество: прямой, не требует дифференцирования; даёт глобальный результат сразу. Ограничение: применим при xi≥0x_i\ge0xi ≥0 и когда цель имеет мультипликативную структуру.
Короткое сравнение методов
- Лагранж: системный, универсальный для дифференцируемых целей и нескольких ограничений; требует проверки сущности стационарной точки (но в нашем случае вогнутость решает).
- Исключение переменной/неявные производные: проще при одном ограничении, даёт те же уравнения; по сути частный случай Лагранжа.
- AM–GM: наиболее короткий и интуитивный для этой задачи; даёт явную глобальную оценку и условие равенства без анализа границ и вторых производных; но применим лишь для неотрицательных переменных и специальной формы задачи.
Замечание по условиям: если допускаются отрицательные xix_ixi , поведение продукта и оптимум меняются (может не существовать конечного максимума), поэтому обычно ставят xi≥0x_i\ge0xi ≥0 или xi>0x_i>0xi >0.