Метод симпсона численного интегрирования. Как вычислить определенный интеграл по формуле Симпсона? Примеры приближенного вычисления определенных интегралов методом парабол

Суть метода Симпсона заключается в приближении подынтегральной функции на отрезке интерполяционным многочленом второй степени p2(x), т.е. приближение графика функции на отрезке параболой. Для интерполирования подынтегральной функции используются три точки.

Рассмотрим произвольный интеграл. Воспользуемся заменой переменной таким образом, чтобы границы отрезка интегрирования вместо стали [-1,1]. Для этого введем переменную z:

Рассмотрим задачу интерполирования подынтегральной функции, используя в качестве узлов три равноудаленные узловые точки z = -1, z = 0, z = +1 (шаг равен 1, длина отрезка интегрирования равна 2). Обозначим соответствующие значения подынтегральной функции в узлах интерполяции:

Система уравнений для нахождения коэффициентов полинома, проходящего через три точки (-1, f-1), (0, f0) и(1, f-+1) примет вид:

Коэффициенты легко могут быть получены:

Вычислим теперь значение интеграла от интерполяционного многочлена:

Путем обратной замены переменной вернемся к исходному интегралу. Учтем, что:

соответствует

Получим формулу Симпсона для произвольного интервала интегрирования:

Полученное значение совпадает с площадью криволинейной трапеции, ограниченной осью x, прямыми x = x0, x = x2 и параболой, проходящей через точки

При необходимости, исходный отрезок интегрирования может быть разбит на N сдвоенных отрезков, к каждому из которых применяется формула Симпсона. Шаг интерполирования при этом составит:

Для первого отрезка интегрирования узлами интерполирования будут являться точки a, a+h, a+2h, для второго a+2h, a+3h, a+4h, третьего a+4h, a+5h, a+6h и т.д. Приближенное значение интеграла получается суммированием N площадей:

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

В данную сумму входят одинаковые слагаемые (для внутренних узлов с четным значением индекса - 2i). Поэтому можно перегруппировать слагаемые в этой сумме таким образом:

Приняв во внимание то, что получаем:

Оценим теперь погрешность интегрирования по формуле Симпсона. Будем считать, что у функции на отрезкесуществуют непрерывные производные. Составим разность:

Применяя к этой разнице последовательно теорему о среднем и дифференцируя R(h) получаем погрешность метода Симпсона:

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

Преимущества и недостатки

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

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

Текст работы размещён без изображений и формул.
Полная версия работы доступна во вкладке "Файлы работы" в формате PDF

Введение

Уже в 10 классе я задумываюсь о том, мне нужно будет сдавать профильный ЕГЭ по математике. Решая задания ЕГЭ, я столкнулся с заданиями на нахождение объема многогранников и тел вращения, хотя это задания из программы 11 класса. Заинтересовавшись этим вопросом, я узнал, что в связи с многообразием геометрических фигур тел существует огромное количество формул для нахождения площадей и объёма (на каждуюфигуруи каждое тело приходится своя формула). Рассматривая формулы по геометрии, я убедился, что огромное количество формул связано с площадями и объемами фигур. Таких формул более двенадцати по площадям плоских фигур и более десяти по объемам пространственных тел.

И я задался вопросом : а существует ли такая универсальная формула для нахождения площади и объёма геометрических фигур и тел?

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

Проблема

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

Гипотеза

В XYIII веке английский математик Томас Симпсон вывел формулу для нахождения некоторых площадей плоских фигур и объемов пространственных тел через вычисление площадей нижнего, верхнего и среднего основания.

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

Цель работы: доказать, что универсальная формула Симпсона может заменить все изучаемые формулы площадей и объемов в школьном курсе геометрии и ей можно пользоваться не только на практике, но и на экзаменах, в том числе и на ЕГЭ.

Задачи работы:

Изучить основные характеристики геометрических тел стереометрии: призмы, пирамиды, конуса, цилиндра, шара;

Изучить имеющуюся литературу по данной теме.

Используя универсальную формулу, вывести формулы площадей и объемов для всех фигур и тел.

Сравнить полученные формулы с формулами, предлагаемыми в учебнике.

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

Практическая значимость моей работы: Результаты данной работы могут иметь применение в школьной практике, а именно использоваться на занятиях по геометрии и алгебре, при подготовке и сдаче ЕГЭ.

Глава 1 Краткие характеристики свойств геометрических тел

Школьный курс геометрии делится на планиметрию и стереометрию. С 7 по 9 класс я изучал свойства фигур на плоскости, в том числе и формулы для нахождения их площадей (Приложение 1-2).

В курсе 10 класса я начал изучать раздел геометрии-стереометрия, в котором изучаются свойства фигур в пространстве. При написании работы, я рассмотрел геометрические тела и их поверхности. Объёмные геометрические тела делятся на многогранники и тела вращения.

Многогранник - поверхность, составленная из многоугольников и ограничивающая некоторое геометрическое тело.

Тела вращения - геометрические тела, полученные путём вращения вокруг своей оси. Тела вращения: цилиндр, конус, шар.

Многогранники бывают выпуклые и невыпуклые. Выпуклые многогранники - расположены по одну сторону от плоскости каждой грани. Невыпуклые многогранники - расположены по обе стороны от плоскости хотя бы одной грани.

Пирамида

Параллелепипед

Глава 2. Формула Симпсона

Томас Симпсон (20 августа1710 - 14 мая1761) - английскийматематик. В 1746 году Симпсон избран в членыЛондонского королевского общества, а ранее - в члены основанного в 1717 году в Лондоне Математического общества. В 1758 избран иностранным членомШведской королевской академии наук. Назначенный профессором вКоролевскую военную академиювВулидже, Симпсон составил учебники поэлементарной математике. В особых отделахгеометриирассматриваются задачи о наибольших и наименьших величинах, решаемые с помощью элементарной геометрии,правильные многогранники, измерение поверхностей, объёмы тел и, наконец, смешанные задачи.

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

где b 1 - площадь (длина) нижнего основания

b 2 - площадь (длина) среднего основания

b 3 - площадь (длина) верхнего основания

2.1 Применение формулы Симпсона для вывода формул площадей плоских фигур.

Наша универсальная формула.b 1 = b 2 =b 3 , тогда получаем:

Ответ: S= hb 1

Вывод. Действительно, площадь параллелограмма равна произведению основания на высоту.

Универсальная формула.

Так какАВСД-трапеция, то b 2 -ее средняя линия, значит

Тогда получаем:

Вывод. Действительно, площадь трапеции равна половине произведения двух оснований на высоту.

Проведя аналогичные доказательства (Приложение 3-4) для формул площадей треугольника, прямоугольника, квадрата и ромба, я пришел к выводу, что универсальная формула Симпсона подошла для вычисления площадей таких плоских фигур как: параллелограмм, трапеция, треугольник, квадрат, ромб, прямоугольник.

2.2. Применение формулы Симпсона для вывода формул объемов пространственных тел.

Так какb 1 =b 2 =b 3 , тогда получаем:

Ответ: V=b 1 h

Доказательство, предложенное в учебнике геометрии авт. Л.С.Атанасяна в Приложении 6.

Вывод. Действительно, объем призмы равен произведению площади основания на высоту. Аналогично проводится доказательство выведения формулы объема цилиндра (Приложение 5)

Решение:Так как b 1 =0, а,то тогда получаем:

Доказательство, предложенное в учебнике геометрии авт. Л.С.Атанасяна в Приложении 9.

Вывод. Действительно, объем конуса равен одной трети произведения площади основания на высоту.Аналогично проводится доказательство выведения формулы объема пирамиды (Приложение 5)

Тогда получаем:

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

Задача 6. Объем шара.

Дано: шар

b 3 - площадь верхнего основании

Найти: Vшара.

(Рис. 11.Шар)

Так какb 1 =b 3 =0, h=2R

Тогда получаем:

Доказательство, предложенное в учебнике геометрии авт. Л.С.Атанасяна в Приложении 10

Вывод: Формулы объемов всех пространственных тел, изучаемые в 11-м классе, также легко выводятся с помощью универсальной формулы Симпсона.

2.3 Практическое применение формулы

Следующим этапом моего исследования является практическое применение (см.Приложение 11-12)

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

Я располагаю формулой, по которой можно приближенно вычислить объем ствола дерева, не задаваясь вопросом о том, на какое геометрическое тело оно похоже: на цилиндр, на полный конус или на усеченный конус. Зная плотности различных пород древесины, можно вычислить вес дерева на корню. Я решил эту задачу с помощью вычисления объема ствола, как объем цилиндра, диаметр основания которого равен диаметру ствола посредине длины: при этом результат получается, однако, преуменьшенный, иногда на 12 %. Без большой ошибки можно принимать объем дерева на корню половину объема цилиндра той же высоты с диаметром, равным поперечнику дерева на высоте груди.

Проделав расчеты, по известным нам ранее формулам, я вычислил объем ствола дерева на корню (см. Приложение 13)

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

Глава 3. Анкетирование учащихся

3.1 Исследование и опрос

Среди учащихся 11-х классов я провел исследование (см.Приложение 13).

Цель исследования: определение количества формул, которые учащиеся могут воспроизвести без повторения за 10 минут, т.е. объема «остаточных» формул.

Результаты оказались следующими (см.Приложение 14):

Наибольшее количество воспроизведенных формул - 41, наименьшее - 5. Учитывая то, что количество формул могло достигать 500 за неограниченное время, я пришел к выводу, что огромное количество формул, изучаемых в школе, учащиеся не помнят. Воспроизведенные формулы составляют лишь 8,2 % от общего количества изученных формул. Чаще всего учащиеся воспроизводили формулы по алгебре (формулы тригонометрии, логарифмические формулы, формулы сокращенного умножения, формула корней квадратного уравнения, производные); по геометрии (формулы площадей плоских фигур, некоторые объемы пространственных тел); несколько формул по физике (формула кинетической энергии, силы тяжести, силы трения и МКТ); по информатике () Это было естественно, т.к. в математике формул больше, чем в любой другой науке.

Увидев полученные результаты, я решил определить причины столь низкого результата. Мною был проведен опрос (см. приложение 14-15) учащихся 11-х классов, в котором предлагалось ответить на следующие вопросы:

Вопросы анкеты.

Как Вы считаете, сколько примерно формул должен знать выпускник школы?

А) зазубривание

Б) понимание

В) метод ассоциаций

Г) другое

Результаты оказались следующими (см.Приложение 15).

Вопрос 1. От 60 до 250 формул

Вопрос 2 . Из полученных ответов можно сделать вывод, что учащиеся 11-х классов при заучивании формул стараются их понять или применяют зазубривание.

Вопрос 3. Мнение учащихся по данному вопросу разошлись, хотя по диаграмме видно, что в основном отвечали «да», т.е. учащиеся считают, что количество формул для запоминания соответствуют уровню памяти среднего ученика.

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

3.2 Тестирование

Теперь я знаю, что формула Симпсона действительно универсальна, и её вполне можно применять в жизни. Но действительно ли она так необходима? Чтобы ответить на этот вопрос, я представил формулу на уроке 11 классу, после чего провел тестирование (см. приложение 16-17), и получил следующие результаты:

Тест № 1

23% признались, что им трудно запомнить все формулы.

17% сказали, что выучить все формулы им не составляет труда, в том числе и формулу Симпсона.

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

Тест № 2

100% утверждают, что формула Симпсона запоминается им легко.

0% признались, что испытывают некоторые трудности в её запоминании.

Тест № 3

76% будут применять эту формулу в дальнейшем.

24% признались, что она им вряд ли понадобится.

Тест № 4

82% считают, что формулу Симпсона стоит включить в школьную программу.

0% считают, что формулу не стоит включать в школьную программу.

18% утверждают, что формулу стоит включить в школьную программу, но только в профильных классах.

Тест № 5

35% считают, что помнить одну формулу для определения объёма сразу нескольких геометрических тел гораздо проще.

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

6% считают, что достаточно помнить только формулы, включённые в школьную программу.

Эту формулу так же можно применить в решении задач, в том числе и на ЕГЭ. Приведу примеры задач, которые были даны в 11 классе, и которые были решены учениками без труда:

Задача1 Правильная шестиугольная призма с высотой 18см вписана в цилиндр, с радиусом основания 4см. Найдите объём призмы.

Задача2 Правильная четырехугольная пирамида, с высотой 24см и стороной основания 5см, вписана в цилиндр. Найдите объём цилиндра.

Вывод:

Заключение

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

Я предположил, что формула английского математика Томаса Симпсона

позволит заменить формулы площадей фигур и объемов тел одной формулой.

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

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

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

Теперь прогуливаясь по лесу, вам наверно будет, вероятно, интересно определить объём любого дерева. Вычислить сколько в нём кубических метров древесины, а заодно и взвесить его - узнать, можно ли было бы, например, увезти такой ствол на одной телеге.

Считаю, свою работу полезной, т.к. мною были выведены все формулы площадей и объемов изучаемых в школе.

Из результатов анкетирования я убедился, что формула Симпсона достаточно проста для запоминания, и её стоит включить в школьную программу.

Эту формулу так же можно применять на экзаменах, включая ЕГЭ.

Список использованной литературы:

Я.И.Перельман. Занимательная алгебра. Занимательная геометрия. - М., «АСТ»,1999.

CD-ROM. Большая энциклопедия Кирилла и Мефодия, 2002.

Л.С. Атанасян и др. Геометрия 10-11 . Учебник для общеобразовательных учреждений,- М., «Просвещение», 2002.

https://ru.wikipedia.org/wiki

https://studfiles.net/preview/5433881/page:10/

https://studopedia.ru/6_126004_formula-simpsona.html

https://vuzlit.ru/940376/vyvod_formuly_simpsona

Приложение 1

Краткие характеристики свойств геометрических тел

Треугольник

Приложение 2

Прямоугольник

Приложение 3

b 3 =0, так как верхнее основание является точкой.

Так как b 2 - является в треугольнике средней линией, то, тогда получаем:

Вывод. Действительно, площадь треугольника равна половине произведения основания на высоту.

Решение: - универсальная формула.

Так как АВСД- квадрат, то b 1 =b 2 =b 3 =h, тогда получаем

Приложение 4

Вывод. Действительно, площадь квадрата равна квадрату его стороны.

Решение: - универсальная формула.

Так как АВСД - прямоугольник, то b 1 =b 2 =b 3 , тогда получаем:

Ответ: S=hb 1 .

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

Решение: - универсальная формула.

b 1 =b 2 =b 3 , тогда получаем:

Приложение 5

Задача 2. Объем цилиндра.

Дано: Цилиндр

b 1 - площадь нижнего основания:

b 2 -площадь среднего сечения:

b 3 - площадь верхнего основания.

Найти: Vцилиндра

(Рис. 22. Цилиндр)

Т.к. b 1 =b 2 =b 3 , тогда получаем:

Ответ: V=b 1 h

Доказательство, предложенное в учебнике геометрии авт. Л.С.Атанасяна в Приложении 7.

Вывод. Действительно, объем цилиндра равен произведению площади основания на высоту.

Решение:Так как b 3 =0, а, то тогда получаем:

Ответ: Доказательство, предложенное в учебнике геометрии авт. Л.С.Атанасяна в Приложении 8.

Приложение 6

Приложение 7.

Приложение 8

Приложение 9.

Приложение 10

Приложение 11

Задача № 1. Вычисляем объём модели куба по обычной формуле. Для этого измеряем ребро модели куба: а = 10,5 см. V=a 3 = 1157,625 cм 3

Задача № 2. Вычисляем объём модели правильной шестиугольной пирамиды по обычной формуле. Для этого измеряем высоту модели h = 17,2 см и сторону основания а = 6,5 см.

Задача № 3. Вычисляем объём модели цилиндра по обычной формуле. Для этого измеряем высоту модели h = 20,4 см и радиус основания R = 14 см.

Приложение 12

	Вычисляем S = π R 2 = 3,14 14 2 см 2 , V =Sh = 3,14196*20,4 = 12554,976 cм 3 Вычисляем объем модели по формуле Симпсона V = h/6(S нижнего основания + S верхнего основания + 4S среднего сечения):
Площади верхнего, нижнего основания и среднего сечения равны между собой S = π R 2 = 3,14 14 2 = 615,44см 2 , h= 20,4 см. V =20,4/6*(20,4+20,4)=12554,976 см 3
Задача № 4. Вычисляем объём модели конуса по обычной формуле. Для этого измеряем высоту модели h = 21 см и радиус основания R = 6 см.
Задача № 5. Вычисляем объём модели шара по обычной формуле. Для этого измеряем радиус шара R = 7 см.

Приложение 13

Расчёт для берёзы :

Расчёт для осины.

Расчёт для сосны.

Приложение 14

Результаты исследования «Определение объема «остаточных» формул»

Диаграмма 1. Определение количества «остаточных» формул.

Диаграмма 2. Предметы, по которым указаны формулы.

Приложение 15

Какой способ для запоминания формул Вы используете?

А) зазубривание

Б) понимание

В) метод ассоциаций

Г) другое

Диаграмма 3. Методы запоминания формул

Считаете ли Вы, что количество формул для заучивания соответствует уровню памяти среднего ученика?

Диаграмма 4. Соответствие количества формул уровню памяти среднего ученика

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

Диаграмма 5. Необходимость применения универсальной формулы

Приложение 16

Приложение 17

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

Напомним простейшие формулы численного интегрирования.

Вычислим приближенное численное значение . Интервал интегрирования [а, b] разобьем на п равных частей точками деления
, называемыми узлами квадратурной формулы. Пусть в узлах известны значения
:

Величина

называется интервалом интегрирования или шагом. Отметим, что в практике -вычислений число я выбирают небольшим, обычно оно не больше 10-20.На частичном интервале

подынтегральную функцию заменяют интерполяционным многочленом

который на рассматриваемом интервале приближенно представляет функцию f (х).

а) Удержим в интерполяционном многочлене только один первый член, тогда

Полученная квадратная формула

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

б) Удержим в интерполяционном многочлене два первых члена, тогда

(2)

Формула (2) называется формулой трапеций.

в) Интервал интегрирования
разобьем на четное число 2n равных частей, при этом шаг интегрирования h будет равен. На интервале
длиной 2h подынтегральную функцию заменим интерполяционным многочленом второй степени, т. е. удержим в многочлене три первых члена:

Полученная квадратурная формула называется формулой Симпсона

(3)

Формулы (1), (2) и (3) имеют простой геометрический смысл. В формуле прямоугольников подынтегральная функция f(х) на интервале
заменяется отрезком прямой у = ук, параллельной оси абсцисс, а в формуле трапеций - отрезком прямой
и вычисляется соответственно площадь прямоугольника и прямолинейной трапеции, которые затем суммируются. В формуле Симпсона функция f(х) на интервале
длиной 2h заменяется квадратным трехчленом - параболой
вычисляется площадь криволинейной параболической трапеции, затем площади суммируются.

ЗАКЛЮЧЕНИЕ

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

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

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

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

Интеграл от любой рациональной функции может быть выражен через элементарные функции в конечном виде, а именно:

через логарифмы- в случаях простейших дробей 1 типа;

через рациональные функции- в случае простейших дробей 2 типа

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

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

Формула Ньютона – Лейбница представляет собой общий подход к нахождению определенных интегралов.

Что касается приемов вычисления определенных интегралов, то они практически ничем не отличаются от всех тех приемов и методов.

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

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

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

Если необходимо получить наиболее точный результат, идеально подходит метод Симпсона .

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

Кафедра «Высшей математики»

Выполнил: Матвеев Ф.И.

Проверила: Бурлова Л.В.

Улан-Удэ.2002

1.Численные методы интегрирования

2.Вывод формулы Симпсона

3.Геометрическая иллюстрация

4.Выбор шага интегрирования

5.Примеры

1. Численные методы интегрирования

Задача численного интегрирования заключается в вычислении интеграла

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

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

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

Численные методы условно можно сгруппировать по способу аппроксимации подынтегральной функции.

Методы Ньютона-Котеса основаны на аппроксимации функции

полиномом степени . Алгоритм этого класса отличается только степенью полинома. Как правило, узлы аппроксимирующего полинома – равноотносящие.

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

сплайном-кусочным полиномом.

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

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

суммарная погрешность погрешность усечения

погрешность округления

Независимо от выбранного метода в процессе численного интегрирования необходимо вычислить приближенное значение интеграла и оценить погрешность. Погрешность уменьшается при увеличении n-количества

разбиений отрезка

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

за счет суммирования значений интегралов, вычисленных на частичных отрезках.

Погрешность усечения зависит от свойств подынтегральной функции и длины

частичного отрезка.

2. Вывод формулы Симпсона

Если для каждой пары отрезков

