International Journal of Fuzzy Logic and Intelligent Systems 2020; 20(1): 1-7
Published online March 25, 2020
https://doi.org/10.5391/IJFIS.2020.20.1.1
© The Korean Institute of Intelligent Systems
Van Suong Nguyen
Faculty of Navigation, Vietnam Maritime University, 484 Lachtray Str., Hai Phong City, Vietnam
Correspondence to :
Van Suong Nguyen (nguyenvansuong@vimaru.edu.vn)
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
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.
Keywords: Celestial navigation, Astronomical ship positioning, Single body, Secant method, Numerical method
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 (
Generally,
If
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
From
Substituting
Placing function (
The
The iterative steps are performed until the difference
The procedure of computation process of fixing the ship position by single celestial body is conducted as following steps: First step: The ship position
Second step: The coordinate of the celestial body is also determined at the time
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
At time
By applying
The solution of latitude by Secant technique was obtained after three iterations (
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
At time UTC = 10:51:22, the coordinate of DUBHE was determined as (
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.
Table 1. 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 ≤ |
Table 2. 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 |
Table 3. 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
International Journal of Fuzzy Logic and Intelligent Systems 2020; 20(1): 1-7
Published online March 25, 2020 https://doi.org/10.5391/IJFIS.2020.20.1.1
Copyright © The Korean Institute of Intelligent Systems.
Van Suong Nguyen
Faculty of Navigation, Vietnam Maritime University, 484 Lachtray Str., Hai Phong City, Vietnam
Correspondence to:Van Suong Nguyen (nguyenvansuong@vimaru.edu.vn)
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
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.
Keywords: Celestial navigation, Astronomical ship positioning, Single body, Secant method, Numerical method
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 (
Generally,
If
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
From
Substituting
Placing function (
The
The iterative steps are performed until the difference
The procedure of computation process of fixing the ship position by single celestial body is conducted as following steps: First step: The ship position
Second step: The coordinate of the celestial body is also determined at the time
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
At time
By applying
The solution of latitude by Secant technique was obtained after three iterations (
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
At time UTC = 10:51:22, the coordinate of DUBHE was determined as (
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.
Table 1 . 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 ≤ |
Table 2 . 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 |
Table 3 . 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 |
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.