Лекції: Наближене розв`язування інтегралів

1. Постановка задачі

2. Формули прямокутників

· лівих

· правих

· середніх

3. Формула трапеції

4. Формула Сімпсона

5. Формула Монте - Карло

6. Графічне інтегрування

Нехай функція f(x) задана на деякому відрізку [a;b]. Розглянемо задачу обчислення її визначеного інтеграла . Якщо для f(x) відома первісна F(x), то інтеграл може бути обчислений точно по основній формулі інтегрального числення – формулі Ньютона-Лейбніца

Однак первісна лише для вузького классу функцій виражається через елементарні функції, причому нерідко її відшукання зв`язане з дуже громісткими обчисленнями. Крім того, можлива ситуація, коли підінтегральна функція задана не аналітично, а таблично чи графічно. Тому для обчислення визначеного інтеграла часто приходиться застосовувати різні наближені формули. Досить просто ці формули можна одержати, виходячи з геометричного змісту визначеного інтегралу: якщо f(x) на [a;b], то - площа криволінійної трапеції, обмеженої відрізком [a;b] осі ох, кривою у=f(x) і прямими х=а, х=b. При наближеному обчисленні криволінійну трапецію заміняють фігурою, обмежену тим же відрізком [a;b], площа якої обчислюється значно простіше. Звідси одержують наближену формулу =S , яку використовують для обчислення лише тоді, коли можна оцінити її похибку.

В дальнішому розглядається функція f(x) для х [a;b], але одержані результати будуть вірні для будь-якої інтегруючої на [a;b] функції.

Ідея формули прямокутників заключається в тому, що на малому відрізку [x ;x +h] площа криволінійної трапеції приблизно рівна площі прямокутника з основою [x ;x +h] і висотою, що рівна ординаті в якійсь точці [x ;x +h], тобто

.

В залежності від того, яку точку відрізка [x ;x +h] вибирають в якості , і одержимо різновиди формули прямокутників.

Розіб`ємо відрізок [a;b] на n різних частин точками де h=

На кожному частинному відрізку [ ] замінимо відповідну криволінійну трапецію на прямокутник, висоту якого можна визначити по різному.

Формули лівих прямокутників. Якщо вибрати в якості висоти прямокутника на кожному з відрізків [ ] (і= ) ординату в лівому кінці, тобто , то криволінійна трапеція заміниться на ступіньчасту фігуру, площа якої можна прийняти за площу трапеції. Відповідно,

S , =S =h .

Одержану формулу називають формулою лівих прямокутників.

Ідея формули прямокутників заключається в тому, що на малому відрізку [x ;x +h] площа криволінійної трапеції приблизно рівна площі прямокутника з основою [x ;x +h] і висотою, що рівна ординаті в якійсь точці [x ;x +h], тобто

.

В залежності від того, яку точку відрізка [x ;x +h] вибирають в якості , і одержимо різновиди формули прямокутників.

Розіб`ємо відрізок [a;b] на n різних частин точками де h=

На кожному частинному відрізку [ ] замінимо відповідну криволінійну трапецію на прямокутник, висоту якого можна визначити по різному.

Формула правих прямокутників. Якщо в якості висоти прямокутника вибрати ординату правого кінця частинного відрізка, тобто висота прямокутника на відрізку [ ] рівна , i=( ), то для наближеного обчислення інтеграла одержимо формулу правих прямокутників: S , =S =h .

Ідея формули прямокутників заключається в тому, що на малому відрізку [x ;x +h] площа криволінійної трапеції приблизно рівна площі прямокутника з основою [x ;x +h] і висотою, що рівна ординаті в якійсь точці [x ;x +h], тобто

.

В залежності від того, яку точку відрізка [x ;x +h] вибирають в якості , і одержимо різновиди формули прямокутників.

Розіб`ємо відрізок [a;b] на n різних частин точками де h=

На кожному частинному відрізку [ ] замінимо відповідну криволінійну трапецію на прямокутник, висоту якого можна визначити по різному.

Формула середніх прямокутників. Висотою прямокутника, побудованого на частинному відрізку, вважають ординату, взяту в середній точці відрізка, тобто . Тоді S , =S =h .

Формула трапеції заснована на тому, що на відрізку [x ;x +h] дугу кривої у=f(x) заміняють хордою, що стягує кінці цієї дуги, тобто проводять лінійне інтерполювання функції у=f(x). При цьому площу криволінійної трапеції заміняють площою трапеції з основами f(x ) і f(x +h) і висотою h, відповідно, Звідси, щоб обчислити , розіб`ємо відрізок [a;b] на n рівних частинних відрізки, так що крок .

