Методы решения нелинейных уравнений
Автор: Wearo Куколдов • Ноябрь 14, 2020 • Лабораторная работа • 421 Слов (2 Страниц) • 420 Просмотры
Московский Технический Университет Связи и Информатики
Кафедра информатики
[pic 1]
Лабораторная работа №1 по теме:
«Методы решения нелинейных уравнений»
Выполнил: Кожемякин Н.А. УБСС1803
Вариант: №9
Москва 2019
Задание.
Решите уравнение.
[pic 2].
Отделение корней с использованием MathCAD.
Отделение корней производим графическим методом с обязательным подтверждением результата аналитически.
[pic 3]
На отрезке [1,1; 1,2] функция f(x) меняет знак, т.е. существует, по крайней мере, один корень. Поскольку знак первой производной f1(x) на выбранном отрезке остается постоянным, то можно сказать, что функция на этом отрезке монотонна. Знакопостоянство второй производной f2(x)<0 на выбранном отрезке является необходимым условием применения метода Ньютона и метода хорд. Следовательно, заданное уравнение имеет единственный корень на отрезке [1,1; 1,2].
Уточнение корня с использованием MathCAD
Метод половинного деления
- Исследование задания.
- Метод половинного деления сходится, если на выбранном отрезке отделен один корень. Так как на отрезке [1,1; 1,2] функция [pic 4] меняет знак ([pic 5]) и монотонна (f’(x)>0), то условие сходимости выполняется.
- Выберем за начальное приближение середину отрезка [pic 6].
- Условие окончания процесса уточнения корня. Для оценки погрешности метода половинного деления справедливо условие |bn–an|<ε , т.е. длина отрезка, полученного на n-ом шаге должна быть меньше заданной точности - [pic 7].
- Результаты «ручного расчета» четырех итераций.
[pic 8]
[pic 9]
Результаты вычислений представлены в форме табл. 1.
Таблица 1. Результаты вычислений
n | a | b | f(a) | f(b) | (a+b)/2 | f((a+b)/2) | b-a |
1 | 1,1 | 1,2 | -0,219 | 0,528 | 1,15 | 0,15856 | 0,1 |
2 | 1,1 | 1,15 | -0,219 | 0,159 | 1,125 | -0,02916 | 0,05 |
3 | 1,125 | 1,15 | -0,029 | 0,159 | 1,1375 | 0,06495 | 0,025 |
4 | 1,125 | 1,1375 | -0,029 | 0,065 | 1,13125 | 0,01796 | 0,0125 |
5 | 1,125 | 1,13125 | -0,029 | 0,018 | 1,12816 | -0,00558 | 0,00625 |
После трех итераций приближение к корню – середина отрезка x3=1,13125.
Вычислим значение корня с помощью встроенной функции root.
...