
- •1. Источники и классификация погрешностей
- •2. Постановка задачи приближённого решения алгебраических уравнен
- •3. Метод половинного деления.
- •4. Метод хорд
- •5. Метод касательной
- •6. Комбинированный метод хорд и касательных
- •7. Решение систем линейных алгебраических уравнений методом Гаусса
- •8. Вычисление определителя и обратной матрицы методом Гаусса
- •9. Решение систем уравнений методом простой итерации
- •10. Метод Зейделя
- •11. Постановка задачи аппроксимации функций
- •12. Интерполяционный многочлен Лагранжа
- •13. Интерполяционные многочлены Ньютона
- •14.Интерполяция сплайнами.
- •15.Квадратурные формулы Ньютона - Котеса. Формула трапеций.
- •16.Квадратурные формулы Ньютона - Котеса. Формула Симпсона.
- •Формула Симпсона
- •17.Квадратурные формулы Гаусса.
- •19.Метод Эйлера.
- •20.Модификации метода Эйлера. Метод Эйлера - Коши.
- •21.Модификации метода Эйлера. Метод серединных точек.
- •22.Метод Рунге - Кутта.
- •23.Методы минимизации функции одной переменной: метод Фибоначчи, метод золотого сечения.
- •24.Методы минимизации функции одной переменной: оптимальный пассивный поиск, метод деления отрезка пополам.
- •25.Многомерные методы оптимизации.
4. Метод хорд
Пусть выполняются следующие условия: f(a)*f(b)<0, функция f имеет непрерывные производные f ΄ ≠ 0 и f ̀̀̀ ̀≠ 0, которые сохраняют постоянный знак на отрезке [a,b] Пусть f ΄ > 0 и f ̀̀̀ ̀> 0 на [a,b], тогда f(a)<0 и f(b)>0.
Пусть
=а
начальное приближение через точки
и В проведем хорду. Ее уравнение
=
Пусть у=0 и х=
,
тогда
=
-
*f(
).
Через точку
(
)
и B
(b;f(b))
проведем хорду. Из ее уравнения при у=0
получим
=
-
*f(
)
и т.д. Таким образом получим итерационную
последовательность, которая вычисляется
рекуррентной формулой.
=
-
*f(
).
N=0,1…
=a.
Эта формула верна и для случая, когда f
̀<0 и f
̀̀ ̀ <0 причем
=a.
Если f
̀и f
̀̀ ̀ различны, то за
=b
и рекуррентная формула имеет вид:
=
-
*f(
).
Погрешность вычисляется по формуле: |
t
-
|
((M-m)/m)*|
-
|
n=0,1…,
где 0
m
|
f
΄(x)
|
M
для x
[a;b]
5. Метод касательной
Пусть
выполняются следующие условия:
f(a)*f(b)<0,
функция f
имеет непрерывные производные f
΄ ≠ 0 и f
̀̀̀ ̀≠ 0, которые сохраняют постоянный
знак на отрезке [a,b]
t[a;b].
Рекуррентная формула имеет вид:
=
-
,
n=0,1,2…
Формула справедлива для всех случаев
относительных знаков f
(x)
и f
''(x)
на отрезке [a,b].
Правило выбора начального приближения:
1. Если f
' и f
'' одного знака, то
=b.
2. Если f
' и f
'' разных знаков, то
=a.
Погрешность
вычисляется по формуле:
|t-
|
(
)^2.
0<m
; 0<|f
''(x)|
для [a,b].
6. Комбинированный метод хорд и касательных
Обозначим
приближения метода хорд через
,
а приближения метода касательных через
или
,
причем
t
и
при n
.
Процесс уточнения до
>0
можно остановить как только окажется
|
и взяв в качестве приближенного корня
середину отрезка между
и
,
т.е. t
=
,
когда вычисления ведутся без заданной
степени точности, то абсолютную
погрешность вычисляют по формуле:
=
,
Комбинированный метод дает возможность
ускорить процесс уточнения корней. Для
этого при вычислении
метод хорд вместо соответствующего
неподвижного конца отрезка [a;b]
используют приближение
найденное
методом касательных. Вычисление пары
чисел (
)
начинают с
,
которое определяется формулой метода
касательных
=
n=0,1,2…
при соответствующем начальном приближении
.
Затем находят
по формуле метода хорд (во всех случаях
относительно знаков производных)
=
-
(n=0,1,2…).
Если f
' и f
'' одного знака, то
=a
и
=b,
если разных, то
=b
и
=a.
7. Решение систем линейных алгебраических уравнений методом Гаусса
Система n линейных уравнений с n независимыми имеет вид:
+
+…+
=
;
+
+…+
=
;
+
+…+
=
.
В матричной форме система имеет вид:
A*X=B;
…
матрица
коэф. системы
матрица – столбец
матрица-столбец
…
Х=
переменных B=
свободных членов
…
Метод
Гаусса – метод последовательного
исключения переменных. Суть метода:
преобразовать исходную систему к
равносильной ей системе с треугольной
матрицей (прямой ход), из которой затем
последовательно записывают значения
всех переменных (обратный ход). Метод
Гаусса относится к точным методам, т.к.
погрешность самого метода равна нулю.
Из-за вычислительных погрешностей (при
округлении) получить точное решение
практически невозможно. Раздел А:
коэффициенты исходной системы и свободные
члены. Для исключения случайных ошибок
предусмотрен текущий контроль, для
этого используют столбец контроля сумм
()
и столбец строчных сумм (S).
Контроль в прямом ходе: после того как
в раздел А внесены коэффициенты и
свободные члены исходной системы находят
контрольные суммы – складывают
коэффициенты и свободные члены по
строкам, результат вносят в раздел ∑.
Над контрольными суммами производят
те же операции, что и над свободными
членами. После выполнения каждого
преобразования находят строчную сумму
результатов и записывают ее в столбец
S(
).
При отсутствии случайных ошибок числа
в столбцах ∑ и S
должны практически совпадать. Пусть b
– вычисляемый элемент нового раздела,
a
– соответствующий элемент предыдущего
раздела,
- горизонтальная проекция элемента b,
- вертикальная проекция элемента b.
Тогда b=a-
.
Контроль в обратном ходе раздел В: в
столбце свободных членов получают
значения переменных, числа в столбце ∑
должны быть ровно на единицу больше
самих переменных. Значения в столбце ∑
можно получить, если при вычислении
переменных заменить свободные члены и
значения уже найденных переменных на
соответствующие числа из столбца ∑.
Вычислим погрешность: для этого подставим
найденные значения в исходную систему
и найдем разности между полученными
значениями и свободными членами
полученной системы. Эти разности
называются невязками.