
- •Элементы теории погрешностей.
- •1.1.Определение абсолютной и относительной погрешности численного результата.
- •1.2.Основные составляющие абсолютной погрешности.
- •1.3.Формулы для оценки абсолютной и относительной погрешности для значения функции и переменных.
- •Решение уравнений с одной неизвестной.
- •2.1.Отделение корней уравнения и отделение корней алгебраического уравнения.
- •2.2.Понятие кратного корня.
- •2.3.Методы уточнения корней простой итерации.
- •2.4.Метод Ньютона.
- •2.5.Метод хорд.
- •2.6.Достаточное условие сходимости метода простой итерации.
- •Вывод условия сходимости:
- •2.7.Оценка погрешности n-приближения.
- •Решение систем линейных уравнений.
- •3.1.2.Классификация слау
- •3.1.3.Обусловленность слау
- •3.1.Метод Гаусса.
- •3.2.Метод прогонки.
- •3.3.Метой простой итерации. Приведение системы к итерационной форме
- •Выполнение итерации
- •Проверка условия окончания решения
- •3.4.Достаточное условие сходимости простой итерации.
- •3.5.Оценка погрешности n-приближения.
- •4.1.Метод простых итераций.
- •1) Функции 1(X,y) и 2(X,y) определены и непрерывно дифференцируемы в r;
- •3) В r выполнены неравенства
- •4.2.Метод Ньютона.
- •Интерполяция функций.
- •5.1.Постановка задачи.
- •5.2.Формула Лагранжа.
- •5.3.Формула Ньютона.
- •5.4.Оценка погрешности интерполяции.
- •5.5.Интерполяция сплайнами.
- •Определение кубического сплайна.
- •Формулировка системы уравнений для коэффициентов кубического сплайна.
- •Аппроксимация функций.
- •6.1.Постановка задачи.
- •6.2.Метод наименьших квадратов.
- •6.3.Матричная формула для коэффициентов многочлена аппроксимации.
- •Численное интегрирование.
- •8.1.Формулы прямоугольников, трапеций, парабол (Симпсона). Оценки погрешностей. Методы прямоугольников
- •Метод трапеций
- •Метод парабол
- •8.2.Метод Рунге двойного счета для оценки погрешности.
- •Методы решения задачи Коши для обыкновенного дифференциального уравнения.
- •9.1.Интегральное представление задачи Коши.
- •10.1.Обусловленность задачи поиска минимума.
- •10.3.Метод золотого сечения.
- •10.4.Метод координатного спуска.
- •10.5.Метод наискорейшего спуска.
3.5.Оценка погрешности n-приближения.
Количество шагов при точности Е из условия: ||A||k+1||F||/(1-||A||)<E, где F-вектор начальных приближений. Решение считается найденным, если два последовательных приближения оказываются достаточно близкими, т.е. норма вектора их разности сравнима с заданной точностью
3.6.Итерации Зейделя.
Этот
метод является модификацией метода
простой итерации. Стратегия этого метода
отличается тем, что при расчете элементов
вектора
наряду с элементами предыдущего
приближения
учитываются
также и элементы, которые уже определены
на данной итерации
.
Расчетная формула метода
3.7.Оценка обусловленности системы линейных уравнений.
Таким образом, малые погрешности вычислений или исходных данных могут привести к существенным погрешностям в решении. Такие системы уравнений называются плохо обусловленными.
Для оценки полученного решения используют следующее соотношение:
3.8.Число обусловленности.
Произведение
называется числом обусловленности
матрицыA
и обозначается как cond(A).
3.9.Оценка погрешности решения.
Решение считается найденным, если два последовательных приближения оказываются достаточно близкими, т.е. норма вектора их разности сравнима с заданной точностью
Решение систем нелинейных уравнений.
Многие практические задачи сводятся к решению системы нелинейных уравнений.
Пусть для вычисления неизвестных x1, x2, … , xn требуется решить систему n нелинейных уравнений:
F1(x1, x2, … , xn) = 0
F2(x1, x2, … , xn) = 0
. . . . .
Fn(x1, x2, … , xn) = 0
В отличие от решения СЛАУ не существует прямых методов решения систем нелинейных уравнений. Лишь в отдельных случаях эту систему можно решить непосредственно. Например, для случая двух неизвестных иногда удается выразить одно неизвестное через другое и таким образом свести задачу к решению одного нелинейного уравнения относительно другого.
В общем случае для решения систем нелинейных уравнений обычно используются итерационные методы.
4.1.Метод простых итераций.
Пусть требуется найти действительные решения системы двух уравнений
F1(x,y) = 0
F2(x,y) = 0
с заданной точностью.
Для этого перепишем исходную систему в приведенном (итерационном) виде:
x = 1(x,y)
y = 2(x,y)
Пусть x0 и y0 - начальные приближения корней, полученные графическим или каким-либо другим способом. Подставив эти значения в правые части приведенной системы уравнений можно получить
x1 = 1(x0,y0)
y1 = 2(x0,y0)
Аналогично можно получить второе приближение
x2 = 1(x1,y1)
y2 = 2(x1,y1)
В общем случае
xn = 1(xn-1,yn-1)
yn = 2(xn-1,yn-1)
Если функции 1(x,y) и 2(x,y) непрерывны и последовательности x1, x2, … , xn, … и y1, y2, … , yn, … сходятся, то пределы их дают решение приведенной, а следовательно, и исходной системы.
Теорема об условиях сходимости итерационного процесса в этом случае выглядит следующим образом:
Пусть в некоторой замкнутой окрестности R (a x A, b y B) имеется одно и только одно решение x = x* и y = y* приведенной системы.
Тогда если:
1) Функции 1(X,y) и 2(X,y) определены и непрерывно дифференцируемы в r;
2) начальные приближения x0, y0 и все последующие приближения xn, yn (n = 1,2, …) принадлежат R;