Реферат: Решение дифференциальных уравнений

Общие сведения. Требуется найти функцию, удовлетворяющую уравнению

(30)

и принимающую при x=x0заданное значение

(31)

При этом будем для определенности считать, что решение надо получить при .

Из курса дифференциальных уравнений известно, что решение задачи (30)-(31) существует, единственно и является гладкой функцией, если правая часть уравнения (30), являющаяся функцией двух переменных x, Y, удовлетворяет некоторым условиям гладкости. Будем считать, что эти условия выполнены и существует единственное гладкое решение .

Численное решение задачи Коши (30)-(31) состоит в том, чтобы получить искомое решение в виде таблицы его приближенных значений для заданных значений аргумента x на некотором отрезке [a,b]:

Для решения задачи Коши (30)-(31) будем использовать разностные методы. Введем последовательность точек и шаги. Точки называют узлами, а множество этих точек называют сеткой. В каждой точке вместо значений функции вводятся числа, аппроксимирующие точное решение на данном множестве точек. Функцию Y, заданную в виде таблицы, называют сеточной функцией.

Далее, заменяя значение производной в уравнении (30) отношением конечных разностей, осуществляем переход от дифференциальной задачи (30), (31) относительно функции y:

(32)

. (33)

Здесь разностное уравнение (32) записано в общем виде, а конкретное выражение его правой части зависит от способа аппроксимации производной. Для каждого численного метода получается свой вид уравнения (32).

Одношаговые методы. Простейшим численным методом решения задачи Коши для обыкновенного дифференциального уравнения является метод Эйлера. Он основан на разложении искомой функции Y(x) в ряд Тейлора в окрестностях узлов в котором отбрасываются все члены, содержащие производные второго и более высокого порядков. Запишем это разложение в виде

(34)

Заменим значения функции Y в узлах значениями сеточной функции yi. Кроме того, используя уравнение (30), полагаем

Будем считать для простоты узлы равноотстоящими, т.е.. Учитывая введенные обозначения и пренебрегая членами порядка, из равенства (34) получаем

(35)

Полагая i=0, с помощью соотношения (35) находим значение сеточной функции y1 при x=x1:

Требуемое здесь значение y0задано начальным условием (31), т.е. Аналогично могут быть найдены значения сеточной функции в других узлах:

Построенный алгоритм называется методом Эйлера. Разностная схема этого алгоритма представлена соотношениями (35). Они имеют вид рекуррентных формул, с помощью которых значение сеточной функции в узле вычисляется по её значению в предыдущем узле .

Рассмотрим вопрос о погрешности метода Эйлера. Погрешность eiв точке xi равна разности между значением сеточной функции yi и точным значением искомой функции Эта погрешность состоит из двух частей: определяется погрешностью начального значения Как правило, начальное значение задается точно, т.е. и тогда и, следовательно, равна нулю та часть погрешности решения, которая связана с Погрешность обусловлена отброшенными членами в разложении в ряд Тейлора (34). На каждом шаге эта погрешность имеет порядок, так как именно члены такого порядка отброшены в (34).

Пример. Найти решение задачи Коши методом Эйлера на отрезке [0; 0.4] с шагом h=0.1.

Решение. Здесь,,,.

Используя рекуррентные формулы

,

последовательно находим

при:

при:

при:

при:

Решение задачи Коши представлено в таблице

Существуют и другие явные одношаговые методы. Наиболее распространенным из них является метод Рунге–Кутта. На его основе могут быть построены разностные схемы разного порядка точности. Приведем схему Рунге–Кутта четвертого порядка. Вычисление приближенного значения в точке производится по следующим формулам

,

(36)

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

Все вычисления удобно производить по следующей схеме.

Порядок заполнения таблицы.

1. Записать в первой строке таблицы данные значения, .

2. Вычислить и записать в таблицу в качестве .

3. Во второй строке таблицы записать, .

4. Вычислить и записать в таблицу в качестве .

5. В третьей строке таблицы записать, .

6. Вычислить и записать в таблицу в качестве .

7. В четвертую строку таблицы записать, .

8. Вычислить и записать в таблицу в качестве .

9. В столбец записать числа,,,, просуммировать эти числа, полученную сумму разделить на 6 и записать в таблицу в качестве .

10. Вычислить. Затем все вычисления продолжаем в том же порядке, принимая за начальную точку. Результаты вычислений записать в следующий блок таблицы.

Метод Эйлера может рассматриваться как метод Рунге–Кутта первого порядка. Метод Рунге–Кутта (36) требует большего объема вычислений, однако окупается повышенной точностью, что дает возможность проводить счет с большим шагом. Другими словами, для получения результатов с одинаковой точностью в методе Эйлера потребуется значительно меньший шаг, чем в методе Рунге–Кутта.

Пример. Методом Рунге – Кутта найти решение задачи Коши на отрезке, приняв шаг h=0.1.

Решение. Покажем начало процесса.

Вычисление Последовательно имеем

Отсюда и, следовательно,

Аналогично вычисляются дальнейшие приближения. Результаты вычислений приведены в таблице 5.

Таблица 5

Решение задачи Коши представлено в таблице

Задание 6

Найти решение задачи Коши методом Эйлера на отрезке [а; b] с шагом h

Библиографический список

1 Омельченко В.П. Математика: учеб.пособие / В.П. Омелченко, Э.В. Курбатова.- Ростов н/Д: Феникс, 2008.- 380с.

2 Пантина И.В. вычислительная математика / И.В.Пантина, А.В.Синчуков.-М.: Маркет ДС, 2010.- 176с.

2 Амосов А.А., Вычислительные методы для инженеров / А.А.Амосов, Ю.А.Дубинский, Н.В. Копченова.- М.: Высш. шк.,1994.- 544 с.

3 Турчак Л.И. Основы численных методов / Л.И. Турчак. — М.: Наука, 1987. – 320 с.

4 Демидович Б.Н., Основы вычислительной математики / Б.Н.Демидович, И.А. Марон.- М.: Высш. шк., 1994. — 172 с.

5 Боглаев Ю.П. Вычислительная математика и программирование / Ю.П. Боглаев.- М.: Высш. шк., 1990. — 544 с.

6 Семенов М.П., Катрахова А.А. Жучкова В.В. Основы численных методов / М.П.Семенов, А.А.Катрахова, В.В.Жучкова.- Воронеж: Изд-во ВГТУ, 1997. — 62 с.

7 Ушаков Д.М.Введение в математические основы САПР / Д.М. Ушаков, В.П.Корячко, И.П.Норенков.-

8 Федорков Е.Д.Численные методы: учеб. пособие / Е. Д. Федорков, А. И. Бобров. — Воронеж: ВГТУ, 2004. — 164 с. — 33-00.

9 Федорков, Е.Д. Вычислительная математика: Учеб.пособие / Е. Д. Федорков, А. И. Бобров. — Воронеж: ГОУВПО «Воронежский государственный технический университет», 2006. — 168 с. — 35-00.

СОДЕРЖАНИЕ

МЕТОДИЧЕСКИЕ УКАЗАНИЯ

к выполнению контрольной работы

по дисциплине «Математическое обеспечение САПР»

для студентов специальности

230100.62 «Информатика и вычислительная техника», профиль «Системы автоматизированного проектирования»,

заочной формы обучения

Составитель

Пак Алла Анатольевна

Компьютерный набор А.А.Пак

Подписано к изданию. .2014.

Уч.-изд. л.. «С»

ФГБОУ ВПО «Воронежский государственный технический университет»

394026 Воронеж, Московский просп., 14