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

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

הפרשים קדמיים

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

$\ f'(x_{i})={1 \over h}\Delta f(x_{i})+\left\{-{1 \over h}{\Delta ^{2}f(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 ^{3}f(x_{i})+...\right\}={1 \over 2h}(f(x_{i+1})-f(x_{i-1}))+\left\{-{1 \over 6h}\mu \delta ^{3}f(x_{i})+...\right\}$

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

$\ {1 \over 6h}\mu \delta ^{3}f(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 ^{2}hf'(c)={1 \over 6}\mu \delta \left[f'(c+{h \over 2})-f'(c-{h \over 2})\right]=$

$\ ={h \over 6}\mu \delta f''({\hat {c}})={h \over 6}{1 \over 2}[f''({\hat {c}}+h)-f''({\hat {c}}-h)]={h^{2} \over 6}f'''({\bar {c}})$

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

אם כן, השגיאה פרופורציונית ל-h², או בכתיב מתמטי: $\ O(h^{2})$. לכן אם $\ h\to 0$ אז השגיאה לפי הפרשים מרכזיים קטנה יותר.

נשים לב כי קודם הפעלנו את $\ \Delta ^{2}$ וקיבלנו $\ h^{2}f''$ ואילו כאן הפעלנו את $\ \delta ^{3}$ וקיבלנו $\ h^{3}f'''$.

הכללה

כללית: אם אבר השגיאה הוא $\ \Gamma ^{n}f(x_{i})$, כאשר גאמא הוא אופרטור הפרשים סופיים, אז השגיאה תהיה $\ h^{n}f^{(n)}(c)$.