3.3 Метод Трапеций
Теоретические сведения
Пусть нам требуется вычислить определенный интеграл, где y = f(x) непрерывна на отрезке [a; b]. Разобьем отрезок [a; b] на n равных интервалов длины h точками. В этом случае шаг разбиения определяется так же как в методе парабол. Рассмотрим функцию на элементарных отрезках [x(i-1);x(i)]. Возможны четыре случая (на рисунке показаны простейшие из них, к которым все сводится при увеличении n):
На каждом отрезке [x(i-1);x(i)] заменим функцию y = f(x) отрезком прямой, проходящей через точки и . На рисунке показаны синими линиями:
В качестве приближенного значения интеграла возьмем выражение
Давайте выясним, что означает в геометрическом смысле записанное приближенное равенство. В первом случае площадь криволинейной трапеции приближенно равна площади трапеции с основаниями f(x(i-1)), f(x(i)) и высотой h, в последнем случае определенный интеграл приближенно равен площади трапеции с основаниями -f(x(i-1)), -f(x(i))и высотой h, взятой со знаком минус. Во втором и третьем случаях приближенное значение определенного интеграла равно разности площадей красной и синей областей, изображенных на рисунке ниже.
Рис. 11
Рис. 12
Блок схема
Текст задачи
Const aa=1; bb=2; ee=0.00001; nn=1;
var a,b,s,ss,e,deltax,x:real;
k,n:integer;
function f(x:real):real;
begin
f:=exp(-x)-sqr(x-1)+1;
end;
begin
a:=aa; b:=bb; n:=nn;
repeat ss:=s;
k:=k+1;
s:=0;
n:=n*2;
DELTAx:=(b-a)/n;
x:=a;
while x<b do begin
s:=s+deltax*(f(x)+f(x+deltax))/2; x:=x+deltax;
end;
e:=abs((ss-s)/s);
writeln(k=,k, ,s=,s:0:7, ,e=,e:0:7);
until e<ee;
end.
Результаты
k |
s |
e |
|
1 |
0.8623688 |
1.0000000 |
|
2 |
0.8900041 |
0.0310508 |
|
3 |
0.8969094 |
0.0076990 |
|
4 |
0.8986355 |
0.0019208 |
|
5 |
0.8990670 |
0.0004800 |
|
6 |
0.8991749 |
0.0001200 |
|
7 |
0.8992018 |
0.0000300 |
|
8 |
0.8992086 |
0.0000075 |
Вывод
В данном методе небольшое количество итераций, средняя точность относительно методов Симпсона и прямоугольников.
4. Вывод о трех методах
Каждый способ приближённого решения определённого интеграла имеет свои преимущества и недостатки, в зависимости от поставленной задачи следует использовать конкретные методы. Если необходимо быстро получить решение, но нет необходимости в большой точности ответа, следует воспользоваться одним из методов прямоугольника. Если же необходимо получить наиболее точный результат, идеально подходит метод Симпсона. Метод трапеций даёт ответ более точный, чем метод прямоугольников, но методу Симпсона он сильно уступает, этот метод можно назвать «золотой серединой» между двумя другими.
- 1. Операции с файлами
- 2. Решение иррациональных уравнений
- 2.1 Метод хорд
- 2.2 Метод половинного деления
- 3. Вычисление определенного интеграла
- 3.2 Метод прямоугольников
- 3.3 Метод Трапеций
- Решение систем линейных алгебраических уравнений
- 4.1 Метод Гаусса
- 4.2 Метод Ньютона
- 4.3. Метод Зейделя
- 5. Решение дифференциальных уравнений
- 5.1 Метод Эйлера
- 5.2 Метод Рунге-Кутта
- 6. Ряды Фурье
- Список литературы
- 1.1. Математическое моделирование и численные методы
- Модуль 3. Основные численные методы. Роль численных методов
- 9. Математические модели и численные методы
- Численные методы построения математических моделей
- 6.1 Математические модели и численные методы решения задач в различных предметных областях
- 19. Численные методы решения научно-технических задач
- Тема 1. Математическое моделирование и численные методы.