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

Конспект лекций Высшая математика (Басканова)

.pdf
Скачиваний:
151
Добавлен:
17.01.2018
Размер:
5.08 Mб
Скачать

429

где z – новая неизвестная функция, а α(x) и β(x) – произвольные

непрерывно дифференцируемые функции от x, причем α(x) 0 в ( a,b). Действительно, так как

y=α′(x)z+α(x)z+β′(x)

то после преобразования получим

α′(x)z+α(x)z+β′(x)+p(x)[α(x) z + β(x)]=q(x)

или

z′+

α′(x) + p(x)α(x) z = q(x) −β′(x) p(x)β(x)

,

 

α(x)

α(x)

 

(4.3.12)

(4.3.13)

т. е. опять линейное уравнение, у которого коэффициент при искомой функции и правая часть непрерывны в (a, b).

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

Покажем, что линейное уравнение всегда интегрируется в квадратурах. Рассмотрим сначала однородное линейное уравнение (4.3.3)

dydx + p(x)y = 0

где функция p(x) непрерывна в интервале (a, b). Переписав это уравнение в виде

dy + p(x)ydx = 0

(4.3.14)

получаем уравнение с разделяющимися переменными. Разделяя переменные, получим:

 

dy

+ p(x)dx = 0

(y = 0?) ,

 

 

y

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

y = Ce

p(x)dx

,

(4.3.15)

 

 

 

 

где C – произвольная постоянная.

Все решения уравнения (4.3.3) содержатся в формуле (4.3.15), так как разделяя переменные, мы могли потерять лишь очевидное нулевое решение y = 0, но и оно содержится в (4.3.15) при C = 0. Из формулы (4.3.15) видно, что всякое решение уравнения (4.3.3) определено во всем интервале (a, b).

Покажем, что функция (4.3.15) является общим решением уравнения (4.3.3) в области (4.3.7), т. е. во всей области задания уравнения (4.3.3).

В самом деле (4.3.15) разрешимо относительно C в области (4.3.7), так что мы имеем:

C = yep(x)dx ,

где функция справа определена в области (4.3.7). Кроме того, по

430

построению функция (4.3.15) является решением уравнения (4.3.3) в интервале (a, b) при всех значениях произвольной постоянной C. А это и означает, что функция (4.3.15) есть общее решение уравнения (4.3.3) в области (4.3.7).

Пример 1. Рассмотрим уравнение

y

+

 

 

 

x

 

 

 

 

y = 0.

 

 

 

 

 

 

 

 

 

1

x2

 

Здесь коэффициент p(x) =

 

 

 

 

x

 

 

 

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

 

 

 

 

 

 

 

 

 

 

1 x

2

 

 

 

 

 

 

 

 

непрерывная в интервале (1,+1). Пользуясь формулой (4.3.15), найдем: y = Ce1x2 .

Это есть общее решение рассматриваемого уравнения в области:

1<x<+1, −∞<y<+∞.

Пример 2. Найти решение уравнения: y+2xy=0.

В этом случае p(x)=2x не имеет точек разрыва. Поэтому всякое решение определено при всех x. Действительно, интегрируя данное уравнение, получим:

y = Cex2 ,

откуда и вытекает наше утверждение.

Решения однородного линейного уравнения обладают следующими двумя характерными для этого уравнения свойствами.

1. Если y1 есть частное решение уравнения (4.3.3), т. е. имеет место тождество

y1′ + p(x)y1 0 (a<x<b),

(4.3.16)

то функция

 

y=Cy1,

(4.3.17)

где C – произвольная постоянная, тоже является решением этого

уравнения.

 

Действительно, полагая в левой части уравнения (4.3.3) y = Cy1 и

принимая во внимание тождество (4.3.16), получим:

 

(Сy1 )′ + p(x)(Сy1 ) = С[y1′ + p(x)y1 ] 0

(a<x<b).

Следовательно, y = Cyl есть решение уравнения (4.3.3).

2. Если y1 – ненулевое частное решение уравнения (4.3.3), то формула (4.3.17), где C – произвольная постоянная, дает общее решение уравнения

(4.3.3) в области (4.3.7).

В самом деле, уравнение (4.3.17) разрешимо в о ластиб ( 4.3.7) относительно C:

431

 

С =

y

(4.3.18)

y

 

 

 

1

 

и, как показано выше, функция (4.3.17) является решением уравнения (4.3.3) при всех значениях C. Следовательно, функция (4.3.17) есть общее решение уравнения (4.3.3) в области (4.3.7).

