Трапезоиден метод во интеграли
Во применетата математика, често се среќаваме со проблемот на пресметување на површината под крива. Теоретски, ова може да се реши со користење на дефинирани интеграли. Меѓутоа, во пракса, не сите функции лесно се интегрираат аналитички. Некои функции се сложени, податоците се достапни само во табеларна форма или на моделите им недостасуваат едноставни антидеривати. Во вакви ситуации, нумеричките методи стануваат неопходни. Еден од најпознатите и најчесто користените нумерички методи е трапезоидното правило, техника на апроксимација за интеграли што ја дели површината под крива на трапезоидни форми.
Основен концепт на определен интеграл
Геометриски, дефинитивниот интеграл (\int_a^bf(x)\,dx\) може да се разбере како површина ограничена со кривата \(y=f(x)\), оската \(x\) и вертикалните линии \(x=a\) и \(x=b\). Ако \(f(x)\ge 0\), интегралот е позитивна површина. Ако функцијата има негативни вредности на одреден интервал, интегралот дава „означена површина“.
Главниот проблем се јавува кога:
1. Функцијата \(f(x)\) нема антидериват што може да се изрази со елементарни функции.
2. Вредноста на \(f(x)\) е позната само во одредени точки (експериментални податоци).
3. Симболичкото пресметување е неефикасно или невозможно.
Тука се користи трапезоидниот метод за приближна пресметка на интегралната вредност со користење на вредностите на функцијата во одредени точки.
Идеја за трапезоиден метод
Трапезоидниот метод започнува со едноставна идеја: приближете крива со права линија на секој подинтервал. Додека методот на Риманово збирно пресметување користи нормални страни и правоаголна површина, трапезоидниот метод користи две крајни точки на подинтервалот и ги поврзува со права линија. Површината под правата линија формира трапез, а не правоаголник.
Да претпоставиме дека сакаме да пресметаме:
\[
\int_a^bf(x)\,dx
\]
Го делиме интервалот \([a,b]\) на \(n\) подинтервали со еднаква должина. Должината на секој подинтервал е:
\[
h=\frac{ba}{n}
\]
Поделбени точки:
\[
x_0=a,\;x_1=a+h,\;x_2=a+2h,\;\точки,\;x_n=b
\]
Со вредности на функцијата:
\[
f(x_0), f(x_1), \точки, f(x_n)
\]
На секој подинтервал \([x_{i-1}, x_i]\), површината под кривата е приближно пресметана со површината на трапезоидот:
\[
A_i приближно \frac{h}{2}\лево[f(x_{i-1}) + f(x_i)\десно]
\]
Потоа вкупната површина (интегрална апроксимација) се добива со собирање на сите трапезоиди:
\[
\int_a^bf(x)\,dx \approx \sum_{i=1}^{n}\frac{h}{2}\left[f(x_{i-1}) + f(x_i)\right]
\]
Ако се среди, формулата за сложено трапезоидно правило станува:
\[
\int_a^b f(x)\,dx \approx \frac{h}{2}\left[f(x_0)+2f(x_1)+2f(x_2)+\cdots+2f(x_{n-1})+f(x_n)\right]
\]
Забележете го моделот на коефициенти: крајните точки \(x_0\) и \(x_n\) се множат со 1, додека точките во средината се множат со 2.
Едноставен пример за пресметка
Да претпоставиме дека сакаме да пристапиме кон:
\[
\int_0^2 x^2\,dx
\]
Аналитички, интегралот е \(\left.\frac{x^3}{3}\right|_0^2=\frac{8}{3}\approx 2{,}6667\). Сепак, ќе го користиме трапезоидниот метод со \(n=4\).
1. Интервал \([0,2]\), потоа:
\[
h=\frac{2-0}{4}=0{,}5
\]
2. Titik-titik: \(x_0=0\), \(x_1=0{,}5\), \(x_2=1\), \(x_3=1{,}5\), \(x_4=2\)
3. Вредност на функцијата:
\(f(0)=0\)
\(f(0{,}5)=0{,}25\)
\(f(1)=1\)
\(f(1{,}5)=2{,}25\)
\(f(2)=4\)
Користете ја формулата за сложен трапез:
\[
T=\frac{h}{2}\left[f(x_0)+2f(x_1)+2f(x_2)+2f(x_3)+f(x_4)\right]
\]
\[
T=\frac{0{,}5}{2}\left[0+2(0{,}25)+2(1)+2(2{,}25)+4\right]
\]
\[
T=0{,}25\left[0{,}5+2+4{,}5+4\right]=0{,}25(11)=2{,}75
\]
Приближниот резултат од 2,75 е доста блиску до точната вредност од 2,6667, со грешка од околу 0,0833.
Грешка во трапезоидниот метод
Трапезоидниот метод е апроксимациски метод, што значи дека секогаш постои разлика помеѓу приближната вредност и вистинската интегрална вредност. Големината на грешката е под влијание на:
1. Број на подинтервали (n): колку е поголем (n), толку е помал (h) и генерално толку се поточни резултатите.
2. Закривеност на функцијата: ако функцијата е многу закривена (има голем втор извод), трапезоидниот метод може да биде помалку точен од методите од повисок ред.
За доволно мазни функции (кои имаат континуирани втори изводи), грешката на сложениот трапезоиден метод може да се процени:
\[
E_T = -\frac{(ba)}{12}h^2 f”(\xi)
\]
за некој \(\xi\) помеѓу \(a\) и \(b\). Од оваа формула може да се види дека грешката е пропорционална на \(h^2\). Тоа е, ако го преполовиме \(h\) (т.е. го дуплираме бројот на подинтервали), грешката се намалува за приближно една четвртина.
Предности на трапезоидниот метод
Трапезоидниот метод е популарен од неколку причини:
1. Едноставно и брзо
Пресметката е релативно лесна, па дури може да се направи и рачно за мали случаи.
2. Погодно за табеларни податоци
Ако вредноста на \(f(x)\) е позната само во дискретни точки (на пр. резултати од мерењето), овој метод може да се примени директно.
3. Подобра точност од правоаголниот метод
Бидејќи користи две крајни точки и ја апроксимира функцијата со права линија, генерално е попрецизна од правоаголната апроксимација за ист број на партиции.
4. Стабилен за многу инженерски апликации
Многу физички и инженерски проблеми (на пр. пресметување на работа, енергија, површина на пресек, празнење и така натаму) може да се решат со овој метод.
Ограничувања на трапезоидниот метод
Иако е корисен, овој метод има ограничувања:
1. Помалку прецизно за високо закривени функции
Бидејќи користи само линеарна апроксимација на секој подинтервал, функциите со брзи промени или остри криви бараат голем \(n\).
2. Не е метод од висок ред
Грешки од редот на \(h^2\); други методи како што е Симпсоновиот може да дадат помали грешки во ист број на подинтервали (под услов функцијата да е доволно мазна).
3. Чувствителен на избор на интервал
Ако интервалот е преширок, а \(n\) е мал, резултатите можат значително да отстапуваат.
Затворање
Трапезоидниот метод на интеграција е една од најфундаменталните техники во нумеричката интеграција. Со делење на интервал на мали делови и замена на кривата со права линија на секој дел, можеме да го приближиме интегралот без да мора аналитички да го наоѓаме антидериватот. Формулата е едноставна, лесна за имплементација и особено корисна кога податоците се дискретни или функцијата е тешка за интегрирање.
Сепак, корисниците на трапезоидниот метод треба да обрнат внимание на бројот на подинтервали и однесувањето на функцијата што се интегрира. За многу закривени функции или кога е потребна голема прецизност, трапезоидниот метод сè уште може да се користи со зголемување на \(n\), или може да се земат предвид други нумерички методи како што се Симпсоновите или Ромберговите методи. Со добро разбирање, трапезоидниот метод станува ефикасна и практична алатка за решавање на разни интегрални проблеми во науката, инженерството и економијата.