построить многочлен второй степени, затем проинтегрировать его и воспользоваться свойством аддитивности интеграла, то получим формулу Симпсона. Рассмотрим подынтегральную функцию на отрезке . Заменим эту подынтегральную функцию интерполяционным многочленом Лагранжа второй степени, совпадающим с в точках :

Проинтегрируем

и называется формулой Симпсона.

Полученное для интеграла

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

Оценим теперь погрешность интегрирования по формуле Симпсона. Будем считать, что у

на отрезке существуют непрерывные производные . Составим разность

К каждому из этих двух интегралов уже можно применить теорему о среднем, поскольку

непрерывна на и функция неотрицательна на первом интервале интегрирования и неположительна на втором (то есть не меняет знака на каждом из этих интервалов). Поэтому:

(мы воспользовались теоремой о среднем, поскольку

- непрерывная функция; ).

Дифференцируя

дважды и применяя затем теорему о среднем, получим для другое выражение: , где

Из обеих оценок для

следует, что формула Симпсона является точной для многочленов степени не выше третьей. Запишем формулу Симпсона, напрмер, в виде: , .

Если отрезок

интегрирования слишком велик, то его разбивают на равных частей (полагая ), после чего к каждой паре соседних отрезков , ,..., применяют формулу Симпсона, именно:

Запишем формулу Симпсона в общем виде.

В этом методе предлагается подынтегральную функцию на частичном отрезке аппроксимировать параболой, проходящей через точки
(x j , f (x j )), где j = i -1; i -0.5; i , то есть подынтегральную функцию аппроксимируем интерполяционным многочленом Лагранжа второй степени:

Проведя интегрирование, получим:

Это и есть формула Симпсона или формула парабол. На отрезке
[a, b ] формула Симпсона примет вид

Графическое представление метода Симпсона показано на рис. 2.4.

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

Избавимся в выражении (2.16) от дробных индексов, переобозначив переменные:

Тогда формула Симпсона примет вид

Погрешность формулы (2.18) оценивается следующим выражением:

где h·n = b - a , . Таким образом, погрешность формулы Симпсона пропорциональна O (h 4 ).

Замечание. Следует отметить, что в формуле Симпсона отрезок интегрирования обязательно разбивается на четное число интервалов.

10.5. Вычисление определенных интегралов методами
Монте–Карло

Рассматриваемые ранее методы называются детерминированными , то есть лишенными элемента случайности.

Методы Монте–Карло (ММК) – это численные методы решения математических задач с помощью моделирования случайных величин. ММК позволяют успешно решать математические задачи, обусловленные вероятностными процессами. Более того, при решении задач, не связанных с какими-либо вероятностями, можно искусственно придумать вероятностную модель (и даже не одну), позволяющую решать эти задачи. Рассмотрим вычисление определенного интеграла

При вычислении этого интеграла по формуле прямоугольников интервал [a, b ] разбиваем на N одинаковых интервалов, в серединах которых вычислялись значения подынтегральной функции. Вычисляя значения функции в случайных узлах, можно получить более точный результат:

Здесь γ i - случайное число, равномерно распределенное на интервале
. Погрешность вычисления интеграла ММК ~ , что значительно больше, чем у ранее изученных детерминированных методов.

На рис. 2.5 представлена графическая реализация метода Монте-Карло вычисления однократного интеграла со случайными узлами (2.21) и (2.22).

(2.23)

Рис. 10.6. Интегрирование методом Монте-Карло (2-й случай)

Как видно на рис. 2.6, интегральная кривая лежит в единичном квадрате, и если мы сумеем получать пары случайных чисел, равномерно распределенных на интервале , то полученные значения (γ 1, γ 2) можно интерпретировать как координаты точки в единичном квадрате. Тогда, если этих пар чисел получено достаточно много, можно приблизительно считать, что
. Здесь S – число пар точек, попавших под кривую, а N – общее число пар чисел.

Пример 2.1. Вычислить следующий интеграл:

Поставленная задача была решена различными методами. Полученные результаты сведены в табл. 2.1.

Таблица 2.1

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

11 ПРИБЛИЖЕННОЕ РЕШЕНИЕ НЕЛИНЕЙНЫХ
И ТРАНСЦЕНДЕНТНЫХ УРАВНЕНИЙ