Файл: Методические указания по выполнению курсовой работы по дисциплине (модулю) Введение в математический анализ.docx

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 30.11.2023

Просмотров: 135

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.


24
R0 Ri

i1

hf xi

i1

24 f

xdx. (3.12)


К интегралу в формуле (3.12) мы перешли, «используя» метод средних

прямоугольников для функции fx.

Формула (3.12) представляет собой теоретическую оценку погрешности вычисления интеграла методом средних прямоугольников, данная оценка является априорной, так как не требует знания значения вычисляемого интеграла. Степень, в которую возводится шаг h , называется порядком метода интегрирования. Метод средних прямоугольников имеет второй порядок точности. Аналогично можно получить априорные оценки погрешностей других рассмотренных ранее методов.

Оценим погрешность метода левых прямоугольников. Погрешность

интегрирования Riна отрезке xi1, xi равняется разности между точным

значением интеграла и его приближенным значением hfxi1 :

Ri

xi

f

xi1

xdx hf

h2



xi1 2

fxi1 .


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

n hn

hb




2
R0 Ri

i1

hf

i1

xi1

f xdx.

a


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

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

Так как значение интеграла на отрезке xi1, xi вычисляется по формуле

h fx


fx

, то погрешность может быть рассчитана по формуле:

2 i i1
Ri
xi

f

xi1
xdx

h

2
fi1

h3


fi
12
xi.


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

a,b может быть рассчитана следующим образом:


n  h2 n

  h2 b



a





12
R0 Ri

i1

hf xi

i1

12 f

xdx



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

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


2 f



b

0
R h xdx.

a

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


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

h2 b


24 f


R0

a


a
  h2 b

xdx.


R0


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

12 f

xdx.

h4 b

IV


a
R0 180 f

xdx.

Методы левых и правых прямоугольников являются методами первого порядка точности. Методы средних прямоугольников и трапеций имеют второй порядок точности, при этом метод трапеций обладает вдвое большей по абсолютной величине погрешностью. Поэтому, если подынтегральная функция задана аналитически, то предпочтительнее из методов второго порядка применять метод средних прямоугольников вследствие его меньшей погрешности. Метод Симпсона имеет четвертый порядок точности с очень малым числовым коэффициентом. Формула Симпсона позволяет получить очень высокую точность, если четвертая производная подынтегральной функции не слишком велика. В противном случае, методы второго порядка точности могут дать большую точность, чем метод Симпсона.
    1. 1   ...   7   8   9   10   11   12   13   14   15

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


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

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

Ih

вычисленное с шагом hприближенное значение интеграла I. Если затем

вычислить приближенное значение

Ih2

с шагом

, то в качестве оценки


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

величину

Ih 2 I

Ih 2 Ih .


При необходимости вычислить результат с заданной точностью

расчеты повторяют с последовательно уменьшающимся (обычно вдвое)

шагом hдо тех пор, пока не выполнится условие Можно применить указанное правило для

Ih 2 Ih

контроля

  . погрешности на


каждом частичном отрезке xi1, xi,


i 1, n. При этом длина очередного шага

hi xi xi1 , посредством последовательного уменьшения начальной длины

вдвое, устанавливается такой, чтобы выполнялось неравенство

Ih 2 Ih hi
. (3.13)



ii b a

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

n

hi

i1 ,

b a

то есть не будет превосходить заданного уровня погрешности.

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

без потери точности вычисления.

Одним из подходов к экономии вычислительных ресурсов ЭВМ при