- •Слайд-лекции по дисциплине «Численные методы»
- •Лекция 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 – Геометрическая интерпретация метода Рунге-Кутты
- •Соотношения модифицированного метода Эйлера можно записать в форме схемы Рунге—Кутты:
- •Контрольные вопросы
- •Литература
В качестве x0 берут обычно один из
концов отрезка [а, b]; если окажется, что'(x) < 0, то может оказаться, что x1 будет
вне отрезка. Это значит, что искомый корень ближе ко второму концу отрезка и именно его следует взять в качестве начального приближения.
Рис. 7 — Графическая интерпретация метода итераций (случай (x) — убывает)
3. Метод хорд
Рассмотрим более быстрый способ нахождения корня на интервале [a,b], в предположении, что f a f b 0
|
|
0 |
f x 0 |
f x |
|
|
f (b) f b 0 |
f (a) f a 0 |
Рис.8а |
Рис. 8б |
Рассмотрим рис.8а. Проведем через точки А и В хорду. Уравнение хорды
|
y f a |
f b f a |
x a |
|||||||
|
|
|
||||||||
|
|
|
|
|
|
b a |
|
|
||
В точке |
x x1 , y 0 |
, в результате получим |
||||||||
первое приближение корня |
|
|
||||||||
|
|
|
|
f |
|
|
|
|
||
|
x1 a |
|
|
|
a. |
|
b |
a |
||
(1) |
f b |
f a |
||||||||
|
|
|
Проверяем условия
а) |
f x1 f b 0 |
б) |
f x1 f a 0 |
Если выполняется условие (а), то в формуле (1) заменяя точку а на х1,,
получим |
|
f x1 |
|
|
b x1 |
|
x2 x1 |
|
|
|
|
||
f b f x1 |
|
|
||||
|
|
|
|
|||
Окончательная |
|
формула |
|
|
для |
n-го |
приближения: |
|
f xn 1 |
|
|
|
|
xn xn 1 |
|
|
|
b xn 1 |
|
|
f b f xn 1 |
|
|
||||
|
|
|
Здесь подвижен конец а, то есть
f xi f b 0
Аналогичная ситуация на рис 9а.
Рассмотрим случай, когда неподвижен конец а.
|
|
0 |
|
|
|
f x |
|
0 |
|||
|
|
|
f x |
|
|
f (b) f b 0 |
f (a) f a 0 |
Рис.9а |
Рис.9б |
На рис 8б, 9б выполняется |
|
|
||||
|
|
f xi f a 0 |
|
|
||
Затем вводим |
|
b1 x1 |
|
(в формуле (1) |
||
точку b заменяем на x1), получим |
||||||
x1 a |
|
f a |
|
x1 a |
||
|
f xn 1 f a |
|||||
Продолжая процесс, придем к формуле |
||||||
xn a |
|
|
f a |
xn 1 a |
||
|
f xn 1 f a |
|
Выполняется до тех пор, пока
xn xn 1 ; xn
- корень уравнения.
Рис. 10
На рис. 10 f x меняет знак, поэтому подвижными будут оба конца.
Теорема. |
|
|
||
Пусть |
|
задана |
непрерывная: |
дважды |
дифференцируемая функция f x |
на [a,b] |
|||
и пусть |
|
f a f b 0 |
|
|
|
|
|
|
|
f x |
|
f x |
сохраняют свои знаки на |
|
а |
|
и |
||
[a,b] |
|
|
|
|
(см. рис 8а, 8б и рис 9а, 9б). Тогда итерационный процесс метода хорд сходится к корню с любой наперед заданной точностью .
Лекция 3. Решение систем линейных алгебраических уравнений.
Содержание
1.Постановка задачи. Правило Крамера.
2.Метод обратной матрицы.
3.Метод Гаусса.
4.Итерационный метод Гаусса- Зейделя.
1. Постановка задачи. Правило Крамера.
Методы их решения разбиваются на две группы:
точные (прямые), позволяющие получить решение системы в точном виде за конечное число арифметических операций. В их числе: Правило Крамера,
метод Гаусса, метод обратной матрицы.
приближенные (итерационные). В них задаются начальные значения и с помощью некоторого алгоритма проводится цикл уточняющих корни вычислений (итерация). Пример: Метод Гаусса—Зейделя.