Next Article in Journal
Poisson Bracket Filter for the Effective Lagrangians
Previous Article in Journal
The Zakharov–Shabat Spectral Problem for Complexification and Perturbation of the Korteweg–de Vries Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geometric Properties of Planar and Spherical Interception Curves

School of IT and Engineering, ADA University, Ahmadbey Aghaoglu Str. 61, Baku AZ1008, Azerbaijan
Axioms 2023, 12(7), 704; https://doi.org/10.3390/axioms12070704
Submission received: 15 June 2023 / Revised: 11 July 2023 / Accepted: 13 July 2023 / Published: 20 July 2023
(This article belongs to the Special Issue Advances in Mathematics and Its Applications II)

Abstract

:
In this paper, some geometric properties of the plane interception curve defined by a nonlinear ordinary differential equation are discussed. Its parametric representation is used to find the limits of some triangle elements associated with the curve. These limits have some connections with the lemniscate constants A , B and Gauss’s constant G, which are used to compare with the classical pursuit curve. The analogous spherical geometry problem is solved using a spherical curve defined by the Gudermannian function. It is shown that the results agree with the angle-preserving property of Mercator and Stereographic projections. The Mercator and Stereographic projections also reveal the symmetry of this curve with respect to Spherical and Logarithmic Spirals. The geometric properties of the spherical curve are proved in two ways, analytically and using a lemma about spherical angles. A similar lemma for the planar case is also mentioned. The paper shows symmetry/asymmetry between the spherical and planar cases and the derivation of properties of these curves as limiting cases of some plane and spherical geometry results.

1. Introduction

The study of curves and their geometric properties is one of the oldest topics in mathematics. There are many interesting curves on a plane or sphere that are defined kinematically or geometrically and can be studied using the methods of calculus, differential equations, and differential geometry. For example, on a plane, one can mention tractrix, pursuit curve, lemniscate, brachistochrone, cycloid, conchoid of Nicomedes, limaçon of Pascal, etc., and on sphere great circle, spherical spiral, which is a special case of loxodrome also known as rhumb line, etc. [1]. The study of curves played a crucial role in the development of mathematics and its applications. Various attempts to describe the curves led to the discovery of different coordinate systems, such as cartesian, polar, barycentric, etc. The study of geometric properties such as lengths, areas, angles, intersections, etc., led to the development of branches of mathematics such as analytic geometry, calculus, differential geometry, differential equations, topology, etc. In antiquity, the curves were used for the solution of construction problems such as trisection of an angle, squaring the circle, and doubling the cube. In the middle ages and renaissance, curves were used for the solution of various algebraic and geometric problems, for the description of trajectories of projectiles and planets, etc. See [2,3,4] for the history of curves. In modern times, many problems of science, technology, engineering, art, and mathematics require the study of various curves on different spaces and their geometric properties. Applications of curves include, but are not bound to, aeronautics, architecture, design, sea and air navigation, music, etc. With the dawn of computer graphics, it became easy to create and represent curves and use them in many aspects of life.
In Section 2 of the current paper, the curve, which is named the interception curve, was studied. This curve can be obtained if one point moves uniformly along a straight line, and another point, initially one unit apart from both the line and the first point on this line, moves with the same speed so that it always stays on the line passing through the first point and the initial position of the second point. This plane curve appears in problems related to the interception of high-speed targets by beam rider missiles (hence the name interception curve) [5,6]. The curve was also studied in [7,8,9]. In [10], in Section 1.460 (polar coordinates) and Section 1.507 (cartesian coordinates), some methods were proposed to determine a formula for this curve. In both cases, the obtained results were probably complicated enough to discourage further study of this curve. In the present paper, the equation in cartesian coordinates was used to prove some geometric properties of the curve. The cartesian equation was studied further to find a simple parametrization among other possible parametrizations. This parametrization helped us to prove some more geometric characteristics, which show its connections with Lemniscate constants and Gauss’s constant [11,12]. The use of mathematical software to perform experiments and computations and to draw accurate and colorful graphs was very important for the results of the current paper. For example, the connection between the limit of | P Q | below and the lemniscate constants was suggested by Wolfram Alpha. Most of the graphs and pictures in the current paper were created using the website https://www.geogebra.org/ (accessed on 17 February 2023).
Some of the curves appear naturally in pursuit-evasion problems. The interception curve can also be studied in this context. The comparison of the well-known pursuit curve with the interception curve is performed at the end of Section 2.
In Section 3, a similar question for a sphere was asked. In the case of the sphere, it was shown that the curve defined by the Gudermannian function has similar geometric properties, and these properties are proved first directly using analytic methods. At the end of Section 3, the connection of these characteristics with Mercator and stereographic projections was discussed. These conformal projections also revealed the symmetry of this curve in connection with spherical and logarithmic spirals.
In Section 4, it was shown how to use a lemma from spherical geometry to prove the mentioned geometric properties of the spherical curve indirectly. The planar case also accepts simple geometric methods, and these were mentioned at the end of Section 4. There is a remarkable symmetry/asymmetry between the planar and spherical cases (see Appendix A, Table A1), and going back and forth between these two cases was very helpful in writing more of these geometric characteristics.

2. Planar Curve

Let us start with discussing the following question and setting the notation that is used throughout the paper.
Question 1. Suppose that two points, P ( x , y ) and Q, initially at O ( 0 , 0 ) and A ( 1 , 0 ) , respectively, move with constant and equal velocities so that Q is on the line x = 1 and P is on the ray O Q . What curve is defined by point P?
Similar questions were discussed in [7,9]. Polar coordinates can be used to describe the resulting curve. Let us denote r = | O P | and A O Q = θ .
Since the speeds of points P and Q are equal, the length of curve O P and the length of line segment A Q , which is tan θ , are equal for each θ . By using the well-known formula for the length of a curve r = r ( θ ) given in polar coordinates, we find that
0 θ r ( t ) 2 + ( r ( t ) ) 2 d t = tan θ .
By taking the derivative of both sides of Equation (1) and simplifying, we obtain an ODE,
r ( θ ) 2 + ( r ( θ ) ) 2 = 1 cos 4 θ ,
with initial condition r ( 0 ) = 0 . This nonlinear equation appears in problems related to the interception of high-speed targets by beam rider missiles [5]. Therefore, in this paper, the curve defined by Equation (2) is called the interception curve. See [13] for a similar problem with r ( 0 ) = 1 .
One can find approximate solutions of Equation (2) using power series. If we substitute θ = 0 in Equation (2), we find that r ( 0 ) = ± 1 . By taking the derivative of Equation (2) and substituting θ = 0 , we can find r ( 0 ) = 0 . By continuing this process, we can find r ( 0 ) = ± 1 , etc. Therefore, the two solutions are r 1 ( θ ) = θ + 1 6 θ 3 + 7 40 θ 5 + 43 720 θ 7 + O ( θ 9 ) and r 2 ( θ ) = r 1 ( θ ) . Maple 2022 software was used to calculate these four nonzero terms of the series, which were then used in GeoGebra to draw Figure 1.
Equation (2) can be solved by quadratures [8] (see also [13]). For example, one can follow the steps given in [10] to find r 1 ( θ ) .
  • Use substitution r = r cot u to obtain r cos 2 θ = ± a sin u . By differentiating and excluding r and r , we obtain u + 2 tan u tan θ = 1 ([10], Section 1.460);
  • Use substitutions η ( ξ ) = tan u , ξ = tan θ to obtain Abel’s equation ( ξ 2 + 1 ) η = ( η 2 + 1 ) ( 1 2 ξ η ) ([10], Section 1.81);
  • Use substitution ξ 4 η ( ξ ) = ( ξ 2 + 1 ) z + ξ 3 to obtain one more Abel’s equation ξ 7 z + 2 ( ξ 2 + 1 ) z 3 + 5 ξ 3 z 2 = 0 ([10], Section 1.151);
  • Use substitution v = 1 z to obtain ξ 7 v v = 2 ( ξ 2 + 1 ) + 5 ξ 3 v ([10], Section 1.185);
  • Use substitution ξ w = ξ 3 v + 1 to obtain linear equation d ξ d w ξ w 2 ( w 2 + 1 ) + 1 2 ( w 2 + 1 ) = 0 , which can be solved by quadratures ([10], Section 1.185).
