Отыскание корня уравнения методом половинного деления

Содержание

 

1. Индивидуальное задание

2. Постановка задачи и формализация

3. Выбор, обоснование, краткое описание методов

3.1 Численное интегрирование

3.1.1 Постановка задачи

3.1.2 Выбор и описание метода

3.2 Отыскание корня уравнения

3.2.1 Постановка задачи

3.2.2 Выбор и описание метода (половинное деление)

4. Проверка условий сходимости методов

5. Тестирование программных модулей

5.1 Тестирование модуля численного интегрирования

5.1.1 Схема алгоритма тестирующей программы

5.1.2 Код тестирующей программы

5.1.3 Результат тестирования

5.2 Тестирование модуля отыскания корня уравнения методом половинного деления

5.2.1 Схема алгоритма тестирующей программы

5.2.2 Код тестирующей программы

5.2.3 Результат тестирования

5.3 Прогонка программы

5.3.1 Схема алгоритма программы при прогонке

5.3.2 Код программы при прогонке

5.3.3 Результаты работы программы при прогонке

6. Детализированная схема алгоритма

7. Код программы

8. Полученные результаты

9. Проверка результатов в MathCAD

10. Основные выводы

Список литературы

модуль корень половинный деление

1. Индивидуальное задание

 

Решить уравнение  на отрезке x℮[0;2р]

2. Постановка задачи и формализация

Задача заключается в поиске корня уравнения f(x)=0 численным методом на отрезке неопределённости [0; 2р], где

Интегрирование проводится численным методом.

Для решения поставленной задачи необходимо разработать следующие модули:

- главный модуль, вводящий исходные данные (требуемую точность и концы отрезка неопределённости) и выводящий конечный результат (решение уравнения)

- модуль, задающий подынтегральное выражение

- модуль, выполняющий численное интегрирование и вычитающий р/2

- модуль, решающий нелинейное уравнение f(x)=0, где f(x) – значение функции, полученное в предыдущем модуле

Укрупнённый алгоритм решения задачи:

3. Выбор, обоснование, краткое описание методов

3.1 Численное интегрирование

 

3.1.1 Постановка задачи

Если функция f(x) непрерывна на отрезке [a;b] и дифференцируема, то определённый интеграл от этой функции в пределах от a до b существует и может быть вычислен по формуле Ньютона-Лейбница:

Причём

Задача численного интегрирования заключается в нахождении значения определённого интеграла через ряд значений подынтегральной функции y_i=f(x_i), заданной в точках x_i (i=0,1,…,n), причём x₀=a, x_n=b. Чащё всего интервал разбивают на подынтервалы длиной h=x_i+1 – x_i

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

Замена функции f(x) интерполирующим полиномом приводит к образованию погрешности вычисления значения интеграла

Здесь I₁ – точное значение интеграла, I – значение, вычисленное численным методом, R- погрешность расчёта численным методом.

3.1.2 Выбор и описание метода

Выбор метода:

Находить значение интеграла можно многими способами, среди которых:

1)  формула прямоугольников

2)  формула трапеций

3)  формула Симпсона

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

Наибольшую точность, т.к. подынтегральная функция аппроксимируется полиномом 2 порядка.

Описание метода:

Если для каждой пары отрезков [x_i;x_i+2] построить многочлен второй степени, затем проинтегрировать его и воспользоваться свойством аддитивности интеграла, то получим формулу Симпсона:

 n=2*m – чётное число

Геометрическая интерпретация формулы Симпсона:

На отрезке [x_i;x_i+2] длиной 2h строится парабола, проходящая через три точки (x_i;y_i), (x_i+1;y_i+1), (x_i+2;y_i+2). Площадь под параболой, заключённой между осью абсцисс и прямыми x=x_i, x=x_i+2, принимают равной интегралу

3.2 Поиск корня нелинейного уравнения

 

3.2.1 Постановка задачи

Пусть требуется найти решение уравнения f(x)=0. f(x) – непрерывная функция в конечном или бесконечном интервале. Если f(x) представляет собой многочлен, то уравнение называют алгебраическим, в противном случае – трансцендентным.

Всякое значение x=x^*, обращающее f(x) в ноль, называется корнем этого уравнения.

Решение задачи отыскания изолированных корней состоит из двух этапов: отделение корней, уточнение корней. При отыскании действительных корней этап отделения производится либо графически, либо аналитически, основываясь на теореме: если f(x) принимает на разных концах отрезка [a;b] разные знаки, то на [a;b] существует по меньшей мере один корень уравнения f(x)=0.

Корень будет единственным на отрезке [a;b], если производная f(x) существует и сохраняет знак внутри [a;b].