Решение краевой задачи для обыкновенного дифференциального уравнения с заданной точностью
Курсовой проект - Математика и статистика
Другие курсовые по предмету Математика и статистика
?ая погрешность меньше заданной точности, следовательно, проверка точности выполнена и полученная точность удовлетворяет заданной.
Таблица 1
XY(X)E
1) 1.30 2.200000 0.00e+000
2) 1.32 2.029192 4.53e-010
) 1.34 1.859824 8.27e-010
) 1.36 1.692046 1.13e-009
) 1.38 1.525999 1.36e-009
) 1.40 1.361814 1.53e-009
) 1.42 1.199613 1.65e-009
) 1.44 1.039511 1.72e-009
) 1.46 0.881613 1.75e-009
) 1.48 0.726016 1.74e-009
) 1.50 0.572810 1.70e-009
) 1.52 0.422075 1.63e-009
) 1.54 0.273886 1.54e-009
) 1.56 0.128307 1.42e-009
) 1.58 -0.014601 1.30e-009
) 1.60 -0.154787 1.17e-009
) 1.62 -0.292205 1.03e-009
) 1.64 -0.426817 8.85e-010
) 1.66 -0.558591 7.44e-010
) 1.68 -0.687500 6.07e-010
) 1.70 -0.813522 4.77e-010
22) 1.72 -0.936644 3.55e-010
) 1.74 -1.056855 2.45e-010
) 1.76 -1.174151 1.48e-010
) 1.78 -1.288531 6.56e-011
) 1.80 -1.400000 4.44e-016
3. Метод конечных разностей
.1 Описание метода
Решим дифференциальное уравнение (1.1) с краевыми условиями (1.2) методом конечных разностей. Для того чтобы получить систему конечно-разностных уравнений, разобьем отрезок [a,b] на n равных частей длины как показано на рисунке 1. Точки разбиения имеют абсциссы:
(i=0,1,2,…,n), .
Рис. 1
Значения в точках искомой функции y=y(x) и ее производных обозначим соответственно:
Введем также обозначения:
.
Заменяя производные симметричными конечно-разностными отношениями для внутренних точек отрезка [a, b] будем иметь:
(i=1,2,…,n-1) (3.1)
где - погрешность формулы порядка .
Подставим конечно-разностные выражения (3.1) в исходное дифференциальное уравнение (1.1). В результате подстановки получим:
(3.2)
Или запишем выражение (3.2) в виде:
,(i=1,2,...,n-1), (3.3)
Так как всего n+1 неизвестных на отрезке [a, b], а уравнений n-1, то не хватает еще двух уравнений для нахождения всех (i=0,1,…,n). Эти недостающие уравнения дают краевые условия (1.2). По условию курсовой работы , , поэтому система (1.2) принимает вид:
С учетом этого запишем систему линейных алгебраических уравнений:
,
где (i=1,2,...,n-1) (3.4)
Таким образом, получена линейная система (3.4) из n+1 уравнений с n+1 неизвестными , представляющими собой значения искомой функции y=y(x) в точках .
Систему уравнений (3.4) можно записать в матричном виде, при этом получается трехдиагональная матрица коэффициентов
Вычислить решение системы можно методом алгебраической прогонки. Метод прогонки состоит из двух частей: прямой и обратный ход.
Приведем матрицу (3.5) к двухдиагональному виду. Предположим, что из (3.4) исключена неизвестная . Тогда это уравнение примет вид:
(3.7)
где (i=1,2,...,n-1) некоторые коэффициенты.
В прямом ходе метода прогонки определяются коэффициенты зная Найдем эти коэффициенты, для этого запишем
Подставим это выражение в (3.4) и, выразив , будем иметь:
(3.8)
Сравнивая формулы (3.7) и (3.8), получим рекуррентные формулы для определения и :
(i=1,2,...,n-1)(3.9)
Определим теперь и . Из (3.7) при i=0 получим:
(3.10)
Сравнивая (3.10) и первое краевое условие из (1.2), находим:
При обратном ходе, используя формулу (3.7) и условие , последовательно находим
Точность в методе конечных разностей достигается аналогичным способом, приведенным в методе сведения краевой задачи к задаче Коши, так как полученная точность зависит от шага разбиения отрезка [a, b].
3.2 Описание результатов
При решении данного дифференциального уравнения второго порядка с заданными краевыми условиями (1.3) методом конечных разностей, получены следующие результаты представленные в таблице 2. В столбце Х приведено разбиение отрезка [1.3; 1.8] с шагом h = 0.02, в столбце Y(X) - значение функции (n=1,…,26) в соответствующих точках , в столбце E - значения найденных абсолютных погрешностей.
В результате работы программы, листинг которой приведен в приложении 2, точность была достигнута при шаге 0.01, максимальная погрешность как видно из таблицы 2 при n=14.
Для достижения заданной точности Е = шаг h = 0.02 уменьшили в 2 раза, поскольку при шаге большем, чем h = 0.01 полученная точность не удовлетворяла заданной.
Легко заметить, что полученная максимальная погрешность меньше заданной точности, следовательно, проверка точности выполнена и полученная точность удовлетворяет заданной.
Таблица 2
XY(X)E
1) 1.30 2.200000 0.00e+000
2) 1.32 2.029191 6.88e-007
) 1.34 1.859823 1.33e-006
) 1.36 1.692044 1.92e-006
) 1.38 1.525996 2.45e-006
) 1.40 1.361810 2.94e-006
) 1.42 1.199608 3.36e-006
) 1.44 1.039506 3.73e-006
) 1.46 0.881607 4.04e-006
) 1.48 0.726010 4.29e-006
) 1.50 0.572804 4.48e-006
) 1.52 0.422069 4.61e-006
) 1.54 0.273879 4.68e-006
) 1.56 0.128301 4.69e-006
) 1.58 -0.014607 4.62e-006
) 1.60 -0.154793 4.50e-006
) 1.62 -0.292211 4.32e-006
) 1.64 -0.426823 4.07e-006
) 1.66 -0.558596 3.77e-006
) 1.68 -0.687504 3.40e-006
) 1.70 -0.813526 2.98e-006
) 1.72 -0.936648 2.50e-006
) 1.74 -1.056858 1.95e-006
) 1.76 -1.174153 1.36e-006
) 1.78 -1.288532 7.06e-007
) 1.80 -1.400000 0.00e+000
4. Метод Галёркина
4.1 Описание метода
В данном методе решения дифференциального уравнения (1.1) с краевыми условиями (1.2) для удобства описания метода введем:
линейный дифференциальный оператор
,
линейные операторы
Тогда постановка задачи для дифференциального уравнения (1.1) с краевыми условиями (1.2) примет вид:
(4.1)
Решение краевой задачи (4.1) будем искать в виде суммы:
(4.2)
где ( i = 1, 2,..., n ) - конечная система базисных функций. Базисные функции должны составлять часть полного класса функций, то есть нет ни одной функции, которую