Matematika_Praktikum

5.Каковы характерные особенности применения метода замены переменной к вычислению определенного интеграла?

6.Каковы характерные особенности применения метода интегрирования по частям к вычислению определенного интеграла?

7.В чем состоит геометрический смысл определенного интеграла?

8.Приведите примеры применения интегрального исчисления в биологии и медицине.

Глава 9. Численное интегрирование

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

9.1. Формула прямоугольников

Если известны значения функции f(x) в некоторых точках x0, x1, … , xm, то в качестве функции «близкой» к f(x) можно взять многочлен Р(х) степени не выше m, значения которого в выбранных точках равны значениям функции f(x) в этих точках.

b b

f (x)dx P(x)dx

a a

Разобьем отрезок интегрирования [a; b] на n равных частей.

Обозначим длину каждой части x b a . При этом y0 = f(x0), y1 =

n

f(x1), …., yn = f(xn).

Составим суммы: y0 x + y1 x + … + yn-1 x

81

a

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

формулой прямоугольников.

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

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

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

Рис. 9.1 Очевидно, что чем больше взять точек n разбиения

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

y1 x + y2 x + … + yn x

Это соответственно нижняя и верхняя интегральные суммы. Первая соответствует вписанной ломаной, вторая – описанной.

Тогда

	b	b a
	f (x)dx		(y0 y1 ... yn 1)
		n
	a

		b	b a
	или	f (x)dx		(y1 y2 ... yn)	(9.1)
			n

y0 y1

yn 1 yn

x;

y1 y2

x; ...

,

x

2

2

b

2

y0 y1

y1 y2

yn 1 yn

f (x)dx

x

x ...

x

2

2

2

a

После приведения подобных слагаемых получаем формулу трапеций:

82

b

b a y

y

f (x)dx

0

n

y

y

... y

(9.2)

n 1

n

2

1

2

a

9.3. Метод средних

Пусть для определенности a = x0, b = xn. Обозначим через x1 ,

									2
x3 ,	x5 , … середины интервалов (x0,					x1),	(x1, x2),	(x2,	x3), …
2	2
Полагаем f(x1		) = y1	, f(x3	) = y3 , f(x5 ) = y5			, … . Тогда
	2	2	2	2	2	2
		b
		f (x)dx bna[y1			y3	... y2n 1		]	(9.3)
		a		2	2		2

Выражение (9.3)

(риc.

a	x1/2	x1 x3/2
b	x

определяет площадь				ступенчатой фигуры
				8.3.1).
				В большинстве
				случаев при данном n
				формула (9.3) точнее,
				чем (9.1) и (9.2). С
				увеличением n
				точность формул (9.1),
				(9.2), (9.3)
				неограниченно
				возрастает.
x2 x5/2

Рис. 9.2

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

(формула парабол или квадратурная формула)

Разделим отрезок интегрирования [a, b] на четное число отрезков п=2m. Площадь криволинейной трапеции, ограниченной графиком функции f(x) заменим на площадь криволинейной трапеции, ограниченной параболой второй степени с осью

83

симметрии, параллельной оси Оу и проходящей через точки кривой, со значениями f(x0), f(x1), f(x2).

Для каждой пары отрезков построим такую параболу (рис. 8.4.1).

у

Уравнения

этих

парабол

имеют вид Ax2 + Bx + C, где

коэффициенты А, В, С могут

быть легко найдены по трем

точкам пересечения параболы

с исходной кривой.

y0 Ax02 Bx0 C

y Ax2 Bx C

1

1

1

0 х0 х1 х2 х3

х4

х

y2 Ax22 Bx2 C

Рис. 9.3

(9.4)

Обозначим 2h x2

x0 .

x2

x3

x2

x2

(Ax2

S

Bx C)dx A

B

Cx

3

2

x0

x0

Если принять х0 = –h, x1 = 0, x2 = h, то

S

h

(2Ah2 6C)

(9.5)

3

