שיטות נומריות לחישוב אינטגרלים מסוימים

בעוד שחלק מהאינטגרלים ניתן לחשב באמצעות שיטות אנליטיות הרי יש אינטגרלים אותם לא ניתן לחשב בצורה כזאת אלא רק עם אנליזה נומרית, כלומר התוצאה המתקבלת היא מספר מסוים המדויק בכל רמת דיוק שנרצה.

להלן רשימה חלקית של שיטות לביצוע תהליך זה:

שימוש בהגדרת האינטגרל לפי רימן

לפי הגדרת האינטגרל לפי רימן מתקיים: $$\underset{a}{\overset{b}{\int }}f(x)\mathrm{d}x=\underset{n\to \mathrm{\infty }}{\lim }{\sigma }_n(f)=\underset{n\to \mathrm{\infty }}{\lim }\sum _{i=1}^{n}f({\xi }_i)\cdot (x_i-x_{i-1})$$ כאשר ${\xi }_i\in [x_{i-1},x_i]$ וההפרש בין נקודות סמוכות שואף לאפס כאשר $n$ שואף לאינסוף.

באנליזה נומרית לא ניתן אמנם להגיע לטור אינסופי אבל ניתן להשתמש בנוסחת הסכום על מנת להגיע לדיוק הנדרש. דבר זה נעשה על ידי שימוש בסכום: $\sum _{i=1}^n{w}_i\cdot f({\xi }_i)\cdot (x_i-x_{i-1})$ (כאשר $w_i$ הוא המשקל של הנקודה) כקירוב לאינטגרל. ניתן להוכיח כי בתנאים מסוימים (השונים במקצת בין השיטות אולם בעיקרם קשורים לרציפות ולאינטגרביליות של הפונקציה) הגדלת מספר הנקודות מביא לקירוב טוב יותר (מידת הקירוב המדויקת תלויה בשיטה).

את הנקודות ואת המשקלים קובעים על פי השיטה המדויקת, כאשר ישנו מספר רב של שיטות מסוג זה כמו שיטת הטרפז, המקרבת את השטח הכלוא מתחת לגרף הפונקציה לאוסף של טרפזים (ראו הדוגמה להלן), שיטת סימפסון (המקרבת את גרף הפונקציה לאוסף של פרבולות), שיטת רומברג (שיטה העושה שימוש בשיטת הטרפז עם מרווחים הולכים וקטנים), שיטת גאוס (עושה שימוש במקדמים המחושבים על פי פולינומים אורתוגונליים). כמו כן השיטות מתחלקות לשיטות פתוחות (בלי התחשבות בנקודות הקצה) ולשיטות סגורות (המתחשבות בכל הנקודות).

שיטת הטרפז עושה שימוש במשקלות של $w_i=\frac{1}{2}$ כאשר המרחקים קבועים בין הנקודות, כלומר: $$\underset{a}{\overset{b}{\int }}f(x)\mathrm{d}x\approx \frac{b-a}{n}(\frac{f(a)+f(b)}{2}+\sum _{i=1}^{n-1}f(a+\frac{b-a}{n}i))$$ כאשר ניתן לשפר את הדיוק באמצעות הגדלת מספר הנקודות.

שיטת מונטה קרלו

תבנית:הפניה לערך מורחב שיטת מונטה קרלו היא שיטת חישוב סטטיסטית שניתן לתאר בקיצור כך: נצייר את צורת הגרף על משטח בעל גודל ידוע ונסמן עליו "פגיעות" באופן אקראי. באמצעות חישוב מספר הפגיעות בתוך התחום החסום על ידי הגרף חלקי מספר הפגיעות מחוץ לו ניתן לחשב את גודל אותו התחום. באמצעות מחשב ניתן לחשב כך תוצאות מדויקות באופן מפתיע.

שימוש באינטגרלים ידועים

ישנם מספר אינטגרלים ידועים הממומשים לרוב כפונקציות מיוחדות, הנפוצה בהן היא ההתפלגות הנורמלית בעזרתה ניתן לחשב אינטגרלים אחרים.

דוגמה:

$\underset{-2x}{\overset{x}{\int }}{y}^2{e}^{-\frac{y^2}{2}}\mathrm{d}y=\sqrt{2\pi }(\mathrm{\Phi }(x)-\mathrm{\Phi }(-2x))-e^{-\frac{x^2}{2}}x+2e^{-2x^2}x$

דבר זה מאפשר לחשב את האינטגרל תוך שימוש בפונקציה שלה יש קירובים רציונליים מוכרים.

ראו גם

• שיטות אינטגרציה
• שיטות למציאת אינטגרלים לא מסוימים
• שיטות אנליטיות לחישוב אינטגרלים מסוימים

קישורים חיצוניים

• תבנית:MathWorld