Файл: Методические указания по выполнению курсовой работы по дисциплине (модулю) Введение в математический анализ.docx
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 10.11.2023
Просмотров: 141
Скачиваний: 1
ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.
частичных отрезках не линейными функциями, а полиномами более высоких порядков, то ошибка соответствующего метода численного интегрирования должна уменьшиться. В методе Симпсона в качестве функции, с помощью
которой на частичных отрезках интегрирования заменяется исходная
функция
fx, выбрана парабола.
Пояснение. Более высокая точность вычисления интегралов обеспечивается при использовании параболической аппроксимации (полиномом второй
степени)
y ax2 bx c
по трем соседним точкам. Для нахождения
коэффициентов a, bи c полинома, проходящего через точки x0 , y0 ,
x1 , y1 , x2 , y2 , требуется найти решение следующей системы линейных уравнений
y0 ax2 bx
-
c,
0 0
1
y1 ax2 bx1 c,
y2 ax2 bx c,
2 2
относительно неизвестных a, b, c.
Разделим отрезок интегрирования a,b
на четное число nравных от-
резков с шагами h. На каждом отрезке xi1 , xi1
длиной
2h, содержащем
три узла xi1, xi, xi1 , заменим подынтегральную функцию полиномом вто- рой степени ix ax2 bx c. Пример представлен на рис. 5.
Рис. 5. Пример замены функции
y fx параболой.
Для рассматриваемого примера значения коэффициентов a, bи cмо- гут быть вычислены следующим образом:
ah2 bh c
fi1 ,
c f,
i
ah2 bh c
fi1 ,
откуда c
fi,
2ah2 2 fi
fi1
fi1 , a
fi1 2 fi
2h2
fi1 , а b
fi1
2h
fi1 .
В результате рассматриваемый полином второй степени примет вид
fx i
x
fi
fi1
2h
fi1 x
fi1 2 fi
2h2
fi1 x2 .
Интегрируя приведенное выражение на отрезке xi1 , xi1 , получим
xi1
xi1
f f
f 2 f f 2 h
ixdx
fii1 i1 xi1 i i1 x
dx
fi1 4 fi fi1 .
(3.7)
xi1
xi1 2h
2h2 3
Приближенное значение интеграла на исходном отрезке интегрирова-
ния a,b
получим суммированием частичных интегралов (3.7) по всем от-
резкам xi1 , xi1 :
b
hn21
I f
a
xdx
3 i0
f2i
4 f2i1
f2i2
(3.8)
h f
3 0
4 f1
2 f2
4 f3
2 f4
... 2 f
n2
-
4 f
n1
fn.
Полученное соотношение называется формулой Симпсона или
формулой парабол, в соответствии с которой искомый определенный
интеграл вычисляется как суммарная площадь параболических сегментов на всех частичных отрезках интегрирования xi1 , xi1 .
Если подынтегральную функцию
fx
заменять полиномом второй
степени на отрезках xi1, xi,
i 1, n, с привлечением дополнительных точек
x xi xi1 ,
i 2
i 1, n
– середин данных отрезков, то число отрезков разбиения
nможет быть любым (не обязательно четным), а формула (4.8) будет иметь вид:
b hn
f xdx 6 fi1 4 fi
fi
a i1
h f 4 f 2 f 4 f
... 2 f
4 f f
(3.9)
6 0 1 1 2
n1 n n
Напомним, что
fi
fxi,
fi1
fxi1 ,
fi
fxi, i 1, n.
Замечание 1. Формула (3.8) может быть использована для вычисления интегралов от функций, заданных как аналитическим выражением, так и таблично, тогда как формула (3.9) применима только в тех случаях, когда подынтегральная функция задана аналитически.
Замечание 2. Рассмотренные формулы, используемые для приближенного вычисления интегралов, называются квадратурными формулами. Нетрудно заметить, что все рассмотренные выше формулы имеют следующую структуру:
b n
fxdx Aifi,
a i0
где
fi – значения подынтегральной функции в узловых точках
xi, а
Ai –
весовые коэффициенты.
- 1 ... 4 5 6 7 8 9 10 11 ... 15