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

Дифференциальные уравнения. Учебное пособие-1

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

Решение. 1) Перейдем от оригиналов к изображениям. Воспользуемся свойством 5° и таблицей изображений:

 

 

 

 

 

 

 

p

2

X( p) p x(0)

 

2

X( p) p 1, cos 2t

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2 4 .

x(t) X ( p) ,

x( )t

 

 

x (0) p

 

 

 

Получили

алгебраическое

 

уравнение

для

изображений:

p2 X ( p) p 1 4X ( p)

p

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

p2 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) Решим полученное уравнение:

 

 

 

 

 

 

 

 

 

 

 

 

 

p

2

X ( p) 4X ( p)

 

p

 

 

p 1

 

X ( p)

 

 

p

 

p

 

 

1

.

 

p2 4

 

p2 4 2

p2 4

p2 4

3) Воспользуемся таблицей изображений и найдем оригинал для функции

X ( p) : x(t) 1 t sin 2t cos 2t 1 sin 2t .

 

4

2

Таким

образом,

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

x(t)

1 t sin 2t cos 2t 1 sin 2t .

 

4

2

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

Пример 2.15. Решить задачу Коши:

x x 2 y 9t,

x(0) 1,

y 2x y 4et ,

y(0) 2..

Решение. 1) Перейдем от оригиналов к изображениям:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x(t)

 

 

 

 

 

 

pX ( p) x(0)

pX ( p) 1,

 

 

X ( p) , y(t) Y ( p) , x (t)

 

 

 

 

 

 

 

 

 

pY ( p) y(0) pY ( p)

2 ,

 

 

 

 

 

 

9

 

4e

t

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y( )t

 

9t

 

 

 

p2 ,

 

 

 

 

 

p 1.

 

 

Запишем систему уравнений для изображений:

 

 

 

 

 

 

 

 

 

 

 

pX ( p) 1 X ( p) 2Y ( p)

 

9

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pY ( p) 2 2X ( p) Y ( p)

4

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) Решим систему уравнений для изображений.

 

 

 

 

 

 

 

 

 

 

 

Из первого уравнения выразим

 

X ( p) 2Y ( p)

 

1

 

 

9

 

 

и подставим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p 1

 

 

 

 

p 1

 

 

 

p2 ( p 1)

во второе: pY ( p) 2 4Y ( p)

 

 

2

 

 

 

 

18

 

 

Y ( p)

 

4

 

.

 

 

p 1

 

p2 ( p 1)

 

p 1

 

 

 

 

 

 

 

p 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

51

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выразим Y ( p) :

Y ( p)

 

2 p3

4 p2

18

. Разложим полученную дробь на

p2 ( p 1)( p 3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

простейшие: Y ( p)

 

4

 

 

 

6

 

 

 

 

4

 

 

 

 

2

 

 

.

 

 

 

 

 

 

 

 

p

p2

 

p 1

p 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставим Y ( p) в выражение для X ( p) :

 

 

 

 

 

 

 

 

 

 

X ( p)

 

8

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

4

 

1

 

 

9

 

 

p( p 1)

p2 ( p 1)

 

( p 1)(p 1)

( p 3)(p 1)

p 1

p2 ( p 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

1

 

 

 

 

3

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

4

 

.

 

 

 

 

 

 

p2 1

 

p 1

 

p2 ( p 1)

 

 

p( p 1)

 

( p 3)(p 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Последние три слагаемых разложим на простейшие и упростим выражение

для

X ( p) : X ( p)

5

 

3

 

6

 

 

2

 

8

 

.

p

p2

p 1

p 3

p2 1

 

 

 

 

 

 

 

3) С помощью таблицы изображений найдем оригиналы для функций

Y ( p) и X ( p) : y(t) 4 6t 4e t 2e3t ,

x(t) 5 3t 6et 2e3t 8sh t

или,

учитывая, что

x(t) 5 3t 2et 2e3t 4e t .

 

 

Таким образом, решением

задачи Коши

x(t) 5 3t 2et 2e3t 4e t

и y(t) 4 6t 4e t

sh t et e t , получим

2

являются функции

2e3t .

Применение формулы Дюамеля

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

an x(n) (t) an 1x(n 1) (t) ... a1x (t) a0 x(t) f (t) ,

(1.43)

при нулевых начальных условиях:

 

(n 1)

 

 

 

... x

(0) 0 .

(1.44)

x(0) x (0)

 

 

Рассмотрим уравнение с такой же левой частью и правой частью,

равной 1, при тех же начальных условиях:

 

 

 

anh(n) (t) an 1h(n 1) (t) ... a1h (t) a0h(t) 1,

(1.45)

 

... h

(n 1)

(0) 0.

 

h(0) h (0)

 

 

 

Решая задачу (1.45), найдем функцию h(t) . Тогда решение x(t)

исходной

задачи можно выразить через h(t) и f (t)

с помощью формулы Дюамеля.

С учетом нулевых начальных условий получим:

52

t

 

 

x(t) h ( ) f (t )d ,

 

0

 

t

 

 

 

x(t) h (t ) f ( )d ,

(1.46)

0

t

x(t) h(t ) f ( )d f (0)h(t),

0

t

x(t) h( ) f (t )d f (0)h(t).

0

Замечание. Метод решения дифференциальных уравнений, основанный на формуле Дюамеля, применяют, как правило, в тех случаях, когда возникают трудности при нахождении изображения правой части f (t) уравнения (1.43), а также при необходимости многократного решения

задачи (1.43)-(1.44) для различных функций f (t) .

Пример 2.16. Решить задачу Коши:

 

 

 

x

 

3x

 

2x

2e

3t

,

 

x(0) 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

(0) 0 .

 

h (0) 0 и

Решение. Составим уравнение h 3h 2h 1,

 

 

h(0) 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решим его:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

2

H ( p) ,

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1) h(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p ;

 

 

H ( p) , h( )t pH ( p) ,

 

h( )t

 

 

 

 

 

 

 

 

p2H ( p) 3pH( p) 2H ( p)

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) Выразим H ( p) из уравнения: H ( p)

 

1

 

 

 

 

 

.

 

 

p( p 1)( p 2)

 

 

3) Разложим

 

полученную

 

 

 

 

 

дробь

 

 

 

 

на

простейшие

 

 

1

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

H ( p)

 

 

 

 

 

 

 

 

 

 

 

 

 

и с помощью таблицы

p( p 1)( p 2)

2 p

p 1

2 p 2

изображений

 

найдем

оригиналы для

 

 

 

функции

H ( p) :

h(t) 12 et 12 e2t .

4)

 

 

t

e

2t

и воспользуемся одной из формул (1.46):

Найдем h (t) e

 

 

 

t

 

 

 

t

 

 

2

 

3(t )

t

3t 2

 

3t

 

 

 

 

 

 

 

e

2e

d 2e

2e

d

x(t) h ( ) f (t )d

e

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

0

 

 

 

 

e3t 2 2e3t

t

et 2e2t

e3t .

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

53

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задания для типового расчета

Задания по подразделам 1.1 – 2.5 рекомендуется выполнять студентам из [13]. По подразделу 4.4 задания в сборнике [13] не представлены.

В качестве заданий по подразделу 4.4 предлагаются задания 23, 24, 26 раздела 1 сборника [29].

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

[13], [29].

1. Приближенный метод решения с помощью рядов (подраздел 2.6)

При указанных начальных условиях найти три первых члена разложения в степенной ряд функции y(x) , являющейся решением

заданного дифференциального уравнения.

1.

y x 2 y 2 e x , y(0) 0 .

2.

y 2x2 3y cos x 2 , y(0) 0 .

3.y 2x 3 ln y y , y(0) 1.

4.

y x3 y 2 e x ,

y(0) 1.

5.

y x 2 y e y x ,

y(0) 0 .

6.y sin 2x xy , y(0) 1.

7. y e 2 x y 2 ,

y(0) 0 .

8.y cos x x e y , y(0) 0 .

9.y xe x y 2 1, y(0) 0 .

10.y tgx xy 2 e x , y(0) 0 .

11.y 2x3 y 2 2x , y(0) 1.

12.y ln(x 1) e y , y(0) 0 .

13.y x 2 sin y 1, y(0) 0 .

14.y xe x ln y , y(0) 1.

15.y sin 2x cos y , y(0) 0 .

16.y x 2 y ln y y , y(0) 1.

17.y e 2 x 3y 2 , y(0) 0 .

18.y cos 2x x y 2 , y(0) 0 .

19.y cos x y 2 , y(0) 1.

20.y y y 2 , y(0) 3 .

21.y 2e y xy , y(0) 0 .

54

22. y x x 2 y 2 , y(0) 5 .

23.y sin x y 2 , y(0) 1.

24.y x 2 y 2 , y(0) 2.

25.y 2e y xy , y(0) 0 .

2.Системы дифференциальных уравнений (подраздел 3)

Найти общее решение систем уравнений методом Эйлера:

x 2x 19 y 30z 1. y 5 y 12z

z 2 y 5z

x x 2 y 12z 3. y 4 y 3z

z 5 y 6z

x x 8 y 23z

5. y 5 y 7zz 3y z

x x y 16z 7. y y 15z

z y 3z

x 6x 9 y 7z 9. y 3y 2z

z 4 y zx 4x 3y

11. y x

z 2x y

x y

13. y 3x 4 yz 2x y 2z

x 5x 6 y 3z 15. y x z

z x 2 y z

x 3x 2 y 2. y 16x y

z 15x 7 y 14zx 5x 7 y

4. y 3x y

z 12x 6 y 3z

x 4x 5z

6. y 7x 2 y 9zz 3x 6z

x 3x 11y 7z 8. y 5 y 4z

z 8 y zx 5x 21z

10. y 21x 2 y 16zz x z

x 5x 2 y 12. y 3x

z 6x 3y

x 4x 5 y 7z

14.y x 4 y 9zz 4x 5zx 7x

16.y 10x 19 y 10zz 12x 24 y 13z

55

 

x 7 y 4z

17.

y y

 

 

 

z x 13y

 

x 2x y 2z

19.

y 4 y 3z

 

 

 

z y

 

x 5x 4z

21.

y 9x 4 y z

 

 

 

z 7x 5 y 4z

 

x 13x 24 y 12z

23.

y 10x 19 y 10z 24.

 

 

 

z 7z

 

x 2x y 3z

 

 

25.

y 2 y

 

 

 

z x 4z.

x 5x 2 y 3z 18. y 4x 5 y 4zz 6x 4 y 4z

x z

20. y 2x 2 y zz 3x 4zx 4x 7 y 5z

22. y 4x 5 yz x 9 y 4zx 7x

y 12x 13y 24zz 10x 10 y 19z

Следующие системы дифференциальных уравнений решить методом исключения:

1. x 2x y 2ty 2x 4 y

x x y 2 sin t 3. y 2x 3y

x 5x 4 y 9

5. y 2x 11y 9 3et

12x y 2x 0 7. x y 2t

x x y

9.y x y 2 sin t

11.x 5x 2 y 2ety 10x 4 yx x y

13.y 3x 3x 2e 4t

 

 

2

y

 

2.

x 3x

3

 

 

 

 

 

 

 

 

2

 

y 3x 5 y t

 

x x y

4.y 17x 9 y cos t

x 3x 2 y 2

6.y 2x y 3t

8.36x y 2x e 2tx y 0

10.x x 3y 3 cos t

y 2x 6 y

12.x 4x 10 y 3 e3ty 2x 5 y

14.x 4x 2 y

y x 2 y sin t

56

 

x

5x 3y cos t

 

x 3x 2 y

15.

 

 

2

 

16.

 

 

x 3y

 

 

y

3

 

y x y t 2

 

 

 

 

 

 

17.

x 9x 17 y

18. x 11x 2 y 1 t

 

y x y 2t

 

y 4x 5 y

19.

x x 2 y 1

20. y x 3t

 

y 2x 3y 4t

 

12 y x 2 y 0

21.

y

x 0

22.

x y x e3t

 

 

 

 

 

 

36 y x 2 y 2 sin 2t

 

y y x

23.

x 6x 2 y

24. x 4x 10 y

 

y 3x y 2e 4t

 

y 2x 5 y sin t

25.x 5x 2 y

y 10x 4 y 2 cos3t.

3. Метод малого параметра (подраздел 2.6)

Построить два первых члена разложения для решений следующих д. у. ( 1; если уравнение не содержит , то y 1):

1.

y y y 3 2

14.

y y e y

1

2.

y 9 y y 2

15.

y y y 2

3.

y e y

16.

y 9 y e y 1

4.

y 4 y y 2 1

17.

y 4 y y 3

5.

y (2 y 2 y 4 )

18.

y y 2

y 4

6.

y sin y

19.

y 2(e

y 2

1)

7.

y cos y

 

20.

y 2 cos y

8.

y e 2 y 3

21.

y y e 2 y 1

9.

y e3 y

22.

y y y 3

10.

y y cos( x) 0

 

y y cos 2 ( x) 0

11.

y y(1 sin( x)) 0

23.

12.

y yy 1

24.

y y sin y

13.

y sin(x y)

25.

y y cos y .

57

РАЗДЕЛ 2. ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ С ЧАСТНЫМИ ПРОИЗВОДНЫМИ

ТЕОРЕТИЧЕСКИЙ МАТЕРИАЛ

1.Основные понятия

Определение. Уравнение, связывающее неизвестную функцию u x1,..., xn , независимые переменные x1,..., xn и частные производные от

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

