菲涅爾轉換

**菲涅爾轉換**是線性標準轉換的一個特例，廣泛應用在光學領域上，主要在描述電磁波於近場區域在空氣中傳播的狀況，他也可以藉由克希荷夫繞射公式去做近似而得到。若是一個遠場傳遞的電磁波，則比較類似夫朗和斐繞射公式。對於一個近場的電磁波，我們可以用菲涅爾數去觀測其性質。菲涅爾數只要遠遠大於一的時候，就代表此繞射波是處於近場。菲涅爾繞射積分式我們可以使用下述的近似去推導(後面章節也有較為完整的描述)，可以寫成以下式子：

: \frac{F\theta^2}{4} \ll 1,

這裡\theta 是用 \theta \approx \frac{a}{L}所描述的極大值角度, a是孔徑的尺寸，L是孔徑與觀察屏之間的距離。

## 早期針對菲涅爾轉換的研究
追溯菲涅爾轉換為人知的時期，要回到十七世紀的時候，由一位義大利人Francesco Maria Grimaldi在他的專題論文-"light"所提出的。接著，Richard C. MacLaurin 藉由觀測光的傳遞、一條遠距離光源所產生的光束遇到一個有孔或狹長的障礙物時如何被影響，去解釋菲涅爾繞射現象。他使用了惠更斯－菲涅耳原理去做推導，發現光的波前從障礙物的狹縫中穿出且射到一個距離狹縫有一段距離的屏幕上的效果，跟光的波前直接穿過去是差不多的，看起來就好像沒有經過任何的物理障礙一樣。這個結果表示如果一個縫隙如果非常的狹長，具有明亮中心的繞射寬帶則可以被觀測到；如果這個縫隙相對來說比前一個狀況寬一些，則會看見明暗交錯的繞射寬帶。若這個隙縫繼續的變大，則明暗交錯的繞射寬帶會漸漸地變的不明顯，直到繞射的現象無法再被觀測到為止。

MacLaurin沒有提到關於光如果從一個小孔照耀出去的話，繞射的中心會有機率受到影響而看不見。但是他指出相反的情況-如果一個圓形物體的影子區域，可能經由繞射效果而使的此區域變明亮，這也是我們俗稱的阿拉戈點。

在MacLaurin的著作-"光學"中，Francis Weston Sears提供了一個菲涅爾所提出的數學近似法，僅僅用一些簡單的數學，就能拿來預測繞射現象的主要特色。藉由觀測有洞的障礙物到偵測屏幕的直線距離，以及知道入射波的波長，我們是有可能去計算出菲涅爾區域。這個區域的特性是，最內部區域是一個圓，接著會有一層一層的同心環所包圍。如果最內部的圓之直徑在屏幕上已經足夠顯示出第一個菲涅爾區域且偵測屏幕沒有被屏蔽的話，光在屏幕中心點的振幅就會比原本的大兩倍。如果最內部的圓之直徑在屏幕上已經足夠顯示出前兩個菲涅爾區域，則光在屏幕中心點的振幅會趨近於零。這表示菲涅爾繞射還是有機會出現一個暗帶中心。上述所提到的特性都可以藉由實驗觀測到以及量測到，而且可以對應到充足的理論，讓我們可以去計算出這些現象發生時，各個條件所對應到的值。

## 菲涅爾轉換的積分
要推導菲涅爾轉換，我們先考慮一個電場的繞射模式，給定x,y,z座標，其公式如下：

: E(x,y,z)={1 \over {i \lambda}} \iint_{-\infty}^{+\infty}{ E(x',y',0) \frac{e^{ikr}}{r}}dx'dy' 

其中， E(x',y',0)  是電磁波在原點的電場分布， r則是兩個座標點位置的距離，其值為\sqrt{(x-x')^2+(y-y')^2+z^2} ， i \,則為虛數的符號。
除非是最簡單的繞射形式，否則我們無法直接從此積分式得到分析解，所以在實作上我們通常是以數值方法去計算其值。

### 菲涅爾轉換的近似型
原積分式較為難處理的部分便是r變數，因此我們先對他採取代數轉換，引入一個新變數\rho，代換掉x,y維度的距離：

:\rho^2 = (x-x')^2+(y-y')^2 \,

因此r就可以被改寫成以下形式，並且用泰勒展開式去近似成以z和\rho表示

: r=  \sqrt{(x-x')^2+(y-y')^2+z^2} = \sqrt{\rho^2+z^2} = z \sqrt{ 1 + \frac{\rho^2}{z^2} } 
:: = z \left[ 1 + \frac{\rho^2}{2 z^2} - \frac{1}{8} \left( \frac{\rho^2}{z^2} \right)^2 + \cdots \right] 
:: = z + \frac{\rho^2}{2 z} - \frac{\rho^4}{8z^3}  + \cdots

因為我們有以下的性質：

:\sqrt{1+k} = (1+k)^{1/2} = 1 + \frac{k}{2} - \frac{k^2}{8} + \cdots，在此情況下k即為\frac{\rho^2}{z^2}

接著，我們考慮上述展開式的第三項以後，我們希望那之後的項都非常小，小到可以忽略，若要達成這項條件，因為r在原電場關係式出現於複數exponential的次方以及分母，在\rho \ll z的情況下，分母部分第三項以後可以直接省略；而對於複數exponential的次方部分，其必須遠小於一個週期，也就是2\pi 。

這表示我們希望有以下性質：

: k \frac{\rho^4}{8 z^3}  \ll 2 \pi ，又 k = { 2 \pi \over \lambda } \,，故可得 \frac{\rho^4}{z^3 \lambda} \ll 8 

將x,y,x',y'代換回來，我們可將不等式滿足的條件改寫成下式，這表示只要積分式中的x,y,x',y'符合下列不等式，我們便能近似r。

: \frac{[(x-x')^2+(y-y')^2]^2}{z^3 \lambda} \ll 8

近似結果為泰勒展開式的前兩項：

: r \approx z + \frac{\rho^2}{2 z} = z + \frac{(x-x')^2 +(y-y')^2}{2 z} 

同時我們將分母的r簡化，由於第一項會遠比第二項來的大，因此我們直接將r近似成z，並將原本的電場公式改寫如下：

: E(x,y,z)=\frac{e^{ikz}}{i \lambda z} \iint_{-\infty}^{+\infty} E(x',y',0)e^ \cdot e^{i  \pi \frac{d}{b} u^{2}} \int_{-\infty}^\infty e^{-i 2 \pi \frac{1}{b} ut}e^{i  \pi \frac{a}{b} t^2} x(t) \; dt \, ,   || 當 *b*&nbsp;≠&nbsp;0,
|-
| X_{(a,0,c,d)}(u) = \sqrt{d} \cdot e^{i  \pi cdu^{2}} x(du) \, , || 當 *b*&nbsp;=&nbsp;0.
|}
且ad-bc=1
若將線性標準轉換的a,b,c,d四個參數寫成一個二乘二的矩陣以及一個行列式限制條件來表示，如同矩陣\begin{bmatrix}
a & b \\
c & d
\end{bmatrix}
，其條件為ad - bc = 1，可以變換出許多衍生的形式。

根據上述公式，菲涅爾轉換為線性標準轉換的一個特例，表示形式為：

U_o(x,y)=-\ \frac{i}{\lambda} \frac{e^{ikz}}{z} \int_{-\infty}^{\infty}\int_{-\infty}^{\infty} U_i(x',y') e^{i\frac{k}{2z}[(x-x')^2+(y-y')^2]}\ dx'dy' 

其中，\lambda代表波長，k = \frac{2\pi}{\lambda}，z為電磁波傳遞方向的距離。

我們亦可把菲涅爾轉換當成是兩個一維的線性標準轉換的合成，表示如下：

U_o(x,y)=\ e^{ikz} \sqrt{\frac{1}{j\lambda z}}\int_{-\infty}^{\infty}e^{i\frac{k}{2z}(y-y')^2}\sqrt{\frac{1}{j\lambda z}}\int_{-\infty}^{\infty} U_i(x',y') e^{i\frac{k}{2z}(x-x')^2}\ dx'dy' 

若以矩陣形式來看，則為\begin{bmatrix}
a & b \\
c & d
\end{bmatrix}=\begin{bmatrix}
1 & \lambda z \\
0 & 1
\end{bmatrix}


## 菲涅爾轉換的其他形式
### 卷積
菲涅爾轉換的積分式也可以有其他種的表示方式，以達到使用其數學性質的目的，如果我們定義以下的函數：

