- •Слайд-лекции по дисциплине «Численные методы»
- •Лекция 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 – Геометрическая интерпретация метода Рунге-Кутты
- •Соотношения модифицированного метода Эйлера можно записать в форме схемы Рунге—Кутты:
- •Контрольные вопросы
- •Литература
Постановка задачи.
Пусть дана система линейных алгебраических уравнений
а11 x1 |
+ a12 x2 + … + a1n xn = b1 |
а21 x1 |
+ a22 x2 + … + a2n xn = b2 |
…………………………………… |
|
аn1 x1 |
+ an2 x2 + … + ann xn = bn |
Используя понятия матрицы:
a11
A ...
an1
и вектор столбца:
b1 B ...
bn
a12
...
an2
è
...
...
...
X
a1n
...
ann
x1
...
xn
систему уравнений можно переписать в виде:
AX = B.
Правило Крамера состоит в вычислении главного определителя матрицы А:
|
a11 |
a12 |
a13 |
... |
a1n |
D |
... |
|
... ... ... |
||
|
an1 |
an2 |
an3 |
... |
ann |
и побочных определителей:
|
b1 |
a12 |
... |
a1n |
|
a11 |
b1 |
a13 |
... |
a1n |
|
D1 |
... |
... |
... |
... |
, D2 |
... |
|
... |
... |
... |
,... |
|
bn |
an2 |
... |
ann |
|
an1 |
bn |
an3 |
... |
ann |
|
И тогда вычисление корней:
x1 = D1 /D, x2 = D2 /D, … xn = Dn /D.
Число требуемых для вычисления определителей операций:
N n(n! + 1).
2. Метод обратной матрицы
Отталкиваясь от АХ = В, умножим обе части на обратную матрицу A-1 :
A-1 А Х = A-1 В
или
Х = A-1 В.
Таким образом, решение системы сводится к вычислению A-1, но её вычисление также трудоемко.
3. Метод Гаусса
Рассмотрим на примере системы трех уравнений:
а11 x1 + a12 x2 + a13 x3 = b1 а21 x1 + a22 x2 + a23 x3 = b2 а31 x1 + a32 x2 + a33 x3 = b3
Прямой ход метода Гаусса состоит в следующем:
Умножим первое уравнение на – a21 /a11 и сложим со вторым:
a21x1 a21a12 / a11x2 a21a13 / a11x3 a21b1 / a11 a21x1 a22 x2 a23x3 b2
_______________________________________
(a22 a21a12 / a11)x2 (a23 a21a13 / a11)x3 b2 a21b1 / a11
Если обозначить первую скобку за a’22 , а вторую за a’23 и полученную правую часть за b’2 , то получим:
a11 x1 a12 x2 a13 x3 b1
a x |
2 |
a x |
b |
|
22 |
23 |
3 |
2 |
Аналогично: умножим первое уравнение на –a31/a11 , сложим с третьим уравнением
и обозначим полученное уравнение так:
a‘32 x2 + a‘33 x3 = b‘3 .
Обобщенная формула для «штриховых» коэффициентов.
a' |
a |
ai1a1 j |
; |
b' b ai1b1 |
, |
i, j 2,3 |
|
|
|||||||
ij |
ij |
a11 |
i i |
a11 |
|
|
|
|
|
|
|
|
Исключим теперь x2 из 3-го уравнения. Таким образом, система сведется к виду:
a11 x1 a12 x2 a13 x3 b1
|
|
|
|
a22 x2 |
a23 x3 |
b2 |
|
a x |
b |
|
|
33 |
3 |
3 |
|
Для «дважды штриховых» коэффициентов общие формулы имеют вид:
a'' |
a' |
|
a'i2a'2 j |
|
b'' b' |
|
a' b' |
|
||
|
|
; |
|
i2 2 |
, |
i, j 2,3 |
||||
|
|
|||||||||
ij |
ij |
|
a'22 |
|
i |
i |
|
a'22 |
|
|
|
|
|
|
|
|
|
|
Обратный ход метода Гаусса состоит в следующем:
Находим x3 из 3-го уравнения:
x3 = b'‘3 /a'‘33
Потом находим
x2 = (b‘2 – a‘23 x3 )/a‘22 ,
x1 = (b1 – a12 x2 – a13 x3 )/a11 .
Аналогично строится процесс метода Гаусса для системы из N уравнений. Обобщим формулы «штриховых» коэффициентов.
a(k) a(k 1) |
|
a(k 1)ika(k 1)kj ; |
b(k) b(k 1) |
|
a(k 1)ikb(k 1)k , |
|||
ij |
ij |
|
|
a(k 1)kk |
|
i |
i |
a(k 1)kk |
|
|
|
|
|
|
|||
i 2,3,..., n, |
j i 1,i,..., n 1, k |
i 1,i 2,..., n. |
Пример. Методом Гаусса решить систему:
|
10x1 7x2 7 |
|
||
|
|
|
|
|
3x1 2x2 6x3 4 |
||||
|
5x |
x |
5x |
6 |
|
1 |
2 |
3 |
|
4. Итерационный метод Гаусса—Зейделя
Рассмотрим метод Гаусса—Зейделя (на примере трех уравнений)
a11 x1a21 x1a31x1
a12 x2 a13 x3 b1
a22 x2 a23 x3 b2
a32 x2 a33 x3 b3
Выразим x1, x2 и x3 соответственно из 1- го, 2-го и 3-го уравнений:
x1 1 a11 (b1 a12 x2 a13 x3 ) |
||||||||
|
|
1/ a22 (b2 a21x1 a23 x3 ) |
||||||
x2 |
||||||||
x 1/ a (b a x a x |
2 |
) |
||||||
|
3 |
33 |
3 |
31 |
1 |
32 |
|