Реферат: Некоторые дополнительные вычислительные методы
в) Метод конечных разностей (метод сеток) …………………………………………….. 24
г) Разностные схемы для решения уравнения теплопроводности ……………………… 25
д) Разностные схемы для решения уравнения колебания струны ……………………… 26
7. Список литературы ………………………………………………………………………… 27
1. Решение систем линейных уравнений
Системы линейных уравнений (СЛУ) имеют в вычислениях очень большое значение, так как к ним может быть приведено приближенное решение широкого круга задач. Так, основными источниками возникновения СЛУ являются теория электрических цепей, уравнения балансов и сохранения в механике, гидравлике и т.д. Существует несколько способов решения таких систем, которые в основном делятся на два типа: 1) точные методы , представляющие собой конечные алгоритмы для вычисления корней системы, 2) итерационные методы , позволяющие получать корни системы с заданной точностью путем сходящихся бесконечных процессов. Заметим, что даже результаты точных методов являются приближенными из-за неизбежных округлений. Для итерационных процессов также добавляется погрешность метода.
Пример системы линейных уравнений:
Или в матричном виде: ,
где матрица коэффициентов системы;
- вектор неизвестных; - вектор свободных членов.
Схема Халецкого
Запишем систему линейных уравнений в матричном виде: ,
где A=[aij ] – квадратная матрица порядка nи
, - векторы-столбцы.
Представим матрицу Aв виде произведения нижней треугольной матрицы B=[bij ] и верхней треугольной матрицы C=[cij ] с единичной диагональю , где
и .
Тогда элементы bij и cij определяются по формулам
и
Отсюда искомый вектор xможет быть вычислен из уравнений и .
Так как матрицы Bи C – треугольные, то системы легко решаются:
и
Из этих двух формул видно, что числа yi выгодно вычислять вместе с коэффициентами cij . Этот метод получил название схемы Халецкого . В схеме применяется обычный контроль с помощью сумм. Если матрица A – симметрическая aij =aji , то
Пример. Решить систему
Решение.
В первый раздел таблицы впишем матрицу коэффициентов системы, ее свободные члены и контрольные суммы. Далее так как , то первый столбец из раздела 1 переносится в первый столбец раздела II. Чтобы получить первую строку раздела II, делим все элементы первой строки раздела I на элемент, в нашем случае на 3.
Имеем: ; ; ; ; .
Переходим к заполнению второго столбца раздела II, начиная со второй строки. Пользуясь формулами, определяем : ; ; .
Далее определяя по формулам, заполняем вторую сетку для раздела II: