Метод деления отрезка пополам — КиберПедия 

Индивидуальные очистные сооружения: К классу индивидуальных очистных сооружений относят сооружения, пропускная способность которых...

Эмиссия газов от очистных сооружений канализации: В последние годы внимание мирового сообщества сосредоточено на экологических проблемах...

Метод деления отрезка пополам

2017-11-22 332
Метод деления отрезка пополам 0.00 из 5.00 0 оценок
Заказать работу

В отличие от предыдущего метода метод гарантированно дает корень, если он имеется на заданном отрезке [ a, b ]. Суть метода показана на рис. 1.68. В качестве исходных данных задается отрезок [ a, b ] с корнем внутри его. Вычисляется функция fa = f (a). Отрезок [ a, b ] делится пополам точкой c = (a + b) / 2. Вычисляется функция f c = f (c). Далее проверяются знаки fa и f c. Если эти величины имеют одинаковые знаки, то точка a переносится в точку c:

a = с, fa = fc

.Если f a и fc имеют разные знаки, то точка b переносится в точку с: b = c. Далее процесс сужения отрезка [ a, b ] продолжается до тех пор, пока b – a не станет меньше заданной погрешности. Корень при этом в точке с.

Рис. 1.68. Суть метода деления отрезка пополам

Метод Ньютон

Суть метода показана на рис. 1.69. В качестве начального приближения в методе задается значение x 0. Метод используется тогда, когда производная функции f (x) определена в аналитическом виде.

В точке x 0 вычисляется значение функции f 0 = f (x 0). Далее по формуле

x 1 = x 0f 0 / f ’(x 0) (2)

вычисляется точка x 1 пересечения касательной, проведенной из точки f 0. Если | x 1 x 0|>E (E —заданная погрешность), то полагается x 0 = x 1, вычисляется f 0 = f (x 0) и итерационный процесс по формуле (2) продолжается до тех пор, пока не окажется |x1 – x0| < E. Корень при этом в x 1. Метод Ньютона не гарантирует сходимость алгоритма к корню.

Рис. 1.69. Суть метода ньютона

Метод Вегстейна

Метод является модификацией метода Ньютона и используется для функций, у которых неизвестна в аналитическом виде производная f ’(x 0). Суть метода показана на рис. 1.70.

В качестве начального приближения в методе задается два начальных приближения x 0 и x 1. Желательно эти начальные приближения задавать по разную сторону от корня. Далее вычисляются значения функции f 0 и f 1 и через точки f 0 и f 1 проводится секущая до пересечения ее с осью Х в точке x 2:

x 2 = x 0 – (f 0 (x 0x 1) / (f 0f 1)) (3)

Далее вычисляется значение функции f 2. Если |x 2x 0| > E (E —погрешность), то полагается

x 0 = x 1 , x 1 = x 2, f 0 = f 1, f 1 = f 2

и итерационной процесс по формуле (3) продолжается до тех пор, пока не окажется | x 2x 0| < E. Метод Вегстейна практически всегда сходится к корню.

Рис. 1.70. Суть метода Вегстейна

Численное вычисление определенного интеграла

Во всех методах исходный отрезок [ a, b ] разбивается на более мелкие интервалы длиной h. Величина шага h определяется точностью вычисления интеграла.

Метод средних

Суть метода показана на рис. 1.71.

В этом методе площадь, ограниченная кривой f (x) заменяется площадью прямоугольника fi * h, где fi = f (xi + h / 2). В общем случае интеграл по формуле средних определяется формулой:

Рис. 1.71. Суть метода средних

Метод трапеций

Суть метода показана на рис. 1.72.

В этом методе площадь, ограниченная кривой f (x) заменяется площадью трапеции . В общем случае интеграл по формуле трапеций определяется по формуле:

Рис. 1.72. Суть метода трапеций

Метод Симпсона

Суть метода показана на рис. 1.73.

В этом методе площадь, ограниченная кривой f (x) заменяется площадью, ограниченной параболой, где xi+1/2=(xi + xi +1) / 2. В общем случае интеграл по формуле Симпсона определяется формулой:

Рис. 1.73. Суть метода Симпсона

Аппроксимация функций

В качестве исходных данных аппроксимации задается таблично заданная функция fi для дискретного набора аргументов xi, . При аппроксимации известна кривая, которая должна покрывать набор табличных значений. У этой кривой f (a, x) неизвестны только параметры a, которые обеспечивают оптимальное расположение кривой f (a, x) относительно табличных данных. Подбор параметров а осуществляется по методу наименьших квадратов, в соответствии с которым:

Минимум этой функции достигается, если производные по всем параметрам аj вектора а равны нулю:

Очень часто используется линейная аппроксимация:

В этом случае нелинейная система уравнений преобразуется к линейной

где:

Для решения данной задачи необходимо решить систему линейных алгебраических уровнений относительно aj .В случае представления f (a, x) набором ортогональных функций матрица сводится к диагональному виду, так как =0 для ik. Для ортогональных величина ak определяется выражением:

Очень часто аппроксимируемая функция f (a, x) называется линией тренда. В частности, в Excel существует сервис получения линий тренда для фиксированного набора функций f (a, x). Приведем пример получениялинии тренда в среде Excel. Прежде всего необходимо задать табличные значения функции с х в 1-й строке и f во второй строке (рис. 1.74).

Рис. 1.74. Исходная функциональная зависимость

По этим даны далее необходимо построить точечный график. Для этого выделяем исходные данные (рис. 1.74) в выбираем меню Вставка | Диаграмма. В процессе диалога Мастера (рис. 1.75) выбираем точечную диаграмму и далее (рис. 1.76) заполняем поля надписей диаграммы.

Рис. 1.75. Окно Мастера для выбора типа диаграммы

В конце диалогового режима с Мастером будет построен точечный график (рис. 1.77).

Рис. 1.76. Заполнение полей надписей диаграммы

Рис. 1.77. Построенный точечный график

Для построения линии тренда необходимо подвести указатель мыши на любую точку графика, что бы высветилась всплывающая подсказка (рис. 1.78) и далее правой клавиши мыши вывести контекстное меню (рис. 1.79). В контекстном меню выбрать пункт Добавить линию тренда. Из предложенных вариантов линий тренда выбрать полиномиальную степени 2 (рис. 1.80). В результате на имеющемся уже точечном графике будет построена оптимальная для этих точек линия тренда (рис. 1.81).

Рис. 1.78. Всплывающая подсказка

Рис. 1.79. Контекстное меню

Рис. 1.80. Из предложенных линий выбираем полиномиальную

Рис. 1.81. Построенная линия тренда

Интерполяция функции

В качестве исходных данных при интерполяции задается таблица (fi, xi). Задачей интерполяции является восстановление значений функции при любом значении x по узлам табличных значений. Наиболее простым методом является линейная интерполяция, суть которой демонстрирует рис. 1.82.

Рис. 1.82. Суть линейной интерполяции

Несложно показать геометрически, что при линейной интерполяции значение f в точке x определяется выражением:

при условии, что xixxi +1.

Более точной является квадратичная интерполяция, суть которой продемонстрирована на рис. 1.83.

Рис. 1.83. Суть квадратичной интерполяции

Несложно показать геометрически, что при квадратичной интерполяции значение f в точке x определяется выражением:

при условии, что | x-xi |=min для всех .


Поделиться с друзьями:

Автоматическое растормаживание колес: Тормозные устройства колес предназначены для уменьше­ния длины пробега и улучшения маневрирования ВС при...

История создания датчика движения: Первый прибор для обнаружения движения был изобретен немецким физиком Генрихом Герцем...

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

Таксономические единицы (категории) растений: Каждая система классификации состоит из определённых соподчиненных друг другу...



© cyberpedia.su 2017-2024 - Не является автором материалов. Исключительное право сохранено за автором текста.
Если вы не хотите, чтобы данный материал был у нас на сайте, перейдите по ссылке: Нарушение авторских прав. Мы поможем в написании вашей работы!

0.027 с.