Ray Propagation Based on the Gaussian Imaging Equation in Ideal Lens Systems …(IG)
First, we derive the equations relating to imaging. Then, we show the matrix representation of geometric ray propagation by an ideal lens, based on the Gaussian imaging equation.
*This page deals with a ray propagation model based on geometrical optics. It does not include the effects of wave-optical diffraction, such as that of Gaussian beams.
The flow of explanation on this page is based on Toshiro Kishikawa's "Introduction to Optics for User Engineers" [1] . However, the derivations and descriptions have been independently restructured by the author.
table of contents:
This page will proceed with the discussion based on the following rules.
Rule A: Sign convention
Regarding dimension notation, in accordance with Appendix: Definitions of Direction and Sign in One-Dimensional Vectors and Angles , the up-down direction is considered positive when pointing upwards, the left-right direction is considered positive when pointing to the right, and the rotational direction is considered positive when turning counterclockwise.
Rule B: Definition of the paraxial region
All rays discussed on this page are assumed to be paraxial rays with sufficiently small angles. That is, when the angle is θ, the following approximation is applied:
sin θ = θ, tan θ = θ, cos θ = 1
Rule C: Parameter notation (color coding)
For the color coding of each parameter, quantities given by the settings (in the recurrence relation model) are shown in blue, and quantities given by the calculation formula are shown in black.
・Derivation of the imaging relationship formula
The manner in which a ray of light from an object point passes through an optical element with a focal length f and is focused onto an image point is illustrated below.
(i) When f > 0, z < 0,
* Apply Rule A (Symbol Convention) and Rule B (Definition of Paraxial Region).
(ii) When f > 0, 0 < z < f,
* Apply Rule A (Symbol Convention) and Rule B (Definition of Paraxial Region).
(iii) When f > 0, f < z,
* Apply Rule A (Symbol Convention) and Rule B (Definition of Paraxial Region).
(iv) When f < 0, z < f
* Apply Rule A (Symbol Convention) and Rule B (Definition of Paraxial Region).
(v) When f < 0, f < z < 0,
* Apply Rule A (Symbol Convention) and Rule B (Definition of Paraxial Region).
(vi) When f < 0, z > 0 ,
* Apply Rule A (Symbol Convention) and Rule B (Definition of Paraxial Region).
In the coding convention of Rule A, (i) to (vi) are all forms of the same model, and if the magnification is β, then from the relationship of similarity of triangles, the following can be said in common for all of (i) to (vi):
Next, let's consider the following major rays.
* Apply Rule A (Symbol Convention) and Rule B (Definition of Paraxial Region).
From the figure, the following can be said:
From (1)-①', (1)-②a, and (1)-②b,
From Lagrange invariant,
From (1)-④, (1)-②abc, and (1)-①',
From (1)-① and (1)-⑤,
Next, if we define the longitudinal magnification α as follows, then from (1)-①, (1)-①', and (1)-⑥,
It should be noted that this only holds true when z is a small quantity.
*When the medium is air:
Since n=n'=1, the following can be said.
From (1)-①, (1)-①', and (1)-⑧, we obtain the following equations, including Newton's equations (or Newton's relations).
Next, by applying (1)-(3) to (1)-(8), the following Gaussian lens formula is obtained.
Next, the object-to-image distance will be described.
If the object-image distance is defined as a one-dimensional vector with the object point as the start point and the image point as the end point, it can be expressed as follows:
Next, for the vertical magnification, by substituting n=n'=1 into (1)-⑦, the following formula is obtained.
・Light propagation based on Gaussian imaging with an ideal lens
1. Relationship between the conversion of incident and emitted ray by optical elements.
Consider a system in which a ray of light passes sequentially from the 0th plane (object plane), the 1st plane, the 2nd plane, ..., the j-th plane, ..., to the nth plane (image plane). In this case, the parameters of the j-th plane are illustrated as follows.
* Apply Rule A (Signing Convention), Rule B (Definition of Paraxial Region), and Rule C (Color Distinction of Parameter).
• Distance from the front focal point to the front principal point of the optical element: fj
Distance from the rear principal point of the optical element to the rear focal point: fj'
• Distance from the front principal point of the optical element to the object-side aerial image: sj
Distance from the rear principal point of the optical element to the image-side aerial image: s'j
• Height of the incident light ray from the object side at the optical element position: hj
• Height of the light rays emitted towards the image side at the optical element position: h'j
- Angle of light rays on the object: uj
• Ray angle on the image side: u'j
Use air as the medium and apply (1)-⑧.
Here, from the diagram,
From (1)-⑩ (Gauss's imaging formula),
From (2)-①a,b,c, (2)-②,
(2)-①a and (2)-②' can be expressed as matrices as follows.
The height of the ray remains unchanged before and after the front and rear principal points; only the angle changes due to refraction. The matrix representing this refraction relationship is denoted by Sj and is called the interface matrix.
2. Theoretical description of light ray propagation between optical elements
Next, we will describe the propagation of light rays between the j-th and (j+1)-th planes. The parameters for this propagation are shown in the following diagram.
* Apply Rule A (Signing Convention), Rule B (Definition of Paraxial Region), and Rule C (Color Distinction of Parameter).
• Distance between the principal points of the optical element: HH'j
• Distance from the rear principal point of one optical element to the front principal point of the next optical element: dj
• Distance from the origin (first object plane) to the front principal point of the optical element: Lj
At this point, the following can be said from the diagram.
(2)-④b and c can be expressed as matrices as follows:
A ray passing through the rear principal point of the j-th plane propagates while maintaining its angle until it reaches the front principal point of the next (j+1)-th plane, with only its height changing. The matrix representing this propagation process is denoted by Tj and is called the spatial propagation matrix.
3. Matrix representation of ray propagation
The propagation of light rays is tracked using the height and angle of the ray as state variables and described in matrix form. This allows for a unified and concise treatment of the behavior of the optical system.
The product of the interface matrix ((2)-③) and the spatial propagation matrix ((2)-⑤) is defined as follows.
Applying this, (2)-③ and (2)-⑤ can be expressed as follows.
By using the recurrence relation (2)-⑤', we can trace the propagation of light rays from the first surface to the final surface.
Furthermore, the inverse matrices of (2)-③a, (2)-⑤a, and (2)-⑥a are as follows:
Applying this, reverse ray propagation can be represented as follows:
From the above, we have shown that the propagation of light rays in both the forward and reverse directions can be described uniformly using matrix form. Based on these representations, the next chapter will describe a method for evaluating the characteristics of the entire optical system.
・Derivation of optical system parameters in light ray propagation using an ideal lens
For an optical system consisting of k ideal lenses, the optical system parameters are determined using the tracing results of each ray shown in the figure below.
* Apply Rule A (Signing Convention), Rule B (Definition of Paraxial Region), and Rule C (Color Distinction of Parameter).
In the paraxial region, "angle" and "ray height" are typically treated as infinitesimal quantities. However, in the calculation process of optical system parameters, these do not necessarily need to be limited to infinitesimal quantities and can be extended to any scale.
Therefore, in on-axis paraxial tracking, the ray height at the principal point is normalized to 1, and in principal ray paraxial tracking, the angle between the object height and the principal point is normalized to -1, enabling a simple and unified analysis.
Here, we will distinguish each ray by changing the subscript of each parameter (the part marked with "*" below).
-
From (2)-③,⑤', (2)-③',⑤'',
-
From (2)-①a,b,c,
Parameters without subscripts are considered system-specific parameters that characterize the entire optical system.
Furthermore, the total lens length TL in the figure can be expressed as follows using Lj defined in (2)-④a.
(i) Parallel light incident from the front
Here, we determine the combined focal length f, the distance of the rear principal point HD', and the image length BFinf for parallel light incidence.
* Apply Rule A (Signing Convention), Rule B (Definition of Paraxial Region), and Rule C (Color Distinction of Parameter).
From (3)-①a,b,
Since uF',1 = 0 and hF',1 = 1,
From (3)-②a, the following image distance for parallel light incidence is obtained.
Therefore, by applying (3)-④a' and (3)-④b to the relationship obtained from the figure, we can obtain the focal length f and the rear principal point distance HD'.
(ii) Parallel light incident from the rear (reverse direction)
Here, we determine the combined focal length f, the front principal point distance HD, and the object distance FFinf for parallel light incidence.
* Apply Rule A (Signing Convention), Rule B (Definition of Paraxial Region), and Rule C (Color Distinction of Parameter).
From (3)-①a',b',
From u'F,k = 0 and h'F,k = -1,
From (3)-②a, the object distance for parallel light incidence is obtained as follows.
Therefore, by applying (3)-⑤a' and (3)-⑤b to the relationship obtained from the figure, we can obtain the focal length f and the front principal point distance HD.
Next, we calculate the distance HH' between the principal points as defined in the figure below.
* Apply Rule A (Signing Convention), Rule B (Definition of Paraxial Region), and Rule C (Color Distinction of Parameter).
From the figure, the distance between principal points HH' can be calculated as follows using TL, HD, and HD' obtained so far.
(iii) On-axis paraxial tracking
Here, the image distance BF and magnification β are determined using the object distance FF and the front principal point distance HD.
* Apply Rule A (Signing Convention), Rule B (Definition of Paraxial Region), and Rule C (Color Distinction of Parameter).
From the diagram,
Applying these, from (3)-①a,b,
Therefore, by applying (3)-⑥a,b,c to the following equation obtained from (3)-②a,b, we can obtain the image distance BF and magnification β.
(iv) Primal ray paraxial stop
Here, the exit pupil position XP is determined using the object distance FF, the front principal point position HD, and the entrance pupil position EP.
* Apply Rule A (Signing Convention), Rule B (Definition of Paraxial Region), and Rule C (Color Distinction of Parameter).
From the diagram,
Applying these, from (3)-①a,b,
Therefore, by applying (3)-⑦a,b,c to the following equation obtained from (3)-②a, the exit pupil position XP can be obtained.
Furthermore, although the magnification β was found in (3)-⑥e, another method for deriving that magnification β is shown below.
First, the following can be said from the diagram.
Therefore, the magnification factor β can be calculated as follows.
-
Calculation format
As an example of implementing the optical system parameter derivation method described above, a corresponding Excel format is provided as an appendix. This format is designed to allow systematic calculation of various parameters (focal length, principal point position, back focus, etc.) based on matrix optics, and can be used for both theoretical studies and practical analyses.
This format also includes the calculations for pupil magnification and effective F-number derived on a separate page (FN) .
・Specific example: Derivation of the combined focal length
This section derives the combined focal length when two optical elements with focal lengths f1 and f2 are placed at a distance d apart. When this system is represented using a paraxial ray propagation model with an ideal lens, it takes the configuration shown in the figure below.
* Apply Rule A (Signing Convention), Rule B (Definition of Paraxial Region), and Rule C (Color Distinction of Parameter).
Applying the parameter values ( d1 = d ) shown in the figure to (3)-4a', we obtain the following:
Calculating (4)-① yields the following result.
Therefore, applying (4)-①' to the following equation obtained from (3)-④b and c, we obtain the following relationship between the image distance BFinf and the combined focal length f.
-
Related Topics
For wave-optical propagation models, including diffraction, such as Gaussian beams, please refer to the following link.
We have discussed ray propagation in an ideal lens system, but for its relationship to ray propagation in an actual spherical lens, please refer to the paraxial calculation theory shown below.
Based on the optical system presented here, the pupil magnification and effective F-number are derived as follows.
-
References
[1] Kishikawa, Toshiro. Yūzā Enjinia no Tame no Kōgaku Nyūmon (An Introduction to Optics for User Engineers). Optronics Co.
-
Update History
-
2026-07: Comprehensive revision
-
2026-01: Comprehensive revision
-
2025-10: Comprehensive revision
-
2025-06: Newly released