Метод множителей лагранжа примеры. Метод множителей Лагранжа
Решение . Найдем экстремум функции F(X) = 9 x 1 +x 1 2 +6 x 2 +x 2 2 , используя функцию Лагранжа:
где
- целевая функция вектора .
- ограничения в неявном виде (i=1..n)
В качестве целевой функции, подлежащей оптимизации, в этой задаче выступает функция:
F(X) = 9 x 1 +x 1 2 +6 x 2 +x 2 2
Перепишем ограничение задачи в неявном виде:
Составим вспомогательную функцию Лагранжа:
= 9 x 1 +x 1 2 +6 x 2 +x 2 2 + λ(x 1 +x 2 -150)
Необходимым условием экстремума функции Лагранжа является равенство нулю ее частных производных по переменным х i и неопределенному множителю λ.
Составим систему:
∂L/∂x 1 = 2 x 1 +λ+9 = 0
∂L/∂x 2 = λ+2 x 2 +6 = 0
∂F/∂λ = x 1 +x 2 -150= 0
Систему решаем с помощью метода Гаусса или используя формулы Крамера .
Запишем систему в виде:
Для удобства вычислений поменяем строки местами:
Добавим 2-ую строку к 1-ой:
Умножим 2-ую строку на (2). Умножим 3-ую строку на (-1). Добавим 3-ую строку к 2-ой:
Умножим 2-ую строку на (-1). Добавим 2-ую строку к 1-ой:
Из 1-ой строки выражаем x 3
Из 2-ой строки выражаем x 2
Из 3-ой строки выражаем x 1
Таким образом, чтобы общие издержки производства были минимальны, необходимо производить x 1 = 74.25; x 2 = 75.75.
Задание . По плану производства продукции предприятию необходимо изготовить 50 изделий. Эти изделия могут быть изготовлены 2-мя технологическими способами. При производстве x 1 - изделий 1-ым способом затраты равны 3x 1 +x 1 2 (т. руб.), а при изготовлении x 2 - изделий 2-ым способом они составят 5x 2 +x 2 2 (т. руб.). Определить сколько изделий каждым из способов необходимо изготовить, чтобы общие затраты на производство были минимальные.
Решение: составляем целевую функцию и ограничения:
F(X) = 3x 1 +x 1 2 + 5x 2 +x 2 2 → min
x 1 +x 2 = 50
an(t)z(n)(t) + an − 1(t)z(n − 1)(t) + ... + a1(t)z"(t) + a0(t)z(t) = f(t)
состоит в замене произвольных постоянных ck в общем решении
z(t) = c1z1(t) + c2z2(t) + ...
Cnzn(t)
соответствующего однородного уравнения
an(t)z(n)(t) + an − 1(t)z(n − 1)(t) + ... + a1(t)z"(t) + a0(t)z(t) = 0
на вспомогательные функции ck(t), производные которых удовлетворяют линейной алгебраической системе
Определителем системы (1) служит вронскиан функций z1,z2,...,zn, что обеспечивает её однозначную разрешимость относительно .
Если - первообразные для , взятые при фиксированных значениях постоянных интегрирования, то функция
является решением исходного линейного неоднородного дифференциального уравнения. Интегрирование неоднородного уравнения при наличии общего решения соответствующего однородного уравнения сводится, таким образом, к квадратурам.
Метод Лагранжа (метод вариации произвольных постоянных)
Метод для получения общего решения неоднородного уравнения, зная общее решение однородного уравнения без нахождения частного решения.
Для линейного однородного дифференциального уравнения n-го порядка
y(n) + a1(x) y(n-1) + ... + an-1 (x) y" + an(x) y = 0,
где y = y(x) - неизвестная функция, a1(x), a2(x), ..., an-1(x), an(x) - известные, непрерывные, справедливо: 1) существуют n линейно независимых решений уравнения y1(x), y2(x), ..., yn(x); 2) при любых значениях констант c1, c2, ..., cn функция y(x)= c1 y1(x) + c2 y2(x) + ... + cn yn(x) является решением уравнения; 3) для любых начальных значений x0, y0, y0,1, ..., y0,n-1 существуют такие значения c*1, c*n, ..., c*n, что решение y*(x)=c*1 y1(x) + c*2 y2(x) + ... + c*n yn (x) удовлетворяет при x = x0 начальным условиям y*(x0)=y0, (y*)"(x0)=y0,1 , ...,(y*)(n-1)(x0)=y0,n-1.
Выражение y(x)= c1 y1(x) + c2 y2(x) + ... + cn yn(x) называется общим решением линейного однородного дифференциального уравнения n-го порядка.
Совокупность n линейно независимых решений линейного однородного дифференциального уравнения n-го порядка y1(x), y2(x), ..., yn(x) называется фундаментальной системой решений уравнения.
Для линейного однородного дифференциального уравнения с постоянными коэффициентами существует простой алгоритм построения фундаментальной системы решений. Будем искать решение уравнения в виде y(x) = exp(lx): exp(lx)(n) + a1exp(lx)(n-1) + ... + an-1exp(lx)" + anexp(lx)= = (ln + a1ln-1 + ... + an-1l + an)exp(lx) = 0, т.е. число l является корнем характеристического уравнения ln + a1ln-1 + ... + an-1l + an = 0. Левая часть характеристического уравнения называется характеристическим многочленом линейного дифференциального уравнения: P(l) = ln + a1ln-1 + ... + an-1l + an. Таким образом, задача о решении линейного однородного уравнения n -го порядка с постоянными коэффициентами сводится к решению алгебраического уравнения.
Если характеристическое уравнение имеет n различных действительных корней l1№ l2 № ... № ln, то фундаментальная система решений состоит из функций y1(x) = exp(l1x), y2(x) = exp(l2x), ..., yn(x) = exp(lnx), и общее решение однородного уравнения имеет вид: y(x)= c1 exp(l1x) + c2 exp(l2x) + ... + cn exp(lnx).
ундаментальная система решений и общее решение для случая простых действительных корней.
Если какой-либо из действительных корней характеристического уравнения повторяется r раз (r-кратный корень), то в фундаментальной системе решений ему отвечают r функций; если lk=lk+1 = ... = lk+r-1, то в фундаментальную систему решений уравнения входят r функций: yk(x) = exp(lkx), yk+1(x) = xexp(lkx), yk+2(x) = x2exp(lkx), ..., yk+r-1(x) =xr-1 exp(lnx).
ПРИМЕР 2. Фундаментальная система решений и общее решение для случая кратных действительных корней.
Если характеристическое уравнение имеет комплексные корни, то каждой паре простых (имеющих кратность 1) комплексных корней lk,k+1=ak ± ibk в фундаментальной системе решений отвечает пара функций yk(x) = exp(akx)cos(bkx), yk+1(x) = exp(akx)sin(bkx).
ПРИМЕР 4. Фундаментальная система решений и общее решение для случая простых комплексных корней. Мнимые корни.
Если же комплексная пара корней имеет кратность r, то такой паре lk=lk+1 = ... = l2k+2r-1=ak ± ibk, в фундаментальной системе решений отвечают функции exp(akx)cos(bkx), exp(akx)sin(bkx), xexp(akx)cos(bkx), xexp(akx)sin(bkx), x2exp(akx)cos(bkx), x2exp(akx)sin(bkx), ................ xr-1exp(akx)cos(bkx), xr-1exp(akx)sin(bkx).
ПРИМЕР 5. Фундаментальная система решений и общее решение для случая кратных комплексных корней.
Таким образом, для отыскания общего решения линейного однородного дифференциального уравнения с постоянными коэффициентами следует: записать характеристическое уравнение; найти все корни характеристического уравнения l1, l2, ... , ln; записать фундаментальную систему решений y1(x), y2(x), ..., yn(x); записать выражение для общего решения y(x)= c1 y1(x) + c2 y2(x) + ... + cn yn(x). Для решения задачи Коши нужно подставить выражение для общего решения в начальные условия и определить значения постоянных c1,..., cn, которые являются решениями системы линейных алгебраических уравнений c1 y1(x0) + c2 y2(x0) + ... + cn yn(x0) = y0, c1 y"1(x0) + c2 y"2(x0) + ... + cn y"n(x0) =y0,1, ......... , c1 y1(n-1)(x0) + c2 y2(n-1)(x0) + ... + cn yn(n-1)(x0) = y0,n-1
Для линейного неоднородного дифференциального уравнения n-го порядка
y(n) + a1(x) y(n-1) + ... + an-1 (x) y" + an(x) y = f(x),
где y = y(x) - неизвестная функция, a1(x), a2(x), ..., an-1(x), an(x), f(x) - известные, непрерывные, справедливо: 1) если y1(x) и y2(x) - два решения неоднородного уравнения, то функция y(x) = y1(x) - y2(x) - решение соответствующего однородного уравнения; 2) если y1(x) решение неоднородного уравнения, а y2(x) - решение соответствующего однородного уравнения, то функция y(x) = y1(x) + y2(x) - решение неоднородного уравнения; 3) если y1(x), y2(x), ..., yn(x) - n линейно независимых решений однородного уравнения, а yч(x) - произвольное решение неоднородного уравнения, то для любых начальных значений x0, y0, y0,1, ..., y0,n-1 существуют такие значения c*1, c*n, ..., c*n, что решение y*(x)=c*1 y1(x) + c*2 y2(x) + ... + c*n yn (x) + yч(x) удовлетворяет при x = x0 начальным условиям y*(x0)=y0, (y*)"(x0)=y0,1 , ...,(y*)(n-1)(x0)=y0,n-1.
Выражение y(x)= c1 y1(x) + c2 y2(x) + ... + cn yn(x) + yч(x) называется общим решением линейного неоднородного дифференциального уравнения n-го порядка.
Для отыскания частных решений неоднородных дифференциальных уравнений с постоянными коэффициентами с правыми частями вида: Pk(x)exp(ax)cos(bx) + Qm(x)exp(ax)sin(bx), где Pk(x), Qm(x) - многочлены степени k и m соответственно, существует простой алгоритм построения частного решения, называемый методом подбора.
Метод подбора, или метод неопределенных коэффициентов, состоит в следующем. Искомое решение уравнения записывается в виде: (Pr(x)exp(ax)cos(bx) + Qr(x)exp(ax)sin(bx))xs, где Pr(x), Qr(x) - многочлены степени r = max(k, m) с неизвестными коэффициентами pr , pr-1, ..., p1, p0, qr, qr-1, ..., q1, q0. Сомножитель xs называют резонансным сомножителем. Резонанс имеет место в случаях, когда среди корней характеристического уравнения есть корень l =a ± ib кратности s. Т.е. если среди корней характеристического уравнения соответствующего однородного уравнения есть такой, что его действительная часть совпадает с коэффициентом в показателе степени экспоненты, а мнимая - с коэффициентом в аргументе тригонометрической функции в правой части уравнения, и кратность этого корня s, то в искомом частном решении присутствует резонансный сомножитель xs. Если же такого совпадения нет (s=0), то резонансный сомножитель отсутствует.
Подставив выражение для частного решения в левую часть уравнения, получим обобщенный многочлен того же вида, что и многочлен в правой части уравнения, коэффициенты которого неизвестны.
Два обобщенных многочлена равны тогда и только тогда, когда равны коэффициенты при сомножителях вида xtexp(ax)sin(bx), xtexp(ax)cos(bx) с одинаковыми степенями t. Приравняв коэффициенты при таких сомножителях, получим систему 2(r+1) линейных алгебраических уравнений относительно 2(r+1) неизвестных. Можно показать, что такая система совместна и имеет единственное решение.
|
Метод Лагранжа ─ это метод решения задачи условной оптимизации, при котором ограничения, записываемые как неявные функции, объединяются с целевой функцией в форме нового уравнения, называемого лагранжианом .
Рассмотрим частный случай общей задачи нелинейного программирования:
Дана система нелинейных уравнений (1):
(1) gi(x1,x2,…,xn)=bi (i=1..m),
Найти наименьшее (или наибольшее) значение функции (2)
(2) f (х1,х2,…,хn),
если отсутствуют условия неотрицательности переменных и f(х1,х2,…,хn) и gi(x1,x2,…,xn) ─ функции, непрерывные вместе со своими частными производными.
Чтобы найти решение этой задачи, можно применить следующий метод: 1. Вводят набор переменных λ1, λ2,…, λm, называемых множителями Лагранжа, составляют функцию Лагранжа (3)
(3) F(х1,х2,…,хn , λ1,λ2,…,λm) = f(х1,х2,…,хn)+ λi .
2. Находят частные производные от функции Лагранжа по переменным xi и λi и приравнивают их нулю.
3. Решая систему уравнений, находят точки, в которых целевая функция задачи может иметь экстремум.
4.Среди точек, подозрительных не экстремум, находят такие, в которыхдостигается экстремум, и вычисляют значения функции в этих точках.
4. Сравнить полученные значения функции f и выбрать наилучшее.
По плану производства продукции предприятию необходимо изготовить 180 изделий. Эти изделия могут быть изготовлены двумя технологическими способами. При производстве х1 изделия I способом затраты равны 4*х1+х1^2 руб., а при изготовлении х2 изделий II способом они составляют 8*х2+х2^2 руб. Определить, сколько изделий каждым из способов следует изготовить, так чтобы общие затраты на производство продукции были минимальными.
Решение: Математическая постановка задачи состоит в определении наименьшего значения функции двух переменных:
f = 4*x1+x1^2 +8*x2 +x2^2, при условии x1 +x2 = 180.
Составим функцию Лагранжа:
F(x1,x2,λ) = 4*x1+x1^2+8*x2+x2^2+λ*(180-x1-x2).
Вычислим ее частные производные по х1,х2, λ и приравняем их к 0:
Перенесем в правые части первых двух уравнений λ и приравняем их левые части, получим 4 + 2*x1 = 8 + 2*x2, или x1 − x2 = 2.
Решая последнее уравнение совместно с уравнением x1 + x2 = 180, находим x1 = 91, x2 = 89, то есть получили решение, удовлетворяющее условиям:
Найдем значение целевой функции f при этих значениях переменных:
F(x1, x2) = 17278
Эта точка является подозрительной на экстремум. Используя вторые частные производные, можно показать, что в точке (91,89) функция f имеет минимум.
Для начала рассмотрим случай функции двух переменных. Условным экстремумом функции $z=f(x,y)$ в точке $M_0(x_0;y_0)$ называется экстремум этой функции, достигнутый при условии, что переменные $x$ и $y$ в окрестности данной точки удовлетворяют уравнению связи $\varphi (x,y)=0$.
Название «условный» экстремум связано с тем, что на переменные наложено дополнительное условие $\varphi(x,y)=0$. Если из уравнения связи можно выразить одну переменную через другую, то задача определения условного экстремума сводится к задаче на обычный экстремум функции одной переменной. Например, если из уравнения связи следует $y=\psi(x)$, то подставив $y=\psi(x)$ в $z=f(x,y)$, получим функцию одной переменной $z=f\left(x,\psi(x)\right)$. В общем случае, однако, такой метод малопригоден, поэтому требуется введение нового алгоритма.
Метод множителей Лагранжа для функций двух переменных.
Метод множителей Лагранжа состоит в том, что для отыскания условного экстремума составляют функцию Лагранжа: $F(x,y)=f(x,y)+\lambda\varphi(x,y)$ (параметр $\lambda$ называют множителем Лагранжа). Необходимые условия экстремума задаются системой уравнений, из которой определяются стационарные точки:
$$ \left \{ \begin{aligned} & \frac{\partial F}{\partial x}=0;\\ & \frac{\partial F}{\partial y}=0;\\ & \varphi (x,y)=0. \end{aligned} \right. $$
Достаточным условием, из которого можно выяснить характер экстремума, служит знак $d^2 F=F_{xx}^{""}dx^2+2F_{xy}^{""}dxdy+F_{yy}^{""}dy^2$. Если в стационарной точке $d^2F > 0$, то функция $z=f(x,y)$ имеет в данной точке условный минимум, если же $d^2F < 0$, то условный максимум.
Есть и другой способ для определения характера экстремума. Из уравнения связи получаем: $\varphi_{x}^{"}dx+\varphi_{y}^{"}dy=0$, $dy=-\frac{\varphi_{x}^{"}}{\varphi_{y}^{"}}dx$, поэтому в любой стационарной точке имеем:
$$d^2 F=F_{xx}^{""}dx^2+2F_{xy}^{""}dxdy+F_{yy}^{""}dy^2=F_{xx}^{""}dx^2+2F_{xy}^{""}dx\left(-\frac{\varphi_{x}^{"}}{\varphi_{y}^{"}}dx\right)+F_{yy}^{""}\left(-\frac{\varphi_{x}^{"}}{\varphi_{y}^{"}}dx\right)^2=\\ =-\frac{dx^2}{\left(\varphi_{y}^{"} \right)^2}\cdot\left(-(\varphi_{y}^{"})^2 F_{xx}^{""}+2\varphi_{x}^{"}\varphi_{y}^{"}F_{xy}^{""}-(\varphi_{x}^{"})^2 F_{yy}^{""} \right)$$
Второй сомножитель (расположенный в скобке) можно представить в такой форме:
Красным цветом выделены элементы определителя $\left| \begin{array} {cc} F_{xx}^{""} & F_{xy}^{""} \\ F_{xy}^{""} & F_{yy}^{""} \end{array} \right|$, который является гессианом функции Лагранжа. Если $H > 0$, то $d^2F < 0$, что указывает на условный максимум. Аналогично, при $H < 0$ имеем $d^2F > 0$, т.е. имеем условный минимум функции $z=f(x,y)$.
Примечание относительно формы записи определителя $H$. показать\скрыть
$$ H=-\left|\begin{array} {ccc} 0 & \varphi_{x}^{"} & \varphi_{y}^{"}\\ \varphi_{x}^{"} & F_{xx}^{""} & F_{xy}^{""} \\ \varphi_{y}^{"} & F_{xy}^{""} & F_{yy}^{""} \end{array} \right| $$
В этой ситуации сформулированное выше правило изменится следующим образом: если $H > 0$, то функция имеет условный минимум, а при $H < 0$ получим условный максимум функции $z=f(x,y)$. При решении задач следует учитывать такие нюансы.
Алгоритм исследования функции двух переменных на условный экстремум
- Составить функцию Лагранжа $F(x,y)=f(x,y)+\lambda\varphi(x,y)$
- Решить систему $ \left \{ \begin{aligned} & \frac{\partial F}{\partial x}=0;\\ & \frac{\partial F}{\partial y}=0;\\ & \varphi (x,y)=0. \end{aligned} \right.$
- Определить характер экстремума в каждой из найденных в предыдущем пункте стационарных точек. Для этого применить любой из указанных способов:
- Составить определитель $H$ и выяснить его знак
- С учетом уравнения связи вычислить знак $d^2F$
Метод множителей Лагранжа для функций n переменных
Допустим, мы имеем функцию $n$ переменных $z=f(x_1,x_2,\ldots,x_n)$ и $m$ уравнений связи ($n > m$):
$$\varphi_1(x_1,x_2,\ldots,x_n)=0; \; \varphi_2(x_1,x_2,\ldots,x_n)=0,\ldots,\varphi_m(x_1,x_2,\ldots,x_n)=0.$$
Обозначив множители Лагранжа как $\lambda_1,\lambda_2,\ldots,\lambda_m$, составим функцию Лагранжа:
$$F(x_1,x_2,\ldots,x_n,\lambda_1,\lambda_2,\ldots,\lambda_m)=f+\lambda_1\varphi_1+\lambda_2\varphi_2+\ldots+\lambda_m\varphi_m$$
Необходимые условия наличия условного экстремума задаются системой уравнений, из которой находятся координаты стационарных точек и значения множителей Лагранжа:
$$\left\{\begin{aligned} & \frac{\partial F}{\partial x_i}=0; (i=\overline{1,n})\\ & \varphi_j=0; (j=\overline{1,m}) \end{aligned} \right.$$
Выяснить, условный минимум или условный максимум имеет функция в найденной точке, можно, как и ранее, посредством знака $d^2F$. Если в найденной точке $d^2F > 0$, то функция имеет условный минимум, если же $d^2F < 0$, - то условный максимум. Можно пойти иным путем, рассмотрев следующую матрицу:
Определитель матрицы $\left| \begin{array} {ccccc} \frac{\partial^2F}{\partial x_{1}^{2}} & \frac{\partial^2F}{\partial x_{1}\partial x_{2}} & \frac{\partial^2F}{\partial x_{1}\partial x_{3}} &\ldots & \frac{\partial^2F}{\partial x_{1}\partial x_{n}}\\ \frac{\partial^2F}{\partial x_{2}\partial x_1} & \frac{\partial^2F}{\partial x_{2}^{2}} & \frac{\partial^2F}{\partial x_{2}\partial x_{3}} &\ldots & \frac{\partial^2F}{\partial x_{2}\partial x_{n}}\\ \frac{\partial^2F}{\partial x_{3} \partial x_{1}} & \frac{\partial^2F}{\partial x_{3}\partial x_{2}} & \frac{\partial^2F}{\partial x_{3}^{2}} &\ldots & \frac{\partial^2F}{\partial x_{3}\partial x_{n}}\\ \ldots & \ldots & \ldots &\ldots & \ldots\\ \frac{\partial^2F}{\partial x_{n}\partial x_{1}} & \frac{\partial^2F}{\partial x_{n}\partial x_{2}} & \frac{\partial^2F}{\partial x_{n}\partial x_{3}} &\ldots & \frac{\partial^2F}{\partial x_{n}^{2}}\\ \end{array} \right|$, выделенной в матрице $L$ красным цветом, есть гессиан функции Лагранжа. Используем следующее правило:
- Если знаки угловых миноров $H_{2m+1},\; H_{2m+2},\ldots,H_{m+n}$ матрицы $L$ совпадают с знаком $(-1)^m$, то исследуемая стационарная точка является точкой условного минимума функции $z=f(x_1,x_2,x_3,\ldots,x_n)$.
- Если знаки угловых миноров $H_{2m+1},\; H_{2m+2},\ldots,H_{m+n}$ чередуются, причём знак минора $H_{2m+1}$ совпадает с знаком числа $(-1)^{m+1}$, то исследуемая стационарная точка является точкой условного максимума функции $z=f(x_1,x_2,x_3,\ldots,x_n)$.
Пример №1
Найти условный экстремум функции $z(x,y)=x+3y$ при условии $x^2+y^2=10$.
Геометрическая интерпретация данной задачи такова: требуется найти наибольшее и наименьшее значение аппликаты плоскости $z=x+3y$ для точек ее пересечения с цилиндром $x^2+y^2=10$.
Выразить одну переменную через другую из уравнения связи и подставить ее в функцию $z(x,y)=x+3y$ несколько затруднительно, поэтому будем использовать метод Лагранжа.
Обозначив $\varphi(x,y)=x^2+y^2-10$, составим функцию Лагранжа:
$$ F(x,y)=z(x,y)+\lambda \varphi(x,y)=x+3y+\lambda(x^2+y^2-10);\\ \frac{\partial F}{\partial x}=1+2\lambda x; \frac{\partial F}{\partial y}=3+2\lambda y. $$
Запишем систему уравнений для определения стационарных точек функции Лагранжа:
$$ \left \{ \begin{aligned} & 1+2\lambda x=0;\\ & 3+2\lambda y=0;\\ & x^2+y^2-10=0. \end{aligned} \right. $$
Если предположить $\lambda=0$, то первое уравнение станет таким: $1=0$. Полученное противоречие говорит о том, что $\lambda\neq 0$. При условии $\lambda\neq 0$ из первого и второго уравнений имеем: $x=-\frac{1}{2\lambda}$, $y=-\frac{3}{2\lambda}$. Подставляя полученные значения в третье уравнение, получим:
$$ \left(-\frac{1}{2\lambda} \right)^2+\left(-\frac{3}{2\lambda} \right)^2-10=0;\\ \frac{1}{4\lambda^2}+\frac{9}{4\lambda^2}=10; \lambda^2=\frac{1}{4}; \left[ \begin{aligned} & \lambda_1=-\frac{1}{2};\\ & \lambda_2=\frac{1}{2}. \end{aligned} \right.\\ \begin{aligned} & \lambda_1=-\frac{1}{2}; \; x_1=-\frac{1}{2\lambda_1}=1; \; y_1=-\frac{3}{2\lambda_1}=3;\\ & \lambda_2=\frac{1}{2}; \; x_2=-\frac{1}{2\lambda_2}=-1; \; y_2=-\frac{3}{2\lambda_2}=-3.\end{aligned} $$
Итак, система имеет два решения: $x_1=1;\; y_1=3;\; \lambda_1=-\frac{1}{2}$ и $x_2=-1;\; y_2=-3;\; \lambda_2=\frac{1}{2}$. Выясним характер экстремума в каждой стационарной точке: $M_1(1;3)$ и $M_2(-1;-3)$. Для этого вычислим определитель $H$ в каждой из точек.
$$ \varphi_{x}^{"}=2x;\; \varphi_{y}^{"}=2y;\; F_{xx}^{""}=2\lambda;\; F_{xy}^{""}=0;\; F_{yy}^{""}=2\lambda.\\ H=\left| \begin{array} {ccc} 0 & \varphi_{x}^{"} & \varphi_{y}^{"}\\ \varphi_{x}^{"} & F_{xx}^{""} & F_{xy}^{""} \\ \varphi_{y}^{"} & F_{xy}^{""} & F_{yy}^{""} \end{array} \right|= \left| \begin{array} {ccc} 0 & 2x & 2y\\ 2x & 2\lambda & 0 \\ 2y & 0 & 2\lambda \end{array} \right|= 8\cdot\left| \begin{array} {ccc} 0 & x & y\\ x & \lambda & 0 \\ y & 0 & \lambda \end{array} \right| $$
В точке $M_1(1;3)$ получим: $H=8\cdot\left| \begin{array} {ccc} 0 & x & y\\ x & \lambda & 0 \\ y & 0 & \lambda \end{array} \right|= 8\cdot\left| \begin{array} {ccc} 0 & 1 & 3\\ 1 & -1/2 & 0 \\ 3 & 0 & -1/2 \end{array} \right|=40 > 0$, поэтому в точке $M_1(1;3)$ функция $z(x,y)=x+3y$ имеет условный максимум, $z_{\max}=z(1;3)=10$.
Аналогично, в точке $M_2(-1;-3)$ найдем: $H=8\cdot\left| \begin{array} {ccc} 0 & x & y\\ x & \lambda & 0 \\ y & 0 & \lambda \end{array} \right|= 8\cdot\left| \begin{array} {ccc} 0 & -1 & -3\\ -1 & 1/2 & 0 \\ -3 & 0 & 1/2 \end{array} \right|=-40$. Так как $H < 0$, то в точке $M_2(-1;-3)$ имеем условный минимум функции $z(x,y)=x+3y$, а именно: $z_{\min}=z(-1;-3)=-10$.
Отмечу, что вместо вычисления значения определителя $H$ в каждой точке, гораздо удобнее раскрыть его в общем виде. Дабы не загромождать текст подробностями, этот способ скрою под примечание.
Запись определителя $H$ в общем виде. показать\скрыть
$$ H=8\cdot\left|\begin{array}{ccc}0&x&y\\x&\lambda&0\\y&0&\lambda\end{array}\right| =8\cdot\left(-\lambda{y^2}-\lambda{x^2}\right) =-8\lambda\cdot\left(y^2+x^2\right). $$
В принципе, уже очевидно, какой знак имеет $H$. Так как ни одна из точек $M_1$ или $M_2$ не совпадает с началом координат, то $y^2+x^2>0$. Следовательно, знак $H$ противоположен знаку $\lambda$. Можно и довести вычисления до конца:
$$ \begin{aligned} &H(M_1)=-8\cdot\left(-\frac{1}{2}\right)\cdot\left(3^2+1^2\right)=40;\\ &H(M_2)=-8\cdot\frac{1}{2}\cdot\left((-3)^2+(-1)^2\right)=-40. \end{aligned} $$
Вопрос о характере экстремума в стационарных точках $M_1(1;3)$ и $M_2(-1;-3)$ можно решить и без использования определителя $H$. Найдем знак $d^2F$ в каждой стационарной точке:
$$ d^2 F=F_{xx}^{""}dx^2+2F_{xy}^{""}dxdy+F_{yy}^{""}dy^2=2\lambda \left(dx^2+dy^2\right) $$
Отмечу, что запись $dx^2$ означает именно $dx$, возведённый в вторую степень, т.е. $\left(dx \right)^2$. Отсюда имеем: $dx^2+dy^2>0$, посему при $\lambda_1=-\frac{1}{2}$ получим $d^2F < 0$. Следовательно, функция имеет в точке $M_1(1;3)$ условный максимум. Аналогично, в точке $M_2(-1;-3)$ получим условный минимум функции $z(x,y)=x+3y$. Отметим, что для определения знака $d^2F$ не пришлось учитывать связь между $dx$ и $dy$, ибо знак $d^2F$ очевиден без дополнительных преобразований. В следующем примере для определения знака $d^2F$ уже будет необходимо учесть связь между $dx$ и $dy$.
Ответ : в точке $(-1;-3)$ функция имеет условный минимум, $z_{\min}=-10$. В точке $(1;3)$ функция имеет условный максимум, $z_{\max}=10$
Пример №2
Найти условный экстремум функции $z(x,y)=3y^3+4x^2-xy$ при условии $x+y=0$.
Первый способ (метод множителей Лагранжа)
Обозначив $\varphi(x,y)=x+y$ составим функцию Лагранжа: $F(x,y)=z(x,y)+\lambda \varphi(x,y)=3y^3+4x^2-xy+\lambda(x+y)$.
$$ \frac{\partial F}{\partial x}=8x-y+\lambda; \; \frac{\partial F}{\partial y}=9y^2-x+\lambda.\\ \left \{ \begin{aligned} & 8x-y+\lambda=0;\\ & 9y^2-x+\lambda=0; \\ & x+y=0. \end{aligned} \right. $$
Решив систему, получим: $x_1=0$, $y_1=0$, $\lambda_1=0$ и $x_2=\frac{10}{9}$, $y_2=-\frac{10}{9}$, $\lambda_2=-10$. Имеем две стационарные точки: $M_1(0;0)$ и $M_2 \left(\frac{10}{9};-\frac{10}{9} \right)$. Выясним характер экстремума в каждой стационарной точке с использованием определителя $H$.
$$ H=\left| \begin{array} {ccc} 0 & \varphi_{x}^{"} & \varphi_{y}^{"}\\ \varphi_{x}^{"} & F_{xx}^{""} & F_{xy}^{""} \\ \varphi_{y}^{"} & F_{xy}^{""} & F_{yy}^{""} \end{array} \right|= \left| \begin{array} {ccc} 0 & 1 & 1\\ 1 & 8 & -1 \\ 1 & -1 & 18y \end{array} \right|=-10-18y $$
В точке $M_1(0;0)$ $H=-10-18\cdot 0=-10 < 0$, поэтому $M_1(0;0)$ есть точка условного минимума функции $z(x,y)=3y^3+4x^2-xy$, $z_{\min}=0$. В точке $M_2\left(\frac{10}{9};-\frac{10}{9}\right)$ $H=10 > 0$, посему в данной точке функция имеет условный максимум, $z_{\max}=\frac{500}{243}$.
Исследуем характер экстремума в каждой из точек иным методом, основываясь на знаке $d^2F$:
$$ d^2 F=F_{xx}^{""}dx^2+2F_{xy}^{""}dxdy+F_{yy}^{""}dy^2=8dx^2-2dxdy+18ydy^2 $$
Из уравнения связи $x+y=0$ имеем: $d(x+y)=0$, $dx+dy=0$, $dy=-dx$.
$$ d^2 F=8dx^2-2dxdy+18ydy^2=8dx^2-2dx(-dx)+18y(-dx)^2=(10+18y)dx^2 $$
Так как $ d^2F \Bigr|_{M_1}=10 dx^2 > 0$, то $M_1(0;0)$ является точкой условного минимума функции $z(x,y)=3y^3+4x^2-xy$. Аналогично, $d^2F \Bigr|_{M_2}=-10 dx^2 < 0$, т.е. $M_2\left(\frac{10}{9}; -\frac{10}{9} \right)$ - точка условного максимума.
Второй способ
Из уравнения связи $x+y=0$ получим: $y=-x$. Подставив $y=-x$ в функцию $z(x,y)=3y^3+4x^2-xy$, получим некоторую функцию переменной $x$. Обозначим эту функцию как $u(x)$:
$$ u(x)=z(x,-x)=3\cdot(-x)^3+4x^2-x\cdot(-x)=-3x^3+5x^2. $$
Таким образом задачу о нахождении условного экстремума функции двух переменных мы свели к задаче определения экстремума функции одной переменной.
$$ u_{x}^{"}=-9x^2+10x;\\ -9x^2+10x=0; \; x\cdot(-9x+10)=0;\\ x_1=0; \; y_1=-x_1=0;\\ x_2=\frac{10}{9}; \; y_2=-x_2=-\frac{10}{9}. $$
Получили точки $M_1(0;0)$ и $M_2\left(\frac{10}{9}; -\frac{10}{9}\right)$. Дальнейшее исследование известно из курса дифференциального исчисления функций одной переменой. Исследуя знак $u_{xx}^{""}$ в каждой стационарной точке или проверяя смену знака $u_{x}^{"}$ в найденных точках, получим те же выводы, что и при решении первым способом. Например, проверим знак $u_{xx}^{""}$:
$$u_{xx}^{""}=-18x+10;\\ u_{xx}^{""}(M_1)=10;\;u_{xx}^{""}(M_2)=-10.$$
Так как $u_{xx}^{""}(M_1)>0$, то $M_1$ - точка минимума функции $u(x)$, при этом $u_{\min}=u(0)=0$. Так как $u_{xx}^{""}(M_2)<0$, то $M_2$ - точка максимума функции $u(x)$, причём $u_{\max}=u\left(\frac{10}{9}\right)=\frac{500}{243}$.
Значения функции $u(x)$ при заданном условии связи совпадают с значениями функции $z(x,y)$, т.е. найденные экстремумы функции $u(x)$ и есть искомые условные экстремумы функции $z(x,y)$.
Ответ : в точке $(0;0)$ функция имеет условный минимум, $z_{\min}=0$. В точке $\left(\frac{10}{9}; -\frac{10}{9} \right)$ функция имеет условный максимум, $z_{\max}=\frac{500}{243}$.
Рассмотрим еще один пример, в котором характер экстремума выясним посредством определения знака $d^2F$.
Пример №3
Найти наибольшее и наименьшее значения функции $z=5xy-4$, если переменные $x$ и $y$ положительны и удовлетворяют уравнению связи $\frac{x^2}{8}+\frac{y^2}{2}-1=0$.
Составим функцию Лагранжа: $F=5xy-4+\lambda \left(\frac{x^2}{8}+\frac{y^2}{2}-1 \right)$. Найдем стационарные точки функции Лагранжа:
$$ F_{x}^{"}=5y+\frac{\lambda x}{4}; \; F_{y}^{"}=5x+\lambda y.\\ \left \{ \begin{aligned} & 5y+\frac{\lambda x}{4}=0;\\ & 5x+\lambda y=0;\\ & \frac{x^2}{8}+\frac{y^2}{2}-1=0;\\ & x > 0; \; y > 0. \end{aligned} \right. $$
Все дальнейшие преобразования осуществляются с учетом $x > 0; \; y > 0$ (это оговорено в условии задачи). Из второго уравнения выразим $\lambda=-\frac{5x}{y}$ и подставим найденное значение в первое уравнение: $5y-\frac{5x}{y}\cdot \frac{x}{4}=0$, $4y^2-x^2=0$, $x=2y$. Подставляя $x=2y$ в третье уравнение, получим: $\frac{4y^2}{8}+\frac{y^2}{2}-1=0$, $y^2=1$, $y=1$.
Так как $y=1$, то $x=2$, $\lambda=-10$. Характер экстремума в точке $(2;1)$ определим, исходя из знака $d^2F$.
$$ F_{xx}^{""}=\frac{\lambda}{4}; \; F_{xy}^{""}=5; \; F_{yy}^{""}=\lambda. $$
Так как $\frac{x^2}{8}+\frac{y^2}{2}-1=0$, то:
$$ d\left(\frac{x^2}{8}+\frac{y^2}{2}-1\right)=0; \; d\left(\frac{x^2}{8} \right)+d\left(\frac{y^2}{2} \right)=0; \; \frac{x}{4}dx+ydy=0; \; dy=-\frac{xdx}{4y}. $$
В принципе, здесь можно сразу подставить координаты стационарной точки $x=2$, $y=1$ и параметра $\lambda=-10$, получив при этом:
$$ F_{xx}^{""}=\frac{-5}{2}; \; F_{xy}^{""}=-10; \; dy=-\frac{dx}{2}.\\ d^2 F=F_{xx}^{""}dx^2+2F_{xy}^{""}dxdy+F_{yy}^{""}dy^2=-\frac{5}{2}dx^2+10dx\cdot \left(-\frac{dx}{2} \right)-10\cdot \left(-\frac{dx}{2} \right)^2=\\ =-\frac{5}{2}dx^2-5dx^2-\frac{5}{2}dx^2=-10dx^2. $$
Однако в других задачах на условный экстремум стационарных точек может быть несколько. В таких случаях лучше $d^2F$ представить в общем виде, а потом подставлять в полученное выражение координаты каждой из найденных стационарных точек:
$$ d^2 F=F_{xx}^{""}dx^2+2F_{xy}^{""}dxdy+F_{yy}^{""}dy^2=\frac{\lambda}{4}dx^2+10\cdot dx\cdot \frac{-xdx}{4y} +\lambda\cdot \left(-\frac{xdx}{4y} \right)^2=\\ =\frac{\lambda}{4}dx^2-\frac{5x}{2y}dx^2+\lambda \cdot \frac{x^2dx^2}{16y^2}=\left(\frac{\lambda}{4}-\frac{5x}{2y}+\frac{\lambda \cdot x^2}{16y^2} \right)\cdot dx^2 $$
Подставляя $x=2$, $y=1$, $\lambda=-10$, получим:
$$ d^2 F=\left(\frac{-10}{4}-\frac{10}{2}-\frac{10 \cdot 4}{16} \right)\cdot dx^2=-10dx^2. $$
Так как $d^2F=-10\cdot dx^2 < 0$, то точка $(2;1)$ есть точкой условного максимума функции $z=5xy-4$, причём $z_{\max}=10-4=6$.
Ответ : в точке $(2;1)$ функция имеет условный максимум, $z_{\max}=6$.
В следующей части рассмотрим применение метода Лагранжа для функций большего количества переменных.
Краткая теория
Метод множителей Лагранжа является классическим методом решения задач математического программирования (в частности выпуклого). К сожалению, при практическом применении метода могут встретиться значительные вычислительные трудности, сужающие область его использования. Мы рассматриваем здесь метод Лагранжа главным образом потому, что он является аппаратом, активно используемым для обоснования различных современных численных методов, широко применяемых на практике. Что же касается функции Лагранжа и множителей Лагранжа, то они играют самостоятельную и исключительно важную роль в теории и приложениях не только математического программирования.
Рассмотрим классическую задачу оптимизации:
Среди ограничений этой задачи нет неравенств, нет условий неотрицательности переменных, их дискретности, и функции и непрерывны и имеют частные производные по крайней мере второго порядка.
Классический подход к решению задачи дает систему уравнений (необходимые условия), которым должна удовлетворять точка , доставляющая функции локальный экстремум на множестве точек, удовлетворяющих ограничениям (для задачи выпуклого программирования найденная точка будет одновременно и точкой глобального экстремума).
Предположим, что в точке функция (1) имеет локальный условный экстремум и ранг матрицы равен . Тогда необходимые условия запишутся в виде:
есть функция Лагранжа; – множители Лагранжа.
Существуют также и достаточные условия, при выполнении которых решение системы уравнений (3) определяет точку экстремума функции . Этот вопрос решается на основании исследования знака второго дифференциала функции Лагранжа. Однако достаточные условия представляют главным образом теоретический интерес.
Можно указать следующий порядок решения задачи (1), (2) методом множителей Лагранжа:
1) составить функцию Лагранжа (4);
2) найти частные производные функции Лагранжа по всем переменным и приравнять их
нулю. Тем самым будет получена система (3, состоящая из уравнений. Решить полученную систему (если это окажется возможным!) и найти таким образом все стационарные точки функции Лагранжа;
3) из стационарных точек, взятых без координат выбрать точки, в которых функция имеет условные локальные экстремумы при наличии ограничений (2). Этот выбор осуществляется, например, с применением достаточных условий локального экстремума. Часто исследование упрощается, если использовать конкретные условия задачи.
Пример решения задачи
Условие задачи
Фирма производит товар двух видов в количествах и . Функция полезных издержек определена соотношением . Цены этих товаров на рынке равны и соответственно.
Определить, при каких объемах выпуска достигается максимальная прибыль и чему она равна, если полные издержки не превосходят
Испытываете сложности с пониманием хода решения? На сайте действует услуга Решение задач по методам оптимальных решений на заказ
Решение задачи
Экономико-математическая модель задачи
Функция прибыли:
Ограничения на издержки:
Получаем следующую экономико-математическую модель:
Кроме того, по смыслу задачи
Метод множителей Лагранжа
Составим функцию Лагранжа:
Находим частные производные 1-го порядка:
Составим и решим систему уравнений:
Так как , то
Максимальная прибыль:
Ответ
Таким
образом необходимо выпускать
ед. товара 1-го вида и
ед. товара 2-го вида. При этом прибыль будет
максимальной и составит 270.
Приведен образец решения задачи квадратичного выпуклого программирования графическим методом.
Решение линейной задачи графическим методом
Рассмотрен графический метод решения задачи линейного программирования (ЗЛП) с двумя переменными. На примере задачи приведено подробное описание построения чертежа и нахождения решения.
Модель управления запасами Уилсона
На примере решения задачи рассмотрена основная модель управления запасами (модель Уилсона). Вычислены такие показатели модели как оптимальный размер партии заказа, годовые затраты на хранение, интервал между поставками и точка размещения заказа.
Матрица коэффициентов прямых затрат и матрица "Затраты-выпуск"
На примере решения задачи рассмотрена межотраслевая модель Леонтьева. Показано вычисление матрицы коэффициентов прямых материальных затрат, матрицы «затраты-выпуск», матрицы коэффициентов косвенных затрат, векторов конечного потребления и валового выпуска.
Самые разрушительные цунами современности
Инопланетянин, родившийся на земле
Основатель франции. Конституция франции. Энергетика и добыча полезных ископаемых