OSA's Digital Library

Optics Express

Optics Express

  • Editor: Michael Duncan
  • Vol. 12, Iss. 19 — Sep. 20, 2004
  • pp: 4347–4352
« Show journal navigation

Redistribution of output weighting coefficients for complex multiplexed phase-diffractive elements

Hua Liu, Zhenwu Lu, and Fengyou Li  »View Author Affiliations


Optics Express, Vol. 12, Issue 19, pp. 4347-4352 (2004)
http://dx.doi.org/10.1364/OPEX.12.004347


View Full Text Article

Acrobat PDF (146 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

The formation of multiplexed phase-only holograms with more weighted phase functions creates spurious cross terms and nonlinear scaling. We extend previously reported work [Appl. Opt. 25, 3767 (1986)] by proposing a normal method to analyze multiplexed holograms mathematically. We show that the output of holograms with any number weighted phase function can be written as a new linear combination for the original phase function with new weights. The relationship between the original weights and the new weights is developed for real-time optimization of hologram performance. We focus on the analysis of two and three multiplexed holograms to demonstrate the effectiveness of this approach.

© 2004 Optical Society of America

1. Introduction

We extend previously reported results [8

8. D. P. Casasent and W. A. Rozzi, “Computer-generated and phase-only synthetic discriminant function filters,” Appl. Opt. 25, 3767–3772 (1986). [CrossRef] [PubMed]

] by proposing a normal method to analyze mathematically multiplexed holograms with N-weighted phase functions. We show that the output can be written as a new linear combination for the original phase functions with new weights. The relationship between the original weights and the new weights is developed for real-time optimization of hologram performance. We use a fast and efficient way to fit polynomials to obtain the weights for multiplexed holograms that correspond to the desired output and efficiency. In particular we analyze a trifocal hologram with our method. In addition, a bifocal hologram can be used in lieu of a trifocal hologram under some conditions, and the results obtained with the trifocal hologram are identical with those reported in Ref. 8

8. D. P. Casasent and W. A. Rozzi, “Computer-generated and phase-only synthetic discriminant function filters,” Appl. Opt. 25, 3767–3772 (1986). [CrossRef] [PubMed]

, which corroborates our method. It is believed that holograms with less than N phase functions can be used without the need for analysis when we analyze a hologram with N-weighted phase functions.

2. Theory

We consider a linear combination of N phase functions exp(n ) (ϕn , where n = 1,2…N are two-dimensional functions) with real weights An we assume that i=1nAi2=1,, which defines a new function M exp(ia), where M is the amplitude and a is the phase:

Mexp(ia)=n=1NAnexp(iϕn)
(1)

We need to perform some simple transformations in Eq. (1) to obtain a useful expression:

exp(ia)=n=1NAnMexp(iϕn)
(2)

where

M=[A12+A22++AN2+2A1A2cos(ϕ1ϕ2)+2A1A3cos(ϕ1ϕ3)++2A1ANcos(ϕ1ϕN)
+2A2A3cos((ϕ1ϕ3)(ϕ1ϕ2))++2A2ANcos((ϕ1ϕN)(ϕ1ϕ2))
+⋯⋯
+2AN1ANcos((ϕ1ϕN)(ϕ1ϕN1))]12
(3)

For M > 0, a key point in our development is that 1/M can be considered a periodic function of (ϕ 1 -ϕ 2), (ϕ 1 -ϕ 3)⋯(ϕ 1 - ϕ N)with period 2π. This allows a Fourier series expansion to be performed:

M(β1β2βN1)=1M=m1mN1am1m2mN1exp(im1β1+im2β2++imN1βN1)
β1=ϕ1ϕ2βN1=ϕ1ϕN1
(4)

Here m 1, m 2m n-1 are integral numbers, and weights a m1m2mN-1 can be written as

am1m2mN1=1(2π)N102π02πM(im1β1imN1βN1)dβ1N1
(5)

Here two main properties of weights a m1m2mN-1can be obtained. First, weights a m1m2mN-1 are real, because, from Eq. (3), M′(β 1,β 2β N-1)=M′(-β 1,-β 2⋯-β N-1). Second, from Eqs. (3) and (5) a m1m2mN-1 depends only on An and n and is completely independent of ϕn . Inserting Eq. (4) into Eq. (2) and using the definition of β 1β N-1 1 given in Eq. (4) finally yield

exp(ia)=m1mN1am1mN1{A1exp[i(m1++mN1+1)ϕ1im1ϕ2imN1ϕN]+A2exp[i(m1++mN1)ϕ1i(m11)ϕ2imN1ϕN]++ANexp[i(m1++mN1)ϕ1im1ϕ2i£mN11£©ϕN]}
(6)

This new function leads to two consequences. First, the phase-only diffractive element corresponding to N-multiplexed phase functions can be written as a new linear combination of the original functions affected by the new weights an . Second, the phase-only operation can introduce spurious cross terms. Since the new weights of the output depend only on the original weight of the input, one can perform an excellent polynomial fit on their dependence.

3. Validation

If N = 3, according to Eq. (6) the phase-only diffractive element that corresponds to three multiplexed phase functions can be written as

exp(ia)=+(a00A1+a10A2+a01A3)exp(iϕ1)+(a00A2+a10A1+a11A3)exp(iϕ2)
+(a00A3+a01A1+a11A2)exp(iϕ3)+
(7)

It is clear that the output can be written as a new linear combination of the original functions and spurious cross terms introduced. The new weights depend only on the original weights. We define the output ratios as x 1 = (a 2/a 1)2 and x 2 = (a 3/a 1)2, the input ratios as y 1 = (A 2/A 1)2 and y 2 = (A 3/A 1)2. Then the relationship between the input ratios and the output ratios shown in Fig. 1 can be obtained quickly and efficiently by use of a fast Fourier transform (FFT) operation. The diffraction efficiency, defined as a12 + a22 + a32 , is plotted in Fig. 2 as a function of the desired output ratios.

Fig. 1. Input ratio versus the desired output ratio of the three multiplexed holograms.
Fig. 2. Diffraction efficiency versus output ratio of the three multiplexed holograms.

Note that the relationship between the input ratios and the output ratios can be written as polynomials:

y2=i=0kj=0imijx1jx2ij,y1=i=0kj=0inijx1jx2ij
(8)

where k is the degree of these polynomials. The polynomial coefficients can be determined by performing a least-squares fit of the polynomial on the dependence of input and output ratios. Here we chose k = 7 and the polynomials are

y1=0.00005096072781+3.94783x10.00873x210.19480x12+⋯⋯+11.96649x12x25+
4.82817x1x264.11001x27
(9)
y2=0.000050960727820.00873x1+3.94783x2+0.44625x12+⋯⋯+8.09256x12x25+
6.93143x1x26+2.49845x27
(10)

The maximum absolute error of y 1 and y 2 is less than 0.002. We could achieve higher accuracy by increasing the degree of the polynomial and the sampling number. Now we can control the weights of the output in real time by adjusting the weights of the input.

It should be noted that the mathematical result of the two multiplexed holograms can be obtained from that of the three multiplexed holograms with y 1 = 0 (results in x 1 = 0) or y 2 = 0 (results in x 2 = 0). So the relationship between the input ratio and the output ratio of the two multiplexed holograms shown in Fig. 3 can be obtained quickly from Fig. 1 [with x 2 = 0 in Fig. 1(a) or with x 1 = 0 in Fig. 1(b)]. The diffraction efficiency shown in Fig. 3(b) can be obtained from Fig. 2 with x 1 = 0 or x 2 = 0. It is clear that the numerical results are the same as those reported in Ref. [8

8. D. P. Casasent and W. A. Rozzi, “Computer-generated and phase-only synthetic discriminant function filters,” Appl. Opt. 25, 3767–3772 (1986). [CrossRef] [PubMed]

].

Fig. 3 (a) Input ratio and (b) diffraction efficiency versus the desired output ratio of two multiplexed holograms.

4. Simulation

As an example, we created a binary phase-only hologram with N = 3 by quantifying the two-dimensional phase function with our computer program. Many possibilities can be considered. We consider the particular case of three Fresnel lens with the same focal length f one equivalent to a center lens and each of the other two shifted from the origin by an amount a and b in the x direction. If x 1 = 0.5 and x 2 = 0.6, with Eqs. (9) and (10) we obtain y 1 = 0.5383 and y 2 = 0.6564. This implies that A 1 = 0.7623, A 2 = 0.4038, and A 3 = 0.5004. In such a case ϕ 1 = (k/2f)(x 2+y 2), ϕ 2 = (k/2f)⌊(x-a)2+y 2⌋ , and ϕ 3 = (k/2f)⌊(x+b)2+y 2⌋, where k is the wave number. Substituting the phase functions above into Eq. (6), we can write the resulting trifocal lens as

exp(ia)=m1m2am1m2{A1exp{ik2f[[x+(m1am2b)]2+y2]}exp[ik2f((m1am2b)2m1a2m2b2)]+A2exp{ik2f[[x+(m1am2ba)]2+y2]}exp[ik2f((m1am2ba)2m1a2m2b2+a2)]+A3exp{ik2f[[x+(m1am2b+b)]2+y2]}exp[ik2f((m1am2b+b)2m1a2m2b2+b2)]}
(11)

Here, substituting some pairs of numbers m 1 and m 2 such as m 1 = 0, m 2 = 0; m 1 = 1, m 2 = 0; m 1 = 0, m 2 = 1; m 1 = -1, m 2 = 0; m 1 = 0, m 2 = -1; m 1 = -1, m 2 = 1; m 1 = 1, m 2 = -1, into Eq. (11), we obtain the desired output orders that posit at x = -b, x = 0 and x = +a as shown in Eq. (7). At the same time, we also obtain the spurious adjustments that appear at positions x = (m 1 a-m 2 b), x = (m 1 a-m 2 b - a) or x = (m 1 a-m 2 b + b). Note that some spurious orders could coincide with the desired Fresnel lenses but with different phase functions. This could create a difference between the real weight of the output and the calculated weights. However, the problem can be compensated by taking these unwanted orders into account, when we calculate the relationship between the input ratios and the output ratios according to Eqs. (5) and (7).

A binary representation of a trifocal hologram is shown in Fig. 4(a) with a = 450 μm, b = 525 μm, and f= 1.138 m for a wavelength of 632.8 nm (He-Ne laser). The reconstruction of the hologram calculated by a Fourier transform of the hologram is shown in Fig. 4(b), in which the results are in agreement with the expected ratios. In addition, the locations of the peaks are at x = -525 μm, x= 0, and x = 450 μm, in agreement with Eq. (11). Some spurious orders appear at positions x = (m 1 a-m 2 b), x = (m 1 a-m 2 b - a), or x = (m 1 a-m 2 b + b).

Fig. 4. (a) Binary representation of a trifocal lens desired for x 1 = 0.5 and x 2 = 0.6 and (b) reconstruction.

5. Conclusion

In conclusion, the phase-diffractive element constructed from a linear combination of N-weighted phase functions has been studied mathematically. We have shown that the final diffractive element contains a linear combination of the original phase functions affected by new weights and that some undesirable additional terms were created. These new weights depend only on the original weights of the input. Polynomial fitting is used to control the performance of the holograms quickly and efficiently. As an example, we studied a phase-diffractive element constructed from a linear combination of three weighted phase functions and the codification of three multiplexed Fresnel lenses to discover their performance characteristics. This method is suitable for any complex multiplexed phase-diffractive element except for M = 0; however, M = 0 can be calculated by use of the polynomial fitting reported in this paper.

Acknowledgments

This study is supported by the National Natural Science Foundation (60078006) and the State Key Laboratory of Applied Optics of the Chinese Academy of Science.

References and links

1.

D. P. Casasent, “Unified synthetic discriminant function computational formulation,” Appl. Opt. 23, 1620–1627 (1984). [CrossRef] [PubMed]

2.

J. A. Davis, S. W. Connely, G. W. Bach, R. A. Lilly, and D. M. Cottrell, “Programmable optical interconnections with large fun-out capability using the magneto-optic spatial light modulator,” Opt. Lett. 14, 102–104 (1989). [CrossRef] [PubMed]

3.

J. A. Davis, D. M. Cottrell, R. A. Lilly, and S. W. Connely, “Multiplexed phase-encoded lenses written on spatial light modulator,” Opt. Lett. 14, 420–422 (1989). [CrossRef] [PubMed]

4.

S. Reichelt and H. J. Tiziani, “Twin-CGHs for absolute calibration in wavefront testing interferometry,” Opt. Commun. 220, 23–32 (2003). [CrossRef]

5.

M. Beyerlein, N. Lindlein, and J. Schwider, “Dual-wave-front computer-generated holograms for quasi-absolute testing of aspherics,” Appl. Opt. 41, 2440–2447 (2002). [CrossRef] [PubMed]

6.

J. A. Davis, E. A. Merrill, D. M. Cottrell, and R. M. Bunch, “Effects of sampling and binarization in the output of the joint Fourier transform correlator,” Opt. Eng. 29, 1094–1100 (1990). [CrossRef]

7.

J. L. Horner and P. D. Gianino, “Applying the phase-only filter concept to the synthetic discriminant function correlation filter,” Appl. Opt. 24, 851–855 (1985). [CrossRef] [PubMed]

8.

D. P. Casasent and W. A. Rozzi, “Computer-generated and phase-only synthetic discriminant function filters,” Appl. Opt. 25, 3767–3772 (1986). [CrossRef] [PubMed]

9.

R. R Kallman, “Optimal low noise phase-only and binary phase-only optical correlation filters for threshold detectors,” Appl. Opt. 25, 4216–4217 (1986). [CrossRef] [PubMed]

10.

E. Carcole, M. S. Millan, and J. Campos, “Derivation of weighting coefficients for multiplexed phase-diffractive elements,” Opt. Lett. 20, 2360–2362 (1995). [CrossRef] [PubMed]

OCIS Codes
(090.4220) Holography : Multiplex holography
(230.4110) Optical devices : Modulators

ToC Category:
Research Papers

History
Original Manuscript: August 12, 2004
Revised Manuscript: August 26, 2004
Published: September 20, 2004

Citation
Hua Liu, Zhenwu Lu, and Fengyou Li, "Redistribution of output weighting coefficients for complex multiplexed phase-diffractive elements," Opt. Express 12, 4347-4352 (2004)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-12-19-4347


Sort:  Journal  |  Reset  

References

  1. D. P. Casasent, �??Unified synthetic discriminant function computational formulation,�?? Appl. Opt. 23, 1620�??1627 (1984). [CrossRef] [PubMed]
  2. J. A. Davis, S. W. Connely, G. W. Bach, R. A. Lilly, and D. M. Cottrell, �??Programmable optical interconnections with large fun-out capability using the magneto-optic spatial light modulator,�?? Opt. Lett. 14, 102�??104 (1989). [CrossRef] [PubMed]
  3. J. A. Davis, D. M. Cottrell, R. A. Lilly, and S. W. Connely, �??Multiplexed phase-encoded lenses written on spatial light modulator,�?? Opt. Lett. 14, 420�??422 (1989). [CrossRef] [PubMed]
  4. S. Reichelt and H. J. Tiziani, �??Twin-CGHs for absolute calibration in wavefront testing interferometry,�?? Opt. Commun. 220, 23�??32 (2003). [CrossRef]
  5. M. Beyerlein, N. Lindlein, and J. Schwider, �??Dual-wave-front computer-generated holograms for quasi-absolute testing of aspherics,�?? Appl. Opt. 41, 2440�??2447 (2002). [CrossRef] [PubMed]
  6. J. A. Davis, E. A. Merrill, D. M. Cottrell, and R. M. Bunch, �??Effects of sampling and binarization in the output of the joint Fourier transform correlator,�?? Opt. Eng. 29, 1094�??1100 (1990). [CrossRef]
  7. J. L. Horner and P. D. Gianino, �??Applying the phase-only filter concept to the synthetic discriminant function correlation filter,�?? Appl. Opt. 24, 851�??855 (1985). [CrossRef] [PubMed]
  8. D. P. Casasent and W. A. Rozzi, �??Computer-generated and phase-only synthetic discriminant function filters,�?? Appl. Opt. 25, 3767�??3772 (1986). [CrossRef] [PubMed]
  9. R. R Kallman, �??Optimal low noise phase-only and binary phase-only optical correlation filters for threshold detectors,�?? Appl. Opt. 25, 4216�??4217 (1986). [CrossRef] [PubMed]
  10. E. Carcole, M. S. Millan, and J. Campos, �??Derivation of weighting coefficients for multiplexed phase-diffractive elements,�?? Opt. Lett. 20, 2360�??2362 (1995). [CrossRef] [PubMed]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Figures

Fig. 1. Fig. 2. Fig. 3
 
Fig. 4.
 

Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited