Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Математика для инженеров(теория)I том

.pdf
Скачиваний:
128
Добавлен:
18.02.2016
Размер:
4.09 Mб
Скачать

 

 

 

 

 

 

 

 

 

 

 

2 z

=

 

 

(4x3

+ 4xy2 )=12x2 + 4y2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 z

=

 

 

 

(4x2 y + 4y3 )= 4x2 +12y2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 z

 

=

 

 

(4x3 + 4xy2 )= 8xy ,

 

2 z

 

 

=

 

 

(4x2 y + 4y3 )= 8xy .

 

 

 

 

 

xy

y

y x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

Сравнивая значения смешанных частных производных

 

2 z

 

и

 

2 z

 

,

 

 

xy

 

y x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

видим, что они совпадают.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема

1.

 

Если

 

в

 

 

некоторой

 

окрестности

точки

 

(x0 ; y0 ) функция z = f (x, y)

 

имеет смешанные частные про-

 

изводные

 

2 z

,

 

 

 

2 z

 

,

 

и

 

эти производные

непрерывны в

 

xy

 

 

y x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

самой точке (x0; y0 ) ,

то они в этой точке равны.

 

 

 

 

 

 

 

Доказательство.

 

 

Пусть

x

 

и

 

y приращения аргу-

 

ментов такие, что точка

(x0 +

x; y0 +

 

 

y) не выходит из об-

 

ласти

определения

 

 

 

 

функции

 

 

 

 

 

 

f (x, y) .

 

Составим

 

вспомогательное выражение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A =

é f

(x + D x, y

0

 

+ D y) - f (x + Dx, y

0

)ù

- é f

(x , y + D y)

- f (x , y

)ù.

 

 

 

ë

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

û ë

 

0 0

 

 

 

 

 

0 0

û

 

Введя обозначение

F (x) = f (x, y0 +

y) f (x, y0 ) ,

будем

 

и

 

 

 

 

 

м

 

 

 

 

 

 

 

 

 

 

 

 

 

 

е

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

ь

 

 

 

 

 

 

 

 

 

 

 

 

 

A = F (x0 +

x) F (x0 ) .

 

 

 

 

 

 

(1)

 

 

Применив к правой части выражения (1) теорему Ла-

 

гранжа,

получим

 

 

 

 

A = F

(x

 

 

x)

 

 

x

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A = ( fx(x0 1

 

 

x, y0 + y) fx(x0 1 x, y0 ))

 

x ,

 

(2)

 

 

где

0 <θ1 <1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x

В равенстве (2) в скобках стоит приращение функции

 

x, y) при изменении у от y

 

 

до y

+

y .

 

 

 

 

 

 

 

 

x 0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

Применив еще раз теорему Лагранжа к правой части выражения (2),

 

получим

 

 

 

 

A = f ′′

 

(x

 

 

x, y

 

 

y)

 

 

x

y ,

 

 

 

 

 

(3)

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xy

 

 

0

 

1

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

0 <θ1 <1,

 

0 <θ2 <1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

483

Обозначив

Φ(y) = f (x0 +

x, y) − f (x0 , y) ,

представим

вспомогательное выражение A в виде

 

 

A = Φ( y0 +

y) − Φ( y0 ) .

(4)

Действуя таким же образом, как и раньше, после дву- кратного применения теоремы Лагранжа, равенство (4)

приведем к виду

 

A = f ′′

(x

3

x, y

0

4

y)

x y ,

 

(5)

 

 

yx

0

 

 

 

 

 

 

 

 

 

 

где 0

3 <1, 0 <θ4 <1.

 

 

 

 

 

 

 

 

 

 

 

Сравнивая выражения (3) и (4) приходим к равенству

f ′′

(x x, y

2

 

y) = f

′′

(x

3

x, y

4

y) .

xy

0

1

0

 

 

 

 

yx

 

0

 

0

 

Отсюда при

x → 0,

y → 0 , с учетом непрерывности смешанных

производных в точке (x0; y0 ) , будем иметь

 

 

 

 

 

fxy′′ (x0 , y0 ) = fyx′′ (x0 , y0 ) .

 

 

Следствие. Если все последовательно вычисляемые производные являются непрерывными функциями в дан- н о й о б л а с т и , т о а н а л о г и ч н о п о л у ч а ю т с я р а в е н с т в а :

f ′′′

= f ′′′

= f ′′′

;

f ′′′

= f ′′′

= f ′′′

и

т . д .

xxy

xyx

yxx

 

xyy

yxy

yyx

 

 

Итак, при многократном дифференцировании функ-

ции многих переменных порядок дифференцирования не существенен, если получающиеся при этом производные являются непрерывными функциями. Таким образом, число различаю-

щихся частных производных любого порядка k (k n) со- кращается и любая частная производная порядка k может

k z

быть записана в виде xmykm , где m одно из чисел

0,1,2,K,k . Следовательно, число различных частных про-

изводных порядка n у функции двух переменных со всеми

н е п р е р ы в н ы м и в некоторой области частными производными, равно не 2n ,

а n +1 .

§ 12. Полные дифференциалы высших порядков

Полный дифференциал dz функции двух переменных z = f (x, y) есть, в свою очередь, функция переменных х и у.

484

Следовательно, можно определить полный дифференциал этой функции. Он называется полным дифференциалом

второго порядка функции z = f (x, y) и обозначается d 2 z .

Полный дифференциал от полного дифференциала второго порядка называют полным дифференциалом третьего по-

рядка d3z и т. д. Аналогично определяются и полные дифференциалы функций большего числа переменных.

Пусть функция z = f (x, y) имеет непрерывные первые

и вторые частные производные на некотором открытом плоском множестве G . Из непрерывности частных произ-

водных fx и fy следует дифференцируемость самой функции

f (x, y) в каждой точке этого множества G . Таким образом, для в с е х т о ч е к (x; y) G о п р е д е л е н д и ф ф е р е н ц и а л

 

 

 

 

dz = f (x, y) dx +

f (x, y) dy .

(1)

 

 

 

 

x

y

 

Согласно сделанным предположениям, частные про-

изводные

z

и

z

имеют на открытом множестве непрерывные частные

x

y

 

 

 

 

 

п

р

о

 

 

 

и

 

 

з

 

в

 

 

¶ æ

z ö

=

2 z

,

æ

z ö

=

2 z

,

 

 

 

ç

÷

 

2

 

ç

÷

 

 

 

x

 

xy

 

 

x è

x ø

 

 

y è

x ø

 

 

 

о

 

 

д

 

 

 

н

 

 

ы

 

 

е

æ

z

ö

=

2

z

,

æ

z

ö

=

2

z .

 

ç

÷

 

ç

÷

 

 

 

 

xy

 

 

 

 

x è

y ø

 

 

y è

y ø

 

y2

Поэтому, в силу утверждения 2.2,

z

и

z

также

 

x

 

y

 

дифференцируемы на множестве G . Таким образом, диф- ференциал dz , рассматриваемый как функция только пе- ременных x и y , в свою очередь, является дифференци-

руемой на множестве G функцией. Вычислим дифферен- циал от первого дифференциала dz , считая dx и dy фикси-

рованными, а точку (x; y) − принадлежащей области G ,

при этом новое дифференцирование обозначим символом

δ :

485

δ (dz)

æ

ç

è

= δ æ z çè x

2 z δ x

xy

 

 

z

ö

æ

δ

z

ö

æ

δ

z

ö

dx +

 

dy ÷

= ç

 

 

÷dx + ç

 

 

 

÷

y

x

y

 

 

ø

è

 

 

ø

è

 

 

ø

 

2

 

ö

 

 

2

 

 

 

 

 

2

z

 

 

+

 

z δ y ÷dy =

 

z dxδ x +

 

 

 

 

 

 

 

 

 

 

 

 

 

y

2

÷

 

 

x

2

 

 

 

xy

 

ø

 

 

 

 

 

 

dy = æç 2 z çè x2

(dxδ y + δ

 

2

z

ö

δ x +

 

δ y ÷dx +

 

 

 

 

yx

÷

 

ø

xdy) +

2 z

dyδ y.

y2

 

 

 

Отметим, что непрерывность вторых производных была использована не только для того, чтобы во всех рас-

