ВЫЧИСЛЕНИЕ ИНТЕГРАЛОВ И РЕШЕНИЕ УРАВНЕНИЙ
Методы вычисления определенных интегралов
Приближённое вычисление определенного интеграла основано на геометрическом смысле интеграла и сводится к приближенному вычислению площади, ограниченной графиком подынтегральной функции f(x), прямыми x = a = x0, x = b = xn и осью OX (рис. 19.1).
Интервал [a, b] делится на n равных частей длиной . Рис. 19.1. График подынтегральной функции
Тогда значениям xi = xi-1 + h, i = 1,2, ..., n соответствуют значения yi = f(xi).
Метод прямоугольников. Согласно методу левых прямоугольников, искомая площадь вычисляется как сумма площадей прямоугольников, основание которых равно h, а высота равна соответственно y0 для первого прямоугольника, y1 – для второго и т.д. вплоть до последнего с высотой yn-1. Тогда
Для метода правых прямоугольников аналогично
Метод трапеций. По методу трапеций определяется сумма площадей трапеций, основаниями которых являются ординаты y0, y1 и т. д., а высоты равны h:
Погрешность метода оценивается как , где М – максимальное значение второй производной f ( x ) на отрезке [ a , b ].Используя это соотношение, можно определить количество точек, на которое делится отрезок, исходя из заданной погрешности.
Рассмотрим алгоритм метода трапеций(рис. 19.2):
|
|
1. Ввод a, b, n.
2. Вычисление , x = a + h,
s = 0.
3. Расчет s = s + f(x), x = x + h.
4. Если x > (b – h), то переход к пункту 5, иначе – переход к п. 3.
5. Вычисление значения интеграла
6. Вывод z.
Рис. 8.2. Схема алгоритма метода трапеций
2. Содержание задания
1. С помощью микрокалькулятора вычислить приближенно значение определенного интеграла из табл.1 для n =4. Номер варианта определяет преподаватель. Для нечетных по номеру вариантов использовать метод трапеций, для четных – метод парабол.
2. Написать программу, реализующую вычисление определённого интеграла соответствующим методом. Для всех вариантов принять n =20.
3. Выполнить вычисления в пакете MathCAD. Результаты сравнить между собой.
Таблица 1
ПРИБЛИЖЕННОЕ РЕШЕНИЕ УРАВНЕНИЙ апеций |
|
|
i = 1, 2, 3 …, 2n,
Алгоритм метода парабол:
1. Ввод a, b, n.
2. Вычисление , x = a + 2h, s1=0, s2=0, i=1.
3. Расчет s2 = s2 + f(x) , x = x + h, s1 = s1 + f(x), x = x + h, i = i + 1.
4. Если i < n – 1, то переход к п. 3, иначе – переход к следующему пункту.
5. Вычисление значения интеграла:
6. Вывод z.
Здесь s1 = y3 + y5 + … + y2n-1, а s2 = y2 + y4 + … + y2n-2.
Дата добавления: 2018-11-24; просмотров: 181; Мы поможем в написании вашей работы! |
Мы поможем в написании ваших работ!