Главная » Курсовые » Численные методы

Лабараторная работа 2. Интегральное исчисление

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

 

  Если функции u = j(x) и v = y(x) непрерывны на отрезке [a, b], а также непрерывны на этом отрезке их производные, то справедлива формула интегрирования по частям:

Приближенное вычисление определенного интеграла.

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

Формула прямоугольников.

   Если известны значения функции f(x) в некоторых точках x0, x1, … , xm, то в качестве функции “близкой” к f(x) можно взять многочлен Р(х) степени не выше m, значения которого в выбранных точках равны значениям функции f(x) в этих точках.

Если разбить отрезок интегрирования на n равных частей . При этом:

y0 = f(x0), y1 = f(x1), …. , yn = f(xn).

Составим суммы: y0Dx + y1Dx + … + yn-1Dx

 y1Dx + y2Dx + … + ynDx

Это соответственно нижняя и верхняя интегральные суммы. Первая соответствует вписанной ломаной, вторая – описанной.

  Тогда  или

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

 

Формула трапеций.

 

  Эта формула является более точной по

 у сравнению с формулой прямоугольников.

  Подинтегральная функция в этом случае

  заменяется на вписанную ломаную.

 

  Геометрически площадь криволинейной трапеции заменяется суммой площадей вписанных трапеций. Очевидно, что чем больше взять точек n разбиения интервала, тем с большей точностью будет вычислен интеграл.

 

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

После приведения подобных слагаемых получаем формулу трапеций:

 

Формула парабол (формула Симпсона или квадратурная формула).

 (Томас Симпсон (1710-1761)- английский математик)

  Разделим отрезок интегрирования [a, b] на четное число отрезков (2m). Площадь криволинейной трапеции, ограниченной графиком функции f(x) заменим на площадь криволинейной трапеции, ограниченной параболой второй степени с осью симметрии, параллельной оси Оу и проходящей через точки кривой, со значениями f(x0), f(x1), f(x2).

  Для каждой пары отрезков построим такую параболу.

 у

   Уравнения этих парабол имеют вид Ax2 + Bx + C, где коэффициенты А, В, С могут быть легко найдены по трем точкам пересечения параболы с исходной кривой.

   (1)

Обозначим .

Если принять х0 = -h, x1 = 0, x2 = h, то  (2)

Тогда уравнения значений функции (1) имеют вид:

C учетом этого: .

Отсюда уравнение (2) примет вид: 

Тогда

 

  Складывая эти выражения, получаем формулу Симпсона:

Чем больше взять число m, тем более точное значение интеграла будет получено.

 

  Пример. Вычислить приближенное значение определенного интеграла

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

 

  По формуле Симпсона получим:

  m

  0

  1

  2

  3

  4

  5

  6

  7

  8

  9

  10

  x

  -2

  -1

  0

  1

  2

  3

  4

  5

  6

  7

  8

  f(x)

2.828

3.873

  4

4.123

4.899

6.557

8.944

11.874

15.232

18.947

22.978

 

 

  Точное значение этого интеграла – 91.173.

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

 

  Для сравнения применим к этой же задаче формулу трапеций.

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

 

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

  Принцип этого метода состоит в том, чтобы заменить подинтегральную функцию по формуле Тейлора и почленно проинтегрировать полученную сумму.

 

Пример. С точностью до 0,001 вычислить интеграл

 

 

Т.к. интегрирование производится в окрестности точки х=0, то можно воспользоваться для разложения подинтегральной функции формулой Маклорена.

  Разложение функции cosx имеет вид:

 

 

Зная разложение функции cosх легко найти функцию 1 – cosx:

 

В этой формуле суммирование производится по п от 1 до бесконечности, а в предыдущей – от 0 до бесконечности. Это – не ошибка, так получается в результате преобразования.

 

Теперь представим в виде ряда подинтегральное выражение.

Теперь представим наш интеграл в виде:

В следующем действии будет применена теорема о почленном интегрировании ряда. (Т.е. интеграл от суммы будет представлен в виде суммы интегралов членов ряда).

Вообще говоря, со строго теоретической точки зрения для применения этой теоремы надо доказать, что ряд сходится и, более того, сходится равномерно на отрезке интегрирования [0, 0,5]. Эти вопросы будут подробно рассмотрены позже (См. Действия со степенными рядами.) Отметим лишь, что в нашем случае подобное действие справедливо хотя бы по свойствам определенного интеграла (интеграл от суммы равен сумме интегралов).

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

Для справки: Точное (вернее – более точное) значение этого интеграла: 0,2482725418…

levitra paypal kaufen viagra spanien bestellen what is the cost of generic viagra in canada cialis generique canada
levitra paypal kaufen viagra spanien bestellen what is the cost of generic viagra in canada cialis generique canada