לדלג לתוכן

אנליזה נומרית/אינטרפולציה: הפרשים סופיים

מתוך ויקיספר, אוסף הספרים והמדריכים החופשי

אינטרפולציה היא מציאת קירוב לערך הפונקציה בין שתי נקודות נתונות. מעתה והלאה נניח, כמו מקודם, כי המרחקים על ציר x בין כל שתי נקודות סמוכות הוא h. לכן כל נקודה ניתן לייצוג על ידי  xi+θh, כאשר  0≤θ≤1.

אינטרפולציה באמצעות הפרשים קידמיים

[עריכה]

נשתמש באופרטורים אשר פיתחנו קודם על מנת לייצג נקודות שאין לנו מידע עבורן. ניעזר בקשר  E=1+Δ ובמקרה הפרטי של בינום ניוטון:

 f(x+θh)=Eθf(x)=(1+Δ)θf(x)=[1+θΔ+θ(θ−1)2!Δ2+θ(θ−1)(θ−2)3!Δ3+...]f(x)

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

קירוב לינארי

[עריכה]
 f(x+θh)≈[1+θΔ]f(x)=f(x)+θ[f(x+h)−f(x)]

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

קירוב ריבועי

[עריכה]
 f(x+θh)≈[1+θΔ+θ(θ−1)2!Δ2]f(x)=f(x)+θ[f(x+h)−f(x)]+θ(θ−1)2![f(x+2h)−2f(x+h)+f(x)]

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

נשים לב כי עבור θ=0,1 נקבל חזרה את ערכי הפונקציה המקוריים. כללית: ע"י k נקודות ידועות נוכל לבנות שיטה מסדר k אשר תתן תוצאות מדויקות באותן נקודות ידועות, אך מכיוון שאנו מאלצים על פולינום מדרגה גבוהה לעבור דרך נקודות נתונות, הדיוק עלול להינזק.

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

 Df(xi)=f′(xi)=f(xi+1)−f(xi)h=Δhf(xi)⇒ f‴(xi)=(Δh)3f(xi)

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

 R=(θ(θ−1)(θ−2)3!Δ3)f(xi)≈h3θ(θ−1)(θ−2)3!f‴(c)

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

אינטרפולציה באמצעות הפרשים אחוריים

[עריכה]

ניעזר בקשר  E=(1+∇)−1 ובמקרה הפרטי של בינום ניוטון:

 f(x+θh)=Eθf(x)=(1−∇)−θf(x)=[1+θ∇+θ(θ+1)2!∇2+θ(θ+1)(θ+2)3!∇3+...]f(x)

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

קירוב לינארי

[עריכה]
 f(x+θh)≈[1+θ∇]f(x)=f(x)+θ[f(x)−f(x−h)]

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

קירוב ריבועי

[עריכה]

פרק זה לוקה בחסר. אתם מוזמנים לתרום לוויקיספר ולהשלים אותו. ראו פירוט בדף השיחה.



אינטרפולציה באמצעות הפרשים מרכזיים

[עריכה]

(להשלים)


הפרק הקודם:
אופרטורים של הפרשים סופיים
אינטרפולציה: הפרשים סופיים הפרק הבא:
אינטרפולציה: מינימום ריבועים