אינטגרל גאוסיאני

אינטגרל גאוסיאני הוא אינטגרל מסוים על פונקציית צפיפות של התפלגות נורמלית, כלומר:
<math display="block">\ \int \frac{1}{\sqrt{2\pi \sigma^2}} e^{-\frac{ (x - \mu)^2 }{2 \sigma^2}} \,\mathrm{d}x</math>
והכללותיו. באופן עקרוני, החלק הלא טריוויאלי באינטגרל זה הוא האינטגרל:
<math display="block">\ \int e^{-x^2} \,\mathrm{d}x</math>
וברגע שיודעים לפתור אותו קל לפתור גם את האינטגרל הכללי יותר שמופיע למעלה.

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

== גאוסיאן במשתנה אחד ==
הנוסחה הכללית עבור גאוסיאן במשתנה אחד היא כדלהלן:
<math display="block">\ \int_{-\infty}^{+\infty} e^{-ax^2 - bx - c}\, \mathrm{d}x = \sqrt{\frac{\pi}{a}} e^{\frac{b^2 - 4ac}{4a}} = \sqrt{\frac{\pi}{a}} \exp{\left(\frac{b^2 - 4ac}{4a}\right)}, \quad a > 0 </math>

=== הוכחת הנוסחה ===
את הנוסחה מוכיחים בשלושה שלבים:
# מחשבים את <math>\int_{-\infty}^{+\infty} e^{-x^2}\,\mathrm{d}x</math>
# מחשבים את <math>\int_{-\infty}^{+\infty} e^{-ax^2}\,\mathrm{d}x</math> באמצעות החלפת משתנים.
# מחשבים את <math>\int_{-\infty}^{+\infty} e^{-ax^2 - bx - c}\, \mathrm{d}x</math> באמצעות השלמה לריבוע.

==== שלב 1 ====
קודם נראה שהאינטגרל מתכנס:
<math display="block">\begin{align}
0 < \int_{-\infty}^{\infty} e^{-x^2} \mathrm{d}x & \leq \int_{-1}^{1} e^{-x^2}\,\mathrm{d}x + \int_{-\infty}^{-1} -e^{-x^2}\,\mathrm{d}x + \int_{1}^{+\infty} e^{-x^2}\,\mathrm{d}x 
\\
& \le \int_{-1}^{1} e^{-x^2}\,\mathrm{d}x + \int_{-\infty}^{-1} -xe^{-x^2}\,\mathrm{d}x + \int_{1}^{+\infty} xe^{-x^2}\,\mathrm{d}x < \infty
\end{align}</math>
שכן עבור <math>| x | > 1</math> מתקיים <math>e^{-x^2} < |x| e^{-x^2}</math> והאינטגרלים <math>\int_1^{\infty} x e^{-x^2} \mathrm{d}x = \left[ -\frac{e^{-x^2}}{2} \right]_1^{\infty} = \frac{e^{-1}}{2} = \frac{1}{2e}</math> והמקביל לו מתכנסים.

שנית, נציב <math>I = \int_{-\infty}^{+\infty} e^{-x^2}\,\mathrm{d}x</math>, נכפול את <math>\ I</math> באותו אינטגרל:
<math> I^2 = \int_{-\infty}^{+\infty} e^{-x^2}\,\mathrm{d}x \int_{-\infty}^{+\infty} e^{-y^2}\,\mathrm{d}y </math>.

מחוקי אינטגרל כפול (משפט פוביני) וכפל אקספוננטים מקבלים ש
<math>\ I^2 = \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} e^{-(x^2 + y^2)}\,\mathrm{d}x \, \mathrm{d}y</math>.
נשים לב שזוהי אינטגרציה על כל המישור <math>\ x-y</math>. נבצע החלפת משתנים לתיאור האינטגרציה בקואורדינטות קוטביות <math>\ r, \phi</math>
כאשר <math>\ r^2 = x^2 + y^2</math>
ו-<math>\ \phi</math> היא הזווית בין <math>\ r</math> לציר ה-<math>\ x</math>.

את אלמנט הנפח האינפיניטסימלי מחשבים באמצעות היעקוביאן של הטרנספורמציה <math>\ \mathrm{d}x \, \mathrm{d}y = r \mathrm{d}r \, \mathrm{d}\phi</math>, ומקבלים ש-<math>\ I^2 = \int_{0}^{2 \pi} \int_{0}^{\infty} e^{-r^2} r \,\mathrm{d}r \,\mathrm{d}\phi</math>.

את החלק הזוויתי קל לחשב, שכן האינטגרנד לא תלוי בזווית. מקבלים ש-<math display="block">\ I^2 = \int_{0}^{\infty} e^{-r^2} r \,\mathrm{d}r \int_{0}^{2 \pi}\, \mathrm{d}\phi =
2 \pi \int_{0}^{\infty} e^{-r^2} r \,\mathrm{d}r = \pi \int_{0}^{\infty}e^{-r^2} 2r \,\mathrm{d}r </math>.
כעת נשים לב ש-<math>\ (r^2)' = 2r</math>, ולכן <math>\ \frac{\mathrm{d}}{\mathrm{d}r} e^{-r^2} = -e^{-r^2} 2r </math>, ולכן <math>\ I^2 = \pi \left[ -e^{-r^2} \right]_{0}^{\infty} = \pi ( 0 - (-1)) = \pi</math>.

לוקחים שורש ריבועי חיובי (כי האינטגרנד חיובי), ומקבלים <math>\ I = \sqrt{\pi}</math>
או לסיכום:
<math display="block">I = \int_{-\infty}^{+\infty} e^{-x^2} \,\mathrm{d} x = \sqrt{\pi}</math>

==== שלב 2 ====
מבצעים את החלפת המשתנים הבאה: <math>\ y^2 = a x^2 \ , \ y = \sqrt{a} x \ , \ \mathrm{d}x = \frac{\mathrm{d}y}{\sqrt{a}}</math> 

ואז: <math>\int_{-\infty}^{+\infty} e^{-ax^2}\,\mathrm{d}x = \frac{1}{\sqrt{a}} \int_{-\infty}^{\infty} e^{-y^2}\, \mathrm{d}y = \sqrt{\frac{\pi}{a}} </math>

==== שלב 3 ====
ההשלמה לריבוע:
: <math>\ -a x^2 - b x - c = -a \left( x + \frac{b}{2a} \right)^2 + \frac{b^2}{4a} - c</math>
האינטגרציה על הריבוע נותנת <math>\ \sqrt{\pi / a}</math> ואילו הגורם הקבוע באקספוננט כופל אותו.

=== הערות נוספות ===
יש לציין שהפונקציה עליה מחושב האינטגרל היא פונקציה זוגית ולכן
:<math>\ \int_{-\infty}^{\infty} e^{-x^2}\, \mathrm{d}x = 2 \int_{0}^{\infty} e^{-x^2}\, \mathrm{d}x </math>

ניתן לחשב את האינטגרל גם כאשר האינטגרנד נכפל בחזקה של <math>\ x</math>, 
:<math>\begin{align}
&\int_0^\infty x^{2n} e^{-x^2/a^2}\,\mathrm{d}x
 &&= \sqrt{\pi} \frac{(2n)!}{n!} \left(\frac{a}{2}\right)^{2n+1} \\
&\int_0^\infty x^{2n+1}e^{-x^2/a^2}\,\mathrm{d}x
 &&= \frac{n!}{2} a^{2n+2}
\end{align}</math>

== אינטגרל גאוסיאני במספר משתנים ==
אינטגרל של n משתנים בתבנית ביליניארית:
:<math>
\int \exp\left( - \frac 1 2 A_{ij} x^i x^j \right) \,\mathrm{d}^nx
=
\sqrt{\frac{(2\pi)^n}{\det A}}
</math>
כאשר A היא מטריצה סימטרית חיובית לחלוטין.

== אינטגרל גאוסיאני עם מקדם מדומה ==
כאשר המקדם של <math>x^2</math> הוא <math>\pm i = \sqrt{-1}</math> עדיין אפשר לחשב את האינטגרל. השיטה הנכונה היא לבצע זאת באמצעות מסילה במישור המרוכב. מקבלים:
*<math>\int_{-\infty}^{\infty} e^{-ix^2} \,\mathrm{d}x = (1-i) \sqrt{ \frac{\pi}{2} }</math>
*<math>\int_{-\infty}^{\infty} e^{+ix^2} \,\mathrm{d}x = (1+i) \sqrt{ \frac{\pi}{2} }</math>

== קישורים חיצוניים ==
* סרטון עם פתרון של [https://www.youtube.com/watch?v=LXL_djg0Q2Y אינטגרל גאוסיאני] בעברית מאת אלון סלע
* 
* 

גאוסיאני
קטגוריה:קרל פרידריך גאוס