LSI系統的空域和頻域的表示

在前面的文章中提到過關於LSI系統在空域和頻域的表示方法。在頻域中,可以把傳遞函數看成是對輸入函數傅裏葉變換

F(xi) 的調製。這一過程被稱為線性濾波(linear filters),對輸入的頻譜特性進行濾波,從而得到一個理想的輸出。

可以把傳遞函數寫成: [H(xi ){	ext{ = }}{A_H}(xi ){e^{ - j{Phi _H}(xi )}}]

其中 [{A_H}(xi )] 是濾波的振幅譜, [{e^{ - j{Phi _H}(xi )}}] 是濾波的相位譜。振幅和相位譜的功能不同的,本文接下來會闡述。

1. 4-f成像系統

從本章開始會闡述衍射和「傅裏葉光學」這個領域成像的原理。在開始的時候,我們考慮一個特殊的系統,該系統可以在成像過程中討論濾波特性。

透鏡的焦距為f,u1在前焦面,u2在後焦面。u1和u2之間存在傅裏葉變換的關係。

如圖所示,透鏡的焦距為f,光波長為 lambda 。後焦面與前焦面之間的光波場的關係為:

[{u_2}({x_2},{y_2}) propto int {int_{ - infty }^infty  {{u_1}({x_1},{y_1}){e^{ - jfrac{{2pi }}{{lambda f}}({x_1}{x_2} + {y_1}{y_2})}}d{x_1}d{y_1}} } ]

令:

[xi  = frac{{{x_2}}}{{lambda f}}]

[eta  = frac{{{y_2}}}{{lambda f}}]

那麼有:

[{u_2}(xi ,eta ) propto int {int_{ - infty }^infty  {{u_1}({x_1},{y_1}){e^{ - j{	ext{2}}pi ({x_1}xi  + {y_1}eta )}}d{x_1}d{y_1}} } ]

這便是二維傅裏葉變換的形式,這個結果告訴我們,透鏡的後焦面上的光場正比於前焦面光場的傅裏葉變換。

既然瞭解了透鏡的傅裏葉變換特性,那麼我們接下來把兩個透鏡放在一起,形成一個4-f系統,如下圖所示。

4-f系統示意圖

第一個透鏡對 u_1 經行了一次傅裏葉變換得到了 u_2 ,第二個透鏡對 u_2 進行了一次傅裏葉變換得到了 u_3 ,也就是說 u_3u_1 傅裏葉變換的傅裏葉變換。那麼根據傅裏葉變換的性質,可以得到:

[egin{gathered}   {u_{	ext{3}}}({x_3},{y_3}) propto int {int_{ - infty }^infty  {{u_2}({x_2},{y_2}){e^{ - jfrac{{2pi }}{{lambda f}}({x_2}{x_3} + {y_2}{y_3})}}d{x_2}d{y_2}} }  hfill \    propto int {int_{ - infty }^infty  {{U_1}(frac{{{x_2}}}{{lambda f}},frac{{{y_2}}}{{lambda f}}){e^{ - jfrac{{2pi }}{{lambda f}}({x_2}{x_3} + {y_2}{y_3})}}dxi deta } }  hfill \    propto mathcal{F}left{ {{U_1}(frac{{{x_2}}}{{lambda f}},frac{{{y_2}}}{{lambda f}})} 
ight}left| {_{xi  = frac{{{x_3}}}{{lambda f}},eta  = frac{{{y_3}}}{{lambda f}}}} 
ight. hfill \    propto {u_1}( - lambda fxi , - lambda feta ) = {u_1}( - {x_3}, - {y_3}) hfill \  end{gathered} ]

根據推導的結果,4-f系統成像的倍率為-1。

  • 舉例——衍射光柵

假設一束平面波通過一個方形的光柵,通過後光場的分佈為:

[{u_1}({x_1},{y_1}) = rect(frac{x}{{5lambda }})*comb(frac{x}{{10lambda }})]

是週期為 10lambda 的週期函數。根據傅裏葉變換的性質,週期函數可以寫成:

[{u_1}({x_1},{y_1}) = sumlimits_{n =  - infty }^infty  {{c_n}{e^{jfrac{{2pi }}{{10lambda }}{x_1}}}} ]

[{U_1}(xi ,eta ) = sumlimits_{n =  - infty }^infty  {{c_n}delta (xi  - frac{n}{{10lambda }})delta (eta )} ]

那麼在中間平面 u_2

[{u_2}({x_2}) propto sumlimits_{n =  - infty }^infty  {{c_n}delta (frac{{{x_2}}}{{lambda f}} - frac{n}{{10lambda }})delta (frac{{{y_2}}}{{lambda f}})} ]