Оно имеет вид

 

 

 

 

 

 

 

u

 

u

 

 

u

 

 

 

F x

, x

 

,..., x

n

, u,

,...,

,

 

 

0 ,

(2.1)

 

 

 

xk1

... xkn

 

1

 

2

 

 

x1

xn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

n

 

 

 

где F – заданная функция своих аргументов.

Порядок старшей частной производной, входящей в уравнение (2.1),

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

Общий вид уравнения с частными производными первого порядка с двумя независимыми переменными x и y

 

u

, u

 

(2.2)

F x, y, u,

0 .

 

x

y

 

 

 

 

 

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

 

u

 

u

 

2u

 

2u

 

2u

 

 

F x, y, u,

 

,

 

,

 

2

,

 

,

 

2

 

0 .

(2.3)

 

 

 

 

 

 

x

 

y

 

x

 

x y

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Например, уравнение

A x, y,u,ux ,uy

2u

B x, y,u,ux ,uy

2u

 

x2

x y

 

 

 

 

(2.4)

C x, y,u,ux ,uy 2u f x, y,u,ux ,uy 0

y2

есть квазилинейное уравнение второго порядка.

58

Замечание. Для сокращения записи иногда частные производные

будем

обозначать

индексами

снизу,

например,

u

x

 

u

,

u

xy

 

2u

,

u

yy

 

2u

, и т. д.

x

x y

y2

 

 

 

 

 

 

 

 

 

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

Например, уравнение

A x, y 2u 2B x, y

 

2u

C x, y 2u

 

 

x y

x2

 

 

y2

(2.5)

D x, y u

 

 

 

 

E x, y

u G x, y u F x, y

x

 

y

 

 

есть линейное уравнение второго порядка относительно неизвестной функции u x, y .

Определение. Решением уравнения с частными производными (2.1) называется всякая функция u u x1,..., xn , которая, будучи подставлена в

уравнение, обращает это уравнение в тождество.

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

1) при изучении различных видов волн – упругих, звуковых, электромагнитных, а также других колебательных явлений мы приходим к

волновому уравнению

2u

c2

 

2u

 

2u

 

2u

,

(2.6)

 

2

 

 

2

 

2

 

2

 

t

 

 

x

 

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где с – скорость распространения волны в данной среде; 2) процессы распространения тепла в однородном изотропном теле,

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

u

a2

 

2u

 

2u

 

2u

;

(2.7)

 

 

 

2

 

2

 

2

 

t

 

 

x

 

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

3) при рассмотрении установившегося теплового состояния в однородном изотропном теле мы приходим к уравнению Пуассона

2u

 

2u

 

2u

f (x, y, z) .

(2.8)

x2

y2

z2

 

 

 

 

59

При отсутствии источников тепла внутри тела уравнение (2.8) переходит в уравнение Лапласа, которое описывает стационарное распределение температур

2u

 

2u

 

2u

0 .

(2.9)

x2

y2

z2

 

 

 

 

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

Уравнения (2.6)–(2.9) называют основными уравнениями математической физики.

Каждое из уравнений (2.6)–(2.9) имеет бесчисленное множество частных решений. При решении конкретной физической задачи необходимо из всех этих решений выбрать то, которое удовлетворяет некоторым дополнительным условиям, вытекающим из ее физического смысла. Такими дополнительными условиями чаще всего являются так называемые граничные условия, т. е. условия, заданные на границе рассматриваемой среды, и начальные условия, относящиеся к одному какому-нибудь моменту времени, с которого начинается изучение данного физического явления.

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

1) решение должно существовать, 2) решение должно быть единственным и 3) решение должно быть устойчивым. Это значит, что малые изменения любого из данных задачи должны вызывать соответственно малые изменения решения. Задача, удовлетворяющая всем трем требованиям, называется корректно поставленной задачей.

2.Классификация уравнений второго порядка

2.1.Классификация уравнений с частными производными 2-го порядка

Рассмотрим уравнение второго порядка

 

 

 

 

 

 

 

 

 

 

n

 

 

 

2

u

 

 

 

 

u

 

 

u

 

 

 

 

aij x1,..., xn

 

 

 

 

 

 

 

 

0 .

 

x

x

 

 

 

x

,...,

x

 

 

(2.10)

 

f x1,..., xn , u,

 

n

 

i, j 1

 

 

i

 

 

j

 

1

 

 

 

D

 

Коэффициенты aij

 

заданные

функции

в области

 

пространства

x1,..., xn , причем

aij

a ji .

Все

функции

и

независимые

переменные

считаем вещественными.

60

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