Автор: Пользователь скрыл имя, 10 Января 2012 в 17:41, реферат
Обычный прием построения квадратурных формул состоит в замене подынтегральной функции f(x) на отрезке [a,b] интерполирующей или аппроксимирующей функцией g(x) сравнительно простого вида, например полиномом, с последующим аналитическим интегрированием. Это приводит к следующему представлению интеграла:
Министерство образования и науки Российской Федерации
ФГАОУ ВПО
«Уральский Федеральный Университет
имени первого
Президента России Б.Н. Ельцина»
Институт радиоэлектроники и информационных технологий - РТФ
Кафедра
«Автоматика и информационные технологии»
Реферат
Тема: Целостный педагогический
процесс: структура,
движущие силы, принципы и закономерности
по дисциплине
«Психология и педагогика»
Группа:
Студент:
Преподаватель: Баландина
Т.Ю.
Екатеринбург
2011
МЕТОДЫ
ВЫЧИСЛЕНИЯ ОПРЕДЕЛЕННОГО
ИНТЕГРАЛА
Задача численного интегрирования функций заключается в вычислении приближенного значения определенного интеграла
на основе ряда значений подынтегральной функции f(x):{ f(xi) = yi, i=0,n}.
Формулы численного вычисления однократного интеграла называются квадратурными формулами, двойного и более кратного – кубатурными.
Обычный прием построения квадратурных формул состоит в замене подынтегральной функции f(x) на отрезке [a,b] интерполирующей или аппроксимирующей функцией g(x) сравнительно простого вида, например полиномом, с последующим аналитическим интегрированием. Это приводит к следующему представлению интеграла:
Пренебрегая остаточным членом R[f], получаем приближенную формулу для вычисления интеграла:
Обозначим через yi = f(xi) значение подынтегральной функции в различных точках на [a,b]. Квадратурные формулы являются формулами замкнутого типа, если x0 = a , xn = b.
В качестве приближенной функции g(x) рассмотрим интерполяционный полином на в форме полинома Лагранжа:
где , при этом , здесь – остаточный член интерполяционной формулы Лагранжа.
Формула (1) дает
где
В
формуле (2) величины {
} называются узлами, {
} – весами,
− погрешностью квадратурной формулы.
Если веса {
} квадратурной формулы вычислены по
формуле (3), то соответствующую квадратурную
формулу называют квадратурной формулой
интерполяционного типа.
Подведем итог:
Рассмотрим частные случаи формул (2) и (3): метод прямоугольников, трапеций, парабол (метод Симпсона). Названия этих методов обусловлены геометрической интерпретацией соответствующих формул.
Определенный интеграл от функции f(x) численно равен площади криволинейной трапеции, ограниченной кривыми у = 0, x = a,
x = b, y = f(x) (рис. 1).
Рис.
1. Площадь под кривой y
= f(x)
Для вычисления этой площади весь отрезок интегрирования [a,b] разбивается на n равных подынтервалов длины h=(b-a)/n. Площадь под подынтегральной кривой приближенно заменяется на сумму площадей прямоугольников, как это показано на рис. 2.
Рис.
2. Площадь под кривой y
= f(x), аппроксимирующаяся суммой площадей
прямоугольников
Сумма площадей всех прямоугольников вычисляется по формуле
Метод, представленный формулой (4), называется методом левых прямоугольников, а метод, представленный формулой (5), – методом правых прямоугольников:
Погрешность
вычисления интеграла определяется
величиной шага интегрирования h.
Чем меньше шаг интегрирования, тем точнее
интегральная сумма S аппроксимирует
значение интеграла I. Исходя из этого
строится алгоритм для вычисления интеграла
с заданной точностью. Считается, что интегральная
сумма S представляет значение интеграла
I c точностью eps, если разница по
абсолютной величине между интегральными
суммами
и
, вычисленными с шагом h и h/2
соответственно, не превышает eps.
Метод
средних прямоугольников
Для
нахождения определенного интеграла
методом средних
Рис. 3. Площадь под кривой y=f(x), аппроксимирующаяся суммой площадей прямоугольников
где ;
n – количество разбиений отрезка [a,b].
Для нахождения определенного интеграла методом трапеций площадь криволинейной трапеции разбивается на n прямоугольных трапеций с высотами h и основаниями у1, у2, у3, ..уn, где n − номер прямоугольной трапеции. Интеграл будет численно равен сумме площадей прямоугольных трапеций (рис. 4).
Рис.
4. Площадь под кривой y
= f(x), аппроксимирующаяся суммой площадей
прямоугольных трапеций
Вычисление определенного
интеграла идет по формуле (6):
где ; n – количество разбиений отрезка [a,b].
Погрешность формулы трапеций оценивается числом
Погрешность формулы трапеций с ростом уменьшается быстрее, чем погрешность формулы прямоугольников. Следовательно, формула трапеций позволяет получить большую точность вычисления определенного интеграла, чем метод прямоугольников.
Если для каждой пары отрезков построить многочлен второй степени, затем проинтегрировать его на отрезке и воспользоваться свойством аддитивности интеграла, то получим формулу Симпсона.
В
методе Симпсона для вычисления определенного
интеграла весь интервал интегрирования
[a,b] разбивается на подынтервалы
равной длины h=(b−a)/n. Число отрезков
разбиения является четным числом.
Затем на каждой паре соседних подынтервалов подынтегральная функция f(x) заменяется многочленом Лагранжа второй степени (рис. 5).
Рис.
5. Функция y = f(x) на отрезке
, заменяемая многочленом 2-го порядка
Рассмотрим
подынтегральную функцию
на отрезке
. Заменим эту подынтегральную функцию
интерполяционным многочленом Лагранжа
второй степени, совпадающим с y=
в точках
:
Проинтегрируем на отрезке :
Введем замену переменных:
Учитывая формулы замены,
Выполнив интегрирование, получим формулу Симпсона:
Полученное для интеграла значение совпадает с площадью криволинейной трапеции, ограниченной осью , прямыми , и параболой, проходящей через точки На отрезке формула Симпсона будет иметь вид
В
формуле параболы значение функции
f(x) в нечетных точках разбиения х1,
х3, ..., х2n-1 имеет
коэффициент 4, в четных точках х2,
х4, ..., х2n-2
− коэффициент 2 и в двух граничных точках
х0 = а, хn
= b − коэффициент 1.
Геометрический смысл формулы Симпсона: площадь криволинейной трапеции под графиком функции f(x) на отрезке [a, b] приближенно заменяется суммой площадей фигур, лежащих под параболами.
Если
функция f(x) имеет на [a,
b] непрерывную производную четвертого
порядка, то абсолютная величина погрешности
формулы Симпсона не больше, чем
, где М − наибольшее значение
на отрезке [a, b]. Так как n4
растет быстрее, чем n2, то погрешность
формулы Симпсона с ростом n уменьшается
значительно быстрее, чем погрешность
формулы трапеций.
f(x)=x
f(x)=x2
f(x)= x3
f(x)= x4
на отрезке [0, 1] с шагом , , ;
Метод | шаг | 1 | x | x2 | x3 | x4 |
Метод левых прямоугольников | 0,05 | 1 | 0,475 | 0,3088 | 0,2256 | 0,1758 |
Метод средних прямоугольников | 0,05 | 1 | 0,5 | 0,3331 | 0,2497 | 0,1996 |
Метод трапеции | 0,05 | 1 | 0,5 | 0,3338 | 0,2506 | 0,2008 |
Формула Симпсона | 0,05 | 1 | 0,5 | 0,3333 | 0,25 | 0,2 |
Метод левых прямоугольников | 0,1 | 1 | 0,45 | 0,285 | 0,2025 | 0,1533 |
Метод средних прямоугольников | 0,1 | 1 | 0,5 | 0,3325 | 0,2488 | 0,1983 |
Метод трапеции | 0,1 | 1 | 0,5 | 0,335 | 0,2525 | 0,2033 |
Формула Симпсона | 0,1 | 1 | 0,5 | 0,3333 | 0,25 | 0,2 |
Метод левых прямоугольников | 0,2 | 1 | 0,4 | 0,24 | 0,16 | 0,1133 |
Метод средних прямоугольников | 0,2 | 1 | 0,5 | 0,33 | 0,245 | 0,1934 |
Метод трапеции | 0,2 | 1 | 0,5 | 0,34 | 0,26 | 0,2133 |
Формула Симпсона | 0,2 | 0,8 | 0,333 | 0,1947 | 0,1349 | 0,1051 |