- •5 Семестр Калуга, 2008
- •Содержание
- •1. Домашнее задание №1: «Методы решения систем линейных алгебраических уравнений»
- •1.1. Задание
- •1.2. Пример выполнения задания Постановка задачи
- •Метод Гаусса с выбором главного элемента.
- •Результат работы алгоритма:
- •Результат работы алгоритма:
- •3. Метод вращений решения.
- •Результат работы алгоритма:
- •4. Метод простых итераций.
- •Пример выполнения программы:
- •5. Метод Зейделя.
- •Результат выполнения программы:
- •Выводы:
- •Требования к оформлению
- •1.4. Список литературы
- •Домашнее задание №2: «Методы решения нелинейных алгебраических уравнений»
- •2.1. Задание
- •2.2. Пример выполнения задания
- •1. Метод биекций.
- •Результат работы программы:
- •2. Метод простых итераций.
- •3. Метод Ньютона.
- •Требования к оформлению
- •2.4. Список литературы
2. Метод простых итераций.
Приведём данное уравнение к виду удобному для итерации
В качестве начального приближения берём левую границу отрезка локализации корня x0= a0 и вычисляем значение функции в этой точке, таким образом, x1= (x0) . Каждое последующее приближение будем искать как xn+1= (xn) .
Рисунок иллюстрирует правило для вычисления следующего приближения по методу простых итераций. На этом рисунке точка М имеет координаты (x1, j (x1)), точка N – координаты (x2, j ( x1).
Условием окончания приближений будет выполнение неравенства:
,
где q берётся исходя из условия:
Скорость сходимости метода – геометрическая.
Выберем отрезок локализации корня [-3;-2] и исследуем изменение значений производной на этом отрезке. Из графика видно, что на выбранном отрезке локализации, следовательно не выполняется одно из условий сходимости метода, а значит в данном случае МПИ применять нельзя.
3. Метод Ньютона.
Пусть задано уравнение
f(x)=0.
Запишем его в виде
, где
и
.
Пусть хк – некоторое приближение к корню х*. Для ускорения сходимости итераций желательно, чтобы был как можно меньше. Положим , то есть
Отсюда находим, что
.
Подставляя в исходное уравнение, получаем рекуррентную формулу:
.
Это и есть итерационная процедура Ньютона. Скорость сходимости метода – квадратичная.
Геометрически метод выглядит так:
Теорема.
Пусть функция удовлетворяет условиям
Тогда, если члены последовательности , определяемые методом Ньютона, при любом фиксированном принадлежат отрезку [a,b] и эта последовательность сходится на [a,b] к корню уравнения, то справедливы неравенства:
Упрощённый метод Ньютона. Итерационный процесс осуществляется без вычислений значений матрицы Якоби на каждом шаге, т.е.: , при этом скорость сходимости метода снижается.
Практическая часть (листинг программы MatLab)
Вывод
Скорость сходимости МПИ – геометрическая, скорость сходимости метода Ньютона – квадратичная, следовательно, метод Ньютона сходится к решению за меньшее число итераций чем в методе Биекций или МПИ. Метод биекций имеет самую низкую сходимость, но прост в реализации.
Полученные результаты:
|
Метод биекций |
МПИ |
Метод Ньютона |
Решение |
-2.1265 |
- |
|
итераций |
12 |
- |
|
Провести сравнительный анализ экспериментальных данных не представляется возможным, так как для данного уравнения МПИ не сходится.
Требования к оформлению
Отчет должен содержать:
Постановку задачи
Алгоритм решения уравнения методом биекций.
Алгоритм решения уравнения методом простых итераций.
Алгоритм решения уравнения методом Ньютона.
Программное обеспечение в среде MatLab и найти решение системы линейных уравнений указанными методами.
Сравнительный анализ методов решений СЛАУ.
Выводы.