
- •Глава 1. Интерполирование
- •Постановка задачи интерполирования
- •1.2. Единственность интерполяционного многочлена
- •. Интерполяционный многочлен Лагранжа
- •1.4. Погрешность интерполяционной формулы Лагранжа
- •. Конечные разности и их свойства
- •1.6. Первая интерполяционная формула Ньютона
- •1.7. Вторая интерполяционная формула Ньютона
- •1.8. Интерполяционная формула Гаусса
- •1.10. Численное дифференцирование
- •Глава 2. Численное интегрирование
- •2.1. Общие замечания
- •Глава 3. Численные методы решения алгебраических и трансцендентных уравнений
- •3.1. Вводные замечания
- •3.2. Отделение корней
- •3.3. Метод половинного деления
- •3.5. Метод Ньютона (метод касательных)
- •3.6. Комбинированный метод
- •3.7. Метод итерации (метод последовательных приближений)
- •Глава 4. Решение систем линейных уравнений
- •4.1. Метод Гаусса (схема единственного деления)
- •4.2. Вычисление определителей
- •4.3. Вычисление обратной матрицы
- •4.4. Некоторые сведения из линейной алгебры
- •Рассмотрим систему линейных алгебраических уравнений
- •Вычитая второе равенство из первого, получаем
- •4.6. Приведение линейной системы к виду,
- •4.7. Метод Зейделя Пусть система линейных алгебраических уравнений приведена к виду
- •Глава 5. Обработка результатов наблюдений
- •5.1. Постановка задачи
- •5.2. Приближение функции, заданной таблично,
- •Глава 6. Численные методы решения
- •6.1. Методы решения задачи Коши. Вводные замечания
- •6.2. Решение дифференциальных уравнений
- •6.3. Метод Эйлера Пусть дано дифференциальное уравнение
- •6.4. Уточненный метод Эйлера Рассмотрим дифференциальное уравнение
- •6.6. Методы Рунге – Кутта Пусть дано дифференциальное уравнение
- •6.7. Метод Милна
- •6.8. Приближенное решение систем дифференциальных
- •6.9.Краевые задачи для обыкновенных дифференциальных
- •6.10. Решение краевых задач
- •6.11. Метод прогонки Рассмотрим линейную систему уравнений (6.54). Преобразуя первые (n-1) уравнения:
- •Глава 7. Численное решение дифференциальных уравнений в частных производных
- •7.1. Классификация дифференциальных уравнений
- •7.2. Метод сеток решения краевых задач
- •7.3. Погрешность аппроксимации дифференциальных уравнений
- •7.4. Аппроксимация граничных условий
- •7.6. Метод сеток для уравнений
4.6. Приведение линейной системы к виду,
удобному для итерации
Задана
система уравнений
,
где
.
Требуется заданную систему уравнений
привести к виду, удобному для итерации.
О
чень
просто преобразуется к нужному виду
система в случае, когда диагональные
элементы матрицы А значительно преобладают
над остальными элементами. Запишем
систему
в развернутом виде:
(4.36)
…………………………………..
Пусть
………………………………
Делим
каждое уравнение системы (4.36) на
диагональный элемент
и находим
:
(4.37)
……………………………………………..
Легко
проверить, что полученная система (4.37)
удовлетворяет первому из достаточных
условий (4.35).
Рассмотрим теперь систему в случае, когда диагональные элементы матрицы А не преобладают над остальными элементами. Из заданной системы выделим уравнения с коэффициентами, модули которых больше суммы модулей остальных коэффициентов уравнения. Каждое выделенное уравнение выписываем в строку новой системы так, чтобы наибольший по модулю коэффициент оказался диагональным. Из оставшихся неиспользованных и выделенных уравнений системы составляем линейные комбинации так, чтобы получить недостающие уравнения, причес диагональные коэффициенты по модулю должны быть больше суммы модулей всех остальных коэффициентов. Необходимо, чтобы при составлении новой системы использовалось каждое уравнение исходной системы. Полученная таким образом система имеет матрицу с преобладающими диагональными элементами, и к ней применимы вышеописанные преобразования.
Естественно, заданную систему можно привести к виду, удобному для итерации, любым другим способом, лишь бы выполнялось одно из условий (4.36).
Пример. Методом простой итерации решить систему
;
;
.
Замечаем, что в третьем уравнении заданной системы коэффициент при х2 по модулю больше суммы модулей коэффициентов при х1 и х2%
.
Ставим уравнение вторым в новой системе. Для получения первого уравнения новой системы сложим первые два уравнения заданной системы:
.
Для получения третьего уравнения новой системы умножим первое уравнение заданной системы - на –2, второе уравнение на 1, третье уравнение - на 2 и все это сложим. Получим
.
Новая система уравнений, эквивалентная заданной, имеет вид
преобразуем ее к виду
.
Полученная система удовлетворяет первому из достаточных условий (4.36).
Возьмем в качестве нулевого приближения х(0) столбец свободных членов. Все вычисления оформим в виде табл.4.5.
Таблица 4.5
|
|
|
|
|
0.6 |
0.6 |
1.875 |
|
0.855 |
-0.3 |
1.95 |
|
0.996 |
-0.009 |
1.981876 |
|
0.998175 |
0.00645 |
1.9995 |
|
0.99861 |
-0.000165 |
1.999771 |
Если остановится на четвертой итерации, то, округляя полученные значения до трех знаков после запятой, имеем ответ:
;
;
.
Отметим, что точным решением системы являются значения:
;
;
.