גרעין דיריכלה

באנליזה מתמטית, גרעין דיריכלה (באנגלית: Dirichlet kernel), הנקרא על שם המתמטיקאי הגרמני יוהאן פטר גוסטב לז'ן דיריכלה, הוא שמן המשותף של קבוצת פונקציות המוגדרות באופן הבא: D n ( x ) = 1 2 π ∑ k = − n n e i k x = 1 2 π ( 1 + 2 ∑ k = 1 n cos ⁡ ( k x ) ) = sin ⁡ ( ( n + 1 / 2 ) x ) 2 π sin ⁡ ( x / 2 ) , {\displaystyle D_{n}(x)={\frac {1}{2\pi }}\sum _{k=-n}^{n}e^{ikx}={\frac {1}{2\pi }}\left(1+2\sum _{k=1}^{n}\cos(kx)\right)={\frac {\sin \left(\left(n+1/2\right)x\right)}{2\pi \sin(x/2)}},} כאשר n מספר שלם אי-שלילי. פונקציות הגרעין הן מחזוריות עם מחזור 2 π {\displaystyle 2\pi } . חשיבותו של גרעין דיריכלה נובעת מהקשר שלו לטורי פורייה. הקונבולוציה של גרעין דיריכלה ה-n עם כל פונקציה f בעלת מחזור 2 π {\displaystyle 2\pi } היא קירוב פורייה ה-n של f, דהיינו: ( D n ∗ f ) ( x ) = ∫ − π π f ( y ) D n ( x − y ) d y = ∑ k = − n n f ^ ( k ) e i k x , {\displaystyle (D_{n}*f)(x)=\int _{-\pi }^{\pi }f(y)D_{n}(x-y)\,dy=\sum _{k=-n}^{n}{\hat {f}}(k)e^{ikx},} כאשר f ^ ( k ) = 1 2 π ∫ − π π f ( x ) e − i k x d x {\displaystyle {\widehat {f}}(k)={\frac {1}{2\pi }}\int _{-\pi }^{\pi }f(x)e^{-ikx}\,dx} הוא מקדם פורייה ה-k של f. מעובדות אלו עולה כי כדי לחקור התכנסות של טורי פורייה מספיק לחקור תכונות של גרעין דיריכלה. דיריכלה הציג את הסכומים הטריגונומטריים הללו במאמר מ-1829 שעסק בביסוס התאוריה העוסקת בהתכנסותם של טורי פורייה, ובמסגרתו הוא עשה שימוש בגרעין זה כדי להוכיח את משפט פורייה על התכנסות נקודתית של טורי פורייה. הניסוח המודרני של משפט פורייה הוא כדלקמן: "תהי f פונקציה רציפה למקוטעין מחזורית בעלת מחזור 2 π {\displaystyle 2\pi } ותהי x 0 {\displaystyle x_{0}} נקודה שבה קיימות הנגזרות החד-צדדיות של f. אז טור פורייה של f בנקודה x 0 {\displaystyle x_{0}} מתכנס ל-: f ( x 0 + ) + f ( x 0 − ) 2 {\displaystyle {\frac {f(x_{0}^{+})+f(x_{0}^{-})}{2}}} ." כאשר f ( x 0 + ) , f ( x 0 − ) {\displaystyle f(x_{0}^{+}),f(x_{0}^{-})} הם הגבולות החד-צדדיים של f מימין ומשמאל לנקודה x 0 {\displaystyle x_{0}} . נסמן ב- S n ( f ; x 0 ) {\displaystyle S_{n}(f;x_{0})} את קירוב פורייה מסדר n של f בנקודה x 0 {\displaystyle x_{0}} ; נראה ש-: lim n → ∞ ( S n ( f ; x 0 ) − f ( x 0 + ) + f ( x 0 − ) 2 ) = 0 {\displaystyle \lim _{n\to \infty }(S_{n}(f;x_{0})-{\frac {f(x_{0}^{+})+f(x_{0}^{-})}{2}})=0} לשם כך נבסס תחילה שתי תכונות חשובות של גרעין דיריכלה D n ( u ) {\displaystyle D_{n}(u)} : D n ( u ) {\displaystyle D_{n}(u)} היא פונקציה זוגית - נובע מיידית מהזוגיות של הקוסינוסים. ∫ − π π D n ( u ) d u = 2 ∫ 0 π D n ( u ) d u = 1 {\displaystyle \int _{-\pi }^{\pi }D_{n}(u)du=2\int _{0}^{\pi }D_{n}(u)du=1} - נובע מכך שרק האיבר הקבוע 1 2 π {\displaystyle {\frac {1}{2\pi }}} תורם לערך האינטגרל, שכן הממוצע של כל אחד מהקוסינוסים בגרעין דיריכלה על פני הקטע [ − π , π ] {\displaystyle [-\pi ,\pi ]} הוא אפס. לפיכך, אם ניעזר בעובדה שהקונבולוציה של (x)f עם גרעין דיריכלה ה-n נותנת את S n ( f ; x 0 ) {\displaystyle S_{n}(f;x_{0})} , נוכל לרשום מחדש את הגבול אותו ניסינו לחשב כ-: S n ( f ; x 0 ) − f ( x 0 + ) + f ( x 0 − ) 2 = ∫ − π 0 ( f ( x 0 + u ) − f ( x 0 − ) ) D n ( u ) d u + ∫ 0 π ( f ( x 0 + u ) − f ( x 0 + ) ) D n ( u ) d u {\displaystyle S_{n}(f;x_{0})-{\frac {f(x_{0}^{+})+f(x_{0}^{-})}{2}}=\int _{-\pi }^{0}(f(x_{0}+u)-f(x_{0}^{-}))D_{n}(u)du+\int _{0}^{\pi }(f(x_{0}+u)-f(x_{0}^{+}))D_{n}(u)du} כעת, נוכיח שכל אחד מהאינטגרלים המחוברים שואף לאפס כאשר n שואף לאינסוף. תחילה נרשום את המחובר השני (ההוכחה עבור המחובר הראשון זהה) בעזרת הצורה הסופית של גרעין דיריכלה: ∫ 0 π f ( x 0 + u ) − f ( x 0 + ) 2 π s i n ( u / 2 ) s i n ( ( n + 1 2 ) u ) d u {\displaystyle \int _{0}^{\pi }{\frac {f(x_{0}+u)-f(x_{0}^{+})}{2\pi \mathbb {sin} (u/2)}}\mathbb {sin} ((n+{\frac {1}{2}})u)du} על פי למת רימן-לבג, לכל פונקציה רציפה למקוטעין (g(u מקדמי פורייה שלה שואפים לאפס כאשר n שואף לאינסוף, ולכן מספיק שנוכיח כי g ( u ) = f ( x 0 + u ) − f ( x 0 + ) 2 s i n ( u / 2 ) {\displaystyle g(u)={\frac {f(x_{0}+u)-f(x_{0}^{+})}{2\mathbb {sin} (u/2)}}} רציפה למקוטעין בקטע [ − π , π ] {\displaystyle [-\pi ,\pi ]} . ברור שלכל u ≠ 0 {\displaystyle u\neq 0} המכנה של (g(u שונה מאפס ולכן הגבולות החד-צדדיים שלה קיימים (תחת ההנחה של המשפט ש-f סופית בכל מקום). הנקודה "הבעייתית" היחידה היא x 0 {\displaystyle x_{0}} עצמה, כלומר כאשר u שווה לאפס. לשם חישוב הגבול של g כאשר u שואף לאפס, נסתייע בגבול הידוע lim u → 0 s i n ( u ) u = 1 {\displaystyle \lim _{u\to 0}{\frac {\mathbb {sin} (u)}{u}}=1} ונקבל: lim u → 0 + f ( x 0 + u ) − f ( x 0 ) 2 s i n ( u / 2 ) = ( lim u → 0 + f ( x 0 + u ) − f ( x 0 ) u ) ⋅ ( lim u → 0 + u / 2 s i n ( u / 2 ) ) = f x 0 + ′ ⋅ 1 = f x 0 + ′ {\displaystyle \lim _{u\to 0^{+}}{\frac {f(x_{0}+u)-f(x_{0})}{2\mathbb {sin} (u/2)}}=(\lim _{u\to 0^{+}}{\frac {f(x_{0}+u)-f(x_{0})}{u}})\cdot (\lim _{u\to 0^{+}}{\frac {u/2}{\mathbb {sin} (u/2)}})=f'_{x_{0}^{+}}\cdot 1=f'_{x_{0}^{+}}} על פי הנחות המשפט הנגזרת החד-צדדית קיימת, ולכן הוכחנו את מה שנדרש. בעלת חשיבות רבה היא העובדה שנורמת L1 של Dn בקטע [ 0 , 2 π ] {\displaystyle [0,2\pi ]} מתבדרת לאינסוף כאשר n שואף לאינסוף. ניתן להראות (ראו למטה) ש-: ‖ D n ‖ L 1 = Ω ( log ⁡ n ) . {\displaystyle \|D_{n}\|_{L^{1}}=\Omega (\log n).} באמצעות טיעון מדויק יותר העושה שימוש בסכומי רימן כדי להעריך את תרומת הסביבה הגדולה ביותר של אפס שבה D n {\displaystyle D_{n}} חיובי, ושימוש באי שוויון ינסן כדי להעריך את התרומה של שאר הקטע לנורמה, ניתן להראות גם ש-: ‖ D n ‖ L 1 ≥ 4 Si ⁡ ( π ) + 8 π log ⁡ n . {\displaystyle \|D_{n}\|_{L^{1}}\geq 4\operatorname {Si} (\pi )+{\frac {8}{\pi }}\log n.} כאשר Si {\displaystyle \operatorname {Si} } היא הפונקציה הקדומה של sinc. היעדרות זו של אינטגרביליות אחידה עומדת מאחורי תופעות מתמטיות רבות הקשורות להתבדרותם של טורי פורייה. הוכחה מדויקת של התוצאה הראשונה שהוזכרה, ש- ‖ D n ‖ L 1 [ 0 , 2 π ] = Ω ( log ⁡ n ) {\displaystyle \|D_{n}\|_{L^{1}[0,2\pi ]}=\Omega (\log n)} , היא: ∫ 0 2 π | D n ( x ) | d x ≥ ∫ 0 π | sin ⁡ [ ( 2 n + 1 ) x ] | x d x ≥ ∑ k = 0 2 n ∫ k π ( k + 1 ) π | sin ⁡ ( s ) | s d s ≥ | ∑ k = 0 2 n ∫ 0 π sin ⁡ ( s ) ( k + 1 ) π d s | = 2 π H 2 n + 1 ≥ 2 π log ⁡ ( 2 n + 1 ) , {\displaystyle {\begin{aligned}\int _{0}^{2\pi }|D_{n}(x)|\,dx&\geq \int _{0}^{\pi }{\frac {|\sin[(2n+1)x]|}{x}}\,dx\\[5pt]&\geq \sum _{k=0}^{2n}\int _{k\pi }^{(k+1)\pi }{\frac {|\sin(s)|}{s}}\,ds\\[5pt]&\geq \left|\sum _{k=0}^{2n}\int _{0}^{\pi }{\frac {\sin(s)}{(k+1)\pi }}\,ds\right|\\[5pt]&={\frac {2}{\pi }}H_{2n+1}\\[5pt]&\geq {\frac {2}{\pi }}\log(2n+1),\end{aligned}}} כאשר נעזרנו בגבול ידוע כדי להניח ש- 2 / x ≤ 1 / | sin ⁡ ( x / 2 ) | {\displaystyle 2/x\leq 1/|\sin(x/2)|} , ואילו H n {\displaystyle H_{n}} הוא המספר ההרמוני ה-n. ניקח את מסרק דיראק, שאמנם אינו פונקציה במובן המקובל אלא סוג של "פונקציה מוכללת" המכונה גם "התפלגות", ונכפיל אותו ב- 2 π {\displaystyle 2\pi } . כך מקבלים את איבר הזהות עבור קונבולוציות עם פונקציות בעלות מחזור 2 π {\displaystyle 2\pi } . במילים אחרות: f ∗ ( 2 π δ ) = f {\displaystyle f*(2\pi \delta )=f} בעבור כל פונקציה f בעלת מחזור 2 π {\displaystyle 2\pi } . הצגת פורייה של פונקציית דיראק מחזורית זאת היא: 2 π δ ( x ) ∼ ∑ k = − ∞ ∞ e i k x = ( 1 + 2 ∑ k = 1 ∞ cos ⁡ ( k x ) ) . {\displaystyle 2\pi \delta (x)\sim \sum _{k=-\infty }^{\infty }e^{ikx}=\left(1+2\sum _{k=1}^{\infty }\cos(kx)\right).} הזהות הטריגונומטרית ∑ k = − n n e i k x = sin ⁡ ( ( n + 1 / 2 ) x ) sin ⁡ ( x / 2 ) {\displaystyle \sum _{k=-n}^{n}e^{ikx}={\frac {\sin((n+1/2)x)}{\sin(x/2)}}} ניתנת לקבלה בקלות אם עושים שימוש בהצגה המרוכבת של טור הקוסינוסים, זאת שכן ההצגה המרוכבת היא למעשה טור הנדסי סופי עם מנה e i x {\displaystyle e^{ix}} . לפיכך, מכיוון שהסיכום של טור הנדסי סופי הוא ∑ k = 0 n a r k = a 1 − r n + 1 1 − r . {\displaystyle \sum _{k=0}^{n}ar^{k}=a{\frac {1-r^{n+1}}{1-r}}.} ובפרט, מתקיים ∑ k = − n n r k = r − n ⋅ 1 − r 2 n + 1 1 − r . {\displaystyle \sum _{k=-n}^{n}r^{k}=r^{-n}\cdot {\frac {1-r^{2n+1}}{1-r}}.} אז לאחר הכפלת המונה והמכנה ב- r − 1 / 2 {\displaystyle r^{-1/2}} , נקבל: r − n − 1 / 2 r − 1 / 2 ⋅ 1 − r 2 n + 1 1 − r = r − n − 1 / 2 − r n + 1 / 2 r − 1 / 2 − r 1 / 2 . {\displaystyle {\frac {r^{-n-1/2}}{r^{-1/2}}}\cdot {\frac {1-r^{2n+1}}{1-r}}={\frac {r^{-n-1/2}-r^{n+1/2}}{r^{-1/2}-r^{1/2}}}.} במקרה זה r = e i x {\displaystyle r=e^{ix}} ולכן: ∑ k = − n n e i k x = e − ( n + 1 / 2 ) i x − e ( n + 1 / 2 ) i x e − i x / 2 − e i x / 2 = − 2 i sin ⁡ ( ( n + 1 / 2 ) x ) − 2 i sin ⁡ ( x / 2 ) = sin ⁡ ( ( n + 1 / 2 ) x ) sin ⁡ ( x / 2 ) {\displaystyle \sum _{k=-n}^{n}e^{ikx}={\frac {e^{-(n+1/2)ix}-e^{(n+1/2)ix}}{e^{-ix/2}-e^{ix/2}}}={\frac {-2i\sin((n+1/2)x)}{-2i\sin(x/2)}}={\frac {\sin((n+1/2)x)}{\sin(x/2)}}} כנדרש. PGL Dirichlet, "On the convergence of trigonometric series which serve to represent an arbitrary function between given limits" J. für Math., 4 (1829) pp. 157–169