Мегаобучалка Главная | О нас | Обратная связь


Алгоритм схемы Халецкого



2015-11-07 2159 Обсуждений (0)
Алгоритм схемы Халецкого 4.67 из 5.00 6 оценок




· Вводим матрицу A и вектор f.

· Обнуляем вспомогательные матрицы B и C.

· Задаем первый столбец B и первую строку C по формулам (3.7) и (3.8).

· Вычисляем остальные элементы матриц B и C.

Þ Перебор строк (цикл i от 1 до N)

Þ Перебор столбцов (цикл j от 1 до N)

Þ Если , то вычисляем по формуле (3.7).

Þ Если , положим .

Þ Если , то вычисляем по формуле (3.8).

Þ Конец перебора строк и столбцов.

· Проверяем правильность вычисления матриц B и C. Для этого достаточно вычислить произведение и сравнить с A.

· Определим вектор y по формуле (3.9).

· Вычислим решение x по формуле (3.10).

· Выведем матрицу A и вектор правой части f, решение x и невязку Ax-f на экран.

Вычисление невязки решения

Для проверки правильности решения системы линейных алгебраических уравнений необходимо вычислить невязку решения. Невязка v является вектором и вычисляется по формуле

, (3.11)

где - приближенное решение системы уравнений Ax=f.

Если решение является точным решением, то невязка v будет равна нулевому вектору. При решении задачи на компьютере, точное решение получить практически невозможно из-за ошибок округления. Приближенное решение является приемлемым, если его невязка не будет превосходить заранее заданной погрешности Eps. Обычно погрешность Eps равна 0.001, но это значение можно изменять в зависимости от поставленной задачи и требуемой точности решения.

Итерационные методы

Итерационные методы основаны на построении итерационной последовательности , сходящейся к искомому решению системы .

Каждый такой метод характеризуется своей итерационной формулой, позволяющей вычислять очередное приближение по ранее найденным.

В простейшем случае при вычислении используется только одно приближение . Такие методы называют одношаговыми или двухслойными. Итерационную формулу для одношаговых методов принято записывать в стандартной канонической форме.

. (3.12)

где - итерационный параметр ( >0), - вспомогательная невырожденная матрица.

Метод называется стационарным, если для и .

Стационарные методы - наиболее простые с точки зрения организации вычислительного процесса. Однако нестационарные методы имеют другие преимущества, связанные с выбором и для ускорения сходимости итерационного процесса.

Определение очередного приближения с помощью итерационной процедуры (3.12) требует решения уравнения

, (3.13)

где

.

Такую процедуру приходится выполнять на каждом шаге. Поэтому итерационный метод можно применять при условии, что определение отдельных членов итерационной последовательности требует существенно меньшего объема вычислений, чем прямое решение исходной системы. Это накладывает ограничения на выбор матрицы .

Самый простой случай, когда . Тогда формула (3.13) приобретает вид

. (3.14)

Итерационные методы, со схемой вычислений (3.14), называются явными.

Среди неявных методов ( ) наибольшее распространение получили методы с треугольной матрицей .

Разность между приближенным решением и точным решением называется погрешностью .

Погрешность удовлетворяет итерационной формуле:

.

Итерационный процесс сходится , если при .

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

Матрица называется симметричной ( ), если для всех .

Матрица называется положительно определенной ( ), если для имеем .

Матрица называется матрицей с диагональным преобладанием, если ее элементы удовлетворяют неравенству:

для любого .

Примеры:

, и .

- симметричная матрица, и - несимметричные.

- не положительно определенная, т.к. для имеем .

- положительно определенная, т.к. для получаем .

и - без диагонального преобладания, - с диагональным преобладанием.



2015-11-07 2159 Обсуждений (0)
Алгоритм схемы Халецкого 4.67 из 5.00 6 оценок









Обсуждение в статье: Алгоритм схемы Халецкого

Обсуждений еще не было, будьте первым... ↓↓↓

Отправить сообщение

Популярное:
Личность ребенка как объект и субъект в образовательной технологии: В настоящее время в России идет становление новой системы образования, ориентированного на вхождение...
Почему люди поддаются рекламе?: Только не надо искать ответы в качестве или количестве рекламы...
Как построить свою речь (словесное оформление): При подготовке публичного выступления перед оратором возникает вопрос, как лучше словесно оформить свою...



©2015-2020 megaobuchalka.ru Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. (2159)

Почему 1285321 студент выбрали МегаОбучалку...

Система поиска информации

Мобильная версия сайта

Удобная навигация

Нет шокирующей рекламы



(0.008 сек.)