- •Вместо введения: о погрешностях при решении прикладных задач
- •Глава I. Численные методы решения уравнений
- •§ 1. Задача локализации корней
- •Ограничение корней
- •Локализация корней
- •Простейший (грубейший) алгоритм локализации корней:
- •§ 2. Понятие об итерационных методах уточнения корней
- •Метод деления пополам (метод вилки)
- •§ 3. Методы хорд и касательных
- •Метод хорд для монотонных выпукло-вогнутых функций
- •Метод касательных для монотонных выпукло-вогнутых функций
- •§ 4. Метрические и банаховы пространства. Теорема о неподвижной точке
- •Матричные нормы
- •§ 5. Метод простой итерации
- •§ 6. Применение метода простой итерации к решению
- •Условие h([; ]) [; ] :
- •Глава II. Вычисления в линейной алгебре
- •§ 1. Метод Гаусса и его улучшения для повышения точности решения
- •§ 2. Метод простой итерации и метод Зейделя
- •§ 3. Подготовка к применению метода простой итерации
- •§ 4. Проблема собственных значений
- •Глава III. Численное интегрирование
- •§ 1. Метод прямоугольников
- •§ 2. Метод трапеций
- •§ 3. Метод Симпсона (параболическое интерполирование)
- •Глава IV. Некоторые методы аппроксимации функций
- •§ 1. Интерполяционный многочлен Лагранжа
- •§ 2. Интерполяционный многочлен Ньютона
- •§ 3. Метод наименьших квадратов
- •Глава V. Некоторые методы численного решения дифференциальных уравнений
- •Приложение: Сводка характеристик численных методов
- •Характеристики метода:
- •Характеристики метода:
- •Характеристики метода:
- •Характеристики метода:
- •Характеристики метода:
Глава I. Численные методы решения уравнений
В этой главе будут изучены некоторые методы приближённого решения уравнений вида f(x) = 0 , где f : R R – заданная непрерывная функция действительного переменного. Требуется найти все корни этого уравнения. Обычно процесс решения предполагает следующие этапы:
локализация корней;
нахождение приближённых значений корней с заданной погрешностью .
§ 1. Задача локализации корней
Процесс локализации корней подразумевает нахождение в области определения функции конечной системы непересекающихся отрезков [i ; i ] (1 i k), на каждом из которых функция имеет ровно один корень, и любой корень принадлежит одному из рассматриваемых отрезков. Как правило, вначале находят ограниченный интервал [; ], в котором содержатся все корни, т.е. решают задачу ограничения корней, а затем, рассматривая функцию f : [; ] R на этом интервале, уже локализуют её корни.
Ограничение корней
Задача ограничения корней разрешима не всегда, т.к. у функции может быть бесконечное количество корней, которые не возможно заключить ни в один ограниченный отрезок [; ]: простейший пример: f(x) = sin x с корнями ·k (k Z). Даже в случае, когда подобная ситуация не возникает, общих методов решения задачи ограничения корней для непрерывной функции нет. Оценки корней известны лишь для простейших элементарных функций. Например, лемма о модуле старшего члена, изучаемая в курсе алгебры, приводит к следующему результату:
Теорема (о
границах корней многочлена).
Все корни
любого многочлена f(z) = anzn
+ an–1zn–1
+ … + akzk
+ a0
С[z]
степени n,
где a1
= … … = ak–1
= 0, ak
0, лежат в кольце
, где
.
Для произвольной
функции задача ограничения корней
требует аналитической работы. При этом
можно пытаться применять различные
методы решения неравенств f(x)
0 , а также
пользоваться следующим очевидным
результатом:
Теорема (достаточное условие для ограничения корней). Пусть задано уравнение f(x) = 0 с дифференцируемой функцией f : R R .
(1) Если числа , R таковы, что f() 0, f() 0 и f(x) > 0 при x (–∞ ; ) (; +∞), то все корни рассматриваемого уравнения принадлежат отрезку [; ].
(2) Если числа , R таковы, что f() 0, f() 0 и f(x) < 0 при x (–∞ ; ) (; +∞), то все корни рассматриваемого уравнения принадлежат отрезку [; ].
Доказательство. (1) Функция f(x) возрастает на (–∞ ; ) (; +∞), поэтому для x (–∞ ; ) верно f(x) < f() 0, а для x (; ∞) верно f(x) > f() 0. Утверждение (2) доказывается аналогично.
Теорема доказана.
Примеры: 1. Ограничить корни уравнения x – cos 2·x = 0.
Здесь f(x) = x – cos 2·x, и ввиду |cos 2·x| 1 равенство x = cos 2·x возможно лишь при |x| 1. Таким образом, все корни рассматриваемого уравнения принадлежат отрезку [–1; 1].
2. Ограничить корни уравнения ex – 1000·ln x = 0.
Функция f(x)
= ex
– 1000·ln
x
определена при x
> 0, причём
справедливо неравенство
,
т.к. числитель дроби x·ex
– 1000 строго
возрастает при x
> 0 и
= +∞ . Поэтому
можно найти такое
> 0, для
которого f()
> 0: например,
воспользовавшись оценкой x·ex
– 1000 > 2x
– 1000 > 0 при
x
> 9. Кроме
того,
f(9) = e9 – 1000·ln 9 > 2,79 – 1000·3 7625,597484987 – 3000 > 0.
Таким образом, все корни рассматриваемого уравнения принадлежат отрезку [0; 9].
Можно ещё заметить, что при x (0; 1] значения рассматриваемой функции положительны, поскольку ln x 0. Так что найденный интервал [0; 9] можно сузить до [1; 9].
3. Ограничить корни уравнения ln(2·x + 1) – 0,1·x2.
Для функции f(x)
= ln(2·x
+ 1) – 0,1·x2,
определённой при x
> –0,5, имеем
.
Решая
неравенство 2·x2
+ x
– 10 > 0,
находим, что при x
> –0,5
справедливо f(x)
> 0
x
(–0,5; 2). При
этом функция f(x)
возрастает
на (–0,5; 2),
убывает на
(2; +∞), и
,
f(2)
= ln
5 – 0,4 > 0,
.
Ясно, что f(x)
имеет два
корня (один на отрезке (–0,5;
2), второй –
на (2; +∞)).
Поскольку
f(10)
= ln
21 – 10 < 4 – 10 = –6 < 0,
то корни принадлежат (–0,5;
10].
Можно заметить (?!), что один корень – число 0, а второй принадлежит отрезку (2; 5].
