אנליזה נומרית/גזירה נומרית/הערכת השגיאה

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

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

[עריכה] הפרשים קדמיים

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

$\ f'(x_i)={1\over h}\Delta f(x_i)+ \left\{ -{1\over h}{\Delta^2f(x_i)\over 2} \right\}= {1\over h}(f(x_{i+1})-f(x_i)) + \left\{ -{\Delta\over 2h}(f(x_{i+1})-f(x_i)) \right\}$

נביט באיבר השגיאה, ונשתמש במשפט לגראנז':

$\ {\Delta\over 2h}(f(x_{i+1})-f(x_i))= {\Delta\over 2h}[hf'(c)]= {1\over 2}[f'(c+h)-f'(c)]={h\over 2}f''(\hat{c})$

$\ \Rightarrow f'(x_i)={1\over h}(f(x_{i+1})-f(x_i))- \left\{ {h\over 2}f''(c) \right\}$

אם כן, השגיאה פרופורציונית ל-h, או בכתיב מתמטי: $\ O(h)$ .

[עריכה] הפרשים מרכזיים

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

$\ f'(x_i)= {1\over h}\mu\delta f(x_i)+ \left\{ -{1\over 6h}\mu\delta^3f(x_i)+... \right\}= {1\over 2h}(f(x_{i+1})-f(x_{i-1}))+ \left\{ -{1\over 6h}\mu\delta^3f(x_i)+... \right\}$

נביט באיבר השגיאה ונשתמש במשפט לגראנז':

$\ {1\over 6h}\mu\delta^3f(x_i)= {1\over 6h}\mu\delta^2 \left[f(x_{i+{1\over 2}})-f(x_{i-{1\over 2}}) \right]= {1\over 6h}\mu\delta^2hf'(c)={1\over 6}\mu\delta \left[ f'(c+{h\over 2})-f'(c-{h\over 2}) \right]=$