- •Оглавление
- •Предисловие
- •Лекция 1. Понятие погрешности решения систем алгебраических уравнений
- •1.1. Точные и приближенные числа, источники погрешностей, классификация погрешностей
- •Источники погрешностей
- •1.2. Абсолютные и относительные погрешности
- •1.3. Понятие о системе линейных уравнений
- •2.1. Формулы Крамера для решения системы линейных уравнений
- •2.2. Решение системы линейных уравнений методом последовательного исключения неизвестных (метод Гаусса)
- •Пусть дана система
- •1.4. Матричные уравнения
- •Лекция 2. Точные методы решения систем линейных алгебраических уравнений (слау)
- •Обозначим
- •Лекция 3. Приближенные методы решения слау
- •3.1. Понятие предела для векторов и матриц Пусть дана последовательность векторов
- •Аналогично, если имеется последовательность квадратных матриц
- •3.2. Метод простой итерации Пусть дана система линейных уравнений:
- •В матричном виде:
- •Условия сходимости итерационного процесса
- •3.4. Оценка погрешности приближенного процесса методом итерации
- •Лекция 4. Метод зейделя
- •4.1. Условия сходимости процесса Зейделя
- •Пусть дана линейная система
- •4.2. Оценка погрешности процесса Зейделя
- •4.3. Приведение системы линейных уравнений к виду, удобному для итераций
- •Лекция 5. Методы решения нелинейных уравнений
- •5.1. Итерационные методы решения нелинейных уравнений а). Метод простой итерации.
- •Число итераций при использовании этого метода
- •Если корень нас не устраивает, то мы находим
- •Г). Метод Ньютона.
- •5.2. Метод Рыбакова
- •5.3. Метод наискорейшего спуска
- •Лекция 6. Интерполирование и экстраполирование
- •6.1. Интерполяционный многочлен Лагранжа
- •Перепишем этот многочлен в другой форме
- •Окончательно получаем
- •Лекция 7. Конечные разности
- •Рассмотрим некоторые свойства конечных разностей
- •7.1. Первая интерполяционная формула Ньютона для равноотстоящих узлов интерполяции
- •7.2. Вторая интерполяционная формула Ньютона
- •7.3. Оценка погрешностей интерполяционных формул Ньютона
- •Лекция 8. Линейное интерполирование
- •8.1. Интерполирование по Эйткину
- •8.2. Интерполяция и приближение сплайном
- •Приближение линейными сплайнами
- •Лекция 9. Численное интегрирование функций
- •9.1. Постановка задачи
- •9.2. Формула прямоугольников
- •Лекция 10. Численные методы решения задачи коши для обыкновенных дифференциальных уравнений
- •10.1. Задача Коши. Общие замечания. Постановка задачи
- •10.2. Метод Эйлера
- •10.4. Метод Рунге-Кутта
- •Лекция 11. Разностный метод решения уравнений математической физики (метод сеток)
- •11.1. Метод сеток для уравнения параболического типа
- •Лекция 12. Метод сеток для уравнения гиперболического типа
- •12.2. Метод сеток для уравнений Пуассона и Лапласа
- •Образцы выполнения контрольных работ Контрольная работа № 5
- •Контрольная работа № 6
- •Решение уравнения методом Эйлера
- •Решение уравнения модифицированным методом Эйлера
- •Решение уравнения методом Рунге-Кутта
- •Аналитическое решение заданного уравнения
- •Сравнение точного решения и приближенных решений исходного дифференциального уравнения
- •Контрольная работа №7
- •Контрольная работа №8
- •Решение дифференциального уравнения методом Рунге-Кутта
- •Аналитические решения заданного уравнения
- •Примеры заданий для самостоятельной работы Задание 1
- •Задание 2
- •Задание 3
- •Задание 4
- •Задание 5
- •Варианты задания
- •Задание 6
- •Варианты задания
- •Задание 7
- •Варианты задания
- •Задание 8
- •Варианты задания
- •Задача для самостоятельного решения
- •Вопросы для самостоятельной подготовки
- •Библиографический список
5.2. Метод Рыбакова
Можно
рассматривать этот метод как модификацию
метода Ньютона. При замене
некоторым числом
,
где
– значение х
на [a,
b]
, при котором производная
максимальна.
При
сходимость не нарушается, но замедляется.
Метод Рыбакова удобен для поиска всех корней уравнения f(x) = 0 на [a,b] .
1. Задаем начальные значения х = х0 = а .
2. Для каждой последовательной итерации ( n = 0, 1, 2, …) вычисляем
и
проверяем условие xn
< b,
если оно выполняется, то, значит, найдены
все корни, в противном случае проверяем
выполнение условия
.
Если оно не выполняется, то повторяем
цикл с пункта 2 и переходим к пункту 3.
3.
Задаем начальное приближение
и снова идем на пункт 2.
5.3. Метод наискорейшего спуска
Распространенным методом минимизации функций большого числа переменных является метод градиентного спуска. Последующее приближение получается из предыдущего смещением в направлении, противоположном градиенту функции F(x). Каждое следующее приближение ищется в виде
.
Приведенное
описание не определяет алгоритм
однозначно, т.к. ничего не сказано о
выборе параметра
. Его можно определять из условия минимума
величины
.
В этом случае рассматриваемый метод называют методом наискорейшего градиентного спуска.
Для
функции
,
соответствующей системе линейных
уравнений с матрицей
,
задача нахождения минимума решается
в явном виде
так как
и
.
Обозначим
через
, т.е.
.
Предположим,
что
. Учитывая, что
,
вычислим
:
,
.
Лекция 6. Интерполирование и экстраполирование
Математическая постановка задачи интерполирования.
В экономике и технике, часто приходится сталкиватьcя с необходимостью вычисления значений функции у = f(х) в точках, отличных от значения аргумента, фиксированных в таблице.
Подобные задачи практики формализуются как математические задачи интерполирования.
Пусть
на отрезке [a,
b]
задана
функция у
= f(х)
своими n
+ 1
значениями
;
;
… ;
в точках x0,
x1,
…, xn
, которые
назовем узлами интерполяции.
Требуется найти аналитическое выражение табулированной функции F(х), совпадающее в узлах интерполяции со значениями заданной функции,
y0 |
y1 |
… |
yn-1 |
yn |
x0 |
x1 |
… |
xn-1 |
xn |
т.е.
;
;
… ;
.
Процесс вычисления значений функций в точках х, отличных от узлов интерполяции, называется интерполированием функции f(х).
Если аргумент х находится за пределами отрезка интерполирования [x0 , xn], то задача определения значения функции в точке х называется экстраполированием.
Задача становится однозначной, если в качестве интерполирующей функции f(х) для функции у = f(х), заданной своими n+1 значениями, выбрать многочлен Fn (x) степени не выше n, такой, что
;
;
… ;
.
Многочлен Fn(x) удовлетворяющий этим условиям, называют интерполяционным многочленом, а соответствующие формулы – интерполяционными.
В случае, когда F(x) выбирается в классе степенных функций, интерполяция называется параболической.
При интерполировании возникает ряд задач:
Выбор наиболее удобного способа построения интерполяционной функции для каждого конкретного случая.
Оценка погрешности при замене f(x) интерполирующей функцией F(x) на отрезке [a, b].
Оптимальный выбор узлов интерполяции для получения минимальной погрешности.