Проведемо через точки поділу ординати і замінимо дану криву ламаною, відрізки якої з`єднують кінці двух сусідніх ординат. Криволінійна трапеція при цьому заміниться на фігуру, що складається з n трапецій, висоти яких h, а основи – ординати , i=( ) в точках поділу. Площа такої фігури виражається формулою , а для обчислення інтеграла маємо формулу =S =h , яку називають формулою трапеції.

Метод Сімпсона .

Спосіб наближеного обчислення визначеного інтеграла по формулі Сімпсона заснований на тому, що на відрізку [x ;x +2h] дугу кривої у=f(x) заміняють дугою квадратичної параболи, що проходить через точки А(х ;f(x )), В(х +h;f(x +h)), C(х +2h;f(x +2h)), тобто проводять квадратичне інтерполювання функції у=f(x). Тоді за наближене значення площі криволінійної трапеції приймають площу параболічної трапеції, яка має ту саму основу [x ;x +2h] і обмежена зверху дугою параболи. Для того, щоб скласти рівняння цієї параболи, використаємо першу інтерполяційну формулу Н`ютона. Складемо многочлен другої степені по трьох вузлах інтерполяції: х , х =x +h, х =x +2h. Одержимо .

Перетворимо цей вираз. Так як х =x +h, то (х-x )*(х-х )=(х-x )*((х-x )-h)=(х-x ) -(х-x )*h i . Площа параболічної трапеції S

. Якщо ввести позначення - ордината початку відрізка, - ордината кінця відрізка, то одержана формула матиме вигляд S ) (1), де .

Розділимо тепер відрізок [a;b] на n их частин, причому ввазатимемо n-парним числом, тобто n=2m, тоді .

Нехай - точки поділу. Прведемо ординати в цих точках .З`єднаємо кінці кожних трьох сусідніх ординат дугами парабол, тобто на відрізках криву дугами парабол. Застосуємо до кожного з відрізків формулу (1). Тоді

, звівши подібні доданки, одержимо

(2)

Формулу (2) називають параболічною формулою чи формулою Сімпсона. Розглянемо задачу відшукання наближеного значення по формулі (2) і оцінимо похибку одержаного наближення. Для контролю задамо число Е>0 – точність шуканого наближення. Метод розв`язання задачі полягає в наступному. Значення S послідовно обчислюються по фомулі (2) при n=2,4,8,...,2 ,..., тобто обчислюються наближення . Після наступного обчислення (n>2) визначається похибка одержаного наближення: . Якщо , то наближення задовільняє заданій точності. Обчислення закінчене і = - наближене значення визначеного інтеграла, при заданому Е>0. В інакшому випадку обчислення повторюються.

Метод Монте-Карло.

Постановка задачі: ми маємо криволінійну трапецію, яка обмежена y=f(x), x=a, x=b і відрізком [a;b]. Необхідно знайти таку фігуру, яка повністю вміщає в себе криволінійну трапецію і достатньо просто знайти площу цієї фігури. Позначимо через S – площу фігури, яка включає криволінійну трапецію. Через S - площу криволінійної трапеції. Через N –кількість точок, які ми будемо вкидати у фігуру площі S. А через N - кількість точок, які попадуть у криволінійну трапецію. Тоді справедливе співвідношення . Звідси .

Зауваження 1: В оптимальному випадку відшуковується max значення функції y=f(x) на відрізку [a;b] і обчислюється площа прямокутника зі сторонами :

1) 1) сторона b-a,

2) 2) сторона max(f(x)).

Зауваження 2: Результат обчислення не залежить від вигляду і площі фігури, яка включає криволінійну трапецію.

Зауваження 3: Для складання програми необхідно використати генератор випадкових чисел. Перед цим необхідно активізувати його за допомогою стрічечки рандомайс.

Графічного інтегрування.

Наближене обчислення інтеграла методом графічного інтегрування застосовується тоді, коли підінтегральна функція задана графічно. Нехай на [a;b] задана неперервна крива, рівняння якої y=f(x). На основі теореми про середнє для визначеного інтеграла існує така точка [a;b], що

. Геометрично це означає, що площа криволінійної трапеції чисельно рівна площі прямокутника з основою [a;b] і висотою f( ). Розглянемо криволінійну трапецію і проведемо горизонтальну пряму приблизно так, щоб одержати потрібний прямокутник. Абсцисами точок перетину прямої і кривої будуть ті точки , про яких згадується у теоремі про середнє. Відкладемо на осі ох зліва від початку координат одиничний відрізок ОР і продовжимо проведену горизонтальну пряму до перетину з віссю ординат (якщо а<0, то краще спочатку зліва від а провести вертикальну пряму і при дальніших діях замінити вісь оу цією прямою). Нехай пряма перетне вісь оу в точці Q, тоді ОQ =f( ). З`єднаємо точки Р і Q і з точки а проведемо пряму аМ, паралельно РQ, до перетину в точці М з вертикаллю, проведеною з точки b. Покажемо, що bМ= , тобто величина побудованого відрізка чисельно рівна значенню визначеного інтеграла. Дійсно, . Звідси

Зауваження: Функція f(x)>0. Однак одержаний результат має місце для будь-якої неперервної на [a;b] функції y=f(x). Наприклад, функція f(x) міняє знак на [a;b]. Заштриховані площі на рисунку 2 приблизно рівні. Оскільки площа частини криволінійної трапеції, розміщеної нище осі ох, більша, то буде від`єним. Проводячи попередні побудови, одержимо відрізок bM, величина якого від`ємна і тут bM= .

На основі проведеної побудови і проводиться графічне інтегрування.