Вычисление определенных интегралов методом прямоугольников

Размещено на http://www.allbest.ru/

Содержание

Введение

Глава 1. Суть метода прямоугольников

1.1 Оценка абсолютной погрешности метода средних прямоугольников

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

Глава 2. Примеры применения метода прямоугольников при приближенном вычислении определенных интегралов

2.1 Разберем каждый случай

2.2 Замечание

2.3 Подведем итог

2.4 Программа вычисления по методу левых прямоугольников

Заключение и выводы

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



Введение

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

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

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



Глава 1. Суть метода прямоугольников

Пусть функция y = f(x) непрерывна на отрезке [a; b]. Нам требуется вычислить определенный интеграл

Обратимся к понятию определенного интеграла. Разобьем отрезок [a;b] на n частей точками .Внутри каждого отрезка выберем точку .

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

.

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

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

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

Если отрезок интегрирования [a;b] разбить на РАВНЫЕ части длины h точками

(то есть ) и в качестве точек выбрать СЕРЕДИНЫ элементарных отрезков (то есть ), то приближенное равенство

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

называют шагом разбиения отрезка [a;b].

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

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

С геометрической точки зрения для неотрицательной функции y=f(x) на отрезке [a;b] точное значение определенного интеграла представляет собой площадь криволинейной трапеции, а приближенное значение по методу прямоугольников - площадь ступенчатой фигуры.

1.1 Оценка абсолютной погрешности метода средних прямоугольников

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

На каждом отрезке имеем приближенное равенство.

Абсолютную погрешность метода прямоугольников на i-ом отрезке вычисляем как разность между точным и приближенным значением определенного интеграла:

.

Так как



есть некоторое число и

, то выражение

в силу четвертого свойства определенного интеграла можно записать как

.

Тогда абсолютная погрешность формулы прямоугольников на i-ом элементарном отрезке будет иметь следующий вид

Если считать, что функция y = f(x) имеет в точкеи некоторой ее окрестности производные до второго порядка включительно, то функцию y = f(x) можно разложить в ряд Тейлора по степеням с остаточным членом в форме Лагранжа:

По свойствам определенного интеграла равенства можно интегрировать почленно:

Абсолютная погрешность формулы прямоугольников на отрезке

[a; b] равна сумме погрешностей на каждом элементарном интервале, поэтому

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

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

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

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

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

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

Абсолютная погрешность методов левых и правых прямоугольников оценивается как



Глава 2. Примеры применения метода прямоугольников при приближенном вычислении определенных интегралов

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

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

Формулировки задач примерно следующие:

· вычислить приближенно определенный интеграл методом прямоугольников, разбив отрезок интегрирования на n частей;

· Методом прямоугольников найти приближенное значение определенного интеграла с точностью до одной сотой (одной тысячной и т.п.).

2.1 Разберем каждый случай

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

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

Решение.

В нашем примере a = 4, b = 9, n = 10,.

Внимательно посмотрим на формулу прямоугольников

.

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

Вычислим шаг:

.

, то

.

Для i = 1 имеем

.

Находим соответствующее значение функции

.

Для i = 2 имеем

.

Находим соответствующее значение функции

.

И так продолжаем вычисления до i = 10. Для удобства представим результаты в виде таблицы.

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

Значение исходного определенного интеграла можно вычислить по формуле Ньютона-Лейбница:

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

Как видите, точное значение определенного интеграла отличается от значения, полученного по методу прямоугольников для n = 10, менее чем на шесть сотых долей единицы.

Графическая иллюстрация.

По условию имеем

a = 1, b = 2,.

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

Нам известно, что

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

,

то будет достигнута требуемая степень точности.

Найдем - наибольшее значение модуля первой производной подынтегральной функции на отрезке [1; 2]. В нашем примере это сделать достаточно просто.

Графиком функции производной подынтегральной функции является парабола, ветви которой направлены вниз, на отрезке [1; 2] ее график монотонно убывает. Поэтому достаточно вычислить модули значения производной на концах отрезка и выбрать наибольшее:

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

Таким образом:

Число n не может быть дробным (так как n - натуральное число - количество отрезков разбиения интервала интегрирования). Поэтому, для достижения точности 0.01 по методу правых или левых прямоугольников, мы можем брать любое n = 9, 10, 11, …Для удобства расчетов возьмем n = 10.

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

,

а правых прямоугольников

.

Для их применения нам требуется найти h и для n = 10.

Итак,

Точки разбиения отрезка [a; b] определяются как .

Для i = 0 имеем

.

Для i = 1 имеем и

.

И так далее до i = 10.

Полученные результаты удобно представлять в виде таблицы:

Подставляем в формулу левых прямоугольников:

Подставляем в формулу правых прямоугольников:

