Разностные аппроксимации - (реферат)
Разностные аппроксимации - (реферат)
Дата добавления: март 2006г.
Разностные аппроксимации
1. Примеры разностных аппроксимаций.
Различные способы приближенной замены одномерных дифференциальных уравнений разностными изучались ранее. Напомним примеры разностных аппроксимаций и введем необходимые обозначения. Будем рассматривать равномерную сетку с шагомh, т. е. множество точек wh={xi=ih, i=0, ±1, ±2, …}.
Пусть u(x) – достаточно гладкая функция, заданная на отрезке [xi-1, xi+1]. Обозначим
Разностные отношения
называются соответственно правой, левой и центральной разностными производными функцииu(x) в точке xi , т. е. при фиксированном xi и при h®0 (тем самым при i®Ґ) пределом этих отношений является u’(xi). Проводя разложение по формуле Тейлора, получим
ux, i – u’(xi) = 0, 5hu’’(xi) + O(h2),
ux, i – u’(xi) = -0, 5hu’’(xi) + O(h2),
ux, i – u’(xi) = O(h2),
Отсюда видно, что левая и правая разностные производные аппроксимируют u’(x) с первым порядком по h, а центральная разностная производная – со вторым порядком. Нетрудно показать, что вторая разностная производная
аппроксимирует u’’(xi) со вторым порядком по h, причем справедливо разложение
Рассмотрим дифференциальное выражение
(1)
с переменным коэффициентом k(x). Заменим выражение (1) разностным отношением
(2)
где a=a(x) – функция, определенная на сетке wh. Найдем условия, которым должна удовлетворять функция a(x) для того, чтобы отношение (aux)x, i аппроксимировало (ku’)’ в точке xi со вторым порядком по h. Подставляя в (2) разложения
где ui’ = u’(xi), получим
С другой стороны, Lu = (ku’)’ = ku’’ + k’u’,
т. е.
Отсюда видно, что Lhu–Lu = O(h2), если выполнены условия
(3)
Условия (3) называются достаточными условиями второго порядка аппроксимации. При их выводе предполагалось, что функция u(x) имеет непрерывную четвертую производную и k(x)–дифференцируемая функция. Нетрудно показать, что условиям (3) удовлетворяют, например, следующие функции:
Заметим, что если положить ai = k(xi), то получим только первый порядок аппроксимации. В качестве следующего примера рассмотрим разностную аппроксимацию оператора Лапласа
(4)
Введем на плоскости (x1, x2) прямоугольную сетку с шагом h1 по направлению x1 и с шагом h2 по направлению x2, т. е. множество точек
wh = xi1 = ih1, xj2 = jh2; i, j = 0, ±1, ±2, …,
и обозначим
Из предыдущих рассуждений следует, что разностное выражение
(5)
аппроксимирует дифференциальное выражение (4) со вторым порядком, т. е. Lhuij – Lu(xi1, xj2) = O(h21) + O(h22). Более того, для функций u(x1, x2), обладающих непрерывными шестыми производными, справедливо разложение
Разностное выражение (5) называется пятиточечным разностным оператором Лапласа, так как оно содержит значения функции u(x1, x2) в пяти точках сетки, а именно в точках (x1i, x2j), (x1i±1, x2j), (x1i, x2 j±1). Указанное множество точек называется шаблоном разностного оператора. Возможны разностные аппроксимации оператора Лапласа и на шаблонах, содержащих большее число точек.
2. Исследование аппроксимации и сходимости
2. 1. Аппроксимация дифференциального уравнения. Ранее рассматривалась краевая задача (k(x) u’(x))’ – q(x) u(x) + f(x) = 0, 0 < x < l, (1)
– k(0) u’(0) + bu(0) = m1, u(l) = m2, (2)
k(x) і c1 > 0, b і 0,
для которой интегро-интерполяционным методом была построена разностная схема
(3)
(4)
где
(5)
(6)
Обозначим через Lu(x) левую часть уравнения (1) и через Lhyi – левую часть уравнения (3), т. е.
Пусть u(x) – достаточно гладкая функция и u(xi) – ее значение в точке xi сетки wh = {xi = ih, i = 0, 1, …, N, hN = l} (7)
Говорят, что разностный оператор Lh аппроксимирует дифференциальный оператор L в точке x=xi, если разность Lhui – Lhu(xi) стремится к нулю при h®0. В этом случае говорят также, что разностное уравнение (3) аппроксимирует дифференциальное уравнение (1).
Чтобы установить наличие аппроксимации, достаточно разложить по формуле Тейлора в точкеx=xi значения ui±1 = u(xi ± h), входящие в разностное выражение Lhui. Большая часть этой работы проделана в предыдущей главе, где показано, что при условиях
(8)
выполняется соотношение
Если кроме того, докажем, что
di = q(xi) + O(h2), ji = f(xi) + O(h2) (9)
то тем самым будет установлено, что оператор Lh аппроксимирует L со вторым порядком по h, т. е.
Lhui – Lu(xi) = O(h2), i = 1, 2, …, N–1 (10)
Итак, доказательство второго порядка аппроксимации сводится к проверке сводится к проверке условий (8), (9) для коэффициентов (5), (6). Проверим сначала выполнение условий (8). Обозначаяp(x) = k-1(x), получим
следовательно,
Аналогично
Отсюда получим
т. е. условия (8) выполнены. Условия (9) выполнены в силу того, что замена интегралов (6) значениямиqi, fiсоответствует приближенному вычислению этих интегралов по формуле прямоугольников с узлом в середине отрезка интегрирования.
2. 2. Аппроксимация граничного условия. Исследуем погрешность аппроксимации разностного граничного условия (4). Обозначимlhu(0) = –a1ux, 0 + bu0. Если u(x) – произвольная достаточно гладкая функция, то очевидно lhu(0) = –k(0) u’(0) + bu(0) + O(h),
т. е. имеет место аппроксимация первого порядка по h. Однако если u=u(x) –решение задачи (1), (2), то разностное граничное условие (4) имеет второй порядок аппроксимации, т. е.
Докажем последнее утверждение. Используя разложение
ux, 0 = (u1 – u0)/h = u’(x1/2) + O(h2), x1/2 = 0, 5h,
a1 = k1/2 + O(h2)
получим
Отсюда имеем
Учитывая граничное условие (2), получаем
lhu(0) = 0, 5h [– (ku’)’(0) + d0u0 – j0] + O(h2).
Выражение, стоящее в квадратных скобках, преобразуем, учитывая уравнение (1), к виду
– (ku’)’(0) + d0u0 – j0 = – (ku’)’(0) + q(0)u(0) – f(0) +
+ (d0 – q(0))u0 – (f(0) – j0) = (d0 – q(0))u0 – (f(0) – j0).
Из соотношений
получаем
что и требовалось доказать.
Таким образом, при достаточной гладкости коэффициентов k(x), q(x), f(x) и решения u(x)разностная схема (10) аппроксимирует исходную задачу (2) со вторым порядком поh. При практическом использовании разностной схемы для нахождения ее коэффициентов не обязательно вычислять интегралы (4), (6) точно. Можно воспользоваться коэффициентами, полученными путем замены этих интегралов квадратурными формулами, имеющими точностьO(h2)и выше. Например, в результате применения формулы прямоугольников получим следующие коэффициенты: ai = k(xi – 0, 5h), di = q(xi), ji = f(xi). Применяя формулу трапеций, получим
Представление коэффициентов разностной схемы в виде интегралов (4), (6) оказывается полезным при исследовании сходимости в случае разрывных функцийk(x), q(x), f(x).
2. 3. Уравнение для погрешности. Решение yi = y(xi) разностной задачи (3), (4) зависит от шага h сетки, y(xi) = yh(xi). По существу, мы имеем семейство решений {yh(xi)}, зависящее от параметра h. Говорят, что решение yh(x) разностной задачи сходится к решению u(x) исходной дифференциальной задачи, если при h®0 погрешность yh(xi) – u(xi), i = 0, 1, …, N, стремится к нулю в некоторой норме. В настоящем параграфе в качестве такой нормы будем брать норму в сеточном пространствеC(wh), т. е. положим
Говорят, что разностная схема имеет m-й порядок точности (или сходится с порядком m), если
где m>0, M>0 – константы, не зависящие от h.
Выше было установлено, что схема (3), (4) имеет второй порядок аппроксимации. Докажем теперь, что эта схема имеет и второй порядок точности. Для этого прежде всего выпишем уравнение, которому удовлетворяет погрешностьzi = yi – u(xi). Поставим yi = zi + u(xi) в уравнения (3), (4). Тогда получим уравнения
(11)
(12)
где обозначено
Функция yi, входящая в правую часть уравнения (11), называется погрешностью аппроксимации дифференциального уравнения (1) разностным уравнением (3) на решении задачи (1), (2). В п. 1 было доказано, чтоyi = O(h2) при h®0, i=1, 2, …, N–1. Аналогично, величина n1является по определению погрешностью аппроксимации краевого условия (2) разностным краевым условием (4) на решении задачи (1), (2), причемn1=O(h2). Таким образом, структура уравнений для погрешности (11), (12) та же, что и у разностной схемы (3), (4), отличаются только правые части.
Чтобы доказать сходимость разностной схемы, оценим решение задачи (11), (12) через правые частиyi, n1, т. е. получим неравенство вида
(13)
с константой M1, не зависящей от h. Из этого неравенства и будет следовать, что
Отметим, что неравенства вида (13), называемые априорными оценками, нашли широкое применение в теории разностных схем. Поскольку структура для погрешности (11), (12) та же, что и у разностной схемы (3), (4), а отличаются только правые части, то оценка (13) выполняется одновременно с аналогичной оценкой
