绕射光学元件definitionofdiffractiveopticalelementsdoe.pptVIP

設計繞射光學元件之模擬退火法的研究 林正峰及邱華楠 南台科技大學光電系 2005.12.10 繞射光學元件 Definition of Diffractive Optical Elements (DOE) 1D or 2D 表面起伏型週期性相位光柵 Several μm ~ several tens ofμm d =~ 1m m Laser beam (plane wave) Period: M x N phase elements DOE Transmittance function g(x,y)=exp( i f (x,y)) Intensity | F{g(x,y)} |2 Lens Diffraction pattern References: [1]. ”Introduction to Fourier Optics”, Joseph W. Goodman , McGRAW-HELL (1996) DOE 成像系統示意圖 我們所定義的繞射元件是一般純相位的元件，其會因元件表面的起伏變化或內部的折射率的變化，而對於射入光的波前產生相位調變。 start initialization T0 ,L0 T1 ,L1 Tn ,Ln End A sequence of Markov chain 1. Input diffraction pixel matrix (N 2) 2. Set design pattern 3. Set parameters : a , b , g0 , m . 4. Generate a transmittance function (g(x,y)) with random phase 5. Compute initial cost function (Cini ) 6. Generate initial temperature (T0) 7. Algorithm stop criteria. 模擬退火法 模擬退火法的基本流程 References: [2].”A quantitative analysis of the simulated annealing algorithm: A case study for the traveling salesman problem.”, E. H. L. Aarts, J. H. M. Korst, P. J. M. van Laarhoven, Journ. of Statistical Physics 50, p189-206. (1988) [3]. ”計算方法叢書－非數值並行算法（第一冊）模擬退火算法”, 康立山等著, 科學出版社. (1997) The cost function is globally minimized. If the elapsed time is infinite, simulated annealing can obtain global minimum value. The real simulated annealing can only obtain near global minimum value . 使用模擬退火法的目的 1.Randomly select a pixel and randomly change its phase to a different value 2.Compute new cost function Probability P(ΔC)=exp(-ΔC/T) Accept the change Cold = C new Unacceptable change End A Markov chain N = P(ΔC ) N P(ΔC ) The N is a uniformly distributed random number in [0,1] The cost function is reconstruction error Yes No The iteration to reach length of Markov chain start DC 0 DC 0 References: [2].”A quantitative analysis of the simulated annealing algorithm: A case study for the traveling salesman problem.”, E. H. L. Aarts, J. H. M. Korst, P. J. M. va