- •1. Метод Данилевского
- •1.1 Метод Гаусса.
- •1.2 Метод Гаусса с выбором главного элемента. ( )
- •1.3 Метод оптимального исключения ( )
- •1.4 Метод прогонки.
- •II. Методы, основанные на разложении матрицы коэффициентов.
- •2.1. Метод квадратного корня (случай эрмитовой матрицы)
- •2.2 Частный случай метода квадратного корня (случай симметричной матрицы)
- •2.3 Схема Халецкого (случай матрицы с отличными от 0 глав. Минорами).
- •Приближенные (итерационные) методы.
- •2.1 Метод простой итерации.
- •2.2 Метод Якоби.
- •2.3 Метод Зейделя.
- •2. Решение полной проблемы собственных чисел.
- •1. Метод Данилевского
- •2. Итерационный метод вращений.
- •Решение нелинейных уравнений.
- •1.2.1 Метод половинного деления (метод дихотомии).
- •1.2.2 Метод Ньютона или метод касательных
- •1) Существование производных 1-го и 2-го порядков;
- •1.2.3 Метод хорд. Метод секущих.
- •1.2.4 Метод простой итерации.
- •III. Решение дифференциальных уравнений.
- •Задачи с начальными условиями для оду (задача Коши).
- •Численные методы решения дифференциальных уравнений Одношаговые методы. Метод Эйлера и его модификации
- •Одношаговые методы. Метод Эйлера и его модификации.
- •1 Метод Эйлера
- •2. Метод Рунге-Кутта.
- •Многошаговые методы
- •1. Формулировка методов
- •Постановка линейных краевых задач для оду 2-го порядка.
- •Конечно-разностный метод решения краевой задачи.
1.2.3 Метод хорд. Метод секущих.
В методах хорд и секущих функция должна удовлетворять тем же условиям, что и в методе касательных. Сходимость и погрешность методов определяется так же, как и в п.1.2.2.
Идея метода хорд состоит в замене кривой хордами, проходящими через концы отрезков, в которых имеет противоположные знаки. В методе хорд требуется, чтобы один конец отрезка, на котором ищется корень, был неподвижен. Заменяя в алгоритме Ньютона (12) производную приближенно отношением: , получим алгоритм метода хорд с неподвижным правым концом (рис.5а):
, (20)
или с неподвижным левым концом (рис.5б):
. (21)
В качестве неподвижного конца ( или ) выбирают тот конец отрезка, для которого знак совпадает со знаком второй производной .
Вычисление приближенных значений корней уравнения следует вести до тех пор, пока не перестанут изменяться те десятичные знаки, которые требуется сохранить в ответе (т.е. пока не будет достигнута заданная степень точности).
Рис. 5.
На рис.5 стрелкой показано направление сходимости.
Для оценки погрешности приближений можно использовать формулу (18) или
, (22)
где
Идея метода секущих состоит в замене кривой прямой, проходящей через точки и . В качестве следующего приближения к корню принимается абсцисса точки пересечения этой прямой с осью (рис.6).
Уравнение прямой
.
Из условия получаем расчетную формулу
. (23)
Формула (23) определяет метод как двухшаговый (результат -го шага зависит от результатов -го и -го шага).
Теорема 5. Метод секущих имеет порядок по крайней мере .
Рис. 6.
Метод секущих является одним из наиболее эффективных итерационных методов решения уравнений (1), так как имеет высокий порядок скорости сходимости в сочетании с минимальными вычислительными затратами.
При применении метода секущих выбор начальной точки нужно осуществлять по тому же принципу, что и в методе касательных (касательная – предельное положение секущей), вторую из начальных точек , требуемую в двухшаговом методе (23), желательно выбирать между и искомым корнем .
Окончание счета по методу секущих, учитывая его быструю сходимость, можно контролировать по формуле (19).
1.2.4 Метод простой итерации.
Применение метода простой итерации требует предварительного приведения уравнения ( ) к виду
, . (24)
Построим график обеих частей уравнения (24): для левой части -- прямая линия – биссектриса первого координатного угла, для правой части – некоторая линия с уравнением (обозначим буквой ). Решением уравнения является абсцисса точки пересечение и биссектрисы.
Допустим, что для имеется начальное приближение . В простейшем варианте метода все дальнейшие приближения строятся по формуле
, (25)
где , например, .
Геометрический смысл процесса вычислений понятен из рис. 7. По находится на точка . Через нее проводится прямая, параллельная оси , и находится точка ее пересечения с биссектрисой. Абсцисса этой точки принимается за следующее приближение к и т.д.
Рис.7.
1. Когда , погрешность и приближение будет отстоять от дальше, чем . Решение будет “точкой отталкивания” для приближений , близких к нему, и в этом случае не будет сходимости последовательности к .
2. Если , то погрешность . Последовательность , если взято достаточно близким к , будет сходиться к приблизительно со скоростью геометрической прогрессии со знаменателем .
При погрешности и будут иметь одинаковые знаки, и сходимость к будет монотонной. Если , то погрешности и имеют разные знаки, и приближения будут сходиться к , колеблясь около .
3. Случай требует специального рассмотрения, так как будет малой величиной высшего порядка по сравнению с . Поэтому, если взято достаточно близким к , то будут весьма быстро сходиться к при возрастании ,
а погрешность будет стремиться к нулю со скоростью, превосходящей сходимость геометрической прогрессии со сколь угодно малым знаменателем. Это часто используют для ускорения сходимости последовательности к путем преобразования заданного уравнения (24) к новому , имеющему то же решение , но такому, что .
Теорема 6. (Достаточные условия сходимости метода простой итерации)
Пусть функция в уравнении (24) определена и дифференцируема на отрезке . Тогда, если существует число такое, что
(26)
на отрезке , то последовательность (25) сходится к единственному корню уравнения (24) при любом начальном приближении .
На практике итерационный процесс останавливают при выполнении условия .