Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Кравченко Н.И. УП_М_Ч2.pdf
Скачиваний:
11
Добавлен:
13.03.2016
Размер:
1.11 Mб
Скачать

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

90

 

 

 

 

 

 

 

d

 

d 2

 

 

 

d 3

 

 

 

 

 

 

 

 

 

r

r

r

 

 

 

 

 

 

 

 

 

 

 

2 ×

 

 

 

3

 

1

 

dt

 

 

 

dt

 

 

dt

 

 

 

 

 

 

 

Если

r

=

r

(t) , то

 

= −

 

 

 

 

 

 

 

 

.

Т

d

 

 

 

d 2

 

2

 

 

 

 

 

 

r

×

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

dt

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 28. Вычислить кручение винтовой линии

 

 

 

 

=

 

 

аcos t +

 

аsin t +

 

 

аmt

 

 

 

 

 

 

 

 

 

 

 

 

r

i

j

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

аsin t

 

 

аcos t

am

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d r

d 2 r

 

d 3 r

 

 

= а3m .

Решение.

×

=

аcos t

 

аsin t

0n

 

 

 

 

 

 

 

2

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

dt

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

аsin t

 

аcos t

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

× d

2

 

 

2

= а4 (1 + m2 ); Т = − а

4

(1 + m

2

)= −

а(1 + m

2

).

 

 

 

d r

 

 

r

 

 

 

 

 

dt

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

3

n

 

 

 

m

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13 ФУНКЦИИ НЕСКОЛЬКИХ ПЕРЕМЕННЫХ

13.1 Основные понятия

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

Пример 29. Площадь S прямоугольника со сторонами, длины которых равны х и у вычисляется по формуле: S = х · у,

где S – является функцией двух переменных т.к. каждой паре значений х и у соответствует определенное значение площади S.

Пример 30. Объем V прямоугольного параллелепипеда с ребрами, длины которых равны х, у, z определяется по формуле V = x y z .

Здесь V – функция трех переменных x y z.

Определение. Функцией n переменных х1, х2,…,хn, где (х1, х2,…,хn) D Rn будем называть правило или закон, по которому каждому

набору переменных (х1, х2,…,хn) D ставится в соответст-

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

образом: у = f(х1, х2,…,хn) .

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

Определение. Функцией двух переменных х, у будем называть правило или закон, по которому каждой паре чисел (х, у) D ставится в соответствие единственное число z Е.

91

Тот факт, что задана функция двух переменных, будем записывать в виде: z = f(x, у). При этом х и у будем называть независимыми переменными (аргументами), а z – зависимой переменной (функцией).

Множество D(z) называется областью определения функции. Множество Е значений, принимаемых z в области определения, называется областью изменения функции.

Функцию двух переменных z = f(х, у), где (х, у) D можно рассматривать как функцию точки М(х, у) координатной плоскости Оху.

Значение функции z = f(х, у) в точке М0(х0, у0) называют частным значением

функции и обозначают одним из способов: z(х0, у0), z x = x0 , f(х0, у0).

y = y0

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

Функция двух переменных, как и функция одной переменной, может быть задана разными способами.

1)Аналитический способ состоит в том, что функция z представлена

спомощью формулы. Если при этом область определения D(z) не указана, то под ней понимают множество таких пар значений (х, у), при которых заданная формула имеет смысл. Областью определения может быть вся плоскость или её часть, ограниченная некоторыми линиями. Линию, ограничивающую область, называют границей области. Точки области, не лежащие

у

 

х 2 - у2< 4

на границе, называются внутренними.

 

 

Область, состоящая только из внутренних то-

 

 

 

 

 

 

чек, называется открытой, область с присое-

 

 

 

 

диненной к ней границей называется замкну-

0

 

2 х

 

той областью. Например, функция

 

 

 

 

 

 

 

 

z = ln(4 – х2 у2) имеет областью определения

Рисунок 62

внутреннюю часть круга х2 + у2 < 4 (открытая

область, рисунок 62).

2) Табличный способ. В первой строке таблицы выписываются возможные значения переменной х: х1, х2, …, хn, в первом столбце – значения переменной у: у1, у2, …, уn, на пересечении строки и столбца указывается соответствующее значение функции, например: z11 = z(x1, y1), zij = z(xi, yj) (i = 1, 2, …, n; j = 1, 2, …, n).

 

 

 

 

 

 

 

у / х

х1

х2

 

xi

 

xn

y1

z11

z21

 

zm1

 

zn1

y2

z12

z22

 

zm2

 

zn2

 

 

 

уi

z1j

z2j

 

zij

 

znj

 

 

 

 

 

 

 

yn

z1n

z2n

 

zin

 

znn

92

3) Графический способ. Примем z за аппликату некоторой точки P(х, у, z) в пространстве. Всей области D(z) соответствует множество точек Р, образующее в пространстве некоторую поверхность. Графиком функции z = f(х, у) в прямоугольной декартовой системе координат в пространстве является поверхность, каждая точка которой Р(х, у, f(х, у)).

Если поверхность является графиком функции двух переменных, то уравнение, определяющее эту функцию, является уравнением поверхности.

Например, функция z = 9 х2 у2 имеет областью определения

круг х2 + у2 9 (рисунок 63) и изображается верхней полусферой с центром в точке О(0, 0, 0) и радиусом R = 3 (рисунок 64).

z

 

 

z

 

 

 

 

3

 

Д(Z)

 

 

 

3

0

3

х

0

у

 

 

 

 

 

 

х

 

Рисунок 63

 

Рисунок 64

 

13.3 Предел функции

Понятия предела функции двух (и более) переменных и непрерывности вводится аналогично понятию предела и непрерывности функции одной переменной.

Определение. Расстоянием от точки М1(х1, у1) до точки М2(х2, у2) назовем

число: ρ = ρ (М

1

, М

2

) =

(х

2

х )2

+ ( у

2

у )2 .

 

 

 

 

1

 

1

Определение. δ-окрестностью точки М0(х0, у0) назовем множество всех точек М(х, у) плоскости, таких, что ρ(М, М0) < δ.

Геометрически δ-окрестность точки М0

у

 

 

 

 

 

 

 

 

 

 

это все внутренние точки круга с центром М0 и

 

 

δ

радиусом δ (рисунок 65).

 

 

М М0

Пусть функция z = f(х, у) определена в не-

 

 

которой окрестности точки М0(х0, у0), кроме быть

 

 

 

 

0

 

х

может, самой этой точки.

 

 

Рисунок 65 Определение. Число b называется пределом функции z = f(х, у) в точке

М0(х0, у0), если для любого сколь угодно малого положительного числа ε > 0 существует δ > 0 такое, что для всех точек М(х, у) D(z), отличных от точки М0 и удовлетворяющих неравенству ρ(М, М0) < δ выполняется неравенство

f(х, у) – b< ε .

 

93

Записывается так b = lim

f (М) или b = lim f (х, у) .

ММ0

хх0

 

уу0

Из определения следует, что если предел существует, то он не зависит от пути, по которому точка М стремится к точке М0.

Геометрический смысл предела функции двух переменных состоит в следующем. Каково бы ни было число ε > 0, найдется δ-окрестность точки М0, что во всех точках М(х, у), отличных от М0, аппликаты соответствующих точек на поверхности z = f(х, у) отличаются от числа b по модулю меньше, чем на ε .

Отметим, что практически все теоремы о пределах, приведенные для функций одной переменной остаются справедливыми и для функции нескольких переменных.

Пример 31. Найти предел lim

 

х

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2х + у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Будем приближаться к О(0, 0): а) по оси Ох

lim

 

х

 

 

= 1 ;

 

 

2х + у

 

 

у

 

 

 

 

 

 

 

 

 

 

 

 

у=0

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х0

 

 

 

 

 

 

 

 

 

 

 

б) по оси Оу:

lim

 

х

= 0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

 

 

 

 

 

х=0 2х + у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с) по прямой у = kх:

 

 

 

х

=

 

1

 

 

Рисунок 66

lim

 

 

 

.

 

 

+ kx

 

2 + k

 

 

 

 

 

 

 

 

 

 

 

 

у=kx 2х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х0

 

 

 

 

 

 

 

 

Отсюда следует, что функция у =

 

х

 

в точке О(0, 0) предела не име-

 

2х + у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ет, т.к. результат зависит от способа приближения (х, у) к началу координат (рисунок 66).

13.4 Непрерывность функции двух переменных

Определение. Функция z = f(х, у) называется непрерывной в точке М0(х0, у0) D(z), если её предел в этой точке совпадает со значением функции в данной точке,

т.е. lim f (х, у) = f (x0 , у0 ) .

хх0 уу0

Определение. Будем говорить, что функция z = f(х, у) непрерывна на множестве, если она непрерывна в каждой точке этого множества.

Определение. Точки, в которых нарушается непрерывность (т.е. не выполняется хотя бы одно из условий непрерывности функции в точке) называются точками разрыва этой функции.

94

Точки разрыва функции z = f(х, у) могут образовывать целые линии разрыва. Например, функция у = х +х у имеет линию разрыва у = – х.

Можно доказать, что арифметические операции над непрерывными функциями двух переменных приводят к непрерывным функциям – подобная теорема имеет место и для функции одной переменной.

13.5 Частные и полное приращения функции двух переменных

Пусть задана функция z = f(х, у). Так как х и у – независимые переменные, то одна из них может изменяться, а другая сохранять постоянное значение. Дадим переменной х приращение ∆х, сохраняя значение переменной у неизменным. Тогда функция z получит приращение, которое назовем частным приращением z по х и обозначим ∆х z: ∆х z = f(x + ∆x, y) – f(х, у).

Аналогично получаем частное приращение z по у:

у z = f(x, у + ∆y) – f(х, у).

Наконец, если аргументу х дать приращение ∆х, а аргументу у – приращение ∆у, то получим полное приращение функции z: ∆z=f(x+∆x, y+∆у)–f(х, у).

Надо заметить, что полное приращение функции, вообще говоря, не равно сумме частных её приращений, т.е. ∆z ≠ ∆х z + у z.

Геометрически полное приращение функции ∆z равно приращению аппликаты графика функции z = f(х, у) при переходе от точки М(х, у) к точ-

ке М1 (х + ∆х, у + ∆у)

(рисунок 67).

 

 

 

 

z

z(M1)

z=f(x, y)

 

z(M)

 

 

 

 

 

z

 

 

0

у

 

у+ у

 

х

 

 

у

х+ х

М(х, у)

 

 

 

 

 

х, у+ у)

х

М1( х+

Рисунок 67

Пример 32. Найти частные и полное приращения функции z = ху.

Решение. х z = (x + ∆x)·у – ху = ху + у·∆х – ху = у·∆х. у z = х(у + ∆у) – ху = ху + х·∆у – ху = х·∆у.

z=(x + ∆x) ( y + ∆у) – ху=ху + у·∆х + х·∆у+∆х·∆у – ху= у·∆х + х·∆у + ∆х·∆у.

Получили, что ∆z≠ х z + у z.

95

13.6Частные производные первого порядка функции двух переменных и их геометрическая интерпретация

Определение. Частной производной по переменной х от функции z = f(х, у) называется предел отношения частного приращения функции ∆х z по переменной х к приращению ∆х

при стремлении ∆х к нулю (если этот предел существует). Частная производная по х от функции z = f(х, у) обозначается одним из сим-

волов:

z

,

f

, z

х

,

f

(х, у) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

х

 

 

х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

хz

 

 

 

f (x + х,

у) f (х,

у)

 

Тогда по определению:

=

lim

=

lim

.

х

х

 

х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х0

 

 

х0

 

 

 

 

Аналогично частная производная по у от функции z = f(х, у) определя-

ется как предел отношения частного приращения функции ∆у z по пере-

 

 

менной у к вызвавшему его приращению ∆у при ∆у → 0.

 

 

 

Обозначения:

Z

,

f

, Z

,

 

f (х, у) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у

 

у

 

у

 

у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом,

z

= lim

 

уz

=

lim

 

f

(x, у + у) f

(х, у)

.

 

 

у

 

у

 

 

 

у

 

 

 

 

 

 

 

 

 

 

 

 

 

у0

 

 

у0

 

 

 

 

 

 

 

Заметим, что частная производная функции двух переменных определяется, как производная функции одной из переменных при условии постоянства значения другой независимой переменной. Поэтому частные производные функции z = f(х, у) находят по формулам и правилам вычисления производных функции одной переменной и только требуется каждый раз помнить, по какой переменной ищется производная.

Выясним геометрический смысл частных производных функции двух переменных. Известно, что графиком функции z = f(х, у) является некоторая поверхность (рисунок 68). Положив у = у0, получим функцию z = f(х, у0), график которой есть линия пересечения поверхности z = f(х, у) с плоско-

стью у = у0.

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

переменной заключаем, что f (x0 , у0 ) = tgα , где α – угол между положи-

x

тельным направлением оси Ох и касательной, проведенной к линии пересечения поверхности z = f(х, у) и плоскости у = у0 в точке Р0(х0, у0, z0).

Аналогично, f(х0 0 ) = tgβ, где β – угол наклона касательной, про-

у

веденной в точке Р0(х0, у0, z0) к линии пересечения поверхности z = f(х, у) и плоскости х = х0, который она образует с положительным направлением оси Оу (рисунок 68).

z
х z
у

 

96

z

z=f(x, y0 )

z=f(x, y )

 

 

 

Р0

 

 

 

0

у0

у

 

 

β

 

 

 

х

x0

 

М0

 

α

 

 

 

 

 

 

Рисунок 68

Пример 33. Найти частные производные функции z = x2 · sin(2x-3y).

Решение. = (х2 sin(2x 3y))х = 2х sin(2x 3y) + 2x2сos(2x 3y) = (х2 sin(2x 3y))у = −3x2сos(2x 3y) .

13.7 Полный дифференциал функции z = f(х, у)

Пусть функция z = f(х, у) определена в некоторой окрестности точки М(х, у). Составим полное приращение функции в точке М:

z = f(x + ∆x, y + ∆у) – f(х, у).

Определение. Функция z = f(х, у) называется дифференцируемой (т.е. имеет полный дифференциал) в точке М(х, у), если её полное приращение можно представить в виде:

∆z = A · ∆x + B · ∆y + α · ∆x + β · ∆y ,

(13.1)

где α = α(∆х, ∆у) → 0 и β = β(∆х, ∆у) → 0 при ∆х → 0,

∆у → 0.

Сумма первых двух слагаемых в равенстве (13.1) есть выражение линейное относительно ∆х и ∆у и представляет собой главную часть приращения функции. Сумма двух последних слагаемых правой части равенства (13.1) является бесконечно малой высшего порядка относительно

ρ = х2 +

у2 , т.е. α х 0 или

ρ 0 .

 

ρ

 

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

Главная часть полного приращения функции z = f(х, у),

 

линейная относительно ∆х и ∆у, называется полным

 

дифференциалом этой функции и обозначается симво-

лом dz: dz = А · ∆х + В · ∆у. (13.2)

Выражения А · ∆х и В · ∆у называют частными дифференциалами.

Для независимых переменных ∆х = dх

и ∆у = dу.

Поэтому равенство (13.2) можно записать в виде:

dz = А · dх + В · dу.

(13.3)

97

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

Если функция z = f(х, у) имеет в точке М(х, у) полный дифференциал, то она

имеет в этой точке частные производные и при этом

z

= А,

z

= В.

х

у

 

 

 

Доказательство. Так как функция дифференцируема в точке М, то справед-

ливо равенство (13.1). Отсюда вытекает, что

lim

z = 0 .

 

 

 

 

 

 

 

 

х0

 

 

 

 

 

 

 

 

 

 

 

у0

 

 

 

 

Положив ∆у = 0, ∆х ≠ 0 в равенстве (13.1), получим:

хz = A · ∆x + α · ∆х.

Отсюда находим, что

хz

 

= А+ α . Переходя в этом равенстве к пределу

х

 

 

 

хz

 

z

 

 

 

 

 

 

 

 

при ∆х → 0, получим

lim

 

= А, т.е.

=

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

 

х

х

 

х0

 

 

 

 

 

z

 

 

 

 

что в точке М существует частная производная

 

= В.

 

у

 

 

 

 

 

 

 

 

 

 

 

 

 

Равенство (13.1) можно записать в виде

z =

z

 

х +

z

 

у + γ , где

х

 

у

 

 

 

 

 

 

 

 

 

 

 

γ = α х + β у 0

при ∆х → 0, ∆у → 0.

 

 

 

 

 

 

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

Из теоремы следует формула для вычисления полного дифференциа-

ла. Формула (13.3) принимает вид:

 

 

 

dz =

z

+

z

.

(13.4)

 

 

 

х

у

 

Теорема 10 (достаточное условие дифференцируемости функции).

 

Если функция z = f(х, у) имеет в точке непрерывные частные производные

 

z

и

z

, то она дифференцируема в этой точке и её полный дифференци-

 

 

 

 

х

у

 

ал выражается формулой (13.4).

 

Примем эту теорему без доказательства. Правила и формулы вычисления дифференциалов функции одной переменной справедливы и для дифференциалов функции двух (и большего числа) переменных.

13.8Применение полного дифференциала к приближенным вычислениям

Из определения дифференциала функции z = f(х, у) следует, что при достаточно малых |∆х| и |∆у| имеет место приближенное равенство

∆z dz.

(13.5)

98

Так как полное приращение функции равно

z = f(x + ∆x, y + ∆у) – f(х, у), то равенство (13.5) можно переписать в виде

f(x + ∆x, y + ∆у) f(х, у) +

z(х, у)

х +

z(х, у)

у. (13.6)

х

у

 

 

 

Полученной формулой пользуются для приближенных вычислений

значений функции двух переменных в точке М(х +

х; у +

у), близкой к

точке М(х; у), если известны значения функции и ее производных в точке М(х; у).

Пример 34. Вычислить приближенно а =

3,972 + 3,072 .

 

 

 

 

 

 

 

Решение. Рассмотрим функцию z =

х2 + у2 . Данное число а есть нара-

щенное значение функции в точке М0(4, 3) при ∆х = -0,03, ∆у = 0,07.

 

 

z(М0) = 42 + 32 = 5.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

х

2

+ у

2

 

 

 

 

х

 

 

Найдем частные производные функции:

х

=

 

 

 

х =

 

 

 

 

 

,

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

+ у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

х

2

+ у

2

 

 

у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у

=

 

 

у =

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

+ у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Значение полного дифференциала в точке М0(2, 1) равно:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

 

 

 

 

 

у

 

 

 

 

 

= 2(0,03) + 1 0,07 = 0,002

 

dz

 

х

=2

=

 

 

х

+

 

 

у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

х2 + у2

 

 

 

 

х2 + у2

 

 

х0

=2

 

 

5

 

 

 

5

 

 

 

 

 

 

 

у0

=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у

 

=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х = -0,03

 

 

 

 

 

 

 

 

 

 

 

х

= -0,03

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у = 0,07

 

 

 

 

 

 

 

 

 

 

 

у = 0,07

 

 

 

 

 

 

 

 

 

 

 

По формуле (60) получаем приблизительное значение величины а:

аf( 4,3 ) + dz М0 = 5 + 0,002 ≈ 5,002.

13.9Производная неявной функции

Напомним, что функция одной переменной называется неявной, если она задается уравнением:

F(х, у) = 0. (13.7)

Предположим, что функция F(х, у) и её частные производные

F , F непрерывны в некоторой области Д, содержащей точку М(х, у).

х у

Дадим переменной х приращение ∆х. Функция у получит приращение ∆у. В силу равенства F(х, у) = 0 будем иметь: F(х + ∆х; у + ∆у) = 0.

Следовательно, F(х + ∆х; у + ∆у) - F(х, у) = 0.

Левую часть равенства, являющуюся полным приращением функции двух переменных, по формуле (13.1) можно переписать так:

 

 

 

 

 

99

 

 

 

 

 

 

 

 

 

 

 

 

F(х + ∆х; у + ∆у) - F(х, у) =

F

х +

F

 

у + α х + β у,

 

 

х

у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α → 0 и β → 0

при ∆х → 0, ∆у → 0.

 

 

 

 

 

 

 

 

 

Так как левая часть полученного равенства равна 0, то

 

 

 

 

 

 

F

х +

F

у + α х + β у = 0.

 

 

 

 

 

 

 

 

 

 

 

х

у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у

 

у

 

 

F

+ α

 

Разделим это равенство на ∆х и выразим из него

:

= −

х

.

х

х

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ β

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у

Устремим ∆х → 0, тогда α → 0

и β → 0. Полагая, что

F

0, в пре-

у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

деле получим:

 

 

ух = −

х

 

 

 

 

 

 

 

 

(13.8)

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у

Пример 35. Найти ddух, если х3 + у3 = 3х + 3у. Решение. В данном случае F(х, у) = х3 + у3 - 3х - 3у.

F = 3 + у3 - 3х - 3у)'х = 2 – 3 = 3(х2 – 1).

х

F = 3 + у3 - 3х - 3у)'у = 2 – 3 = 3(у2 – 1).

у

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

ух =

 

х2

1

.

 

 

 

 

 

1

у2

 

 

 

13.10Частные производные и полные дифференциалы высших порядков

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

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

Частные производные второго порядка определяются и обозначаются следующим образом:

z

=

2 z

= z′′хх

- вторая частная производная по переменной х;

 

 

 

х2

 

х

х

 

 

 

100

z

= ∂х у

z =

у х

z

= ∂у у

2 z

 

= z′′ух

- вторая частная смешанная производная по у и по х;

ух

 

 

 

 

2 z

 

= z′′ху

- вторая частная смешанная производная по х и по у;

ху

 

 

 

 

2 z

 

= z′′уу

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

у2

 

 

 

 

Аналогично определяются частные производные более высоких по-

 

3 z

 

 

2 z

 

4 z

 

 

 

3 z

 

 

рядков. Так

 

 

 

=

 

 

 

 

 

,

 

 

 

=

 

 

 

 

 

и т.д.

 

2

 

 

 

2

 

2

 

 

 

2

 

х

у

 

 

 

х

 

 

ух

х

 

 

ух

 

 

 

 

 

у

 

 

 

 

х

 

 

 

Определение. Частная производная второго и более высокого порядка, взятая по различным переменным, называется смешанной частной производной.

Имеет место следующая теорема.

Теорема 11. Если смешанные частные производные одного порядка, отличающиеся лишь порядком дифференцирования непрерывны, то они равны между собой.

В частности для функции z = f(х, у) имеем: 2 z = 2 z .

ху ух

Определение. Полным дифференциалом второго порядка называется полный дифференциал от полного дифференциала: d 2 z = d (dz) .

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

от полного дифференциала (n – 1)-го порядка:

d n z = d(d n1z) .

 

 

Пользуясь формулой (13.4), получим

 

z dx +

z

 

=

d 2 z = d

 

 

х

у

 

 

 

 

 

 

 

 

z

 

z

 

 

z

 

z

 

2 z

dx 2 +

2 z

 

=

 

 

 

dx +

 

dу dx +

 

 

 

dx +

 

dу dу =

 

 

dydx +

 

 

 

 

 

 

 

 

 

 

 

х

 

у

 

 

 

х

 

у

 

х2

 

ух

 

 

х

 

 

у

 

 

 

 

+

2 z

dxdy + 2 z 2

= 2 z dx2

+ 2

2 z

 

dxdy +

2 z 2 .

 

ху

ху

 

 

у2

х2

 

 

 

у2

 

 

 

Отсюда,

d 2 z = 2 z 2

+ 2

2 z

dxdy + 2 z 2 .

 

 

 

 

 

 

 

 

 

 

 

х2

 

ху

 

у2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

Символически это можно записать так: d 2 z =

dx +

z .

х

у

 

 

 

 

 

 

 

 

 

 

 

 

 

Методом математической индукции можно показать, что

101

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

d n z =

 

dx +

 

 

dу z .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

 

 

у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 36. Найти d 2 z , если z = x3 · sin2y.

 

 

 

 

Решение. Найдем частные производные первого порядка

 

 

 

z

=

(х3 sin 2

 

у)х = 3х2 sin 2 у.

 

 

 

 

 

 

х

 

 

 

 

 

 

 

 

 

(х3 sin 2

 

у)у = х3 2 sin у cos y = x3 sin 2 y .

 

 

 

 

z

 

=

 

 

 

 

 

у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдем частные производные второго порядка:

 

 

2 z

= (3х2

sin 2 у)

х = 6х sin 2 у;

2 z

 

= (х3

sin 2 у)

х = 3х2

sin 2 у;

х2

ух

 

 

2 z

= (х3 sin 2 у)у = 2х3 соs2 у.

 

 

 

 

 

 

у2

 

 

 

 

 

 

 

 

 

 

 

 

Окончательно имеем:

d 2 z = 6хsin 2 уdx2 + 6х2 sin 2 ydxdy + 2x3сos2 ydy 2 .

13.11 Скалярное поле. Производная по направлению. Градиент

Предположим, что в каждой точке Р некоторой области D нам задано значение скалярной физической величины u, т.е. такой величины, которая полностью характеризуется своим числовым значением. Например, это может быть температура точек неравномерно нагретого тела, плотность распределения электрических зарядов в изолированном наэлектризованном теле, потенциал электрического поля и т.д. При этом u называется скалярной функцией точки; записывается это так: u =u (P ).

Область D, в которой определена функция u (P ), может совпадать со

всем пространством, а может являться некоторой его частью. Определение. Если в области D задана скалярная функция точки u (P ), то

говорят, что в этой области задано скалярное поле.

Мы будем считать, что скалярное поле стационарное, т.е. что величина u (P ) не зависит от времени t. Если величина u будет зависеть не только

от точи Р, но и от времени t, то будем иметь нестационарное поле. Если физическая величина векторная, то ей будет соответствовать

векторное поле, например силовое поле, электрическое поле напряжённости, магнитное поле и др.

Если скалярное поле отнесено к системе координат Охуz, то задание точки Р равносильно заданию её координат х, у, z; и тогда функцию u (P )

можно записать в обычном виде функции трёх переменных: u (х,у,z ).

102

Определение. Поверхностью уровня скалярного поля называется геометрическое место точек, в которых функция u принимает постоянное значение, т.е. u (х,у,z )=С.

В курсе физики при рассмотрении поля потенциала поверхности уровня называют обычно эквипотенциальными поверхностями (т.е. поверхностями равного потенциала).

Если в частном случае скалярное поле плоское, т.е. мы изучаем распределение значений физической величины в какой-то плоской области, то функция u зависит от двух переменных, например, х и у. Линиями уровня этого поля будут линии уровня функции u ( х,у), т.е. u (х,у)=С.

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

Производная по направлению.

Важной характеристикой скалярного поля является скорость изменения поля в заданном направлении.

Пусть задано скалярное поле, т.е. задана

 

функция u (х,у,z ). Возьмём точку Р(х,у,z )

 

и какой-нибудь луч λ , из неё выходящий.

 

Направление этого луча зададим углами

 

α ,β ,γ , которые образуются с направлениями

 

осей Ох, Оу, Оz (рисунок 69).

Рисунок 69

Если еλ - единичный вектор, направленный по лучу λ , то его проекциями будут направляющие косинусы

еλ ={cosα ,cosβ ,cosγ }.

Пусть точка P1 (х1 1 ,z1 ) лежит на луче λ ; расстояние РР1 обозна-

чим через ρ . Проекции вектора РР1 на оси координат будут, с одной стороны, равны ρcosα , ρcosβ , ρcosγ , а с другой стороны, - разностями х1 х,

у1 у и z1 z . Следовательно,

х1 = хcosα , y1 = ycosβ , z1 = zcosγ .

Рассмотрим теперь приращение функции u при переходе из точки Р в точку Р1 :

u (P1 )u (P )=u ( хcosα , ycosβ , zcosγ )u (x ,y ,z ).

Если точка Р1 будет изменять своё положение на луче λ , то в выражении для разности u (P1 )u (P ) будет меняться только величина ρ . Составим

отношение u (P1 )u (P ) и перейдём к пределу при ρ→0 , предполагая, что

ρ

этот предел существует.

 

 

 

 

 

103

 

 

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

 

 

 

lim

u (P1

)u (P )

= lim

u ( хcosα , ycosβ , zcosγ )u ( x ,y ,z )

(13.9)

 

ρ

 

ρ

 

ρ→0

 

ρ→0

 

называется производной от функции u (х,у,z ) по направлению λ в точке Р. Этот предел будем обозначать символом λu или u ( x ,y ,z ). Величина его зависит от выбранной точки Р(х,у,z ) и от направления

луча

λ

, т.е. от α ,β и γ .

 

 

Если точка Р фиксирована, то величина производной u

будет зави-

 

 

λ

 

сеть только от направления луча λ .

Из определения производной по направлению следует, что если направление λ совпадает с положительным направлением оси Ох, т.е. α = 0 ,

β= γ =

π , то предел (13.9) будет просто равен частной производной от функ-

 

2

 

 

 

 

 

ции u (х,у,z ) по х:

u (х, y , z )u (x ,y ,z )

 

 

 

 

lim

=

u

.

 

ρ

 

x

 

ρ→0

 

 

 

Аналогичную картину мы получим, если направление λ будет совпа-

дать с направлениями осей Оу и Оz.

 

 

 

 

Подобно тому как частные производные u х

,u у и u z характеризуют

скорость изменения функции u в направлении осей координат, так и производная по направлению u λбудет являться скоростью изменения функции

u ( х,у,z ) в точке Р по направлению луча λ . Абсолютная величина производной u λпо направлению λ определяет величину скорости, а знак произ-

водной – характер изменения функции u (возрастание или убывание). Вычисление производной по направлению производится при помощи

следующей теоремы.

Теорема. Если функция u (х,у,z ) дифференцируема, то её производная u λпо любому направлению λ существует и равна

u

=

u cosα+

u cosβ+

u cosγ ,

(13.10)

∂λ

 

x

y

z

 

где cosα , cosβ и cosγ - направляющие косинусы луча λ .

Доказательство. Так как функция u (х,у,z ) дифференцируема, то её полное приращение можно записать в виде

u =u (х+ x ,у+ y ,z+ z )u ( x ,y ,z )=

u

x+

u

y+

u

z,

x

y

z

 

 

 

 

104

где ε - бесконечно малая величина более высокого порядка, чем

ρ = x 2 + y 2 +

z 2 . Полагая

xcosα ,

 

 

y cosβ ,

z cosγ , пред-

ставим разность u (P1 )u (P ) в виде

 

 

 

 

u

 

 

 

 

 

 

 

 

 

u (P )u (P )=

u

ρcosα+

u

ρcosβ+

ρcosγ +ε ,

 

 

 

 

 

 

 

 

 

1

 

 

x

 

y

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(P )

 

 

 

 

причём

ε

0 при ρ→0 . Отношение

 

u (P1 )u

будет равно

ρ

 

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u (P1 )u (P )

= u cosα+ u cosβ+ u cosγ +

ε

.

 

 

 

 

 

 

 

 

 

 

 

ρ

 

 

x

y

 

z

 

u

 

ρ

 

Так как значения частных производных

 

u

,

 

и

u

в точке Р, а так-

 

 

x

y

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

же α ,β и γ от ρ не зависят, то, переходя к пределу при ρ→0, получим

u

= lim

u (P1 )u (P )

=

u cosα+

u cosβ+

u cosγ , что и требовалось

∂λ

ρ

ρ→0

 

x

y

z

доказать.

Из этой формулы непосредственно следует сделанное выше замечание, что если направление λ совпадает с положительным направлением одной из осей координат, то производная по этому направлению равна соот-

ветствующей частной производной, например, если α =0 , β = γ = π2 , то

u = u . ∂λ ∂x

Из формулы (13.10) видно, что производная по направлению λ 0 , противоположному направлению λ , равна производной по направлению λ ,

взятой с обратным знаком. Действительно, при перемене направления углы

α ,β и γ изменятся на π

и

 

 

 

u

=

u cos (α+π )+

u cos (β+π )+

u cos (γ+π )= −

u .

 

∂λ ′

 

 

x

y

z

∂λ

Это означает, что при перемене направления на противоположное абсолютная величина скорости изменения функции u не меняется, а изменяется только характер её изменения; если, например, в направлении λ функция

возрастает, то в направлении λ 0 она убывает, и наоборот.

Если поле плоское, то направление луча λ вполне определяется углом α его наклона к оси абсцисс. Формулу для производной по направлению в

случае плоского поля можно получить из общей формулы, положив γ =

π

и

 

π

 

 

 

 

2

 

β =

−α . Тогда

 

 

 

 

 

 

2

u

 

u cosα+

u sinα .

 

 

 

 

=

 

 

 

 

∂λ

 

x

y

 

 

 

 

 

 

 

105

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если α = 0, то

u

=

u

, а если α = π , то

 

u

 

=

u

.

 

 

 

∂λ

 

 

∂λ

 

 

 

 

x

2

 

 

 

 

y

 

 

 

 

 

Градиент

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим снова формулу для производной по направлению

u = u cosα+

u cosβ+

u cosγ .

 

 

 

 

 

 

 

∂λ

x

y

z

 

 

 

 

 

 

 

 

 

 

 

Вторые множители в каждом слагаемом являются, как мы уже отме-

чали, проекциями единичного вектора

 

 

 

 

 

 

 

е

λ , направленного по лучу

λ

:

еλ ={cosα ,cosβ ,cosγ }.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Возьмём теперь вектор, проекциями которого на оси координат будут

служить значения частных производных

 

u

,

u

 

и

u

в выбранной точке

 

 

y

 

 

 

 

 

 

 

 

 

 

x

 

 

z

P ( х,у,z ). Назовём этот вектор градиентом функции u ( х,у,z ) и будем обозначать его символами grad u или u .

Определение. Градиентом функции u (х,у,z ) называется вектор, проекция-

ми которого служат значения частных производных этой функции, т.е. grad u = ux i+ uy j+ uz k .

Подчеркнем, что проекции градиента зависят от выбора точки

P ( х,у,z ) и изменяются с изменением координат этой точки. Таким обра-

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

u ( х,у,z ), соответствует определенный вектор - градиент этой функции. Отметим, что градиент линейной функции u = ax+by+cz+d есть постоянный вектор: grad u = ai+bj+ck .

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

u

= grad u

 

λ .

(13.11)

e

λ

 

 

 

 

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

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

Производная функции по данному направлению равна проекции градиента функции на направление дифференцирования, т. е.

u

=

 

grad u

 

cosϕ ,

Рисунок 70

 

 

∂λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

ϕ - угол между

вектором grad u и лучом

λ

(рисунок 70).

106

Отсюда сразу следует, что производная по направлению достигает наибольшего значения, когда cosϕ=1, т. е. при ϕ= 0 . Это наибольшее

значение равно grad u .

Итак, grad u есть наибольшее возможное значение производной u λв данной точке Р, а направление grad u совпадает с направлением луча, вы-

ходящего из точки Р, вдоль которого функция меняется быстрее всего, т. е. направление градиента есть направление наискорейшего возрастания функции. Ясно, что в противоположном направлении функция u будет быстрее всего убывать.

Докажем теперь теорему, устанавливающую связь между направлением градиента функции и поверхностями уровня скалярного поля.

Теорема. Направление градиента функции u ( х,у,z ) в каждой точке

совпадает с направлением нормали к

поверхности уровня скалярно-

го поля, проходящей через эту точку

(рисунок 71).

Доказательство. Выберем произвольную точку P0 (х0 0 ,z 0 ). Уравнение поверхности уровня, проходящей через точку P0 =u 0 ,

где u 0 =u (х0 0 ,z 0 ).

Напомним, что прямая, перпендикулярная к касательной плоскости в точке касания, называется нормалью к поверхности в этой точке.

На основании известных из аналитической геометрии правил уравнение нормали к поверхности U (x ,y ,z )=0 в её точке

P0 (x 0 ,y 0 ,0 ) имеет вид:

 

 

 

 

хх0

 

=

 

 

yy 0

 

=

 

zz 0

 

.

Рисунок 71

 

 

 

u

(P )0

 

 

u

 

(P )

0

 

u

(P )0

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

Отсюда и следует, что направляющий вектор нормали, имеющий про-

екции

 

 

 

u

(Р

0

),

 

 

u

(Р

0

),

u

(Р

0

), является градиентом функции u ( х,у,z )

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

у

 

 

 

 

z

 

 

 

в точке P0 , что и требовалось доказать.

Таким образом, градиент в каждой точке перпендикулярен касательной плоскости к поверхности уровня, проходящей через данную точку, т. е. его проекция на эту плоскость равна нулю. Следовательно:

Производная по любому направлению, касательному к поверхности уровня, проходящей через данную точку, равна нулю.

Укажем теперь некоторые свойства градиента функции, часто облегчающие его вычисление.

107

1)grad (u1 +u 2 )= grad u1 + grad u 2 .

2)grad Cu1 =C gradu1, где С – постоянная.

3)grad u1 u 2 =u 2 gradu1 +u1 gradu 2 .

4)grad f (u )= f (u )grad u .

Перечисленные свойства градиента показывают, что правила его отыскания совпадают с правилами отыскания производной функции.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x

 

 

 

 

 

)= x 2 +

x 2

Пример 37.

Пусть в

R 2 задана функция f

1

; x

2

2

. Найти ли-

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

нии уровня

f (x1 ; x 2 )=С и изобразить их в плоскости Ох1 х2 .

Решение. При C < 0 её линия уровня

f ( x )=C представляет собой пустое

множество, так как при всех x1 , x 2

верно неравенство

 

 

 

 

 

 

 

 

 

 

 

 

 

(x

 

 

 

)= x 2 +

 

x 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

1

; x

2

 

2

 

0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

При C =0 линия уровня состоит из одной точки О(0; 0) , так как равенство

 

 

 

 

 

 

 

 

 

 

f (x1 ; x 2

)= x 2 +

x 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

= 0

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

возможно лишь при x1 =0 , x 2 =0.

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При C >0 уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x

1

; x

2

)= x 2

+

 

x 22

=C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

О(0; 0) и полуосями a =

C и b= 2

 

C ,

 

 

 

 

 

 

 

 

 

поскольку при C >0 уравнение можно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 2

 

 

 

x 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

записать в виде

 

1

+

 

 

2

 

=1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При увеличении С полуоси а и b увеличиваются

 

Рисунок 72

пропорционально друг другу, так что эксцентриситеты всех эллипсов совпадают.

Пример 38. Найти производную функции f ( x ; y )=3 x 2 +2 xy в точке M 0 (1; 2 ) по направлению , составляющему угол 30° с направлением

оси Ох.

Решение. Из условия следует, что направляющие косинусы оси равны cosα = cos30 0 = 23 ; cosβ = cos60 0 = 12 .

Частные производные функции f равняются f x′ =6 x+ 2 y ; f y′ = 2 x .

Их значение в точке M 0 (1; 2 ) таково:

f x(M 0 )= 6+2 (2 )= 2; f y(M 0 )= 2 1= 2 .

 

 

 

 

 

 

 

 

108

 

 

Для нахождения производной по направлению применяем формулу

 

 

 

f

(М 0

)=

f

(М 0 )cosα+

f

(М 0 )cosβ .

 

 

 

 

х

y

 

 

 

 

 

 

 

Получаем:

f

(М 0

)= 2

3

+2 1

= 3+1.

 

 

 

2

 

 

 

 

 

2

 

 

 

Пример 39. Найти производную функции f ( x ; y ;z )= xy 2 z+3 x 2 yz 3

в точке M 0 (1;1; 2 ) по направлению , заданному вектором a (=1; 2; 2 ).

Решение. Находим единичный направляющий вектор τ оси : τ =

 

1

 

a .

 

a

 

 

 

 

 

Поскольку a = 12 +(2 )2 +2 2 =3 , получаем, что

τ= 13(1; 2; 2 )= 13; 23; 23 .

Найдём теперь градиент: поскольку

 

 

 

f

= у 2 z+6 xyz 3

;

f

= 2 xyz+3 x 2 z 3 ;

f

 

= 2

+9 x 2 yz 2

 

 

 

 

х

y

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

и

 

f

(M 0 )= 20;

 

f

(M 0 )=35,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

z

 

 

 

 

 

 

 

 

 

то (grad f

)(M 0 )=(46; 20; 35 ).

 

 

 

 

 

 

 

 

Теперь находим производную по направлению:

 

 

 

 

 

f

(M 0 )=( grad f )(M 0

)τ = −46 1

+20

2

+35 2

= −16 .

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

3

 

 

3

3

 

 

Пример 40. Найти уравнение касательной плоскости к поверхности S, за-

данной уравнением

х2

у 2 z = 0

, в точке M 0 (2; 1;1).

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

Решение.

Уравнение

касательной

 

 

плоскости

к

поверхности

 

S ={ f (x ; y ;z )= 0} в точке M 0 (x 0

; y 0 ;z 0 ) S

имеет вид

 

 

f x(x 0 ;y 0 ;z 0 )(xx 0 )+ f y(x 0 ; y 0 ;z 0 )( yy 0 )+ f z(x 0 ; y 0 ;z 0 )(zz 0 )=0

Находим частные производные и их значения в точке M 0 : f x′ = x ; f y′ = −2 y ; f z′ = −1;

f x(x 0 ; y 0 ;z 0 )= 2; f y(x 0 ; y 0 ;z 0 )= 2; f z(x 0 ; y 0 ;z 0 )= −1.

Поэтому искомое уравнение касательной плоскости имеет вид

2 (x2 )+2 ( y+1)(z1)= 0 , или 2 x+2 yz1=0 .

109

Пример 41. Найти уравнения касательной и нормали, проведённых к ли-

нии уровня C =3 функции

f ( x ; y )= 2 x 2 y 3 + xy 4 в точке M 0

(1;1).

Решение. Линия уровня задаётся уравнением 2 x 2 y 3 + xy 4 =3.

Касательная к этой линии уровня имеет вид

 

 

f

(М 0

)(xx 0 )+

f

(М 0 )( yy 0 )= 0 ,

 

 

 

 

 

 

x

 

y

 

а нормаль (напомним, что нормаль - это прямая, проходящая через точку касания перпендикулярно касательно) можно задать уравнением вида

 

 

xx 0

 

=

 

yy 0

 

 

 

 

 

 

 

 

.

 

 

f

(М 0

)

 

f

(М 0 )

 

 

 

 

 

 

x

 

 

 

y

Вычисляем значения частных производных:

f x′ = 4 3 + y 4 ; f y′ = 6 x 2 у 2 +4 xy 3 ; f x(M 0 )= 4+1=5; f y(M 0 )=6+4=10 ,

откуда получаем уравнение касательной:

5(x1)+10 ( y1)=0, или x+2 y3= 0 , и нормали:

x1= y1, или 2 xy1=0. 5 10

Пример 42. Найти уравнения касательной плоскости и нормали, проведён-

ных к поверхности уровня функции f (x; y; z) = x2 + y2 + 4z2 , проходящей через точку M 0 (2;1;1).

Решение. Найдём уравнение поверхности уровня: поскольку

f (M 0 ) = 4 + 1 + 4 = 9 , поверхность уровня задаётся уравнением

x2 + y

2 + 4z2 = 9, или

x2

+

 

y2

+

z2

= 1.

 

 

 

 

 

 

 

 

9

 

 

 

 

 

 

 

 

 

 

 

 

 

9

 

 

 

4

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это уравнение задаёт эллипсоид вращения (вокруг оси Оz) с центром в точ-

ке О(0;0;0) и полуосями а = 3, b = 3 и c

= 3 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

f

в точке M 0 :

Найдём значения частных производных функции

f

(M

 

) = 2x

 

 

= −4 ; f (M

 

 

) = 2 y

 

 

 

= 2 ; f

(M

 

) = 8z

 

= −8 .

 

 

 

 

 

 

 

 

 

 

x

 

0

 

 

M 0

 

y

0

 

 

 

 

M 0

z

 

0

 

 

M 0

 

 

 

 

 

 

 

 

 

 

 

Записываем уравнение касательной плоскости в виде

 

 

 

 

f x(M 0 )(x x0 ) + f y(M 0 )(y y0 ) + f z(M 0 )(z z0 ) = 0 ,

 

то есть 4(x + 2)+ 2(y 1)8(z + 1) = 0 , или 2x y + 4z + 9

= 0 .

Уравнения нормальной прямой записываем в виде

 

 

 

 

 

x x0

y y0

z z0

x + 2

y 1

 

z + 1

 

 

 

=

 

=

 

, то есть

 

=

 

 

=

 

.

 

f x(M 0 )

f y(M 0 )

f z(M 0 )

4

2

 

1

 

 

 

 

 

110

 

 

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

 

Пусть функция z = f(х, у) определена в некоторой области D и точка

М0(х0, у0) D (внутренняя точка области).

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

 

Точка М0(х0, у0) называется точкой максимума (минимума)

 

 

 

функции z = f(х, у), если в достаточно малой δ-окрестности

 

 

 

точки М0

для каждой точки М(х, у) отличной от точки М0(х0,

 

 

 

у0), выполняется неравенство:

 

 

 

f(x, y) < f(x0, y0)

(f(x, y) > f(x0, y0)).

 

z

Р0

 

На рисунке 73 М0 – точка макси-

 

 

Р

 

мума, а точка М1 – точка мини-

 

 

 

Р Р1

мума функции z = f(х, у).

 

 

 

 

 

 

 

 

Значение функции в точке макси-

 

 

 

 

 

 

 

 

 

у

мума (минимума) называется максиму-

 

 

 

 

мом (минимумом) функ-

х

М0

М

ММ1

ции. Максимум и минимум функ-

ции называется её экстремумом.

 

 

 

Рисунок 73 Замечание. Согласно определению, точка экстремума функции является внутренней точкой области определения функции. Максимум и минимум имеют локальный (местный) характер: значения функции в точке М0(х0, у0) сравниваются с её значениями в точках достаточно близких к М0(х0, у0). В области D функция может иметь несколько экстремумов или не иметь ни одного.

13.13Необходимые и достаточные условия существования экстремума функции z = f(х, у)

Теорема 12 (необходимые условия существования экстремума).

Если функция z = f(х, у) дифференцируема в точке М0(х0, у0) и имеет в ней экстремум, то частные производные функции в этой точке рав-

ны нулю:

z

 

М

0

= 0,

z

 

М

0

= 0.

 

 

х

у

 

 

 

 

 

 

 

Доказательство. Зафиксируем одну из переменных, например, положим у = у0. Тогда функция f(x, y0) зависит от одной переменной и имеет экстремум в точке х = х0. Согласно необходимому условию экстремума функции одной переменной в этой точке производная функции равна нулю, т.е.

f 'x (х0, y0) = 0.

111

Рассуждая аналогично, можно показать, что f 'у (х0, y0) = 0. Следовательно, если функция z = f(х, у) имеет в точке М0(х0, у0) экстремум, то

 

z

 

 

= 0

 

 

 

х

 

М0

 

.

 

z

 

 

 

 

 

 

 

М0

= 0

 

у

 

 

 

 

 

 

 

 

 

Следствие 1. Функция z = f(х, у) может иметь экстремум только в тех точках, где частные производные либо равны нулю, либо не существуют.

Определение. Точки, в которых частные производные функции z = f(х, у) либо обращаются в нуль, либо не существуют, называются критическими точками функции.

Итак, функция может иметь экстремум только в критических точках. Однако не всякая критическая точка является точкой экстремума. Например, функция z = x2 – y2 имеет частные производные и точка О(0; 0) является для неё критической. Но эта точка не является точкой экстремума, т.к. в

достаточно малой окрестности точки О(0; 0) найдутся точки, для которых z > 0 (z(x, 0) = x2 > 0) и z < 0 (z(0, y) = -y2 < 0).

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

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

Теорема 13 (достаточные условия существования экстремума).

Пусть функция z = f(х, у) в критической точке М0(х0, у0) и некоторой её окрестности имеет непрерывные частные производные второго порядка. Положим

А=

2 z

 

М0 ,

В =

2 z

 

 

М0 ,

С =

2 z

 

М0 ,

=

 

А В

 

= АС В

2

 

 

 

 

 

 

 

 

 

 

 

х2

 

 

 

 

у2

 

 

В С

 

 

 

ху

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда: 1) если ∆ > 0, то функция z = f(х, у) имеет экстремум в точке М0(х0, у0): максимум, если А< 0 и минимум, если А > 0;

2)если ∆<0, то функция z=f(х, у) в точке М0(х0,, у0) экстремума не имеет;

3)в случае, если ∆ = 0, то в точке М0(х0, у0) экстремум может быть, может не быть. Необходимы дополнительные исследования.

Пример 43. Исследовать на экстремум функцию z = 3xy – x2y – xy2 . Решение. Областью определения функции является множество точек на

плоскости хОу, т.е. (х, у) R2. Найдем частные производные функции:

112

хz = (3ху х2 у ху2 )х = 3у 2ху у2 ; уz = (3ху х2 у ху2 )у = 3х х2 2ху.

Точки, в которых частные производные не существуют, отсутствуют. Следовательно, критические точки найдем, решая систему уравнений:

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

= 0 у(3 2х у) = 0

у = 0

или

 

у = 0

3у 2ху у

 

 

 

2

2ху = 0

х(3

х 2 у) = 0

х = 0

 

3 х 2 у

3х х

 

 

 

 

 

 

 

 

 

х = 0

 

или

3 2х у = 0

.

 

 

 

 

 

2х у = 0

 

2 у = 0

 

 

 

 

3

 

3 х

 

 

 

 

 

Отсюда получим точки: М1(0, 0), М2(3, 0), М3(0, 3), М4(1, 1).

Найдем частные производные второго порядка данной функции:

2 z

= (3у 2ху у2 )х = −2 у,

2 z

= (3х х2 2ху)х = 3

х2

ху

= (3х х2 2ху)у = −2х.

 

2 z

 

 

у2

 

 

 

= 0 или

2х 2 у,

Вточке М1(0, 0) имеем: А1 = 0, В1 = 3, С1 = 0, ∆1 = АС – В2 = –9 < 0.

Т.к. ∆1 < 0, то в точке М1(0, 0) функция экстремума не имеет.

Вточке М2(3, 0) имеем: А2 = 0, В2 = –3, С2 = –6, ∆2 = –9 < 0. Также экстремума нет.

Вточке М3(0, 3) имеем: А3 = –6, В2 = –3, С3 = 0, ∆3 = –9 < 0. Функция экстремума не имеет.

Вточке М4(1, 1) имеем: А4 = –2, В4 = –1, С4 = –2, ∆4 = 3, т.е. ∆4 > 0.

Т.к. А4<0, то в точке М4 функция имеет максимум: zmax=z(1; 1)=3 –1 – 1 =1.

13.14 Метод наименьших квадратов

На практике часто требуется по результатам наблюдений двух величин установить приближенную аналитическую зависимость между ними.

 

х

х1

х2

хn

 

у

у1

y2

yn

 

Один из способов решения этой задачи состоит в следующем: вид ис-

комой функции у = f(х, а, b, с, …), где , а, b, с, … - параметры, устанавливается или из теоретических соображений, или на основании характера расположенных на координатной плоскости точек, соответствующих экспериментальным данным (хi, уi). Неизвестные параметры следует подобрать так,

чтобы разность между экспериментальными значениями функции у2э и расчетными значениями у2р была минимальной.

Таким образом, требуется подобрать а, b, с … так, чтобы ∆уj =

уiэ уiр было минимальным. Удобнее минимизировать функцию, равную сумме квадратов указанных разностей, которые принимают за меру откло-

113

нения экспериментальных и расчетных значений, т.е.

n

 

n

2

уi2

(уiэ уiр ) .

S =

=

i=1

 

i=1

 

Sа

Для этого достаточно составить систему уравнений: S

bSс

=0

=0 и, решив её,

=0

найти значения параметров.

Пусть, например, экспериментальные точки расположены на координатной плоскости так, как изображено на рисунке 74.

Естественно предположить, что искомую функцию можно искать в виде линейной функции у = ах + b, где а и b – параметры. В этом слу-

чае имеем: S(a, b) = n (уi ахi b)2

min .

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

Найдем частные производные функции S(a, b) по неизвестным а и b:

S

= n 2(уi ахi b) (хi ),

S

= n 2(уi

 

ахi

b) (1) и после упро-

а

i=1

b

 

i=1

 

 

 

 

 

 

 

 

щений получим систему уравнений:

 

 

 

 

 

 

 

 

 

 

n

2

n

х

 

n

х

 

у

 

 

 

 

а х

+ b

 

=

 

 

 

 

i=1

i

i=1

 

i

i=1

 

i

 

i

(13.12)

 

 

 

n

 

 

 

n

 

 

 

.

 

 

а хi + b

n = уi

 

 

 

 

 

i=1

 

 

i=1

 

 

 

 

Решая полученную систему уравнений, например, по формулам Крамера, найдем значения параметров а и b.

Если же экспериментальные точки расположены так, как на рисунке 75, то искомую функцию естественно искать в виде у = ах2 + bx + c.

n ( 2 )2

В этом случае S(a, b, c) = уi ахi i с .

i=1

Система уравнений относительно параметров а, b, с имеет вид:

 

n

 

n

 

n

n

 

а хi4

+ b хi3

+ с хi2

= хi2 уi

 

i=1

 

i=1

 

i=1

i=1

 

 

n

 

n

 

n

n

уi .

а хi3

+ b хi2

+

с хi

= хi

 

i=1

 

i=1

 

i=1

i=1

 

 

n

хi2

n

 

 

n

 

а

+ b

хi + с

n = уi

 

 

i=1

 

i=1

 

 

i=1

 

114

у

у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 х1 х2 х3 х4

х

0 х1 х2 х

3 х4 х5 х6

х7

х

Рисунок 74

 

 

Рисунок 75

Пример 44. Исследуя истечение жидкости через щель, получены следующие данные зависимости коэффициента расхода μ от напора Н:

μ

0,448

0,432

0,421

0,417

0,414

0,412

Н

0,164

0,328

0,636

0,984

1,312

1,640

Найти формулу зависимости в виде μ = На + b .

Решение. Перепишем искомую зависимость в виде μ = a Н1 + b . Для нахождения неизвестных коэффициентов составим следующую таблицу:

Нi

 

1

 

μi

 

1

 

μi

1

 

 

 

Нi2

 

 

Нi

 

H i

 

 

 

 

 

 

 

0,164

6,098

0,448

37,186

2,732

 

0,328

3,049

0,432

9,296

1,317

 

0,636

1,572

0,421

2,471

0,662

 

0,984

1,016

0,417

1,032

0,424

 

1,312

0,762

0,414

0,581

0,315

 

1,640

0,610

0,412

0,372

0,251

 

Σ

13,107

2,544

50,938

5,701

 

 

 

6

1

6

1

 

6

1

 

a

+ b

=

μi

 

 

 

Hi

Составим систему уравнений вида (63):

 

i=1Hi2

i=1Hi

i=1

 

6

1

 

6

 

 

 

 

а

 

+ 6b = μi

 

 

 

 

 

 

 

 

i=1Hi

 

i=1

 

 

 

 

 

 

 

50 ,938 а+13 ,107 b=5 ,701. По формулам Крамера имеем:

13 ,107 а+6 b=2 ,544

=50,938 13,107 = 305,628 171,793 = 133,835

13,107 6

 

 

 

 

 

 

 

 

115

а =

 

 

5,701 13,107

 

 

 

= 34,206 33,344 = 0,862

 

 

 

 

 

 

 

 

2,544

6

 

 

 

 

 

 

b =

 

 

50,938

5,701

 

= 129,586 74,723 = 54,863

 

 

 

 

 

 

13,107

2,544

 

 

 

 

а =

 

0,862

= 0,006 , b =

54,863

= 0,410 .

133,835

 

 

 

 

 

133,835

 

Таким образом, искомая формула имеет вид: μ = 0,Н006 + 0,410 .

ЛИТЕРАТУРА

1.Щипачев, В.С. Высшая математика. Полный курс. Учебник для бакалавров / В.С. Щипачев. – М.: Юрайт, 2012.

2.Пискунов, Н.С. Дифференциальное и интегральное исчисление. В 2 т. Т.1: Учеб. пособие для студ. техн. вузов. - М.: Интеграл-Пресс, 2005.

3.Поспелов, А.С. Сборник задач по высшей математике Учебное пособие для бакалавров. Часть 1 / А.С. Поспелов. – М.: Юрайт, 2011.

116

Учебное издание

Кравченко Наталья Ивановна

М А Т Е М А Т И К А

ВВЕДЕНИЕ В МАТЕМАТИЧЕСКИЙ АНАЛИЗ, ДИФФЕРЕНЦИАЛЬНОЕ ИСЧИСЛЕНИЕ ФУНКЦИЙ ОДНОЙ И НЕСКОЛЬКИХ ПЕРЕМЕННЫХ, ПРИЛОЖЕНИЯ ДИФФЕРЕНЦИАЛЬНОГО ИСЧИСЛЕНИЯ

КИССЛЕДОВАНИЮ ФУНКЦИЙ ОДНОЙ

ИНЕСКОЛЬКИХ ПЕРЕМЕННЫХ

Учебное пособие В шести частях Часть 2

Подписано к печати

Тираж

Формат 60×84 1/16

Объем

Заказ

Отдел оперативной полиграфии ФГБОУ ВПО НГМА, 346428, г. Новочеркасск, ул. Пушкинская, 111