Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Моделирование процессов средствами дифференциальных уравнений. Учебное пособие.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
808 Кб
Скачать

Глава 4

КРАЕВЫЕ ЗАДАЧИ И СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

1.КРАЕВЫЕ ЗАДАЧИ

1.1.Как отмечалось выше, общее решение линейного дифференциального уравнения второго порядка содержит две произвольных постоянных (две «степени свободы»). Частное решение может быть выделено из общего путем задания двух начальных условий. Например, в механике это может быть задание положения движущегося объекта в начальный момент и начальной скорости. Однако, может быть задано также положение объекта в два различных момента времени. Например, процесс механических колебаний объекта массы m относительно положения равновесия описывается уравнением

my′′+ hy′+ ky = f (t) ,

(19)

где y = y(t) – отклонение в момент t точки от положения равновесия;

h – коэффициент трения; k – коэффициент упругости восстанавливающей силы; f (t) – внешняя сила. Если задать положения α и β объекта в мо-

менты соответственно τ и τ

y (τ) = α

(20)y (τ ) ,

то приходим к так называемой краевой задаче (19)–(20). Подобная математическая модель возникает и в задаче об электрических колебаниях и др.

1.2. В общем случае краевая задача для ЛНУ y′′+ py′+qy = f (t)

с условиями (20) может быть решена путём выделения из общего решения

y = C1y1(t) +C2 y2 (t) + yч(t)

(21)

того частного решения, который удовлетворяет (20). Подставляя t = τ и t = τ в соотношение (21), получаем систему алгебраических линейных уравнений

C1 y1(τ) +C2 y2 (τ) = α− yч(τ)C1 y1(τ ) +C2 y2 (τ ) =β− yч(τ )

88

для нахождения соответствующих значений постоянных C1 и C2 . Най-

денные значения следует затем подставить в (21). Тем самым будет получено решение краевой задачи.

Пример 1. Механические колебания материальной точки описываются уравнением

у′′+8у′+17 у = е4t ,

причём положение точки в начальный момент и в момент t =1 заданы: y (0) = 0, y (1) = 0 .

Определить отклонение y (t) точки от положения равновесия в лю-

бой момент времени t.

Решение. На первом этапе мы ищем общее решение линейного неоднородного уравнения. Для этого рассмотрим соответствующее ЛОУ. Запишем характеристическое уравнение

λ2 +8λу +17 = 0

и найдём его корни λ1,2 =−4 ±i . Поскольку корни оказались комплексносопряжёнными, то

у

0

= е4t (C

cost +C

2

sin t) .

 

 

1

 

 

 

 

Поскольку контрольное число S =−4 не совпадает ни с одним из кор-

ней λ1,2 , то частное решение неоднородного уравнения имеет вид

 

 

 

уч

= Ме4t .

 

Подставляя производные

 

 

 

 

 

 

 

= −4Ме

4t

,

′′

 

 

4t

уч

 

 

уч =16 Ме

 

в неоднородное уравнение, получаем

16Me4t 32Me4t +17Me4t = e4t ,

откуда М =1 так что уч =е4t . Теперь мы можем записать общее решение ЛНУ

у = у0 + уч = е4t (1+C1 cost +C2 sin t) .

На втором этапе решения подбираем значения постоянных C1 и C2 , так чтобы соответствующее решение удовлетворяло краевым условиям.

89

Для этого значения t = 0 и

y = 0, t =1 и

 

y = 0 подставим в найденное

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

 

 

 

 

 

 

 

 

 

 

1

+С

+0

= 0

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

4

(1

+С

 

cos1+C

 

sin1) = 0,

е

 

1

2

 

 

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

C

= −1

 

 

 

 

 

 

1

 

1+cos1

 

 

 

 

 

 

 

 

 

C2

=

 

 

 

 

.

 

 

 

sin1

 

 

 

 

 

 

 

 

 

Окончательно, найдено решение

у(t) = e4t 1cost + cos11sin t

sin1

которое и определяет искомое отклонение в любой момент времени t. Вычисляя приближённо (с точностью до 0,1) постоянный коэффициент

cos11 , получим sin1

 

 

 

 

 

у(t) e4t

 

(1cost 0,6sin t) .

 

 

 

Пример 2. Процесс колебания материальной точки массой m под

действием силы упругости

Fу = −ky ,

силы

сопротивления

среды

Fc

= −hy

и внешней силы F(t) , где t – время, а y(t) – отклонение от со-

стояния

 

равновесия y = 0, может быть

описан уравнением

вида

 

′′

 

 

 

h

 

 

k

 

 

 

F(t)

 

 

y

+ py

+ qy = f (t) . Здесь p = m , q = m ,

f (t)

=

 

m . Найти закон движе-

 

 

 

ния точки, если известны значения p = 2, q = 10,

f (t) = sin 3t , а также ко-

ординаты точки в моменты времени t0 = 0 и t1 = 5: y(0) = 0, y(5) = 0.

 

 

 

Решение. Требуется решить уравнение

 

 

 

 

 

у′′+2у′+10у = sin 3t ,

причём положение точки в начальный момент и в момент t1 =1 заданы: y (0) = 0, y (5) = 0 .

Найдём общее решение ЛНУ. Составляем характеристическое уравнение

λ2 + 2λ+10 = 0 .

90

Это уравнение для соответствующего ЛОУ имеет комплексные корни λ1,2 =−1±3i , поэтому общее решение ЛОУ получаем в виде

у

0

= еt (C cos3t +C

2

sin 3t) .

 

1

 

Поскольку контрольное число S =3i не совпадает ни с одним из корней, то частное решение ЛНУ ищем в виде

уч = М cos3t + N sin 3t .

Находя

уч = −3М sin 3t +3N cos3t , уч′ = −9M cos3t 9N sin 3t

и подставляя результаты в ЛНУ, получаем

9M cos3t 9N sin3t +2(3M sin3t +3N cos3t) +10(M cos3t + N sin3t) =sin3t; cos3t(9M +6N +10M ) +sin3t(9N 6M +10N) =sin3t;

cos3t(M +6N) +sin3t(6M + N) =sin3t.

Приравнивая коэффициенты при одинаковых функциях в левой и правой частях, получаем

 

 

 

 

 

 

 

 

M +6N = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6M + N

 

 

 

 

 

 

откуда

М = −

6

, N =

1

так что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

37

 

 

 

37

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

уч =−

6

cos3t +

 

 

1

sin 3t .

 

 

 

 

 

 

 

 

 

 

 

 

37

37

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, общее решение ЛНУ

 

 

 

 

 

 

 

 

 

у = у + у

ч

=et (C cos3t +C

2

sin 3t)

6

cos3t +

1

sin 3t .

 

 

 

 

0

 

 

1

 

 

 

 

 

 

37

 

 

37

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь подставим краевые условия: t = 0 и y = 0, t = 5 и y = 0 :

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С1

 

 

= 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

37

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

е

(С1 cos15 +C2

sin15)

 

 

 

cos15 +

 

sin15 = 0.

 

37

37

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

91

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