Difference between revisions of "User:Tohline/Appendix/CGH/ParallelApertures2D"

From VistrailsWiki
Jump to navigation Jump to search
Line 103: Line 103:
   <td align="left">
   <td align="left">
<math>~
<math>~
[Z^2 + y_1^2  + x_1^2]^{1 / 2} \, .
L \biggl[1 - \frac{2(x_1 X_j + y_1 Y_k ) }{L^2} \biggr]^{1 / 2}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~
L \biggl[1 - \frac{(x_1 X_j + y_1 Y_k ) }{L^2}  \biggr] \, .
</math>
  </td>
</tr>
</table>
 
Hence, the double-summation expression for the amplitude at screen location (x<sub>1</sub>, y<sub>1</sub>) becomes,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~A(x_1, y_1)</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~ A_0 \sum_j \sum_k
a_{jk} e^{i\phi_{jk} } \cdot \exp\biggl\{ -i \biggl[ \frac{2\pi(x_1 X_j + y_1Y_k}{\lambda L} \biggr] \biggr\}
\, ,
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
where,
<div align="center">
<math>~A_0 \equiv e^{i(2\pi L/\lambda)}</math>.
</div>
When written in this form, it should be apparent why discrete Fourier transform techniques &#8212; specifically, 2D-FFT techniques &#8212; are useful tools for evaluation of the complex amplitude, A(x<sub>1</sub>, y<sub>1</sub>).
<font color="red"><b>NOTE:</b></font>  If x<sub>1</sub> and/or y<sub>1</sub> are ever comparable in size to Z &#8212; which may be the case for large apertures or for apertures tilted by nearly 90&deg; to the image screen &#8212; then the variation of L with image screen position cannot be ignored and, accordingly, the coefficient A<sub>0</sub> cannot be moved outside of the double summation.


=See Also=
=See Also=

Revision as of 22:36, 7 February 2020

CGH: 2D Rectangular Appertures that are Parallel to the Image Screen

This chapter is intended primarily to replicate §I.B from the online class notes — see also an updated Table of Contents — that I developed in conjunction with a course that I taught in 1999 on the topic of Computer Generated Holography (CGH) for a subset of LSU physics majors who were interested in computational science. This discussion parallels the somewhat more detailed one presented in §I.A on the one-dimensional aperture oriented parallel to the image screen.


Whitworth's (1981) Isothermal Free-Energy Surface
|   Tiled Menu   |   Tables of Content   |  Banner Video   |  Tohline Home Page   |


Utility of FFT Techniques

Consider the amplitude (and phase) of light that is incident at a location (x1, y1) on an image screen that is located a distance Z from a rectangular aperture of width w and height h. By analogy with our accompanying discussion in the context of 1D apertures, the complex number, A, representing the light amplitude and phase at (x1, y1) will be,

<math>~A(x_1, y_1)</math>

<math>~=</math>

<math>~ \sum_j \sum_k a_{jk} e^{i(2\pi D_{jk} /\lambda + \phi_{jk})} \, , </math>

where, here, the summations are taken over all "j,k" elements of light across the entire 2D aperture, and now the distance Djk is given by the expression,

<math>~D^2_{jk}</math>

<math>~\equiv</math>

<math>~ (X_j - x_1)^2 + (Y_k - y_1)^2 + Z^2 </math>

 

<math>~=</math>

<math>~ Z^2 + y_1^2 - 2y_1 Y_k + Y_k^2 + x_1^2 - 2x_1 X_j + X_j^2 </math>

 

<math>~=</math>

<math>~ L^2 \biggl[1 - \frac{2(x_1 X_j + y_1 Y_k ) }{L^2} + \frac{X_j^2 + Y_k^2}{L^2} \biggr] \, , </math>

and,

<math>~L</math>

<math>~\equiv</math>

<math>~ [Z^2 + y_1^2 + x_1^2]^{1 / 2} \, . </math>

If <math>~|X_j/L| \ll 1</math> and <math>~|Y_k/L| \ll 1</math> we can drop the quadratic terms in favor of the linear ones in the expression for Djk and deduce that,

<math>~D_{jk}</math>

<math>~\approx</math>

<math>~ L \biggl[1 - \frac{2(x_1 X_j + y_1 Y_k ) }{L^2} \biggr]^{1 / 2} </math>

 

<math>~\approx</math>

<math>~ L \biggl[1 - \frac{(x_1 X_j + y_1 Y_k ) }{L^2} \biggr] \, . </math>

Hence, the double-summation expression for the amplitude at screen location (x1, y1) becomes,

<math>~A(x_1, y_1)</math>

<math>~\approx</math>

<math>~ A_0 \sum_j \sum_k a_{jk} e^{i\phi_{jk} } \cdot \exp\biggl\{ -i \biggl[ \frac{2\pi(x_1 X_j + y_1Y_k}{\lambda L} \biggr] \biggr\} \, , </math>

where,

<math>~A_0 \equiv e^{i(2\pi L/\lambda)}</math>.

When written in this form, it should be apparent why discrete Fourier transform techniques — specifically, 2D-FFT techniques — are useful tools for evaluation of the complex amplitude, A(x1, y1).

NOTE: If x1 and/or y1 are ever comparable in size to Z — which may be the case for large apertures or for apertures tilted by nearly 90° to the image screen — then the variation of L with image screen position cannot be ignored and, accordingly, the coefficient A0 cannot be moved outside of the double summation.

See Also


Whitworth's (1981) Isothermal Free-Energy Surface

© 2014 - 2021 by Joel E. Tohline
|   H_Book Home   |   YouTube   |
Appendices: | Equations | Variables | References | Ramblings | Images | myphys.lsu | ADS |
Recommended citation:   Tohline, Joel E. (2021), The Structure, Stability, & Dynamics of Self-Gravitating Fluids, a (MediaWiki-based) Vistrails.org publication, https://www.vistrails.org/index.php/User:Tohline/citation