Тогда уравнения значений функции (9.4) имеют вид:

y0 Ah2 Bh C

y1 C

y2 Ah2 Bh C

C учетом этого: y

0

4y

y

2

2Ah2

6C .

1

h

Отсюда уравнение (8.5.1) примет вид:S

(y

4y y

2

)

Тогда

3

0

1

x2

h

f (x)dx

(y0 4y1 y2 )

x0

3

x4

h

f (x)dx

(y2 4y3 y4 )

x2

3

...............................................

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

84

b

b a

f (x)dx 6m y0 y2m 2(y2 y4 ... y2m 2) 4(y1 y3 ... y2m 1)

a

Чем больше взять число m, тем более точное значение интеграла будет получено.

9.5. Примеры

№1. Найти приближенное значение определенного интеграла

8

x3 16dx с помощью формулы Симпсона и формулы трапеций,

2

разбив отрезок интегрирования на 10 частей.

Решение.

Так как п=2т, то в нашем примере т=5. По формуле Симпсона получим:

	8		dx	b a	(y0 y10 2 y2 y	4 y6 y8 4(y1 y3 y5 y7
		x3 16
	2			6m

8 ( 2)

y9 )) (y( 2) y(8) 2(y(0) y(2) y(4) y(6)) 4(y( 1) 6 5

y(1) y(3) y(5) y(7)))

Все дальнейшие расчеты приведены в таблице:

		i	0		1			2	3		4	5	6	7	8	9	10
	xi		–2			–1		0	1		2	3	4	5	6	7	8
y(xi)			2,83		3,87			4	4,12	4,9		6,56	8,94	11,87	15,23	18,95	22,
																	98
			Окончательно получим,
8				dx		8 2
		x3 16					(2,83 22,98 2(4 4,9 8,94 15,23) 4(3,87 4,12 6,56
2						6 5

11,87 18,95)) 91,14

Точное значение этого интеграла – 91.173.

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

Абсолютная погрешность равна 91,14 91,173 0,033.

Относительная погрешность 0,033 100% 0,036%

91,173

85

Для сравнения применим к этому же интегралу формулу трапеций.

	8			b a y				y						8 2		2,83 22,98
		x3 16dx					0		n	y1 y2	... yn 1
								2						10		2
	2			n

3,87 4 4,12 4,9 6,56 8,94 11,87 15,23 18,95) 91,35

Абсолютная погрешность равна 91,35 91,173 0,177.

Относительная погрешность 0,177 100% 0,194%

91,173

Формула трапеций дала менее точный результат по сравнению с формулой Симпсона.

1	dt
№2. Вычислить определенный интеграл			с помощью
	1 t	2
0

формулы прямоугольников, если п=10.
	Решение.
	По												формуле							прямоугольников					получим:
1	dt					b a						(y0 y1			y2		... y9 )			1 0	(y(0) y(0,1) y(0,2) ... y(0,9)).
	2																			10
0	1 t						n
	Резльтаты вычислений поместим в таблицу:
	i				0					1					2			3		4		5	6	7		8	9
	ti				0					0,1					0,2			0,3		0,4		0,5	0,6	0,7		0,8	0,9
y(ti)					1					0,99					0,96			0,9		0,86		0,80	0,74	0,67		0,61	0,55
																		2
	Таким образом,
	1			dt								1	(1 0,99 0,96 0,92 0,86 0,80 0,74 0,67 0,61 0,55) 0,81.
				2
	0 1 t									10
Точное значение этого интеграла – 0,79.
				Найдем точное значение интеграла
	1				dt										1
										acrtg(t)							arctg(1) arctg(0) /4 0.785
							2
	01 t
															0

Абсолютная погрешность равна 0.785 0.79 0,005.

Относительная погрешность 0,005 100% 0,64%

0.785

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

86

меньше 1%).

№3. Найти площадь криволинейной трапеции, ограниченной

графиком функции y = 1 и прямыми: y = 0, x = a = 1, x = b =11

