Решение краевой задачи для ОДУ методом конечных разностей

МИНИСТЕРСТВО НАУКИ И ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ

МАРИЙСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ

Кафедра ПиПЭВС

Курсовая работа

по дисциплине «Моделирование систем»

Решение краевой задачи для ОДУ методом конечных разностей.

Йошкар-Ола

2010

Теоретическая часть.

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

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

Конечные разности.

Исчисление конечных разностей связано с изучением свойств и применений разностей между соседними членами какой-нибудь последовательности или между значениями функции в точках, расположенных с постоянным интервалом в некотором пространстве. Слово «конечные» используется здесь в несколько устаревшем смысле «не бесконечно малые», т.е. не связанные с предельными переходами. Поскольку дифференциальное исчисление занимается изучением пределов разностей, а исчисление конечных разностей – самими разностями, то естественно, что между этими двумя теориями существуют много параллелей. Исчисления конечных разностей используются при интерполяции в математических таблицах, при суммировании числовых рядов, при вычислении интегралов и дифференцировании функций. Разности встречаются также в любой ситуации, когда надо описать поведение объекта, который испытывает воздействие меняющихся условий на определенном расстоянии (во времени и в пространстве). Например, термостату требуется значительное время, чтобы отреагировать на изменение температуры, поэтому он реагирует не на текущую температуру, а на ту, что была минуту назад. Другой пример: автомашиной управляет водитель, которому требуется какое-то время, чтобы отреагировать на возникшую на дороге ситуацию.

Под конечной разностью первого порядка функции f (x) принято понимать величину

где d – некоторая постоянная, которую часто, но не всегда, принимают равной 1. Разность второго порядка обозначается D2f и представляет собой разность разностей, т.е.

Продолжив этот процесс, мы получим разности более высоких порядков D3f (x), D4f (x), ј .

У истоков теории.

Хотя исследование свойств и использование конечных разностей приходится на современный период развития математики, Птолемей (ок. 150 н.э.) ввел в Альмагесте таблицу разностей первого порядка, чтобы облегчить расчеты в таблице длин хорд. Разности второго порядка использовал при вычислении своих таблиц логарифмов в 1624 Г.Бриггс. Теория интерполяции берет начало со знаменитой пятой леммы из 3-й книги Математических начал (1687) И.Ньютона, в которой впервые была приведена формула, носящая ныне его имя. Частный случай формулы Ньютона, открытый также независимо его современником Дж.Грегори (1638–1675), приведен ниже (см. формулу (7)). В общей формуле интерполяции Ньютона использовались разделенные разности, хотя этот термин, по-видимому, был введен О.де Морганом (1806–1871) в 1848. Первое применение исчисления конечных разностей к задачам теории вероятностей принято связывать с именами П.де Монтмора (1678–1719) и А.де Муавра (1667–1754).

Хотя Л.Эйлер (1707–1783) в своих работах по дифференциальному исчислению использовал предельные переходы в конечных разностях, основания современной теории конечных разностей были заложены в основном Ж.Лагранжем (1736–1813) и П.Лапласом (1749–1827). Первый из них ввел в исчисление конечных разностей символические методы, второй сделал конечные разности главным инструментом в своей Аналитической теории вероятностей (1812).

Под влиянием этих работ математики 19 в. принялись интенсивно разрабатывать предмет, и в 1860 Дж.Буль выпустил свой классический Трактат об исчислении конечных разностей. С тех пор это исчисление и круг его приложений существенно расширились. Одно из наиболее важных приложений конечные разности нашли в статистике. Особенно полезными они оказались в теории сериальной корреляции, в анализе случайных последовательностей и статистических временных рядов.

Алгоритм метода конечных разностей.

1. Разобьем отрезок [a,b] на n равных частей длины, или шага

Точки разбиения называются узлами, а их совокупность – сеткойна отрезке [a,b]

2. Вводим обозначения

3. Заменим производные односторонними конечно-разностными отношениями:

Эти формулы приближенно выражают значения производных во внутренних точках интервала [a,b].

Кроме того, краевые условия дополнительно дают еще два уравнения:

Составляем матрицу уравнений, которую можно решить любым численным методом.

Практическая часть.

Внимание, отключите Adblock

Вы посетили наш сайт со включенным блокировщиком рекламы!
Ссылка для скачивания станет доступной сразу после отключения Adblock!

Скачать