Таким образом, для построения общего решения однородного линейного уравнения достаточно найти какое-нибудь одно ненулевое частное решение.

Рассмотрим теперь неоднородное уравнение (4.3.1)

dydx + p(x)y = q(x) .

Предположим, что нам известно некоторое решение y1 этого уравнения, т. е. имеем тождество:

y1′ + p(x)y1 q(x)

(a<x<b),

(4.3.19)

Введем новую неизвестную функцию z по формуле

 

y = y1 + z

 

(4.3.20)

Подставляя (4.3.20) в (4.3.1) имеем:

 

 

y1′ + z′ + p(x)y1 + p(x)z = q(x)

(4.3.21)

Отсюда, в силу тождества (4.3.19), получаем:

 

 

z′ + p(x)z = 0

 

(4.3.22)

Мы получили для определения z однородное линейное уравнение, левая часть которого имеет тот же вид, что и левая часть уравнения (4.3.3).

Уравнение (4.3.22) называется однородным линейным уравнением, соответствующим неоднородному линейному уравнению (4.3.1). Общее решение уравнения (4.3.22) имеет вид:

 

z = Ce

p(x)dx

(4.3.23)

 

 

где C – произвольная постоянная. Подставляя найденное значение z в

(4.3.20), получим:

 

 

 

y = y

+ Cep(x)dx

(4.3.24)

1

 

 

Все решения уравнения (4.3.1) содержатся в формуле (4.3.24). Эта формула представляет собою общее решение уравнения (4.3.1) в полосе

a< x <b, −∞<y< +(4.3.25)

т. е. во всей области задания уравнения (4.3.1).

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

Теорема. Если y1 есть частное решение неоднородного линейного уравнения (4.3.1),

y′+ p(x)y = q(x)

то общее решение этого уравнения дается формулой (4.3.20), y=y1+z,

432

где

z = Cep(x)dx

есть общее решение соответствующего однородного линейного уравнения (4.3.22),

z′ + p(x)z = 0 .

Из этой теоремы следует, что знание одного частного решения неоднородного уравнения (4.3.1) дает возможность получить общее решение при помощи одной квадратуры.

Если мы знаем не одно, а два частных решения y1 и y2 неоднородного уравнения (4.3.1), то общее решение можно получить вовсе без квадратур, а именно общим решением будет

y=y1+C(y2y1). (4.3.26)

Действительно, из формулы (4.3.20) имеем: z = y y1; заменяя здесь y на y2, видим, что y2y1 есть частное решение однородного уравнения (4.3.22).

Тогда общее решение уравнения (4.3.22) дается формулой z = С(y 2y1) и, согласно доказанной теореме, формула (4.3.26) дает общее решение уравнения (4.3.1).

Выяснив структуру общего решения уравнения (4.3.1), укажем один общий способ фактического построения общего решения

4.3.2. Метод вариации произвольной постоянной (метод Лагранжа)

Будем искать решение уравнения (4.3.1) в том же виде, что и общее решение (4.3.23) соответствующего однородного уравнения (4.3.20), но будем считать C не постоянной, а некоторой непрерывно дифференцируемой функцией от x, т. е. положим

 

 

 

 

y = C(x)e

p(x)dx

 

 

 

 

(4.3.27)

 

 

 

 

 

 

 

 

 

 

 

и выберем функцию C(x) так, чтобы (4.3.27) удовлетворяло уравнению

(4.3.1). Подставляем (4.3.27) в (4.3.1):

 

 

 

 

 

 

 

 

 

 

p(x)dx

C(x) p(x)e

p(x)dx

+ p(x)C(x)e

p(x)dx

= q(x)

C (x)e

 

 

 

 

 

 

 

откуда:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= q(x)e

p(x)dx

 

 

 

 

 

 

 

 

C (x)

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C(x) = q(x)ep(x)dx dx + C

 

 

 

(4.3.28)

где C – произвольная постоянная. Подставляя это значение C(x) в

формулу (4.3.27), получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = e

p(x)dx

q(x)e

p(x)dx

dx + C

 

 

 

(4.3.29)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

433

Это есть решение уравнения (4.3.1) по построению и притом общее в полосе (4.3.25), так как она имеет структуру (4.3.20).

Таким образом, общее решение неоднородного линейного уравнения (4.3.1) всегда может быть найдено двумя квадратурами.

Замечание. Общее решение неоднородного линейного уравнения (4.3.1) можно найти также следующим методом, принадлежащим Эйлеру.

Умножим обе части уравнения (4.3.1) на функцию

μ(x) = e

 

p(x)dx

.

 

 

 

 

 

 

 

(4.3.30)

 

 

 

 

 

 

 

 

 

 

 

 

Получим уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

′ ∫ p(x)dx

 

+ p(x)ye

p(x)dx

= q(x)e

p(x)dx

,

(4.3.31)

y e

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

ye

p(x)dx

,

 

 

 

 

 

 

(4.3.32)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

так что мы можем переписать это уравнение в виде

 

 

 

 

ye

p(x)dx

= q(x)e

p(x)dx

 

 

 

(4.3.33)

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yep(x)dx

= q(x)ep(x)dx dx + C

 

 

 

(4.3.34)

и, следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = e

p(x)dx

q(x)e

p(x)dx

dx + C

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Мы получили тот же вид общего решения, что и при применении метода Лагранжа.

Функция (4.3.30) называется интегрирующим множителем

линейного уравнения (4.3.1), а изложенный метод Эйлера – методом интегрирующего множителя.

Пример 1. Рассмотрим уравнение

y′ − 2x y = x.

Найдем его решение методом вариации произвольной постоянной.

Соответствующее однородное уравнение

z′ −

2

z =

0 имеет общее решение

x

z=Cx2.

 

 

 

 

 

 

 

Ищем общее решение данного неоднородного уравнения в виде

 

 

y=C(x)x2.

 

 

 

 

Подставляя его в уравнение имеем:

2

 

 

 

 

2

 

 

 

2

 

С (x)x

 

+ C(x)2x x C(x)x

 

= x

или

434

C(x) = 1x ,

откуда

C(x) = ln x +C .

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

y = x2 (ln x + C).

Это и есть общее решение уравнение.

Проинтегрируем это уравнение методом интегрирующего множителя. Имеем:

µ(x) = e

2dx

=

1

.

 

x

 

 

1

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

Умножая обе части уравнения на

 

 

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

x2

 

 

1

 

 

1

 

 

 

 

 

 

 

y =

 

,

 

 

2

 

x

 

x

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

y = x2 (ln

 

x

 

+ C).

 

 

 

 

Пример 2. Найти общее решение уравнения xy+2x2y=1, пользуясь формулой общего решения (4.3.29).

Имеем:

p(x)=2x, q(x)=1/x.

Подставляя в формулу (4.3.29) получим:

y = e

 

 

 

 

C +

1

e

 

 

dx

 

 

2xdx

 

 

 

2xdx

 

 

или

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

2

 

e

x2

 

 

 

y = e

 

 

dx

 

 

 

 

 

x

 

 

 

C +

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.3.3. Уравнение Бернулли

Определение 3. Уравнение вида

y+p(x) y = q(x)ym (4.3.35)

где m любое вещественное число, называется уравнением Бернулли. Будем считать, что m отлично от0 и 1,ибо в этих случаях уравнение Бернулли вырождается в линейное. Относительно функций p(x) и q(x) будем

предполагать, что они непрерывны в интервале (a, b).

Уравнение Бернулли всегда может быть сведено к линейному

435

уравнению. В самом деле, преобразуем сначала правую часть уравнения Бернулли к виду правой части линейного уравнения. Для этого разделим обе части уравнения (4.3.35) на ym:

ym y+p(x) y1m = q(x) (ym=0?)

(4.3.36)

Введем теперь новую неизвестную функцию z, положив

 

1m

 

1

 

 

1m

(4.3.37)

y

= z y = z

.

Тогда

 

 

 

(1m)ymy=z.

(4.3.38)

 

Поэтому, умножая обе части уравнения (4.3.36) на 1 m и выполняя

подстановку (4.3.37), приходим к линейному уравнению

 

z+ (1m)p(x)z = (1m)q(x)

(4.3.39)

Интегрируя это уравнение и возвращаясь к переменной y, получим

общее решение уравнения Бернулли в виде

1

 

 

 

 

y ={e(m1) p(x)dx C + (1 m)q(x)e(m1) p(x)dx dx }1m

(4.3.40)

Деля уравнение (4.3.35) на ym, теряем решение y = 0. Очевидно, что это могло случиться лишь при m > 0. Поэтому решение y = 0 также является решением уравнения (4.3.35), которое не содержится в формуле общего решения (4.3.40).

