Loading [MathJax]/jax/element/mml/optable/BasicLatin.js
Skip to main content
Library homepage
 
Global

2.1: סדרות פורייה ואינטגרלים, פונקציית דיראק

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

δN(x)=(1/2π)(sin(N+12)x)/sin12x

כלומר,

fN(θ)=xxδN(θθ)f(θ)dθ

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

אנו ממשיכים לטרנספורמציה פורייה, שבה פונקציה בקו האינסופי באה לידי ביטוי כאינטגרל על פני רצף של סינוסים וקוסינוסים (או אקספוננציאלים שווים). eikx מתברר כי טיעונים מקבילים לאלה שהובילו δN(x) כעת לתת פונקציה δ(x) כזו

f(x)=δ(xx)f(x)dx

מול זה, אפשר בהחלט לתהות מה הטעם בפונקציה δ(x) שעל פיתול עם f(x) מחזירה את אותה פונקציה. f(x) הרלוונטיות של δ(x) תתברר בהמשך הקורס, כאשר מצבים של חלקיק קוונטי מיוצגים על ידי פונקציות גל בקו האינסופי, כמוf(x), והפעולות עליהם כוללות אופרטורים אינטגרליים הדומים לפיתול לעיל. עבודה עם פעולות על פונקציות אלה היא הכללת רצף של מטריצות הפועלות על וקטורים במרחב סופי-ממדי, והיא הייצוג האינסופי ממדי של δ(x) מטריצת היחידה. בדיוק כמו באלגברה של מטריקס המצבים העצמיים של מטריצת היחידה הם קבוצה של וקטורים המשתרעים על החלל, ורכיבי מטריצת היחידה קובעים את קבוצת תוצרי הנקודות של וקטורי בסיס אלה, פונקציית הדלתא קובעת את התוצר הפנימי הכללי של בסיס רצף של מצבים. זה ממלא תפקיד חיוני בפורמליזם הסטנדרטי למצבי רצף, ואתה צריך להכיר אותו!

סדרת פורייה

π<θπניתן להרחיב כל פונקציה אמיתית חלקה למדי f(θ) המוגדרת במרווח בסדרת פורייה,

f(θ)=A02+n=1(Ancosnθ+Bnsinnθ)

שם ניתן למצוא את המקדמים באמצעות מצב האורתוגונליות,

ππcosmθcosnθdθ=πδm,n

ואותו תנאי לתת sinnθ של:

An=1πππf(θ)cosnθdθ

Bn=1πππf(θ)sinnθdθ

שים לב כי עבור פונקציה זוגית, רק An הם nonzero, עבור פונקציה מוזרה רק Bn הם nonzero.

כמה חלק הוא "חלק באופן סביר"?

מספר המונחים של הסדרה הדרושים בכדי לתת קירוב טוב לפונקציה תלוי במהירות שבה הפונקציה משתנה. כדי לקבל מושג מה משתבש כאשר פונקציה אינה "חלקה", מאלף למצוא את סדרת הסינוס פורייה עבור פונקציית הצעד f(θ)=1forπ<θ0,f(θ)=1for0<θπ

באמצעות הביטוי Bn לעיל קל למצוא:

f(θ)=4π(sinθ+sin3θ3+sin5θ5+...)

לקיחת חצי תריסר הקדנציות הראשונות בסדרה נותנת:

התופעה של גיבס

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

כיצד הסכום על פני N תנאים קשור לפונקציה השלמה

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

אז החלפת ערכי המקדמים (משוואה\ ref {2.1.6} ו-\ ref {2.1.7})

An=1πππf(θ)cosnθdθ

Bn=1πππf(θ)sinnθdθ

בסדרה (משוואה\ ref {2.1.4})

fN(θ)=A02+Nn=1(Ancosnθ+Bnsinnθ)

נותן

fN(θ)=12πππf(θ)dθ+1πNn=1ππ(cosnθcosnθ+sinnθsinnθ)f(θ)dθ

=12πππf(θ)dθ+1πNn=1ππcosn(θθ)f(θ)dθ

כעת אנו יכולים להשתמש בזהות הטריגונומטרית

Nn=1cosnx=sin(N+12)x2sin12x12

כדי למצוא

fN(θ)=12πππsin(N+12)(θθ)sin12(θθ)f(θ)dθ=ππδN(θθ)f(θ)dθ

היכן

δN(x)=12πsin(N+12)xsin12x

(שים לב שהוכחת הזהות הטריגונומטרית היא פשוטה: כתוב z=eixcosnx=12(zn+zn), כך וסכם את ההתקדמות הגיאומטרית.)

הולכים אחורה לרגע וכותבים

δN(x)=1π(Nn=1cosnx+12)

קל לבדוק את זה ππδN(x)dx=1

כדי לעזור לדמייןδN(θ), הנהN=20:

image063.jpg

זה עתה קבענו כי השטח הכולל מתחת לעקומה = 1, וברור מהתרשים שכמעט כל השטח הזה נמצא מתחת לפסגה המרכזית, מכיוון שהאזורים המרוחקים מהמרכז חיוביים ושליליים כמעט באותה מידה. רוחב הפסגה המרכזית הואπ/(N+12), גובהו(N+12)/π.

תרגיל: בגדולN, בערך כמה רחוק למטה הוא טובל בתנודה הראשונה? (N/π2)

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

fN(θ)=ππδN(θθ)f(θ)dθ

ייתן fN(θ) קרובf(θ), ועבור פונקציות fN(θ) אלה נוטים f(θ) ככל ש - N יגדל.

ברור גם מדוע פיתול עקומה זו עם פונקציית צעד נותן חריגה ותנודות. נניח שהפונקציה f(θ) היא צעד, קופץ מ -0 ל -1 ב θ=0 מהצורה המתפתלת של האינטגרל, אתה אמור להיות מסוגל לשכנע את עצמך שהערך של fN(θ) בנקודה θ הוא השטח הכולל מתחת δN(θ) לעקומה משמאל לנקודה זו (אזור מתחת לאפס - כלומר מתחת לציר x - כמובן ספירה שלילית). בשביל θ=0 זה חייב להיות בדיוק 0.5 (שכן כל השטח תחת δN(θ) מוסיף 1). אבל אם אנחנו רוצים את הערך של fN(θ) at θ=π/(2N+1) (כלומר, הנקודה הראשונה מימין למקור שבו העקומה חותכת את ציר x), עלינו להוסיף את כל השטח משמאל לθ=π/(2N+1), אשר למעשה מסתכם בסך הכל שטח גדול מאחד, שכן שטח השאריות מימין לנקודה זו הוא שלילי בסך הכל. זה נותן את ההשתלטות.

סדרת פורייה במכניקת הקוונטים: אלקטרון בקופסה

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

תרגיל 2.1.1

להוכיח כי עבור פונקציהf(θ)=n=aneinθ, עם מורכבות כללית, an

12πππ|f(θ)|2dθ=n=|an|2

הרלוונטיות הפיזית של תוצאה זו היא כדלקמן: עבור אלקטרון מוגבל להיקף של טבעת של רדיוס יחידה, θ הוא המיקום של האלקטרון. בסיס אורתונורמלי של מצבי האלקטרון בטבעת זו הוא מכלול הפונקציות (1/2π)einθ עם מספר n שלם, סופרפוזיציה מנורמלת נכונה של מצבים אלה חייבת להיותn=|an|2=1, כך שההסתברות הכוללת למצוא את האלקטרון במצב כלשהו היא אחדות. אבל זה חייב גם אומר שההסתברות הכוללת למצוא את האלקטרון בכל מקום על הטבעת היא אחדות - וזה הצד השמאלי של המשוואה 2π הנ"ל - ביטול.

סדרת פורייה אקספוננציאלית

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

הצעד הראשון הוא טריוויאלי: עלינו להכליל מפונקציות אמיתיות לפונקציות מורכבות, כדי לכלול פונקציות גל בעלות זרם שאינו נעלם. ניתן לכתוב פונקציה מורכבת חלקה בסדרת פורייה פשוט על ידי מתן אפשרות An Bn להיות מורכבת, אך במקרה זה הרחבה טבעית יותר תהיה בסמכויות של eiθ,eiθ

אנו כותבים:

f(θ)=n=aneinθwithan=12πππeinθf(θ)dθ

וחזרה על השלבים לעיל fN(θ)=12πππNn=Nein(θθ)f(θ)dθ=12πππf(θ)dθ+1πππNn=1cosn(θθ)f(θ)dθ

בדיוק אותו ביטוי כמו קודם, ולכן נותן את אותו הדברδN(θ). זה לא מפתיע, כי השימוש N במונחים cosnθ=12(einθ+einθ),sinnθ=12i(einθeinθ) הראשונים An,Bn ניתן פשוט לארגן מחדש לסכום על פני einθ תנאים עבורNnN.

אלקטרון מחוץ לקופסה: טרנספורמציית פורייה

כדי לפרק חבילת גל לרכיבי גל המישור שלה, עלינו להרחיב את טווח האינטגרציה (π,π) מהמשמש לעיל ל(,). אנו עושים זאת על ידי שינוי גודל תחילה מ- (π,π) עד (L/2,L/2) ואז לקיחת הגבול. L

קנה המידה של המרווח מ- 2π עד L (בייצוג המורכב) נותן:

f(x)=n=ane2πinx/Lwherean=1LL/2L/2f(x)e2πinx/Ldx

הסכום n בלהיות מעל כל המספרים השלמים. זהו ביטוי עבור f(x) במונחים של גלי מישור eikx שבהם המותרים k נמצאים2πn/L, עם n=0,±1,±2,

אם נחזור על השלבים לעיל בגזירת הפונקציהδN(x), אנו מוצאים את הפונקציה המקבילה δLN(x)=1L(1+2Nn=1cos2πnxL)=sin((2N+1)πx/L)Lsin(πx/L)

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

לכן עלינו לקחת את הגבול N שהולך לאינסוף לפני שלוקחים את L לאינסוף.

זה מה שאנחנו עושים בשאר החלק הזה.

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

נניח שיש לנו חבילת גל כזו, נניח על אורךL1, לפיה אנו מתכוונים שהגל הוא בדיוק אפס מחוץ למתיחה של ציר האורךL1. מדוע לא פשוט לבטא זאת במונחים של סדרת N פורייה אינסופית המבוססת על מרווח גדול כלשהו (L/2,L/2) בתנאי שאורך חבילת הגל L1 נמצא לגמרי בתוך המרווח הזה? העניין הוא שניתוח כזה אכן ישחזר במדויק את חבילת הגלים בתוך המרווח, אך אותו סכום של גלי מישור שהוערך על פני כל ציר ה- x יגלה מחרוזת אינסופית של מנות גל זהות L זו מזו! זה לא מה שאנחנו רוצים.

כמקדים לנטילה L לאינסוף, הבה נכתוב את מונחי גל המישור האקספוננציאלי בסימון הסטנדרטיk,

e2πinx/L=eiknx

אז אנו מסכמים קבוצה (אינסופיתN) של גלי מישור בעלי ערכי מספר גלים

kn=2πn/L,n=0,±1,±2,,

קבוצה של מרווחים שווים עם k הפרדה. Δk=2π/L

שקול עכשיו מה קורה אם נכפיל את המרווח הבסיסי מ- (L/2,L/2) עד(L,L).

kהערכים החדשים המותרים הםkn=πn/L,n=0,±1,±2,, כך שההפרדה היא עכשיוΔk=π/L, חצי ממה שהיה קודם. ניכר שככל שאנו גדליםL, המרווח בין kn ערכים עוקבים הולך ופחות.

חוזרים למרווח האורךL, כותבים שיש Lan=a(kn),kn=2πn/L לנו

f(x)=n=ane2πinx/L=1Ln=a(kn)eiknx

נזכיר כי ניתן להגדיר את אינטגרל רימן על ידי

f(k)dk=lim

עםk_n=n\Delta k \; ,\; n = 0, \pm 1, \pm 2, ….

הביטוי בצד ימין של המשוואה עבור הוא f(x) בעל אותה צורה כמו הצד הימני של ההגדרה האינטגרלית של רימן, וכאן. \Delta k=2\pi/L

זאת אומרת, 2\pi f(x)=\dfrac{2\pi}{L}\sum_{n=-\infty}^{\infty} a(k_n) e^{ik_nx}=\sum_{n=-\infty}^{\infty} a(k_n) e^{ik_nx}\Delta k=\int\limits_{-\infty}^{\infty} a(k)e^{ikx}dk \label{2.1.26}

בגבול\Delta k\rightarrow 0, או באופן שווהL\rightarrow\infty. אנו כמובן מניחים כאן שהפונקציהa(k_n), שהגדרנו רק (עבור נתוןL) על קבוצת הנקודותk_n, נוטה לפונקציה a(k) רציפה בגבולL\rightarrow\infty.

מכאן נובע L שבגבול האינסופי, יש לנו את משוואות הטרנספורמציה של פורייה:

f(x)=\dfrac{1}{2\pi}\int\limits_{-\infty}^{\infty} a(k)e^{ikx}dk \;\; and \;\; a(k)=\lim_{L\rightarrow \infty} La_n=\lim_{L\rightarrow \infty}\int\limits_{-L/2}^{L/2} f(x)e^{-2\pi inx/L}dx=\int\limits_{-\infty}^{\infty} f(x)e^{-ikx}dx \label{2.1.27}

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

עכשיו לקחנו את שניהם N L ולאינסוף, מה קרה לתפקוד שלנו\delta_N^L(x)? זכור כי הנוהל שלנו למציאת f_N(\theta) במונחים של f(\theta) נתן את המשוואה

f_N(\theta)=\dfrac{1}{2\pi} \int\limits_{-\pi}^{\pi} f(\theta') d\theta' +\dfrac{1}{\pi} \int\limits_{-\pi}^{\pi} \sum_{n=1}^{N} \cos n(\theta-\theta')f(\theta')d\theta' \label{2.1.10}

ומתוך זה מצאנו\delta_N(\theta).

בעקבות אותו הליך פורמלי עם (L=\infty) טרנספורמציות פורייה, אנו נאלצים לקחת N אינסוף (לזכור שההליך N היה הגיוני רק אם נלקח לאינסוף לפני כןL), כך שבמקום משוואה עבור f_N(\theta) במונחים שלf(\theta), אנו מקבלים משוואה עבור f (x) במונחים של עצמו! בואו נכתוב את זה קודם ונחשוב אחר כך:

f(x)=\dfrac{1}{2\pi}\int\limits_{-\infty}^{\infty} \left( \int\limits_{-\infty}^{\infty} f(x')e^{-ikx'}dx' \right) e^{ikx}dk=\int\limits_{-\infty}^{\infty} \left( \int\limits_{-\infty}^{\infty} \dfrac{dk}{2\pi}e^{ik(x-x')} \right) f(x')dx' \label{2.1.28}

במילים אחרות, f(x)= \int\limits_{-\infty}^{\infty} \delta(x-x')f(x')dx' \label{2.1.29}

איפה \delta(x)=\int\limits_{-\infty}^{+\infty} \dfrac{dk}{2\pi}e^{ikx} \label{2.1.30}

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

מהטיעונים שלנו לעיל, עלינו להיות מסוגלים להתאושש \delta(x) כגבול של \delta_N^L(x) על ידי לקיחת תחילה N לאינסוף, לאחר מכןL. זאת אומרת,

\delta(x)=\lim_{L\rightarrow\infty}\left(\lim_{N\rightarrow\infty}\delta_N^L(x)\right) =\lim_{L\rightarrow\infty}\left( \lim_{N\rightarrow\infty}\dfrac{\sin((2N+1)\pi x/L)}{L\sin(\pi x/L)} \right) \label{2.1.31}

דרך להבין את הגבול הזה היא לכתוב M=(2N+1)\pi/L M ולשחרר לאינסוף לפני כןL. (פירוש הדבר שכאשר אנו לוקחים L בגדול בדרך לאינסוף, אנו לוקחים N הרבה יותר גדול!)

אז המונה הוא פשוט\sin Mx. בגבול האינסופיL, לכל סופי x המכנה צודק\pi x, שכן \sin\theta =\theta בגבול הקטן. \theta

מתוך זה,

\delta(x)=\lim_{M\rightarrow\infty} \dfrac{\sin Mx}{\pi x} \label{2.1.32}

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

כדי לראות איך זה קשור ל (גם לא מוגדר)\delta(x)=\int\limits_{-\infty}^{+\infty}(dk/2\pi)e^{ikx}, הזיכרון \delta_N^L(x) הגיע מהסדרה

\delta_N^L(x)=\dfrac{\sin(2N+1)\pi x/L}{L\sin\pi x/L}=\dfrac{1}{L}\left( 1+2\sum_{n=1}^N \cos\dfrac{2\pi nx}{L} \right) \label{2.1.33}

ביטוי הקוסינוס במונחים של אקספוננציאלים, ואז החלפת הסכום באינטגרל N בגבול הגדול, באותו אופן שעשינו קודם לכן, כתיבהk_n=2\pi n/L, כך שהמרווח בין עוקבים k_n הוא\Delta k=2\pi/L, כך: \int f(k)dk\ \cong (2\pi/L)\sum f(k_n)

\delta_N^L(x)=\dfrac{1}{L}\sum_{n=-N}^N e^{2\pi in/L} \cong \int\limits_{-2\pi N/L}^{2\pi N/L}\dfrac{dk}{2\pi}e^{ikx} =\dfrac{\sin (2\pi Nx/L)}{\pi x} \label{2.1.34}

אז ברור שאנחנו מגדירים את \delta(x) כגבול של האינטגרל\int\limits_{-2\pi N/L}^{2\pi N/L}(dk/2\pi)e^{ikx}, שמנותק בפתאומיות בערכים הגדולים\pm (2\pi N/L). למעשה, זה לא מאוד פיזי: תרחיש הרבה יותר מציאותי עבור חבילת גל אמיתית יהיה הפחתה הדרגתית בתרומות ממצבי תדר גבוה (או אורך גל קצר) - כלומר ניתוק עדין באינטגרל מעל k זה שימש להחלפת הסכום. n לדוגמה, הליך ניתוק סביר יהיה להכפיל את האינטגרנד ב-exp(-\Delta^2k^2), ואז לקחת את הגבול של קטן. \Delta

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

\delta(x)=\lim_{\Delta\rightarrow 0}\int\limits_{-\infty}^{\infty}\dfrac{dk}{2\pi}e^{ikx}e^{-\Delta^2 k^2}=\lim_{\Delta\rightarrow 0} \dfrac{1}{(4\pi\Delta^2)^{1/2}}e^{-x^2/4\Delta^2}=\lim_{\Delta\rightarrow 0} \delta_{\Delta}(x) \label{2.1.35}

כלומר, ניתן להגדיר את פונקציית הדלתא כ"גבול הצר "של חבילת גל גאוסית בשטח כולל 1. בניגוד לפונקציה\delta_N(\theta), אין \delta_{\Delta}(x) לה קווי צד מתנדנדים, הודות להחלקה שלנו מחוץ לניתוק הרווח העליון של k, כך שהפסקות צעד אינן יוצרות את תופעת הצילום של גיבס - במקום זאת, צעד יוחלק על פני מרחק של סדר. \Delta

מאפייני פונקציית הדלתא

פשוט לאמת את המאפיינים הבאים מההגדרה כמגבלה של חבילת גל גאוסית:

\int \delta(x)dx=1, \; \delta(x)=0 \; for \; x\neq 0 \label{2.1.36}

\delta(x)=\delta(-x), \; \delta(ax)=\dfrac{1}{|a|}\delta(x) \label{2.1.37}

\int \delta(a-x)\delta(x-b)dx=\delta(a-b) \label{2.1.38}

הגדרה נוספת, וקשר עם אינטגרל הערך העיקרי

אין דרך ייחודית להגדיר את פונקציית הדלתא, ונהלי ניתוק אחרים יכולים לתת תובנות שימושיות. לדוגמה, ניתן לפצל את אינטגרל ה- k -space לשניים ולחתוך אקספוננציאלי פשוט על שני החצאים, כלומר נוכל לקחת את ההגדרה \delta(x)= \lim_{\varepsilon\rightarrow 0}\left( \int\limits_{-\infty}^{0}\dfrac{dk}{2\pi}e^{ikx}e^{\varepsilon k}+\int\limits_{0}^{\infty}\dfrac{dk}{2\pi}e^{ikx}e^{-\varepsilon k} \right) \label{2.1.39}

הערכת האינטגרלים,

\delta(x)=\lim_{\varepsilon\rightarrow 0}\dfrac{1}{2\pi} \left(\dfrac{1}{ix+\varepsilon}-\dfrac{1}{ix-\varepsilon}\right)=\lim_{\varepsilon\rightarrow 0}\dfrac{1}{\pi}\left(\dfrac{\varepsilon}{x^2+\varepsilon^2}\right) \label{2.1.40}

קל לבדוק שפונקציה זו מנורמלת כראוי על ידי ביצוע שינוי המשתנה x=\varepsilon\tan\theta ושילוב מ- -\pi/2 to\pi/2. ייצוג זה של פונקציית הדלתא יתגלה כשימושי בהמשך. שים לב שנחשב כפונקציה של משתנה מורכב, לפונקציית הדלתא יש שני קטבים על הציר הדמיוני הטהור בz=\pm i\varepsilon.

ההגדרה הסטנדרטית של אינטגרל הערך העיקרי היא:

\int\limits_{-D}^{D}f(x)\dfrac{P}{x}dx=\lim_{z\rightarrow 0}\left(\int\limits_{-D}^{-z}\dfrac{f(x)}{x}dx+\int\limits_{z}^{D}\dfrac{f(x)}{x}dx\right) \label{2.1.41}

לא קשה לראות שעבור פונקציה ניתנת להבחנה רציפה f(x) זה שווה ערך ל

\int\limits_{-D}^{D}f(x)\dfrac{P}{x}dx=\lim_{z\rightarrow 0}\int\limits_{-D}^{D}f(x)\dfrac{x}{x^2+\varepsilon^2}dx \label{2.1.42}

לכן ניתן לכתוב את מפעיל הערך העיקרי באופן סמלי:

\dfrac{P}{x}=\lim_{z\rightarrow 0}\dfrac{x}{x^2+\varepsilon^2}=\lim_{z\rightarrow 0}\dfrac{1}{2}\left(\dfrac{1}{x+i\varepsilon}+\dfrac{1}{x-i\varepsilon}\right) \label{2.1.43}

אם מחברים את זה יחד עם הייצוג הדומה של פונקציית הדלתא לעיל, ולוקחים את הגבול להבנה, יש לנו את התוצאה השימושית: \varepsilon\rightarrow 0

\dfrac{1}{x\pm i\varepsilon}=\dfrac{P}{x}\mp i\pi\delta(x) \label{2.1.44}

תרגילים

1. להוכיח את משפט פרסבל: If \;\; f(x)=\int_{-\infty}^{\infty} \dfrac{dk}{2\pi}a(k)e^{ikx}, \;\; then\;\; \int_{-\infty}^{\infty}|f(x)|^2dx=\int_{-\infty}^{\infty} \dfrac{dk}{2\pi}|a(k)|^2. \label{2.1.45}

2. הוכיח את הכלל עבור טרנספורמציית פורייה של פיתול של שתי פונקציות:

If \;\; f(x)=\int_{-\infty}^{\infty} \dfrac{dk}{2\pi}a(k)e^{ikx},\;\; g(x)=\int_{-\infty}^{\infty} \dfrac{dk'}{2\pi}b(k')e^{ik'x},\;\; then\;\; \int_{-\infty}^{\infty} f(x-x')g(x')dx'=\int_{-\infty}^{\infty} \dfrac{dk}{2\pi}a(k)b(k)e^{ikx} \label{2.1.46}