x

методами:

а) прямоугольников; б) трапеций; в) Симпсона;

г) аналитическим с использованием формулы НьютонаЛейбница.

Построить график заданной функции с разбиением отрезка

a;b на n = 10 подынтервалов и график функции s 0x				1			dt на
					t
отрезке x a;b .
Решение.	Составим	таблицу	разбиения		отрезка

интегрирования на n = 10 равных участков с длинами интервалов

x	=		b a					11 1	1	(табл. 11.1). Во			второй		строке		таблицы
					n		10
представлены увеличенные в 10 раз значения															yk	f (xk )		1
																				xk
(k =					).
		1;11
								Таблица 9.1 Данные для численных методов
xk				1			2		3	4	5	6	7	8		9 10		11
10yk				10			7,07		5,77 5,00		4,47	4,08	3,78	3,54		3,33	3,16	3,02
	а)				Используя					формулу		прямоугольников				с	высотами,

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

87

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

n

Sпр yk x (1 0,707 0,577 0,5 0,333 0,316) 1 5,020

k 1

Та же формула прямоугольников, но с подстановкой в нее высот, равных правым значениям функции на концах подынтервалов, дает значение интеграла, равного площади ограниченных пунктиром прямоугольников:

n 1

Sпр yk x (0,707 0,577 0,5 ... 0,316 0,302) 1 4,671.

k 2

б) По формуле трапеции получим

			y1		y n 1			n				1			0,302
S	тр							y k x										0,707			... 0,316		1 4,671
		2													2
								k 2
Следует обратить внимание на очевидное равенство
							S			1							S			.
										2
												S		пр					пр
в) По формуле Симпсона (n=5–количество спаренных
подынтервалов)
		b a										n					n
S сим					y1			y2 n 1		4 y2 k						2 y2 k 1
			6n
							4 0,707				k 1						k 2				0,316
	11 1 1			0,302							0,5				0,408				0,354				4,637
				0,577												0,333
	6 5							0,447		0,378
			2

г) Определим точное значение интеграла, являющегося

табличным:

S		1	dxx			2	x
		11

11

1

		2		11	1			4,633 .

Найдем относительные ошибки определения площадей различными использованными методами численного

88

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

			S			S		5 ,02	4 ,633
				пр
	пр											0 ,08 .
				S				4 ,633
Аналогично:
пр	0,07	7		10 2 ; тр				7	10 3		сим 9		4
										;			10

9.6. Варианты заданий

Для всех вариантов выполнить:

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

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

Используемые формулы: Sприближенное Sточное

100%

Sточное

в). Сделать выводы о точности полученных результатов.

	0
2.	(x 5)exdx, n 8;
	4
	4			dx
3.					, n 12	;
			2
	2 x			4x 5

1

4.arctg2xdx, n 10;

0

2,4

5.sin2 xcosxdx, n 8;

	0
	0,8	x
6.				dx, n 8;
		2
	0	cos	x

89

1,2

7.x2 ln(x 1,5)dx, n 10;

0,2

8cosx

8.2 2 sin xdx, n 12;

1

9.ex sin xdx, n 10;

	0
	4
				1 ln x						dx, n 10;
10.
	1					x
11.	1				x	2		dx, n 10;
	0		1 x
	1
12.	x(2 x2 )2 dx,										n 10;
	0

1x2

13.2 2 x dx, n 10;

2xdx

14., n 12;

4	2		x
		2
1
15. (x2		1)3 dx,			n 10.
0

9.7. Контрольные вопросы

1.Что такое первообразная? Обладают ли первообразные одной функции свойством единственности?

2.Дайте определение, в том числе в виде математического выражения, неопределенного интеграла.

3.Что такое подынтегральная функция? подынтегральное выражение?

4.В чем заключается метод замены переменной в определенном интеграле? Метод интегрирования по частям?

5.Что называется определенным интегралом?

6.Запишите формулу Ньютона-Лейбница.

90