37809

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

Лабораторная работа

Информатика, кибернетика и программирование

4 Формула Симпсона Формула Симпсона записывается так: . Погрешность формулы Симпсона прямо пропорциональна в четвертой степени. На практике как и в случае метода трапеций расчеты ведут на сгущающихся сетках и оценку погрешности формулы Симпсона осуществляют по формуле 5. Критерием завершения процесса вычисления определенного интеграла с заданной точностью методом Симпсона на сгущающихся сетках служит условие .

Русский

2013-09-25

248 KB

10 чел.

Лабораторная работа № 5

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

ЦЕЛЬ РАБОТЫ: изучить и программно реализовать на языке высокого уровня широко применяемые на практике численные методы вычисления одномерных и кратных определенных интегралов, исследовать их на тестовых задачах.

Элементы теории

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

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

Квадратурные формулы в общем виде записываются так:

 

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

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

Формула трапеций имеет следующий вид:

. (5.1)

Она базируется на двух положениях: интервал интегрирования

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

Погрешность формулы трапеций

, (5.2)

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

Вычисление определенного интеграла по формуле (5.1) в условиях ошибок округления сопровождается также  вычислительной ошибкой:

,

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

Использование для оценки погрешности формулы (5.2) вызывает определенные трудности вследствие необходимости вычисления . Поэтому на практике привлекают прием вычисления интеграла на сгущающихся сетках с шагом  и , где  и формулу Рунге для оценки главной составляющей погрешности:

, (5.3)

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

. (5.4)

Формула Симпсона

Формула Симпсона записывается так:

. (5.5)

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

Погрешность формулы Симпсона

прямо пропорциональна  в четвертой степени.

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

. (5.6)

Кубатурная формула Симпсона

Вычисление двойных интегралов

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

где  обозначает  .  

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

Тогда

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

Задание

  1.  Разработать, программно реализовать и исследовать на задачах, предложенных преподавателем (см. табл. 5.1), алгоритмы интегрирования функций одной переменной методами трапеций, Симпсона на сгущающихся сетках с критерием завершения вычислительного процесса в виде (5.4) и (5.6) для ; двух переменных – методом Симпсона на заданной преподавателем сетке.

Содержание электронного  отчета

1 Алгоритмы.

  1.  Тексты программ.
  2.  Задачи, результаты их решения.

Задачи

Таблица 5.1

Подинтегральная функция

либо

Интервал

Интервал

1

[0.8; 1.762]

2

[1.3; 2.621]

3

[0.6; 1.724]

4

[3.0; 4.254]

Продолжение табл. 5.1

Подинтегральная функция

либо

Интервал

Интервал

5

[0; 1.234]

6

[0; 1.047]

7

[1.2; 2.471]

8

[1.0; 2.835]

9

[1.0; 2.631]

10

[2.0; 3.104]

11

[0; 1.075]

12

[0; 4.0]

13

[0; π/2]

14

[0; π/4]

15

[0; 1.0]

16

[3.0; 29.0]

17

[0; ln5]

18

[1.0; 4.0]

19

[0; π]

20

[0; π/2]

21

[-1.0; 1.0)

22

[-1.0; 1.0)

23

[0; 1.0]

24

[0; 1.0]

25

[0; 1.0]

26

[0; 1.0]

Окончание табл. 5.1

Подинтегральная функция

либо

Интервал

Интервал

27

[0; 1.0]

28

[0; 1.0]

29

[0; 4.0]

[1.0; 2.0]

30

[3.0; 4.0]

[1.0; 2.0]

31

[0; 2.0]

[0; 1.0]

32

[-1.0; 1.0]

[-1.0; 1.0]

33

[0; π/2]

[0; π/4]

34

[0; 1.0]

[1.0; 2.0]

35

[0; 2.0]

[0.5; 1.5]

ЛИТЕРАТУРА

1. Мулярчик С. Г.  Численные методы. Мн., 2001.

2. Вержбицкий В. М. Численные методы. М., 2000.

37


 

А также другие работы, которые могут Вас заинтересовать

46356. ОСНВАНИЯ И ФУНДАМЕНТЫ 885.28 KB
  При выборе варианта учитывается: Стоимость Долговечность Срок возведения конструкции фундамента Возможность производства работ в зимнее время Сохранение естественной структуры грунта КОНСТРУКЦИИ ФУНДАМЕНТОВ Основные типы фундаментов на естественном основании Рис1.верхняя плоскость фундамента. нижняя плоскость фундамента. hф высота фундамента.
46357. КОНСПЕКТ ЛЕКЦИЙ ПО КУРСУ «КОТЕЛЬНЫЕ УСТАНОВКИ ПРОМЫШЛЕННЫХ ПРЕДПРИЯТИЙ» 1.85 MB
  Соответственно большой роли энергетики в промышленном производстве современные промышленные предприятия имеют сложные и многообразные энергетические системы состоящие из комплексов установок и устройств предназначенных для сжигания топлива и производства транспорта распределения и потребления электроэнергии теплоты сжатого воздуха газа кислорода. Широкое применение пара для производства электроэнергии в технологических процессах и в быту определяет использование в котлах более 25 вcегo добываемого топлива.2 ...
46359. Опытно-экспериментальная работа по развитию интереса дошкольников к математике на основе использования игровых занимательных задач 276 KB
  Задачи математического развития детей дошкольного возраста Ерофеева и другие которые разработали методику формирования математических представлений детей дошкольников. В математике заложены огромные возможности для развития мышления детей в процессе обучения с самого раннего возраста.
46360. ИНФОРМАЦИОННЫЕ СИСТЕМЫ С БАЗАМИ ЗНАНИЙ 553 KB
  Менеджмент организаций специализации Информационные системы в менеджменте Моделирование знаний в информационных системах Интеллектуальные расчетнологические и экспертные системы Математические концептуальные методологии проектирования систем Базовые модели математической теории...
46361. Транспортная энергетика. Методические указания 481.5 KB
  Выполнение работы прививает определенные навыки практического применения основных методик термодинамического анализа для расчетов термодинамических процессов расширяет знания об идеализированных циклах реальных тепловых машин позволяет на конкретных примерах усвоить методику энергетического анализа и расчета важнейших процессов и циклов. Паровые циклы циклы паросиловых установок ПТУ их рабочим телом является вещество которое в течение цикла находится то в парообразном то в жидком состоянии. Характерные особенности идеального цикла...
46362. АНАЛИТИЧЕСКИЕ И ИМИТАЦИОННЫЕ МОДЕЛИ 7.07 MB
  Получаемая информация о состояниях процесса используется для определения тех характеристик процесса которые нужно получить в результате моделирования. Содержательное описание в словесном выражении концентрирует сведения о физической природе и характеристиках элементарных явлений исследуемого объекта или процесса о степени и характере взаимодействия между ними месте и значении каждого элементарного явления в общем процессе функционирования рассматриваемой реальной системы. Тщательное изучение объекта сводится к наблюдению и фиксации...
46363. Проект и изготовление модели-макета фонтанной арматуры крестового типа 1.5 MB
  Фонтанирование нефтяной скважины это процесс движения нефти от её забоя к устью происходящий под действием пластовой энергии. Естественное фонтанирование нефтяной скважины возможно лишь при условии если ее забойное давление больше гидростатического давления создаваемого на забой весом столба газонефтяной смеси поднимающейся к устью. Фонтанная елка предназначена для контроля и регулирования фонтанной струи направления ее по тому или иному выкиду и в случае необходимости для закрытия скважины. С его помощью поддерживают заданный режим...
46364. Активность сетей Петри. Задача о чтении/записи 1.96 MB
  Причиной рассмотрения сохранения в сети Петри было распределение ресурсов в операционной системе ЭВМ. Сеть Петри на рис. Тупик в сети Петри это переход или множество переходов которые не могут быть запущены.