Пусть требуется найти определенный интеграл где непрерывна на отрезке .
В большинстве случаев решать эту задачу приходится приближенно численными методами, которые основаны на замене подынтегральной функции интерполяционным многочленом. Для достижения необходимой точности вычислений интервал интегрирования делится на более мелкие части и на каждом из этих участков производится соответствующая замена. Тогда весь интеграл представляется в виде линейной комбинации нескольких значений
При фиксированном n
или
(1)
Формулы вида (1) называются квадратурными.
называется остаточным членом.
Значения называются узлами или сеткой, n– число разбиений интервала интегрирования.
Если делится на равные части, то называется шагом сетки.
Приближенное значение интеграла будет зависеть от h, поэтому будем обозначать его J(h).
Функция на участке заменяется многочленом нулевой степени, т. е. постоянной для всех величиной тогда
Геометрически это означает замену графика
на участке прямой .
Интеграл, как площадь криволинейной трапеции, будет приближенно равен площади прямоугольника.
Формула прямоугольников:
. (2)
Остаточный член .
На каждом из участков разбиения заменяется многочленом первой степени. Геометрическая кривая заменяется секущей и интеграл берется приближенно равным площади трапеции с основаниями и высотой
Формула трапеции:
(3)
Остаточный член .
Интервал интегрирования делится на четное число участков, на каждой паре участков с узлами функция заменяется многочленом второй степени, т. е. кривая заменяется квадратичной параболой.
Формула метода Симпсона:
(4)
Остаточный член .
Контрольная работа
Задание №1 Решить уравнение методом половинного деления, хорд с точностью .
Задание №2 Найти решение данной системы. Метод Гаусса.
№ варианта | Коэффициенты при неизвестных | Свободные члены | |||
0,11270 | -2,39990 | 8,95146 | 0,75000 | 8,60527 | |
9,58778 | -3,45350 | 0,24300 | 1,46840 | 16,40216 | |
0,86400 | 4,23700 | -2,50200 | -1,72927 | -15,88846 | |
-0,28427 | -4,58674 | -1,85970 | 0,14940 | 10,90588 | |
1,11270 | -3,02270 | -10,91328 | 1,06140 | 11,56420 | |
8,40446 | -3,45350 | 0,12430 | 0,84560 | 5,25400 | |
-0,33640 | 5,11230 | -1,83880 | 16,03250 | -11,79026 | |
-0,28427 | 5,85754 | -2,48250 | -0,16200 | -13,67224 | |
1,42410 | -2,71130 | 9,60540 | 0,43860 | 6,30236 | |
0,33853 | -5,34326 | -2,17110 | -0,16200 | 12,83405 | |
-0,02500 | 5,11230 | -2,46160 | -16,71758 | -11,58650 | |
8,40446 | -2,83070 | 0,43570 | 1,15700 | 15,77090 | |
0,28640 | 5,11230 | -2,15020 | 16,60758 | -12,52887 | |
0,80130 | -2,39990 | -8,29752 | 0,75000 | 7,078579 | |
8,52378 | -2,83070 | -0,18710 | 1,46840 | -2,20182 | |
0,33853 | 4,72046 | -1,85970 | -0,16200 | -11,78629 | |
0,11270 | -2,71130 | -9,60540 | 0,75000 | 8,93943 | |
-8,99612 | -3,45350 | 0,12430 | 1,15700 | 1,07023 | |
0,02500 | 5,11230 | -2,15020 | 16,03250 | -11,77124 | |
-0,28427 | 5,23474 | -2,17110 | -0,16200 | -12,58937 | |
0,80130 | -2,71130 | 9,60540 | 1,06140 | 6,16237 | |
8,52378 | -3,14210 | -0,18710 | 1,15700 | 16,18665 | |
0,02500 | 8,00900 | -1,83880 | -14,66234 | -10,15728 | |
0,02713 | -5,34326 | -2,17110 | -0,47340 | 14,18018 | |
0,86400 | 4,80090 | -2,46160 | 16,60758 | -12,88453 | |
1,42410 | -2,39990 | -8,95146 | 0,43860 | 6,53240 | |
-10,17944 | -3,45350 | 0,3570 | 1,46840 | -0,61624 | |
-0,28427 | 5,23474 | -1,85970 | -0,47340 | -12,05482 | |
0,80130 | -3,02270 | 9,60540 | 0,75000 | 5,53137 | |
-0,28427 | -5,85754 | -2,48250 | -0,16200 | 15,60785 | |
-0,33640 | 5,11230 | -2,15020 | -16,71758 | -13,11164 | |
8,52378 | -3,45350 | -0,18710 | 0,84560 | 15,88634 | |
-0,33640 | 5,42370 | -2,46160 | -10,08774 | -14,95126 | |
1,42410 | -3,02270 | 10,25934 | 0,43860 | 4,97590 | |
8,99612 | -3,45350 | 0,43570 | 8,45600 | 15,15486 | |
-0,28427 | -5,83234 | -2,48250 | 0,14940 | 13,79060 | |
8,01300 | -2,71130 | -8,95146 | 0,75000 | 9,11636 | |
0,28427 | 5,20954 | -2,17110 | 0,14940 | -13,29494 | |
0,02300 | 5,42370 | -2,15020 | 16,71758 | -10,78791 | |
-9,11544 | -3,45350 | -0,18710 | 1,15700 | 1,72450 |
Задание №3 Вычислить определённый интеграл с точностью методом Симпсона.
№ | интеграл | № | интеграл | ||
0,001 | 0,0001 | ||||
0,0001 | 0,01 | ||||
0,01 | 0,001 | ||||
0,001 | 0,01 | ||||
0,0001 | 0,0001 | ||||
0,01 | 0,01 | ||||
0,001 | 0,0001 |
Литература
1. Демидович Б.П. Численные методы анализа / Б.П.Демидович, И.А. марон, Э.З. Шувалова. – м.: Наука, 1976. – 368с.
2. Соболь Б.В. Практикум по вычислительной математике / Б.В. Соболь, Б.Ч. Месхи, И.М. Пешхоев. – Ростов-на-Дону: Феникс, 2008. – 344 с.