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

Учебное пособие 800446

.pdf
Скачиваний:
7
Добавлен:
01.05.2022
Размер:
2.77 Mб
Скачать

y 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p(x)dx,

 

ln

 

y

 

ln

 

C

 

 

 

p(x)dx,

(C 0),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p(x)dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

(C 0),

 

 

 

y Ce

 

 

 

 

 

 

 

y Cexp p(x)dx ,

 

 

C .

(2.44)

y

y0 M0

0 a x0 b x

Рис. 6. Полоса единственности решения линейного уравнения

В формуле (2.44) содержатся все решения однородного уравнения (2.42). Заменим в формуле (2.44) неопределенный интеграл определенным интегралом

с фиксированным нижним пределом

x0 (a,b)

и переменным верхним

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

Получим

 

 

 

 

x

 

 

 

 

 

 

y Cexp p(x)dx .

(2.45)

 

x0

 

 

 

 

 

Обозначим y(x0) y0

и положим в этой формуле x x0 . Тогда получим

C y0. Подставив это значение C в формулу (2.45), получим:

 

 

 

x

 

 

 

 

 

 

 

y y0 exp

p(x)dx .

(2.46)

 

 

x0

 

 

 

 

 

 

Если y0 - произвольно, то формула (2.46) является общим решением линейного однородного уравнения (2.42) в форме Коши. Если же y0 - фиксировано, то (2.46) представляет собой решение уравнения (2.42) с начальными данными x0, y0.

Из формулы общего решения (2.46) однородного уравнения видно, что если y0 0, то решение не обращается в нуль ни в одной точке.

30

Пример.

Найти общее решение уравнения

 

 

 

y

 

x

 

 

y 0.

 

 

 

 

 

 

 

 

 

1 x2

 

 

 

 

 

 

 

x

 

 

Решение.

В этом случае

p(x)

 

 

 

 

- функция непрерывная в

 

 

 

 

 

 

 

1 x2

 

 

 

 

 

 

 

 

интервале ( 1,1). Воспользуемся формулой (2.44), предварительно вычислив

p(x)dx

 

x

 

 

 

t 1 x2

 

 

1

 

dt

 

 

 

 

 

;

 

 

dx

 

 

 

 

 

1 x2

 

 

 

 

t

 

 

 

 

 

 

 

 

1 x2

 

 

 

 

 

 

dt 2xdx

2

 

 

t

 

 

 

 

 

 

 

 

y Ce

1 x2

,

C .

Полученная формула дает общее решение уравнения в области:

1 x 1,

y .

Пример. Найти решение задачи Коши для линейного однородного уравнения:

y ycosx 0,

y(0) 1.

Решение. Воспользовавшись формулой (2.46), получаем:

 

 

x

 

 

 

sin x

 

 

 

 

( cosx)dx

 

e

.

 

 

 

y 1 exp

 

 

 

 

 

0

 

 

 

 

 

Рассмотрим свойства решений однородного линейного уравнения.

1. Если y1(x) есть частное решение уравнения (2.42), то функция y Cy1

, где C - произвольная

постоянная,

тоже

является

решением

этого

уравнения.

как y1 - решение

 

 

 

 

Действительно, так

уравнения,

то имеет

место

тождество

y1 p(x)y1 0

(a x b).

 

 

 

 

 

Подставим в уравнение (2.42) y Cy1 и примем во внимание записанное

выше тождество. Имеем:

 

 

 

(a x b).

 

(Cy1) p(x)Cy1 C(y1 p(x)y1) 0

 

Следовательно, функция y Cy1 является решением уравнения (2.42).

2. Если y1 - частное решение уравнения (2.42), то формула y Cy1, где C - произвольная постоянная, дает общее решение уравнения (2.42) в области a x b, y .

Действительно, функция y Cy1 является решением уравнения (2.42) (как показано выше) для всех значений параметра C. Рассматривая равенство y Cy1 как уравнение относительно C, видим, что оно разрешимо в области

31

a x b,

y при всех C:C y y1.

Поэтому формула y Cy1 дает

общее решение.

4. В этом разделе рассматривается важнейшее свойство решения неоднородного уравнения – теорема о структуре общего решения линейного неоднородного уравнения.

Теорема. Если yчн - частное решение неоднородного уравнения (2.41),

то общее решение дается формулой

 

yон yоо yчн,

где yоо - общее решение однородного уравнения (2.42)

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

yоо Ce p(x)dx.

Пусть нам известно некоторое частное решение yчн

неоднородного уравнения (2.41), т.е.

 

 

 

q(x)

a x b.

 

yчн p(x)yчн

Введем новую

неизвестную

функцию z по формуле y yчн z.

Подставляя функцию y в уравнение (2.41), получим, требуя, чтобы функция y

удовлетворяла уравнению:

yчн p(x)yчн z p(x)z q(x).

Таким образом, чтобы функция y yчн z являлась решением (2.41) достаточно, чтобы функция z удовлетворяла однородному линейному уравнению, и соответствующему уравнению (2.41):

z p(x)z 0.

Общее решение этого уравнения, как показано в п. 2 настоящего раздела, имеет вид (2.45). Итак, функция z должна являться общим решением однородного уравнения: z yoo. Таким образом, получим:

p(x)dx

( C ).

y yчн Ce

Покажем, что эта формула представляет общее решение уравнения (2.41).

Запишем это решение в виде

 

x

 

p(x)dx

y yчн Ce x0

(2.47)

где, как и в формуле (2.45), неопределенный интеграл заменим определенным

интегралом

с переменным верхним

пределом. Покажем, что

при любых

x0 (a,b)

и

y0 ( , )

частное решение уравнения (2.41), удовлетворяющее

условию

y(x0) y0 ,

принадлежит

семейству (2.47) при

однозначно

определяемой постоянной C.

 

 

Полагая в формуле (2.47) x x0,

y y0, получим:

 

 

 

 

C y0 yчн (x0).

 

 

 

 

 

32

 

Таким образом, формула (2.47) дает общее решение неоднородного линейного уравнения:

yон yчн уоо,

yоо Ce p(x)dx .

(2.48)

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

5. Методы построения общего решения неоднородного линейного уравнения.

В настоящем разделе будут рассмотрены два метода решения уравнения (2.41), первый из которых - метод Лагранжа, применяется к решению не только линейных уравнений или сводящимся к ним.

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

Для нахождения решения уравнения (2.41) применим следующий метод: будем искать решение в том же виде, что и решение соответствующего однородного уравнения (2.44), но будем считать C не постоянной, а неизвестной непрерывно дифференцируемой функцией от x. Таким образом, положим

 

 

y C(x)e p(x)dx.

 

 

 

(2.49)

Требуем, чтобы функция (2.49) удовлетворяла уравнению (2.41):

 

p(x)dx

C(x)p(x)e

p(x)dx

p(x)C(x)e

p(x)dx

q(x).

C (x)e

 

 

 

 

Отсюда получаем следующее дифференциальное уравнение, которому

должна удовлетворять функция C(x):

 

 

 

 

 

 

 

 

 

dC(x)

q(x)e p(x)dx.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

Это уравнение – простейшее дифференциальное уравнение с

разделяющимися переменными. Интегрируя его, имеем:

 

 

 

dC(x) q(x)e p(x)dxdx C,

C(x) C q(x)e p(x)dxdx,

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

( C ).

Подставляя эту функцию в

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

 

 

 

 

 

 

 

 

y e p(x)dx C q(x)e p(x)dxdx,

( C ).

 

(2.50)

Полученная функция по построению является решением уравнения (2.41).

Более того, она является общим решением в полосе

a x b,

y ,

поскольку имеет структуру (2.48).

 

 

 

 

 

 

 

Действительно, записав (2.50) в виде

 

 

 

 

 

 

 

y Ce p(x)dx e p(x)dx q(x)e p(x)dxdx,

 

(2.51)

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

33