сматриваемых точках существовали дифференциалы δ xz и δ yz , но и для того,

чтобы в процессе вычислений не обращать внимания на порядок

дифференцирования.

Полагая теперь δ x = dx, δ y = dy, получим

d 2 z =

2 z

dx2

+ 2

2 z

dxdy +

2 z

dy2

,

(2)

x

2

xy

y

2

 

 

 

 

 

 

 

 

так называемый,

второй

дифференциал

функции

z = f (x, y).

Аналогичным образом, при непрерывности частных производных соответствующего порядка, можно вычис- лить:

d3z = d (d2 z)

= 3z dx3

+ 3

 

 

3z

 

dx2dy + 3

 

3z

dxdy2

+ 3z dy3,

 

 

 

 

 

 

xy2

 

 

 

 

 

 

x3

 

 

x2y

 

 

 

 

 

 

 

 

y3

 

 

 

 

 

 

d 4 z =

4 z

 

dx4

+ 4

 

4 z

dx3dy +

 

 

 

 

 

 

 

 

 

 

 

 

x

4

 

x3y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+6

4 z

 

dx2dy2

+ 4

 

4 z

 

dxdy3

+

4 z

dy4

 

 

 

 

(3)

 

 

 

x2y2

 

xy3

y4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и т.д.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По индукции можно показать, что

 

 

 

 

 

 

 

 

 

 

n

 

n z

 

n

 

1 n z

 

 

 

 

 

n−1

 

 

2

 

n z

 

 

n−2

 

2

 

d

 

z =

 

dx

 

+ C

 

 

 

 

 

 

 

dx

 

 

dy + C

 

 

 

 

 

dx

 

 

dy

 

+K+

 

xn

 

n xn−1y

 

 

n xn−2y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+Ck

 

 

n z

 

 

 

dxnk dyk +K +

n z dy n .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n xnkyk

 

 

 

 

 

 

 

 

 

 

 

yn

 

 

 

 

 

 

 

486

Полученное выражение символически можно запи-

сать в виде равенства

 

 

 

ön

 

æ

 

 

d n z = ç

 

dx +

 

dy ÷

z, nÎ .

(4)

x

y

è

 

ø

 

 

Эту запись следует понимать так: сумму, стоящую в скобках, надо возвести в степень n по формуле бинома

Ньютона, после чего показатель степеней у x и y счи-

тать указателями порядка производных по х и по у от функции z = f (x, y) .

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

функции

u = f (x, y, z)

ее n-й дифференциал

 

d nu символиче-

ски можно записать в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ön

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d nu =

æ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ç

 

 

dx +

 

 

 

 

dy +

 

 

 

 

dz ÷ u .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

 

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим теперь дифференциалы высших порядков

сложной функции z = f (x, y) ,

где x = x(u,v) ,

y = y(u,v).

 

 

 

 

 

В силу инвариантности формы первого дифференциа-

ла имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

z dx +

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz =

dy ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где dx и dy не постоянные, а функции переменных u

и v:

 

 

 

 

 

 

 

x

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx =

 

d u +

dv ;

 

dy =

 

d u +

dv .

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

v

 

 

 

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

z

 

z

ö

 

 

 

æ z ö

 

 

 

 

z

 

(dx)

 

 

æ

z

ö

 

 

z

d (dy) =

d 2 z = d ç

 

dx +

 

 

 

dy ÷

= d

ç ÷dx +

 

 

 

 

d

+ d

ç

 

÷dy +

 

 

 

x

y

x

 

y

 

 

è

 

ø

 

 

 

è ¶x ø

 

 

 

 

 

 

 

 

 

 

 

è

y

ø

 

 

 

=

2 z

dx

2

+

2 z

 

dxdy +

z

d

2

x

+

 

2 z

 

dxdy +

2 z

 

dy

2

+

z

d

2

y =

x2

 

y x

x

 

xy

y2

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=æç ¶x dx + y dy ö÷2 z + xz d 2 x + yz d 2 y. è ø

487

Сравнивая это выражение с формулой (2) для второго дифферен- циала функции z = f (x, y) , замечаем, что второй дифферен-

