הצגת המקור של הדף "אופטיקת פורייה"

'''אופטיקת פורייה''' היא ענף ב[[אופטיקה]] שהחל להתפתח בסוף ה[[המאה ה-19|מאה ה-19]]. תורה זו מתארת את מהלך האור על ידי שימוש ב[[גל|מכניקת הגלים]] ובכלים [[מתמטיקה|מתמטיים]], הלקוחים מ[[הנדסת חשמל]] ומ[[תורת האינפורמציה]]. היא מהווה כלי חישובי ותאורטי פופולרי בקרב פיזיקאים ומהנדסים כיום והיא שימושית מאוד באופטיקה, [[תקשורת]], [[מכ"ם]] ועוד. קרויה על שמו של ה[[מתמטיקאי]] וה[[פיזיקאי]] ה[[צרפת]]י [[ז'אן-בטיסט ז'וזף פורייה]] ועל שם ההתמרה האינטגרלית הנושאת את שמו - [[התמרת פורייה]].
[[קובץ:Analog_Signal_Processing.pdf|שמאל|350px|ממוזער|מערכת אופטית של שתי עדשות מדמה את תמונתו של ז'אן-בטיסט פורייה. עדשה ראשונה מייצרת את התמרת פורייה של האובייקט ב"מישור פורייה". עיבוד האות מתבצע במישור פורייה. לאחר מכן, עדשה שנייה מבצעת התמרת-פורייה נוספת ומתקבלת התמונה ביציאה]]

==רקע היסטורי==
===אופטיקה גלית===
עד [[המאה ה-19]] נהגו פיזיקאים לתאר את האור כזרם של חלקיקים. גישה זו נקבעה במאה ה-17 על ידי [[אייזק ניוטון]] והולידה את [[אופטיקה גאומטרית|האופטיקה הגאומטרית]] אשר נלמדת עד היום ומצליחה להסביר בצורה מדויקת למדי הרבה תופעות אופטיות "יומיומיות". למרות הצלחתה של האופטיקה הגאומטרית, במאה ה-19 גילו פיזיקאים כמו [[אוגוסטן ז'אן פרנל]], [[יוזף פראונהופר]], [[תומאס יאנג]] ואחרים כי האור עובר תופעות גליות של [[התאבכות]] ו[[עקיפה]], אשר אינן ניתנות להסבר על ידי המודל החלקיקי והאופטיקה הגאומטרית. לעומת זאת, תופעות אלה בהחלט ניתנות להסבר על ידי אימוץ התפיסה כי האור הוא [[גל]] המתפשט במרחב. בכך נזנח המודל הניוטוני של האופטיקה הגאומטרית לטובת האופטיקה הגלית. גישה פיזיקלית חדשה זו יכולה הייתה להסביר הן את התופעות המוכרות של האופטיקה הגאומטרית והן את התופעות, שבאותה עת נחשבו חדשות, של עקיפה והתאבכות. תורה זו הגיעה לבשלות בשנת [[1864]], עת הציג הפיזיקאי הסקוטי [[ג'יימס קלרק מקסוול]] בפני [[החברה המלכותית]] סט של ארבע [[משוואות דיפרנציאליות]] המתארות את התנהגות הגלים האלקטרומגנטיים. משוואות אלו מכונות היום "[[משוואות מקסוול]]" והן מהוות את הבסיס התאורטי לכל התיאור הגלי של האור.

על אף יעילותה ורב-גוניותה של האופטיקה הגלית, היו עוד כמה תופעות שלא ניתנות היו להסבר גם על ידה. אחת מהן הייתה הקטסטרופה באולטרה-סגול של [[קרינת גוף שחור]] והשנייה הייתה [[האפקט הפוטואלקטרי]]. ארבעים ואחת שנים מאוחר יותר, בשנת [[1905]], הציע [[אלברט איינשטיין]] הסבר לאפקט הפוטואלקטרי על ידי אימוץ מודל חלקיקי חדש לאור, אולם טען כי אין בכך כוונה לזנוח את ההנחה כי האור הוא גם גל. גישה זו, המקבלת שני מודלים שונים ונפרדים (המודל החלקיקי והמודל הגלי) בו-זמנית, נקראת היום "[[דואליות גל-חלקיק|הדואליות של האור]]" והיא מניחה כי לאור יש גם תכונות גליות אשר מתוארות היטב על ידי המודל הגלי וגם תכונות חלקיקיות המתוארות היטב על ידי המודל החלקיקי. גישה זו הולידה בתחילת המאה ה-20 את [[מכניקת הקוונטים|המכניקה הקוונטית]], אשר הצמיחה את [[אופטיקה קוונטית|האופטיקה הקוונטית]], המתארת את התנהגות אור במערכות מסדר גודל קטן מאוד - גודל שבו נכשלת האופטיקה הגלית.

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

===אופטיקה, תקשורת ואינפורמציה===
מאז שנות השלושים של המאה ה-20, הייתה התקרבות של פיזיקאים רבים העוסקים באופטיקה לענפי האינפורמציה והתקשורת מתחום הנדסת החשמל. מגמה זו מובנת לאור העובדה כי קיים דמיון רב בין מערכות תקשורת ומערכות הדמיה, בשני המקרים מטרתן היא לאסוף ולהוביל מידע. כלומר, שני סוגי המערכות מקבלים מידע ([[קלט]]) ופולטים מידע ([[פלט]]). עבור מערכות תקשורת המידע הוא בדרך כלל פונקציה זמנית ([[מתח חשמלי|מתח]] או [[זרם חשמלי|זרם]] התלויים בזמן) ועבור מערכות אופטיות המידע הוא פונקציה מרחבית (התפלגות מרחבית של עצמת האור או של אמפליטודת השדה האלקטרומגנטי). יתרה מזאת, מערכות תקשורת ומערכות אופטיות חולקות תכונות נוספות כמו ליניאריות ואינווריאנטיות. כל מערכת בעלת שתי תכונות אלו ניתנת לתיאור מתמטי על ידי [[ניתוח תדירויות]] ו[[אנליזת פורייה]].

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

==רקע מתמטי==
[[קובץ:Fourier-series01.png|שמאל|ממוזער|150px|בניית גל ריבועי באמצעות סכום של סינוסים.]]
הבסיס התאורטי לאופטיקת פורייה נובע מענף במתמטיקה הנקרא "[[טור פורייה|טוריי פורייה]]". לפי תורה זו, בהתקיימם של תנאים מסוימים, ניתן לרשום כל פונקציה [[טור (מתמטיקה)|כטור]] של [[פונקציה הרמונית|פונקציות הרמוניות]] (סינוסים וקוסינוסים). פונקציות הרמוניות אלה מקיימות את משוואת הגלים ולכן גם כל [[צירוף ליניארי]] שלהן מקיים אותה. לפיכך, אם נסתכל על אות כלשהו שמציית ל[[משוואת הגלים]] (גל אור, אות חשמלי), מורכב ככל שיהיה, ניתן יהיה לפרק אותו לסדרה של רכיבים - סינוסים וקוסינוסים פשוטים הנקראים "רכיבי פורייה". על רכיבים אלה ניתן לעשות מניפולציות בצורה פשוטה יחסית - פעולה זאת נקראת "עיבוד האות". לאחר עיבוד האות מרכיבים את כל הרכיבים מחדש ליצירת הפלט (תמונה במקרה של [[מערכת אופטית]] או אות חשמלי במקרה של מערכת תקשורת).

הפעולה המתמטית שמפרקת את הקלט לרכיביו נקראת [[התמרת פורייה]] והפעולה שמחברת את רכיבי פורייה מחדש נקראת "התמרת פורייה ההפוכה". עדשה מבצעת את התמרת פורייה על גל אור שנכנס אליה.
[[קובץ:Ряд Фурье.gif|250px|מרכז|המחשת טור פורייה]]
==התמרת פורייה דו-ממדית==
{{לשכתב|הפרק משתמש בטון דידקטי ולא אנציקלופדי|נושא=מדעי הטבע}}
לאורך הערך כולו נשתמש בהגדרה הבאה עבור [[התמרת פורייה]] של פונקציה דו־ממדית <math>g_{(x,y)}</math>:
{{ש}}
<math>\mathcal{F}\left\{g_{(x,y)} \right\}=G_{(f_x,f_y)}=\int_{-\infty}^{+\infty}{\int_{-\infty}^{+\infty}{g_{(x,y)}\mbox{exp}\left[-j2\pi\left( f_xx+f_yy \right) \right]}\,dx}\,dy</math>
{{ש}}
כאשר <math>j</math> היא [[מספר מדומה|היחידה המדומה]], <math>\left(f_x,f_y\right)</math> נקראות התדירויות המרחביות ויש להן יחידות הפוכות ליחידות של <math>\left(x,y\right)</math>.
{{ש}}
את ההתמרה ההפוכה נגדיר באופן הבא:
{{ש}}
<math>\mathcal{F}^{-1}\left\{G_{(f_x,f_y)} \right\}=g_{(x,y)}=\int_{-\infty}^{+\infty}{\int_{-\infty}^{+\infty}{G_{(f_x,f_y)}\mbox{exp}\left[j2\pi\left( f_xx+f_yy \right) \right]}\,df_x}\,df_y</math>
{{ש}}
לשתי התמרות אינטגרליות אלה תכונות שימושיות כמו ליניאריות, [[הפרדת משתנים]] וסימטריה והן מצייתות למספר משפטים מתמטיים שכדאי לקורא להכיר כמו משפט פרסבל, משפט ה[[קונבולוציה]], משפט אינטגרל פורייה ועוד.
{{ש}}
בנוסף, נציג את התמרת הנקל מסדר אפס, שהיא התמרת אנלוגית להתמרת פורייה המתאימה לפונקציות דו־ממדיות עם סימטריה רדיאלית
<math> g_{(r,\theta)}=g_{R(r)} </math>:
{{ש}}
<math> G_{(\rho)}=2\pi \int_{0}^{\infty}{rg_{R(r)}J_0 \left( 2\pi r \rho \right)}\,dr </math>
{{ש}}
כאשר <math>J_0</math> היא [[פונקציית בסל]] מסדר אפס ונתונה ע"י:
{{ש}}
<math>J_0\left(a\right)=\frac{1}{2\pi}\int_{0}^{2\pi}{\mbox{exp}\left[ -ja\cos\left( \theta-\phi \right) \right]}\,d\theta </math>
{{ש}}
בשלב זה נמליץ לקורא ללמוד ולהכיר התמרות פורייה/הנקל של פונקציות שימושיות ([[פונקציית מדרגה|מדרגת היחידה]], מפתח ריבועי, מפתח מעגלי, [[פונקציית דלתא של דיראק]], פונקציית "אוהל", [[מסרק דיראק|רכבת הלמים]] וכו'), אשר לא הובאו כאן כדי לא להכביד על קריאת הערך.

==מערכות ליניאריות==
מערכת מוגדרת כ[[מיפוי (מתמטיקה)|מיפוי]] של סט נתוני קלט לסט נתוני פלט. עבור המקרה של מערכות חשמליות, הקלט והפלט הוא פונקציית חד־ממדיות של הזרם או המתח התלויות בזמן; עבור המקרה של מערכות אופטיות, הקלט והפלט הן פונקציות דו־ממדיות ממשיות (עוצמה) או מרוכבות (השדה) המשתנות במרחב. נסמן את פעולת המערכת על ידי [[אופרטור]] מתמטי <math>\mathcal{S}\left\{\cdot \right\}</math> הפועל על פונקציית הקלט <math>g_1\left(x_1,y_1\right)</math> וממפה אותה לפונקציה <math>g_2\left(x_2,y_2\right)</math>
באופן הבא:
<math>g_2\left(x_2,y_2\right)=\mathcal{S}\left\{g_1\left(x_1,y_1\right)\right\}</math>
{{ש}}
[[קובץ:LinearSys.pdf|300px|מרכז|ממוזער|תיאור סכמטי של פעולת אופרטור]]

מתוך קבוצת כל האופרטורים נתמקד רק [[העתקה ליניארית|באופרטורים הליניאריים]]. ליניאריות של אופרטור אומרת שהפלט מאופרטור הפועל על [[צירוף ליניארי|קומבינציה ליניארית]] של שתי פונקציות (או יותר) הוא קומבינציה ליניארית של פלט האופרטור לאחר שפעל על כל פונקציה בנפרד. או בכתיב מתמטי - עבור כל שתי פונקציות
<math> p\left(x_1,y_1\right),q\left(x_1,y_1\right) </math>
ועבור כל שני סקלרים
<math> \alpha, \beta </math>
מתקיים:
{{ש}}
<math> \mathcal{S}\left\{ \alpha p\left(x_1,y_1\right)+\beta q\left(x_1,y_1\right) \right\}=\alpha\mathcal{S}\left\{p\left(x_1,y_1\right) \right\}+\beta\mathcal{S}\left\{ q\left(x_1,y_1\right)\right\} </math> .}{{ש}}
[[קובץ:LinearOp.pdf|מרכז|300px|ממוזער|תיאור סכימטי של פעולת אופרטור ליניארי]]

==ליניאריות ואינטגרל הסופרפוזיציה==
נסתכל על פעולת אופרטור ליניארי הפועל על פונקציית כניסה <math>g_1\left(x_1,y_1\right)</math>. אך מעט לפני כן, נשתמש בתכונה ידועה של פונקציית-דלתא כדי לכתוב את פונקציית הכניסה בצורה הבאה:
<math> g_1\left(x_1,y_1\right)=\int{\int{g_1\left(\xi,\eta\right)\delta\left(x_1-\xi,y_1-\eta\right)}\,d\xi}\,d\eta </math> .
{{ש}}
למציאת התגובה של המערכת יש להפעיל את האופרטור על הקלט
{{ש}}
<math> g_2\left(x_2,y_2\right)=\mathcal{S}\left\{\int{\int{g_1\left(\xi,\eta\right)\delta\left(x_1-\xi,y_1-\eta\right)}\,d\xi}\,d\eta \right\} </math> .
{{ש}}
משום שהאופרטור לא פועל על פונקציית משתני-הדמה שב[[אינטגרל]], ניתן להכניס אותו לתוך האינטגרל ולהפעילו על פונקציית הדלתא בלבד
<math> g_2\left(x_2,y_2\right)=\int{\int{g_1\left(\xi,\eta\right)\mathcal{S}\left\{\delta\left(x_1-\xi,y_1-\eta\right)\right\}}\,d\xi}\,d\eta </math> .
{{ש}}
בשלב הסופי, נגדיר את פונקציית [[תגובת הלם|התגובה להלם]] של המערכת
<math> h\left(x_1,y_1;\xi,\eta\right)=\mathcal{S}\left\{\delta\left(x_1-\xi,y_1-\eta\right)\right\} </math> .
{{ש}}
כעת, ניתן לרשום את הקשר בין הקלט והפלט של המערכת הליניארית באופן הבא
<math> g_2\left( x_2,y_2 \right)=\int{\int{g_1\left( \xi,\eta\right)h\left( x_1,y_1;\xi,\eta \right)}\,d\xi}\,d\eta </math> .
{{ש}} אינטגרל זה נקרא "אינטגרל הסופרפוזיציה" וחשוב לשים לב אינטגרל זה הוא בעצם [[קונבולוציה]]. ובכך אנו מקבלים עקרון חשוב מאוד באופטיקת פורייה - '''הפלט של מערכת אופטית שווה לקונבולוציה של פונקציית הכניסה עם פונקציית התגובה להלם של המערכת''' <math> g_2=g_1*h </math> .
{{ש}}
כאן ראוי להזכיר את משפט הקונבולוציה הטוען כי התמרת פורייה של קונבולוציה בין שתי פונקציות שווה למכפלת התמרות פורייה של כל פונקציה בנפרד (בהנחה שתנאי המשפט מתקיימים) - <math> G_2=G_1\cdot H </math>. כאשר
<math> G_2 </math>
היא התמרת פורייה של פונקציית הפלט
<math> g_2 </math>
והפונקציה
<math> G_1 </math>
היא התמרת פורייה של פונקציית הקלט
<math> g_1 </math>
והפונקציה
<math> H </math>
היא התמרת פורייה של פונקציית התגובה להלם של המערכת והיא נקראת גם [[פונקציית תמסורת|פונקציית התמסורת]] של המערכת.

==מעבר ממשוואות מקסוול הווקטוריות למשוואת הלמהולץ הסקלרית==
כאמור, הבסיס לכל המודל הגלי של האור מסוכם בצורה אלגנטית ופשוטה על ידי ארבעת [[משוואות מקסוול]]. משוואות אלו הן [[משוואה דיפרנציאלית|משוואות דיפרניציאליות]] וקטורית והן מצמדות בין [[שדה חשמלי|השדה החשמלי]] [[שדה מגנטי|והשדה המגנטי]]. עבור תווך ליניארי, [[איזוטרופיות|איזוטרופי]], [[הומוגניות (פיזיקה)|הומוגני]], לא [[נפיצה|נפיץ]] ולא-מגנטי משוואות מקסוול יכולות להכתב בצורה הבאה ([[מערכת היחידות הבין-לאומית|במערכת יחידות MKS]])
{{ש}}
<math> \nabla\times\vec{E}=-\mu\frac{\partial\vec{H}}{\partial t} </math>
{{ש}}
<math> \nabla\times\vec{H}=-\epsilon\frac{\partial\vec{E}}{\partial t} </math>
{{ש}}
<math> \nabla\cdot\epsilon\vec{E}=0 </math>
{{ש}}
<math> \nabla\cdot\mu\vec{H}=0 </math> .
{{ש}}
לפי [[אנליזה וקטורית|זהות וקטורית]] שימושית מתקיים: {{ש}}
<math> \nabla\times\left(\nabla\times\vec{E} \right)=\nabla\left( \nabla\cdot\vec{E}\right)-\nabla^2\vec{E} </math> .
{{ש}}
אם נציב בזהות זו את המשוואת השלישית ונציב אותה מחדש במשוואת מקסוול הראשונה, נקבל את [[משוואת הגלים]] עבור השדה החשמלי
<math> \nabla^2\vec{E}=\frac{n^2}{c^2}\frac{\partial^2 \vec{E}}{\partial t^2} </math> .{{ש}}
כש- <math> n=\sqrt{\frac{\epsilon}{\epsilon_0}} </math> הוא [[מקדם שבירה|מקדם השבירה]] של החומר ו- <math> c=\frac{1}{\sqrt{\mu_0 \epsilon_0}} </math> היא [[מהירות האור]] בריק.{{ש}}
[[משוואה]] וקטורית זו ניתנת להפרדה לרכיבים - כל רכיב של השדה החשמלי, <math> E_x, E_y, E_z </math>, מציית למשוואה זו בנפרד ולכן ניתן לפרק אותה לסט של שלוש משוואות סקלריות <math> \nabla^2E_i=\frac{n^2}{c^2}\frac{\partial^2 E_i}{\partial t^2} </math> כאשר
<math>i=x,y,z</math>.
{{ש}}
משום שכל גל ניתן לכתיבה על ידי סכום של [[גל מישורי|גלים מישוריים]] מהצורה
<math> \vec{u}\left(\vec{r};t \right)=\vec{u}_0\mbox{exp}\left[\pm j\left(\vec{k}\cdot\vec{r}-\omega t \right) \right] </math>
בעלי וקטורי-גל
<math> \vec{k} </math> שונים, יש טעם להציב את הגל המישורי במשוואת הגלים.{{ש}}
לשם כך יש לחשב את הנגזרת הזמנית ולהציבה במשוואת הגלים :
<math> \frac{\partial^2\vec{u}}{\partial t^2}=-\omega^2\vec{u} </math> .{{ש}}
לאחר הכנסת [[יחס נפיצה|יחס הדיספרסיה]] המתאים לתווך איזוטרופי והומוגני <math> \omega=vk</math> מקבלים את המשוואה
<math> \nabla^2\vec{u}+k^2\vec{u}=0 </math> .
{{ש}}
(תהליך דומה ניתן לעשות גם עבור השדה המגנטי).{{ש}}
משוואה זו נקראת [[משוואת הלמהולץ]] והיא הייתה מוכרת למהנדסי החשמל שעסקו ב[[קו תמסורת|קווי תמסורת]] ובתקשורת [[טלגרף]].
היתרון של עבודה עם משוואה זו במקום עם משוואת הגלים הוא שכעת אין פותרים משוואה עם תלות זמנית ומרחבית, אלא משוואה מרחבית בלבד.

==דיפרקציה במרחב חופשי==
הבעיה הבסיסית באופטיקה היא זו - נתונה פונקציית ה[[התפלגות]] של האמפליטודת השדה האלקטרומגנטי במישור הכניסה למערכת. נניח כי הגל הוא [[מונוכרומטיות|מונוכרומטי]] (בעל אורך גל בודד
<math> \lambda </math>
) ומתקדם בכיוון ציר z. פונקציית הכניסה היא
<math> u_0\left(x_0,y_0,z=0\right) </math>.
אנו מחפשים את פונקציית התפלגות השדה ביציאה מהמערכת, שבמקרה הנוכחי היא מרחב חופשי,
<math> u_1\left(x_1,y_1,z\right) </math>
, כאשר השדה מציית למשוואת הלמהולץ.
{{ש}}
לפי [[עקרון הויגנס|עקרון הויגינס-פרנל]], כל נקודה בפונקציית הכניסה היא מקור של [[גל כדורי]], כאשר
<math> u_0 </math>
נותנת את המשקל של כל מקור כזה. לפיכך, פתרון הבעיה הוא [[סופרפוזיציה]] של גלים כדוריים, אשר לכל אחד מהם מקבל משקל שונה:
{{ש}}
<math> u_1\left( x_1,y_1,z \right)=\frac{1}{j\lambda z}\int{ \int{u_0\left( x_0,y_0,z \right) \mbox{exp}\left[ j\frac{2\pi}{\lambda}\sqrt{ \left( x_1-x_0 \right)^2+\left( y_1-y_0 \right)^2+z^2 } \right] }\,dx_0 }\,dy_0 </math>.

==קירוב פרנל==
כאשר גודל המפתח של פונקציית הכניסה קטן מספיק ממרחק ההתקדמות של הגל, ניתן לקרב את השורש שבאקספוננט לפי [[טור טיילור]] באופן הבא:

<math> \sqrt{ \left( x_1-x_0 \right)^2+\left( y_1-y_0 \right)^2+z^2 }\approx z\left( 1+\frac{\left( x_1-x_0 \right)^2}{2z^2}+\frac{\left( y_1-y_0 \right)^2}{2z^2} \right) </math>
{{ש}}
לאחר הכנסת קירוב זה לאקספוננט שבאינטגרל מקבלים:
{{ש}}
<math> u_1\left( x_1,y_1,z \right)\approx \frac{1}{j\lambda z} \mbox{e}^{j\frac{2\pi}{\lambda}z} \mbox{e}^{j\frac{\pi}{\lambda z} \left( x_1^2+y_1^2 \right)} \int{ \int{u_0\left( x_0,y_0,z \right) \mbox{e}^{j \frac{\pi}{\lambda z} \left( x_0^2+y_0^2 \right) } \mbox{e}^{-j \frac{2\pi}{\lambda z} \left(x_0 x_1+y_0 y_1 \right) } }\,dx_0 }\,dy_0 </math>.
{{ש}}
בשלב זה ניתן לסמן
<math> f_x\equiv\frac{x_1}{\lambda z}, f_y\equiv\frac{y_1}{\lambda z} </math>.{{ש}}
בנוסף, אפשר להבחין כי האינטגרל הנ"ל, המכונה גם אינטגרל ה[[עקיפת פרנל|דיפרקציה של פרנל-כירכהוף]], הוא למעשה התמרת פורייה של הפונקציה
<math> u_0\mbox{e}^{j\frac{\pi}{\lambda z}\left( x_0^2+y_0^2 \right) } </math>. {{ש}}
נרשום זאת כך:
{{ש}}
<math> u_1\left( x_1,y_1,z \right)\approx \frac{1}{j\lambda z} \mbox{e}^{j\frac{2\pi}{\lambda}z}\mbox{e}^{j\frac{\pi}{\lambda z} \left( x_1^2+y_1^2 \right)} \mathcal{F} \left\{ u_0\mbox{e}^{j\frac{\pi}{\lambda z}\left( x_0^2+y_0^2 \right) } \right\} </math>.

==קירוב פראנהופר==
כאשר מרחק ההתקדמות [[גדול מספיק]] כדי לקרב את האקספוננט
<math> \mbox{e}^{j\frac{\pi}{\lambda z} \left( x_0^2+y_0^2 \right)} </math>
ל-1, אנו נותרים בלעדיו באינטרגרנד ואינטגרל פרנל-כירכהוף הופך להיות התמרת פורייה של
<math> u_0 </math>. כך מקבלים:
{{ש}}
<math> u_1\left( f_x,f_y,z \right)=\Omega\mathcal{F}\left\{ u_0\left( x_0,y_0 \right) \right\} </math>.
{{ש}}
כאשר:
{{ש}}
<math> \Omega=\frac{1}{j\lambda z}\mbox{e}^{jkz}\mbox{e}^{j\frac{k}{2z}\left(x_1^2+y_1^2 \right)} </math>
{{ש}}
<math> k=\frac{2\pi}{\lambda} </math>.
{{ש}}
משום שהעצמה תלויה בריבוע ערכה המוחלט של האמפליטודה ניתן לכתוב את התמונה ביציאה כך:
{{ש}}
<math> I_{out}=\left| u_1 \right| ^2 = \left| \Omega \right| ^2 \left| \mathcal{F}\left\{ u_0 \right\} \right| ^2 </math>.
{{ש}}
קירוב זה נקרא [[עקיפת פראונהופר|קירוב פראנהופר]] ובו טמון העיקרון החשוב ביותר באופטיקת פורייה והסיבה לשמה:
'''פונקציית היציאה מהמערכת היא התמרת פורייה של פונקציית הכניסה'''.

===דוגמאות===
====עקיפה ממפתח מלבני====
[[קובץ:Rect_aperture.jpg|שמאל|300px|ממוזער|מפתח מלבני והתמרת פורייה שלו]]
במשטח מישורי מחומר אטום לאור, קיים חור מלבני בעל רוחב a וגובה b המואר בתאורה אחידה. נחשב כאן את תמונת המפתח על מסך המרוחק מרחק z מן המפתח, לפי קירוב פראנהופר:
{{ש}}
ראשית, נרשום את פונקציית הכניסה -
<math> u_{in}\left( x_0,y_0,0 \right)=\mbox{rect}\left( \frac{x_0}{a}\right) \mbox{rect}\left( \frac{y_0}{b} \right) </math>.
לפי קירוב פראנהופר, כדי לחשב את פונקציית היציאה צריך למצוא את התמרת פורייה של הכניסה. במקרה זה היא:
{{ש}}
<math> u_1=\Omega\mathcal{F}\left\{u_0\right\}=\Omega ab\mbox{sinc}\left( af_x \right)\mbox{sinc}\left( bf_y \right) </math>.
כאשר:
<math> f_x=\frac{x_1}{\lambda z} f_y=\frac{y_1}{\lambda z} </math>.
{{ש}}
והעוצמה תהיה:
{{ש}}
<math> I_{out}=\frac{a^2b^2}{\lambda^2z^2}\mbox{sinc}^2 \left( \frac{ax_1}{\lambda z} \right)\mbox{sinc}^2\left( \frac{by_1}{\lambda z} \right) </math>

====עקיפה ממפתח עגול====
[[קובץ:Round_aperture.jpg|שמאל|300px|ממוזער|מפתח עגול והתמרת פורייה שלו]]
נחשב את התמונת ביציאה מהמערכת כשהמפתח הוא מעגלי בקוטר l. פונקציית הכניסה הפעם היא:
{{ש}}
<math> u_{in}=\mbox{circ}\left( \frac{r_0}{\frac{l}{2}} \right) </math>.
{{ש}}
כאשר
<math> r_0=\sqrt{x_0^2+y_0^2} </math>
{{ש}}
הפעם נשתמש בהתמרת הנקל לחישוב פונקציית היציאה.
{{ש}}
<math> u_{out}=2\Omega A \frac{J_1\left( \pi\rho l \right)}{\pi\rho l} </math> כש-A הוא שטח המעגל ו-
<math> \rho=\frac{r_1}{\lambda z} </math> עם
<math> r=\sqrt{x_1^2+y_1^2} </math>.
{{ש}}
עצמת השדה ביציאה במקרה זה תהיה:
{{ש}}
<math> I_{out}=4\frac{A^2}{\lambda^2 z^2}\left[ \frac{J_1\left( \pi\rho l \right)}{\pi\rho l} \right]^2 </math>

====עקיפה מסריג אמפליטודה סינוסואידי====
[[קובץ:sin_aperture.jpg|שמאל|300px|ממוזער|מפתח ריבועי עם שינוי אמפליטודה סינוסי והתמרת פורייה שלו]]
במקרה זה פונקציית הכניסה היא סינוס עם תדירות מרחבית
<math>f_0</math> לאורך ציר x ועומק מודולציה
<math> m </math>
ה"רוכב" על עוצמה קבועה בת יחידה אחת:
{{ש}}
<math> u_{in}=1+m\mbox{cos}\left( 2\pi f_0 x_0 \right) </math>.
במקרה זה התמרת פורייה של הכניסה תיתן:
{{ש}}
<math> u_{out}=\Omega\mathcal{F}\left\{ u_{in} \right\}=\Omega \left[ \delta\left(f_x,f_y\right)+\frac{m}{2}\delta\left( f_x-f_0,f_y \right)+\frac{m}{2}\delta\left( f_x+f_0,f_y \right) \right] </math>.
{{ש}}

זוהי תוצאה שימושית מאוד - "הולדתן" של שתי פונקציות הדלתא נותנת שניתן להשתמש בסריג אמפליטודה דק כמעין [[מפצל קרניים]] - קרן נכנסת מתפצלת לשלוש קרנים: קרן אחת ממשיכה באותה דרך, קרן שנייה המוסטת ימינה וקרן שלישית מוסטת שמאלה. כדי לממש רכיב בעל פעולה דומה ההמסתמך רק על אופטיקה גאומטרית יש לתכנן רכיב עבה, כבד ומסורבל.{{ש}}
באופן דומה, ניתן לתכנן רכיב בעל פונקציית שן-מסור שיעשה עבודה דומה (אם כי שם נקבל יותר משתי קרניים) או שן-מסור א-סימטרי, שבאמצעותו ניתן לקבל קרן חדשה אחת בלבד.
זו דוגמה טובה ליתרונותיה של אופטיקת פורייה בתכנון רכיבים אופטיים, ובעיקר בתכנון רכיבים דיפרקטיביים שלא ניתן היה לתכנן באמצעות אופטיקה גאומטרית.

==רכיבים דקים==
בחלק הקודם עסקנו במערכת "מרחב חפשי". בחלק זה נביא רכיבים אופטיים נוספים - אל כולם נתייחס כאל רכיבים "דקים". כלומר, אין בתוכם תופעות דיפרקציה משניות. בנוסף, לשם פשטות החישובים, נתעלם מהחזרות פנימיות ומתופעות לא ליניאריות (אם ישנן). {{ש}}
בדומה לסימונים בחישובי מערכות ליניאריות, מטרתנו היא לבטא את פעולת הרכיב כאופרטור t הפועל על פונקציית הכניסה: {{ש}}
<math> u_{out}=t\left(x,y,z\right)u_{in} </math>.

===לוח דיאלקטרי===
[[קובץ:Dielectric_slab.pdf|שמאל|200px]]
נסתכל על לוח [[חומר דיאלקטרי|דיאלקטרי]] שקוף, בעל מקדם שבירה n ועובי d.{{ש}}
ללא המעבר דרך הלוח, היה האור צובר איבר פאזה
<math> \mbox{e}^{jkd} </math>, אך משום שהאור עובר דרך הלוח הוא יקבל איבר פאזה
<math> \mbox{e}^{jknd} </math>.{{ש}}
לכן הפרש הפאזה ביחס לאור שהיה עובר בוואקום מכתיב שאופרטור הרכיב יהיה:{{ש}}
<math> t=\mbox{e}^{jk\left(n-1 \right)d} </math> .{{ש}}
כאשר הגל פוגע בלוח בזווית
<math> \theta </math>,
אז העובי האפקטיבי של הלוח הוא ב[[הקירוב הפראקסיאלי|קירוב הפראקסיאלי]]
<math> d_{eff}\approx\frac{d}{\mbox{cos}\theta} </math>. אם נציב זאת באופרטור הלוח נקבל:{{ש}}
<math> t=\mbox{e}^{jk\left(n-1 \right)\frac{d}{\mbox{cos}\theta}} </math> .{{ש}}

===פריזמה===
[[קובץ:thin_prism.pdf|שמאל|200px]]
נתונה [[מנסרה (אופטיקה)|פריזמה]] דקה (נקראת לפעמים גם "טריז" - wedge), בעלת מקדם שבירה n ובעלת זווית ראש
<math> \alpha </math>. נוכל לכתוב ביטוי עבור עובי הלוחית כפונקציה של y:{{ש}}
<math> d\left(y\right)=\mbox{tan}\left(\alpha\right)y </math>.{{ש}}
לכן אופרטור הפריזמה יהיה:{{ש}}
<math> t=\mbox{e}^{jk\left(n-1 \right)\mbox{tan}\left(\alpha\right)y} </math> .

===עדשה דקה===
[[קובץ:Negative.lens.svg|שמאל|200px]]
נחשב את הביטוי עבור [[עדשה]] דקה בעלת עובי מרבי a, מקדם שבירה n ורדיוסי עקמומיות
<math> R_1, R_2 </math> המצייתים להסכם הסימנים של נוסחת לוטשי-העדשות. לפי [[נוסחת לוטשי העדשות]] לעדשה כזו יהיה אורך מוקד
<math> f=\left[ \left( n-1 \right) \left( \frac{1}{R_1}-\frac{1}{R_2} \right) \right]^{-1} </math>.
{{ש}}
באופן דומה לשני הרכיבים הקודמים, נרשום את עובי העדשה ונציב בביטוי עבור הפרש הפאזה בין אור שעבר בעדשה ובין אור שעבר בוואקום:
{{ש}}
<math> d\left(x,y\right)=a-\frac{x^2+y^2}{2}\left( \frac{1}{R_1}-\frac{1}{R_2} \right) </math> .
{{ש}}
לכן פונקציית ההעברה (או ה"אופרטור") של העדשה תהיה:
{{ש}}
<math> t\left(x,y\right)=\mbox{e}^{j\left(n-1\right)d\left(x,y\right)}=\mbox{e}^{jka\left(n-1\right)}\mbox{e}^{-j\frac{k}{2f}\left(x^2+y^2\right)} </math> .
{{ש}}
בביטוי זה ניתן להבחין בשני איברים: האיבר הראשון - איבר פאזה גלובלית קבוע במרחב והוא פחות מעניין מביניהם. האיבר השני - פאזה ריבועית, היוצרת גל כדורי. זהו האיבר הגורם לגל מישורי להתכנס או להתבדר, בהתאם למוקד העדשה ומרחק ההתפשטות.

==שיטת האופרטורים==
לאחר שמצאנו ביטויים עבור רכיבים אופטיים ועבור התקדמות במרחב חופשי, נוכל כעת להרכיב ביטויים עבור מערכות מורכבות יותר, מרובות רכיבים. כפי שניתן לצפות, הביטוי עבור מערכת מרובת מקטעים של התקדמות במרחב חופשי יכול להיות ארוך ומסורבל וידרוש פתרון אינטגרלים רבים. כדי להתמודד עם אתגר זה, נציג נוטציה שימושית:{{ש}}
ראשית, נגדיר אופרטור מתאים לכל רכיב ולאחר מכן נציג את היחסים ביניהם.
נסמן את הגדלים הבאים:{{ש}}
<math> \rho=\sqrt{x^2+y^2}, k=\frac{2\pi}{\lambda}</math>
===אופרטור קנה-מידה (scaling)===
<math> \nu[a]f(x,y)=f(ax,ay) </math>
===אופרטור ההזזה===
<math> S[\vec{s}]f(x,y)=f(x-s_x,y-s_y) </math>
===אופרטור הפאזה הריבועית===
<math> Q[a]=e^{\frac{k}{2}a\rho^2} </math>
===אופרטור עדשה דקה בעלת מוקד f===
<math> Q[-1/f]=e^{l\frac{k}{2}\left(-\frac{1}{f}\right)\rho^2} </math>
===אופרטור ההתקדמות בתווך חופשי===
<math> R[z]=\frac{e^{jkz}}{jkz}Q[1/z]\nu[1/\lambda z]FQ[1/z]=e^{jkz}F^{-1}Q[-\lambda^2z]F </math> {{ש}}
כאשר F מציין התמרת פורייה ו-<math> F^{-1} </math> מציין את התמרת פורייה ההפוכה.{{ש}}

===יחסי אופרטורים===

{| class="wikitable"
|-
!| ||<math>\nu</math>|| F || Q || R
|- align="center"
| <math>\nu</math> || <math>\nu[t_2]\nu[t_1]=\nu[t_2t_1]</math> || <math>\nu[t]F=F\nu[1/t]</math> || <math>\nu[t]Q[c]=Q[t^2c]\nu[t] </math> || <math>\nu[t]R[d]=R[d/t^2]\nu[t] </math>
|- align="center"
| F || <math> F\nu[t]=\nu[1/t]F </math> || <math> FF=\nu[-1] </math> || <math> FQ[c]=R\left[\frac{-c}{\lambda^2}\right]F </math> || <math> FR[d]=Q[-\lambda^2d]F </math>
|- align="center"
| Q || <math> Q[c]\nu[t]=\nu[t]Q[c/t^2] </math> || <math> Q[c]F=FR[-c/\lambda^2] </math> || <math> Q[c_2]Q[c_1]=Q[c_2+c_1] </math> || <math> Q[c]R[d]=</math>{{ש}}
<math> R\left[(d^{-1}+c)\right]^{-1}\cdot</math>{{ש}}
<math>\nu[1+cd]\cdot </math> {{ש}}
<math> Q\left[ (c^{-1}+d)^{-1}\right] </math>
|- align="center"
| R || <math> R[d]\nu[t]=\nu[t]R[t^2d] </math> || <math> R[d]F=FQ[-\lambda^2d] </math> || <math> R[d]Q[c]=</math> {{ש}}
<math> Q\left[ (c^{-1}+d)^{-1} \right]\cdot</math>{{ש}}
<math>\nu\left[ (1+cd)^{-1} \right]\cdot</math>{{ש}}
<math>R\left[ (d^{-1}+c)^{-1} \right] </math>
|| <math> R[d_2]R[d_1]=R[d_2+d_1] </math>

|}
===דוגמאות===
====עדשה בודדת - מישור פורייה====
נתונה עדשה בודדת בעלת אורך מוקד f. נרצה לבדוק מהי הטרנספורמציה שיעבור גל שפוגע בעדשה ומתקדם מרחק של f.{{ש}}
אופרטור המערכת הוא:{{ש}}
<math> \tau=R[f]Q\left[-\frac{1}{f}\right] </math> .{{ש}}
נרשום את R בצורתו המפורשת:
<math> R[f]=Q \left[ \frac{1}{f} \right] \nu \left[\frac{1}{\lambda f} \right] FQ \left[ \frac{1}{f} \right] </math> ונציב בביטוי עבור אופרטור המערכת:{{ש}}
<math> \tau=Q\left[\frac{1}{f} \right]\nu\left[\frac{1}{\lambda f} \right] FQ\left[\frac{1}{f} \right]Q\left[-\frac{1}{f}\right] </math> .{{ש}}
ניתן כעת לאחד את שני ה־Q-ים העוקבים ולבטלם. נקבל את הביטוי:
<math> \tau=Q\left[\frac{1}{f} \right]\nu\left[\frac{1}{\lambda f} \right] F </math> . {{ש}}
מכאן ניתן ללמוד כי לאחר שהאור עובר בעדשה ומתקדם למישור המוקד הוא מקבל איבר פאזה ריבועית (שאינו בא לידי ביטוי בעוצמה) ועובר כיווץ/הגדלה, אך הדבר החשוב ביותר הוא שהאור עובר התמרת פורייה. משום כך קוראים למישור המוקד של העדשה בשם "מישור פורייה". כפי שנראה בהמשך, עובדת היותה של מערכת זו כרכיב שמייצר טרנספורם פורייה של הגל הנכנס תהיה שימושית מאוד לצורכי עיבוד נתונים אנלוגי, למשל עבור סינון מרחבי.

====הדמיה על ידי עדשה בודדת====
נתונה המערכת הבאה: עדשה בעלת מרחק מוקד f מדמה אובייקט המרוחק ממנה מרחק u לתמונה אשר נמצאת במרחק v. {{ש}}
לפי תנאי הדימות מתקיים
<math> \frac{1}{f}=\frac{1}{u}+\frac{1}{v} </math> . {{ש}}
אופרטור המערכת יהיה:
<math> \tau=R[v]Q\left[-\frac{1}{f}\right]R[u] </math> . {{ש}}
נרשום את האופרטורים R בצורתם המפורשת ונציב בביטוי עבור אופרטור המערכת לקבלת:{{ש}}
<math> \tau=Q\left[ \frac{1}{v}\right]\nu\left[ \frac{1}{\lambda v}\right] FQ\left[\frac{1}{v} \right] Q\left[-\frac{1}{f} \right] Q\left[\frac{1}{u} \right]\nu\left[\frac{1}{\lambda u} \right] FQ[\left[\frac{1}{u} \right] </math> . {{ש}}
נשתמש ביחסי האופרטורים בין שני Q-ים עוקבים ובתנאי הדימות עבור שלושת ה־Q-ים שבמרכז:{{ש}}
<math> Q\left[\frac{1}{v} \right] Q\left[-\frac{1}{f} \right] Q\left[\frac{1}{u} \right]=Q\left[ \frac{1}{v}+\frac{1}{u}-\frac{1}{f}\right]=Q[0]=1 </math>. איבר זה מתבטל. {{ש}}
אופרטור המערכת הוא כעת: {{ש}}
<math> \tau=Q\left[ \frac{1}{v}\right]\nu\left[ \frac{1}{\lambda v}\right] F\nu\left[\frac{1}{\lambda u} \right] FQ[\left[\frac{1}{u} \right] </math> . {{ש}}
ניתן להזיז את ה-F הראשון מקום אחד ימינה לפי היחס
<math>F\nu\left[\frac{1}{\lambda u} \right]=F\nu\left[ \lambda u\right] </math> ולאחר מכן, להשתמש ביחס בין שתי התמרות פורייה עוקבות
<math>FF=\nu[-1] </math> לקבלת הביטוי: {{ש}}
<math> \tau=Q\left[ \frac{1}{v}\right]\nu\left[ \frac{1}{\lambda v} \right]\nu\left[ -1 \right]\nu\left[ \lambda u \right]Q\left[\frac{1}{u} \right] </math>.{{ש}}
את שלושת ה
<math>\nu</math>-ים העוקבים ניתן לצמצם ל-
<math>\nu\left[ -\frac{u}{v} \right] </math> ולהחליף מקום בין <math>\nu</math> ו־Q לקבלת:{{ש}}
<math> \tau=Q\left[\frac{1}{v} \right]Q\left[\frac{u}{v^2}\right]\nu\left[ -\frac{u}{v} \right] </math> .{{ש}}
לבסוף, ניתן לאחד את שני ה־Q-ים לקבלת התוצאה הסופית:{{ש}}
<math> \tau=Q\left[\frac{u+v}{v^2} \right]\nu\left[ -\frac{u}{v}\right] </math>. {{ש}}

מסקנה: מערכת זו, בסה"כ מוסיפה איבר פאזה ריבועית, משנה את קנה-המידה של התמונה (מגדילה או מקטינה) והופכת את התמונה. איבר הפאזה הריבועית לא בא לידי ביטוי בפונקציית העוצמה (העצמה תלויה רק בריבוע [[ערך מוחלט|ערך המוחלט]]).
====עיבוד אותות אנלוגי====
[[קובץ:Analog_sig_pro_scheme.pdf|שמאל|350px|ממוזער|מערכת אופטית מדמה סריג שתי וערב. תמונת הסריג במישור פורייה היא סריג ריבועי של נקודות. במישור פורייה מציבים צמצם אופקי צר החוסם את כל הנקודות בעלות מחזוריות בכיוון ציר y. התמונה המתקבלת במישור היציאה היא התמרת פורייה של סדרת הנקודות שלא נחסמו - סדרת קווים ישרים עם מחזוריות בכיוון ציר x.]]
טרנפורם פורייה הפועל על פונקציית כניסה מפרק אותה לרכיביה הספקטרליים השונים, כשלכל רכיב כזה יש משקל שונה. כפי שראינו, עדשה יכולה למפות גל נכנס לטרנספורם פורייה שלו במישור המוקד שלה ("מישור פורייה"). משום כך, ניתן לעשות מניפולציות שונות על רכיבי פורייה של פונקציית הכניסה במישור פורייה. למשל, ניתן להנמיך את חלקם, או אפילו לחסום לחלוטין.{{ש}}
למשל, נחשוב על רשת שתי וערב של קווים ישרים: במישור פורייה של עדשה, נקבל את טרנספורם פורייה שלו. במקרה זה - סריג ריבועי של נקודות. על ידי הצבת צמצם אופקי נוכל לחסום את כל התדרים השייכים לרכיבי פורייה בעלי מחזוריות בכיוון y ולאפשר רק לחלק מרכיבי פורייה בעלי מחזוריות בכיוון x לעבור. כעת נציב עדשה נוספת. עדשה זו תמפה את התמונה במישור פורייה לטרנספורם פורייה שלה במישור המוקד של העדשה החדשה. כלומר, מפעילים את התמרת פורייה פעמיים וחוזרים לתמונה המקורית, אך הפעם ללא הרכיבים בעלי המחזוריות בציר y, אלא רק בציר x. כך ניתן לקבל את קווי האורך בלבד. באופן דומה, על ידי הוספת צמצם אורכי במישור פורייה, נוכל לסנן את קווי האורך ולקבל את [[קו רוחב|קווי הרוחב]] בלבד.{{ש}}
הטכניקה אינה מוגבלת רק לתמונות מחזוריות. ניתן לסנן באמצעותה רעשים, לחסום תדרים גבוהים (הממוקמים רחוק מהראשית במישור פורייה) או תדרים נמוכים (הממוקמים קרוב לראשית), או להעביר תחום תדרים סופי. לשם כך עלינו רק לחסום פיזית נקודות במישור פורייה.

==מערכות אופטיות בעלות מפתח סופי==
[[קובץ:Finite_lens.pdf||שמאל|300px|ממוזער]]
עד עתה, כל הרכיבים האופטיים נחשבו כרכיבים אינסופיים בשטחם. כאשר לרכיבים יש גודל סופי, למשל עבור עדשה בעלת מפתח מעגלי בקוטר d, גודלם הסופי משפיע על התמונה ביציאה כפי שנראה בדוגמה הבאה.{{ש}}
כאמור, נתבונן בגל הנפלט מאובייקט, הנמצא מרחק u משמאל לעדשה סופית בעלת מוקד f ומפתח סופי המיוצג על ידי הפונקציה
<math> p(x,y)</math>, אשר ממופה לתמונה במישור במרחק v מימין לעדשה. אופרטור המערכת במקרה זה ניתן לביטוי ע"י{{ש}}
<math>\tau=R[v]\cdot p(x,y)Q\left[-\frac{1}{f}\right]R[v]</math> . כאשר לפי תנאי הדימות מתקיים
<math> \frac{1}{f}=\frac{1}{u}+\frac{1}{v} </math> .{{ש}}
לאחר הצבה של הביטוי המפורש של R, תנאי הדימות וקצת אלגברת אופרטורים נקבל: {{ש}}
<math> u_{out}=\tau u_{in}\approx Q\left[\frac{1}{v} \right]\nu\left[\frac{1}{\lambda f} \right]F\left\{ p(x,y)\cdot F\nu\left[\lambda v \right]Q\left[\frac{1}{u} \right]u_{in} \right\} </math> . {{ש}}
לפי חוקי התמרת פורייה ומשפט הקונבולוציה מתקבל הביטוי הסופי:{{ש}}

<math> u_{out}\approx Q\left[\frac{1}{v} \right]\nu\left[\frac{1}{\lambda f} \right]\left\{ P(f_x,f_y)* \nu\left[-\lambda v \right]Q\left[\frac{1}{u} \right]u_{in} \right\} </math> . {{ש}}
כאשר
<math>P(f_x,f_y) </math>
היא התמרת פורייה של פונקציית המפתח p ו-
<math>*</math> הוא סימן הקונבולוציה.{{ש}}
ניתן לראות כי פעולת המערכת נותנת פונקציית יציאה שמורכבת מקונבולוציה של התמרת פורייה של פונקציית המפתח, שבמקרה של עדשה עגולה היא פונקציית בסל, עם פונקציית הכניסה לאחר שעברה כיווץ/הרחבה. ישנם גם איברי פאזה ריבועית שאותם ניתן לבטל בקלות על ידי הוספת עדשות בעלות אורך מוקד של u ו-v במישור הכניסה ובמישור היציאה בהתאמה. אם העדשה הייתה אינסופית אז פונקציית המפתח שלה הייתה מישור אינסופי, שהתמרת הפורייה שלו היא פונקציית דלתא. הקונבולוציה של פונקציית דלתא עם פונקציית הכניסה היא פונקציית הכניסה בעצמה. זו התוצאה שקיבלנו כשחישבנו את פעולת עדשה אינסופית.

==הארה לא קוהרנטית==
עד כה טיפלנו במערכות אופטיות ליניאריות המוארות בתאורה [[קוהרנטיות (פיזיקה)|קוהרנטית]].{{ש}}
קיבלנו שעבור מערכות אלו הפלט הוא קונבולוציה של הקלט עם פונקציית התגובה להלם של המערכת
<math>u_{out}=h*u_{in} </math> ושהעוצמה ביציאה היא
<math>I_{out}^{(coherent)}=\left|h*u_{in} \right|^2 </math>.{{ש}}
מערכת המוארת בתאורה לא קוהרנטית תהיה ליניארית בעוצמה ולא באמפליטודה כמו במקרה הקוהרנטי. העוצמה במוצא תהיה קונבולוציה של עוצמת התגובה להם עם עוצמת פונקציית הכניסה:
<math> I_{out}^{(incoherent)}=|h|^2*|u_{in}|^2 </math> . {{ש}}
לשם נוחות, נהוג להפעיל את משפט הקונבולוציה כדי לעבוד עם מכפלה במקום קונבולוציה:
<math>\mathcal{F}\left\{ I_{out}\right\}=\mathcal{F}\left\{ |h|^2\right\}\cdot \mathcal{F}\left\{ |u_{in}|^2\right\} </math> .{{ש}}
עוד נהוג לנרמל את כל הפונציות לפי ערכן בראשית:{{ש}}
<math> g_{out}=\frac{\mathcal{F}\left\{ I_{out}\right\}}{\mathcal{F}\left\{ I_{out}\right\}_{(0,0)}} </math> {{ש}}
<math> g_{in}=\frac{\mathcal{F}\left\{ I_{in}\right\}}{\mathcal{F}\left\{ I_{in}\right\}_{(0,0)}} </math> {{ש}}
<math> \hat{\mathcal{H}}=\frac{\mathcal{F}\left\{ |h|^2\right\}}{\mathcal{F}\left\{ |h|^2\right\}_{(0,0)}} </math> {{ש}}
ולרשום את הנוסחה הנ"ל בצורה:
<math> g_{out}=\hat{\mathcal{H}}\cdot g_{in}</math> .{{ש}}
הפוקציה
<math>\hat{\mathcal{H}}_{(f_x,f_y)} </math>
נקראת "פונקציית התמסורת הלא קוהרנטית" או "[[פונקציית תמסורת אופטית]]" של המערכת (Optical Tranfer Function - OTF) ויש לה משמעות גאומטרית:{{ש}}
עבור שני מפתחים שאחד נמצא במיקום
<math> \left( \frac{\lambda f_x}{2}v,\frac{\lambda f_y}{2}v \right) </math>
והשני נמצא במקום
<math> \left( -\frac{\lambda f_x}{2}v,-\frac{\lambda f_y}{2}v \right) </math>, כאשר v הוא המרחק למישור המוקד, פונקציית התמסורת האופטית מייצגת את היחס בין שטח החפיפה של שני המפתחים ובין השטח הכולל של שניהם.{{ש}}

==כושר הפרדה - הגבלת רזולוציה אופטית==
כאמור, סופיותו וצורתו של המפתח של המערכת האופטית משפיעים על התמונה ביציאה בכך שהתמונה היא קונבולוציה של התמרת פורייה של המפתח עם פונקציית הכניסה. עובדה זו, יוצרת בעיה של הבחנה בין שתי נקודות שונות במישור הכניסה. כשהעדשה הייתה אינסופית היא הייתה יכולה למפות שתי נקודות (פונקציות דלתא) שונות במישור הכניסה לשתי נקודות (פונקציות דלתא) שונות במישור היציאה. כלומר, הייתה לצופה המתבונן בתמונה אפשרות להבחין בין שתי נקודות שונות ולייחס אותן לשתי נקודות שונות במישור הכניסה. לשם הדוגמה, נחשוב על שני כוכבים רחוקים. מערכת אופטית בעלת מפתח אינסופי הייתה יכולה למפות את שניהם לנקודות שונות על גלאי, כך שצופה היה יכול להבין כי אלו הם שני כוכבים שונים. כאשר למערכת מפתח סופי, שתי נקודות במישור הכניסה ימופו לשתי פונקציות סופיות (למשל, במקרה של מפתח עגול - לשתי פונקציות Airy). משום שלפונקציות אלה יש רוחב, אם הן יהיו "קרובות מדי" הן יכולות להראות כמו פונקציה אחת "מרוחה" על פני השטח והצופה לא יוכל להבדיל אם זהו כתם או שתי נקודות מובדלות. אם נישאר עם דוגמת הכוכבים, אז צופה לא יוכל להבחין כי אלו שני כוכבים נפרדים ויטעה לחשוב כי זהו גוף אחד גדול.{{ש}}
בעיה זו נקראת בעיית הרזולוציה האופטית. הפיזיקאי זוכה [[פרס נובל]] [[ג'ון ויליאם סטראט ריילי|לורד ריילי]] קבע קריטריון לקביעת [[רזולוציה אופטית]] של מערכת אופטית, הנקרא קריטריון ריילי. ריילי קבע כי שני מקורות נקודתיים הם "מופרדים" על ידי מערכת אופטית בעלת מפתח עגול בעל [[קוטר]] D כאשר המרחק בין שני שיאי תבניות Airy של שני המקורות שווה למרחק של האפס הראשון של פונקציית Airy ממרכזה.{{ש}}
===דוגמה 1: הפרדה בין שני מקורות נקודתיים - הארה קוהרנטית===
נתונה מערכת אופטית בעלת מפתח עגול ברדיוס R ושני מקורות נקודתיים וקוהרנטיים המרוחקים במרחק
<math>\Delta </math> זה מזה.{{ש}}
מפתח המערכת מבוטא על ידי הביטוי הבא:
<math>p(r)=\mbox{circ}\left(\frac{r}{R} \right) </math> {{ש}}
והתמרת הפורייה שלה היא הפונקציה:
<math> P(f_x,f_y)=\frac{J_1\left(2\pi R\rho \right)}{R\rho} </math> כש-
<math> \rho=\sqrt{f_x^2+f_y^2}</math> . {{ש}}
התפלגות האמפילטודה בכניסה תבוטא כך:
<math> u_{in}(x,y)=\delta(x,y)+\delta(x-\Delta,y) </math> .{{ש}}
[[קובץ:Bessel func.jpg|שמאל|300px|ממוזער|פונקציות בסל]]
התמונה המתקבלת במישור פורייה של המערכת נתונה ע"י:
<math>u_{out}\propto P\left(\frac{x}{\lambda v},\frac{y}{\lambda v} \right)*u_{in}\left(-\frac{x}{M},-\frac{y}{M} \right) </math> {{ש}}
כש- v הוא מרחק העדשה ממישור התמונה, u הוא המרחק בין האובייקט והעדשה ו-M היא ההגדלה הליניארית
<math>M=\frac{v}{u} </math> .{{ש}}

חישוב הקונבולוציה נותן:{{ש}}
<math> u_{out}\propto \frac{J_1\left(2\pi R\sqrt{x^2+y^2}/\lambda v \right)}{\sqrt{x^2+y^2}/\lambda v}+ \frac{J_1\left(2\pi R\sqrt{(x-M\Delta)^2+y^2}/\lambda v \right)}{\sqrt{(x-M\Delta)^2+y^2}/\lambda v}</math>{{ש}}
כש-
<math>J_1</math> הוא [[פולינום]] בסל מסדר ראשון. {{ש}}

(מן הגרף המצורף, ניתן לראות כי האפס הראשון של פונקציית בסל מסדר 1 הוא בערך
<math>x_0\approx 2.41</math> . {{ש}}
העוצמה ביציאת מהמערכת תהיה אם כן-{{ש}}
<math> I_{out}^{(coherent)}\propto\left| \frac{J_1\left(2\pi R\sqrt{x^2+y^2}/\lambda v \right)}{\sqrt{x^2+y^2}/\lambda v}+ \frac{J_1\left(2\pi R\sqrt{(x-M\Delta)^2+y^2}/\lambda v \right)}{\sqrt{(x-M\Delta)^2+y^2}/\lambda v} \right|^2</math> {{ש}}
כאשר עובדים בתאורה קוהרנטית ישנה אפשרות להגביר את כושר ההפרדה. אם ניתן להוסיף איבר פאזה שרירותי לאחר המקורות, אז פונקציית העצמה תהיה: {{ש}}
<math> I_{out}^{(coherent)}\propto\left| \frac{J_1\left(2\pi R\sqrt{x^2+y^2}/\lambda v \right)}{\sqrt{x^2+y^2}/\lambda v}+e^{j\phi} \frac{J_1\left(2\pi R\sqrt{(x-M\Delta)^2+y^2}/\lambda v \right)}{\sqrt{(x-M\Delta)^2+y^2}/\lambda v} \right|^2</math> .{{ש}}
כפי שנראה בקרוב, במקרה של הפרש פאזה
<math> \phi=\pi</math> ניתן להגביר את כושר ההפרדה של המערכת - השדה במרכז תמיד מתאפס ולכן תמיד ניתן יהיה להבחין בין שני המקורות.

===דוגמה 2: הפרדה בין שני מקורות נקודתיים - הארה לא-קוהרנטית===
בהינתן אות המערכת ואותה הבעיה כמו בדוגמה הקודמת אך הפעם - עם תאורה לא-קוהרנטית, פונקציית העצמה ביציאה מהמערכת היא סכום ריבועי תבניות איירי: {{ש}}
<math> I_{out}^{(incoherent)}\propto \left| \frac{J_1\left( 2\pi R \sqrt{x^2+y^2}/\lambda v \right)}{\sqrt{x^2+y^2}/\lambda v} \right|^2 + \left| \frac{J_1\left( 2\pi R \sqrt{(x-M\Delta)^2+y^2}/\lambda v \right)}{\sqrt{(x-M\Delta)^2+y^2}/\lambda v} \right|^2 </math>

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

הגרף התחתון מתאר שני מקורות קוהרנטיים עם הפרש פאזה של
<math> \pi</math>
בניהם. ניתן ללמוד מגרף זה שהשדה בין שני המקורות תמיד מתאפס ולכן תמיד ניתן להפריד בניהם, אולם, היתרון בכושר ההפרדה במקרה זה בא עם החיסרון של חוסר דיוק במיקום.
[[קובץ:Two sources.gif|מרכז|ממוזער|500px|חתך העוצמה של תמונה של שני מקורות נקודתיים מופרדים. הנקודות האדומות מציינות את מיקומי המקורות כפי שהיו נראים במערכת אידיאלית.]]

== לקריאה נוספת ==
<div class="mw-content-ltr">
* Goodman, Joseph W. "Introduction to Fourier optics." 1968
</div>

[[קטגוריה:אופטיקה]]
[[קטגוריה:עיבוד אותות]]
[[קטגוריה:התמרת פורייה]]