As one can understand from these steps, the final result will not be a simple expression. In particular, by using substitution w = tan s , the solution of the linear differential equation in the last step can be reduced to the integral sec s d s = 2 F ( x 2 | 2 ) + const , where F ( x | m ) is the elliptic integral of the first kind with parameter m = k 2 (cf. [6]).
Note also that r = r 1 ( θ ) and x = 1 do not intersect for 0 < θ < π / 2 , because otherwise, if the points P and Q coincided, the length of the curve O P would be equal to A Q , which is impossible.
Answer to Question 1. Now, let us try using cartesian coordinates to find the solution of Equation (2) as a parametric curve. First, note that in the cartesian coordinates, Equation (1) can be written as
0 x 1 + ( y ( t ) ) 2 d t = y x .
By taking the derivative of both sides of Equation (3) and simplifying, we obtain
x 2 1 + ( y ( x ) ) 2 = y x y ,
which in particular agrees with y ( 0 ) = 0 . A parametrization of the solution can be written by solving Equation (4) for y to obtain y = y x x 2 ( y ) 2 + 1 = p x x 2 p 2 + 1 , where y = p 0 . By noting that d y = p d x and d y = p d x + x d p 2 x p 2 + 1 d x x 2 p p 2 + 1 d p , a linear differential equation d x d p x p 2 p 2 + 1 = 1 2 p 2 + 1 is obtained. By solving this equation, we obtain the parametrization (cf. [10], Section 1.507, where the roles of x and y are interchanged):
x ( p ) = 1 p 2 + 1 4 0 p d t 2 t 2 + 1 4 , y ( p ) = p x ( p ) ( x ( p ) ) 2 p 2 + 1 ( p 0 ) .
Of course, there are also other possible parametrizations of this curve (see Appendix B). This curve also has some interesting geometric properties. Some of these properties (Theorems 1 and 2) can be taken as its definition.
Theorem 1. 
Suppose that the tangent line of the curve Equation (4) at the point P intersects x-axis, y-axis, and the line x = 1 at points F, U, and T, respectively (see Figure 1). Then,
1 
| U P | = | O U | + | T Q | ,
2 
( 1 x ) · | U P | = | T Q | ,
3 
x · | P T | = | T Q | ,
4 
sin Q P T = | O P | · sin 2 T Q P ,
where x is the abscissa of point P ( x , y ) .
Proof. 
The equation of line O Q is Y y X x = y x . By substituting X = 1 in this equation, we can find the coordinates of point Q 1 , y x . The equation of tangent line U T is Y y X x = y . By substituting Y = 0 , X = 0 , and X = 1 in this equation, we can find the coordinates of points F, U, and T, respectively:
F x y y , 0 , U 0 , y x y , T 1 , y + ( 1 x ) y .
We can now find
| U P | = x 1 + ( y ) 2 , | O U | = y + x y ,
| T Q | = y + y x y y x = ( 1 x ) y y x , | P T | = ( 1 x ) 1 + ( y ) 2 .
It remains only to put these formulas in the required equalities to check that each of the first three equalities is equivalent to Equation (4). The fourth part follows from the previous equalities and the sine theorem for P Q T . □
It is also possible to prove these geometric properties using plane geometry methods, which we discuss in Section 4. Using these plane geometry methods, one can also prove the following geometric property.
Theorem 2. 
The radius of the circle through O and tangent to line U T at point P is equal to the radius of the circle through O and tangent to line A T at point Q (see Figure 1).
Proof. 
The center of the circle through O and tangent to line U T at point P is the intersection of perpendicular bisector Y y 2 = x y X x 2 of line segment O P and line Y y = 1 y X x perpendicular to line P T at point P. We find that X = x ( x 2 + y 2 ) y 2 ( x y y ) and Y = y + x 2 + y 2 2 ( x y y ) . Therefore, the radius of this circle is
X x 2 + Y y 2 = x 2 + y 2 2 ( x y y ) 1 + ( y ) 2 .
Similarly, the center of the circle through O and tangent to the line A T at point Q is the intersection of perpendicular bisector Y y 2 x = x y X 1 2 of line segment O Q , and line Y = y x perpendicular to the line A T at point Q. Using this, we find that the radius of this circle is 1 X = x 2 + y 2 2 x 2 . It remains only to check that by Equation (4), these radii are equal. □
The following result shows that there is a connection between the interception curve and the lemniscate functions.
Theorem 3. 
The length of side P Q and the difference of lengths of the other two sides of P Q T approach the same limit B 2 as x 1 :
lim x 1 | P Q | = lim x 1 | P T | | T Q | = Γ 3 4 4 2 π = B 2 0.3588850048 ,
where B is the second lemniscate constant.
Proof. 
First, we calculate that | P Q | = ( 1 x ) 2 + ( 1 x ) y x 2 . Since x ( p ) 1 as p + , lim x 1 | P Q | = lim p + ( 1 x ) y x . On the other hand, using parametrization Equation (5), we can write
( 1 x ( p ) ) y ( p ) x ( p ) = ( 1 x ( p ) ) 2 p 2 + 1 1 x ( p ) p 2 + 1 + p .
The subtrahend 1 x ( p ) p 2 + 1 + p in the last expression approches 0 as p + . Therefore, it remains only to find
lim p + ( 1 x ( p ) ) 2 p 2 + 1 = lim p + p 2 + 1 4 0 p d t 2 t 2 + 1 4 2 .
Since p 2 + 1 4 = 1 + 0 p t d t 2 ( t 2 + 1 ) 3 4 ,
lim p + ( 1 x ( p ) ) 2 p 2 + 1 = 1 0 + t 2 + 1 t d t 2 ( t 2 + 1 ) 3 4 2 .
By substituting t = tan s , where 0 s < π 2 , in the improper integral, after simplifications we obtain
0 + t 2 + 1 t d t 2 ( t 2 + 1 ) 3 4 = 0 π 2 cos s d s 2 ( 1 + sin s ) .
It is known that [11,12,14,15]
0 π 2 cos s d s 2 = 0 1 x 2 d x 1 x 4 = Γ 3 4 2 2 π = B .
We only need to show that
0 π 2 cos s d s 2 ( 1 + sin s ) = 1 B ,
which follows from
0 π 2 cos s d s 2 ( 1 + sin s ) + 0 π 2 cos s d s 2 =
0 π 2 cos s 2 1 1 + sin s + 1 d s = cos s 2 4 sin s 2 sin s 2 + cos s 2 2 | 0 π 2 = 1 .
Similarly,
| P T | | T Q | = ( 1 x ) 1 + ( y ) 2 y + y x .
Therefore,
lim x 1 | P T | | T Q | = lim p + ( 1 x ( p ) ) p 2 + 1 p + y ( p ) x ( p ) =
= lim p + 1 x ( p ) p 2 + 1 + p + lim p + ( 1 x ( p ) ) y ( p ) x ( p ) = 0 + B 2 = B 2 .
Note that B 2 = 1 4 G 2 , where G is Gauss’s constant defined by the arithmetic–geometric mean [11,12,14,15].
There are other strategies for making point P follow point Q so that lim x 1 | P Q | is smaller. However, it is possible to show that there is no optimal strategy. Consider, for example, the following scenario. Suppose that there is a barrier along line x = 1 and an evader is at point ( 1 , 0 ) . A pursuer at the origin starts to chase with a constant speed, and the evader moves with the same speed along line x = 1 . It would be natural to think that the pursuer should always move toward the current position of the evader. This gives rise to a trajectory for the pursuer, commonly known as pursuit curve. However, the pursuer can also move towards a future position of the evader instead of the evader’s current position. One of the many possible ways for the pursuer to do this is to stay on the line connecting the origin with the evader. The optimal strategy for the evader is to run along barrier x = 1 by never changing its direction. However, there is no optimal trajectory for the pursuer, so any candidate for the optimal curve can be improved by replacing one of its curved parts with a line segment and, therefore, shortening that part. Because of the presence of barrier x = 1 , the curved part should always exist. Note that any straight line through the origin intersects x = 1 , except line x = 0 , which obviously can not be the optimal curve. One can pose a similar question for a sphere, where the barrier and the evader are on the great circle of the sphere, and the pursuer is at one of the poles of the great circle. This is performed in the following section.
Consider also this trajectory: first, P moves with a constant speed from O to the point with coordinates ( 1 , k ) along line y = k x where k 0 , 0 x 1 , while Q moves with the same speed from A to the point with coordinates ( 1 , k 2 + 1 ) , after which P follows Q along line x = 1 and distance | P Q | = k 2 + 1 k does not change. The value of k 2 + 1 k = 1 k 2 + 1 + k can be made arbitrarily close to zero. On the other hand, for any chosen curve, lim x 1 | P Q | = 0 is not possible. Indeed, | P Q | | ( P Q ) y | , where ( P Q ) y is the projection of P Q onto the y-axis, and | ( P Q ) y | is a non-decreasing function of time (or x). Therefore, for any given curve with a particular positive value of lim x 1 | P Q | , there is always a more efficient curve making the limit even smaller. In view of these observations, it would be interesting to compare the curve in Figure 1 with the mentioned trajectory of P chasing Q, the pursuit curve.
Figure 2 was created using Maple 2022. The Maple 2022 software was also used to simplify the limits and the integrals in the paper.
Suppose now that P starts to move from O with a constant speed along such a curve that the tangent line of the curve at P always passes through Q, which moves as before with the same constant speed from A upwards along line x = 1 . The obtained curve is known as the pursuit curve and defined by differential equation 1 + ( y ) 2 = ( 1 x ) y , with initial conditions y ( 0 ) = y ( 0 ) = 0 . See, for example, [16] (page 256), [3,4] (page 241 of the 1911 edition, page 607 of the 1902 edition) for the solution found by P. Bouguer in 1732:
y = 1 4 ( 1 x ) 2 1 2 ln ( 1 x ) 1 4 .
One can check that for this curve, | P Q | = 1 2 + ( 1 x ) 2 2 , and therefore lim x 1 | P Q | = 1 2 > B 2 , which means that the interception curve overperforms the pursuit curve in this regard.

3. Spherical Curve

We now study analogous questions for spherical geometry.
Question 2. Suppose that two points P and Q, initially at B ( 0 , 0 , 1 ) and A ( 1 , 0 , 0 ) , respectively, move with constant and equal velocity so that Q is on the great circle z = 0 , x 2 + y 2 = 1 of sphere x 2 + y 2 + z 2 = 1 with center O ( 0 , 0 , 0 ) , and P is on the great circle through B and Q of the sphere. What curve is defined by point P? (see Figure 3).
Answer. We use spherical coordinates to describe the resulting interception curve. We denote A O Q = θ and P O B = ϕ . Since ρ = | O P | = 1 , for the coordinates of point P ( x , y , z ) , we can write x = cos θ sin ϕ , y = sin θ sin ϕ , and z = cos ϕ , where we assume that ϕ = ϕ ( θ ) is a function of θ . Therefore, x θ = sin θ sin ϕ + ϕ cos θ cos ϕ , y θ = cos θ sin ϕ + ϕ sin θ cos ϕ , and z θ = ϕ sin ϕ . Since the points P and Q travel equal distances, by using the well-known formula from vector calculus for the length of a parametrically defined curve, we can write
0 θ ( x ( t ) ) 2 + ( y ( t ) ) 2 + ( z ( t ) ) 2 d t = θ ,
where θ on the right-hand side of the last equality is the distance traveled by Q. By taking the derivative of both sides of Equation (6) with respect to θ and simplifying using the formula for x θ , y θ , and z θ , we obtain sin 2 ϕ + ( ϕ ( θ ) ) 2 = 1 with initial condition ϕ ( 0 ) = 0 . By solving this ODE, we obtain
ϕ = 2 tan 1 e θ π 2 .
Note that Equation (7), which can also be written as ϕ = tan 1 sinh θ , is sometimes called the Gudermannian function gd ( x ) . See, for example, [17], Section 4.23(viii) and [18], Section 6.12. The Gudermannian function gd ( x ) is the vertical component of the Mercator projection [19], and we see below that this connection with the Mercator projection is not a coincidence.
Theorem 4. 
The points P and Q are approaching each other and lim θ | P Q | = 0 .
Proof. 
One can observe that
lim θ ϕ ( θ ) = lim θ 2 tan 1 e θ π 2 = 2 · π 2 π 2 = π 2 ,
and therefore the distance between the points P ( 1 , θ , ϕ ) and Q ( 1 , θ , π 2 ) approaches zero as the curve winds around the sphere. □
This curve has other geometric properties. In the following, X Y ^ means the spherical distance between points X and Y on a sphere. For a unit sphere with center O, X Y ^ = X O Y .
Theorem 5. 
If a great circle tangent to the curve ϕ = 2 tan 1 e θ π 2 at point P of unit sphere intersects the great circle on the plane x O y at point T (see Figure 3), then
1 
the sum of the lengths of the arcs P T and T Q is not dependent on θ, and P T ^ + T Q ^ = π 2 ,
2 
as θ increases, T Q ^ increases, P T ^ decreases, and both approach to π 4 , as θ .
3 
spherical angle Q P T is equal to B P ^ ,
4 
spherical angle B P T is equal to P Q ^ + π 2 .
Proof. 
We continue to use the notations introduced in Question 2 and its answer. Since ϕ ( θ ) = cos ϕ , we can write the following for tangent vector v = ( x θ , y θ , z θ ) of curve: Equation (7)
v = ( sin θ sin ϕ + cos θ cos 2 ϕ , cos θ sin ϕ + sin θ cos 2 ϕ , sin ϕ cos ϕ ) .
We can find the normal vector of the plane containing the great circle through P and T, as n 1 = v × r , where r = ( x , y , z ) :
n 1 = ( cos θ sin ϕ cos ϕ + sin θ cos ϕ , sin θ sin ϕ cos ϕ cos θ cos ϕ , sin 2 ϕ ) .
Then, the equation of this plane is
( cos θ sin ϕ cos ϕ + sin θ cos ϕ ) x + ( sin θ sin ϕ cos ϕ cos θ cos ϕ ) y ( sin 2 ϕ ) z = 0 .
By putting z = 0 in this equation, we can find the only unknown coordinate Θ of point T ( 1 , Θ , π 2 ) in spherical coordinates:
tan Θ = y x = cos θ sin ϕ cos ϕ + sin θ cos ϕ sin θ sin ϕ cos ϕ + cos θ cos ϕ = tan θ + sin ϕ 1 tan θ sin ϕ .
By Equation (7), sin ϕ = tanh θ . Therefore,
tan Θ = tan ( θ + tan 1 tanh θ ) ,
which means that Θ = θ + tan 1 tanh θ and consequently Q T ^ = tan 1 tanh θ . Using Spherical Pythagorean Theorem ([20], page 112) for P Q T , we obtain
cos P T ^ = cos P Q ^ cos Q T ^ = sin ϕ cos tan 1 tanh θ = tanh θ cosh θ cosh 2 θ = sinh θ cosh 2 θ .
Therefore, tan P T ^ = coth θ , and, consequently, tan P T ^ tan Q T ^ = coth θ tanh θ = 1 . Since P T ^ , T Q ^ < π 2 , this proves that P T ^ + T Q ^ = π 2 . Using formula P T ^ = tan 1 coth θ and Q T ^ = tan 1 tanh θ , we can now show that as θ increases, T Q ^ increases, but P T ^ decreases, and both approach π 4 . Using the formula for a right angled spherical triange ([20], page 127, 4.35), we obtain
tan Q P T = tan Q T ^ sin P T ^ = tanh θ sin π 2 ϕ = sin ϕ cos ϕ = tan ϕ .
Therefore, Q P T = B P ^ . Similarly, B P T = P Q ^ + π 2 . □
Using similar analytic methods, one can prove the following property of curve Equation (7) for tangent small circles. It is also possible to prove this theorem using Lemma 2.
Theorem 6. 
The spherical radius R of a small circle of the unit sphere, which passes through point B and is tangent to curve Equation (7) at point P (see Figure 3), satisfies tan 2 R = 1 4 sec 4 1 2 B P ^ .
We now discuss the mentioned connection with the Mercator projection. Let us perform a projection of the sphere onto cylinder x 2 + y 2 = 1 using formula x = θ and y = ln tan π 4 + φ 2 (see [21], p. 191). Note that here, x and y are new coordinates over the cylinder [19]. Note also the difference between the spherical coordinate ϕ and latitude φ = π 2 ϕ . Using Equation (7) and the last two equalities, we can write
y = ln π 4 + π 2 2 tan 1 e θ π 2 2 = ln coth θ 2 .
Since x = θ on the cylinder, we determine that the Mercator projection of curve Equation (7) is y = ln coth x 2 , which is an inverse function of itself. This means that after the application of the Mercator projection, curve Equation (7) has a symmetry with respect to line y = x (see Figure 4). Note that y = x is a Helix on the cylinder, and it is the projection of a Spherical Spiral (a special case of loxodrom or rhumb line [1]) intersecting the meridians ϕ = const of the sphere at a constant angle π 4 . Since the Mercator projection is conformal (angle-preserving), the slope y = 1 sinh x on the cylinder and cot B P T on the sphere should be equal. Indeed, by Theorem 5,
cot B P T = cot π ϕ = cot ϕ
= cot 2 tan 1 e θ π 2 = tan 2 tan 1 e θ = 2 e θ 1 e 2 θ = 1 sinh θ .
Note that in the literature, the expression spherical spiral is sometimes used for different curves. In Figure 4, we used parametrization x = cos θ cosh θ , y = sin θ cosh θ , z = tanh θ for this curve (red curve) ([22], Lemma 8.27), the parametrization x = cos θ , y = sin θ , z = θ for y = x (orange curve), which is the Mercator projection of the previous curve onto the cylinder, parametrization x = cos θ , y = sin θ , z = ln coth θ 2 for y = ln coth x 2 (green curve), which is the Mercator projection of interception curve Equation (7) (blue curve). Note also that in Figure 4, the stereographic projection of Equation (7) onto the x y plane, with respect to pole B, is a plane spiral curve (black curve), which in polar coordinates can be written as r = coth θ 2 . Indeed, by Equation (7),
r = 1 tan ϕ 2 = 1 tan tan 1 e θ π 4 = coth θ 2 .
Line y = 2 (yellow line) on the x y plane is the asymptote of this curve. It is known that the stereographic projection of the rhumb line is a logarithmic spiral r = e θ (purple curve) ([22], Lemma 8.27). In particular, since curves y = x and y = ln coth x 2 intersect at a right angle, and stereographic and Mercator projections are both conformal, the curves on the sphere, on the cylinder, and the plane in Figure 4 intersect at right angles at the indicated points. See Appendix A, Table A2 for the list of corresponding curves in Figure 4.

