Бельшина Ю.Н., Вакуленко С.В., Воронова В.Б., Галишев М.А., Кондратьев С.А., Моторыгин Ю.Д., Толстых В.И., Шарапов С.В

ОСНОВЫ МЕТОДИКИ УСТАНОВЛЕНИЯ ОЧАГА ПОЖАРА ПРИ РАССЛЕДОВАНИИ ДЕЛ О ПОЖАРАХ

Противопожарной службы

МЧС России

Метод Рунге - Кутта

Метод Рунге - Кутта является одним из наиболее употребительных методов высокой точности. Метод Эйлера можно рассматривать как простейший вариант метода Рунге - Кутта.

Рассмотрим задачу Коши для дифференциального уравнения

y (t) = f (t, y (t))

с начальным условием y (t ₀ ) = y _0.

Как и в методе Эйлера, выберем шаг h = и построим сетку с системой узлов t_i = t ₀ + ih, i = 0, 1, …, n.

Обозначим через y_i приближенное значение искомого решения в точке t_i.

Приведем расчетные формулы метода Рунге - Кутта четвертого порядка точности:

y_{i+ 1} = y_i + h (k + 2 k + 2 k + k),

k = f (t_i, y_i),

k = f (t_i +, y_i + k), (6.17)

k = f (t_i +, y_i + k),

k = f (t_i + h, y_i + hk),

i = 0, 1, …, n.

Оценка погрешности. Оценка погрешности метода Рунге - Кутта затруднительна. Грубую оценку погрешности дает правило Рунге (см. раздел 6.2). Так как метод Рунге - Кутта имеет четвертый порядок точности, т. е. p = 4, то оценка погрешности (6.6) примет вид

R | y- y |. (6.18)

Используя правило Рунге, можно построить процедуру приближенного вычисления решения задачи Коши методом Рунге - Кутта четвертого порядка точности с заданной точностью. Нужно, начав вычисления с некоторого значения шага h, последовательно уменьшать это значение в два раза, каждый раз вычисляя приближенное значение y, i = 0, 1, …, n. Вычисления прекращаются тогда, когда будет выполнено условие:

R | y- y | <. (6.19)

Приближенным решением будут значения y, i = 0, 1, …, n.

Пример 6.4.

Методом Рунге - Кутта четвертого порядка точности найдем решение на отрезке [0, 1] следующей задачи Коши.

y (t) = 2ty, y (0) = 1. (6.20)

Возьмем шаг h = 0.1. Тогда n = = 10.

В соответствии с (6.17) расчетные формулы примут вид:

y_{i+ 1} = y_i + h (k + 2 k + 2 k + k),

k = 2 t_iy_i,

k = 2(t_i +)(y_i + k), (6.21)

k = 2(t_i +)(y_i + k),

k = 2(t_i + h)(y_i + hk),

i = 0, 1, …, 10.

Задача (6.20) имеет точное решение: y (t) = e, поэтому погрешность определяется как абсолютная величина разности между точными и приближенными значениями _i = | y (t_i) - y_i |.

Найденные по формулам (6.21) приближенные значения решения y_i и их погрешности _i представлены в таблице 6.5:

Таблица 6.5

ti	yi	i	ti	yi	i
0.1 0.2 0.3 0.4 0.5	1.01005 1.04081 1.09417 1.17351 1.28403	10^-9 410^-9 210^-8 610^-8 210^-7	0.6 0.7 0.8 0.9 1.0	1.43333 1.63232 1.89648 2.24790 2.71827	510^-7 210^-6 310^-6 610^-6 210^-5

Список литературы 1. Амосов А., Дубинский Ю. А., Копченова Н.В. Вычислительные методы для инженеров: Учеб. пособие. - М.: Высш. шк., 1994. 2. Бахвалов Н. С. Численные методы. - М.: Наука, 1973. 3. Волков Е. А. Численные методы. - М.: Наука, 1987. 4. Дьяконов В. П. Математическая система Maple V R3/R4/R5. - М.: Изд-во "СОЛОН", 1998. 5. Калиткин Н.Н. Численные методы. - М.: Наука, 1978. 6. Копченова Н.В., Марон И. А. Вычислительная математика в примерах и задачах. - М.: Наука, 1972.

7. Пирумов У.Г. Численные методы.: Учебное пособие. - М.: Изд-во МАИ, 199

Санкт-Петербургский институт Государственной

Санкт-Петербург

Основы методики установления очага пожара при расследовании дел о пожарах: Учебное пособие / Под общей редакцией В.С. Артамонова.- СПб.: Санкт-Петербургский институт Государственной пПротивопожарной сСлужбы МЧС России, 2005. - 132 с.

Авторы: Ю.Н. Бельшина, С.В. Вакуленко, В.Б. Воронова, М.А. Галишев, С.А. Кондратьев, Ю.Д. Моторыгин, В.И. Толстых, С.В. Шарапов.

Учебное пособие «Основы методики установления очага пожара при расследовании дел о пожарах» предназначено для курсантов и слушателей высших учебных заведений пожарно-технического профиля при обучении инженерно-техническим аспектам мероприятий по расследованию пожаров. В работе кратко изложены основы методики осмотра места пожара и составления протокола осмотра; общие методические подходы к исследованию различных конструкционных материалов и веществ, при выявлении очага пожара.

Рецензенты:

Дата добавления: 2014-11-29