циал d 2 z сложной функции свойством инвариантности не обладает. Этот вывод справедлив и для дифференциалов более высоких порядков.

§ 13. Формула Тейлора для функции двух переменных

Допустим, что функция z = f (x, y) в окрестности точки (x0; y0 ) имеет непрерывные частные производные до m-го порядка

включительно. Придадим переменным

x0 и

y0

приращения

x и y так, чтобы точка (x0 + Dx; y0 + Dy)

не выходила за пределы

рассматриваемой окрестности.

 

 

 

Составим функцию

 

 

 

F (t) = f (x0 + t D x, y0 + t D y) ,

 

(1)

которая при постоянных x0 , y0 ,D x

и

y

будет слож-

ной функцией переменной t, определенной в окрестности точки t = 0 и имеющей в этой окрестности непрерывные производные до порядка m включительно. Применив к этой функ- ции формулу Тейлора для функции одной переменной, получим

 

 

¢

 

 

 

F

¢¢

(0)

 

2

 

 

 

F

(m)

(0) n

 

F

(m+1)

(θ ) m+1

F (t) = F (0) + F (0)t +

 

 

 

 

t

 

+K +

 

 

 

 

 

 

 

 

t +

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

2!

 

 

 

m!

 

 

(m +1)!

,

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где 0 <θ <1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из равенства (1) имеем

 

 

f (x + t D x, y + t D y)

 

 

 

 

 

f

 

dx f

 

 

dy

 

 

 

 

 

 

¢

(t) = x

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

0

 

 

 

D x +

 

× dt + y

× dt =

 

 

 

 

 

x

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

f (x0 + t D x, y0 + t D y)

 

 

 

 

æ

 

 

 

 

 

 

ö

 

(x0 + t D x, y0 + t D y)

+

 

 

 

 

 

 

 

D y =

ç

 

D x

+

 

 

Dy

÷ f

 

y

 

 

 

 

 

x

y

 

 

 

 

 

 

 

 

 

 

 

è

 

 

 

 

 

ø

 

 

 

 

 

 

 

.

Далее

488

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

2 f

 

x + t D x, y + t D y

)

 

 

 

 

 

 

 

 

 

 

 

 

F¢¢(t ) =

F¢(t) =

 

 

 

 

 

 

 

