Методы численного интегрирования (нахождение определенных интегралов)

Содержание

Слайд 2

1. Аналитический метод Где F(x) – первообразная функции f(x)

1. Аналитический метод


Где F(x) – первообразная функции f(x)

Слайд 3

Аналитический метод интегрирования не всегда может быть применен на практике. Пример «неберущегося» интеграла:

Аналитический метод интегрирования не всегда может быть применен на практике.
Пример «неберущегося»

интеграла:
Слайд 4

Графическая интерпретация определенного интеграла Линии ограничения: y=0; y=f(x); x=a; x=b.

Графическая интерпретация определенного интеграла

Линии ограничения:
y=0;
y=f(x);
x=a;
x=b.

Слайд 5

2. Численные методы 1. [a, b] разбивается на n равных отрезков

2. Численные методы


1. [a, b] разбивается на n равных отрезков

длиной h.
2. Площадь S разбивается на n полос шириной h.
3. Полоса представляется в виде геометрической фигуры.
4. Рассчитывается площадь каждой полосы.
5. Искомый интеграл есть сумма площадей всех полос.
Слайд 6

1. Метод прямоугольников Отдельно взятая полоса представляется в виде прямоугольника шириной

1. Метод прямоугольников
Отдельно взятая полоса представляется в виде прямоугольника шириной h.
ВОПРОС:

Какая величина принимается за высоту прямоугольника?
Слайд 7

А. Метод левых прямоугольников Высота - значение функции в левой точке

А. Метод левых прямоугольников

Высота - значение функции в левой точке основания

каждой полосы.
Формула расчета интеграла:
Слайд 8

B. Метод правых прямоугольников Высота - значение функции в правой точке

B. Метод правых прямоугольников

Высота - значение функции в правой точке основания

каждой полосы.
Формула расчета интеграла:
Слайд 9

С. Метод средних прямоугольников Высота - значение функции в середине основания каждой полосы. Формула расчета интеграла:

С. Метод средних прямоугольников

Высота - значение функции в середине основания каждой

полосы.
Формула расчета интеграла:
Слайд 10

Блок-схема метода средних прямоугольников

Блок-схема метода средних прямоугольников

Слайд 11

2. Метод трапеций Отдельно взятая полоса представляется в виде перевернутой трапеции

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

Отдельно взятая полоса представляется в виде перевернутой трапеции высотой

h.
Основания трапеции будут равны значениям функции в левой и правой точке высоты трапеции.
Площадь трапеции:
Слайд 12

Гладкая кривая заменяется ломаной линией Интеграл рассчитывается по следующей формуле:

Гладкая кривая заменяется ломаной линией

Интеграл рассчитывается по следующей формуле:

Слайд 13

Блок-схема метода трапеций

Блок-схема метода трапеций

Слайд 14

3. Метод Симпсона Гладкая функция заменяется участками парабол. Через любые 3

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

Гладкая функция заменяется участками парабол.
Через любые 3 точки

на плоскости можно провести одну и только одну параболу.
Парабола проводится через точки пересечения границ 2-х соседних полос с графиком подынтегральной функции.
Слайд 15

Гладкая кривая заменяется участками парабол Каждая парабола заменяет исходную подынтегральную функцию

Гладкая кривая заменяется участками парабол

Каждая парабола заменяет исходную подынтегральную функцию

сразу над двумя полосами. Следовательно, число разбиений должно быть четным !!!
Рассмотрим ситуацию с одной параболой (2-мя полосами) и выведем формулу для расчета интеграла.
Слайд 16

Любая парабола описывается уравнением: y=ax2+bx+c Точки (0, y0), (h, y1), (2h,

Любая парабола описывается уравнением:
y=ax2+bx+c
Точки (0, y0), (h, y1), (2h, y2) лежат

на одной параболе, следовательно, должны удовлетворять одной и той же функции.
Число разбиений должно быть четным !!!
Слайд 17

Подставляем координаты 3-х точек в уравнение для параболы, получаем систему линейных

Подставляем координаты 3-х точек в уравнение для параболы, получаем систему линейных

алгебраических уравнений.
Здесь неизвестные - параметры параболы: a, b, c.
Из 1-го уравнения: y0=c.
Произведя замену, получим новую систему уравнений:
Решаем полученную СЛАУ методом Крамера:
Слайд 18

Выведем формулу для расчета коэффициентов a и b:


Выведем формулу для расчета коэффициентов a и b:

Слайд 19

Площадь под фигуры можно вычислить, проинтегрировав полученную параболическую зависимость: y=ax2+bx+c Подставим

Площадь под фигуры можно вычислить, проинтегрировав полученную параболическую зависимость:
y=ax2+bx+c

Подставим в полученную

формулу значения для коэффициентов параболы a, b и c:
Слайд 20

Получим: Если число разбиений будет не 2, а 4, то формула

Получим:
Если число разбиений будет не 2, а 4, то формула для

вычисления интеграла будет иметь следующий вид:
Слайд 21

В общем виде: Формула Симпсона

В общем виде:
Формула Симпсона

Слайд 22

Блок-схема метода Симпсона

Блок-схема метода Симпсона

Слайд 23

Замечания о погрешности численного интегрирования

Замечания о погрешности численного интегрирования

Слайд 24

Для оценки погрешности численного интегрирования сравним значения интеграла, рассчитанные различными численными

Для оценки погрешности численного интегрирования сравним значения интеграла, рассчитанные различными численными

методами с истинным значением интеграла, рассчитанным аналитически.
Пример:
Истинное значение: S=5
Слайд 25

Слайд 26

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

Из таблицы видно, что погрешность зависит от метода интегрирования и от

количества разбиений интервала интегрирования.