Задача Задача №3593637
А
А
- Итерационный метод и точный метод
- Итерационный метод: строится бесконечная последовательность приближений x^(k), которая стремится к точному решению Ax = b при бесконечном числе шагов; удобно для больших и разреженных систем.
- Точный метод: дает решение в finitely-м количестве арифметических операций (например, Гаусса, формула Крамера, обратная матрица).
- Условие сходимости и роль нормы переходной матрицы
- Сходимость обеспечивается, если для переходной матрицы B в виде x^(k+1) = B x^(k) + d выполняется ρ(B) < 1 (спектральный радиус).
- Роль нормы: достаточное условие — ||B|| < 1 в какой-либо норме; если это выполняется, то сходимость гарантирована (но не обязательно необходима). Диагональное преобладание обеспечивает такое условие через нормы.
- Основная идея метода Якоби
- Каждый неизвестный на шаге k+1 вычисляется из соответствующего уравнения, используя все остальные неизвестные из шага k: x_i^(k+1) = (b_i − ∑{j≠i} a{ij} x_j^(k)) / a_{ii}. Все x_i^(k+1) вычисляются параллельно.
- Основная идея метода Гаусса–Зейделя
- Каждое новое значение x_i^(k+1) вычисляется с использованием недавно полученных значений: x_i^(k+1) = (b_i − ∑{j<i} a_{ij} x_j^(k+1) − ∑_{j>i} a{ij} x_j^(k)) / a_{ii}. Обновления последовательны внутри шага.
- Почему Гаусса–Зейделя быстрее Якоби и память
- GS обычно сходится быстрее, т.к. использует свежие значения сразу же в расчёте последующих компонент.
- Память: GS может работать в одну векторную память (in-place), в то время как Якоби требует две копии векторов (старого и нового).
- Метод релаксации
- Вводится параметр релаксации ω: x^(k+1) = x^(k) + ω (x^(k+1,_GS) − x^(k)).
- Роль ω: ускоряет сходимость при подходящем выборе (часто 1 < ω < 2 для ускорения над Гауссом–Зейделем, но слишком большое ω может привести к расходимости). Выбор ω зависит от задачи.
- Условие диагонального преобладания и его роль
- Условие: для каждой строки i выполняется ∣a_ii∣ > ∑_{j≠i} ∣a_ij∣.
- Обоснование: обеспечивает уменьшение нормы итерационной матрицы (например, ||B||_∞ < 1), следовательно ρ(B) < 1 и сходимость методов Якоби и Гаусса–Зейделя.
- Три основных критерия остановки и их недостатки
- Абсолютная/относительная норма остатка: ||b − Ax^k|| ≤ tol. Недостаток: при плохой обусловленности остаток может плохо отражать реальную погрешность.
- Норма изменения решения: ||x^(k+1) − x^k|| ≤ tol. Недостаток: может «застревать» в малых шагах даже при сущеочной погрешности.
- Максимальное число итераций: k ≤ maxIter. Недостаток: можно досрочно остановиться до достижения нормы; не гарантирует реального решения.
- Когда предпочтительны итерационные или точные методы? Основные критерии выбора
- Предпочтение итерационных: очень крупные или разреженные системы, ограниченная память, требуется приближённое решение за ограниченное время, множество правых частей или необходимость частых пересчётов без полного разложения матрицы.
- Предпочтение точных методов: небольшие или плотные системы, нужна высокая точность за раз, хорошо обусловленная система, возможность заранее разложить LU и многократно использовать результат.