Дифференциальных
Скачать 157.37 Kb.
|
3 РЕШЕНИЕ КРАЕВЫХ ЗАДАЧ МЕТОДОМ КОНЕЧНЫХ РАЗНОСТЕЙНаиболее распространенным и универсальным численным методом решения дифференциальных уравнений является метод конечных разно- стей. Основное содержание метода заключается в следующем. Область непрерывного изменения аргумента (например, отрезок) заменяется дис- кретным множеством точек, называемых узлами. Эти узлы составляют разностную сетку. Искомая функция непрерывного аргумента прибли- женно заменяется функцией дискретного аргумента на заданной сетке. Эта функция называется сеточной. Исходное дифференциальное уравне- ние заменяется разностным уравнением относительно сеточной функции. При этом для входящих в уравнение производных используются соответ- ствующие конечно-разностные соотношения. Такая замена дифференци- ального уравнения разностным называется его аппроксимацией на сетке (или разностной аппроксимацией). Решение дифференциального уравнения сводится к отысканию зна- чений сеточной функции в узлах сетки. Обоснованность замены диффе- ренциального уравнения разностным, точность получаемых решений, устойчивость метода – важнейшие вопросы, которые требуют тщатель- ного изучения. Рассмотрим содержание метода на примере решения дифференци- ального уравнения второго порядка y fx, y, y при заданных граничных условиях (3.1)y0 y0 , y1 y1 . (3.2) Разобьем отрезок [0, 1] на n равных частей точками xi= ih, (i = 0, 1, …, n). Решение краевой задачи (3.1), (3.2) сведем к вычислению значений сеточной функции yiв узловых точках xi. Для этого воспользу- емся уравнением (3.2) для внутренних узлов: y xi fxi, yxi, y xi, i 1, 2,..., n- 1. (3.3)Заменим производные, входящие в эти соотношения, их конечно- разностными аппроксимациями: yx yi1 yi1 Oh2 , i2h |