Optical system for performing complex fourier transforms

Patent metadata

Field	Value
Publication number	US-11073860-B2
Application number	US-201815938994-A
Country	US
Kind code	B2
Filing date	Mar 28, 2018
Priority date	Sep 28, 2015
Publication date	Jul 27, 2021
Grant date	Jul 27, 2021

Abstract

A method of performing a complex Fourier transform of an input function including amplitude and phase information, including decomposing the input function into a plurality of sub-functions, wherein the Fourier transform of each sub-function includes an amplitude function and a phase function in which the phase is constrained to a plurality of possible phase values. The phase function of the Fourier transform of each sub-function is determined with an optical system that measures the amplitude function of an optical Fourier transform of the sub-function and changes in the amplitude function of the optical Fourier transform caused by applying a perturbation function to the sub-function. The determined phase functions and the measured amplitude functions are combined for each of the sub-functions to form the complex Fourier transform of the input function.

First claim

The invention claimed is:

1. A method of performing a complex Fourier transform of an input function comprising amplitude and phase information, comprising: decomposing the input function into a plurality of sub-functions with a processor, wherein the complex Fourier transform of each of the sub-functions comprises an amplitude function and a phase function in which a phase of the phase function is constrained to a plurality of possible phase values; determining, with the processor, the phase function of the Fourier transform of each of the sub-functions with an optical system that measures the amplitude function of an optical Fourier transform of each of the sub-functions and changes in the amplitude function of the complex Fourier transform caused by applying a perturbation function to the sub-function, and combining, with the processor, the determined phase functions and the measured amplitude functions for each of the sub-functions to form the complex Fourier transform of the input function.

2. The method claim 1 , wherein the step of decomposing the input function into the plurality of sub-functions comprises:
determining a real component and an imaginary component of the input function;
decomposing the real component of the input function into a first plurality of sub-functions;
decomposing the imaginary component of the input function into a second plurality of sub-functions.

3. The method claim 1 , wherein the step of measuring changes in the amplitude function comprises:
adding the perturbation function to the sub-function to form a perturbed sub-function;
performing an optical Fourier transform of the perturbed sub-function and detecting an amplitude function of the Fourier transform of the perturbed sub-function; and
comparing the measured amplitude function of the Fourier transform of the sub-function to the amplitude function of the Fourier transform of the perturbed sub-function.

4. The method of claim 1 , wherein the phase of the phase function constrained to two of the possible phase values, and wherein a difference between the two possible phase values for at least one of the sub-functions is (2n+1)π radians, where n is an integer.

5. The method of claim 4 , wherein a first of the plurality of sub-functions is a first even function and a first odd function, and wherein a second of the plurality of sub-functions is a second even function and a second odd function.

6. The method of claim 1 , wherein the complex Fourier transform of the sub-function is a pixelated representation of the complex Fourier transform of the sub-function comprising a spatial array of elements, wherein each of the elements comprise an amplitude value and a phase value.

7. The method of claim 1 , further comprising forming the complex Fourier transform of each of the sub-functions using the determined phase function and measured amplitude function for each of the sub-functions.

8. The method of claim 1 , wherein the step of combining the determined phase functions and measured amplitude functions comprises:
adding together the respective complex Fourier transforms of a first of the plurality of sub-functions to form the Fourier transform of a real component of the input function; and
adding together the respective complex Fourier transforms of a second of the plurality of sub-functions to form the Fourier transform of an imaginary component of the input function,
wherein the Fourier transform of the real component of the input function and the Fourier transform of the imaginary component of the input function form the complex Fourier transform of the input function.

9. The method of claim 1 , wherein the changes in the amplitude function of the optical Fourier transform caused by applying the perturbation function to the sub-function to generate a perturbed sub-function are determined by subtracting the measured amplitude function of the Fourier transform of the sub-function from the amplitude function of the Fourier transform of the perturbed sub-function to form a difference function.

10. The method of claim 9 , further comprising determining a sign of the difference function.

11. The method of claim 1 , wherein the step of measuring the amplitude function of the optical Fourier transform of the sub-function comprises:
representing the sub-function on a spatial light modulator;
illuminating the spatial light modulator to form spatially modulated light;
Fourier transforming the spatially modulated light using a Fourier transform lens; and
detecting a spatial intensity distribution at the Fourier plane of the Fourier transform lens, or
representing the sub-function on the spatial light modulator;
illuminating the spatial light modulator to form the spatially modulated light; and
detecting the spatial intensity distribution in an optical far field.

12. The method of claim 11 , further comprising square-rooting the spatial intensity distribution.

13. The method of claim 11 , wherein the step of detecting the spatial intensity distribution comprises detecting the spatial intensity distribution using a photodetector array.

14. The method of claim 1 , wherein the perturbation function is a delta function or a discrete approximation of the delta function.

15. The method of claim 1 , wherein the Fourier transform of the perturbation function comprises the amplitude function wherein the phase of the phase is one of the plurality of possible phase values.

16. A device for performing a complex Fourier transform of an input function comprising amplitude and phase information, comprising:
a processor arranged to decompose the input function into a plurality of sub-functions, wherein the complex Fourier transform of each of the sub-functions comprise an amplitude function and a phase function in which a phase is constrained to a plurality of possible phase values; and
an optical system arranged to measure, for each of the sub-functions the amplitude function of an optical Fourier transform of the sub-function and changes in the amplitude function of the optical Fourier transform caused by applying a perturbation function to the sub-function;
wherein the processor is further arranged to:
determine, for each of the sub-functions the phase function of the complex Fourier transform of the sub-function from the amplitude function of the optical Fourier transform of the sub-function and the changes in the amplitude function of the optical Fourier transform caused by applying the perturbation function to the sub-function; and
combine the phase functions and the amplitude functions of the optical Fourier transforms of each of the sub-functions to form the complex Fourier transform of the input function.

17. The device of claim 16 , wherein the optical system further comprises:
a spatial light modulator arranged to display each of the sub-functions;
a light source arranged to illuminate the spatial light modulator to form spatially modulated light;
a Fourier transform lens arranged to receive the spatially modulated light and Fourier transform the spatially modulated light, and
a photodetector array arranged to detect a spatial intensity distribution at a Fourier plane of the Fourier transform lens.

18. The device of claim 16 , wherein the optical system further comprises:
a spatial light modulator arranged to display each of the sub-functions;
a light source arranged to illuminate the spatial light modulator to form spatially modulated light;
a photodetector array arranged to detect a spatial intensity distribution in the optical far field.

Assignees

• Cambridge Entpr Ltd

Inventors

• Macfaden Alexander Joseph

Classifications

• G06E3/005

  using electro-optical or opto-electronic means · CPC title

• G06F17/14 Primary

  Fourier, Walsh or analogous domain transformations {, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms (for correlation function computation G06F17/156; spectrum analysers G01R23/16)} · CPC title

• G06E3/003 Primary

  forming integrals of products, e.g. Fourier integrals, Laplace integrals, correlation integrals; for analysis or synthesis of functions using orthogonal functions · CPC title

• G06F17/141

  Discrete Fourier transforms · CPC title