Loading [MathJax]/jax/output/HTML-CSS/jax.js
Skip to main content
LibreTexts - Ukrayinska

5.5: Чисельна інтеграція

Фундаментальна теорема обчислення дає конкретну методику знаходження точного значення певного інтеграла. Ця методика заснована на обчисленні антипохідних. Незважаючи на силу цієї теореми, все ж існують ситуації, коли ми повинні наблизити значення певного інтеграла замість того, щоб знайти його точне значення. Перша ситуація, яку ми досліджуємо, полягає в тому, що ми не можемо обчислити антипохідне цілого. Другий випадок, коли ми насправді не знаємо цілісного, а тільки його значення при оцінці в певних точках.

Елементарна функція - це будь-яка функція, яка являє собою комбінацію поліномів,nth коренів, раціональних, експоненціальних, логарифмічних і тригонометричних функцій. Ми можемо обчислити похідну будь-якої елементарної функції, але є багато елементарних функцій, з яких ми не можемо обчислити антидериватив. Наприклад, такі функції не мають антипохідних, які ми можемо висловити елементарними функціями:

$ $ e^ {-x^2},\ квад\ sin (x^3)\ квад\ текст {і}\ квад\ розрив {\ sin x} {x}.\]

Найпростіший спосіб звернутися до антипохіднихex2 - це просто написатиex2 dx.

У цьому розділі окреслено три загальні методи наближення значення визначених інтегралів. Кожен з них описується як систематичний метод апроксимації площі під кривою. Апроксимуючи цю область точно, ми знаходимо точне наближення відповідного певного інтеграла.

Ми будемо застосовувати методи, які ми вивчаємо в цьому розділі, до наступних визначених інтегралів:

$\ int_0^1 e^ {-x^2}\ dx,\ quad\ int_ {-\ frac {\ pi} {4}} ^ {\ frac {\ pi} {2}}\ sin (x^3)\ dx,\ квад\ текст {і}\ quad\ int_ {0.5} ^ {4\ pi}\ frac {\ sin (x)} {x}\ dx,\]

як зображено на малюнку5.5.1.

альтальтальт

Малюнок5.5.1: Графічно представляють три визначені інтеграли, які не можуть бути оцінені за допомогою антипохідних.

Методи правила лівої та правої руки

У розділі 5.3 ми розглянули задачу оцінки визначених інтегралів шляхом наближення площі під кривою за допомогою прямокутників. Ми переглядаємо ці ідеї тут, перш ніж вводити інші методи наближення певних інтегралів.

Почнемо з огляду позначень. fДозволяти бути безперервна функція на інтервалі[a,b]. Бажаємо наблизитиbaf(x) dx. [a,b]Розбиваємо наn однаково розташовані підінтервали, кожен довжиниdx=ban. Кінцеві точки цих підінтервалів позначені як

$ х_1=а,\ x_2 = а+дх,\ x_3 = a+ 2dx,\\ ldpots,\ x_i = a+ (i-1)\ dx,\\ ldpots,\ x_ {n+1} = б.\]

Ключова ідея 8 стверджує, що для використання правила лівої руки ми використовуємо підсумовуванняni=1f(xi) dx і використовувати правило правої руки, яке ми використовуємоni=1f(xi+1) dx. Ми розглядаємо використання цих правил в контексті прикладів.

Приклад5.5.1: Approximating definite integrals with rectangles

Наближене10ex2 dx використання правил лівої та правої руки з 5 однаково розташованими підінтервалами.

Рішення

Починаємо з поділу інтервалу[0,1] на 5 однаково розташованих інтервалів. У нас єdx=105=1/5=0.2, так $x_1 = 0,\ x_2 = 0,2,\ x_3 = 0,4,\ x_4 = 0,6,\ x_5 = 0,8,\\ текст {і}\ x_6 = 1.\]

Використовуючи правило лівої руки, ми маємо:

ni=1f(xi) dx=(f(x1)+f(x2)+f(x3)+f(x4)+f(x5)) dx=(f(0)+f(0.2)+f(0.4)+f(0.6)+f(0.8)) dx(1+0.961+0.852+0.698+0.527)(0.2)0.808.

Використовуючи правило правої руки, ми маємо:

\ [\ почати {вирівняти}\ сума_ {i = 1} ^n f (x_ {i+1})\ dx &=\ великий (f (x_2) + f (x_3) + f (x_4) + f (x_5) +f (x_6)\ великий)\ dx\\ &=\ великий (f (0.2) + f (0.4) + f (0.4) + f (0.6) + f (0.6) + f (0.6) + f (0.6) + f (0,8) +f (1)\ великий)\ dx\\ &\ приблизно\ великий (0,961 +0,852 + 0,698 + 0,527 + 0,368) (0,2)\\\ &\ приблизно 0,681.
\ end {вирівняти}\]

альтальт

Малюнок5.5.2: Апроксимація10ex2 dx в прикладі5.5.1

5.5.2На малюнку показані прямокутники, які використовуються в кожному методі для наближення певного інтеграла. Ці графіки показують, що в даному конкретному випадку правило лівої руки є надмірним наближенням, а Правило правої руки - під наближенням. Щоб отримати краще наближення, ми могли б використовувати більше прямокутників, як ми це робили в розділі 3.1. Ми також могли б усереднити результати правила лівої та правої руки разом, даючи

$$\ гідророзриву {0.808 + 0.681} {2} = 0.7445.\]

Фактична відповідь, з точністю до 4 знаків після десяткового дробу, становить 0.7468, показуючи наше середнє значення є хорошим наближенням.

Приклад5.5.2: Approximating definite integrals with rectangles

Наближенеπ2π4sin(x3) dx використання правил лівої та правої руки з 10 однаково розташованими підінтервалами.

Рішення

Почнемо з знаходження\ Delta x:

$\ розрив {b-a} {n} =\ гідророзриву {\ pi/2 - (-\ pi/4)} {10} =\ гідророзриву {3\ pi} {40}\ приблизно 0,236.\]

Корисно виписати кінцеві точки підінтервалів в таблицю; в таблиці5.5.1 ми наведемо точні значення кінцевих точок, їх десяткові наближення та десяткові наближенняsin(x3) обчислюваних в цих точках.

Таблиця5.5.1: Таблиця значень, використовуваних дляπ2π4sin(x3) dx наближення в прикладі5.5.1.
xi Точний Прибл. sin(x3i)
\ (x_i\) ">x1 π/4 0.785 \ (\ sin (x_i^3)\) ">0.466
\ (x_i\) ">x2 7π/40 0.550 \ (\ sin (x_i^3)\) ">0.165
\ (x_i\) ">x3 π/10 0.314 \ (\ sin (x_i^3)\) ">0.031
\ (x_i\) ">x4 π/40 0.0785 \ (\ sin (x_i^3)\) ">0
\ (x_i\) ">x5 π/20 0.157 \ (\ sin (x_i^3)\) ">0.004
\ (x_i\) ">x6 π/8 0,393 \ (\ sin (x_i^3)\) ">0.061
\ (x_i\) ">x7 π/5 0,628 \ (\ sin (x_i^3)\) ">0.246
\ (x_i\) ">x8 11π/40 0,864 \ (\ sin (x_i^3)\) ">0.601
\ (x_i\) ">x9 7π/20 1.10 \ (\ sin (x_i^3)\) ">0.971
\ (x_i\) ">x10 17π/40 1.34 \ (\ sin (x_i^3)\) ">0.690
\ (x_i\) ">x11 π/2 1.57 \ (\ sin (x_i^3)\) ">0.670

Після створення цієї таблиці можна просто наблизити певний інтеграл за допомогою Правил лівої та правої руки. (Примітка: саму таблицю легко створити, особливо зі стандартною програмою електронних таблиць на комп'ютері. Останні дві колонки - це все, що потрібно.) Правило лівої руки підсумовує перші 10 значеньsin(x3i) і множить суму наdx; Правило правої руки підсумовує останні 10 значеньsin(x3i) і множить наdx. Тому у нас є:

Правило лівої руки:π2π4sin(x3) dx(1.91)(0.236)=0.451.

Правило правої руки:π2π4sin(x3) dx(1.71)(0.236)=0.404.

Середнє значення правил лівої та правої руки: 0.4275.

альтальт

Малюнок5.5.3:π2π4sin(x3) dx у прикладі5.5.2

Фактична відповідь, з точністю до 3 знаків після десяткового дробу, дорівнює 0,460. Наші наближення знову були досить хорошими. Прямокутники, використані в кожному наближенні, наведені на малюнку5.5.3. З графіків зрозуміло, що використання більшої кількості прямокутників (а значить, і більш вузьких прямокутників) повинно привести до більш точного наближення.

Трапецієподібне правило

У прикладі5.5.1 ми наблизили значення10ex2 dx з 5 прямокутників однакової ширини. 5.5.2На малюнку показані прямокутники, які використовуються в правилах лівої та правої руки. Ці графіки чітко показують, що прямокутники не так добре відповідають формі графіка, і що точні наближення прийдуть лише за допомогою безлічі прямокутників.

Замість того, щоб використовувати прямокутники для наближення площі, ми можемо замість цього використовувати трапеції. На малюнку5.5.4 ми показуємо область підf(x)=ex2 на[0,1] наближену 5 трапеціями однакової ширини; верхні «кути» кожної трапеції лежать на графікуf(x). З цієї цифри зрозуміло, що ці трапеції більш точно наближають площу підf і, отже, повинні давати краще наближення10ex2 dx. (Насправді ці трапеції, здається, дають велике наближення площі!)

альт

Малюнок5.5.4: Апроксимація10ex2 dx за допомогою 5 трапецій однакової ширини.

Формула для площі трапеції наведена на малюнку5.5.5. Ми10ex2 dx наближаємо з цими трапеціями в наступному прикладі.

Приклад5.5.3: Approximating definite integrals using trapezoids

Використовуйте 5 трапецій однакової ширини для наближення10ex2 dx.

Рішення

Для обчислення площ 5 трапецій на малюнку знову буде корисно створити таблицю значень5.5.6, як показано в табл5.5.2.

альт

Малюнок5.5.5: Площа трапеції

Таблиця5.5.2: Таблиця значеньex2.
xi ex2i
\ (x_i\) ">0 \ (e^ {-x_i^2}\) ">1
\ (x_i\) ">0.2 \ (e^ {-x_i^2}\) ">0.961
\ (x_i\) ">0.4 \ (e^ {-x_i^2}\) ">0.852
\ (x_i\) ">0.6 \ (e^ {-x_i^2}\) ">0.698
\ (x_i\) ">0.8 \ (e^ {-x_i^2}\) ">0.527
\ (x_i\) ">1 \ (e^ {-x_i^2}\) ">0.368

Крайня ліва трапеція має ніжки довжиною 1 і 0,961 і висотою 0,2. Таким чином, за нашою формулою площа крайньої лівої трапеції дорівнює:

$\ гідророзриву {1+0.961} {2} (0.2) = 0.1961.\]

Рухаючись вправо, наступна трапеція має ноги довжиною 0,961 і 0,852 і висотою 0,2. Таким чином, його площа становить:

$\ гідророзриву {0.961+0.852} 2 (0,2) = 0.1813.\]

Сума площ всіх 5 трапецій становить:

1+0.9612(0.2)+0.961+0.8522(0.2)+0.852+0.6982(0.2)+0.698+0.5272(0.2)+0.527+0.3682(0.2)=0.7445.

Ми наближаємо10ex2 dx0.7445.

У цьому прикладі можна спостерігати багато речей. Зверніть увагу, як кожен член в остаточному підсумовуванні множився як на 1/2, так і наdx=0.2. Ми можемо врахувати ці коефіцієнти, залишивши більш стисле підсумовування як:

$\ фрак12 (0.2)\ великий [(1+0.961) + (0.961+0.852) + (0.852+0.698) + (0.698+ 0,527) + (0.527 + 0,368)\ Великий].\]

Тепер зверніть увагу, що всі числа, крім першого і останнього, додаються двічі. Тому ми можемо написати підсумовування ще більш лаконічно, як

$\ гідророзрив {0.2} {2}\ Великий [1 + 2 (0.961+0.852+0.698+0.527) + 0,368\ Великий].\]

Це серце трапецієподібного правила, де певний інтегралbaf(x) dx наближається за допомогою трапецій рівної ширини для наближення відповідної площі підf. Використовуючиn однаково розташовані підінтервали з кінцевими точкамиx1x2,,xn+1,, ми знову маємоΔx=ban. Таким чином:

baf(x) dxni=1f(xi)+f(xi+1)2 dx= dx2ni=1(f(xi)+f(xi+1))= dx2[f(x1)+2ni=2f(xi)+f(xn+1)].

Приклад5.5.4: Using the Trapezoidal Rule

Перегляньте приклад5.5.4 та наблизитесьπ2π4sin(x3) dx за допомогою трапецієподібного правила та 10 однаково розташованих підінтервалів.

Рішення

Повертаємося до Таблиці5.5.1 для таблиці значеньsin(x3). Нагадаємо, щоdx=3π/400.236. Таким чином, ми маємо:

π2π4sin(x3) dx0.2362[0.466+2(0.165+(0.031)++0.69)+(0.67)]=0.4275.

Зверніть увагу, як «швидко» трапецієподібне правило може бути реалізовано після створення таблиці значень. Це справедливо для всіх методів, досліджених у цьому розділі; реальною роботою є створення таблиціxi таf(xi) значень. Як тільки це буде завершено, наблизити певний інтеграл не складе труднощів. Знову ж таки, використовувати технологію розумно. Електронні таблиці можуть зробити швидку роботу цих обчислень і полегшити використання безлічі субінтервалів.

Також зверніть увагу на наближення, які дає трапецієподібне правило. Це середнє значення наближень, заданих Правилами лівої та Правої Руки! Це ефективно робить правила лівої та правої руки застарілими. Вони корисні при першому вивченні певних інтегралів, але якщо потрібне реальне наближення, то, як правило, краще використовувати трапецієподібне правило замість правила лівої або правої руки.

Як ми можемо вдосконалити трапецієподібне правило, окрім використання все більшої кількості трапецій? Відповідь зрозуміла, як тільки ми озирнемося назад і розглянемо, що ми насправді зробили досі. Правило лівої руки насправді не стосується використання прямокутників для наближення площі. Натомість він наближає функціюf з постійними функціями на малих підінтервалах, а потім обчислює певний інтеграл цих постійних функцій. Трапецієподібне правило дійсно апроксимує функціюf з лінійною функцією на малому підінтервалі, потім обчислює певний інтеграл цієї лінійної функції. В обох цих випадках певні інтеграли легко обчислити в геометричному вираженні.

Отже, ми маємо прогресію: ми починаємо з наближенняf з постійною функцією, а потім з лінійною функцією. Що далі? Квадратична функція. Апроксимуючи криву функції з безліччю парабол, ми, як правило, отримуємо ще краще наближення певного інтеграла. Ми називаємо цей процес Правилом Сімпсона, названим на честь Томаса Сімпсона (1710-1761), хоча інші використовували це правило аж 100 років до цього.

Правило Сімпсона

З огляду на один момент, ми можемо створити постійну функцію, яка проходить через цю точку. З огляду на дві точки, ми можемо створити лінійну функцію, яка проходить через ці точки. Враховуючи три точки, ми можемо створити квадратичну функцію, яка проходить через ці три точки (враховуючи, що жодна з двох не має однаковогоx —value).

Розглянемо три точки(x1,y1),(x2,y2) і(x3,y3) чиїx —значення однаково розташовані іx1<x2<x3. fДозволяти квадратичної функції, яка проходить через ці три точки. Це не важко показати, що

x3x1f(x) dx=x3x16(y1+4y2+y3).

Розглянемо Малюнок5.5.6. Функціяf проходить через показані 3 точки, а параболаg, яка також проходить через ці точки, графічна пунктирною лінією. Використовуючи наше рівняння зверху, ми точно знаємо, що  int31g(x) dx= frac316 big(3+4(1)+2 big)=3. Оскількиg є хорошим наближенням дляf включення[1,3], ми можемо констатувати, що $$\ int_1^3 f (x)\ dx\ приблизно 3.

альт

Малюнок5.5.6: Графік функціїf та параболи, яка добре її наближає[1,3].

Зверніть увагу, як інтервал[1,3] був розділений на два підінтервали, оскільки нам потрібно було 3 бали. Через це, всякий раз, коли ми використовуємо Правило Сімпсона, нам потрібно розбити інтервал на парну кількість підінтервалів.

Загалом, щоб наблизитиbaf(x) dx за допомогою Правила Сімпсона,[a,b] поділіть наn підінтервали, деn парний і кожен підінтервал має ширину dx=(ba)/n. Миf наближаємосяn/2 параболічними кривими, використовуючи Рівняння5.5.22 для обчислення площі під цими параболами. Складання цих областей дає формулу:

$\ int_a^b f (x)\ dx\ приблизно\ розрив {\ dx} 3\ великий [f (x_1) +4f (x_2) +2f (x_3) +4f (x_4) +\ ldots+2f (x_ {n-1}) +4f (x_n) +f (x_ {n+1})\ великий]\.]

Зверніть увагу, як коефіцієнти членів при підсумовуванні мають візерунок 1, 4, 2, 4, 2, 4, 2, 4, 1.

Давайте продемонструємо Правило Сімпсона на конкретному прикладі.

Приклад5.5.5: Using Simpson's Rule

Приблизно10ex2 dx за допомогою правила Сімпсона та 4 однаково розташованих підінтервалів.

Рішення

Ми починаємо з складання таблиці значень, як ми мали в минулому, як показано в табл5.5.3.

Таблиця5.5.3: Таблиця значень для10ex2 dx наближення разом з графіком функції.
xi ex2i
\ (x_i\) ">0 \ (e^ {-x_i^2}\) ">1
\ (x_i\) ">0,25 \ (e^ {-x_i^2}\) ">0.939
\ (x_i\) ">0.5 \ (e^ {-x_i^2}\) ">0.779
\ (x_i\) ">0,75 \ (e^ {-x_i^2}\) ">0.570
\ (x_i\) ">1 \ (e^ {-x_i^2}\) ">0.368

Правило Сімпсона стверджує, що

$\ int_0^1e^ {-x^2}\ dx\ приблизно\ розрив {0,25} {3}\ Великий [1+4 (0.939) +2 (0.779) +4 (0.570) + 0,368\ Великий] = 0.7468\ overline {3}.\]

Нагадаємо, в прикладі5.5.1 ми заявили, що правильна відповідь, з точністю до 4 знаків після десяткового дробу, склала 0.7468. Наше наближення з правилом Сімпсона, з 4 підінтервалами, краще, ніж наше наближення з трапецієподібним правилом, використовуючи 5!

альт

Малюнок5.5.7: Використання правила Сімпсона зn=4 для наближення10ex2 dx.

Малюнок5.5.7 показуєf(x)=ex2 разом зі своїми апроксимуючими параболами, демонструючи, наскільки добре наше наближення. Апроксимуючі криві майже не відрізняються від фактичної функції.

Приклад5.5.6: Using Simpson's Rule

Приблизноπ2π4sin(x3) dx за допомогою правила Сімпсона та 10 однаково розставлених інтервалів.

Рішення

Таблиця5.5.4 показує таблицю значень, які ми використовували в минулому для цієї проблеми, показана тут знову для зручності. Знову ж таки, dx=(π/2+π/4)/100.236.

Таблиця5.5.4: Таблиця значень, використовуваних дляπ2π4sin(x3) dx наближення в прикладі5.5.6.
xi sin(x3i)
\ (x_i\) ">0.785 \ (\ sin (x_i^3)\) ">0.466
\ (x_i\) ">0.550 \ (\ sin (x_i^3)\) ">0.165
\ (x_i\) ">0.314 \ (\ sin (x_i^3)\) ">0.031
\ (x_i\) ">0.0785 \ (\ sin (x_i^3)\) ">0
\ (x_i\) ">0.157 \ (\ sin (x_i^3)\) ">0.004
\ (x_i\) ">0,393 \ (\ sin (x_i^3)\) ">0.061
\ (x_i\) ">0.628 \ (\ sin (x_i^3)\) ">0.246
\ (x_i\) ">0.864 \ (\ sin (x_i^3)\) ">0.601
\ (x_i\) ">1.10 \ (\ sin (x_i^3)\) ">0.971
\ (x_i\) ">1.34 \ (\ sin (x_i^3)\) ">0.690
\ (x_i\) ">1.57 \ (\ sin (x_i^3)\) ">0.670

Правило Сімпсона стверджує, що

π2π4sin(x3) dx0.2363[(0.466)+4(0.165)+2(0.031)++2(0.971)+4(0.69)+(0.67)]=0.4701

альт

Малюнок5.5.8: Наближенняπ2π4sin(x3) dx у прикладі5.5.6 з правилом Сімпсона та 10 однаково розташованими інтервалами.

Нагадаємо, що фактичне значення, з точністю до 3 знаків після коми, дорівнює 0,460. Наше наближення знаходиться в межах однієї 1/100th від правильного значення. Графік на малюнку5.5.8 показує, наскільки точно параболи відповідають формі графіка.

Резюме та аналіз помилок

Ми узагальнюємо ключові поняття цього розділу на сьогоднішній день у наступній Key Idea.

Ключова ідея 9: Чисельне інтегрування

fДозволяти бути безперервна функція на[a,b], нехайn буде натуральне число, і нехайΔx=ban.

Набірx1=a,x2=a+ dx,,xi=a+(i1) dx,xn+1=b.

Розглянемоbaf(x) dx.

Правило лівої руки:baf(x) dx dx[f(x1)+f(x2)++f(xn)].

Правило правої руки:baf(x) dx dx[f(x2)+f(x3)++f(xn+1)].

Трапецієподібне правило:baf(x) dx dx2[f(x1)+2f(x2)+2f(x3)++2f(xn)+f(xn+1)].

Правило Сімпсона:baf(x) dx dx3[f(x1)+4f(x2)+2f(x3)++4f(xn)+f(xn+1)] (nпарне)}.

У наших прикладах ми наблизили значення певного інтеграла за допомогою заданого методу, а потім порівняли його з «правильною» відповіддю. Це повинно було викликати кілька питань у свідомості читача, таких як:

  1. Як обчислювалася «правильна» відповідь?
  2. Якщо вдасться знайти правильну відповідь, який сенс наближення?
  3. Якщо є значення для наближення, як ми повинні знати, чи є наближення будь-яким хорошим?

Це хороші питання, і відповіді на них навчальні. У прикладах правильну відповідь так і не обчислювався. Швидше за все, було дано наближення з точністю до певної кількості знаків після десяткового числа. У прикладі5.5.1 ми не знаємо точної відповіді, але ми знаємо, що вона починається з 0.7468. Ці більш точні наближення обчислювалися за допомогою числового інтегрування, але з більшою точністю (тобто більше субінтервалів та допомогою комп'ютера).

Оскільки точної відповіді знайти не вдається, наближення все ще має своє місце. Як ми можемо визначити, чи є наближення будь-яким хорошим?

«Проби і помилки» передбачає один спосіб. Використовуючи технологію, зробіть наближення з, скажімо, 10, 100 і 200 субінтервалами. Це, швидше за все, не займе багато часу, і повинна з'явитися тенденція. Якщо тенденція не виникає, спробуйте використовувати ще більше субінтервалів. Майте на увазі, що спроби та помилки ніколи не є безглуздими; ви можете наштовхнутися на проблему, в якій тенденція не з'явиться.

Другий спосіб полягає у використанні аналізу помилок. Хоча деталі виходять за рамки цього тексту, є деякі формули, які дають межі того, наскільки добре буде ваше наближення. Наприклад, формула може стверджувати, що наближення знаходиться в межах 0,1 від правильної відповіді. Якщо наближення дорівнює 1,58, то відомо, що правильна відповідь знаходиться між 1,48 і 1,68. Використовуючи безліч підінтервалів, можна отримати наближення настільки точно, як подобається. Теорема5.5.1 стверджує, що це за межі.

Теорема5.5.1: Error Bounds in the Trapezoidal and Simpson's Rules

  1. ETДозволяти помилка в наближенні заbaf(x) dx допомогою трапецієподібного правила.
    Якщоf має безперервнуnd похідну 2 на[a,b] іM є будь-якою верхньою межею|f"(x)| on[a,b], то

$$ E_T\ лек\ гідророзриву {(b-a) ^3} {12n^2} М.\]

  1. ESДозволяти бути помилкою в наближенніbaf(x) dx за допомогою правила Сімпсона.
    Якщоf має безперервнуth похідну 4 на[a,b] іM є будь-якою верхньою межею|f(4)| on[a,b], то

$E_S\ лек\ гідророзриву {(б-а) ^5} {180n^4} М.\]

Є деякі ключові речі, які слід зазначити щодо цієї теореми.

  1. Чим більше інтервал, тим більше похибка. Це повинно мати сенс інтуїтивно.
  2. Помилка зменшується, оскільки використовується більше субінтервалів (тобто, якn стає більшим).
  3. Помилка в Правилі Сімпсона має термін, що відноситься до 4th похідної відf. Розглянемо кубічний многочлен: його4th похідна дорівнює 0. Отже, похибка наближення певного інтеграла кубічного полінома правилом Сімпсона дорівнює 0 — Правило Сімпсона обчислює точну відповідь!

Ми переглядаємо Приклади5.5.35.5.5 та обчислюємо межі помилок, використовуючи теорему5.5.1 в наступному прикладі.

Приклад5.5.7: Computing error bounds

Знайдіть межі помилок під час апроксимації10ex2 dx за допомогою трапецієподібного правила та 5 підінтервалів та використання правила Сімпсона з 4 підінтервалами.

Рішення

Трапецієподібне правило зn=5:

Почнемо з обчислення2nd похідної відf(x)=ex2:

f"(x)=ex2(4x22).

5.5.8На малюнку зображений графікf"(x) включення[0,1]. Зрозуміло, що найбільше значенняf", в абсолютному значенні, дорівнює 2. Таким чином, ми дозволяємоM=2 і застосовуємо формулу помилки з теореми5.5.1.

альт

Малюнок5.5.9: Графікf"(x) у прикладі,5.5.7 щоб допомогти встановити межі помилок.

$ $ E_T =\ frac {(1-0) ^3} {12\ cdot 5 ^ 2}\ cdot 2 = 0.00\ оверлайн {6}.\]

Наша формула оцінки помилок стверджує, що наше наближення 0.7445, знайдене в прикладі,5.5.3 знаходиться в межах 0,0067 від правильної відповіді, отже, ми знаємо, що

$0,7445-0,0067 = .7378\ leq\ int_0^1e^ {-x^2}\ dx\ лек 0.7512 = 0.7445 + 0,0067.\]

Раніше ми обчислювали точну відповідь, правильну до 4 знаків після коми, щоб бути 0.7468, підтверджуючи дійсність теореми5.5.1.

Правило Сімпсона зn=4:

Почнемо з обчислення4th похідної відf(x)=ex2:

$ $ ф\, ^ {(4)} (х) = e^ {-х^2} (16х^4-48х^2+12).\]

5.5.9На малюнку зображений графікf(4)(x) включення[0,1]. Зрозуміло, що найбільше значенняf(4), в абсолютному значенні, дорівнює 12. Таким чином, ми дозволяємоM=12 і застосовуємо формулу помилки з теореми5.5.1.

$e_s =\ frac {(1-0) ^5} {180\ точка 4 ^ 4}\ точка 12 = 0,00026.\]

альт

Малюнок5.5.10: Графікf(4)(x) у прикладі,5.5.7 щоб допомогти встановити межі помилок.

Наша формула оцінки помилок стверджує, що наше наближення0.7468¯3 знайденого в прикладі5.5.5 знаходиться в межах 0.00026 від правильної відповіді, отже, ми знаємо, що

$0.74683-0.00026 = .74657\ leq\ int_0^1e^ {-x^2}\ dx\ лек 0.74709 = 0.74683 + 0.00026.\]

Ще раз підтверджуємо обґрунтованість теореми5.5.1.

На початку цього розділу ми згадали дві основні ситуації, коли числове інтегрування було бажаним. Ми розглянули випадок, коли антипохідне цілого числа неможливо обчислити. Ми зараз досліджуємо ситуацію, коли цілісність невідома. Це, по суті, найбільш широко використовуване застосування методів числового інтегрування. «Більшу частину часу» ми спостерігаємо поведінку, але не знаємо «функції», яка її описує. Замість цього ми збираємо дані про поведінку і робимо наближення на основі цих даних. Ми демонструємо це на прикладі.

Приклад5.5.8: Approximating distance traveled

Один з авторів відвозив дочку додому зі школи, поки вона записувала їх швидкість кожні 30 секунд. Дані наведені в табл5.5.5. Приблизна відстань, яку вони пройшли.

Рішення

Нагадаємо, що інтегруючи функцію швидкості, ми отримуємо пройдену відстань. У нас є інформація проv(t); ми будемо використовувати Правило Сімпсона для наближенняbav(t) dt.

Найскладнішим аспектом цієї проблеми є перетворення даних у форму, в якій ми повинні бути. Швидкість вимірюється в милі на годину, тоді як час вимірюється з кроком 30 секунд.

Таблиця5.5.5: Дані про швидкість, зібрані з інтервалом 30 секунд, наприклад5.5.8.
0 0
1 25
2 22
3 19
4 39
5 0
6 43
7 59
8 54
9 51
10 43
11 35
12 40
13 43
14 30
15 0
16 0
17 28
18 40
19 42
20 40
21 39
22 40
23 23
24 0

Нам потрібно обчислити dx=(ba)/n. Зрозуміло,n=24. Що такеa іb? Оскільки ми починаємо вчасноt=0, ми маємо цеa=0. Остаточний записаний час настав через 24 періоди 30 секунд, що становить 12 хвилин або 1/5 години. Таким чином, ми маємо

$\ dx =\ гідророзриву {b-a} {n} =\ гідророзриву {1/5-0} {24} =\ фрак1 {120};\ квад\ гідророзриву {\ dx} {3} =\ гідророзриву {1} {360}.\]

Таким чином пройдена відстань становить приблизно:

\ [\ почати {align*}\ int_0^ {0.2} v (t)\ dt &\ приблизно\ розрив {1} {360}\ Великий [f (x_1) +4f (x_2) + 2f (x_3) +\ cdots + 4f (x_n) +f (x_ {n+1})\ великий]\\
&=\ frac {1} {360}\ Великий [0+4\ cdot25+2\ cdot 2 +\ cdot+ 2\ cdot40+4\ dot 23 + 0\ Великий]\\ &\ приблизно 6.2167\\ текст {милі.} \ end {вирівнювати*}\]

Ми приблизні, автор проїхав 6,2 км. (Оскільки ми впевнені, що читач хоче знати, одометр автора записав відстань приблизно 6,05 миль.)

Ми розпочали цю главу, вивчаючи антипохідні та невизначені інтеграли. Потім ми, здавалося, змінити фокус, дивлячись на області між графіком функції іx -вісь. Ці області ми визначили як певний інтеграл функції, використовуючи позначення, дуже схоже на позначення невизначеного інтеграла. Фундаментальна теорема обчислення пов'язала ці два, здавалося б, окремих поняття воєдино: ми можемо знайти ділянки під кривою, тобто ми можемо оцінити певний інтеграл, використовуючи антипохідні.

Ми закінчили главу, зазначивши, що антипохідні іноді знайти більш ніж важко: вони неможливі. Тому ми розробили числові методи, які дали нам хороші наближення визначених інтегралів.

Ми використовували певний інтеграл для обчислення площ, а також для обчислення переміщень і пройдених відстаней. Ми можемо зробити набагато більше, ніж це. У розділі 7 ми побачимо більше додатків певного інтеграла. Перед цим у розділі 6 ми вивчимо передові методи інтеграції, аналогічні правилам навчання, таким як правила диференціації продукту, частки та ланцюга.

Дописувачі та атрибуція