Xreferat.com » Рефераты по математике » РЕШЕНИЕ СИСТЕМ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ПЯТИТОЧЕЧНЫМ МЕТОДОМ АДАМСА – БАШФОРТА

РЕШЕНИЕ СИСТЕМ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ПЯТИТОЧЕЧНЫМ МЕТОДОМ АДАМСА – БАШФОРТА

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ УКРАИНЫ

НАЦИОНАЛЬНЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ

“ ХАРЬКОВСКИЙ ПОЛИТЕХНИЧЕСКИЙ ИНСТИТУТ”
Кафедра “Системы и Процессы Управления”

ОТЧЕТ

о научно-исследовательской курсовой работе

по численным методам

на тему :

« РЕШЕНИЕ СИСТЕМ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ПЯТИТОЧЕЧНЫМ МЕТОДОМ АДАМСА – БАШФОРТА »


Выполнил студент

гр.И-29 Уханов Е.В.

Руководитель работы

Д.т.н. проф Бреславский Д.В.


Харьков 2001

СОДЕРЖАНИЕ


Введение………………………………………………………………………..3

  1. Постановка задачи …………………………………………………………4

  2. Методы решения………………..…………………………………………6

2.1. Метод прогноза и коррекции …………………………………………6

2.2 Модифицированный метод Гаусса ………………………………….12

3. Описание алгоритма ………………………………………………………14

4. Описание программы ……………………………………………………..15

5. Примеры расчетов ………………………………………………………...17

5.1. Решение одного дифференциального уравнения …………………...17

5.2. Решение системы дифференциальных уравнений ………………….19

Заключение ……………………………………………………………………20

Список использованной литературы ………………………………………..21

Приложение 1 …………………………………………………………………22

Приложение 2 …………………………………………………………………23

Приложение 3 …………………………………………………………………24

Приложение 4 …………………………………………………………………25


ВВЕДЕНИЕ


Во многих областях науки и техники , а также отраслях наукоемкой промышленности , таких как : авиационная , космическая , химическая , энергетическая , - являются весьма распространенные задачи прогноза протекания процессов , с дальнейшей их коррекцией .

Решение такого рода задач связано с необходимостью использования численных методов , таких как : метод прогноза и коррекции , метод Адамса-Башфорта , метод Эйлера , метод Рунге-Кута , и др. При этом , стоит задача решения системы линейных дифференциальных уравнений первого порядка одним из методов интегрирования , на произвольном промежутке времени . Одним из оптимальных методов дающих высокую точность результатов – является пяти точечный метод прогноза и коррекции Адамса-Башфорта . Для повышения точности метода используется трех точечный метод прогноза и коррекции с автоматическим выбором шага , что приводит к универсальному методу интегрирования систем дифференциальных уравнений произвольного вида на любом промежутке интегрирования .

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

Целью данной курсовой работы является разработка алгоритма решения систем линейных дифференциальных уравнений первого порядка пяти точечным методом прогноза и коррекции Адамса-Башфорта .


1. ПОСТАНОВКА ЗАДАЧИ


Рассмотрим произвольную систему линейных дифференциальных уравнений первого порядка :


(1.1)


тогда как :


А = (1.2)


где А заданная матрица размером N x N .

- вектор с N координатами , который подлежит определению ;

N – произвольное целое число ;

  • заданные вектора правых частей с N координатами .

С использованием метода прогноза и коррекции Адамса-Башфорта пятого порядка , необходимо получить значения неизвестных для заданных временных интервалов . Для стартования метода необходимо использовать метод прогноза и коррекции третьего порядка с переменным шагом , на заданных временных промежутках ..


2. МЕТОДЫ РЕШЕНИЯ


2.1. Метод прогноза и коррекции


Метод прогноза и коррекции относится к задачам класса Коши , а именно к численным решениям многошаговыми методами .

Рассмотрим задачу Коши :

, (2.1.1)

Подставим в (2.1.1) точное решение y(x) , и проинтегрируем это уравнение на отрезке , тогда получим :

(2.1.2)


где в последнем член предполагаем , что p(x) полином , аппроксимирующий f(x,y(x)) . Чтобы построить этот полином , предположим , что - приближения к решению в точках . Будем считать для начала , что узлы Xi расположены равномерно с шагом h . тогда fi = f(xi,yi), ( i=k,k-1,k-2,…,k-N) есть приближения к f (x,y(x)) в точках и мы в качестве P возьмем интерполяционный полином для выбора данных (xi,fi) ,

( i =k,k-1,k-2,…,k-N) . Таким образом , P – полином степени N , удовлетворяющий условиям P(xi)=fi , ( i = k,k-1,k-2,…,k-N) . В принципе , можем проинтегрировать этот полином явно , что ведет к следующему методу :

(2.1.3)

В простейшем случае , когда N=0 , полином P есть константа , равная fk , и (2.1.3) превращается в обычный метод Эйлера :

(2.1.4)

Если N=1 , то P есть линейная функция , проходящая через точки

(xk-1,fk-1) и (xk,fk) , т.е.

(2.1.5)

интегрируя этот полином от Xk до Xk+1 , получим следующий метод :

(2.1.6)

который является двухшаговым , поскольку использует информацию в двух точках xk и xk-1 . Аналогично , если N=2 , то P - есть кубический интерполяционный полином , а соответствующий метод определяется формулой :

(2.1.7)

Отметим , что метод (2.1.6) – есть метод Адамса-Башфорта второго порядка , (2.1.7) – метод Адамса-Башфорта четвертого порядка .

Для стартования метода (2.1.7) необходимы сведения о четырех предыдущих точках . Соответственно данный метод требует вычисления стартующих данных . Воспользуемся для нахождения второй точки одношаговым методом Эйлера , который имеет вид :

Таким образом , подставляя начальные условия, мы находим вторую точку . Следует заметить , что степень точности совпадает со степенью точности остальных методов , что является существенным фактором в стартовании метода прогноза и коррекции .

Ввиду того , что стартовые методы имеют более низкий порядок , в начале приходится считать с меньшим шагом и с использованием большего промежутка времени . В данном случае метод Эйлера для дальнейшего интегрирования не оправдывает себя . Для этих целей воспользуемся трехшаговым методом прогноза и коррекции с переменным шагом .

Рассуждая также , как для метода Адамса-Башфорта , который излагается в работах : [1],[2],[3] , мы мы приходим к формулам :

Прогноз :

(2.1.8)

Коррекция :

Если Вам нужна помощь с академической работой (курсовая, контрольная, диплом, реферат и т.д.), обратитесь к нашим специалистам. Более 90000 специалистов готовы Вам помочь.
Бесплатные корректировки и доработки. Бесплатная оценка стоимости работы.

Поможем написать работу на аналогичную тему

Получить выполненную работу или консультацию специалиста по вашему учебному проекту
Нужна помощь в написании работы?
Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Пишем статьи РИНЦ, ВАК, Scopus. Помогаем в публикации. Правки вносим бесплатно.

Похожие рефераты: