| ||||||||
МЕНЮ
| Наближене обчислення визначених інтегралів, що не беруться через елементарні функції /Укр./Наближене обчислення визначених інтегралів, що не беруться через елементарні функції /Укр./Наближене обчислення визначених інтегралів, що не беруться через елементарні функції ЗмістВступ. Формули прямокутників і трапеції. Параболічне інтерполювання. Дроблення проміжку. Залишковий член формули прямокутників. Залишковий член формули трапеції. Залишковий член формули Сімпсона. Додаток 1. Додаток 2. Висновки. Література. Вступ.Багато задач науки і техніки приводять до проблеми обчислення інтегралів, але не всі інтеграли піддаються обчисленню. В даній роботі разглядається питання наближеного обчислення визначених інтегралів, що не беруться через елементарні функції. Зокрема, виводяться формули наближеного обчислення прямокутників, формула трапецій а також формула Сімпсона. Формули прямокутників і трапеції. Нехай треба обчислити значення визначеного інтегралу В даній роботі можно ознайомитися з основними із цих методів, в яких наближені формули для інтегралів складаються по деякому числу значень підінтегральної функції, обчислених для ряду (зазвичай рівновіддалених) значень незалежної змінної. Перші формули, які сюди відносяться, простіші всього отримуються із геометричних міркувань. Витлумачуючи визначений інтеграл Перш за все, вдруге використовуючі ту думку, яка привела нас до самого поняття о визначеном інтегралі, можно розбити усю фігуру (мал. 1) на смуги, скажемо однієї і той же ширини
де На практиці зазвичай беруть
Надалі, кажучи про формулу прямокутників, ми будемо мати на увазі якраз цю формулу. Геометричні міркування природньо приводять і до другої, часто використовуваємій
наближеній формулі. Замінивши дану криву вписаною в неї ламаною, з вершинами
у точках
![]() Мал. 2 Додаючи, прийдемо до нової наближеної формули
Це так звана формула трапецій. Можно показати, що при зростанні Параболічне інтерполювання. Для наближеного обчислення інтеграла
і покласти Можно сказати, що тут – при обрахуванні площі – дана «
крива»
Сам вибір інтерполюючуго многочлена При інтерполюванні виходить лінійний, відносно значень В найпростішому випадку, при
Геометрично – площа криволінійної фігури замінюється тут площадью прямокутника з висотою, яка рівна середній її ординаті. При
і, як легко обчислити, Таким чином, тут ми наближено вважаємо На цей раз площа криволінійної фігури замінюється площею трапеції: замість кривої береться хорда, яка зполучає її кінці. Менш тривіальний результат отримаємо взявши
За допомогою легкого обчислення вираховуємо і, аналогічно
Таким чином, приходимо до наближеної формули
Тут площа фігури під даною кривою замінюється площею фігури, яка обмежена звичайною параболою (з вертикальною віссю), що проходить через крайні і середню точки кривої. Збільшуя степінь Дроблення проміжку. При обчисленні інтегралу
в зв' язку з чим, шуканий інтеграл постане у вигляді суми
Тепер же до кожного із цих проміжків застосуємо параболічне інтерполювання, тобто станемо обчислювати інтеграли (9) по одній із наближених формул – (4), (6), (8). Легко збагнути, що виходячи із формул (4) або (6), ми таким шляхом знов отримаємо вже відомі нам формули прямокутників і трапецій, (1) и (2). Застосуємо тепер до інтегралів (9) формулу (8), при цьому для стислості положимо, як і вище,
Ми отримаємо
. . . . . . . . . . . . . . . . . . . . . . . . . .
Зрештою, додаючи почленно ці равенства, прийдемо до формули
Вона носит назву формули Сімпсона (Th. Simpson); цією формулою користуються для наближенного обчислення інтегралів частіші, аніж формулами прямокутников і трапецій, бо она – при тих же затратах – дає зазвичай більш точний результат. Залишковий член формули прямокутників. Почнемо з формули (4). Припустимо, що у проміжку
де Якщо проінтегрувати цю рівність у проміжку від
Таким чином, отримаємо
так, що залишковий член формули (4), який поновлює її точність має вигляд
Позначив через
де
Якщо зараз розділити проміжок
Додавнши ці равенства (при
де вираз і є залишковий член формули прямокутників (1). Так як вираз також знаходиться між Тому остаточно маємо
При зростанні Залишковий член формули трапеції. Займемось тепер формулою (6) при попередніх здогатках відносно функції
Інтегруя цю формули від
так що залишковий член формули (6) буде
Розмірковуючи, як і вище, і користуючись тим, що другий множник підінтегральної функції і тут не змінює знака, знайдемо
Нарешті, для випадку ділення проміжку на
Таким є залишковий член формули трапецій (2). При зростанні Залишковий член формули Сімпсона. Звернемося, нарешті до формули (8). Можно було б, аналогічно тому, як це було зроблено тількі що, знов скористатись формулою Лагранжа з залишковим членом і покласти
Но ми стикаємося тут з таким станом речей, а саме, проінтегрувавши рівність (15), ми не змогли б спростити інтегральний вираз для додаткового члену за допомогою теореми про середне, бо вираз Вираз
яким би не було число
Тепер проінтегрувавши цю равність від так як
Якщо припустити похідну
користуючись тим, що другий множник в підінтергальному виразі не змінює знака, можно підставити в такому вигляді:
Якщо проміжок
При зростанні Текст программи для автоматичного обчислення інтегралів на мові програмування QBASIC: 'Тут описуються сталі e = 2.718281828459045# pi = 3.141592653589793# 'Тут задається від під інтегральної функції DEF fny# (x#) = e^ x# ^ 2 DEF fncoef# (i#) = (i# MOD 2) * 2 + 2 DEF fnxi# (i#) = a# + i# * h# DEF fnxis# (i#) = a# + i# * h# / 2 DEF fnxic# (i#) = a# + i# * h# + h# / 2 DEF fnxir# (i#) = a# + i# * h# + h# / 2 CLS 'Тут вводяться межі інтегрування та 'кількість проміжків INPUT “Введіть нижню межу інтегрування ” a# INPUT “Введіть верхню межу інтегрування ” b# INPUT “Введіть кількість проміжків ” n# 'Тут обчислюється крок h# = (b# - a#) / n# 'Тут обчислюється наближене значення 'інтеграла за методом Сімпсона integ# = 0 FOR i# = 1 TO ((2 * n#) - 1) integ# = integ# + fncoef#(i#) * fny#(fnxis#(i#)) NEXT integ# = integ# + fny#(a#) + fny#(b#) integ# = integ# * (h# / 6) PRINT "Simpson = "; integ# 'Тут обчислюється наближене значення 'інтеграла за методом трапецій integ# = 0 FOR i# = 1 TO (n# - 1) integ# = integ# + fny#(fnxi#(i#)) NEXT integ# = integ# + (fny#(a#) + fny#(b#)) / 2 integ# = integ# * h# PRINT " Trapeze = " ; integ# 'Тут обчислюється наближене значення 'інтеграла за методом лівих прямокутників integ# = 0 FOR i# = 0 TO (n# - 1) integ# = integ# + fny#(fnxi#(i#)) NEXT integ# = integ# * h# PRINT "L Rectangle = "; integ# 'Тут обчислюється наближене значення 'інтеграла за методом центральних прямокутників integ# = 0 FOR i# = 0 TO n# integ# = integ# + fny#(fnxic#(i#)) NEXT integ# = integ# * h# PRINT "C Rectangle = "; integ# 'Тут обчислюється наближене значення 'інтеграла за методом правих прямокутників integ# = 0 FOR i# = 1 TO n# integ# = integ# + fny#(fnxir#(i#)) NEXT integ# = integ# * h# PRINT "R Rectangle = "; integ# Додаток 2.Далі подані результати роботи програми, яка викладена в додатку 1. 1) n=1000 Метод Сімпсона -8.742278155181581D-08 Метод трапецій -8.742270585611512D-08 Метод лівих прямокутників 3.141505318306509D-03 Метод центральних прямокутників -3.14167628761223D-03 Метод правих прямокутників -6.283265152840917D-03 2) n=1000 Метод Сімпсона 2.000000000000067 Метод трапецій 1.999998355065565 Метод лівих прямокутників 1.999998355202888 Метод центральних прямокутників 1.999995887392223 Метод правих прямокутників 1.999990952591778 3)
|