:h(x,y,z) = \frac{e^{i k z}}{i \lambda z} e^{i \frac{k}{2 z} (x^2 + y^2)}

這個積分可以被表示成卷積的形式

:E(x,y,z) = E(x,y,0) * h (x,y,z) \, = \int_{-\infty}^{\infty}\int_{-\infty}^{\infty}  E(x',y',0) * h (x-x',y-y',z)\ dx'dy'

換句話說，我們可以將這個電磁波的傳遞，以線性濾波器的模型來表示。這也是為什麼我們一般會將h (x,y,z)稱為自由空間傳遞時的脈衝響應。

### 傅立葉轉換
另外一種常見的變換形式就是傅立葉轉換，如果在積分式中，我們將k用\lambda來取代，亦即 k = { 2 \pi \over \lambda } \,，並且將x,x',y,y'的關係展開：

:\begin{align}
  (x - x')^2 &= x^2 + x'^2 - 2 x x', \\
  (y - y')^2 &= y^2 + y'^2 - 2 y y',
\end{align}

那麼我們便可以將此積分式變換為一個二維的傅立葉轉換表示形式，接著，我們使用定義G(p,q)為函數g(x,y)的傅立葉轉換：

:G(p,q) = \mathcal{F} \left\{ g(x,y) \right\} \equiv \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} g(x,y) e^{-i 2 \pi (p x + q y)} \,dx\,dy, 

其中p,q所代表的是空間頻率，也就是所謂的波數的意思，如此一來，將g(x,y)= E(x',y',0) e^{i \frac{\pi}{\lambda z} (x'^2 + y'^2)}代入，我們便能將原本的菲涅爾轉換表示成下列形式：

:\begin{align}
  E(x,y,z) &= \frac{e^{i k z}}{i \lambda z} e^{i \frac{\pi}{\lambda z}(x^2 + y^2)}
              \left. \mathcal{F}\left\{ E(x',y',0) e^{i \frac{\pi}{\lambda z} (x'^2 + y'^2)} \right\}
              \right|_{p = \frac{x}{\lambda z},\ q = \frac{y}{\lambda z}} \\
           &=  \left\{ h(x,y) \right\} \cdot G(p,q) \big|_{p = \frac{x}{\lambda z},\ q = \frac{y}{\lambda z}},
\end{align}

換句話說，當我們運用此性質時，先計算g(x,y)的傅立葉轉換，再將p,q代回 \big( \tfrac{x}{\lambda z}, \tfrac{y}{\lambda z} \big)，並乘上
 h(x,y)。這樣子表示的好處是，當有一個程序可以讓我們知道目標函數的傅立葉變換，我們便可以用這個方法直接計算，而且此表示形式也與線性標準轉換較為類似。

## 離散一般化菲涅爾轉換
2006年，由I. Aizenberg和J.T. Astola所發表，是由線性卷積的離散轉換所發想而來。他們使用一個可被廣泛運用的代數方法，使得能在任意的向量空間中做訊號分析。I. Aizenberg和J.T. Astola找到了一個可被運用在離散傅立葉轉換的一般化菲涅爾轉換，並且得出有偶數次的離散傅立葉轉換有兩個一般化的菲涅爾轉換函數、奇數次的離散傅立葉轉換則只有一個一般化的菲涅爾轉換函數。而這些菲涅爾轉換主要是由Walsh轉換與聯合轉換所推導而來。

2015年，由Xing Ouyang等人發表的"Discrete Fresnel Transform and Its Circular Convolution"利用如離散傅立葉轉換以及正弦餘弦函數等離散三角轉換，運用他們的特性，去導出在無線週期的光柵系統下的離散菲涅爾轉換。相較於以前的離散菲涅爾轉換，Xing等人所得出的不會有衰減的現象。同時，菲涅爾轉換的圓形卷積也是第一次在這篇論文中被探討，他可以證明出兩個序列的圓形卷積等同於其中一個訊號與菲涅爾轉換的圓形卷積，這個性質除了可以幫助我們得出托爾伯特圖像的參數外，亦可在光學上有所運用，也可以對於資料處理和菲涅爾轉換的數值分析有幫助。

## 應用
* 影像加密
* 訊號分離
* 重建信號

## 參考