[{u_2}({x_2}) propto sumlimits_{n =  - infty }^infty  {{c_n}delta ({x_2} - frac{{nf}}{{10}})delta ({y_2})} ]

那麼在中間平面的 [{x_2} = nf/10]y_2=0 處,有不連續分佈的光強點,光強大小為 {c_n}^2。在現實的系統中,不存在無限窄的的由δ函數表示的點,會有有限的寬度(比如說愛裏斑),這些點和光柵的衍射級數有關。4-f系統的第二塊透鏡對中間平面再一次進行了傅裏葉變換,輸出平面的光強和輸入平面相同。

2.濾波

假設在4-f系統的中間平面放置一個模板,其傳輸方程為:

[T({x_2},{y_2}) = {A_t}({x_2},{y_2}){e^{ - j{Phi _t}({x_2},{y_2})}}]

加入這個模板後,在 u_3 平面的光波場為:

[{u_3}({x_3},{y_3}) propto mathcal{F}left{ {T(xi ,eta ){U_1}(xi ,eta )} 
ight}]

這個模板相當於是空間濾波器,可以調製原始光波場 u_1 的空間頻率。空間頻率振幅的衰減(或增大)由 [{A_t}(frac{{{x_2}}}{{lambda f}},frac{{{y_2}}}{{lambda f}})] 決定,相位的調製由 [{Phi _t}(frac{{{x_2}}}{{lambda f}},frac{{{y_2}}}{{lambda f}})] 決定。

2.1振幅濾波

假設模板是二值的,在不同的位置非0即1。考慮一種特殊的情況,模板有一條狹縫,傳輸方程為:

[t({x_2},{y_2}) = rect(frac{{{x_2}}}{b})]

這條縫允許通過的空間頻率範圍為:

[ - frac{b}{{2lambda f}} leqslant xi  leqslant frac{b}{{2lambda f}}]

在該範圍之外的頻率均不允許通過。這個傳輸方程如下圖所示,這種濾波被稱為低通濾波(LPF)。

低通濾波

經過4-f系統後,在後焦面的強度分佈如下圖的紅線圖形所示。因為過濾了高頻分量只保留了低頻分量,所以信號強度的分佈不是完美的方波。這個圖形與文章《傅裏葉光學(二)》中,由部分傅裏葉級數疊加而成的方波形狀相類似。

相類似的過程還有高通濾波(HPF,保留高頻分量,過濾低頻分量),帶通濾波(BPF,保留通帶頻率分量,過濾其他頻率)。

原始方波和經過濾波後的強度圖形

2.2相位濾波

在中間平面的模板只改變頻譜的相位而不改變振幅,傳輸方程為:

[T({x_2},{y_2}) = {e^{jPhi ({x_2},{y_2})}}]

假設有一塊相位屏,穿過它的光線會改變 [alambda /2] 的光程長度(optical path length,OPL), a 是一個常數。這就相當於相位屏為:

[{Phi _H}(xi ) = api step(x_2)]

根據不同的常數 a ,經行相位濾波後的圖形如下圖所示。

不同常數a的相位濾波

  • 線性相位濾波

[H(xi ) = {e^{ - j2pi {x_0}xi }}]

[G(xi ) = F(xi ){e^{ - j2pi {x_0}xi }}]

[g(x) = {mathcal{F}^{ - 1}}left{ {F(xi ){e^{ - j2pi {x_0}xi }}} 
ight} = f(x - {x_0})]

  • 弱相位濾波

當相位改變非常小,即 [left| {{Phi _H}} 
ight| ll 1] ,可以用泰勒級數將相位展開:

[H(xi ) = {e^{ - j{Phi _H}(xi )}} approx 1 - j{Phi _H}(xi )]

用傅裏葉逆變換得到衝激響應:

[h(x) approx delta (x) - j{phi _H}(x)]

其中 [{phi _H}(x) leftrightarrow {Phi _H}(xi )]

舉個例子,假設

[{Phi _H}(xi ) = dsin (2pi bxi ),left| d 
ight| ll 1]

這是弱正弦相位,進行傅裏葉逆變換得到:

[{phi _H}(x) =  - jfrac{d}{{2b}}left[ {delta (x - b) - delta (x + b)} 
ight]]

那麼弱相位下近似的衝激響應為:

[h(x) approx delta (x) + frac{d}{{2j}}left[ {delta (x - b) - delta (x + b)} 
ight]]

是正弦衍射光柵的0,+1和-1級。

2.3 同時進行振幅和相位濾波

對一個方波信號同時進行振幅和相位濾波後得到的圖形如下圖綠色的圖形所示。

藍色圖形是原始的方波信號,紅色圖形只對方波進行低通濾波,綠色圖形同時對方波進行振幅和相位濾波

推薦閱讀:

相關文章