Метод бисекции в маткаде

Метод бисекции в маткаде

Отделить корни уравнения

графически и уточнить один из них:

· методом половинного деления;

· методом простой итерации;

с точностью ε=0,001.

Создать функции, реализующие указанные методы, построить графическую иллюстрацию методов, результаты проверить с помощью встроенных функций, оценить точность полученных значений.

Метод бисекции или метод деления отрезка пополам — простейший численный метод для решения нелинейных уравнений вида F(x)=0. Предполагается только непрерывность функции F(x).

Графическое представление метода бисекций

Решим задание в пакете Маткад:

Идея метода состоит в том, что по двум точкам и построить прямую (то есть хорду, соединяющую две точки графика) и взять в качестве следующего приближения абсциссу точки пересечения этой прямой с осью Ox .

Графическое представление метода хорд

Решим задание в пакете Маткад:

3. Метод касательных (Ньютона)

Графическое представление метода касательных

Решим задание в пакете Маткад:

4. Метод секущих

Графическое представление метода секущих

Решим задание в пакете Маткад:

5. Метод простой итерации

Графическое представление метода простой итерации

Решим задание в пакете Маткад:

Решить задачу Коши для дифференциального уравнения на отрезке [a,b] при начальном заданном условии и шаге интегрирования h:

1. методом Эйлера;

2. методом Рунге – Кутта 4 – го порядка точности.

3. проверить решение с помощью встроенных функций пакета MathCAD.

В решении оставлять 5 цифр после запятой.

Решим задание в пакете Маткад:

2. Метод Рунге – Кутта 4 – го порядка точности

Проверим решение с помощью встроенных функций пакета MathCAD

Лабораторная работа №3

РЕШЕНИЕ НЕЛИНЕЙНЫХ УРАВНЕНИЙ

Теоретические сведения

Уравнение – это равенство, которое выполняется лишь при некоторых значениях переменных, которые входят в уравнение.

Нелинейное уравнение – уравнение, в котором неизвестное находится в степени, не равной 1, или является аргументом какой-то функции.

Общий вид нелинейного уравнения:

где – некоторая нелинейная функция.

Если функция непрерывна на отрезке , то уравнение может иметь корни на этом отрезке при выполнении двух условий:

1) функция меняет знак на отрезке . Это легко проверить выполнением неравенства ;

2) первая и вторая производная сохраняют знак на всем отрезке.

Если какая-либо точка преобразовывает уравнение в равенство, то точка называется корнем нелинейного уравнения. Геометрически точка соответствует точке, в которой график функции пересекается с осью .

Приближенное нахождение действительных корней состоит из 2 этапов:

1. Отделение корня (нахождение промежутка, в котором существует корень)

2. Уточнение корней с заданной степенью точности.

Для отделения корней применяют следующие способы:

— формулы Лагранжа для полиномов.

Для уточнения корней применяют следующие численные методы:

— метод половинного деления (дихотомии или бисекции);

— метод Ньютона (касательных);

— метод итераций или последовательных приближений;

Метод половинного деления

Дано: уравнение и интервал .

Находим точку :

Выполняем проверку: если , то тогда , а если , то тогда . Если такую проверку сделать многократно, то мы будем сужать интервал до тех пор пока не найдем такую точку, для которой значение функции будет равно нулю, а модуль не будет превышать заданную точность .

Метод хорд.

Суть метода состоит в лианезации кривой хордой. Решение задачи покажет уравнение прямой, проходящей через 2 точки и :

Выбираем тот из промежутков или , на концах которого знаки противоположны.

С выбранным промежутком повторяем те же действия, до тех пор, пока и модуль не превышает заданной точности.

После многократных повторений получаем конечную формулу:

при .

Метод касательных.

Суть метода состоит в лианезации кривой касательной. Касательную проводят в том конце кривой, где знак второй производной совпадает со знаком функции, то есть:

Решением является точка пересечения касательной с . Решение получают из уравнения касательной, проходящей через соответствующую точку. После многократного повторения получаем конечную формулу:

Действия повторяются до тех пор, пока и модуль не превышает заданной точности.

Итерационный метод.

Для использования итерационного метода разбиваем таким образом, чтобы можно было записать .

Итерационный процесс записывается следующим образом:

Условием сходимости итерационного процесса является:

Решение уравнений средствами Mathcad

Как известно, многие уравнения не имеют аналитических решений. В первую очередь это относится к большинству трансцендентных уравнений. Такие уравнения могут решаться численными методами с заданной точностью (не более значения заданного системной переменной TOL).

Для уравнений вида решение в Mathcad находится с помощью функции

Возвращает значение , принадлежащее отрезку , при котором выражение или функция обращается в 0. Оба аргумента этой функции должны быть скалярами. Функция возвращает скаляр.

Аргументы:

функция, определенная где-либо в рабочем документе, или выражение. Выражение должно возвращать скалярные значения.

имя переменной, которая используется в выражении. Этой переменной перед использованием функции необходимо присвоить числовое значение. Mathcad использует его как начальное приближение при поиске корня.

– необязательны, если используются, то должны быть вещественными числами, причем .

Наиболее распространен графический метод определения начальных приближений. Принимая во внимание, что действительные корни уравнения – это точки пересечения графика функции с осью абсцисс, достаточно построить график функции и отметить точки пересечения с осью , или отметить на оси отрезки, содержащие по одному корню. Построение графиков часто удается сильно упростить, заменив уравнение равносильным ему уравнением:

где функции и – более простые, чем функция . Тогда, построив графики функций и , искомые корни получим как абсциссы точек пересечения этих графиков.

Отсутствие сходимости функции

Если после многих итераций Mathcad не находит подходящего приближения, то появится сообщение (отсутствует сходимость). Эта ошибка может быть вызвана следующими причинами:

— уравнение не имеет корней;

— корни уравнения расположены далеко от начального приближения;

— выражение имеет локальные и между начальным приближением и корнями;

— выражение имеет разрывы между начальными приближениями и корнями;

— выражение имеет комплексный корень, но начальное приближение было вещественным.

Чтобы установить причину ошибки, исследуйте график . Он поможет выяснить наличие корней уравнения и, если они есть, то определить приблизительно их значения. Чем точнее выбрано начальное приближение корня, тем быстрее будет сходиться.

Рекомендации по использованию функции

Для изменения точности, с которой функция ищет корень, нужно изменить значение системной переменной TOL. Если значение TOL увеличивается, функция будет сходиться быстрее, но ответ будет менее точен. Если значение TOL уменьшается, то функция будет сходиться медленнее, но ответ будет более точен. Чтобы изменить значение TOL в определенной точке рабочего документа, используйте определение вида . Чтобы изменить значение TOL для всего рабочего документа, измените его настройки.

Если два корня расположены близко друг от друга, следует уменьшить TOL, чтобы различить их.

Теорема 2. Итерационный процесс половинного деления сходится к искомому корню ξ с любой наперед заданной точностью ε.
Доказательство: Рассмотрим последовательность чисел ξi являющихся приближением корня на i -ом шаге.
ξi=½(bi+ai), i=0,1.
где a=a; b=b; ai;bi — границы подынтервалов, в которых f(ai)f(bi) 0 мы ни задали, всегда можно найти такое n , что ч.т.д.
Графически метод дихотомии выглядит следующим образом

|f(c)|≤δ f(a)f(c) 10 = 1024 ≈ 10 3 раз. За 20 итераций (n=2) уменьшается в 2 20 ≈ 10 6 раз.

Пример №1 . Найти экстремум функции: y=5x 2 -4x+1 методом дихотомии, если ε=0.1, а исходный интервал [0,10].

  • Решение
  • Видео решение

Пример №3 . Методом бисекции найти решение нелинейного уравнения на отрезке [a,b] с точностью ε = 10 -2 . Выбрав полученное решение в качестве начального приближения, найти решение уравнения методом простой итерации с точностью ε = 10 -4 . Для метода простой итерации обосновать сходимость и оценить достаточное для достижения заданной точности число итераций.
sqrt(t)+x 2 = 10, a = 2.6, b = 3

Найдем корни уравнения:
Используем для этого Метод половинного деления (метод дихотомии)..
Считаем, что отделение корней произведено и на интервале [a,b] расположен один корень, который необходимо уточнить с погрешностью ε.
Итак, имеем f(a)f(b) 1 /2(a+b) и вычисляем f(c). Проверяем следующие условия:
1. Если |f(c)| 1 /2 n (b-a)
В качестве корня ξ. возьмем 1 /2(an+bn). Тогда погрешность определения корня будет равна (bn – an)/2. Если выполняется условие:
(bn – an)/2 1 /2(an+bn).
Решение.
Поскольку F(2.6)*F(3) 0, то a=2.8
Итерация 2.
Находим середину отрезка: c = (2.8 + 3)/2 = 2.9
F(x) = 0.113
F(c) = -0.487
Поскольку F(c)•F(x) 0, то a=2.825
Остальные расчеты сведем в таблицу.

Название: Решение нелинейных уравнений
Раздел: Рефераты по математике
Тип: контрольная работа Добавлен 09:13:39 31 декабря 2010 Похожие работы
Просмотров: 3535 Комментариев: 8 Оценило: 7 человек Средний балл: 4.1 Оценка: 4 Скачать
Читайте также:  Как добавить свою музыку в мафию 2
N c a b f(c) f(x)
1 2.6 3 2.8 -1.6275 -0.4867
2 2.8 3 2.9 -0.4867 0.1129
3 2.8 2.9 2.85 0.1129 -0.1893
4 2.8 2.85 2.825 -0.1893 -0.3386
5 2.825 2.85 2.8375 -0.3386 -0.2641
6 2.8375 2.85 2.8438 -0.2641 -0.2267

Ответ: x = 2.8438; F(x) = -0.2267
Решение было получено и оформлено с помощью сервиса Метод Ньютона онлайн

Пример №2 . Локализовать корень нелинейного уравнения f(x) = 0 и найти его методом бисекции с точностью ε1 = 0,01. Выбрав полученное решение в качестве начального приближения, найти решение уравнения методом простой итерации с точностью ε2 = 0,0001. Для метода простой итерации обосновать сходимость и оценить достаточное для достижения заданной точности ε2 число итераций.

Ссылка на основную публикацию
Лучшие производители шкафов купе в москве рейтинг
*Обзор лучших по мнению редакции expertology.ru. О критериях отбора. Данный материал носит субъективный характер, не является рекламой и не служит...
Куда устанавливается visio viewer
Файлы с расширением .vsd и .vsdx предназначены для программы Microsoft Visio Viewer (vviewer.exe) и содержат визуальные объекты, текст, диаграммы и...
Кью ар код программа
Что такое QR-код: QR код «QR - Quick Response - Быстрый Отклик» — это двухмерный штрихкод (бар-код), предоставляющий информацию для...
Магазины дискаунтеры в россии список
Как закончили 2018 год десять крупнейших розничных ретейлеров РФ? Какая сеть растет быстрее всех? Почему гипермаркеты умирают, а «жесткие дискаунтеры»...
Adblock detector