( 0

 

 

 

0

 

 

 

 

D x

2 +

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

2 f

(x + t D x, y + tD y)

 

 

 

 

 

 

2 f (x + t D x, y

0

 

+ t Dy)

 

+2

 

 

 

0

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

D xDy +

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

Dy2 =

 

 

 

 

 

 

 

 

xy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

æ

Dx +

Dy

ö2

 

f (x + t Dx, y + t D y) .

 

 

 

 

 

 

 

 

 

 

 

ç

 

 

 

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

 

 

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично найдем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F¢¢¢(t)

= çæ

D x +

 

Dy ÷ö f

(x0 + t D x, y0 + t D y) .

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è x

 

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По индукции можно получить

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(k )(t) = æ

Dx +

Dy

ök

 

 

 

f

 

(x

 

+ t Dx, y

+ t D y),k =1,2,K,m . (3)

 

 

 

÷

 

 

 

 

 

 

 

 

 

 

ç

x

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

t = 0 , при k =1,2,K,m , находим

Полагая в формуле (3)

 

F¢(0) = çæ

 

Dx +

 

 

Dy

÷ö f

(x0 , y0 ) =

f (x0 , y0 )

Dx +

 

f (x0 , y0 )

Dy = df (x0 , y0 );

 

 

y

 

 

 

 

 

 

 

 

è x

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ö2

 

 

 

 

 

 

 

2 f

(x , y

0

)

 

 

 

 

 

 

 

F

¢¢(0) = ç

 

 

 

 

Dx +

 

 

 

Dy

÷ f (x0 , y0 )

=

 

 

 

0

 

 

 

 

Dx2 +

 

 

 

x

y

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

 

 

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 f (x , y )

 

 

 

 

 

 

 

 

 

2 f (x , y

 

)

 

2 = d2 f (x , y

 

);

 

 

+2

 

 

 

 

 

0 0

 

 

 

Dx Dy +

 

 

 

 

 

 

 

 

 

0 0

 

Dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

 

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. . .

 

. . .

 

. . . . .

 

.

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

(m−1) (0) = d m−1 f (x , y

0

) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

При k = m , из уравнения (3), имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(m) (θ t) =

æ

Dx +

Dy

öm f

 

(x +θ t Dx, y +θ t Dy) =

 

 

 

ç

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

y

 

 

 

÷

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= d m f (x

+θ t Dx, y

 

+θ t Dy) ,

где 0 <θ <1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставив найденные значения производных в фор-

мулу (2)

 

при

t =1

 

и

 

 

 

 

учитывая,

что

 

F (0) = f (x0 , y0 ) ,

F (1) = f (x0 + Dx, y0 + Dy) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получаем

489

f(x0 + Dx, y0 + Dy) = f (x0 , y0 ) + æç x Dx + y Dy ö÷ f (x0 , y0 ) +

èø

 

 

+

 

1

æ

 

Dx +

Dy ö2

f (x , y

 

) +K+

 

 

 

ç

 

 

 

0

 

 

 

 

 

2!

x

 

 

 

 

y

÷

0

 

 

 

 

 

 

 

 

è

 

 

 

 

ø

 

 

 

 

(4)

 

 

 

 

 

1

 

 

æ

 

 

 

 

Dy öm−1

 

 

 

 

+

 

 

 

 

 

 

Dx +

f

(x , y

) +

 

 

 

 

 

 

 

 

ç

 

 

 

 

 

 

(m -1)!

x

 

 

 

y

÷

 

 

0

0

 

 

 

è

 

 

 

ø

 

 

 

 

+

1

æ

 

Dx +

 

Dy öm

f (x +θ D x, y +θ Dy),

 

ç

 

 

 

 

 

 

 

 

 

 

 

y

÷

 

 

 

0

 

 

0

 

 

m!è x

 

 

 

 

ø

 

 

 

 

 

 

 

 

где 0 <θ <1,

Dx = x - x0 , Dy = y - y0 .

 

 

 

 

Формула (4)

называется формулой Тейлора для функ-

ции z = f (x, y)

 

с остаточным членом в форме Лагранжа. В

дифференциальной форме эта формула имеет вид

 

 

Df = f (x + Dx, y + Dy) - f (x , y

0

) = d f (x , y

0

) +

1

d

2 f (x , y

0

) +

 

0

 

0

 

0

 

 

0

2

 

0

(5)

 

 

1

 

 

 

 

1

 

 

 

 

 

 

+K +

 

 

dm−1 f (x , y

 

) +

d m f (x +θ D x, y +θ Dy),

 

 

 

 

(m -1)!

0

 

 

 

0

 

m!

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

0 <θ <1.

 

 

 

 

 

 

 

 

 

 

 

 

 

Как и для функции одной переменной, для получения формулы Маклорена в формуле Тейлора (4) нужно поло- жить x0 = 0, y0 = 0, x = x, y = y .

При m =1 из формулы (4) для функций двух перемен- ных получаем формулу конечных приращений:

f (x + Dx, y + Dy) - f (x, y) = fx¢ (x +θ Dx, y +θ Dy)Dx + + fy¢ (x +θ Dx, y +θ Dy)Dy .

Из этой формулы следует, что если внутри некоторой области непрерывные частные производные первого по- рядка тождественно равны нулю, то функция внутри упо- мянутой области сохраняет постоянное значение.

Замечание 1. Можно показать, что для функции z = f (x, y) остаточный член формулы Тейлора в форме Ла-

гранжа

R(x, y) = m1!d m f (x0 +θ Dx, y0 +θ Dy)

490

является при (x; y) ® (x0; y0 ) бесконечно малой вели-

чиной

 

 

 

 

 

более

высокого

порядка

малости

по

сравнению

с

ρm−1 ((x; y),(x0; y0 )) , где

ρ =

 

, то есть мож-

(x - x0 )2 + ( y - y0 )2

но представить остаточный член в форме Пеано:

R(x, y) = o(ρm−1 ).

Пример 1. Записать формулу Тейлора с остаточным

членом в форме Пеано для

функции

f (x, y) = 2xy в

точке

(1;1) при m = 3 .

 

 

 

Решение. Для любых

х, у из

окрестности

точки

(x0; y0 ) имеет место формула Тейлора в дифференциальной

форме (5):

f (x, y) = f (x0 , y0 ) + d f (x0 , y0 ) + 12 d2 f (x0 , y0 ) + o(ρ2 ).

Вычислим

f (1,1) = 2 , d f (1,1) = fx¢ (1,1)dx + fy¢ (1,1)dy =

=y × 2xy ln 2 (1;1) ×(x -1) + x × 2xy ln 2 (1;1) ×( y -1) =

=2ln 2×(x -1) + 2ln 2×( y -1) ,

d

2

¢¢

 

 

2

¢¢

¢¢

2

=

 

f (1,1) = fxx (1,1)dx

 

+ 2 fxy (1,1)dxdy + fyy (1,1)dy

 

= y2 × 2xy ln2 2

 

×(x -1)2 + 2

(2xy ln 2 + xy × 2xy ln2 2)

 

×(x -1)( y -1) +

 

 

 

 

 

 

(1;1)

 

 

 

 

 

(1;1)

 

 

 

 

 

 

 

 

 

+x2 × 2xy ln2 2 (1;1) ( y -1)2 = 2ln2 2 ×(x -1)2 + 2(2ln 2 + 2ln2 2)(x -1)( y -1) +

+2ln2 2 ×( y -1)2 .

Следовательно,

2xy = 2 + 2ln 2×(x -1) + 2ln 2 ×( y -1) + ln2 2 ×(x -1)2 + +2ln 2(1+ ln 2) (x -1)( y -1) + ln2 2 ×( y -1)2 + o(ρ2 ),

где ρ2 = (x -1)2 + ( y -1)2 .

§ 14. Экстремум функции двух переменных

491

10 . Экстремум функции двух переменных. Рассмот-

рим функцию z = f (x, y) , определенную в некоторой облас- ти D. Точка M0 (x0; y0 ) называется точкой максимума (ми- нимума) функции f (x, y) , если существует такая δ - окрестность этой точки, где f (x, y) определена и непре- рывна, и для всех точек M (x; y) из этой окрестности, от- личных от точки M0 (x0; y0 ) , выполняется неравенство f (M ) £ f (M0 ) ( f (M ) ³ f (M0 )) .

Точка M0 называется точкой глобального максимума

(соответственно минимума) функции f (M )

в области

D, если

f (M0 ) ³ f (M )

(соответственно

f (M0 ) £ f (M ) ) для

всех М

из D.

 

 

f (x, y) в

 

Минимум

(максимум)

функции

точке

M0 (x0; y0 )

называется строгим, если f (M0 ) < f (M ) ( f (M0 ) > f (M )) .

 

Например, функция z = 5 - 2x2 - y2 в точке M0 (0;0) достигает

максимума, равного 5.

 

 

 

 

 

 

20 . Необходимое условие экстремума. Имеет место

 

Теорема 1. Если в точке (x0; y0 ) дифференцируемая

функция f (x, y)

достигает экстремума, то все ее первые

частные производные в этой точке равны нулю.

 

Доказательство. Пусть M0 (x0; y0 ) точка экстремума

функции f (x, y) . Предположим, что M0 (x0; y0 )

 

точка мак-

симума, тогда

f (x0 , y0 ) ³ f (x, y) для всех точек M (x; y) , на-

ходящихся в окрестности точки M0 (x0; y0 ) .

Зафиксируем

y = y0 ,

получим

f (x, y0 ) функцию одной переменной х.

Эта

функция, как видно из

указанного

неравенства, при

x = x0

 

имеет

максимум,

поэтому

 

ее

 

производная

f ¢ (x , y

0

) = 0 .

 

 

 

 

 

 

x

0

 

 

 

¢ (x , y

 

) = 0 . Таким

 

Аналогично доказывается, что f

0

образом,

 

 

y

0

 

 

 

 

 

 

 

в точке максимума M0 (x0; y0 )

получаем,

что

 

 

492