Главная страница

математическое моделирование. 7915726_мат_мод_испр. Используя заданную функцию f(x), которая выбирается из таблицы 1 по числу N


Скачать 191.88 Kb.
НазваниеИспользуя заданную функцию f(x), которая выбирается из таблицы 1 по числу N
Анкорматематическое моделирование
Дата29.01.2022
Размер191.88 Kb.
Формат файлаdocx
Имя файла7915726_мат_мод_испр.docx
ТипДокументы
#345362
страница5 из 5
1   2   3   4   5

Задание 7.


Решить дифференциальное уравнение у' = f(х) + ху при за­данных начальных условиях хо = а, у(хо)= у(а) = 0 в заданных пределах [a, b] с шагом не менее (b - а)/ 10. Метод численного решения дифференциального уравнения выбрать по числу N3+1 из следующего списка методов:

1. Методы Эйлера

у' = ln(x+1)*sin(x) + х*у

хо = 0, у(0) = 0

h = 1
Решение

yi+1 = yi + h·f(xi, yi)

Приведенное выше соотношение является расчётной формулой метода Эйлера численного решения задачи Коши. Вычислив yi , i = 0,1,..,n получим таблицу значений решения в точках xi , i = 0,1,..,n

Погрешность формулы равна o(h2). К сожалению, метод Эйлера накапливает ошибку от шага к шагу. 

Первый шаг

x0 = 0; y0 = 0

y1 = 0+ ln(0+1)*sin(0) + 0*0 = 0

Первый шаг

x1 = 1; y1 = 0

y2 = 0+ ln(1+1)*sin(1) + 1*0 = 0,583

Третий шаг

x2 = 2; y2 = 0,583

y3 = 3+ ln(2+1)*sin(2) + 2*0,583 = 2,165

Остальные шаги сведем в таблицу

i

xi

yi

f(xi,yi)

hf(xi,yi)

1

0

0

0

0

2

1

0,000

0,583

0,583

3

2

0,583

2,165

2,165

4

3

2,749

8,442

8,442

5

4

11,191

43,545

43,545

6

5

54,735

271,958

271,958

7

6

326,693

1959,616

1959,616

8

7

2,286E+03

1,601E+04

1,601E+04

9

8

1,829E+04

1,463E+05

1,463E+05

10

9

1,646E+05

1,482E+06

1,482E+06

11

10

1,646E+06

1,646E+07

1,646E+07


Модифицированный метод Эйлера

y*i+1 = yi+h*f(xi,yi)

yi+1 = yi+h/2*(f(xi,yi) + f(xi+1,y*i+1))

Расчеты представим в таблице

h =

1

 

 

 

i

xi

yi

y*i+1

h/2*(f(xi;yi)+f(xi+1,y*i+1))

1

0

0,000

0,000

0,292

2

1

0,292

1,167

2,103

3

2

2,395

8,184

15,269

4

3

17,664

70,851

167,686

5

4

185,350

925,533

2683,065

6

5

2868,415

17208,771

58796,218

7

6

6,166E+04

4,317E+05

1,696E+06

8

7

1,757E+06

1,406E+07

6,239E+07

9

8

6,415E+07

5,773E+08

2,855E+09

10

9

2,919E+09

2,919E+10

1,591E+11

11

10

1,620E+11

1,782E+12

1,061E+13


В заключение приведем решение методом Рунге-Кутты, реализованное в MathCad



Как видим расхождение между наиболее точным методом Рунге-Кутты и методом Эйлера весьма существенны, особенно в конце интервала.
1   2   3   4   5


написать администратору сайта