Глава II. Дискретные модели § 2.3. Аппроксимация граничных условий

Назад 2.3.2. Первая краевая задача для обыкновенного дифференциального уравнения второго порядкаВперед

Рассмотрим краевую задачу

uўўku = f(x),   x О W = [0,1],(1)

u(0) = a,   u(1) = b.(2)

Эту краевую задачу можно записать в "операторном" виде

Lu = f,(3)

lu = g,(4)

где L дифференциальный оператор, определяемый равенством L = d2/dx2kI, а l: F(W) ® R2 оператор "граничных условий": lu = (u(0), u(1)), g = (a, b). Последняя краевая задача может быть записана в "объединенной" операторной форме

Lu = F,(5)

где L = (L, l), а F = (f, g).

Допустим, мы аппроксимировали дифференциальный оператор L на равномерной сетке wh = {ih}ni=0следующим образом:

Lhuh = (uh)xx+kuh,   (uh О F(wh))

или, что то же,

(Lhuh)(xi) = uh(xi+1) – 2uh(xi) + uh(xi–1)
h2
 kuh(xi).

Тогда уравнение (1) аппроксимируется, например, уравнением

Lhuh = Lh f =def fh.

Граничные же условия переносятся на сеточные функции в данном случае напрямую (lh = l) и, таким образом, могут быть записаны в виде

lhuh = Lh g =def gh.

Таким образом, краевая задача (3)(4) переписывается в виде

Lhuh = fh,

lhuh = gh,

или, в "объединенной" форме,

Lhuh = Fh,

где Lh = (Lh, lh), а Fh = (fh, gh).