Linear canonical transformation

Linear canonical transformation

Paraxial optical systems implemented entirely with thin lenses and propagation through free space and/or graded index (GRIN) media, are Quadratic Phase Systems (QPS). The effect of any arbitrary QPS on an input wavefield can be described using the linear canonical transform (LCT), a unitary, additive, four-parameter class of linear integral transform.

The former appeared a couple of times before Moshinsky and Quesne (1974) called attention to their significance in connection with canonical transformations in quantum mechanics. A particular case of the latter was developed by Segal (1963) and Bargmann(1961) in order to formalized Fok's boson calculus (1928). [K.B. Wolf, "Integral Transforms in Science and Engineering," Ch. 9:Canonical transforms, New York, Plenum Press, 1979.]

The LCT generalizes the Fourier, fractional Fourier, Laplace, Gauss-Weierstrass, Bargmann and the Fresnel transforms as particular cases.

Mathematical Definition

There are maybe several different ways to represent LCT. However, LCT can be viewed as a 2x2 matrix with determinant of the matrix is equal 1.:X_{(a,b,c,d)}(u) = sqrt{-i} cdot e^{-i pi frac{d}{b} u^{2 int_{-infty}^infty e^{-i 2 pi frac{1}{b} ut}e^{i pi frac{a}{b} t^2} x(t) , dt , when b e 0 ,

:X_{(a,0,c,d)}(u) = sqrt{d} cdot e^{-i pi cdu^{2 x(du) , , when b = 0 ,

:ad-bc = 1 , should be satisfied

pecial Cases of LCT

Since Linear Canonical Transform is a general term for other transforms, here are some examples of the special case in LCT.

Fourier Transform

Fourier Transform is a special case of LCT.

when egin{bmatrix} a & b \ c & dend{bmatrix} = egin{bmatrix} 0 & 1 \ -1 & 0end{bmatrix}

Fractional Fourier Transform

Fractional Fourier Transform is a special case of LCT.

when egin{bmatrix} a & b \ c & dend{bmatrix} = egin{bmatrix} cos heta & sin heta \ -sin heta & cos hetaend{bmatrix}

Fresnel Transform

Fresnel transform is equivalent to

when egin{bmatrix} a & b \ c & dend{bmatrix} = egin{bmatrix} 1 & lambda z \ 0 & 1end{bmatrix}z:distance ; lambda:wave length

=Additivity property of the WDF=

If we denote the LCT by O_F^{(a,b,c,d)} ,

i.e., X_{(a,b,c,d)}(u) = O_F^{(a,b,c,d)} [x(t)] ,


:O_F^{(a2,b2,c2,d2)} left { O_F^{(a1,b1,c1,d1)} [x(t)] ight } = O_F^{(a3,b3,c3,d3)} [x(t)] ,

where:egin{bmatrix} a3 & b3 \ c3 & d3end{bmatrix} = egin{bmatrix} a2 & b2 \ c2 & d2end{bmatrix}egin{bmatrix} a1 & b1 \ c1 & d1end{bmatrix}


Canonical transforms provide a fine tool for the analysis of a class of differential equations. These include the diffusion, the Schrödinger free-particle, the linear potential (free-fall), and the attractive and repulsive oscillator equations. It also includes a few others such as the Fokker-Planck equation. Although this class is far from universal, the ease with which solutions and properties are found makes canonical transforms an attractive tool for problems such as these. [K.B. Wolf, "Integral Transforms in Science and Engineering," Ch. 9&10, New York, Plenum Press, 1979.]

Wave propagation travel through air, lens, and dishes are discussed in here. All of the computations can be reduced to 2x2 matrix algebra. This is the spirit of LCT.

Electromagnetic Wave Propagation

If we assume the system look like this, the wave travel from plane xi, yi to the plane of x and y.

We can use Fresnel Transform to describe the Electromagnetic Wave Propagation in the air.: U_0(x,y) = - frac{i}{lambda} frac{e^{ikz{z} int_{-infty}^{infty} int_{-infty}^{infty} e^{j frac{k}{2z} [ (x-x_i)^2 -(y-y_i)^2 ] } U_i(x_i,y_i) , dx_i dy_i :k= frac {2 pi}{lambda}: , wave number; lambda : , wavelength; z : , distance of propagation

This is equivalent to LCT (shearing), when: egin{bmatrix} a & b \ c & d end{bmatrix}= egin{bmatrix} 1 & lambda z \ 0 & 1 end{bmatrix}

When the travel distance (z) is larger, the shearing effect is larger.

pherical lens

With the above lens from the image, and refractive index = n, we get:: U_0(x,y) = e^{ikn Delta} e^{-j frac{x}{2f} [x^2 + y ^2] } U_i(x,y) : f: , focal lenth Delta : , thickness of length

The distortion passing through the lens is similar to LCT, when: egin{bmatrix} a & b \ c & d end{bmatrix}= egin{bmatrix} 1 & 0 \ frac{-1}{lambda f} & 1 end{bmatrix}

This is also a shearing effect, when the focal length is smaller, the shearing effect is larger.

atellite Dish

Dish is equivalent to LCT, when: egin{bmatrix} a & b \ c & d end{bmatrix}= egin{bmatrix} 1 & 0 \ frac{-1}{lambda R} & 1 end{bmatrix}

This is very similar to lens, except focal length is replaced by the radius of the dish. Therefore, if the radius is larger, the shearing effect is larger.


If the system is considered like the following image. Two dishes, one is the emitter and another one is the receiver, and the signal travel through a distance of D.

First, for dish A (emitter), the LCT matrix looks like this::egin{bmatrix} 1 & 0 \ frac{-1}{lambda R_A} & 1 end{bmatrix}

Then, for dish B (receiver), the LCT matrix looks like this::egin{bmatrix} 1 & 0 \ frac{-1}{lambda R_B} & 1 end{bmatrix}

Last, we need to consider the propagation in air, the LCT matrix looks like this::egin{bmatrix} 1 & lambda D \ 0 & 1 end{bmatrix}

If we put all the effects together, the LCT would look like this::egin{bmatrix} a & b \ c & d end{bmatrix}=egin{bmatrix} 1 & 0 \ frac{-1}{lambda R_B} & 1 end{bmatrix}egin{bmatrix} 1 & lambda D \ 0 & 1 end{bmatrix}egin{bmatrix} 1 & 0 \ frac{-1}{lambda R_A} & 1 end{bmatrix}=egin{bmatrix} 1-frac{D}{R_A} & - lambda D \ frac{1}{lambda} (R_A^{-1} + R_B^{-1} - R_A^{-1}R_B^{-1}D) & 1 - frac{D}{R_B} end{bmatrix},

ee also

Other time-frequency transforms:

* Fractional Fourier transform
* continuous Fourier transform
* chirplet transform


* J.J. Ding, "time-frequency analysis and wavelet transform course note," the Department of Electrical Engineering, National Taiwan University (NTU), Taipei, Taiwan, 2007.
* K.B. Wolf, "Integral Transforms in Science and Engineering," Ch. 9&10, New York, Plenum Press, 1979.
* S.A. Collins, "Lens-system diffraction integral written in terms of matrix optics," "J. Opt. Soc. Amer." 60, 1168–1177 (1970).
* M. Moshinsky and C. Quesne, "Linear canonical transformations and their unitary representations," "J. Math. Phys." 12, 8, 1772–1783, (1971).
* B.M. Hennelly and J.T. Sheridan, "Fast Numerical Algorithm for the Linear Canonical Transform", "J. Opt. Soc. Am. A" 22, 5, 928–937 (2005).
* H.M. Ozaktas, A. Koç, I. Sari, and M.A. Kutay, "Efficient computation of quadratic-phase integrals in optics", "Opt. Let." 31, 35–37, (2006).
* Bing-Zhao Li, Ran Tao, Yue Wang, "New sampling formulae related to the linear canonical transform", "Signal Processing" '87', 983–990, (2007).
* A. Koç, H.M. Ozaktas, C. Candan, and M.A. Kutay, "Digital computation of linear canonical transforms", "IEEE Trans. Signal Process.", vol. 56, no. 6, 2383-2394, (2008).

Wikimedia Foundation. 2010.

Нужно решить контрольную?

Look at other dictionaries:

  • Canonical coordinates — In mathematics and classical mechanics, canonical coordinates are particular sets of coordinates on the phase space, or equivalently, on the cotangent manifold of a manifold. Canonical coordinates arise naturally in physics in the study of… …   Wikipedia

  • Canonical — is an adjective derived from . Canon comes from the Greek word kanon , rule (perhaps originally from kanna reed , cognate to cane ), and is used in various meanings. Basic, canonic, canonical : reduced to the simplest and most significant form… …   Wikipedia

  • Linear complex structure — In mathematics, a complex structure on a real vector space V is an automorphism of V that squares to the minus identity, −I. Such a structure on V allows one to define multiplication by complex scalars in a canonical fashion so as to regard V as… …   Wikipedia

  • Bogoliubov transformation — In theoretical physics, the Bogoliubov transformation, named after Nikolay Bogolyubov, is a unitary transformation from a unitary representation of some canonical commutation relation algebra or canonical anticommutation relation algebra into… …   Wikipedia

  • Projection (linear algebra) — Orthogonal projection redirects here. For the technical drawing concept, see orthographic projection. For a concrete discussion of orthogonal projections in finite dimensional linear spaces, see vector projection. The transformation P is the… …   Wikipedia

  • Generalized linear model — In statistics, the generalized linear model (GLM) is a flexible generalization of ordinary least squares regression. It relates the random distribution of the measured variable of the experiment (the distribution function ) to the systematic (non …   Wikipedia

  • Theorems and definitions in linear algebra — This article collects the main theorems and definitions in linear algebra. Vector spaces A vector space( or linear space) V over a number field² F consists of a set on which two operations (called addition and scalar multiplication, respectively) …   Wikipedia

  • Projective transformation — A projective transformation is a transformation used in projective geometry: it is the composition of a pair of perspective projections. It describes what happens to the perceived positions of observed objects when the point of view of the… …   Wikipedia

  • Trace (linear algebra) — In linear algebra, the trace of an n by n square matrix A is defined to be the sum of the elements on the main diagonal (the diagonal from the upper left to the lower right) of A, i.e., where aii represents the entry on the ith row and ith column …   Wikipedia

  • General linear group — Group theory Group theory …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”