решение неоднородного уравнения, так как оно содержится в формуле (2.51) при C 0.

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

Отметим характерное свойство решения линейного уравнения. Из (2.51)

видно, что общее решение линейно зависит от произвольной постоянной. Для нелинейных уравнений такое свойство уже не выполняется.

Другой вид общего решения получается из (2.51) заменой неопределенных интегралов определенными с переменным верхним пределом:

 

x

x

 

s

 

 

 

p(t)dt

p(t)dt x

p(t)dt

 

y Ce x0

e x0

q(s)ex0

ds,

(2.52)

 

 

 

x0

 

 

 

где x0 - произвольное число из (a,b),

x (a,b).

 

 

Пусть требуется решить задачу Коши для уравнения (2.41), то есть

нужно найти

такое решение уравнения,

которое

удовлетворяет условию

y(x0) y0 , где

x0 (a,b),

y0 ( , ). Подставляя

в (2.52)

x x0 , y y0

находим: C y0.

Таким образом, решением задачи Коши является функция

x

x

 

s

 

 

p(t)dt

p(t)dt x

p(t)dt

 

y y0e x0

e x0

q(s)ex0

ds.

(2.53)

 

x0

 

 

З а м е ч а н и е 1. Формула (2.52) показывает, что если функции

p(x) и

q(x) заданы и непрерывны на всей оси

 

( , ),

то и решение с любыми

начальными данными x0 , y0 будет не только непрерывно, но и непрерывно дифференцируемо при всех значениях x. Это означает, что интегральная кривая уравнения, проходящая через любую точку плоскости M0(x0, y0) будет гладкой кривой при x ( , ).

З а м е ч а н и е 2. Если функции p(x) и q(x) непрерывны в (a,b), за исключением отдельных точек, то формула общего решения (2.51) или (2.52) остается справедливой для всех x из (a,b), если интегралы при переходе через точки разрыва не теряют смысл (если интегралы расходятся, то формула

общего решения уже не будет справедливой во всем интервале).

 

 

2

Пример.

Найти

общее решение дифференциального

уравнения

y

y x.

 

 

 

x

 

 

 

 

Решение.

Данное

уравнение является уравнением вида

(2.41), где

 

 

p(x) 2x, q(x) x, то есть линейным неоднородным уравнением. Функция

34

q(x) определена и непрерывна при всех x, а функция p(x) - при x ( ,0) (0, ). Найдем его общее решение методом вариации постоянной.

Решим

соответствующее

ему однородное

z

2

z 0. Разделяя

 

переменные и интегрируя, получим

 

x

 

 

 

 

 

dz

2

dx

,

ln

 

z

 

ln

 

C

 

2ln

 

x

 

 

(C 0),

 

 

 

 

 

 

 

откуда z Cx2

z

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(C 0) и z

0. Общее решение имеет вид

z Cx2, C ( , ).

Общее решение неоднородного уравнения, в соответствии с методом Лагранжа, ищем в виде y C(x)x2.

Находим производную y C (x)x2 2xC(x), которую, вместе с функцией y подставляем в уравнение. Получаем:

 

 

2

 

2

 

 

 

 

 

2

 

 

 

 

 

 

C (x)x

 

2xC(x)

 

 

 

C(x)x

 

x

или

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

C (x)

,

dC(x)

,

 

C(x) ln

 

x

 

C .

 

 

 

 

x

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя найденную функцию в выражение для общего решения,

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

 

y x2(C ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример.

Найти общее

 

решение

уравнения xy 2x2y 1, пользуясь

формулой общего расчета (2.50).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение.

В данном случае p(x) 2x,

 

q(x) 1 x. Подставляя в формулу

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2xdx

1

 

2xdx

 

 

 

x2

 

 

ex

2

 

 

 

 

 

 

 

y e

C

 

 

e

 

dx

 

e

 

C

 

 

 

x

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx .

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

5.2. Метод Бернулли построения общего решения.

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

y u(x)v(x).

(2.54)

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

vdu u dv p(x)uv q(x), dx dx

35

или

 

du

 

dv

 

 

 

v

 

u

 

p(x)v

q(x).

(2.55)

dx

 

 

 

dx

 

 

 

Функцию v(x) выбираем так, чтобы

 

 

 

 

 

dv

p(x)v 0.

 

(2.56)

 

 

 

 

 

 

 

 

dx

 

 

 

Это означает, что на одну из двух функций налагается требование – она должна удовлетворять однородному уравнению, соответствующему данному неоднородному. Таким образом, дифференциальное уравнение (2.55) будет равносильно системе двух дифференциальных уравнений, одно из которых – уравнение (2.56), а другое получается из (2.55) при условии (2.56). Имеем:

dv

p(x)v 0,

 

 

 

 

 

 

 

 

 

 

 

 

(2.57)

dx

 

 

 

du

 

v

 

 

q(x).

 

 

dx

 

Оба уравнения системы (2.57) являются уравнениями с разделяющимися переменными. Для первого из них мы возьмем частное решение

v exp( p(x)dx) (см. формулу (2.48)). Из второго уравнения получаем

du q(x) dx, u q(x)e p(x)dxdx C. v

Подставив найденные функции v(x) и u(x) в формулу (2.54), получим

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

y e p(x)dx

 

 

 

 

 

 

C

 

q(x)e p(x)dxdx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример. Решить уравнение

y sin x ycosx

sin

2 x

.

x2

 

 

 

 

 

Решение. Ищем его в виде y u(x) v(x). Дифференцируя по х, имеем: y u v uv .

Подставим эти выражения в уравнение

 

 

 

sin

2 x

 

u vsin x uv

sin x uvcosx

 

 

(*)

x

2

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

dvsin x vcosx 0. dx

Отсюда

dv

 

cosx

dx,

lnv lnsin x ,

v sin x.

 

 

v

sin x

 

36

Подставим функцию v в уравнение (*) и учитывая равенство нулю подчеркнутой части, получим уравнение для u(x)

 

 

 

du

sin2 x

sin

2 x

.

 

 

 

 

x

2

 

1

 

dx

 

Отсюда u(x)

C,

C

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

 

 

x

 

 

 

 

 

y sin x Csin x. x

Пример. Найти решение задачи Коши

 

y x

2

(2x 1),

x(x 1)y

 

 

 

 

 

 

 

 

 

y(2) 4.

 

 

 

Решение. Ищем общее решение дифференциального уравнения методом Бернулли y uv. Подставляя y и y u v uv в уравнение, получим:

x(x 1)(u v uv ) uv x2(2x 1), x(x 1)u v u x(x 1)v v x2(2x 1).

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

x(x

x(x

Решаем первое уравнение,

1) dv v 0, dx

1)v du x2(2x 1). dx

разделяя переменные

 

dv

 

dx

 

lnv

dx

 

 

,

 

.

v

x(x 1)

x(x 1)

Разлагаем подынтегральную правильную рациональную дробь на сумму простых дробей:

 

1

 

 

A

 

 

B

 

 

 

 

 

 

 

 

 

x(x 1)

 

 

x 1

 

 

 

 

 

 

x

 

 

откуда A(x 1) Bx 1. Полагая x 0, получим A 1;

полагая x 1, получим

B 1. Таким образом, имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lnv

dx

 

 

 

dx

ln

 

x

 

.

 

 

 

 

 

 

 

 

 

x

x 1

 

x 1

 

 

 

Частное решение берем в виде v

 

 

x

. Подставляя эту функцию во

 

 

 

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

 

 

 

 

 

x 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

du

du (2x 1)dx,

u x2 x C.

x2

 

x2(2x 1),

dx

 

 

 

 

 

 

37

 

 

 

 

 

 

 

 

 

 

 

Общее решение имеет вид

 

 

 

 

 

 

 

 

 

 

x

2

 

 

 

2

 

x

y

 

 

(x

 

x C) x

 

C

 

.

 

 

 

 

 

 

 

x 1

 

 

 

 

 

 

x 1

Для построения решения задачи Коши найдем значение C, используя

начальное условие x0 2,

y0 4. Имеем:

 

 

 

 

 

 

4 4 C

2

, C 0.

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

Решением задачи Коши является функция y x2.

Задачи для самостоятельного решения

1.

y ytgx xtgx 1,

y(0) 1.

Ответ: y Csin x x.

 

2.

y sin 2x 2(y cos x).

Ответ:

y ctgx

 

 

1

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cosx

 

3.

y 2xy 3x2 2x4,

y(0) 1.

 

Ответ: y ex2 .

 

4.

dx (x y)dy 0.

 

 

Ответ: x Ce y y 1.

 

5.

y ytgx 2sin x,

 

y(0) 1. Ответ:

 

y

 

 

1

 

(C cos2 x),

y cosx.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cosx

 

6.

ex(y y) 1.

Ответ: y e x(x C).

 

 

 

 

 

 

 

x

2

y

 

2xy ln x 0.

Ответ: y

 

ln x 1

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

7.

 

x

 

x2 .

 

 

 

 

 

 

y y sin x.

Ответ: y Ce

x

1

(cosx sin x).

 

8.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение вида

 

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

 

 

 

 

 

y p(x)y q(x)yv,

 

 

 

 

 

 

 

 

 

(2.58)

где v - любое действительное число, называется уравнением Бернулли. Если v 0, то уравнение Бернулли вырождается в линейное уравнение (2.41); при v 1 уравнение превращается в уравнение с разделяющимися переменными:

y p(x)y q(x)

(v 0);

y (q(x) p(x))y

(v 1).

Уравнение (2.58) также будет являться уравнением с разделяющимися

переменными при p(x) q(x).

 

Ниже мы

рассмотрим случай, когда v 0,

v 1. Относительно функций

p(x) и q(x) будем предполагать, что они

непрерывны при x (a,b) и

 

p(x)

 

 

 

q(x)

 

.

 

 

 

 

 

 

1. Уравнение Бернулли сводится к линейному уравнению подстановкой

38

1

(z y1 v).

 

y z1 v ,

(2.59)

Действительно, применив эту подстановку, найдем производную новой функции

z (1 v)y vy .

Умножив обе части уравнения (2.58) на (1 v)y v (предполагая, что y 0 если v 0), получим дифференциальное уравнение для новой функции

z (1 v)p(x)z (1 v)q(x),

(2.60)

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

1

 

(v 1)

 

p(x)dx

 

(1 v)

 

p(x)dx

 

1 v

 

 

 

C (1 v) q(x)e

dx

 

.

(2.61)

y e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Читателю рекомендуется самостоятельно получить это решение. Таким образом, уравнение Бернулли также интегрируется в квадратурах.

2. Помимо общего решения уравнение Бернулли при v 0 имеет еще очевидное решение y 0 (при v 0 функция y 0 не является решением уравнения). При 0 v 1 решение y 0 не содержится в формуле общего решения. Оно является особым решением, поскольку в каждой точке этого решения нарушается единственность решения задачи Коши.

При v 1 решение y 0 является частным решением,

так как в этом

случае

через точки оси OX не проходит ни одна интегральная

кривая (кроме

y 0,

т.е. самой оси OX). Таким образом, при v 1 во всякой точке оси OX

решение существует и единственно. Формально это решение можно включить в семейство (2.61), положив C .

3. При решении уравнения Бернулли нет необходимости пользоваться готовой формулой (2.61). Проще превратить его в линейное уравнение, используя подстановку (2.59), или же непосредственно использовать метод Бернулли, отыскивая решение уравнения в виде y u v.

Пример. Найти решение уравнения xy y y2 ln x.

Решение. Разделив обе части уравнения на x видим, что это уравнение типа (2.58) при v 2. Разыскиваем его решение в виде произведения двух функций y u v. Подставим в уравнение y u v uv , y u v, получим:

u v uv 1uv u2v2 ln x, x x

39