## 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

There are several applications including optical pattern recognition [1], optical interconnections [2,3], three-dimensional display, and an absolute interferometric test of aspherics [4,5], in which several phase functions with different weights are multiplexed into a single phase-only or a binary phase-only hologram. In many cases, the multiplexed functions are added and coded onto a spatial light modulator (SLM) such as a magneto-optic spatial light modulator, twisted nematic liquid crystals, and transmissive matrix addressed ferroelectric liquid crystals. Nonlinearity occurs in these processes because each pixel of the SLM must encode the sum of these functions. Research has shown that spurious cross terms are formed and the weights of the output have a nonlinear relation with the weights of the encoded linear combination [6], which is a fundamentally different result compared with the desired output. This problem forces one to use a variety of complicated numerical techniques to compensate for the nonlinearities [7–9]. The conventional way to tackle these kinds of problem is to use various computer algorithms such as simulated annealing or iterative gradient approaches. However, these approaches are computer intensive and not intuitive. We propose a new and simple way to analyze the hologram constructed by adding two weighted phase functions [10]. However this research is limited to the analysis of two multiplexed phase-diffractive elements. The particular technique of this method cannot be extended to multiplexed holograms with more than two phase functions. To our knowledge, no previous research mathematically deduced the relationship between the encoded weights and the output weights with a normal method that is suitable for holograms with *N*-weighted phase functions.

We extend previously reported results [8] 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, 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(*iϕ*_{n}
) (*ϕ*_{n}
, where *n* = 1,2…*N* are two-dimensional functions) with real weights *A*_{n}
we assume that $\sum _{i=1}^{n}}{A}_{i}^{2}=1,$, which defines a new function *M* exp(*ia*), where *M* is the amplitude and *a* is the phase:

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

where

$$\phantom{\rule{9.5em}{0ex}}+2{A}_{2}{A}_{3}\mathrm{cos}\left(\left({\varphi}_{1}-{\varphi}_{3}\right)-\left({\varphi}_{1}-{\varphi}_{2}\right))+\cdots +2{A}_{2}{A}_{N}\mathrm{cos}\left(({\varphi}_{1}-{\varphi}_{N}\right)-\left({\varphi}_{1}-{\varphi}_{2}\right)\right)$$

$$\phantom{\rule{9.5em}{0ex}}+\mathrm{\cdots \cdots}\text{}$$

$$\phantom{\rule{9.5em}{0ex}}+{2{A}_{N-1}{A}_{N}\mathrm{cos}\left(\left({\varphi}_{1}-{\varphi}_{N}\right)-\left({\varphi}_{1}-{\varphi}_{N-1}\right)\right)]}^{\frac{1}{2}}$$

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:

$$\phantom{\rule{18em}{0ex}}{\beta}_{1}={\varphi}_{1}-{\varphi}_{2}\cdots {\beta}_{N-1}={\varphi}_{1}-{\varphi}_{N-1}$$

Here *m*
_{1}, *m*
_{2}⋯*m*
_{n-1} are integral numbers, and weights *a*
_{m1m2⋯mN-1} can be written as

Here two main properties of weights *a*
_{m1m2⋯mN-1}can be obtained. First, weights *a*
_{m1m2⋯mN-1} are real, because, from Eq. (3), *M*′(*β*
_{1},*β*
_{2}⋯*β*
_{N-1})=*M*′(-*β*
_{1},-*β*
_{2}⋯-*β*
_{N-1}). Second, from Eqs. (3) and (5) *a*
_{m1m2⋯mN-1} depends only on *A*_{n}
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

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 *a*_{n}
. 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

$$\phantom{\rule{5em}{0ex}}+\left({a}_{00}{A}_{3}+{a}_{0-1}{A}_{1}+{a}_{1-1}{A}_{2}\right)\mathrm{exp}\left(i{\varphi}_{3}\right)+\cdots $$

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 ${a}_{1}^{2}$ + ${a}_{2}^{2}$ + ${a}_{3}^{2}$ , is plotted in Fig. 2 as a function of the desired output ratios.

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

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

$$4.82817{\text{x}}_{1}{{\text{x}}_{2}}^{6}-4.11001{{\text{x}}_{2}}^{7}$$

$$6.93143{\text{x}}_{1}{{\text{x}}_{2}}^{6}+2.49845{{\text{x}}_{2}}^{7}$$

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].

## 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*/2*f*)(*x*
^{2}+*y*
^{2}), *ϕ*
_{2} = (*k*/2*f*)⌊(*x*-*a*)^{2}+*y*
^{2}⌋ , and *ϕ*
_{3} = (*k*/2*f*)⌊(*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

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*).

## 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]