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

в точке 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 +eNx

0

0

5

20

yx3 +eNx

0

1

5

106

Продолжение табл.

Вариант

f (x, y)

x0

y0

xk

 

 

 

 

 

21

y x + eNx

0

1

5

22

2yx +eNy

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 +ey

0

1

0,5

30

cos(xy / N)

0

0

0,5

 

 

 

 

 

107

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