我們提出了一種處理傅里葉變換的方法,其并不需要二次多項(xiàng)式相位項(xiàng)的抽樣,而是用解析的方法處理。我們提出該理論的同時(shí)也給出了幾個(gè)例子證明其潛力。
0BrAgv"3a_ B1i&HoGbz 1.簡(jiǎn)介
<44A*ux d:8c}t2X 物理光學(xué)建模需要頻繁地從空間轉(zhuǎn)換到角頻域,反之亦然。這可以由電場(chǎng)和磁場(chǎng)分量的傅里葉變換得到。所以,快速傅里葉變換(FFT)算法成了快速
物理光學(xué)建模的支柱[1]。FFT技術(shù)的數(shù)值計(jì)算量與場(chǎng)分量復(fù)振幅所需采樣點(diǎn)的數(shù)量近似成線性關(guān)系。在光學(xué)中,我們經(jīng)常處理有強(qiáng)波陣面相位的場(chǎng)分量,例如:球形。但是由于2π模,平滑的波陣面相位的復(fù)抽樣導(dǎo)致了大量的數(shù)值計(jì)算工作,甚至在FFT中也是如此。
GV0-"9uwX~ Ku`u%5< 2.理論
`M[o.t 2.1 場(chǎng)的表征:提取二次相位
~A@HW!*Z@ wfQ^3HL 我們從空間域的符號(hào)開始,在本文中我們使用符號(hào)
![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP510154EW.png)
對(duì)應(yīng)6個(gè)場(chǎng)分量,也就是V = (E, H):
uD4j.% qAW?\*n5N &>A<{J@VL ![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP51016053R.png)
(1)
]x5+v0 在公式1中,我們假設(shè)場(chǎng)
![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP510163bR.png)
有兩部分:
衍射場(chǎng)
![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP5101F3L8.png)
和一個(gè)平滑的波陣面相位exp(iψ(ρ))。對(duì)于得到的結(jié)果,我們從波陣面相位中提取二次相位exp(iψ(ρ))并且將余下的部分認(rèn)為是余項(xiàng)場(chǎng)
![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP5101JJ11.png)
。假設(shè)exp(iψ(ρ))可由其實(shí)數(shù)系數(shù)C和D = (Dx, Dy)給出:
SIZZFihcYh v%@)I_6[P ![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP5101T0K9.png)
(2)
顯然,在強(qiáng)二次相位情況中,全場(chǎng)
![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP510163bR.png)
比余項(xiàng)場(chǎng)需要更多的抽樣量。所以,我們的目標(biāo)是通過(guò)FFT且無(wú)二次相位項(xiàng)exp(iψ(ρ))抽樣的情況下,計(jì)算V(ρ)的傅里葉變換。
qV5lv-p N~|Z@pU" 2.2.半解析傅里葉變換
-]Y@_T.C p6X-P%s 從卷積定理可知:
$*+IsP! +SNjU"x ![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP5101941608.png)
(3)
Yfr4<;% 通常來(lái)說(shuō),項(xiàng)
![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP5102025X4.png)
必須進(jìn)行數(shù)值計(jì)算處理。另一方面,從數(shù)學(xué)
角度[2]我們可知:
MqjdW e+<'=_x { ![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP510205K61.png)
(4)
"'74GY8, 適用于任何復(fù)
![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP5102145527.png)
,只要R{a} ≥ 0且a ≠ 0。
Seh(G 在該數(shù)學(xué)工具的幫助下,項(xiàng)κ[exp(iψ(ρ))]的解析表征可以推導(dǎo)出來(lái):
=/Ph]f9 ![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP5102241111.png)
(5)
其中:
2 9#jKh ![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP510240CE.png)
(6)
其中常數(shù)項(xiàng)
![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP5102430118.png)
。
-dS@l'$ 將公式5帶入公式3,通過(guò)改變卷積和傅里葉變換積分的階次,我們發(fā)現(xiàn)
![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP5102G4C6.png)
可以表示為:
):lH ![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP5102626107.png)
(7)
其中:
~@$RX:p ![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP5102I5318.png)
(8)
這里,
![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP5102Q3462.png)
和坐標(biāo)項(xiàng)
![](http://www.infotek.com.cn/uploads/allimg/190805/1-1ZP5102UR95.png)
。公式7-8是半解析傅里葉變換的數(shù)學(xué)表達(dá)式。它表示全場(chǎng)的FFT可被兩個(gè)余項(xiàng)場(chǎng)的FFT替代。
.$]-::& fj97_Q= 3.數(shù)值
仿真 W1_.wN$,5 1ne3CA= 這些概念在物理光學(xué)建模和設(shè)計(jì)軟件Wyrowski
VirtualLab Fusion[3]中實(shí)現(xiàn)。
hQ (84u .'PS L 3.1.有效性測(cè)試1:純二次相位
CUnBi?