4. Alternative Methods for Spherical and Planar Curves

We can also prove Theorem 5 using the following lemma, which is interesting on its own.
Lemma 1. 
On a unit sphere with center O, a great circle and one of its poles B is drawn. Two great circles through B intersect the first great circle at Q 1 and Q 2 , so that Q 1 Q 2 ^ < π 2 . Another great circle intersects arcs B Q 1 and B Q 2 at points P 1 and P 2 so that P 1 P 2 ^ = Q 1 Q 2 ^ . This great circle also intersects the first great circle at T (See Figure 5). Then,
P 1 T ^ + T Q 2 ^ = P 2 T ^ + T Q 1 ^ = π 2 .
Proof. 
Let us denote P 1 P 2 ^ = Q 1 Q 2 ^ = x . Since the angles at Q 1 and Q 2 are right angles, by the spherical sine theorem for B Q 1 Q 2 , the angle between the two great circles at vertex B is also x ([20], page 30). P 1 Q 1 ^ = a , P 2 Q 2 ^ = b is also denoted (See Figure 6). Then, B P 1 ^ = π 2 a and B P 2 ^ = π 2 b . By the spherical sine theorem for B P 1 P 2 ([20], page 115),
sin B sin P 1 P 2 ^ = sin P 1 sin B P 2 ^ = sin P 2 sin B P 1 ^ .
Since sin B = sin P 1 P 2 ^ = sin x ,
sin B P 1 P 2 = sin B P 2 ^ = cos b , sin B P 2 P 1 sin B P 1 ^ = cos a .
We take another point P 3 on arc P 1 Q 1 such that P 2 P 3 ^ = x . Such a point, P 3 , should exist, because P 2 Q 1 ^ > x . We denote the midpoint of P 1 P 3 ^ by E and let Q 1 P 3 ^ = c . Then, P 1 E ^ = E P 3 ^ = 0.5 ( a c ) . Since spherical P 1 P 2 P 3 is isosceles, P 2 E is its angle bisector and P 2 E is perpendicular to P 1 P 3 . The three angles of the quadrilateral E Q 1 Q 2 P 2 are right angles. Therefore, the fourth E P 2 Q 2 is obtuse and E P 2 B is acute. Since P 1 P 2 B < E P 2 B , P 1 P 2 B is acute, too. Similarly, B P 3 P 2 = P 3 P 1 P 2 is acute. Consequently, B P 1 P 2 is obtuse. Consequently, B P 1 P 2 and B P 2 P 3 are obtuse, B P 2 P 1 and B P 3 P 2 are acute. Therefore, by Equation (8), B P 1 P 2 = π 2 + b and B P 2 P 1 = π 2 a . Similarly, B P 2 P 3 = π 2 + c and B P 3 P 2 = π 2 b . Therefore, P P 2 P 3 = B P 2 P 3 B P 2 P 1 = c + a and P 1 P 2 E = c + a 2 . We denote the intersection of the great circles E P 2 and Q 1 Q 2 by D. We also denote E P 2 = w , P 2 T = z , Q 2 T = y , T D = z 1 . Then, P 2 D = π 2 w and x + y + z 1 = π 2 . It remains only to show that z = z 1 . This follows from the fact that spherical P 2 T D is isosceles. Indeed, D = Q 1 E ^ = a + c 2 and T P 2 D = P 1 P 2 E = a + c 2 . Consequently, D = T P 2 D or z = z 1 . Finally,
P 1 T ^ + T Q 2 ^ = P 2 T ^ + T Q 1 ^ = x + y + z = π 2 .
For the proof of Theorem 5, one can observe that the infinitesimal line element or the length element of curve ϕ = 2 tan 1 e θ π 2 satisfies the conditions of Lemma 1 and, therefore, one just needs to pass to limit x 0 ( Q 1 approaches Q 2 ) which implies P T ^ + T Q ^ = y + z = π 2 . Since lim θ | P Q | = 0 , y and z both approach π 4 . Passing to limit x 0 also implies a b . Therefore, Q P T = lim x 0 P 1 P 2 B = lim x 0 π 2 a = π 2 b = B P ^ . Similarly, B P T = lim x 0 B P 2 T = lim x 0 π P 1 P 2 B = lim x 0 π 2 + a = π 2 + b = P Q ^ + π 2 .
We can prove Theorem 6 using the following lemma, which is a direct consequence of the formula for the radius of the circumscribed circle of a spherical triangle (see, e.g., [23], Formula (4)).
Lemma 2. 
Under conditions of Lemma 1, let a small circle of the unit sphere be circumscribed about B P 1 P 2 (See Figure 5). If the radius of the small circle is R, then
tan 2 R = 1 4 sec 4 1 2 P 1 P 2 ^ sec 4 1 2 B P 1 ^ sec 4 1 2 B P 2 ^ .
It is possible to give similar simple proofs for the equalities in the planar case, too. The equalities in Theorems 1 and 2 can be proved as a limiting case of the following lemma when Q 1 approaches Q 2 .
Lemma 3. 
Lines passing through point O and points Q 1 and Q 2 on a parallel line are drawn. A line intersects these lines at U and T, and the segments O Q 1 and O Q 2 at points P 1 and P 2 , respectively, so that | P 1 P 2 | = | Q 1 Q 2 | (see Figure 7). Then,
1 
| O U | + | T Q 2 | = | U P 1 | and | O U | + | T Q 1 | = | U P 2 | ,
2 
the radii of circles through points O , P 1 , P 2 and O , Q 1 , Q 2 are equal,
3 
if the distance between the parallel lines is 1, and the distances from points P 1 and P 2 to line O U are x 1 and x 2 , respectively, then
sin Q 1 P 1 T · sin Q 2 P 2 T = | O P 1 | · | O P 2 | · sin 2 T Q 1 P 1 · sin 2 T Q 2 P 2 ,
x 1 x 2 · | P 1 T | · | P 2 T | = | T Q 1 | · | T Q 2 | ,
( 1 x 1 ) ( 1 x 2 ) · | U P 1 | · | U P 2 | = | T Q 1 | · | T Q 2 | .
See Appendix A, Table A1 for a complete picture of corresponding lemmas.

5. Conclusions

In the paper, the two interception curves, one on a plane and the other one on a sphere, were studied. The extensive literature about the planar curve and the mentioned connections with lemnicate constants and plane geometry lemmas show that it is an important and interesting topic for mathematics and its applications. The spherical curve and its properties naturally extend these results to a geometry different from planar. The interplay between the plane and spherical geometry results and the geometric properties of curve Equations (5) and (7) were very helpful to find and prove these properties. Evident symmetry between the planar and spherical cases helped us to find more properties and shorter and more elegant proofs for them. Analogous to the planar case, there is a spherical pursuit curve studied in [24] (see also [2], page 78 of Vol. II), and it would be interesting to compare it with the spherical interception curve. It would also be interesting to consider all these questions in the context of hyperbolic geometry [24] and other planar regions [25].

Funding

This work was supported by ADA University Faculty Research and Development Fund.

Data Availability Statement

Not applicable.

Acknowledgments

I thank the reviewers and editors for their helpful suggestions and comments.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A

Table A1. Corresponding theorems in planar and spherical cases, and lemmata in Section 4 that can be used to prove these results.
Table A1. Corresponding theorems in planar and spherical cases, and lemmata in Section 4 that can be used to prove these results.
Plane CurveSpherical Curve(A)symmetryPlane GeometrySpherical Geometry
Theorem 1Theorem 5SymmetryLemma 3, 1. & 3.Lemma 1
Theorem 2Theorem 6AsymmetryLemma 3, part 2.Lemma 2
Theorem 3Theorem 4Asymmetry--
Table A2. Corresponding curves on the sphere, on the cylinder and on the plane in Figure 4.
Table A2. Corresponding curves on the sphere, on the cylinder and on the plane in Figure 4.
Sphere x 2 + y 2 + z 2 = 1 Cylinder x 2 + y 2 = 1 Mercator ProjectionPlane z = 0 Stereographic Projection
ϕ = 2 tan 1 e θ π 2
Interception curve (blue)
y = ln coth x 2
(green)
r = coth θ 2
(black)
cos θ cosh θ , sin θ cosh θ , tanh θ
Spherical spiral (red)
y = x
Helix (orange)
r = e θ
Logarithmic spiral (purple)

Appendix B

In Section 2, we found parametrization Equation (5) for the curve defined by Equation (4). It is possible to write alternative parametrizations for this curve. Let us denote z = y x and z = p 1 . Then, from Equation (4) we obtain z = p 2 1 x p . Since d z = p d x and d z = p p 2 1 d p p d x x d p , we obtain a linear differential equation d x d p + x 2 p = p 2 p 2 1 . Noting the initial condition x = 0 , y = 0 , the solution of Equation (4) can also be written as
x ( p ) = 1 p 1 p t d t 2 t 2 1 , y ( p ) = x ( p ) p 2 1 ( x ( p ) ) 2 p ( p 1 ) .
Another parametrization of the solution can be written if z = p 2 1 x p is solved first for x to obtain x = p 2 1 p z p and then, noting that d x = d z p , a linear differential equation d z d p z 2 p = 1 2 p p 2 1 is obtained. Its solution z ( p ) = p 1 p d t 2 t t t 2 1 offers us
x ( p ) = p 2 1 p 1 p 1 p d t 2 t t t 2 1 , y ( p ) = x ( p ) · p 1 p d t 2 t t t 2 1 ( p 1 ) .

References

  1. Weisstein, E.W. Loxodrome; From MathWorld—A Wolfram Web Resource. Available online: https://mathworld.wolfram.com/Loxodrome.html (accessed on 17 February 2023).
  2. Loria, G. Curve Sghembe Speciali Algebriche e Trascendenti; N. Zanichelli: Bologna, Italy, 1925; Available online: https://quod.lib.umich.edu/u/umhistmath/ABR0683.0002.001/85?rgn=full+text;view=pdf (accessed on 17 February 2023).
  3. Loria, G. Spezielle Algebraische und Transzendente Ebene Kurven: Theorie und Geschichte; B. G. Teubner: Leipzig, Germany, 1911; Volume 2, Available online: https://archive.org/details/speziellealgebr01lorigoog/page/n252/mode/2up (accessed on 17 February 2023).
  4. Loria, G. Spezielle Algebraische und Transzendente Ebene Kurven: Theorie und Geschichte. First Edition 1902. Available online: https://quod.lib.umich.edu/u/umhistmath/ABR0252.0001.001/632?rgn=full+text;view=pdf (accessed on 17 February 2023).
  5. Elnan, O.R.S.; Lo, H. Interception of High-Speed Target by Beam Rider Missile. AIAA J. 1963, 1, 1637–1639. [Google Scholar] [CrossRef]
  6. Vinh, N.X. Comment on ”Interception of High-Speed Target by Beam Rider Missile”. AIAA J. 1964, 2, 409. [Google Scholar] [CrossRef]
  7. Wilder, C.E. A discussion of a differential equation. Am. Math. Mon. 1931, 38, 17–25. [Google Scholar] [CrossRef]
  8. Zbornik, J. Akademie der Wissenschaften in Wien Mathematisch-Naturwissenschaftliche Klasse. IIa Sitzungsberichte 1957, 166, 42. [Google Scholar]
  9. Bailey, H.R. The Hiding Path. Math. Mag. 1994, 67, 40–44. [Google Scholar] [CrossRef]
  10. Kamke, E. Differentialgleichungen Lösungsmethoden und Lösungen; Springer: Wiesbaden, Germany, 1977. [Google Scholar] [CrossRef]
  11. Todd, J. The Lemniscate Constants. Commun. Acm (Cacm) 1975, 18, 14–19. [Google Scholar] [CrossRef]
  12. Todd, J. The Lemniscate Constants. In Pi: A Source Book; Springer: New York, NY, USA, 2000. [Google Scholar] [CrossRef]
  13. Aliyev, Y. From The First Remarkable Limit to a Nonlinear Differential Equation. In Proceedings of the 5th International Conference on Mathematics: An Istanbul Meeting for World Mathematicians ICOM, Istanbul, Turkey, 1–3 December 2021; pp. 293–304. Available online: http://icomath.com/dosyalar/2021%20P%20123.pdf (accessed on 17 February 2023).
  14. Cox, D.A. The Arithmetic-Geometric Mean of Gauss. L’Enseignement Mathématique 1984, 30, 275–330. [Google Scholar] [CrossRef]
  15. Cox, D.A. The Arithmetic-Geometric Mean of Gauss. In Pi: A Source Book; Springer: New York, NY, USA, 2000. [Google Scholar] [CrossRef]
  16. Savelov, A.A. Planar Curves; State Publisher of Physics-Mathematics Literature: Moscow, Russia, 1960. (In Russian) [Google Scholar]
  17. Olver, F.W.J.; Lozier, D.W.; Boisvert, R.F.; Clark, C.W. NIST Handbook of Mathematical Functions; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  18. Zwillinger, D. (Ed.) CRC Standard Mathematical Tables and Formulae; Advances in Applied Mathematics; CRC Press: Boca Raton, FL, USA, 2011. [Google Scholar]
  19. Weisstein, E.W. Mercator Projection; From MathWorld–A Wolfram Web Resource. Available online: https://mathworld.wolfram.com/MercatorProjection.html (accessed on 17 February 2023).
  20. Whittlesey, M.A. Spherical Geometry and Its Applications; Textbooks in Mathematics; CRC Press: Boca Raton, FL, USA, 2020. [Google Scholar] [CrossRef]
  21. Pearson, F. Map Projections: Theory and Applications; CRC Press: Boca Raton, FL, USA; Routledge: London, UK, 2018. [Google Scholar]
  22. Abbena, E.; Salamon, S.; Gray, A. Modern Differential Geometry of Curves and Surfaces with Mathematica; Chapman and Hall, CRC: Boca Raton, FL, USA, 2006. [Google Scholar] [CrossRef]
  23. Odani, K. On a Relation between Sine Formula and Radii of Circumcircles for Spherical Triangles. Bull. Aichi Univ. Educ. 2010, 59, 1–5. Available online: http://hdl.handle.net/10424/2951 (accessed on 17 February 2023).
  24. Roeser, E. Die Verfolgungskurve auf der Kugel; Dissert.: Halle, Germany, 1907. [Google Scholar]
  25. Morley, F.V. A Curve of Pursuit. Am. Math. Mon. 1921, 28, 54–61. Available online: https://www.jstor.org/stable/2973034 (accessed on 17 February 2023). [CrossRef]
Figure 1. Interception curve (blue) on a Plane. Its tangent line (black), the line containing its position vector (black), and the line x = 1 (blue) are also shown.
Figure 1. Interception curve (blue) on a Plane. Its tangent line (black), the line containing its position vector (black), and the line x = 1 (blue) are also shown.
Axioms 12 00704 g001
Figure 2. Comparison of interception (left, red) and pursuit (right, red) curves. Created using parametrization (5) and Bouguer’s formula. The tangent lines (green), the lines containing the position vectors (thin blue) of the curves, and lines x = 1 (thick blue) are also shown.
Figure 2. Comparison of interception (left, red) and pursuit (right, red) curves. Created using parametrization (5) and Bouguer’s formula. The tangent lines (green), the lines containing the position vectors (thin blue) of the curves, and lines x = 1 (thick blue) are also shown.
Axioms 12 00704 g002
Figure 3. Interception curve (light blue) on a Unit Sphere. Its tangent great circle (dark blue), meridian (red) and equator (black) great circles are also shown.
Figure 3. Interception curve (light blue) on a Unit Sphere. Its tangent great circle (dark blue), meridian (red) and equator (black) great circles are also shown.
Axioms 12 00704 g003
Figure 4. Mercator and stereographic projections of the spherical interception curve (blue) and spherical spiral (red). The other colors are explained in the text and in Appendix A, Table A1.
Figure 4. Mercator and stereographic projections of the spherical interception curve (blue) and spherical spiral (red). The other colors are explained in the text and in Appendix A, Table A1.
Axioms 12 00704 g004
Figure 5. Lemma 1 for spherical case. The green small circle is used in Lemma 2.
Figure 5. Lemma 1 for spherical case. The green small circle is used in Lemma 2.
Axioms 12 00704 g005
Figure 6. Proof of Lemma 1.
Figure 6. Proof of Lemma 1.
Axioms 12 00704 g006
Figure 7. Lemma 3 for planar case.
Figure 7. Lemma 3 for planar case.
Axioms 12 00704 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Aliyev, Y.N. Geometric Properties of Planar and Spherical Interception Curves. Axioms 2023, 12, 704. https://doi.org/10.3390/axioms12070704

AMA Style

Aliyev YN. Geometric Properties of Planar and Spherical Interception Curves. Axioms. 2023; 12(7):704. https://doi.org/10.3390/axioms12070704

Chicago/Turabian Style

Aliyev, Yagub N. 2023. "Geometric Properties of Planar and Spherical Interception Curves" Axioms 12, no. 7: 704. https://doi.org/10.3390/axioms12070704

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop