zoukankan      html  css  js  c++  java
  • [傅里叶变换及其应用学习笔记] 十五. 傅里叶变换在衍射上的应用

    这份是本人的学习笔记,课程为网易公开课上的斯坦福大学公开课:傅里叶变换及其应用。

    光的衍射(Diffraction)

    光在传播过程中,遇到障碍物或小孔(窄缝)时,它有离开直线路径绕到障碍物阴影里去的现象,这种现象称为光的衍射。衍射会产生明暗条纹或光环。

    衍射的形成有三个要素:

    1. 光源(source)

    2. 上面分布着小孔的衍射屏(aperture plane),光线只能通过小孔穿透到另一侧

    3. 接收小孔传播过来的光线,并在上面形成图形的的接收屏(image plane)

    image

    衍射屏与接收屏的距离决定了有两种衍射

    1. 衍射屏与接收屏的距离相比光波长较近,称为近场衍射(near-field),又成为菲涅耳衍射(Frenel diffraction)

    2. 衍射屏与接收屏的距离相比光波长较远,称为远场衍射(far-field),又称为夫琅禾费衍射(Fraunhoffer diffraction)

    我们平常所做的衍射实验还有一个附加条件:光源需要离衍射屏足够远,使得光可以平行经过衍射屏上的小孔,并且光波在衍射平面上有相同的相位。

    光的数学表示

    我们这里只研究光的波动性。因此假设光是震荡的电磁场,并且是单色光(单一频率)。

    在衍射屏上有光波表示为(波的复数表示形式)

    $Ee^{2pi i u t}$

    其中$E$表示电场强度,$ u$表示频率,$t$为时间,表明光波随时间波动。

    接收屏上的光波表示

    要求解这个问题需要用到惠更斯原理(Huyghens' Principle)。惠更斯原理讲的是,衍射屏上的每个小孔,都可以视为一个新的光源。

    image

    如上图,则通过$dx$的光波记为

    $E_0e^{2pi i u t}dx$

    现要求$P$点处光波的变化。光波在$P$点的变化有幅度与相位的变化,我们这里分析的是相位变化,因此忽略幅度变化。

    相位与光传播的距离$r$以及光的波长$lambda$相关。小孔$x$与$P$点间有$frac{r}{lambda}$个光波,因此相位(弧度)的变化为$frac{2pi r}{lambda}$。因此经由$x$到$P$点的光波为

    $dE = E_0e^{2pi i u t}e^{2pi ifrac{r}{lambda}}$

    因此,经过整个衍射屏后传播到$P$点的总光波为

    $E=displaystyle{int_{aperture}E_0e^{2pi i u t}e^{2pi ifrac{r}{lambda}}dx }$

    与$x$相关的部分只有$r$,因此

    $E=displaystyle{E_0e^{2pi i u t}int_ {aperture}e^{2pi i frac{r}{lambda}}dx }$

    现引入夫琅禾费近似公式

    image

    假设$rgg x$,则

    $r=r_0-xsin heta$

    把该等式代入上述光波$E$的式子,得

    $egin{align*}E
    &=E_0e^{2pi i u t}int_{aperture}e^{2pi ifrac{1}{lambda}(r_0-xsin heta)}dx\
    &=E_0e^{2pi i u t}e^{2pi ifrac{r_0}{lambda}}int_{aperture}e^{-2pi ifrac{xsin heta}{lambda}}dx
    end{align*}$

    由于我们只关心相位的相关部分,也就是积分内的相关部分,因此可以写成

    $E propto displaystyle{int_{aperture}e^{-2pi ixfrac{sin heta}{lambda}}dx}$

    令$p=frac{sin heta}{lambda}$,

    $E propto displaystyle{int_{aperture}e^{-2pi ixp}dx}$

    对于衍射屏上的小孔,我们可以用孔径函数$A(x)$来记录

    $A(x)=egin{cases}
    1 & ext{ , } xin aperture \
    0 & ext{ , } otherwise
    end{cases}$

    只有小孔才能使光波穿过,其余地方都不透光,即

    $Epropto displaystyle{int_{-infty}^{infty}A(x)e^{-2pi ipx}dx }$

    我们看到,等式右边是一个傅里叶变换,即

    $Epropto mathcal{F}A(p)quad ,quad p=frac{sin heta}{lambda}$

    结论:

    • 光的强度是孔径函数傅里叶变换的幅值(Intensity of the light is the magnitude of the Fourier Transform of aperture function)

    单缝衍射

    单缝衍射的截面的孔径函数是一个$Pi$函数

    image

    $A(x)=Pi_a(x)$

    那么其衍射图像的截面为

    $mathcal{F}A(x)=asincleft( frac{asin heta}{lambda} ight) quad,quad p=frac{sin heta}{lambda}$

    15_1

    我们肉眼观察到的光强度是其绝对值

    $|mathcal{F}A(x)|=left| asincleft( frac{asin heta}{lambda} ight) ight| $

    15_2

    如果光不是从小孔射进,而是一个点(point),即点光源,而接收屏离衍射屏足够远,那么有如下分析:在点上射出的光源用$delta$表示,其在接收屏的光为$mathcal{F}delta=1$,即光会均匀照亮接收屏。

    杨氏双缝实验(Young's double slits)

    image

    $A(x) = Pi_a(x-frac{b}{2})+Pi_a(x+frac{b}{2})$

    $mathcal{F}A(x) = a(sinap)2cos(pi bp)quad,quad p=frac{2pi heta}{lambda}$

  • 相关阅读:
    Kubernetes 集成研发笔记
    Rust 1.44.0 发布
    Rust 1.43.0 发布
    PAT 甲级 1108 Finding Average (20分)
    PAT 甲级 1107 Social Clusters (30分)(并查集)
    PAT 甲级 1106 Lowest Price in Supply Chain (25分) (bfs)
    PAT 甲级 1105 Spiral Matrix (25分)(螺旋矩阵,简单模拟)
    PAT 甲级 1104 Sum of Number Segments (20分)(有坑,int *int 可能会溢出)
    java 多线程 26 : 线程池
    OpenCV_Python —— (4)形态学操作
  • 原文地址:https://www.cnblogs.com/TaigaCon/p/5067082.html
Copyright © 2011-2022 走看看