התפלגות קושי

התפלגות קושי

פונקציית צפיפות ההסתברות
פונקציית ההסתברות המצטברת
מאפיינים
פרמטרים	${x}_0$ החציון, $\gamma$ סקלה
תומך	$\mathbb{ℝ}$
פונקציית צפיפות הסתברות (pdf)	$\frac{1}{\pi }\left[\frac{\gamma }{(x-x_0{)}^2+{\gamma }^2}\right]$
פונקציית ההסתברות המצטברת (cdf)	$\frac{1}{\pi }\arctan \left(\frac{x-x_0}{\gamma }\right)+\frac{1}{2}$
תוחלת	לא מוגדרת
סטיית תקן	לא מוגדרת
חציון	${x}_0$
ערך שכיח	${x}_0$
שונות	לא מוגדרת
אנטרופיה	$\ln (4\pi \gamma )$
פונקציה יוצרת מומנטים (mgf)	לא מוגדרת
פונקציה אופיינית	$\exp (x_{0}it-\gamma |t|)$
צידוד	לא מוגדר
גבנוניות	לא מוגדרת

התפלגות קוֹשִי (Cauchy), על שם המתמטיקאי הצרפתי אוגוסטן לואי קושי, היא התפלגות רציפה בעלת חשיבות במתמטיקה ובמספר תחומים בפיזיקה. בקרב פיזיקאים ההתפלגות מכונה לעיתים פילוג לורנץ (Lorentz), פילוג ברייט-ויגנר (Breit-Wigner) או לורנציאן.

הגדרה

התפלגות קושי מוגדרת כהתפלגות רציפה בעלת פונקציית צפיפות ההסתברות

$\begin{array}{rl}f(x;x_0,\gamma )& =\frac{1}{\pi \gamma [1+{\left(\frac{x-x_0}{\gamma }\right)}^2]}\\ & =\frac{1}{\pi }\left[\frac{\gamma }{(x-x_0{)}^2+{\gamma }^2}\right]\end{array}$

כאשר ${x}_0$ הוא פרמטר מיקום, אשר קובע את החציון של ההתפלגות, ואילו $\gamma$ הוא פרמטר סקלה, אשר קובע את רוחב ההתפלגות ובהתאם את גובה הערכים. בגבול שבו $\gamma \to 0$ נקבל את פונקציית הדלתא של דיראק.

המקרה הפרטי של התפלגות קושי עם פרמטרים $x_0=0$ ו $\gamma =1$ נקרא התפלגות קושי סטנדרטית, עם צפיפות התפלגות

$f(x;0,1)=\frac{1}{\pi (1+x^2)}$ .

תכונות

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

סכום של משתנים מקריים המתפלגים קושי

אם $X_1,X_2,\dots ,X_n$ הם משתנים מקריים בלתי-תלויים ושווי-התפלגות שנדגמו מהתפלגות קושי סטנדרטית, ממוצע המדגם שלהם $\overline{X}=\frac{1}{n}\sum _i{X}_i$ מתפלג קושי סטנדרטית. בפרט, הממוצע אינו מתכנס לתוחלת, ואכן התפלגות קושי אינה מקיימת את חוק המספרים הגדולים. ההוכחה של תכונה זו אפשרית על ידי אינטגרציה של פונקציית צפיפות ההסתברות או על ידי שימוש בפונקציה האופיינית של התפלגות קושי הסטנדרטית:$${\phi }_X(t)=\mathrm{E}\left[e^{iXt}\right]=e^{-|t|}$$עבור סכום הדגימות מקבלים ${\phi }_{\sum _i{X}_i}(t)=e^{-n|t|}$ כלומר $\overline{X}$ הוא משתנה מקרי בעל התפלגות קושי סטנדרטית.

באופן כללי יותר, אם $X_1,X_2,\dots ,X_n$ הם משתנים בלתי-תלויים בעלי התפלגות קושי עם פרמטרי מיקום $x_1,\dots ,x_n$ ופרמטרי סקלה ${\gamma }_1,\dots ,{\gamma }_n$ ו $a_1,\dots ,a_n$ הם מספרים ממשיים אזי $\sum _i{a}_i{X}_i$ מתפלג קושי עם פרמטר מיקום $\sum _i{a}_i{x}_i$ וסקאלה $\sum _i|a_i|{\gamma }_i$. כלומר חוק המספרים הגדולים אינו מתקיים עבור סכום משוקלל של משתני קושי בלתי תלויים.

פונקציה אופיינית

יהי $X$ משתנה מקרי המתפלג קושי. הפונקציה האופיינית של התפלגות קושי ניתנת על ידי

${\phi }_X(t)=\mathrm{E}\left[e^{iXt}\right]=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}f(x;x_0,\gamma )e^{ixt}dx=e^{i{x}_{0}t-\gamma |t|}$

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

$f(x;x_0,\gamma )=\frac{1}{2\pi }\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}{\phi }_X(t;x_0,\gamma )e^{-ixt}dt$ .

המומנט ה-n של ההתפלגות מתקבל מהצבת $t=0$ בנגזרת ה-n של הפונקציה האופיינית. יש לשים לב כי הפונקציה האופיינית אינה גזירה בראשית הצירים. ואכן המומנטים של התפלגות קושי אינם מוגדרים פרט למומנט האפס.

דיברגנץ קולבק-לייבלר

ניתן להביא את דיברגנץ קולבק-לייבלר בין שתי התפלגויות קושי עם פרמטרים $x_{0,1},{\gamma }_1$ ו $x_{0,2},{\gamma }_2$ כנוסחה סגורה סימטרית:

$\mathrm{K}\mathrm{L}(p_{x_{0,1},{\gamma }_1}:p_{x_{0,2},{\gamma }_2})=\log \frac{{({\gamma }_1+{\gamma }_2)}^2+{(x_{0,1}-x_{0,2})}^2}{4{\gamma }_1{\gamma }_2}$ .

אנטרופיה

האנטרופיה הדיפרנציאלית של התפלגות קושי נתונה על ידי

$\begin{array}{rl}H(\gamma )& =-\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}f(x;x_0,\gamma )\log (f(x;x_0,\gamma ))dx\\ & =\log (4\pi \gamma )\end{array}$

הנגזרת של פונקציית השברונים של התפלגות קושי היא

$Q^{\prime }(p;\gamma )=\gamma \pi {\sec }^2\left[\pi (p-\frac{1}{2})\right]$ .

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

$H(\gamma )=\underset{0}{\overset{1}{\int }}\log (Q^{\prime }(p;\gamma ))\mathrm{d}p=\log (4\pi \gamma )$ .

בניות

דגימה מהתפלגות קושי

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

באופן פורמלי, תהי נקודה $(x_0,\gamma )$ במישור x-y. נבחר קו העובר דרך הנקודה, כך שהזווית שלו ביחס לציר $x$ נבחרת באופן באופן אחיד (בין 90°- ל-90°+) באקראי. נקודת החיתוך של הישר עם ציר $x$ הוא מתפלגת קושי עם פרמטר מיקום $x_0$ ופרמטר סקאלה $\gamma$.

הגדרה זו נותנת דרך פשוטה לדגום מהתפלגות קושי סטנדרטית. תהי $u$ דגימה מהתפלגות אחידה רציפה מהקטע $[0,1]$. אזי ניתן ליצור דגימה $x$ מהתפלגות קושי סטנדרטית על ידי

$x=\tan (\pi (u-\frac{1}{2}))$ .

לחלופין, אם $U$ ו $V$ הם שני משתנים מקריים בלתי תלויים המתפלגים נורמלית עם תוחלת 0 ושונות 1 אז המנה $U/V$ מתפלגת התפלגות קושי סטנדרטית. ובאופן כללי, אם ל-$(U,V)$ סימטריה סיבובית במישור סביב ראשית הצירים, אזי המנה $U/V$ מתפלגת התפלגות קושי סטנדרטית.

פונקציית צפיפות התפלגות

פונקציית צפיפות ההסתברות של התפלגות קושי היא

$f(x;x_0,\gamma )=\frac{1}{\pi \gamma [1+{\left(\frac{x-x_0}{\gamma }\right)}^2]}=\frac{1}{\pi }\left[\frac{\gamma }{(x-x_0{)}^2+{\gamma }^2}\right]$

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

המקסימום של פונקציית הצפיפות של התפלגות קושי הוא ב-$x=x_0$ וערכו $\frac{1}{\pi \gamma }$.

פונקציית התפלגות מצטברת

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

$F(x;x_0,\gamma )=\frac{1}{\pi }\arctan \left(\frac{x-x_0}{\gamma }\right)+\frac{1}{2}$

ופונקציית השברונים (הפונקציה ההופכית לה) היא

$Q(p;x_0,\gamma )=x_0+\gamma \tan \left[\pi (p-\frac{1}{2})\right]$ .

הרבעון הראשון והשלישי הם $(x_0-\gamma ,x_0+\gamma )$, ומכאן $2\gamma$ הוא הטווח הבין-רבעוני של ההתפלגות.

הקשר להתפלגות t

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

שיערוך פרמטרים

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

$\overline{x}=\frac{1}{n}\sum _{i=1}^n{x}_i$

למרות שערכי הדגימות $x_i$ יתרכזו סביב $x_0$, ממוצע המדגם לא יתכנס כאשר מגדילים את הערך של n. זאת מפני שההסתברות לדגימה $x$ עם ערך מוחלט גדול, הולכת וגדלה. למעשה, התפלגות ממוצע המדגם תהיה שווה להתפלגות התצפיות עצמן; כלומר, ממוצע המדגם של מדגם גדול אינו טוב יותר (או גרוע יותר) משערוך הערך של $x_0$ מכל תצפית בודדת במדגם. באופן דומה, חישוב שונות המדגם יביא לערכים שיגדלו ככל שיגדל מספר התצפיות.

זו הסיבה שיש צורך בשיטות אחרות לשיערוך פרמטר המיקום $x_0$ ופרמטר קנה המידה $\gamma$. שיטה פשוטה אחת היא לקחת את הערך החציוני של המדגם כאומד של $x_0$, וחצי מהטווח הבין-רבעוני כאומד של $\gamma$. פותחו גם שיטות אחרות, מדויקות וחזקות יותר, שלעיתים נותנות תוצאות מדויקות בהרבה, למשל הממוצע המקוטע של 24% הדגימות המרכזיות מהמדגם, מייצר אומד של $x_0$ מדויק יותר מחציון המדגם.

ניתן להיעזר בשיטת הנראות המקסימלית כדי לאמוד את הפרמטרים $x_0$ ו $\gamma$. עם זאת, זה גישה זו מורכבת, מכיוון שהיא דורשת מציאת השורשים של פולינום מסדר גבוה, ושורשים מרובים ייצגו נקודות מקסימה מקומיות. כמו כן, בעוד שמעריך הנראות המקסימלית יעיל אסימפטוטית, הוא אינו יעיל יחסית עבור מדגמים קטנים. פונקציית הלוג-נראות עבור התפלגות קושי למדגם בגודל $n$ היא:

$\hat{\ell }(x_1,\dots ,x_n\mid x_0,\gamma )=-n\log (\gamma \pi )-\sum _{i=1}^n\log (1+{\left(\frac{x_i-x_0}{\gamma }\right)}^2)$

גזירה של פונקציית הנראות ביחס ל-$x_0$ ול-$\gamma$ והשוואת הנגזרות לאפס נותנות את המשוואות הבאות:

$\frac{d\ell }{d{x}_0}=\sum _{i=1}^n\frac{2(x_i-x_0)}{{\gamma }^2+{(x_i-x_0)}^2}=0$

$\frac{d\ell }{d\gamma }=\sum _{i=1}^n\frac{2{(x_i-x_0)}^2}{\gamma ({\gamma }^2+{(x_i-x_0)}^2)}-\frac{n}{\gamma }=0$

כאשר נשים לב כי במשוואת הנגזרת ביחס ל-$\gamma$ הביטוי

$\sum _{i=1}^n\frac{{(x_i-x_0)}^2}{{\gamma }^2+{(x_i-x_0)}^2}$

היא פונקציה מונוטונית של $\gamma$ וכי הפתרון עבור $\gamma$ חייב לקיים

$\min |x_i-x_0|\le \gamma \le \max |x_i-x_0|$ .

הפתרון עבור $x_0$ דורש מציאת השורשים של פולינום מסדר $2n-1$. הפתרון עבור $\gamma$ דורש מציאת השורשים של פולינום מסדר $2n$. לכן בדרך כלל יהיה צורך להיעזר בשיטות נומריות. היתרון של שיערוך נראות מקסימלית הוא יעילות אסימפטוטית גבוהה יותר, שמתרגמת לשונות נמוכה יותר של האומד; הערכת $x_0$ באמצעות חציון המדגם היא רק כ-81% יעילה מבחינה אסימפטוטית ביחס אומדן $x_0$ לפי הסבירות המקסימלית. ממוצע המדגם הקטוע המבוסס על 24% הדגימות המרכזיות יעיל אסימפטוטית בכ-88% יעיל ביחס למשערך הנראות המקסימלית של $x_0$.

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

מדיה וקבצים בנושא התפלגות קושי בוויקישיתוף   

• התפלגות קושי, באתר אנציקלופדיה למתמטיקה (באנגלית)
• .html התפלגות קושי, באתר MathWorld (באנגלית)   
• התפלגות קושי, באתר אנציקלופדיה בריטניקה (באנגלית)