- •Рекомендуется изучить раздел «Непрерывность функции в точке» пособия Л.И. Магазинников, А.Л. Магазинников Дифференциальное исчисление.
- •Лабораторная работа №4
- •Проверка сходимости числовых рядов
- •Рекомендуется изучить раздел «Числовые ряды» в пособии Л.И. Магазинников Функции комплексного переменного. Ряды. Интегральные преобразования.
- •Рассмотрим пример.
- •Исследовать на сходимость числовой ряд
- •Численное исследование сходимости.
- •Аналитическое исследование сходимости.
- •Численное решение дифференциальных уравнений
- •Рекомендуется изучить раздел «Уравнения высших порядков» в пособии А.А. Ельцов, Т.А. Ельцова Интегральное исчисление. Дифференциальные уравнения.
- •Решение. Составляем характеристическое уравнение
- •Получаем систему двух уравнений
- •Лабораторная работа №9
- •Рекомендуется изучить раздел «Ряды Фурье» в пособии Л.И. Магазинников Функции комплексного переменного. Ряды. Интегральные преобразования.
15
Лабораторная работа №7
Численное решение дифференциальных уравнений
Рекомендуется изучить раздел «Уравнения высших порядков» в пособии А.А. Ельцов, Т.А. Ельцова Интегральное исчисление. Дифференциальные уравнения.
Цель работы: применение численного и аналитического методов решений линейных дифференциальных уравнений.
Рассмотрим пример.
Найдите частное решение дифференциального уравнения, удовлетворяющего указанным условиям:
y′′−5y′+6y = 0, y(0) =1, y′(0) = 0
Решение. Составляем характеристическое уравнение
k2 − 5k + 6 = 0
и находим его корни: |
|
|
|
|
|
|
|
|
|
|
k = 5 ± |
|
= |
5 ±1, k = 2, k |
|
|
|||||
25 −24 |
2 |
= 3. |
||||||||
1,2 |
2 |
|
|
|
2 |
1 |
|
|||
|
|
|
|
|
|
|
||||
Записываем общее решение дифференциального уравнения: |
||||||||||
|
y(x) = C e2 x +C |
e3x . |
|
|
||||||
|
|
|
|
1 |
|
2 |
|
|
|
|
Определяем константы C1 |
и С2. Поскольку y(0) = 1, то: C1 + C2 = 1. |
|||||||||
′ |
|
|
2 x |
+3C2e |
3x |
. Поскольку y'(0) = 0, то: 2C1 + 3C2 = 0. |
||||
Находим производную y (x) = 2C1e |
|
|
Получаем систему двух уравнений
С1 +С2 =1,2С1 +3С2 = 0,
из которой C1 = 3, C2 = −2.
Записываем частное решение дифференциального уравнения: y(x) = 3e2 x −2e3x
Строим график зависимости y(x):
16
Задание.
1.Составить программу для решения линейных дифференциальных уравнений.
2.Используя рассмотренный выше пример, протестировать программу.
3.Построить график получившейся зависимости y(x). Сравнить графики численного и аналитического решений; сделать выводы.
4.Численно и аналитически решить следующие дифференциальные уравнения. Во всех примерах построить графики зависимостей y(x). Сравнить графики численного и аналитического решений; сделать выводы.
ВАРИАНТ 1 |
ВАРИАНТ 2 |
1. |
y′′+6y′+9 y =10sin x, |
1. |
y′′+9 y = 6e3x |
, |
|
|
|
|
|
|
|
||||||||||||||||||
′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|
||||
|
y (0) = y(0) = 0 |
|
|
|
|
y(0) |
|
|
|
|
|
0 |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
= y |
(0) = |
|
|
|
|
|
|
|||||||||||
2. |
y′′−4 y′ = 6x2 +1, |
|
y′′−6y′+9 y = x2 − x +3, |
||||||||||||||||||||||||||
y(0) |
|
|
|
′ |
|
|
|
2. |
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
= 2; y |
(0) = 3 |
|
|
y(0) |
= |
|
|
′ |
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
27 |
|
|
||||||||
|
y′′−4 y′+5y = 2x2ex , |
|
3 , y (0) = |
|
|
||||||||||||||||||||||||
3. |
|
|
′′ |
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|||||||||
|
y(0) |
|
|
|
′ |
|
|
|
|
y |
+3y |
−4 y = (x + 2)e , |
|||||||||||||||||
|
= 2, y |
(0) = 3 |
|
3. |
|
′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
9 |
|
|
|
|
(0) = y(0) = |
|
|
|
|
|
|
||||||||||||
|
y′′+9 y = |
|
, |
|
|
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|||||||
4. |
cos3x |
|
|
y′′+ 4 y = |
|
|
|
|
, |
|
|
|
|
|
|||||||||||||||
|
y(0) |
|
|
|
′ |
|
0 |
|
|
|
sin 2x |
|
|
|
|
|
|||||||||||||
|
=1, y (0) = |
|
4. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
y |
π |
|
|
|
|
|
|
|
π |
= π |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
4 |
= 2, |
y′ |
4 |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
ВАРИАНТ 3 |
|
|
|
|
|
|
|
|
|
|
ВАРИАНТ 4 |
|||||||||||||
1. |
y′′− |
4 y′ = 2xe4 x , |
1. |
y′′−4 y′ = 4xe4 x , |
|
|
|
|
|
||||||||||||||||||||
y(0) |
|
|
′ |
|
|
|
|
y(0) |
|
|
′ |
|
|
|
0 |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
= y (0) = 0 |
|
|
|
|
= y |
(0) = |
|
|
|
|
|
|
||||||||||||||||
2. |
′′ |
|
′ |
−12 y = 8sin 2x, |
2. |
|
′′ |
+ 4 y = sin x, |
|
|
|
|
|
|
|||||||||||||||
y′ + 4 y |
|
y′ |
|
|
|
|
|
|
|||||||||||||||||||||
|
y (0) = y(0) = 0 |
|
|
|
|
y |
(0) = y(0) =1 |
|
|
|
|
|
|
||||||||||||||||
|
y′′−9 y′+18y = |
|
9e3x |
|
y′′−6y′+8y = |
|
|
|
4 |
|
|
, |
|||||||||||||||||
|
|
|
, |
|
|
|
|||||||||||||||||||||||
3. |
|
1+ e−3x |
3. |
2 + e−2 x |
|||||||||||||||||||||||||
|
y(0) = 0, y′(0) = 0 |
|
y(0) =1+3ln 3, y |
′ |
|
|
|
10ln 3 |
|||||||||||||||||||||
|
|
(0) = |
|||||||||||||||||||||||||||
|
y′′+ y = 4ctg x, |
|
|
|
|
y′′+ |
1 |
|
y = |
|
|
|
1 |
|
|
, |
|
||||||||||||
4. |
π |
|
|
4, |
π |
|
|
4. |
π2 |
|
|
π |
2 |
cos |
π |
|
|||||||||||||
|
y |
= |
y′ |
|
= 4 |
|
|
|
|
|
|
|
|
|
|
x |
|
||||||||||||
|
2 |
|
|
|
|
2 |
|
|
|
y(0) |
|
|
|
|
′ |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
= 2, y |
(0) = 0 |
|
|
|
|
17
|
|
|
|
|
ВАРИАНТ 5 |
|
|
|
|
|
ВАРИАНТ 6 |
|
||||||||||||||||||||||||
1. 2 y′′+ y′− y = 3xex , |
|
|
1. |
y′′− y′−6y = xe2 x , |
|
|
|
|||||||||||||||||||||||||||||
|
y(0) |
|
|
′ |
|
|
|
0 |
|
|
|
|
|
|
|
y(0) |
′ |
|
|
0 |
|
|
|
|
|
|
|
|||||||||
|
= y (0) = |
|
|
|
|
|
|
|
|
= y |
(0) = |
|
|
|
|
|
|
|
||||||||||||||||||
2. |
y′′+ y = 2cos x, |
|
|
|
|
|
y′′+ 4 y =12cos2x, |
|
|
|
||||||||||||||||||||||||||
y(0) |
=1, y |
′ |
|
|
|
|
|
|
|
|
|
2. |
|
|
|
|||||||||||||||||||||
|
(0) = 0 |
|
|
y(0) |
′ |
|
|
0 |
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
9e3x |
|
|
|
|
|
= y |
(0) = |
|
|
|
|
|
|
|
|||||||||||||
3. y′′+3y′ = |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
1+e3x |
|
|
|
y′′−6y′+8y = |
|
|
|
|
|
4 |
, |
|
||||||||||||||||||||||||
|
y(0) |
|
|
|
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
−2 x |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3. |
|
|
|
|
|
|
|
|
1+ e |
|
|
|
|||||||||
|
= ln 4, y (0) = 3 −3ln 2 |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
π |
2 |
|
|
|
|
|
|
|
|
y(0) |
|
|
|
|
|
|
|
|
|
′ |
|
|
6ln 2 |
||||
|
|
′′ |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
=1+ 2ln 2, y |
(0) = |
|||||||||||||||||
4. |
y |
+ π |
|
y = sin πx , |
|
|
|
|
|
|
|
|
|
9 |
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
1 |
|
|
|
|
|
|
1 |
|
|
π2 |
|
|
′′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
y |
=1, |
|
|
= |
|
y |
+ |
9 y = sin3x , |
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
y′ |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
2 |
4. |
y |
π |
|
|
|
|
π |
|
= |
3π |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
= 4, y′ |
|
|
|
|
|
2 |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
ВАРИАНТ 7 |
|
|
|
|
|
ВАРИАНТ 8 |
|
||||||||||||||||||||||||
1. y′′−4 y = xe2 x , |
|
|
|
|
|
|
|
1. |
y′′+ 4 y = e−2 x , |
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
y(0) |
|
|
′ |
|
|
|
0 |
|
|
|
|
|
|
|
y(0) |
′ |
|
|
1 |
|
|
|
|
|
|
|
|||||||||
|
= y (0) = |
|
|
|
|
|
|
|
|
= y |
(0) = |
|
|
|
|
|
|
|
||||||||||||||||||
2. |
y′′−4 y′+13y = 26x +5, |
|
y |
′′ |
+ y = x |
3 |
+1, |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
y(0) |
=1, y |
′ |
|
= 0 |
|
|
2. |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
(0) |
|
|
|
|
y(0) |
′ |
|
|
0 |
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= y |
(0) = |
|
|
|
|
|
|
|
||||||||
|
y′′− |
3y′ = |
|
|
9e |
−3x |
|
, |
|
|
|
y′′+6y′+8y = 4e−2 x (2 +e2 x ), |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
3. |
||||||||||||||||||||||||||||
3. |
|
3 + e |
−3x |
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y (0)= 0, y′(0)= 0 |
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y(0) |
|
|
|
|
|
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
= 4ln 4, y |
(0) = 9ln 4 −3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
y′′+ |
4 y = 8ctg 2x, |
|
|
|
y′′+ π2 y |
= |
π2 |
|
, |
|
|
|
|||||||||||||||||||||||
4. |
|
|
4. |
cos |
πx |
|
|
|
||||||||||||||||||||||||||||
|
π |
|
|
|
|
π |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
y(0) |
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
y |
|
= 5, y′ |
|
|
= 4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
4 |
|
|
|
= 3, y |
(0) = 0 |
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
18
Лабораторная работа №8
Применение степенных рядов к решению дифференциальных уравнений
В случаях, когда интегрирование дифференциального уравнения в элементарных функциях невозможно или способ его решения слишком сложен, решение такого уравне-
∞ |
|
ния следует искать в виде ряда Тейлора y = ∑cn (x − x0 )n . Коэффициенты ряда cn |
нахо- |
n=0 |
|
дят подстановкой ряда в уравнение и приравниванием коэффициентов при одинаковых степенях (x − x0 ) в обеих частях полученного равенства. Если удаётся найти все коэффи-
циенты ряда, то полученный ряд служит решением во всей своей области сходимости. Этим способом можно интегрировать линейные дифференциальные уравнения с перемен-
ными коэффициентами. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
′′ |
+ xy |
′ |
+ y = x cos x; |
y(0) = 0, y |
(0) |
=1. |
||||||||||||||||||||||||||||||||||
Пример 10. Решить задачу Коши для уравнения y |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
′ |
|
|||
Сначала разложим правую часть в степенной ряд по степеням x, т.к. у нас x0 |
= 0. |
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
x |
2 |
|
|
x |
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
+ |
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
x cos x = x 1− |
2! |
4! |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Будем искать решение уравнения в виде ряда y = c |
0 |
+ c x + c |
2 |
x2 |
+ c |
3 |
x3 |
+ c |
4 |
x4 |
+ . |
|
|
||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Тогда |
y′ = c1 + 2c2 x +3c3 x2 |
|
+ 4c4 x3 +5c5 x4 |
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
y′′ = 2c2 + 6c3 x +12c4 x2 + 20c5 x3 +30c6 x4 + |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
Из начальных условий находим c0 |
= 0, c1 |
=1. Подставим полученные ряды в исходное |
|||||||||||||||||||||||||||||||||||||||||||||||||||||
уравнение |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
(2c |
2 |
+ |
6c |
3 |
x +12c |
4 |
x2 |
+ 20c |
5 |
x |
3 +30c |
6 |
x4 + ) + x |
(c |
|
+ |
2c |
2 |
x |
+3c |
3 |
x2 |
+ 4c |
4 |
x3 |
+5c |
5 |
x4 + ) + |
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
2 |
|
x |
4 |
|
|
x |
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
+ (c0 |
+ c1 x + c2 x |
2 |
+ c3 x |
3 |
+ c4 x |
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
+ ) |
= x 1− |
2! |
+ |
|
4! |
|
− |
6! |
|
+ . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Приравниваем коэффициенты при одинаковых степенях x:
x0 |
|
2c2 = 0 |
|
||
x1 |
|
6c3 + 2 =1 |
x2 |
|
12c4 +3c2 = 0 |
x3 |
|
20c5 + 4c3 = − 12 |
……………………. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
1 |
|
|
|
|
1 |
|
|
Решая эту систему, находим: c |
2 |
= c |
4 |
= c |
6 |
= = c |
2n |
= 0 , |
c |
3 |
= − |
, c |
5 |
= |
, |
c |
7 |
= − |
|
, |
||||||||
|
5! |
7! |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
3! |
|
|
|
|
|
||||||||||||
Получаем искомое решение в виде ряда |
|
y = x − |
x3 |
+ |
x5 |
− |
x7 |
+ или |
|
y = sin x . |
|
|
|
|||||||||||||||
|
3! |
5! |
|
7! |
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|