- •Глава 1. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ. ПОНЯТИЙНЫЙ АППАРАТ
- •1. МОДЕЛЬ КАК НЕОПРЕДЕЛЯЕМОЕ ПОНЯТИЕ
- •2. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ: ХАРАКТЕРИСТИКИ И ЭТАПЫ МОДЕЛИРОВАНИЯ
- •3. МОДЕЛИРОВАНИЕ СИСТЕМ
- •4. КЛАССИФИКАЦИЯ СИСТЕМ. СИСТЕМНЫЙ ПОДХОД
- •5. ТРЕБОВАНИЯ К МОДЕЛЯМ
- •Глава 2. УРАВНЕНИЯ ПЕРВОГО ПОРЯДКА
- •1. ОСНОВНЫЕ ПОНЯТИЯ
- •2. УРАВНЕНИЯ С РАЗДЕЛЯЮЩИМИСЯ ПЕРЕМЕННЫМИ
- •3. ОДНОРОДНОЕ УРАВНЕНИЕ ПЕРВОГО ПОРЯДКА
- •4. ЛИНЕЙНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ПЕРВОГО ПОРЯДКА
- •5. ТИПОВЫЕ ЗАДАНИЯ ПО ПОСТРОЕНИЮ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ СРЕДСТВАМИ УРАВНЕНИЙ ПЕРВОГО ПОРЯДКА
- •Глава 3. УРАВНЕНИЯ ВТОРОГО И ВЫСШИХ ПОРЯДКОВ
- •1. УРАВНЕНИЯ, ДОПУСКАЮЩИЕ ПОНИЖЕНИЕ ПОРЯДКА
- •3. ЛИНЕЙНЫЕ ОДНОРОДНЫЕ УРАВНЕНИЯ ВТОРОГО ПОРЯДКА
- •6. ТИПОВЫЕ ЗАДАНИЯ ПО РЕШЕНИЮ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ВТОРОГО ПОРЯДКА
- •Глава 4. КРАЕВЫЕ ЗАДАЧИ И СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
- •1. КРАЕВЫЕ ЗАДАЧИ
- •2. СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
- •3. ТИПОВЫЕ ЗАДАНИЯ ПО РЕШЕНИЮ СИСТЕМ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
- •Глава 5. ПРИБЛИЖЁННОЕ РЕШЕНИЕ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
- •1. РАЗЛОЖЕНИЕ РЕШЕНИЯ В СТЕПЕННОЙ РЯД
- •2. МЕТОД ЭЙЛЕРА ПРИБЛИЖЁННОГО РЕШЕНИЯ ЗАДАЧИ КОШИ ДЛЯ ОБЫКНОВЕННОГО ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ ПЕРВОГО ПОРЯДКА
- •3. ТИПОВЫЕ ЗАДАНИЯ ПО ПРИБЛИЖЁННОМУ РЕШЕНИЮ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
- •ЗАКЛЮЧЕНИЕ
в точке xk = 0,5 (вычисления производить с точностью до 0.001). Разбить отрезок [x0 ; xk ] на n = 5 частей.
Решение. Так как x0 = 0, xk = 0,5, |
n = 5 , то h = |
xk − x0 |
= |
0,5 −0 |
= 0,1. |
|
|
5 |
|||||
|
|
n |
|
|
|
|
Из условия задачи y0 =1 , f (x, y) = x2 + y2 . |
|
|
|
|
||
Далее находим x1 = x0 + h = 0,1; |
x2 = x1 + h = 0,2; |
x3 = x2 + h = 0,3; |
||||
x4 = x3 + h = 0,4; x5 = x4 + h = 0,5 . |
|
|
|
|
|
|
Затем находим значения приближённого решения в узлах:
1)y1 = y0 + hf (x0 , y0 ) =1 + 0,1 (02 +12 ) =1,1;
2)y2 = y1 + hf (x1, y1) =1,1 + 0,1 (0,12 +1,12 ) =1,222;
3)y3 = y2 + hf (x2 , y2 ) =1,222 + 0,1 (0,22 +1,2222 ) =1,375;
4)y4 = y3 + hf (x3, y3 ) =1,375 + 0,1 (0,32 +1,3752 ) =1,573;
5)y5 = y4 + hf (x4 , y4 ) =1,573 + 0,1 (0,42 +1,5732 ) =1,837.
Таким образом, получили, что приближённое решение y(xk ) =1,837.
2.3.Погрешность приближённого решения в точке xk равна разности между точным значением искомой функции y(xk ) и значением сеточной функции yk . Чем больше число n отрезков, на которые разбивается исходный отрезок [x0; xk ] , тем более точным является полученное по методу Эйлера решение.
3.ТИПОВЫЕ ЗАДАНИЯ ПО ПРИБЛИЖЁННОМУ РЕШЕНИЮ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
3.1.Типовое задание с решением. Используя метод Эйлера, определить приближённым образом состояние процесса, смоделированного
ввиде задачи Коши
y′ = cos(xy) + y / 5y(0) =1
в момент x = 0,5 . Отрезок [x0 ; x] разбить на n1 = 5 равных частей. Вычисления проводить с точностью до 0,001.
Решение. Имеем x0 = 0, xk = 0,5, n = 5 , шаг |
разбиения, очевидно, |
h = 0,1. Согласно начальному условию y0 =1 , и |
f (x, y) = cos(xy) + y / 5 . |
Теперь находим значения приближённого решения в узлах:
103
1) x1 = x0 + h = 0,1;
y1 = y0 + hf (x0 , y0 ) =1+ 0,1 (cos(0 1) +1/ 5) =1,12;
2) x2 = x1 + h = 0,2;
y2 = y1 + hf (x1, y1) =1,12 + 0,1 (cos(0,1 1,12) +1/1,12) =1,242;
3) x3 = x2 + h = 0,3;
y3 = y2 + hf (x2 , y2 ) =1,242 + 0,1 (cos(0,2 1,242) +1/1,242) =1,364;
4) x4 = x3 + h = 0,4;
y4 = y3 + hf (x3, y3 ) =1,364 + 0,1 (cos(0,3 1,364) +1/1,364) =1,483;
5) x5 = x4 + h = 0,5;
y5 = y4 + hf (x4 , y4 ) =1,483 +0,1 (cos(0,4 1,483) +1/1,483) =1,595.
Следовательно, приближённое решение найдено: у(х) = y5 =1,595.
3.2. Задания для самостоятельного решения.
Задание 1. Найти три первых отличных от нуля члена разложения по формуле Маклорена решения y = y(x) следующей задачи Коши
1.y′ =
y(0)
2.y′ =
y(0)
3.y′ =
y(0)
4. y′ =y(0)
5. y′ =y(0)
y′ =
6.
y(0)
y′ =
7.
y(0)
2x2 + y2 +6x −1
= 0
2e y +4y + x −2
= 0
ex +4xy +cos x −1
= 0
x2 +2ex + y +1 2
= 0
sin2 2x +9sin y +1
=0
x2 +9y +6e y +1
=0
7 cos x +cos y +5x
= 0
|
|
|
|
11+2y |
2 |
+ x |
2 |
+3x |
||||||||
8. |
y′ = |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= 0 |
|
|
|
|
|
|
|
|
|||
|
y(0) |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
x3 |
|
|
|
|
|
|
|
|
|
|
|
9. |
y′ = |
|
|
|
|
+2sin y +3 |
||||||||||
|
6 |
|
||||||||||||||
|
|
|
|
|
= 0 |
|
|
|
|
|
|
|
|
|||
|
y(0) |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
2sin 2x +5e |
y |
|
|||||||||
10. |
y′ = |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= 0 |
|
|
|
|
|
|
|
|
|||
|
y(0) |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
2e |
2 y |
− x |
2 |
+ 4 |
|
||||||
11. |
y′ = |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= 0 |
|
|
|
|
|
|
|
|
|||
|
y(0) |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
= |
2sin x + |
4y |
2 |
−2 |
|||||||||
12. |
y′ |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
=1 |
|
|
|
|
|
|
|
|
|||
|
y(0) |
|
|
|
|
|
|
|
|
|
||||||
13. |
|
′ |
= cos x +e |
x |
+ y |
|
||||||||||
y |
|
|
|
|||||||||||||
|
|
|
|
|
= 3 |
|
|
|
|
|
|
|
|
|||
|
y(0) |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
x |
|
|
|
2 |
|
|
|
2 |
|
|
14. |
y′ = |
|
|
+ x |
|
+ y |
|
|
|
|||||||
2 |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y(0) = 2 |
|
|
|
|
|
|
|
|
|||||||
104
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
15. |
y′ = sin x +cos x +9y |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= −1 |
|
|
|
|
|
|||||
|
y(0) |
|
|
|
|
|
|
||||||
|
|
|
|
2 |
+ |
6y |
2 |
|
|
|
|
||
16. |
y′ = x |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
=1 |
|
|
|
|
|
|
|
||||
|
y(0) |
|
|
|
|
|
|
|
|||||
|
|
7 cos x + |
5x + y |
2 |
|
||||||||
17. |
y′ = |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= 2 |
|
|
|
|
|
|
|
||||
|
y(0) |
|
|
|
|
|
|
|
|||||
|
|
1+cos x + y |
2 |
|
|
||||||||
18. |
y′ = |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
=1 |
|
|
|
|
|
|
|
||||
|
y(0) |
|
|
|
|
|
|
|
|||||
|
|
|
y2 |
|
|
|
|
|
|
|
|||
19. |
y′ = |
|
|
|
|
+2cos x +3 |
|
||||||
2 |
|
|
|
||||||||||
|
|
=1 |
|
|
|
|
|
|
|
||||
|
y(0) |
|
|
|
|
|
|
|
|||||
|
|
|
|
2 |
|
+2x + y |
|
|
|
||||
20. |
y′ = y |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= 2 |
|
|
|
|
|
|
|
||||
|
y(0) |
|
|
|
|
|
|
|
|||||
|
|
7e |
−x |
+5x +cos y |
|
||||||||
21. |
y′ = |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= 0 |
|
|
|
|
|
|
|
||||
|
y(0) |
|
|
|
|
|
|
|
|||||
|
|
3y |
2 |
+ x |
2 |
− x |
|
|
|||||
22. |
y′ = |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= 0 |
|
|
|
|
|
|
|
||||
|
y(0) |
|
|
|
|
|
|
|
|||||
|
|
x2 |
|
|
|
|
|
−x |
||
23. |
y′ = |
|
|
+7 y +e |
|
|
|
|||
2 |
|
|
|
|||||||
|
|
= 0 |
|
|
|
|
|
|||
|
y(0) |
|
|
|
|
|
||||
|
|
y3 |
|
|
|
|
|
|||
24. |
y′ = |
|
|
|
|
+6cos x |
|
|
||
3 |
|
|
|
|
|
|||||
|
|
=1 |
|
|
|
|
|
|||
|
y(0) |
|
|
|
|
|
||||
|
|
4sin x +9y |
2 |
+1 |
||||||
25. |
y′ = |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
= −1 |
|
|
|
|
||||
|
y(0) |
|
|
|
|
|||||
|
|
2x +9y +6e |
y |
|
||||||
26. |
y′ = |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
= 0 |
|
|
|
|
|
|||
|
y(0) |
|
|
|
|
|
||||
|
|
4xy −1+3e |
x |
|
|
|||||
27. |
y′ = |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
= 0 |
|
|
|
|
|
|||
|
y(0) |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
x |
28. |
y′ = cos x +4y +e |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
= 3 |
|
|
|
|
|
|||
|
y(0) |
|
|
|
|
|
||||
|
|
6e |
y |
+9y +e |
x |
|
||||
29. |
y′ = |
|
|
|
|
. |
||||
|
|
|
|
|
|
|
|
|
||
|
|
= 0 |
|
|
|
|
|
|||
|
y(0) |
|
|
|
|
|
||||
|
|
2 −cos x + y |
2 |
|
||||||
30. |
y′ = |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
= 2 |
|
|
|
|
||||
|
y(0) |
|
|
|
|
|||||
В следующих задачах для самостоятельного решения число N от 1 до 20 задаётся преподавателем.
Задание 2. Найти приближённое решение задачи Коши для обыкновенного дифференциального уравнения первого порядка
y′ = f (x, y)
y(x0 ) = y0
в точке x = xk с помощью метода Эйлера. Использовать разбиение отрезка [x0; x1] на n1 = 5 и n2 =10 равных частей.
105
Вариант |
|
f (x, y) |
x0 |
y0 |
xk |
|
|
|
|
|
|
1 |
sin(xy) + Ny2 |
0 |
1 |
0,5 |
|
2 |
sin(x + Ny) + y |
0 |
0 |
0,5 |
|
|
|
|
|
|
|
3 |
Ny2 +ex |
0 |
1 |
0,5 |
|
4 |
N + |
y +sin x |
0 |
4 |
0,5 |
|
|
|
|
|
|
5 |
N |
y +2x2 |
0 |
1 |
1 |
|
|
|
|
|
|
6 |
10 |
y + sin(x / N) |
0 |
1 |
1 |
|
|
|
|
|
|
7 |
cos y y + x2 y / N |
0 |
0.1 |
1 |
|
|
|
|
|
|
|
8 |
y / N + x2 |
0 |
0.1 |
1 |
|
|
|
|
|
|
|
9 |
3x + Ny +ln y |
0 |
1 |
0,5 |
|
|
|
|
|
|
|
10 |
x / N + 2yxln y |
1 |
1 |
1,5 |
|
|
|
|
|
|
|
11 |
xy / N +ln y |
1 |
1 |
1,5 |
|
|
|
|
|
|
|
12 |
3xy +(ln y) / N |
1 |
1 |
1,5 |
|
|
|
|
|
|
|
13 |
e2x y + N |
0 |
0 |
0,5 |
|
|
|
|
|
|
|
14 |
ex + xy / N |
0 |
1 |
1 |
|
|
|
|
|
|
|
15 |
ex y + Nx |
0 |
1 |
0,5 |
|
16 |
ex sin у + Nx |
0 |
1 |
0,5 |
|
|
|
|
|
|
|
17 |
ex cos y + xy / N |
0 |
1 |
1 |
|
|
|
|
|
|
|
18 |
ex (1+ y) + x / N |
0 |
1 |
1 |
|
19 |
yx2 +e−Nx |
0 |
0 |
5 |
|
20 |
yx3 +e−Nx |
0 |
1 |
5 |
|
106
Продолжение табл.
Вариант |
f (x, y) |
x0 |
y0 |
xk |
|
|
|
|
|
21 |
y x + e−Nx |
0 |
1 |
5 |
22 |
2yx +e−Ny |
0 |
0 |
1 |
23 |
5yx +2sin(x / N ) |
0 |
1 |
0,1 |
|
|
|
|
|
24 |
2x2 +cos y / N |
0 |
1 |
0,5 |
|
|
|
|
|
25 |
2ln x +cos( y / N ) |
1 |
1 |
1,5 |
|
|
|
|
|
26 |
x + y + xy / N |
0 |
1 |
0,5 |
|
|
|
|
|
27 |
xy + x + y / N |
0 |
1 |
1 |
|
|
|
|
|
28 |
x + N / y |
0 |
1 |
0,5 |
|
|
|
|
|
29 |
x / N +e−y |
0 |
1 |
0,5 |
30 |
cos(xy / N) |
0 |
0 |
0,5 |
|
|
|
|
|
107
