Студопедия.Орг Главная | Случайная страница | Контакты | Мы поможем в написании вашей работы!  
 

Аппроксимация производных



Для аппроксимации первой производной можно воспользоваться формулами:

– правая разностная производная,

– левая разностная производная,

– центральная разностная производная.

Существует множество способов аппроксимации производной, которые следуют из определения производной:

.

На основе формул для разностной аппроксимации первой производной можно построить разностную аппроксимацию второй производной:

(6.3)

Аналогично можно получить аппроксимации производных более высокого порядка.

Определение. Погрешностью аппроксимации -ой производной называется разность .

Для определения порядка аппроксимации используется разложение в ряд Тейлора.

Рассмотрим правую разностную аппроксимацию первой производной:

Т.е. правая разностная производная имеет первый по порядок аппроксимации. Аналогичные оценки можно сделать для левой разностной производной.

Центральная разностная производная имеет второй порядок аппроксимации.

Аппроксимация второй производной по формуле (6.3) также имеет второй порядок аппроксимации.

Для того чтобы аппроксимировать дифференциальное уравнение, необходимо в нем заменить все производные их аппроксимациями. Рассмотрим задачу (6.1), (6.2) и заменим в (6.1) производные:

.

В результате получим:

(6.4)

Порядок аппроксимации исходной задачи равен 2, т.к. вторая и первая производные аппроксимированы с порядком 2, а остальные – точно.

Итак, вместо дифференциальных уравнений (6.1), (6.2) получена система линейных уравнений вида

(6.5)

для определения в узлах сетки. Матрица данной системы имеет вид:

Данная матрица является трехдиагональной, т.е. все ненулевые элементы расположены на главной диагонали и двух прилегающих к ней диагоналях. Для решения таких СЛАУ имеется экономичный метод прогонки (см. с. 31-32). Решая полученную систему уравнений, мы получим решение исходной задачи.

Для краевой задачи (6.1), (6.2) имеем:

Коэффициенты СЛАУ определяются формулами:

Прямой этап метода прогонки:

, , i =2, 3, … n.

Обратный этап метода прогонки:

.

Условие устойчивости (условие диагонального преобладания) тогда имеет вид:

.

Пусть . Тогда и, следовательно, , , т.е. .

Тогда условие устойчивости имеет вид и, как можно видеть, справедливо всегда.

Решение

Дано: ,

, h= 0,2

.

В результате получим разностную схему:

Или

Оформим в виде таблицы

i a c b f u q(x) e(x) z(x) x
    4,1 2,8 0,32          
  1,2 4,1 2,8 0,30 0,675     3,8 0,2
  1,2 4,1 2,8 0,29 -0,72     3,6 0,4
  1,2 4,1 2,8 0,27 0,87     3,4 0,6
  1,2 4,1 2,8 0,26 -0,868     3,2 0,8
    -1,0              

Где =(2- h), =4+2 , =2+ h, =2 , i=2,….,n-1

=4+2 , =2+ h, =2 -(2- h)φ

=0, =-1

Оформим в виде таблицы

i a c b f u q(x) e(x) z(x) x
    4,1 2,8 0,24          
  1,2 4,1 2,8 0,22 0,675     2,8 1,2
  1,2 4,1 2,8 0,21 -0,72     2,6 1,4
  1,2 4,1 2,8 0,19 0,87     2,4 1,6
  1,2 4,1 2,8 0,18 -0,868     2,2 1,8
    -1,0              
                   

Решая систему, методом прогонки, находим ,….

=0

4,1 +2,8 =0,24

1,2 +4,1 +2,8 =0,22

1,2 +4,1 +2,8 =0,21

1,2 +4,1 +2,8 =0,19

1,2 +4,1 +2,8 =0,18

=1





Дата публикования: 2015-11-01; Прочитано: 336 | Нарушение авторского права страницы | Мы поможем в написании вашей работы!



studopedia.org - Студопедия.Орг - 2014-2024 год. Студопедия не является автором материалов, которые размещены. Но предоставляет возможность бесплатного использования (0.008 с)...