In this article, a new approach to solve an astronomical ship position is proposed based on a single celestial body and secant technique. The mathematical equations relating the ship position to a celestial body are described using an altitude and an azimuth of celestial body. These equations are established using knowledge of spherical trigonometry. Numerical analysis is then employed based on secant method to solve the mathematical equations for the ship position. This technique uses the altitude and the azimuth simultaneously while existing methods need to observe the altitude either of two celestial bodies or of a single body at two different times. Therefore, it is worth to apply this proposed technique when there is only a single celestial body identified. To validate the feasibility of the proposed technique, two numerical examples were performed.
The fixing of ship position at sea is an important task of navigation officers to keep the ship on her desired route. To date, global positioning system (GPS) has proven their capacity in determining the ship position precisely, quickly, and continuously. However, the signal received from GPS satellites is sometimes impossible because of solar flares, magnetic storms or other reasons, and thus the GPS receivers cannot operate in these situations. To overcome the drawbacks of GPS system in such unavailable cases, the traditional navigation method such as celestial navigation technique may be used. Intercept method has been known as most popular technique of celestial navigation to find out the ship position. To handle this method, the altitude of celestial bodies and an assumed position must be determined to establish the parameters of the line of position equation (LOP), the equations for positioning are then solved by calculation program or graphic plotting. Actually, the ship position is an intersection of circles of an equal altitude (COP) instead of LOP. Recently, with the rapid development of computer science, new computing methods have been proposed to solve the COP equations specifically for the ship position. In research of the authors [1] the trigonometric expanding technique was described to solve directly the COP equations. Later on, in [2], the vector calculus was introduced to obtain the COP equation with three unknowns in the Cartesian coordinate system, and the technique of vector expansion was then suggested to solve the ship position.
By using the evolution algorithms in artificial intelligence, studies referred in [3, 4] applied the genetic algorithm and the particle swarm optimization to determine an optimal position of ship. The advantage of these optimal algorithms is to imitate the natural processes of biological evolution in searching the global solution for complex problems without certain equations. As a continuation of the work based on the vector calculus, the authors in [5] suggested the singular value decomposition (SVD) technique integrated the least square method which is more stable than the normal method. Recently, a stand-alone celestial navigation positioning method was presented by [6]. In this study, the fixing problem was investigated for both stationary and moving observers. The advantage of this method is simplicity of the mathematical theory and applicable to portable devices. For obtaining the ship position when the horizon is invisible, the research [7] proposed an azimuth method based on the governing equations between the ship position and the great-circle azimuth of the celestial bodies. Meanwhile, the study [8] used the positioning equations through the spherical trigonometry. To improve the applicable ability of the azimuth method, a method combining the azimuth of celestial body with factors of ship route was investigated in [9].
Although the above-mentioned researches have significantly contributed to the celestial navigation for fixing, there have been still unsolved problems of fixing a position. For example, it is required to measure the altitude or azimuth of either at least two celestial bodies or a single celestial body at two different times. Besides, the existing methods of celestial navigation with more than two celestial bodies require time-consuming to conduct because many steps need to be done before fixing.
In this research, a new approach to obtain astronomical ship position is studied based on a single celestial body and secant method. The mathematical equations relating the ship position to celestial body are described based on an altitude and an azimuth of celestial body. These equations are established using knowledge of spherical trigonometry. Numerical method based on secant technique is then employed to solve the mathematical equations for ship position. The main contributions of this investigation are as follows:
First, our proposed method can perform with a single celestial body at a given time. Instead of more than two celestial bodies, this method requires only a single body. Therefore, it can be useful for some cases where there is only an identified single celestial body.
Secondly, applying secant technique to solve the system is seen as one of new approaches of advanced mathematics in celestial navigation. This technique can be a useful reference for navigators and marine students.
To verify the effectiveness of the proposed technique, two numerical examples were performed. The results showed that this technique has good performance for finding the ship position. This article consists of the parts as follows: mathematical equations for ship position are presented in Section 2. The detail of the proposed technique is shown in Section 3. Results of numerical examples are described in Section 4. The conclusion is summarized in Section 5.
Before establishing the governing mathematical equations for ship position, it is necessary to describe the parameters on a celestial sphere. Additionally, variables and symbols also are interpreted in this section. The Earth is considered as a perfect sphere. The ship position (P), a celestial body (C) and the factors of the Earth such as the meridians, the circles of latitude, and the equator are projected on a virtual sphere that is called as celestial sphere shown as Figure 1. On this celestial sphere, the ship position is defined by latitude (ϕ) and longitude (λ), and the celestial body is denoted by Greenwich hour angle (GHA) and Declination (δ). The details of variables and symbols are represented as in Table 1. From Figure 2, the relation of longitude and hour angles is presented as follows:
Generally,
LHA in
If LHA is more than 180°, then this value should be (LHA – 180°). If LHA is more than 360°, then this value should be (LHA – 360°). Otherwise, LHA is kept originally in computing process.
In this section, the governing formulas are mentioned based on the mathematical laws of spherical triangle to solve the ship position. As illustrated in Figure 3, the spherical triangle P_{N}PC is formed by three great-circle arcs: the observer meridian (P_{N}P), the celestial body meridian (P_{N}C), and the celestial body-through vertical circle. The detail of using formulas is shown as follows:
Applying the cosine formula of spherical trigonometry to the triangle P_{N}PC, we have:
Rearranging the above equation, the following equation is obtained
The four-part law is also employed for the triangle PNPC to have the below formula
Rewriting the above equation, the following form is yielded
In Section 3,
In this section, the novel technique is represented. Let the ship position P (observer) is at latitude and longitude on the surface of the Earth. At the time UTC, the celestial body (C) having the coordinate (GHA and δ) is identified and the altitude (H) and the azimuth (A) of this body are measured simultaneously from the position of observer. The ship position is an intersection of two circles which are circle of equal altitude and great circle with observed azimuth at the ship position, as shown in Figure 4, and the mathematical equations of these circles are applied to solve the ship position.
From
Substituting
Placing function (F) as follows:
The F is a function of latitude (ϕ), and the objective of the problem is how to find out the latitude (ϕ) so that the function F converge to 0. To deal with this optimization problem, the secant technique in numerical method is employed to conduct the computation procedure. Because
The iterative steps are performed until the difference |ϕ_{n} – ϕ_{n}_{−1}| converges to a desired small value. To solve
The procedure of computation process of fixing the ship position by single celestial body is conducted as following steps: First step: The ship position P_{1}(ϕ_{1}, λ_{1}) obtained at closest time is recorded. The latitudes ϕ_{1} and (ϕ_{1}± 1°) are initialized as two initial solutions for the positioning equation. At the time UTC, a celestial body C is identified and selected for fixing.
Second step: The coordinate of the celestial body is also determined at the time UTC. The observer (on the ship) measures the altitude (H) and the azimuth (A) of the celestial body.
Third step: Using
Final step:
To practice the fixing steps with this proposed method, the detail is illustrated in flowchart as shown in Figure 5.
In this section, two numerical examples were carried out to illustrate the content of the proposed approach for celestial navigation method. The gulf of Tonkin of Viet Nam was chosen to take sight the celestial body for these examples. The detail of numerical examples is shown as follows:
On April 1, 2019, a ship travels with speed 14 knots and course 112°. The closest position of the ship before observing the celestial body is at (20°44′6N, 107°06′4E). The star named MIRFAK was chosen to observe the altitude and the azimuth simultaneously. The true position of ship at observation time was recorded at (20°37′5N, 107°16′3E). At time UTC =11:40:07, the altitude and the azimuth of MIRFAK were measured as (H =34°31′9) and (A =317°18).
At time UTC =11:40:07, the coordinate of MIRFAK was determined as (GHA =313°09′5, δ =49°55′7N). Using the two initial solutions ϕ_{1} =21°N and ϕ_{2} =20°N for
By applying
The solution of latitude by Secant technique was obtained after three iterations (ϕ =20.62519° =20°37′5N).
The longitude was found based on
On April 30, 2019, a vessel moves with speed 12.5 knots and course 120°. The closest position of the ship before observing the celestial body is at (20°40′9N, 107°12′8E). The star named DUBHE was chosen to observe the altitude and the azimuth simultaneously. The true position of ship at observation time was recorded at (20°45′2N, 107°20′6E). At time UTC =10:51:22, the altitude and the azimuth of DUBHE were measured as (H =41°27′1) and (A =22°93).
At time UTC = 10:51:22, the coordinate of DUBHE was determined as (GHA =214°42′5, δ =61°39′1N). Using the two initial solutions ϕ_{1} =21°N and ϕ_{2} =20°N for
Similarly, the equation system for the latitude was solved by applying
The longitude was found based on
From above-mentioned examples, it can find that the new approach is possible to apply theorically to numerical problems. In the future, this technique will be investigated with more examples and real expriments to verify it’s feasibility. In addition, the external disturbances such as wind, current, and wave can be also considered for this celestial navigation technique.
In this study, a new approach is proposed for fixing based on a single body and secant technique to obtain the ship position. The conclusion is summarized as follows:
•First, the mathematical equations relating the ship position to celestial body were described based on an altitude and an azimuth of celestial body. These equations were established using knowledge of spherical trigonometry.
•Second, this proposed technique used both parameters which are the altitude and the azimuth of a single celestial body; therefore this technique can perform with a single celestial body at a specified time when there is only a single celestial body identified.
•Finally, applying secant technique to solve the system is seen as one of new approaches of advanced mathematics in celestial navigation. This approach can be a useful reference for navigators and marine students.
Although numerical examples were carried out by the proposed approach, some limitations were still found out in this method. One of them is a large error in observed azimuth. In the future, it is necessary to research the effect of the error in observed azimuth on the ship position.
No potential conflict of interest relevant to this article was reported.
The relation between the ship position and celestial body on celestial sphere.
The relation between longitude and hour angles.
A spherical trigonometry for obtaining the ship position.
A spherical trigonometry for obtaining the ship position.
The flowchart of the proposed fixing method.
Variables and symbols
Symbol | Variable | Interval |
---|---|---|
Greenwich hour angle | 0 ≤ | |
Local hour angle | 0 ≤ | |
Declination | −90° (S)≤ | |
Latitude | −90° (S)≤ | |
Longitude | −180° (W)≤ | |
Azimuth of celestial body | 0 ≤ | |
Altitude of celestial body | 0 ≤ |
Result of obtaining latitude by secant method (Example 1)
Iterative loop | Inputs (deg) | Solution (deg) | Iterative error ( | |
---|---|---|---|---|
Iter1 | 0.00741063 | 20.629904 | 0.62990432 | |
Iter2 | 20.629904 | 9.39195E-05 | 20.625249 | 0.00465578 |
Iter3 | 20.625249 | 1.17796E-06 | 20.625189 | 5.9135E-05 |
Result of obtaining latitude by secant method (Example 2)
Iterative loop | Inputs (deg) | Solution (deg) | Iterative error ( | |
---|---|---|---|---|
Iter1 | 0.0235923 | 20.76205 | 0.762047 | |
Iter2 | 20.76205 | 0.0007643 | 20.75441 | 0.007632 |
Iter3 | 20.75441 | 2.461E-05 | 20.75416 | 0,000254 |
Iter4 | 20.75416 | −8.19E-09 | 20.75416 | 8.44E-08 |
E-mail: nguyenvansuong@vimaru.edu.vn