- •1 Системы линейных алгебраических уравнений
- •1.Решение системы существует и является единственным.
- •2. Система уравнений вообще не имеет решений.
- •3. Система уравнений имеет бесконечное множество решений
- •2 Итерационные методы решения систем линейных уравнений
- •Геометрическое истолкование процесса
- •3 Метод исключения (метод Гаусса)
- •4 Задания для самостоятельной работы
- •5 Задания для практической работы
- •6 Задания к лабораторной работе
- •Варианты задания
- •7 Содержание отчета
- •8 Контрольные вопросы
- •1 Формула прямоугольников
- •1.1 Формула «левых» прямоугольников
- •1.2 Формула «правых» прямоугольников
- •1.3 Формула «средних» прямоугольников
- •1.4 Случай неравноотстоящих узлов
- •1.5 Алгоритм вычисления интеграла функции, заданной
- •2 Формула трапеций
- •2.1 Вывод формулы
- •2.2 Оценка погрешности формул прямоугольников и трапеций
- •3 Формула Симпсона
- •4 Формула Гаусса
- •5 Задания для самостоятельной работы
- •5 Задания для практической работы
- •6 Задания к лабораторной работе
- •7 Содержание отчета
- •8 Контрольные вопросы.
- •1 Понятие о приближении функции
- •2 Точечная аппроксимация. Задача интерполирования
- •3 Интерполяционный многочлен Лагранжа
- •4 Интерполяционная формула Ньютона
- •5 Аппроксимация с помощью кусочных полиномов
- •6 Аппроксимация сплайнами
- •7 Задания для самостоятельной работы
- •8 Задания к лабораторной работе
- •9 Содержание отчета
- •10 Контрольные вопросы
- •1 Численное решение алгебраических и трансцендентных уравнений
- •2 Графическое решение уравнений
- •3 Метод половинного деления
- •4 Метод хорд (пропорциональных частей, ложного положения)
- •5 Метод Ньютона (метод касательных)
- •Внеся эту поправку в (2), найдем следующее по порядку приближение корня:
- •6 Видоизмененный метод Ньютона (метод Рыбакова)
- •7 Метод секущих (комбинированный метод секущих-хорд)
- •8 Комбинированный метод касательных-хорд
- •9 Метод последовательных приближений
- •10 Усовершенствованный метод последовательных приближений
- •11 Метод Монте-Карло
- •12 Задания для самостоятельной работы
- •13 Задания к лабораторной работе
- •14 Содержание отчета
- •15 Контрольные вопросы
- •Вычислительные методы Методические указания по проведению лабораторных работ
- •65044, Украина, Одесса, пр. Шевченко, 1
- •65044, Украина, г.Одесса, пр. Шевченко, 1, корп. 5
2 Графическое решение уравнений
Для изучения методов решения нелинейных уравнений важным является знание приближенных графических методов их решения.
Действительные корни уравнения F(x)=0 приближенно можно определить как абсциссы точек пересечения функции y=F(x) с осью Ох.
Если уравнение F(x)=0 не имеет близких между собой корней, то этим способом, как показано на рис.4.2, его корни отделяются.
Рисунок 4.2 – Приближенное определение корней уравнения
На практике часто бывает выгоднее это уравнение заменить равносильным ему уравнением
x=(x),
где функции x и (x) более просты, чем функция F(x).
Тогда, построив графики функций
y=x
и
y=(x),
искомые корни получим как абсциссы точек пересечения этих графиков.
Например, необходимо графически решить уравнение x lgx=1.
Запишем его в виде:
.
Построив кривые и , как показано на рис.4.3, приближенно найдем единственный корень уравнения х 2,5.
Рисунок 4.3 – Графическое решение уравнения x lgx=1
Нахождение корней уравнения упрощается, если одна из функций y=x или y=(x) линейная, т.е., например, хax+b. Тогда корни находятся, как абсциссы точек пересечения кривой y=(x) и прямой y=ax+b. Именно такой графический способ нам потребуется использовать при изучении, например, метода последовательных приближений.
Хотя графические методы решения уравнений весьма удобны и сравнительно просты, но они, как правило, применимы лишь для грубого определения корней. Особенно неблагоприятным в смысле потери точности является случай, когда линии пересекаются под очень острым углом и практически сливаются по некоторой дуге.
3 Метод половинного деления
Пусть дано уравнение F(x)=0, где функция F(x) непрерывна в области своего существования. Требуется найти корень этого уравнения, что графически представляется точкой пересечения графика функции F(x) и оси Х (рис.4.4).
Рисунок 4.4 – Иллюстрация метода половинного деления
Алгоритм метода состоит из следующих операций:
Сначала вычисляются значения функции в точках, расположенных через равные интервалы на оси Х. Это делается до тех пор, пока не будут найдены два последовательных значения функции F(a) и F(b), имеющие противоположные знаки. Если функция непрерывна, изменение знака указывает на существование корня.
Затем по формуле
вычисляется среднее значение Х в интервале значений [а,b] и находится значение функции F(с). Если знак функции F(c) совпадает со знаком F(a), то в дальнейшем вместо F(a) используется F(c). Если же F(c) имеет знак, противоположный знаку F(a), т.е. ее знак совпадает со знаком F(b), то F(b) заменяется на F(c).
Далее находим (для нашего случая):
и так далее.
В результате интервал, в котором заключено значение корня, сужается. Если F(c) достаточно близко к нулю, процесс заканчивается, в противном случае он продолжается.
Часто применяется другая оценка сходимости:
|cn+1 – cn|<ε
Величина F(а) определяется лишь один раз, поскольку нам нужен только знак функции F(х) на левой границе, и он в процессе итераций не меняется.
Метод половинного деления не обладает высокой вычислительной эффективностью. Однако с увеличением числа итераций он обеспечивает получение все более точного приближенного значения корня. После того, как впервые найден интервал, в котором заключен корень, его ширина после N итераций убывает в 2N раз.