מתאם עצמי

שמאל|ממוזער|400x400 פיקסלים| למעלה: סרטוט של סדרה של 100 מספרים אקראיים המסתירה פונקציית סינוס. למטה: פונקציית הסינוס מתגלה בקורלוגרם  שנוצר על ידי מתאם עצמי.
ממוזער|400x400 פיקסלים| השוואה חזותית של קונבולציה, מתאם צולב ו'''מתאם עצמי'''. עבור הפעולות הכוללות את הפונקציה , ובהנחה שהגובה של  הוא 1.0, ערך התוצאה ב-5 נקודות שונות מסומן על ידי השטח המוצלל מתחת לכל נקודה. כמו כן, הסימטריה של  היא הסיבה לכך ש־<math>g*f</math> ו־<math>f \star g</math> זהים בדוגמה זו.
**מתאם עצמי** או **אוטו־קורלציה** (נכתב גם **אוטוקורלציה**; ), המכונה גם **מתאם סדרתי** (serial correlation) במקרה של זמן בדיד, הוא פונקציה המודדת את המתאם בין אות לבין גרסה מושהית שלו. מתאם עצמי מעריך את הדמיון בין ערכי משתנה אקראי בנקודות זמן שונות. כלי מתמטי זה מסייע בזיהוי דפוסים חוזרים או מחזוריות סמויה בתוך אות שעלולים להיות מוסתרים על ידי רעש. מתאם עצמי משמש בהרחבה בעיבוד אותות, ניתוח סדרות עִתיות וניתוח במישור הזמן כדי לחקור את התנהגות הנתונים לאורך זמן. ובסטטיסטיקה מרחבית כדי לחקור קשרים בין דגימות או אירועים על פני מרחקים שונים.

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

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

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

יהי \{X_t\} תהליך סטוכסטי, כאשר t מייצגת נקודה כלשהי בזמן – ערך שלם עבור תהליך בזמן בדיד או מספר ממשי עבור תהליך בזמן רציף. הערך X_t הוא התוצאה (או המימוש ) של התהליך בזמן t.

בהנחה שלתהליך יש תוחלת \mu_t ושונות \sigma_t^2 בכל זמן t, אזי **פונקציית המתאם העצמי** בין הזמנים t_1 ו־t_2 מוגדרת כך:

כאשר \operatorname{E} הוא אופרטור התוחלת והקו העליון מייצג צמוד מרוכב. התוחלת עשויה שלא להיות מוגדרת היטב בכל המקרים.

אם מחסירים את התוחלת לפני הכפל, התוצאה היא **פונקציית השונות המשותפת העצמית** בין הזמנים t_1 ו־t_2.

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

### הגדרה לתהליך סטוכסטי סטציונרי במובן הרחב
תהליך סטציונרי במובן הרחב \left\{ X_t \right\} הוא תהליך שבו התוחלת \mu והשונות \sigma^2 אינן תלויות בזמן. בנוסף, פונקציית השונות המשותפת העצמית שלו תלויה רק בהפרש הזמנים (ההשהיה) בין שתי נקודות, t_1 ו־t_2, ולא במיקומן המוחלט בזמן. מכאן נובע כי פונקציות השונות המשותפת העצמית והמתאם העצמי ניתנות לייצוג כפונקציות של השהיית הזמן \tau=t_2-t_1, והן פונקציות זוגיות של \tau. כתוצאה מכך, **פונקציית המתאם העצמי** מקבלת את הצורה:

ו**פונקציית השונות המשותפת העצמית** היא:

ובפרט:

:\operatorname{K}_{XX}(0) = \sigma^2

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

מקדם המתאם העצמי של תהליך סטוכסטי מוגדר כך:

:\rho_{XX}(t_1,t_2) = \frac{\operatorname{K}_{XX}(t_1,t_2)}{\sigma_{t_1}\sigma_{t_2}} = \frac{\operatorname{E}\left[(X_{t_1} - \mu_{t_1}) \overline{(X_{t_2} - \mu_{t_2})} \right]}{\sigma_{t_1}\sigma_{t_2}}

כאשר פונקציה זו \rho_{XX} מוגדרת היטב, ערכיה נמצאים בטווח [-1,1], כך ש־1 מציין מתאם מושלם ו־1- מציין אנטי-מתאם מושלם.

עבור תהליך סטציונרי במובן הרחב (WSS), מקדם המתאם העצמי מקבל את הצורה הפשוטה:

:\rho_{XX}(\tau) = \frac{\operatorname{K}_{XX}(\tau)}{\sigma^2} = \frac{\operatorname{E} \left[(X_{t+\tau} - \mu)\overline{(X_{t} - \mu)}\right]}{\sigma^2}.

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

### תכונות
#### סימטריה
פונקציית המתאם העצמי \operatorname{R}_{XX} היא פונקציה זוגית, כלומר:

:\operatorname{R}_{XX}(t_1,t_2) = \overline{\operatorname{R}_{XX}(t_2,t_1)}

ובמקרה של תהליך WSS:

:\operatorname{R}_{XX}(\tau) = \overline{\operatorname{R}_{XX}(-\tau)}

#### מקסימום באפס
עבור תהליך WSS:

:\left|\operatorname{R}_{XX}(\tau)\right| \leq \operatorname{R}_{XX}(0)

בנוסף, \operatorname{R}_{XX}(0) הוא תמיד ממשי.

#### אי-שוויון קושי-שוורץ
לכל תהליך סטוכסטי מתקיים אי-שוויון קושי-שוורץ:

:\left|\operatorname{R}_{XX}(t_1,t_2)\right|^2 \leq \operatorname{E}\left[ |X_{t_1}|^2\right] \operatorname{E}\left[|X_{t_2}|^2\right]

#### מתאם עצמי של רעש לבן
עבור רעש לבן בזמן רציף, לפונקציית המתאם העצמי יהיה שיא חד ב־\tau=0 המיוצג על ידי פונקציית דלתא של דיראק, והיא תהיה אפס עבור כל ערך אחר של \tau.

#### משפט וינר-חינצ'ין
משפט וינר-חינצ'ין  קובע כי פונקציית המתאם העצמי \operatorname{R}_{XX} קשורה לצפיפות ההספק הספקטרית S_{XX} דרך התמרת פורייה:

:\operatorname{R}_{XX}(\tau) = \int_{-\infty}^\infty S_{XX}(f) e^{i 2 \pi f \tau} \, {\rm d}f

:S_{XX}(f) = \int_{-\infty}^\infty \operatorname{R}_{XX}(\tau) e^{- i 2 \pi f \tau} \, {\rm d}\tau 

במקרה של פונקציות ממשיות, \operatorname{R}_{XX} סימטרית ולכן ההתמרה מכילה רק רכיבי קוסינוס ממשיים (כלומר שאינם מוכפלים ביחידה המדומה):

:\operatorname{R}_{XX}(\tau) = \int_{-\infty}^\infty S_{XX}(f) \cos(2 \pi f \tau) \, {\rm d}f

:S_{XX}(f) = \int_{-\infty}^\infty \operatorname{R}_{XX}(\tau) \cos(2 \pi f \tau) \, {\rm d}\tau 

## יישומים
היכולת של מתאם עצמי לזהות תבניות חוזרות בנתונים הובילה ליישומים רבים בתחומים שונים:
* ספקטרוסקופיית מתאם פלואורסצנטי : משמשת לניתוח דיפוזיה מולקולרית ותגובות כימיות.
* מדידת ספקטרום אופטי ופולסי לייזר: אוטוקורלטורים אופטיים  מסייעים במדידת ספקטרומים אופטיים ופולסים אולטרה-קצרים.
* פיזור אור דינמי : קובע התפלגות של גודלי חלקיקים בתרחיפים ננומטריים באמצעות ניתוח תבניות נקודתיות של לייזר.
* עיבוד אותות ב-GPS: מסייע בתיקון השהיית התפשטות (propagation delay) על ידי יישור אות משוחזר של המקלט עם אות הלוויין.
* פיזור קרני רנטגן וניתוח ננו-מבנים: התמרת פורייה של פונקציית המתאם העצמי המרחבי משמש לחקר צפיפות אלקטרונים.
* מדעי השטח  ומיקרוסקופיה: קושר בין מורפולוגיית משטח לתכונות פונקציונליות במיקרוסקופיית סריקה.
* אופטיקה וניתוח קוהרנטיות: מתאמים עצמיים מנורמלים מודדים את הקוהרנטיות של שדות אלקטרומגנטיים.
* אסטרונומיה: קובע תדירויות פולסרים.
* מוזיקה ועיבוד אותות: משמש בזיהוי גובה צליל (לדוגמה, מכווני כלים, Auto-Tune) ובחישוב קצב מוזיקלי.
* עקיפת קרני X: פונקציית פטרסון מסייעת לשחזר מידע חסר של מופעי פורייה (Fourier phases) במבנים אטומיים.
* סטטיסטיקה וניתוח מרחבי: מתאם עצמי מרחבי מסייע בהערכת אי-ודאות במדגמים של אוכלוסיות הטרוגניות.
* ספקטרומטריית מסה: אלגוריתם SEQUEST מדרג ספקטרומים של פפטידים באמצעות מתאם עצמי.
* אסטרופיזיקה: משמש לניתוח התפלגות גלקסיות ולחקר מקורות X באנרגיה נמוכה.
* שרשראות מרקוב מונטה קרלו : מתאם עצמי נלקח בחשבון עבור הערכת שגיאות מדויקת.
* גאופיזיקה וניתוח סייסמי: מתאם עצמי מחושב לצורך ניתוח תכונות סייסמיות בסקרים תת-קרקעיים בתלת-ממד.
* אולטרסאונד רפואי: מתאם עצמי משמש להדמיית זרימת דם.
* אופטימיזציה של תיקי השקעות: מתאם עצמי בתשואות נכסים משפיע על החלטות השקעה.
* מדידת תדר במערכות חשמל: משמש בממסרי מדידה נומריים לאנליזה מדויקת של תדרי מערכת החשמל.

## ראו גם
*מורן I

## לקריאה נוספת
* 
* 
* 
* Solomon W. Golomb, and Guang Gong. Signal design for good correlation: for wireless communication, cryptography, and radar. Cambridge University Press, 2005.
* Klapetek, Petr (2018). *[https://www.elsevier.com/books/quantitative-data-processing-in-scanning-probe-microscopy/klapetek/978-0-12-813347-7 Quantitative Data Processing in Scanning Probe Microscopy: SPM Applications for Nanometrology]* (Second ed.). Elsevier. pp. 108–112 .
* Weisstein, Eric W. "Autocorrelation". MathWorld.

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

## הערות שוליים

קטגוריה:עיבוד אותות
קטגוריה:סטטיסטיקה