В качестве простейших примеров, иллюстрирующих сказанное, могут служить уравнения y= y2 и y′ = 2y . Для первого из них решение y=0 –

частное, для второго – особое.

Пример. Найти решение уравнения

y′ − 1x y = − 1x y2 .

Это – уравнение Бернулли. Деля обе части на y2, имеем: y2 y′ − 1x y1 = − 1x .

Полагая

y1=z,

получим:

z′ + 1x z = 1x

Интегрируя это (линейное) уравнение, находим: z = 1x (C + x).

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

436

y = x +x C .

Решениями будут также полуоси оси Ox: y=0 (x0).

437

ЛЕКЦИЯ 4.4. УРАВНЕНИЯ В ПОЛНЫХ ДИФФЕРЕНЦИАЛАХ. УРАВНЕНИЯ ВЫСШИХ ПОРЯДКОВ. ОСНОВНЫЕ ПОНЯТИЯ И ОПРЕДЕЛЕНИЯ.

4.4.1. Уравнения в полных дифференциалах

 

Определение 1. Уравнение

 

M(x, y) dx + N(x, y) dy = 0,

(4.4.1)

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

некоторой функции U от x и y, т. е.

 

M(x, y) dx + N(x, y) dy = dU

(4.4.2)

называется уравнением в полных дифференциалах.

 

Относительно функций M и N мы будем предполагать, что они непрерывны по обеим переменным в некоторой односвязной области и ни в

одной точке этой области не обращаются одновременно в нуль.

 

Уравнение в полных дифференциалах можно записать так:

 

dU = 0.

(4.4.3)

Поэтому общий интеграл его имеет вид

 

U(x,y)=C.

(4.4.4)

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

Поэтому возникают два вопроса: 1) как узнать по виду уравнения (4.4.1), является ли оно уравнением в полных дифференциалах? 2) В случае положительного ответа на первый вопрос, как построить функцию U и, следовательно, общий интеграл уравнения (4.4.1)?

Предположим, что функции M(x,y) и N(x,y) имеют непрерывные частные производные соответственно по y и по x. Пусть левая часть уравнения (4.4.1) представляет собою полный дифференциал, т. е.

M (x, y)dx + N(x, y)dy dU U dx +

U dy .

 

x

y

Это равносильно тому, что имеют место тождества

 

U = M (x, y),

U = N(x, y).

(4.4.5)

x

y

 

Дифференцируя первое из этих тождеств по y, а второе по x, получаем тождества:

2U

=

M

,

2U

=

N

;

(4.4.6)

xy

y

yx

x

 

 

 

 

 

левые части полученных тождеств равны между собою, а тогда равны и правые, т. е.

438

 

 

 

M

N .

(4.4.7)

y

 

x

 

Условие (4.4.7) является необходимым для того, чтобы левая часть уравнения (4.4.1) была полным дифференциалом. Покажем, что это условие является достаточным, если область односвязна. Действительно, пусть условие (4.4.6) выполнено. Покажем, что тогда существует функция U, удовлетворяющая соотношению (4.4.2) или, равенствам (4.4.5).

Будем исходить из первого из равенств (4.4.5):

U

= M (x, y) .

(4.4.8)

x

 

 

Нетрудно убедиться, что ему удовлетворяет функция

x

 

U (x, y) = M (x, y)d x+ ϕ(y) ,

(4.4.9)

x0

 

где ϕ (y) – произвольная функция от y, которую мы будем считать дифференцируемой и выберем ее так, чтобы функции (4.4.9) удовлетворяла и второму из равенств (4.4.5), т. е. чтобы

U

 

x

 

 

 

 

 

 

M (x, y)d x+ ϕ (y) = N(x, y)

y

 

y x

 

или

 

 

 

0

 

 

 

 

 

 

x

M (x, y)

 

y

 

dx + ϕ (y) = N(x, y)

 

 

 

 

 

x0

 

 

 

 

 

Используя условие (4.4.7), перепишем это равенство так:

x N(x, y) dx + ϕ(y) = N(x, y)

x0 x

Выполняя интегрирование, получаем:

N(x,y) N(x0,y) + ϕ′(y) = N(x,y)

откуда

ϕ′(y) = N(x0,y),

следовательно

y

ϕ(y) = N(x0, y)dy + C

y0

(4.4.10)

(4.4.11)

(4.4.12)

(4.4.13)

где C– уже произвольная постоянная. Подставляя найденное

выражение функции ϕ (y) в формулу (4.4.9), получаем искомую функцию

U(x, y):