- •Слайд-лекции по дисциплине «Численные методы»
- •Лекция 1.
- •1. Постановка задачи
- •Методы решений нелинейных уравнений делят на:
- •Три этапа поиска корня:
- •2. Отделение корня
- •Рис. 1 — Графическое отделение корня
- •Для отделения корня можно активно использовать компьютер.
- •3. Уточнение корня. Метод деления отрезка пополам
- •Пусть на интервале [a,b] расположен один корень.
- •Через n итераций интервал будет равен bn an 21n b a
- •Графически метод половинного деления выглядит следующим образом
- •Рис. 4 — Блок-схема метода деления отрезка пополам
- •Лекция 2. Уточнение корня нелинейного
- •1. Метод Ньютона (метод касательных)
- •Рассмотрим полученный треугольник Aax. Для него:
- •Достоинство метода — быстрая сходимость. Если выполняется условие: F''(х) (вторая производная) сохраняет знак
- •Рис. 6 — Блок-схема метода Ньютона
- •2. Метод итераций (простых)
- •В качестве x0 берут обычно один из
- •3. Метод хорд
- •Рассмотрим рис.8а. Проведем через точки А и В хорду. Уравнение хорды
- •Если выполняется условие (а), то в формуле (1) заменяя точку а на х1,,
- •На рис 8б, 9б выполняется
- •Теорема.
- •Лекция 3. Решение систем линейных алгебраических уравнений.
- •1. Постановка задачи. Правило Крамера.
- •Постановка задачи.
- •Используя понятия матрицы:
- •Правило Крамера состоит в вычислении главного определителя матрицы А:
- •2. Метод обратной матрицы
- •3. Метод Гаусса
- •Для «дважды штриховых» коэффициентов общие формулы имеют вид:
- •Аналогично строится процесс метода Гаусса для системы из N уравнений. Обобщим формулы «штриховых»
- •4. Итерационный метод Гаусса—Зейделя
- •Зададим начальные (нулевые) приближения для неизвестных:
- •Лекция 4. Численное интегрирование.
- •1. Метод прямоугольников.
- •где исходный отрезок [a, b] с помощью точек
- •есть площадь так называемой криволинейной трапеции,
- •Рис. 13 — Блок-схема метода прямоугольников
- •Ошибка формулы прямоугольников:
- •Рис. 14 — Иллюстрация метода прямоугольников на примере четверти круга
- •Программа на языке Паскаль, реализующую метод прямоугольников на этом примере:
- •2. Метод трапеций
- •Значит интеграл по всему отрезку [a, b]:
- •Ошибка метода трапеций:
- •Рис. 16 — Блок-схема алгоритма метода трапеций
- •3. Метод Симпсона
- •Аналогично,
- •Ошибка её:
- •Лекция 5. Численное решение дифференциальных уравнений
- •1. Постановка задачи
- •Решение дифференциальных уравнений ищется не в виде аналитической функции, а в виде набора
- •2. Метод Эйлера
- •В методе Эйлера движение в каждой следующей точке происходит по касательной к кривой,
- •3.Метод Рунге-Кутты
- •Первый
- •Наибольшее распространение в вычислительной практике нашел метод Рунге—Кутты четвертого порядка:
- •Схема
- •Прямая АЕ, проходящая через точку
- •Рис. 17 – Геометрическая интерпретация метода Рунге-Кутты
- •Соотношения модифицированного метода Эйлера можно записать в форме схемы Рунге—Кутты:
- •Контрольные вопросы
- •Литература
Пусть на интервале [a,b] расположен один корень.
Метод деления отрезка пополам заключается в следующем. Имеем:
|
|
|
|
|
f a f b 0 |
c 1 a b |
|
|
Определяем половину отрезка |
|
|||||||
и вычисляем |
f c . |
|
2 |
|
||||
|
|
|
||||||
|
Проверяем следующие условия: |
|
||||||
Если |
|
f c |
|
, то c – корень. Здесь |
|
|||
|
|
|||||||
|
|
|
|
|
|
- заданная точность. |
|
|
|
f a f с 0 |
|
||||||
Если |
, то корень лежит в |
|
||||||
|
|
[a,c] |
. |
|
|
|
||
интервале |
|
|
|
|
|
|||
Если |
f c f b 0 |
, то корень лежит на |
|
|||||
отрезке [c,b] . |
|
|
|
|
Через n итераций интервал будет равен bn an 21n b a
Если выполняется условие
bn an / 2
то процесс поиска заканчивается и корень уравнения будет
12 bn an
Графически метод половинного деления выглядит следующим образом
| f c | |
f a f c 0 |
f b f c 0 |
|
Рис. 3 – Графический метод |
|
Рис. 4 — Блок-схема метода деления отрезка пополам
Лекция 2. Уточнение корня нелинейного
уравнения.
Содержание
1. Метод Ньютона (метод касательных)
2.Метод итераций (простых)
3.Метод хорд
1. Метод Ньютона (метод касательных)
Вметоде касательных проводят касательную
ккривой y = F(x) в одной из точек a или b отрезка, где имеется корень. Затем находят точку пересечения касательной с осью OX
Рассмотрим полученный треугольник Aax. Для него:
tg F(a) /(x a)
Учитывая, что
F' (a) = tg(π – α)= –tg(α).
и обобщив эти формулу, можем записать:
xn 1 xn F(xn ) / F (xn )
Процесс следует продолжать, пока не будет достигнута требуемая точность. Пока не станет
|F(xn) | < ε (заданная точность).
Достоинство метода — быстрая сходимость. Если выполняется условие: F''(х) (вторая производная) сохраняет знак и при этом F(x) и F''(х) имеют одинаковый знак, то процесс сходится так, что на каждой итерации число верных значащих цифр удваивается.
Недостаток — на каждой итерации велик объём вычислений (нужно, в частности, знать выражения для производной и так далее).
Рис. 6 — Блок-схема метода Ньютона
2. Метод итераций (простых)
Состоит в замене исходного уравнения f(x) = 0
эквивалентным ему уравнением
x= (x)
ипостроением последовательности
xn+1 = (xn),
которая при n сходится к точному решению.
Условием сходимости этого процесса является:
| '(x)| < 1.