- •Методические указания
- •Часть 1
- •1. Численное решение нелинейных уравнений.
- •1.1. Метод деления отрезка пополам.
- •1.2. Метод Ньютона (метод касательных).
- •1.3. Метод простой итерации.
- •2. Методы решения систем линейных алгебраических уравнений.
- •2.1. Метод Гаусса.
- •2.2. Метод обратной матрицы.
- •2.3. Метод прогонки.
- •2.4. Метод простой итерации (метод Якоби).
- •2.5. Метод Зейделя.
- •3. Численные методы решения систем нелинейных уравнений.
- •3.1. Метод простой итерации (метод Якоби).
- •3.2. Метод Зейделя.
- •3.3. Метод Ньютона.
- •Литература
1.1. Метод деления отрезка пополам.
Допустим, что мы
нашли отрезок
,
в котором расположено искомое значение
корня
,
т.е.
.
Пусть для
определенности
,
(рис. 1.1). В качестве начального приближения
корня
принимается
середина этого отрезка, т.е.
.
Далее исследуем значение функции
на концах отрезков
и
.
Тот из них, на концах которого
принимает значения разных знаков,
содержит искомый корень. Поэтому его
принимаем в качестве нового отрезка.
Вторую половину отрезка
отбрасываем. В качестве первой итерации
корня принимаем середину нового отрезка
и т. д.
|
|
|
Рис. 1.1 Метод деления отрезка пополам. |
Пример 1.1.
Найти решение уравнения
c точностью
методом деления отрезка пополам.
|
|
|
|
|
Рис. 1.2. Графический метод изоляции корня уравнения |
Начальное
приближение:
,
,
.
;
;
.
1-е приближение:
,
,
.
Погрешность
.
;
;
.
Корень находится
в интервале
.
2-е приближение:
,
,
.
Погрешность
.
;
;
.
Корень находится
в интервале
.
…
7-е приближение:
,
,
.
Погрешность
.
Приближенным
решением данного уравнения является
.
На рис. 1.3 приведена программа решения данного уравнения методом деления отрезка пополам. В качестве исходных данных вводятся границы интервала, содержащего корень и точность вычисления.
|
CLS DEF FNF(X)=X^3+X-1 1 INPUT A, B, E IF FNF(A)*FNF(B)>0 THEN 1 2 X =(A+B)/2 IF FNF(A)*FNF(X)<0 THEN B=X ELSE A=X IF (B-A)>E THEN 2 PRINT X,FNF(X) END |
|
Рис. 1.3. Пример программы нахождения корней уравнения методом деления отрезка пополам на языке QUICK BASIC. |
Пример 1.2.
Найти решение уравнения
c точностью
методом деления отрезка пополам с
помощью программыExcel.
Найдем интервал,
содержащий единственный корень уравнения.
Для этого необходимо построить таблицу
или график функции
.
Введем в ячейки A2, A3, A4, … значения переменной
.Введем в ячейку B2 формулу =–A2^3+A2–1.
Скопируем формулу и вставим в остальные ячейки столбца B.
Найдем соседние ячейки, в которых значения функции имеют разные знаки (рис. 1.4 а). Соответствующие значения переменной
дают границы интервала, содержащего
корень.Для построения графика вызываем мастер диаграмм. Выбираем тип диаграммы «точечная» - точечная диаграмма со значениями, соединенными сглаживающими линиями.
Границы интервала, содержащего корень, соответствуют значениям шкалы, между которыми линия графика пересекает горизонтальную ось (рис. 1.4 б)
|
a) |
б) | |||||||||||||||||||||||||||
|
|
| |||||||||||||||||||||||||||
|
Рис. 1.4. Изоляция корня уравнения в Excel с помощью: а) таблицы; б) графика. Искомый корень находится в интервале [0; 1]. | ||||||||||||||||||||||||||||
-
A
B
C
D
E
F
G
1
a
b
x
b-a
F(a)
F(b)
F(x)
2
0,0000
1,0000
0,5000
1,0000
-1,0000
1,0000
-0,3750
3
0,5000
1,0000
0,7500
0,5000
-0,3750
1,0000
0,1719
4
0,5000
0,7500
0,6250
0,2500
-0,3750
0,1719
-0,1309
5
0,6250
0,7500
0,6875
0,1250
-0,1309
0,1719
0,0125
6
0,6250
0,6875
0,6563
0,0625
-0,1309
0,0125
-0,0611
7
0,6563
0,6875
0,6719
0,0313
-0,0611
0,0125
-0,0248
8
0,6719
0,6875
0,6797
0,0156
-0,0248
0,0125
-0,0063
9
0,6797
0,6875
0,6836
0,0078
-0,0063
0,0125
0,0030
Рис. 1.5. Решение уравнения методом деления отрезка пополам с помощью программы Excel.
Продолжаем решение на новом листе (рис. 1.5).
Ввести в ячейки A1 – G1 заголовки столбцов.
В ячейку A2 – значение левой границы интервала

В ячейку B2 – значение правой границы интервала

В ячейку C2 – формулу середины отрезка
=(A2+B2)/2В ячейку D2 – формулу погрешности =B2–A2
В ячейку E2 – формулу функции =A2^3+A2-1.
Скопировать формулу из E2 в ячейки F2 и G2. Строка 2 теперь содержит результаты начального приближения
В ячейку A3 – формулу =ЕСЛИ(E2*G2<0;A2;C2)
В ячейку B3 – формулу =ЕСЛИ(E2*G2<0;C2;B2)
Выделить ячейки C2:G2 и скопировать формулы в соседние ячейки C3:G3 при помощи маркера заполнения (небольшой черный квадрат в правом нижнем углу выделенного блока). Строка 3 теперь содержит результаты первого приближения
Выделить ячейки A3:G3 и скопировать формулы в соседние ячейки расположенных ниже строк A4:G4, A5:G5, и т.д. при помощи маркера заполнения. Каждая новая строка содержит результаты очередного приближения.
В столбце С найти значение корня, соответствующее заданной точности.
Приближенное
решение данного уравнения
содержится в ячейкеС9
(погрешность
в ячейкеD9).



1
