UML_4256
.pdf
|
|
|
|
g(b) = a , |
|
f (g( y), y) = 0 y W . |
(5) |
||||||
Комментарий. Теорема утверждает, что функциональное |
|||||||||||||
уравнение |
f (x, y) = 0 , то есть система |
|
|
|
|||||||||
|
f |
1(x1, x2,..., xn, y1, y2,..., ym) = 0, |
|
||||||||||
|
|
2(x1, x2,..., xn, y1, y2,..., ym) = 0, |
|
||||||||||
|
f |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
.................................................. |
|
|||||||||||
|
f |
n |
( |
, |
,..., |
xn |
, y , y |
2 |
,..., y |
m |
) = 0 |
|
|
|
|
|
x1 x2 |
|
1 |
|
x = g( y) , |
непрерывно |
|||||
определяет |
неявную |
|
|
функцию |
|
дифференцируемую в окрестности точки b, если ранг матрицы Якоби
|
∂f |
i |
(a,b) |
равен n . |
|
|
|
|
|
||
|
∂x j |
||||
|
|
|
|
|
|
Доказательство теоремы. Рассмотрим отображение F , |
ставящее |
||||||||||||||
в соответствие точке (x, y) Ε Rn + m точку |
( f ((x, y), y) Rn + m , |
то |
|||||||||||||||
есть |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F(x, y) = ( f (x, y), y) + (z, w) , |
|
(6) |
|
где f (x, y) = z Rn , y = w W Rm . |
|
|
|
||||||||||||||
Напомним, что точка (x, y) = (x1, x2,..., xn, y1, y2,..., ym) Rn + m |
и точка |
||||||||||||||||
( f (x, y), y) = ( f 1, f 2,..., f n, y1, y2,..., ym) Rn + m , |
так что отображение |
||||||||||||||||
F(x, y) действует из Rn+m в Rn+m . |
|
|
|
||||||||||||||
|
|
Поскольку f (x, y) – функция непрерывно дифференцируемая, и |
|||||||||||||||
f (a,b) = 0, то |
|
|
|
||||||||||||||
f (a + h,b + l) = f (a,b) + f '(a,b)(h,l) + o((h,l)) = A(h,l) + o((h,l)), |
(7) |
||||||||||||||||
lim |
o |
((h,l)) |
→ 0 при |
|
h |
|
+ |
|
l |
|
→ 0 . |
|
|
|
|||
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|||||||||||||
|
|
h |
+ |
l |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Тогда |
|
|
|
||||||||||||||
|
|
|
|
|
F(a,b) = ( f (a,b),b) = (0,b) = (0,0,...,0,b1,b2,...,bm) . |
|
(8) |
||||||||||
F(a + h,b +l) − F(a,b) = ( f (a + h,b + l),b + l)(8)−(0,b) = |
|
|
|||||||||||||||
= ( f (a + h,b + l),l)(7)=( A(h,l) + 0((h,l),l)) = ( A(h,l),l) + o((h,l),0) . |
|
(9) |
Из (9) видно, что отображение F(x, y) непрерывно дифференцируемое, а полная производная в точке (a,b) записывается так:
F '(a,b) = ( A(h,l),l) . |
(10) |
351 |
|
Покажем, что линейный оператор (10) взаимно-однозначный, следовательно, обратимый. Для этого достаточно убедиться, что прообраз нулевого вектора является только нулевым вектором (см. теорему 1 §17).
Пусть l = 0, тогда из (10) F '(h,0) = ( A(h,0),0) , а из (4) A(h,0) = 0 .
Ù h = 0, если Rang A = n. Это и означает обратимость оператора
F '(a,b) .
Тогда имеет место теорема 2 §17 об обратной функции. То есть
существуют окрестности |
U |
и V точек |
a и b такие, |
что F(x, y) |
|||
взаимно-однозначно отображает U на V . |
|
|
|
||||
Отображение обратное к отображению (6) имеет вид |
|
|
|||||
|
|
x =ϕ(z, w) , |
y = w. |
|
(11) |
||
При этом из f (x, y) = z следует |
|
|
|
|
|||
|
f (ϕ(z, w), w) = z , (z, w) V . |
|
(12) |
||||
Обозначая через |
W окрестность точки |
b такую, что |
(0, w) W , и |
||||
полагая z = 0 , |
из (12) |
получим |
f (ϕ(0, y), y) = 0 , |
y W . |
Если |
||
обозначить ϕ(0, y) = g( y) , |
то получим |
f (g( y), y) = 0 , |
y W . |
Это |
совпадает с (5). То есть мы доказали существование неявной функции x = g( y) в окрестности W точки y = b. Полагая x = a , y = b из (11)
найдем при z = 0 , a =ϕ(0,b) = g(b) , то есть условия (5) полностью
выполняются.
Единственность g( y) следует из взаимной однозначности
отображения F , а дифференцируемость – из дифференцируемости обратной функции. Теорема доказана.
Так как (5) представляет собой тождество, то дифференцируя его по правилу дифференцирования сложной функции, получим формулу для производной неявно заданной функции
f ' = fx′g '( y) + fy′ = 0 . => g '( y) = −( fx′)−1 fy′. |
(13) |
|||||
Здесь fx′ |
и fy′ |
– |
так |
называемые полные частные производные |
||
отображения f (x, y) . |
|
|
|
|
||
Замечание 1. Производные высших порядков от неявной функции |
||||||
можно получить, дифференцируя (13). |
|
|
||||
Замечание |
2. |
Мы |
считали, что |
функциональное |
уравнение |
|
f (x, y) = 0 |
определяет неявно функцию |
x = g( y) . Если, |
наоборот, |
считать, что функциональное уравнение определяет неявно функцию y =ϕ(x) , то в формулах (13) следует поменять местами x и y . По
352
аналогии с формулой для функции с одной переменной запишем формулу (13) в символической форме
|
|
|
|
|
|
|
|
|
|
|
|
y′x = − |
fx′ |
|
|
|
|
|
|
|
|
(13’) |
|||
f (x, y) = f |
|
|
|
|
|
|
|
fy′ |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
( y1, y2,..., yn, x1, x2,..., xm) = 0 . |
|
|
|
|
|
|
|||||||||||||||||||
|
∂f |
|
|
|
∂f |
1 |
|
|
∂f |
1 |
|
|
|
|
|
∂f |
|
∂ |
∂ |
|
|
||||
|
|
1 |
|
|
|
....... |
|
|
|
|
|
|
|
1 |
|
|
|
f 1 |
...... |
f 1 |
|
|
|||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
∂y1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
∂y2 |
|
∂yn |
|
|
|
∂x1 |
|
∂x2 |
∂xm |
|
|||||||||||||
fy′ = ... |
.............. |
..... |
|
|
, |
fx′ = ..... ............. |
..... |
. |
(14) |
||||||||||||||||
|
∂f n |
|
∂f n ........ |
|
∂f n |
|
|
|
|
∂f n |
|
|
|
∂f n ...... |
∂f n |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
∂y1 |
∂y2 |
|
|
|
|
|
|
|
|
∂x1 |
|
∂x2 |
|
|
|
|||||||||
|
|
∂yn |
|
|
|
|
∂xm |
|
|||||||||||||||||
В случае, |
|
когда |
x, y, f (x, y) R1 |
уравнение |
f (x, y) = 0 |
может |
определять только числовую функцию одной переменной y(x) (или x( y) ). Производная неявной функции есть отношение частных производных функции f (x, y) , взятое со знаком минус.
Пример 1. x |
y |
− y |
2 x |
= |
′ |
|
|
|
|
|
|
0 . Найти y (1). |
|
|
|||||
Решение. Проверим выполнение требований теоремы. Во-первых, |
|||||||||
уравнение |
f (1, y0) = 0 |
|
должно иметь |
решение. |
1 = y2 . |
=> y1,2 = ±1. |
|||
Требование |
выполняется. |
Во-вторых, |
функция |
f (x, y) |
должна быть |
непрерывно дифференцируемой в окрестности точки (x0, y0) . Так как
частные производные |
∂f |
= yxy−1 − 2ln y y2 x , |
∂f |
= x y ln x − 2x y2 x−1 |
|
∂x |
|
∂y |
|
непрерывны в окрестности точки (1,1) , то это требование выполняется (см. теорему 1 §18). В-третьих, якобиан в точке (1,1) не должен
обращаться в нуль. В |
нашем |
случае это ∂f (1,1) |
|
= −2 ≠ 0 . |
Итак, все |
||||||
|
|
|
|
|
|
∂x |
|
|
|
|
|
требования теоремы выполняются. Воспользуемся формулой (13’). |
|||||||||||
|
fx′ |
|
∂f |
|
|
yx y −1−2 ln y y2 x |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
||||
y′x = − |
= − |
∂x |
∂f |
= − |
|
= |
. |
|
|||
|
x y ln x−2x y2 x−1 |
|
|
||||||||
|
fy′ |
|
|
|
2 |
|
|
||||
|
|
|
|
∂y |
|
|
|
(1,1) |
|
|
|
Упражнение. Найти y′′2 в точке (1,1) . |
|
|
|
|
|
||||||
|
|
|
|
x |
|
|
|
|
|
|
|
В случае, когда |
x Rn , |
y, f (x, y) R1 , уравнение |
f (x, y) = 0 |
может определять только числовую функцию n переменных, например, 353
y =ϕ(x1, x2,..., xn) . |
Знаменатель |
(якобиан) в |
(13’) |
– это |
частная |
||||||||
производная |
∂f |
, |
а числитель |
– вектор |
f ' |
|
= ( |
∂f |
, |
∂f |
,..., |
∂f |
) . |
∂y |
x |
|
∂x2 |
|
|||||||||
|
|
|
|
|
|
∂x1 |
|
∂xn |
Производная y′x неявной функции – это полная производная числовой
|
|
∂ϕ |
|
∂ϕ |
|
|
∂ϕ |
∂ϕ |
|
∂f |
|
|
функции ϕ(x), |
|
|
|
|
∂xi ∂f . |
|||||||
то есть вектор y′x = |
|
, |
|
,..., |
|
, |
|
= − |
||||
∂x1 |
∂x2 |
|
∂xi |
|||||||||
|
|
|
|
|
∂xn |
|
|
∂y |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Пример 2. Найти частные производные z′x и z′y функции, |
||||||||||||
заданной неявно уравнением z3 +3xyz −1 = 0 . |
|
|
|
|
|
|
|
|||||
Решение. |
Воспользуемся (13’), считая |
все |
требования теоремы |
выполненными в окрестности некоторой точки (x, y) . Неявная функция
z = z(x, y) |
|
|
– |
|
числовая |
|
|
|
функция |
|
двух переменных, |
z ' = (z′x , z′y ), |
|||||||||||||||||||
|
|
|
∂f |
|
|
|
|
|
yz |
|
|
|
|
|
|
|
|
∂f |
|
|
|
|
xz . |
|
|
|
|||||
z′x = − |
∂x |
∂f |
|
= − |
|
, |
|
|
z′y = − |
∂y |
∂f |
= − |
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
z2 + xy |
|
|
|
|
|
|
|
|
|
|
z2 + xy |
|
|
|
|||||||||
|
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
|
|||
|
В случае, |
когда |
x R1, |
|
а |
y, f (x, y) Rn , уравнение |
f (x, y) = 0 |
||||||||||||||||||||||||
определяет |
|
неявно |
|
векторную |
функцию |
скалярного |
аргумента |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
fx′ = |
|
∂f |
|
|
∂f |
2 |
|
|
∂f |
|
|
|
|
|
|
|
|
|||
y(x) = ( y1, y2,..., yn) . |
|
|
1 |
, |
|
,..., |
|
n |
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂x |
|
∂x |
|
∂x |
|
|
|
|
|
|
|||||||
|
|
∂f |
1 |
|
|
|
∂f |
1 |
....... |
|
∂f |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
fy′ |
|
∂y1 |
|
|
∂y2 |
|
|
∂yn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
= ... .............. ..... |
. |
|
|
|
Производная |
неявной |
функции |
||||||||||||||||||||||||
|
|
|
|
|
|
∂f n |
|
|
∂f n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
∂f n |
|
........ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
∂y1 |
|
|
∂y2 |
|
∂yn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
y′x |
= −( fy′)−1 fx′ – вектор. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
Пример 3. Векторная |
|
функция |
|
(x) |
скалярного |
аргумента x |
||||||||||||||||||||||||
|
|
y |
|||||||||||||||||||||||||||||
задана |
неявно |
системой |
уравнений |
|
x2 + y2 + z2 =1, |
y = 2x. |
Найти |
||||||||||||||||||||||||
производную этой функции. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
Решение. |
В |
данном |
|
|
случае |
f (x, |
|
) = ( f 1, f 2) , |
где |
|||||||||||||||||||||
|
|
|
y |
||||||||||||||||||||||||||||
f 1 = x2 + y2 + z2 −1, f 2 = y − 2x . |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
354 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
1 |
|
|
|
|
|
|
достигает |
в точке |
(x0, y0) = |
|
, |
|
|
|
условного |
минимума, |
равного |
|||
3 |
3 |
||||||||||||
|
2 |
|
|
|
|
|
|
|
|
||||
zmin = |
. |
Заметим, |
что безусловный |
минимум |
достигается |
в начале |
|||||||
3 |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
координат и равен нулю.
Рассмотрим задачу в общем случае. Пусть u = f (x, y) ,
(x, y) = (x1, x2,..., xn, y1, y2,..., ym) Ε Rn+m – числовая функция (n + m)
переменных |
непрерывно |
дифференцируема в |
области Ε. |
И пусть |
ϕ(x, y) = 0 – |
уравнения |
связей, а z =ϕ(x, y) |
– векторная |
функция |
(n + m) переменных, отображающая область Ε в Rm .
Условие ϕ(x, y) = 0 выделяет из области Ε некоторую подобласть,
которую условно назовем поверхностью в пространстве Rn+m по аналогии с поверхностью уровня для числовой функции (см. §9).
Пусть |
ϕ(a,b) = 0 , |
то есть точка (a,b) |
лежит на поверхности. |
Говорят, что функция |
f (x, y) достигает в |
точке (a,b) условного |
|
минимума, |
если неравенство f (x, y) ≥ f (a,b) |
выполняется для всех |
точек поверхности, достаточно близких к точке (a,b) . При этом не требуется выполнения неравенства f (x, y) ≥ f (a,b) для точек хотя и близких к точке (a,b) , но не лежащих на поверхности.
Аналогично определяется условный максимум. Условный максимум и минимум называют условным экстремумом.
Рассмотрим два метода нахождения условного экстремума.
Первый метод. |
Пусть |
функция z =ϕ(x, y) непрерывно |
дифференцируема в Ε, |
ϕ(a,b) = 0 |
и матрица Якоби имеет ранг, равный |
m. Тогда функция удовлетворяет теореме о неявной функции (см. §19), то есть в некоторой окрестности точки a Rn существует непрерывно
дифференцируемая векторная функция y = y(x) |
со значениями в Rm , |
|||
при этом ϕ(a) = b . |
Подставим y = y(x) |
в |
функцию |
u = f (x, y) . |
Получим сложную функцию f (x, y(x)) n |
переменных, |
непрерывно |
||
дифференцируемую |
в окрестности точки |
x = a . Ее |
безусловный |
экстремум в точке x = a и будет условным экстремумом для функции f (x, y) в точке (a,b) . Именно этим методом мы решили задачу в
примере 1.
Второй метод (метод множителей Лагранжа).
Первый метод простой в теоретическом отношении может вызвать затруднения на практике, поскольку не всякую неявную функцию
357
|
n |
∂L |
|
|
|
∂L |
= 0 , так как xk – независимые переменные. |
|||
dL = ∑ |
|
|
dxk = 0. => |
|
||||||
∂xk |
∂xk |
|||||||||
k =1 |
|
|
|
|
|
|||||
Присоединяя к необходимым условиям экстремума функции |
||||||||||
Лагранжа уравнения связи, получим систему n + 2m уравнений: |
|
|||||||||
|
∂L |
= 0 , |
∂L |
= 0 , ϕ j(x, y) = 0 , k =1,2,..., n , j =1,2,..., m , |
(7) |
|||||
|
|
|
|
|||||||
|
∂xk |
|
∂y j |
|
|
|
из которых найдем координаты точки (a,b) возможного экстремума и множители Лагранжа λ j .
Из (5) видно, что при наличии связей, то есть при ϕ(x, y) = 0 , функция Лагранжа совпадает с данной f (x, y) . Поэтому безусловный экстремум функции Лагранжа совпадает с условным для функции
f(x, y) . Отсюда следует алгоритм решения задачи:
1)составляем функцию Лагранжа (5);
2)записываем систему (7), из которой находим точку (a,b) возможного условного экстремума;
3)считая все переменные независимыми, находим второй дифференциал функции Лагранжа
|
n |
∂ |
|
|
m |
∂ |
|
2 |
|
d 2 L = |
∑ |
|
dxi |
+ |
∑ |
|
dy |
L(x, y) ; |
(8) |
|
|
||||||||
|
|
∂xi |
|
|
j = 1 ∂y j |
|
j |
|
|
i = 1 |
|
|
|
|
dy j |
||||
4) используя |
(2), |
дифференциалы зависимых переменных |
выражаем через дифференциалы dxk независимых переменных, а затем исследуем знак квадратичной формы (8).
Пример 2. На кривой y = x2 , x + z = 4 найти точку, ближайшую к
началу координат (см. рис.)
Решение. Пусть искомая точка M (x, y, z) .Тогда квадрат расстояния
|
|
|
|
|
|
|
|
|
|
|
Z |
ее от начала |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
координат ρ2(M ,0) = f (x, y, z) = x2 + y2 + z2 . |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Итак, следует исследовать на условный |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
экстремум функцию f (x, y, z) = x2 + y2 + z2 с |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
двумя связями |
ϕ1 = y − x2 |
и ϕ2 = x + z − 4. |
|
|
|
|
|
|
|
|
M |
Составим |
функцию |
Лагранжа |
|||||||
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
L = x2 + y2 + z2 + λ1( y − x2 ) |
+ λ2(x + z − 4) . |
|||||||||
|
|
|
|
0 |
|||||||||||||
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Запишем необходимые условия (7): |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
Y |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L'x = 2x − 2λ1x + λ1 = 0 , |
L'y = 2 y + λ1 = 0 , |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L'z = 2z + λ2 , y − x2 = 0, x + z − 4 = 0. |
|||
X |
|||||||||||||||||
Решая эту систему, найдем единственное |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
359 |
|
|
решение x0 =1, y0 =1, z0 = 3, λ1 = −2 , λ2 = −6 . Запишем второй дифференциал d 2 L = 6dx2 + 2dy2 + 2dz2 .
Поскольку квадратичная форма положительно определенная, то выражать все через дифференциалы независимой переменной нет смысла.
Ответ: M 0(1,1,3), f min =11.
360