Оценка точности вычисления «неберущихся» интегралов




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

Определенные интегралы от таких функций можно вычислить только приближенно. Для оценки точности вычисления в таких случаях используют, например, правило Рунге. В данном случае интеграл вычисляется по выбранной формуле (прямоугольников, трапеций, парабол Симпсона) при числе шагов, равном n, а затем при числе шагов, равном . Погрешность вычисления значения интеграла при числе шагов, равном , вычисляется по формуле Рунге: , для формул прямоугольников и трапеций , а для формулы Сипсона . Таким образом, интеграл вычисляется для последовательных значений числа шагов , ,..., где – начальное число шагов. Процесс вычислений заканчивается, когда для очередного значения будет выполнено условие , где – заданная точность.

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

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

Квадратурная формула прямоугольников.

Вычислим, при каком шаге погрешность будет составлять 0,01:

.

Поскольку , то .

При шаге отрезок разбивается на равностоящих узлов.

Квадратурная формула трапеций.

Вычислим, при каком шаге погрешность составит 0,01:

.

Поскольку , .

При шаге ,отрезок разбивается на равностоящих узлов.

Квадратурная формула Симпсона.

Вычислим, при каком шаге погрешность составит 0,01:

,

, .

 

При шаге , отрезок разбивается на равностоящих узлов.

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

 


Содержание РГР «Приближенные методы вычисления

Определенных интегралов»

Студенту предлагается работа, состоящая из четырех этапов:

1 этап – точное вычисление определенного интеграла.

2 этап – приближенное вычисление определенного интеграла одним из методов: прямоугольников или трапеций.

3 этап – приближенное вычисление определенного интеграла методом парабол.

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

Построение графика подынтегральной функции.

 

Варианты и образец выполнения РГР приведены ниже.

Варианты

№ варианта f(x) a b Шаг h
      0,1
      0,1
      0,1
      0,1
    π 0,1π
      0,1
      0,1
      0,1
      0,1
    π/2 0,05π
      0,1
      0,1
      0,1
      0,1
    π 0,1π
      0,1
      0,1
    π/2 0,05π
      0,1
    π/2 0,05π
      0,1
      0,1
    π/2 0,05π
    π/2 0,05π
    π/2 0,05π
    π/2 0,05π
      0,1
      0,1
    π/2 0,05π
      0,1

 


Образец выполнения РГР

Задание. Вычислить интеграл

1. Точное вычисление:

= 0,40631714.

 

2. Приближенное вычисление с помощью формул прямоугольников:

,

, .

, .

Составим таблицу:

xi yi = f (xi)
     
  0,1 0,010005
  0,2 0,04016
  0,3 0,091207
  0,4 0,165041
  0,5 0,265165
  0,6 0,396981
  0,7 0,567851
  0,8 0,786966
  0,9 1,065081
    1,414214

 

По первой формуле прямоугольников получаем:

≈ 0,1 = 0,1·3,062514 = 0,306251.

По второй формуле прямоугольников получаем:

≈ 0,1 = 0,1· 4,802669 = 0,480267.

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

Вычислим относительную и абсолютную погрешности.

I = 0,40631714, ,

, .

, .

3. Приближенное вычисление по формуле трапеций:

В нашем случае получаем:

≈ 0,1 = =0,1 = 0,1·4,095562 = =0,409556.

Вычислим относительную и абсолютную погрешности.

I = 0,40631714,

, .

4. Приближенное вычисление по формуле Симпсона:

 

 

В нашем случае получаем:

 

=

= = 0,406325.

 

Вычислим относительную и абсолютную погрешности.

I = 0,40631714

, .

 

В действительности, = 0,40631714.

Таким образом, при разбиении отрезка на 10 частей по формуле Симпсона мы получили 5 верных знаков; по формуле трапеций – три верных знака; по формуле прямоугольников мы можем ручаться только за первый знак.


 

Литература:

1. Задачи и упражнения по математическому анализу для втузов: учеб. пособие для втузов / Г. С. Бараненков [и др.]; под ред. Б.П. Демидовича. - М.; Владимир: Астрель: Изд-во АСТ: ВКТ, 2010. - 495 с.

2. Кудрявцев, Л.Д. Краткий курс математического анализа: учеб. для вузов. Т. 1: Дифференциальное и интегральное исчисления функций одной переменной. Ряды / Л. Д. Кудрявцев. - 3-е изд., перераб. - М.: ФИЗМАТЛИТ, 2009. - 399 с.

3. Пискунов, Н.С. Дифференциальное и интегральное исчисления: учеб. пособие для втузов. В 2 т. Т. 1 / Н. С. Пискунов. - Изд. стер. - М.: Интеграл-Пресс, 2004. - 415 с.

4. Шипачев В.С. Высшая математика: учеб. для вузов / В. С. Шипачев. - 7-е изд., стер. - М.: Высш. шк., 2005. - 479 с.

 


 

 

Баранова И.М., Гущин Г.В., Еловиков А.Б., Часова Н.А.

 

МАТЕМАТИКА

 

ПРИБЛИЖЕННЫЕ МЕТОДЫВЫЧИСЛЕНИЯ

ОПРЕДЕЛЕННЫХ ИНТЕГРАЛОВ

 

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

 

 

Формат Объем Тираж Заказ

 

 

Брянск, Станке Димитрова 3, Редакционно-издательский отдел

Отпечатано: Печатный цех БГИТА

 

 



Поделиться:




Поиск по сайту

©2015-2024 poisk-ru.ru
Все права принадлежать их авторам. Данный сайт не претендует на авторства, а предоставляет бесплатное использование.
Дата создания страницы: 2018-01-08 Нарушение авторских прав и Нарушение персональных данных


Поиск по сайту: