Вычисление интеграла методом Симпсона
Министерство образования и науки
Российской Федерации
Воронежский государственный
технический университет
Кафедра полупроводниковой электроники
и нано электроники
Контрольная работа
по дисциплине информатика
На тему: "Вычисление интеграла
методом Симпсона"
Содержание
1. Теоретическая часть
1.1 Метод деления отрезка пополам
1.2 Метод трапеций
2. Алгоритм работы программы
3. Листинг программы
4. Пример работы программы
5. Вывод
Список литературы
1. Теоретическая часть
1.1
Метод деления отрезка пополам
Задача заключается в нахождении корней нелинейного уравнения
. (1)
Для начала итераций необходимо знать интервал 
значений 
, на концах которого функция принимает значения разных знаков:
Из непрерывности функции f и условия (2) следует, что на интервале
существует хотя бы один корень уравнения
(в случае наличия нескольких корней метод приводит к нахождению одного из них)
Выберем точку внутри интервала:
. (3)
Если 
, то корень найден. Если 
разобьём интервал на два: 
и 
. Теперь найдём новый интервал, в котором
функция меняет знак. Пусть 
и соответственно корень находится внутри интервала . Тогда обозначим 
и повторим описанную процедуру до достижения
требуемой точности. За количество итераций N первоначальный отрезок
делится в 2N раз.
1.2
Метод трапеций
Метод численного интегрирования функции одной переменной,
заключающийся в замене на каждом элементарном отрезке подынтегральной функции
на многочлен первой степени, то есть линейную функцию. Площадь под графиком
функции аппроксимируется прямоугольными трапециями. Алгебраический порядок
точности равен 1.
Если отрезок 
является элементарным и не подвергается дальнейшему
разбиению, значение интеграла можно найти по формуле:
. (4)
Это простое применение формулы для площади трапеции - полусумма
оснований, которыми в данном случае являются значения функции в крайних точках
отрезка, на высоту (длину отрезка интегрирования).
Метод трапеций быстро сходится к точному значению интеграла для
периодических функций, поскольку погрешность за период аннулируется. Метод
может быть получен путём вычисления среднего арифметического между результатами
применения формул правых и левых прямоугольников.
2.
Алгоритм работы программы
3. Листинг программы
program task;= 0.0000001;f (x: Real): Real;: = 4
- sqr (sqr (x));;g (x: Real): Real;x < 0 then: = - exp (ln (-x) / 3)if x
> 0 then: = exp (ln (x) / 3): = 0;;func (x: Real): Real;: = f (x) - g
(x);;bissection (a,b: Real): Real;: Real;: Boolean;: = 0;( (func (a) * func
(b)) > 0) then
WriteLn ('Здесь не может быть корня! ');(abs (func
(a)) < E) then: = aif (abs (func (b)) < E) then: = b: = (a + b) / 2;: =
true;t do(abs (func (x)) < E) or (abs (a - b) < E) then: = x;: =
false;(func (x) * func (a) < 0) then: = x: = x;: = (a + b) / 2;;;;;;Integral
(a,b: Real): Real;= 128;, x, res: Real;: = (b - a) / n;: = a;: = 0;: = res +
0.5 * (func (x) + func (x + h)) * h;: = x + h;x > b-h;: = res;;,b: Real;: =
bissection (-3, 0);: = bissection (0,3);('Корень 1 = ', a: 7: 7);('Корень 2 =
', b: 7: 7);('Площадь интеграла = ', Integral (a, b): 7: 7);;
end.
4.
Пример работы программы
Точное значение интеграла, посчитано средствами
математического пакета Mathcad:
.
5.
Вывод
В данной работе решена задача приближённого интегрирования
функции методом трапеций и нахождения пределов интегрирования методом деления
отрезка пополам.
В процессе создания работы разработан алгоритм решения
поставленной задачи. По этому алгоритму на языке Pascal составлена и отлажена
программа.
В ходе тестирования были получены результаты работы методом
деления отрезка пополам и метода трапеций.
Программа является полностью работоспособной, что
подтверждается результатами её тестированием.
Список
литературы
1. Н.С.
Пискунов. Дифференциальное и интегральное исчисление 370 - 373 с.
2. С.А.
Немнюгин. Turbo Pascal. Учебник для вузов.