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

UML_4256

.pdf
Скачиваний:
21
Добавлен:
01.06.2015
Размер:
6.64 Mб
Скачать

градиент функции в точке M (1,1) эллипса. grad z (1,1) = z = (2,6) . В

этой точке он перпендикулярен эллипсу и определяет направление наискорейшего подъема функции (см. §6). Скорость этого подъема

z = 2 10 . Очевидно, в направлении касательной к эллипсу в этой точке функция не изменяется.

Z

M

Y

X

 

 

 

 

 

 

 

 

 

 

Графиком функции

u =

1

x2

+

1

y2

+

1

z2

является множество

2

9

4

 

 

 

 

 

 

 

 

точек четырехмерного пространства (гиперпараболоид). Пересекая его плоскостью u =1, получим поверхность уровня (эллипсоид, см. рис.)

1

x2

+

1

y

2

+

1

z2

=1.

 

3

 

1

 

1

 

 

 

 

 

grad u 1,

 

,1

= 1,

 

,

 

 

перпендикулярен к

2

9

 

4

2

3

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

эллипсу в точке M 1,

 

,1

и определяет направление наибыстрейшего

2

 

 

 

подъема. Скорость этого подъема z = 76 . В направлениях касательной

плоскости в точке M 1, 32 ,1 функция не изменяется.

§ 10. Производные и дифференциалы высших порядков. Независимость от порядка дифференцирования

Пусть

y = f (x)

– числовая функция n переменных. Если

существует

частная

производная от частной производной, то ее

 

 

321

называют

 

частной

производной

второго

 

порядка и обозначают

 

 

 

f

=

2 f

=

f ′′2 ,

 

f

=

 

2 f

= fx′′x

 

,

 

 

 

 

 

 

xi2

 

 

 

 

 

 

 

 

 

 

 

 

 

x jxi

 

 

xi

xi

 

 

xi

 

x j

xi

 

i

j

 

 

 

 

f

=

 

2 f

= fx′′j xi .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xix j

 

 

 

 

 

 

 

 

 

xi

x j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Последние две производные называются смешанными частными производными второго порядка.

Теорема. Если частные производные fxix j и fx′′j xi существуют в

окрестности точки x0 = (x10, x02,..., x0n) и непрерывны в самой точке x0 , то они равны, то есть

fx′′x

j

(x0) = fx′′x (x0).

(1)

i

j i

 

Доказательство. Так как при вычислении указанных в теореме производных изменяются только две переменные xi и x j , то ради

упрощения записи обозначим

их x и

y , а функцию

f (x10, x02,..., xi,..., x j,..., x0n) через f (x, y) . Тогда следует доказать, что

′′

′′

(1’)

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

Найдем частные приращения:

x f (x0, y0) = f (x0 +∆x, y0) f (x0, y0) ,

y f (x0, y0) = f (x0, y0 + ∆y) f (x0, y0) ,

xy f (x0, y0) = ∆x (y f (x0, y0)) = x ( f (x0, y0 + ∆y) f (x0, y0))=

= f (x0 + ∆x, y0 + ∆y) f (x0, y0 + ∆y) ( f (x0 + ∆x, y0) f (x0, y0)).

(2)

yx f (x0, y0) = ∆y (x f (x0, y0)) = ∆y ( f (x0 + ∆x, y0) f (x0, y0))=

 

= f (x0 + ∆x, y0 + ∆y) f (x0 + ∆x, y0) ( f (x0, y0 + ∆y) f (x0, y0)).

(3)

Сравнивая (2) и (3), видим, что xy f (x0, y0) = ∆yx

f (x0, y0) .

(4)

Применяя дважды теорему Лагранжа, найдем

 

 

xy f (x0, y0) = x ( f (x0, y0 + ∆y) f (x0, y0))=

 

 

 

= ∆x ( fy (x

0,η)y) = ( fy (x0

+ ∆x,η) fy (x0,η))y = fyx (ξ,η)yx,

 

ξ (x0, x0 + ∆x) , η ( y0, y0 + ∆y) .

 

(5)

Аналогично, yx f (x0, y0) = ∆y ( f (x0 + ∆x, y0) f (x0, y0))=

 

′′

 

= ∆y ( fx (ξ

1, y0)x) = ( fx (ξ

1, y0 +∆y) fx (ξ1, y0))x = fxy (ξ1,η1)xy ,

 

ξ1 (x0, x0 + ∆x) , η1 ( y0, y0 + ∆y) .

 

(6)

 

 

322

 

 

Поскольку левые части (5) и (6) равны, то имеем

 

′′

′′

(7)

fyx (ξ,η) =

fxy (ξ1,η1) .

Так как по условию теоремы вторая частная производная непрерывна,

то переходя к пределу в (7) при

x 0

и

y 0, получим (1’).

Теорема доказана.

y = f (x)

n

переменных дифференцируема

в

Пусть

функция

некоторой

области.

Ее

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

 

n

f (x)

d xi

является

dy =

xi

 

 

 

 

 

 

 

i=1

 

 

 

функцией

двух

векторных

переменных

x = (x1, x2,..., xn)

и

dx = (dx1, dx2,..., dxn) .

Если

x

независимая

переменная,

то второй

дифференциал функции находят (если он существует) в

предположении,

 

 

что

 

dx = const .

 

Тогда

имеем

 

n

f

 

 

n

f

n n

2 f

 

 

 

d 2y = d(dy) = d

 

d xi

= d(

 

)d xi = ∑ ∑

 

 

dxidx j .

(8)

 

 

 

 

 

 

xi

 

 

 

xi

i=1 j=1

xix j

 

i=1

 

i=1

 

Функцию A(h, h) =

n

 

 

 

 

 

 

 

aijhih j , h = (h1, h2,..., hn)

 

 

(9)

 

 

 

i, j=1

 

 

 

 

 

= a ji , то

называют квадратичной

формой переменной

h . Если aij

квадратичную форму называют симметричной. Будем считать вторые частные производные непрерывными, тогда (8) – симметричная квадратичная форма дифференциала dx независимой переменной.

Формулу (8) для второго дифференциала удобнее переписать в символической (операторной) форме следующим образом:

 

 

 

 

 

 

 

 

n

 

 

2

 

 

 

 

 

 

 

d 2 f (x) =

 

d xi

f (x) .

(8’)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

= 1

xi

 

 

 

По индукции получим формулу для m-го дифференциала

 

 

 

 

 

 

 

 

 

n

 

m

 

 

 

 

 

 

d m f (x) =

 

d xi

f (x) .

(10)

 

 

 

 

 

 

xi

Пример.

 

 

 

i = 1

 

 

 

 

Найти

первый

и

 

второй

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

функции

u = x y2 z3 .

 

 

 

 

 

 

 

 

 

 

 

 

Решение. du = uxdx + uy dy +uz dz = y2z3dx + 2xyzdy +3x y2z2dz .

 

 

 

2

 

 

 

 

 

 

 

d 2u =

 

 

dx +

 

dy +

 

dz u =

 

 

 

 

 

 

 

x

y

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

323

 

 

 

 

 

2u

 

2u

 

2

 

2u

 

2u

2u

=

 

 

dx2 +

 

 

2

dy

 

+

 

 

dz2 + 2

 

dxdy + 2

 

dxdz +

x

2

y

 

z

2

xy

xz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 2xz3 dy2 + 6x y2zdz2 + 4 yz3dxdy + 6 y2z2dxdz +12xyz2dydz

2 2u dydz = ∂yz

.

Рассмотрим сложную функцию n

переменных ϕ(x) = f (u) , u = u(x) ,

u = (u1,u2,...,um) , x = (x1, x2,..., xn) .

Используя инвариантность формы

первого дифференциала, запишем

 

 

 

 

m

 

n

ui

 

df (u) = ∑

f

du , du =

dx .

 

 

 

 

i

i

 

 

k

i = 1 ui

 

k = 1

uk

Теперь при нахождении второго дифференциала мы не можем считать

dui постоянными, так

как

они

зависят

от x .

Найдем второй

дифференциал сложной функции.

 

 

 

 

 

 

 

 

 

m

 

f

 

 

m

 

f

 

 

f

 

 

d

2 f

(u) = d

 

 

dui

=

d

 

 

dui

+

 

d(dui)

=

 

 

 

 

 

 

 

 

 

 

ui

 

 

 

 

ui

 

 

ui

 

 

 

 

 

i = 1

 

 

 

i = 1

 

 

 

=

m

m

2 f

 

duidu j +

m

f

d 2ui .

 

 

 

 

(11)

∑ ∑

 

 

 

 

 

 

 

 

 

 

 

ui

 

 

 

 

 

i = 1 j = 1 duidu j

 

 

 

i = 1

 

 

 

 

 

 

 

 

Сравнивая (11) с (8), видим, что второй дифференциал не обладает инвариантностью формы.

§ 11. Формула Тейлора

Пусть функция f (x) n переменных определена в некоторой окрестности O(x0 ,δ ) Rn в (открытом шаре радиуса δ ). Пусть f (x) имеет в этой окрестности частные производные до (m +1) -го порядка, которые непрерывны в точке x0 . Пусть h – приращение аргумента такое, что x0 + h O(x0 ,δ ) . Тогда в силу выпуклости шара (см. §2) все точки отрезка x(t) = x0 + ht , t [0,1] принадлежат O(x0 ,δ) .

Рассмотрим сложную функцию

 

ϕ(t) = f (x(t)) ,

x(t) = x0 + ht , t [0,1].

(1)

Из непрерывности частных производных функции

f (x) следует ее

дифференцируемость, а из

дифференцируемости

f (x) следует

дифференцируемость сложной функции (1) в точке t = 0.

Запишем формулу Тейлора для функции ϕ(t) одной переменной.

324

m

1

ϕ(k)(0)tk +

 

 

t

m+1

 

+

 

 

 

ϕ(t) =ϕ(0) +

 

 

 

 

 

 

ϕ(m

1)(θ),

θ (0,1) .

 

 

(m +1)!

k =1 k!

 

 

 

 

 

 

 

 

 

 

По правилу дифференцирования сложной функции (1) имеем

n

f (x)

x

=

 

 

n

 

f (x)

h ,

x = x

0

+ ht .

ϕ '(t) = ∑

 

 

k

 

 

 

 

 

 

 

k = 1 xk

dt

k = 1 xk

 

k

 

 

 

 

 

 

 

 

 

 

Из (1) найдем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ(0) = f (x0) , ϕ(1) = f (x) ,

 

 

 

 

 

а из (3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ '(0) =

n

f (x0 )

h

= df (x0 ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xk

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

k = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

Дифференцируя (3) по t , получим при t = 0

 

 

 

 

 

 

 

 

n

n

2 f (x)

hk h j

= d 2 f (x

0

) .

ϕ''(0) = ∑

xkx j

 

 

 

 

k = 1 j = 1

 

 

 

 

 

 

 

 

Дифференцируя (3), по индукции получим

ϕ(k )(0) = d k f (x0 ) , k =1,2,...,m.

ϕ(m+1)(θ) = d m+1 f (x0 +θh) .

Учитывая (4) и (7) из (2) при t =1, получим

 

 

m

1

 

1

 

 

 

f (x) f (x0 ) =

 

d k f (x0 ) +

 

 

d m+1 f (x0

+θh) .

 

(m +1)!

k =1 k!

 

 

(2)

(3)

(4)

(5)

(6)

(7)

(8)

Равенство (8) называют формулой Тейлора с остаточным членом в форме Лагранжа.

Следствие 1. При m = 0 равенство (8) примет вид

f (x) f (x0 ) = df (x0

n

f (x0 +θh)

dxi .

(9)

+θh) =

xi

 

i=1

 

 

Равенство (9) называют формулой конечных приращений Лагранжа. Следствие 2. Формулу Тейлора (8) можно записать с остаточным

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

m

1

d k f (x0 ) + o(hm )

 

f (x) f (x0 ) =

 

(10)

 

k =1 k!

 

(без доказательства).

325

§ 12. Продолжение функции. Частная производная на границе области. Равномерная сходимость по параметру

Пусть функция f (x) n переменных определена и непрерывна на незамкнутом множестве Ε Rn . Если ее можно доопределить во всех точках множества Ε \ Ε таким образом, что полученная функция будет

непрерывной на замыкании Ε, то эту доопределенную функцию называют непрерывным продолжением функции f (x) . Эту новую

функцию обычно обозначают тем же символом f . Продолженную

функцию называют также расширением данной, а данную функцию по отношению к расширенной называют сужением (см. §4).

Теорема 1. Если функция f (x) равномерно непрерывна на незамкнутом множестве Ε, то ее можно продолжить на Ε \ Ε

единственным образом так, что продолжение будет непрерывной на Ε (без доказательства).

Продолжением пользуются, например, в том случае, когда в

некоторых точках границы области определения функции f (x)

нельзя

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

 

 

Y

A

 

 

Пусть, например, дифференцируемая функция

 

 

B

двух переменных f (x, y) задана в

круге

0

 

 

 

 

 

 

 

 

x2 + y2 1 (см. рис.). Ясно, что в точке A(0,1) не

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

имеет смысла f 'x(0,1) . Но если эта функция в

 

 

 

 

 

 

 

 

 

 

 

 

0

0

X

 

 

 

 

 

 

 

 

 

 

круге x2 + y2 <1 имеет равномерно непрерывные

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

Действительно, пусть в нашем случае функция f (x, y) имеет в

точке B(x0, y0) окружности левую частную производную по x .

 

fx(x0, y0) =

lim

f (x0

+ ∆x, y0) f (x0, y0)

=

 

 

x

 

 

x→−0

 

 

= lim fx(x0 +θx, y0)

def

 

 

= fx(x0, y0) .

 

(1)

x→−0

 

 

 

 

 

(Воспользовались теоремой Лагранжа). В правой части (1) стоит обобщенная производная, а в левой – обычная левая, и они равны.

326

Пусть функция f (x, y)

задана на

множестве

Ε = Χ×Υ, Χ Rn ,

Υ Rm .

Если

вторую

переменную

зафиксировать

y = y0 = ( y0, y0,..., y

0 ) , то

получим

функцию

одной переменной

1

2

m

 

 

 

f (x, y0) ,

x = (x1, x2,..., xn) .

При каждом фиксированном y имеем новую

функцию

переменной x .

Переменная

y

в этом случае

называется

параметром,

а функцию

f (x, y) называют функцией переменной x ,

зависящей от параметра y .

 

 

 

 

 

 

 

В частном случае, когда N = Υ R1 , получим счетное множество

функций

f (x, n) = f n(x) , то есть функциональную последовательность

 

 

1

f 1(x), f 2(x),..., f n(x),....

(2)

 

 

 

1

 

y N ,

 

 

 

 

Например,

z =

 

sin(xy) ,

x R ,

 

функция 2-х

переменных

y

эквивалентна функциональной последовательности

 

 

 

 

 

sin x,

1

sin 2x,

1

sin 3x,...,

1

sin nx,....

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

3

 

 

n

 

Фиксируя

в

(2) x , получим числовую

последовательность, которая

будет сходиться или расходиться. Множество G Χ всех значений x , при которых последовательность (2) сходится, называют областью сходимости функциональной последовательности (2). Таким образом, на множестве G задается функция ϕ(x) , к которой сходится

функциональная последовательность (2), то есть

 

ϕ(x) = lim f n(x) , x G .

 

(3)

 

n→∞

 

 

Равенство (3) можно переписать иначе

 

 

 

f n(x) ϕ(x)

 

< ε , n > N (ε, x) .

 

(3’)

 

 

 

Подчеркнем, что N (ε, x) зависит не только от ε ,

но и от

x . Если

неравенство (3’) выполняется одновременно для всех x G ,

то есть

N (ε) не зависит от x , то последовательность

f n(x) называется

равномерно сходящейся к функции ϕ(x) на множестве G. В этом случае

 

f n(x) ϕ(x)

 

< ε , n > N (ε) , x G .

(3’’)

 

 

Неравенство (3’’) можно записать так:

 

sup

 

f n(x) ϕ(x)

 

< ε , n > N (ε) .

(4)

 

 

 

x G

 

 

 

 

 

327

Вернемся теперь к общему случаю. Пусть y0 – предельная точка

множества Υ (возможна бесконечно удаленная). Если существует предел

 

 

 

lim f (x, y) =ϕ(x) , x G ,

 

(5)

 

 

 

yy0

 

 

 

 

 

то говорят, что функция f (x, y) сходится к функции ϕ(x) в точке

y0 .

Равенство (5) можно переписать так:

 

f (x, y) ϕ(x)

 

< ε при

 

 

 

 

 

 

 

 

 

 

y y0

 

< δ(ε, x) .

 

(5’)

 

 

 

 

 

 

Определение.

Если ε > 0 существует окрестность точки

y0

O( y0,δε) такая,

что

 

f (x, y) ϕ(x)

 

< ε

y Υ ∩O( y0,δε) , x G , то

 

 

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

называется равномерно сходящейся к функции ϕ(x) в

точке y0 .

δ

 

 

 

 

 

 

 

 

 

от x . Если функция

f (x, y)

В этом случае

в (5’) не зависит

равномерно сходится к функции ϕ(x)

в каждой точке y0 Υ ,

то она

называется равномерно сходящейся к функции ϕ(x) на множестве Υ.

Заметим,

что

можно наоборот, считать x

– параметром, а

y

аргументом функции f (x, y) .

 

 

 

 

Теорема

2. Пусть функция f (x, y)

задана на

множестве

Ε = Χ×Υ,

а

функции ϕ(x) и ψ ( y) на

множестве

Χ

и

Υ

соответственно. Если сходимость

 

 

 

 

 

 

lim f (x, y) =ϕ(x)

 

 

 

(6)

 

 

yy0

 

 

 

 

или

 

lim f (x, y) =ψ ( y)

 

 

 

(7)

 

 

 

 

 

 

 

xx0

 

 

 

 

равномерная, то существуют равные между собой повторные пределы:

lim lim f (x, y) = lim lim f (x, y) .

(8)

xx0 yy0

yy0 xx0

 

Доказательство. Пусть, например, сходимость (6) равномерная. Тогда согласно определению равномерной сходимости

 

 

 

 

f (x, y) ϕ(x)

 

< ε y Υ ∩O( y0,δε) , x Χ.

(9)

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

Пусть y1, y2 O( y0,δε), тогда имеем

 

 

 

f (x, y1) f (x, y2)

 

=

 

( f (x, y1) ϕ(x)) + (ϕ(x) f (x, y2))

 

 

 

 

 

 

 

328

 

 

 

 

f (x, y1) ϕ(x)

 

+

 

f (x, y2) ϕ(x)

 

<

ε

+

ε

= ε

 

 

 

 

 

 

 

 

2

2

или

 

f (x, y1) f (x, y2)

 

< ε x Χ.

 

 

 

 

 

 

 

 

Переходя к пределу при x x0 из последнего неравенства получим с учетом (7):

ψ ( y1) ψ ( y2)

< ε y1, y2 O( y0,δε) и x O(x0 ,δε ).

(10)

Согласно критерию Коши (10) означает существование предела функцииψ ( y) в точке y0 , то есть

 

 

 

(7)

 

 

 

 

 

 

 

 

 

 

 

limψ ( y) = lim lim f (x, y) = A.

(11)

 

 

yy0

yy0 xx0

 

При y2 y0 из (10) получим

 

 

 

 

 

 

 

 

 

 

 

ψ ( y1) lim ψ ( y2)

(11)

ψ ( y1) A

 

< ε .

(12)

 

=

 

 

 

 

 

y2y0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Фиксируя y = y1 , перепишем (9) и (7) так:

 

 

 

f (x, y1) ϕ(x)

 

< ε

(9’)

 

 

 

 

 

f (x, y1) ψ ( y1)

 

< ε

(7’)

 

 

 

Тогда имеем

ϕ(x) A = (ϕ(x) f (x, y1)) + ( f (x, y1) ψ ( y1)) + (ψ ( y1) A)

 

(7 ')

 

 

(12)

 

(9)

(7 ')

(7)

 

ϕ(x) f (x, y1)

 

 

 

 

ψ ( y1) A

 

+

f (x, y1) ψ ( y1)

+

< ε

+ ε

+ ε

= 3ε .

 

 

 

 

 

 

 

 

 

 

Последнее означает, что

 

 

 

 

 

 

 

 

 

 

 

 

(6)

 

f (x, y) .

 

 

 

 

 

 

A = limϕ(x) = lim lim

 

 

(13)

 

 

 

xx0

 

xx0 yy0

 

 

 

 

 

Сравнивая (13) и (11) видим, что (8) выполняется, и теорема доказана.

§ 13. Экстремум. Необходимые условия

Пусть скалярная функция f (x) определена на множестве Ε Rn . Точка x0 Ε называется точкой максимума функции f (x) , если существует δ – окрестность O(x0 ,δ) такая, что

329

f (x0 ) f (x) для

x O(x0 ,δ ) ∩Ε.

(1)

Аналогично, если f (x0) f (x) для

 

x O(x0 ,δ ) ∩Ε,

(2)

то x0 – точка минимума функции f (x) .

 

Если неравенства (1) и (2) строгие при

x x0 , то говорят, что

функция достигает в точке x0 строгого максимума или минимума. Точки максимума и минимума называют точками экстремума функции.

Замечание 1. Условия (1) и (2) удобнее записать через

приращение функции в точке x0 :

 

 

 

 

 

f (x0) = f (x) f (x0) 0 .=> max; f (x0) 0 . => min.

(1’, 2’)

Теорема 1 (необходимое условие экстремума). Если функция

f (x) имеет частные

производные

в

точке x0 и

достигает в ней

экстремума, то

f (x0)

= 0 , i =1,2,..., n .

 

 

 

 

xi

 

 

 

 

 

 

 

 

 

 

 

 

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

Зафиксируем

все

составляющие

вектора

x = (x1, x2,..., xn) , исключая xi .

Получим функцию одной переменной

f (x10, x02,..., xi, xi0+1,..., x0n) =ϕ(xi) .

Так

как

функция

f (x)

достигает

экстремума в точке x0 , то функция ϕ(xi) достигает экстремума в точке

x0 . Согласно аналогичной теореме для функции одной переменной, ее

i

f (x) по

производная, то есть частная производная функции

переменной xi , обращается в нуль. Теорема доказана.

f (x0 ) = 0 . А

Следствие 1. В условиях теоремы 1, очевидно, grad

поскольку градиент является полной производной функции n переменных, то необходимое условие существования экстремума можно записать в той же форме, что и для функции одной переменной: f '(x0 ) = 0 , где f '(x0 ) = grad f (x0 ) .

Следствие 2. Если функция f (x) дифференцируема в точке x0 , то необходимое условие экстремума можно записать в виде df (x0 ) = 0 .

Замечание 2. Поскольку функция f (x) может достигать

экстремума в точках, в которых не существует полной производной, то возможные точки экстремума следует искать среди точек, в которых f '(x) обращается в нуль или не существует.

330

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]