Вычислим точное значение определенного интеграла по формуле Ньютона-Лейбница:

Очевидно, точность в одну сотую соблюдена.

Графическая иллюстрация.



2.2 Замечание

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

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

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

Берем произвольное n (например, n = 5) и вычисляем приближенное значение интеграла. Далее удваиваем количество отрезков разбиения интервала интегрирования, то есть, берем n = 10, и вновь вычисляем приближенное значение определенного интеграла. Находим разность полученных приближенных значений для n = 5 и n = 10. Если абсолютная величина этой разности не превышает требуемой точности, то в качестве приближенного значения определенного интеграла берем значение при n = 10, предварительно округлив его до порядка точности. Если же абсолютная величина разности превышает требуемую точность, то вновь удваиваем n и сравниваем приближенные значения интегралов для n = 10 и n = 20. И так продолжаем до достижения требуемой точности.

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

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

Не будем подробно останавливаться на вычислениях.

Для n = 5 имеем

,

для n = 10 имеем

.

Так как , тогда берем n = 20. В этом случае

.

Так как , тогда берем n = 40. В этом случае

.

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

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

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

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

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

2.3 Подведем итог

При вычислении определенного интеграла методом прямоугольников (методом средних прямоугольников) пользуемся формулой

и оцениваем абсолютную погрешность как

.

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

и

соответственно. Абсолютную погрешность оцениваем как

.

2.4 Программа вычисления по методу левых прямоугольников

Program levii;{Метод левых прямоугольников}

uses crt;

var i,n:integer; a,b,h,x,xb,s:real;

function f(x:real):real;

begin f:=(1/x)*sin(3.14*x/2); end;

begin

clrscr;

write('Введите нижний предел интегрирования '); readln(a);

write('Введите верхний предел интегрирования '); readln(b);

write('Введите количество отрезков '); readln(n);

h:=(b-a)/n; s:=0; xb:=a;

for i:=0 to n-1 do

begin x:=xb+i*h; s:=s+f(x)*h; end;

writeln('Интеграл равен ',s:12:10); readln;

end.

a=1 b=2 n=10 S= 0,5318196161

a=1 b=2 n=20 S= 0,5067061037

a=1 b=2 n=100 S= 0,4866793606



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



Program pravii; {Метод правых прямоугольников}

uses crt;

var i,n:integer; a,b,h,x,xb,s:real;

function f(x:real):real;

begin f:=(1/x)*sin(3.14*x/2); end;

begin

clrscr;

write('Введите нижний предел интегрирования '); readln(a);

write('Введите верхний предел интегрирования '); readln(b);

write('Введите количество отрезков '); readln(n);

h:=(b-a)/n; s:=0; xb:=a;

for i:=1 to n do

begin x:=xb+i*h; s:=s+f(x)*h; end;

writeln('Интеграл равен ',s:12:10); readln;

end.

a=1 b=2 n=10 S= 0,4318992805

a=1 b=2 n=20 S= 0,4567459359

a=1 b=2 n=100 S= 0,4766873270





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

Program srednii; {Метод средних прямоугольников}

uses crt;

var i, n: integer; a, b, dx, x, s, xb : real;

function f(x : real):real;

begin f:=(1/x)*sin(3.14*x/2); end;

begin

clrscr;

write('Введите нижний предел интегрирования '); readln(a);

write('Введите верхний предел интегрирования '); readln(b);

write('Введите количество отрезков '); readln(n);

dx:=(b-a)/n; xb:=a+dx/2;

for i:=0 to n-1 do

begin x:=xb+i*dx; s:=s+f(x)*dx; end;

write('Интеграл равен ',s:15:10); readln;

end.

a=1 b=2 n=10 S= 0,4815925913

a=1 b=2 n=20 S= 0,4816593368

a=1 b=2 n=100 S= 0,4816806769





Заключение и выводы

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

Чем больше значение n, тем точнее значение интеграла.

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

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



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

1. Волков Е.А. Численные методы: Учебное пособие. М.: Наука, 1982

2. Лапчик М.П. Численные методы: Учебное пособие для вузов. - М.: Acadimia, 2004

3. Турчак Л.И. Основы численных методов: Учебное пособие - М.: Наука, 1987

4. Николаев Н.Н., Костерин В.А. Методы вычислений в задачах

электроаппаратостроения: Лабораторный практикум. Чебоксары: Изд-во Чуваш ун-та, 1992

5. Демидович Б.П., Марон И.А. Основы вычислительной математики. М.: Наука, 1966.

6. Самарский А.А., Гулин А.В. Численные методы. М.: Наука, 1989.

7. Ракитин В.И., Первуш В.Е. Практическое руководство по методам вычислений сприложением программ для ПК - М. Высшая школа, 1998

Размещено на Allbest.ru

...