Solving the geodesics on the ellipsoid as a boundary value problem Research Article
|
|
- Maximilian Atkins
- 5 years ago
- Views:
Transcription
1 31) DOI: /jogs Journal of Geodetic Science Solving the geodesics on the ellipsoid as a boundary value problem Research Article G. Panou, D. Delikaraoglou, R. Korakitis Department of Surveying Engineering, National Technical University of Athens, Zografou Campus, Athens, Greece Abstract: The geodesic between two given points on an ellipsoid is determined as a numerical solution of a boundary value problem. The secondorder ordinary differential equation of the geodesic is formulated by means of the Euler-Lagrange equation of the calculus of variations. Using Taylor s theorem, the boundary value problem with Dirichlet conditions at the end points is replaced by an initial value problem with Dirichlet and Neumann conditions. The Neumann condition is determined iteratively by solving a system of four rst-order differential equations with numerical integration. Once the correct Neumann value has been computed, the solution of the boundary value problem is also obtained. Using a special case of the Euler-Lagrange equation, the Clairaut equation is veri ed and the Clairaut constant is precisely determined. The azimuth at any point along the geodesic is computed by a simple formula. The geodesic distance between two points, as a de nite integral, is computed by numerical integration. The numerical tests are validated by comparison to Vincenty s inverse formulas. Keywords: boundary value problem Clairaut constant geodesics inverse geodesic problem numerical integration Versita sp. z o.o. Received ; accepted Introduction The shortest path between two points on an ellipsoid spheroid) is along a geodesic curve more simply, a geodesic. There are two main problems in geodesy related to the geodesics: Direct problem: Given a point P 0 ϕ 0,λ 0 ) on the ellipsoid with geodetic latitude ϕ 0 and geodetic longitude λ 0, together with the geodesic distance s and the azimuth α 01 to a point P 1 ϕ 1,λ 1 ), determine the geodetic coordinates ϕ 1, λ 1 and the azimuth α 10 at P 1 ϕ 1,λ 1 ). Inverse problem: Given two points P 0 ϕ 0,λ 0 ) and P 1 ϕ 1,λ 1 ) on the ellipsoid, determine the geodesic distance s between them and the azimuths α 01, α 10 at the end points. Today, with Global Navigation Satellite System GNSS) technologies, this problem is more realistic than the direct problem. A historical summary of solution methods for these problems can be found in Rapp 1993), Deakin and Hunter 2010), and Karney 2013). Among these methods, Vincenty s iterative formulas based on series expansions are widely used Vincenty, 1975). Recently, Karney 2013) gave improved series expansions for the direct case of the geodesic as well as a method for solving the inverse case. However, Sjöberg and Shirazian 2012) solved the direct and inverse problem by decomposing the solutions into those on a sphere and the corrections for the ellipsoid. The spherical solutions are given in closed form, while the corrections for an ellipsoid are expressed with elliptic integrals, suitable for numerical integration. A similar approach is followed by Saito 1970). Also, part of the inverse problem is the determination of the Clairaut constant, which was treated by Sjöberg 2007). Today, considering modern computational capabilities, we prefer solution methods that use numerical integration rather than a series expansion approach, because truncated series solution inevitably makes a mathematical approximation. By comparison, numerical integration suffers only from computational errors, which can be addressed with improved geopanou@survey.ntua.gr
2 Journal of Geodetic Science 41 computational systems and require no change in the theoretical background. Kivioja 1971) and Jank and Kivioja 1980) have solved the direct problem by numerical integration of two of the basic differential equations of the geodesic, using the geodesic distance s as the independent variable. Thomas and Featherstone 2005) improved Kivioja s method, showing numerical agreement with Vincenty s direct formulas. Another advantage of their method is that it provides a numerically efficient and convenient approach to plot the geodesic i.e., to produce the coordinates and azimuths at any point along the geodesic). However, as Thomas and Featherstone 2005) point out, the major disadvantage of this method is that it is not possible to invert the formulas to obtain a closed-form solution to the inverse problem. In this work we present a method which solves the geodesic problem. This method is based on the calculus of variations and uses numerical integration techniques. It principally addresses what is traditionally known as the inverse problem, but it can also be used to plot the geodesic and determine the Clairaut constant without using the Clairaut equation between two points on a biaxial ellipsoid. In addition, we do not use conformal mapping with an auxiliary sphere, as do Saito 1970), Vincenty 1975) and Karney 2013). Finally, our approach can be generalized to describe the geodesics on a triaxial ellipsoid, the Clairaut equation does not hold. 2. Geodesics as a boundary value problem We consider a biaxial ellipsoid which, in Cartesian coordinates, is described by x 2 a 2 + y2 a 2 + z2 b 2 = 1 1) a and b are its two semiaxes a > b). From these we can compute the rst eccentricity by e = a 2 b 2 /a. In geodesy, it is well-known that this ellipsoid is described parametrically by x = N cos ϕ cos λ y = N cos ϕ sin λ z = N 1 e 2) sin ϕ 2a) 2b) 2c) π/2 ϕ +π/2 is the geodetic latitude, π < λ +π is the geodetic longitude and a N = 1 e 2 sin 2 ϕ is the radius of curvature in the prime vertical normal section. In this parametrization, the rst fundamental coefficients E, F, and G are Deakin and Hunter, 2008) 3) E = a2 1 e 2) 2 1 e2 sin 2 ϕ ) 3 4) G = F = 0 5) a2 cos 2 ϕ 1 e 2 sin 2 ϕ. 6) In Eq. 5), F = 0 indicates that the ϕ-curves parallels) and λ-curves meridians) are orthogonal. Also, E 0 for all ϕ and G = 0 when ϕ = ± π/2 at the poles). Let us consider that a curve on the ellipsoid is described by ϕ = ϕλ) 7) i.e., that the geodetic latitude depends on the geodetic longitude. The element of distance ds on the ellipsoid is given by Rapp, 1984) ds = Eϕ ) 2 + Gdλ 8) ϕ = dϕ dλ. 9) Hence, the length s of a curve ϕ = ϕλ) from λ = λ 0 to λ = λ 1 λ 0 < λ 1 ) is given by λ1 s = fϕ, ϕ)dλ 10) λ 0 fϕ, ϕ) = Eϕ ) 2 + G. 11) We assume that λ 0 λ 1 since the case λ 0 = λ 1 can be excluded as a trivial one: all meridians on the ellipsoid are geodesics, the azimuths α along the meridian are 0 or π, and the geodesic distance s between two points located at the meridian can be computed by the well-known meridian distance formula Rapp, 1984). Moreover, without loss of generality, we shall always consider that λ 0 < λ 1, which implies that 0 < α < π. Also, the poles are excluded and hence it holds that G 0. The computation of the geodesic between two points P 0 ϕ 0,λ 0 ) and P 1 ϕ 1,λ 1 ) on the ellipsoid entails determining the curve ϕ = ϕλ) with ϕ 0 = ϕλ 0 ) and ϕ 1 = ϕλ 1 ) such that the length in Eq. 10) is minimum. This implies that the smooth) geodesic ϕ = ϕλ) must satisfy the Euler-Lagrange equation of the calculus of variations van Brunt, 2004; Logan, 2006) ) d f f = 0. 12) dλ ϕ ϕ
3 42 Journal of Geodetic Science In our case, using Eq. 11) we obtain 3. From a boundary to an initial value problem and f Eϕ = ϕ Eϕ ) 2 + G 13) The method which we propose reduces the boundary value problem given by Eqs. 18) and 19) to an initial value problem. Subsequently, the initial value problem can be solved by well-known numerical techniques. First, Eq. 18) is written equivalently as a system of two differential equations of the rst-order f ϕ = E ϕ ) 2 + G 2 Eϕ ) 2 + G 14) d dλ ϕ) = f 1 ϕ, ϕ ) 20a) E = de dϕ, G = dg dϕ. 15) d ) ϕ = f 2 ϕ, ϕ ) 20b) dλ By writing out the total derivative in Eq. 12) using the chain rule, the Euler-Lagrange equation becomes 2 f ϕ ϕ ϕ + 2 f ϕ ϕ ϕ + 2 f f λ ϕ ϕ = 0 16) and with f 1 ϕ, ϕ ) = ϕ 21) f 2 ϕ, ϕ ) = h 1 ϕ ) 2 + h2 22) and the term 2 f λ ϕ = λ ϕ = d2 ϕ dλ 2 17) ) f = 0 ϕ in our case. Substituting Eqs. 13) and 14) into Eq. 16) subsequently yields h 1 = G G 1 E 2 E 23) h 2 = 1 G 2 E. 24) The initial values associated with this system are 2EGϕ 2EG E G ) ϕ ) 2 GG = 0 18) D : ϕ 0 = ϕλ 0 ) 25a) which is a non-linear second-order ordinary differential equation. The boundary values associated with this equation are N : ϕ 0 = ϕ λ 0 ) 25b) ϕ 0 = ϕλ 0 ) 19a) ϕ 1 = ϕλ 1 ) 19b) which are known as Dirichlet conditions. Hence, the geodesic between two points on the ellipsoid is described by a two-point boundary value problem. Generally, there are several numerical approaches for solving a two-point boundary value problem, such as shooting methods, nite differences, and collocation or nite element methods see, e.g., Fox, 1990; and Keller, 1992). However, in next section we develop a method based on Taylor s theorem. D and N denote Dirichlet and Neumann conditions, respectively, and ϕ 0 is an unknown value. Obviously, with given initial values Eqs. 25)), the system of Eqs. 20) can be numerically integrated, allowing determination of the solution ϕ = ϕλ) which corresponds to initial values Eqs. 25)). For this reason, the geodesic can be better described by Our aim is to determine ϕ 0 such that ϕ = ϕϕ 0, ϕ 0; λ) 26) ϕ 1 = ϕϕ 0, ϕ 0; λ 1 ). 27)
4 Journal of Geodetic Science 43 We start with an approximate value ϕ 0) 0 28) and we integrate the system of Eqs. 20) using any convenient numerical method on the interval [λ 0, λ 1 ]. Thus, we determine the geodesic see Fig. 1)) ϕ 1 North α 01 Γ Γ 1). Γ 0) P 1 ϕ 1, λ 1 ) P 1) 1 P 0) 1 ) ϕ 1) 1, λ 1 ϕ 0) 1, λ 1 ) with Γ 0) : ϕ = ϕϕ 00), ϕ 0 ; λ) 29) ϕ 0 P 0 ϕ 0, λ 0 ) λ 0 λ 1 ϕϕ 0) 0, ϕ 0 ; λ 1 ) = ϕ 0) 1 ϕ 1. 30) Figure 1. The geodesic on the ellipsoid. Therefore, we search for a correction δϕ 00) such that ϕϕ 00) + δϕ 00), ϕ 0 ; λ 1 ) = ϕ 1. 31) and via numerical integration on the interval [λ 0, λ 1 ] we determine the geodesic Using Taylor s theorem second and higher order terms ignored), Eq. 31) can be written as Γ 1) : ϕ = ϕϕ 01), ϕ 0 ; λ) 37) ) ϕ 1 = ϕϕ 00) ϕ, ϕ 0 ; λ 1 ) + δϕ 0) ) ) with ϕϕ 01), ϕ 0 ; λ 1 ) = ϕ 1 1) ϕ 1. 38) and from Eqs. 30), 31) and 32) we then obtain δϕ 0 0) = ϕ 1 ϕ 1 0) ϕ 0) )1. 33) In Eq. 33) the derivative has an unknown value. In order to solve this problem we apply the chain rule in Eqs. 20) to obtain ) d ϕ = f 1 dλ ϕ ) d ϕ = f 2 dλ ϕ ϕ ϕ + f 1 ϕ ϕ 34a) + f 2 ϕ. 34b) ϕ Hence, we can integrate the system of Eqs. 20) and 34) on the interval [λ 0, λ 1 ] and obtain at λ 1 the values ) ϕ ϕ 0) 1, ) ) which are required in Eq. 33). In other words, by integrating the system of Eqs. 20) and 34), we obtain the geodesic Γ 0) and the value δϕ 00) which is required to start a new iteration. Now, we start with the value Using the results at λ 1 ) ϕ ϕ 1) 1, ϕ 1) ) we compute the new correction δϕ 01). The process is repeated m times until we reach a value such that ϕ m) 0 40) ϕ1 m) ϕ 1 < ε 41) ε > 0 is a user-de ned threshold for the desired accuracy. We observe that, together with the correct value ϕ 0, this computation yields the geodesic Γ as well. The details related to the numerical integration are included in the next section. 4. Numerical integration Introducing the variables ϕ 01) = ϕ 00) + δϕ 0) 0 36) x 1 = ϕ 42a)
5 44 Journal of Geodetic Science x 2 = ϕ 42b) 4.1. Step size and initial conditions x 3 = ϕ x 4 = ϕ 42c) 42d) and using Eqs. 21), 22) and 34), the system of Eqs. 20) and 34) can be rewritten as x 1 = x 2 43a) The step size δλ is calculated by δλ = λ 1 λ 0 n 47) n is the number of steps. Basically, the choice of n represents a compromise between speed small n) and accuracy small δλ). For the variable x 1 the initial condition is always the geodetic latitude ϕ 0. For the variable x 2 the initial condition can be obtained by the spherical case: x 2 = h 1 x 2 ) 2 + h 2 x 3 = x 4 ) x 4 = h 3 x 2 ) 2 + h 4 x 3 + 2h 1 x 2 x 4 using Eqs. 23) and 24), h 3 = h 1 ϕ 43b) 43c) 43d) 44) A is the azimuth and ϕ 00) = cos ϕ 0 cot A 48) ) sin ϕ1 sin ϕ 0 cos ω A = arccos cos ϕ 0 sin ω 49a) ω = arccos [sin ϕ 0 sin ϕ 1 + cos ϕ 0 cos ϕ 1 cos λ 1 λ 0 )]. 49b) and h 4 = h 2 ϕ. 45) We now observe that x 2 and x 4 are expressed in terms of the x 1, x 2 and x 1, x 2, x 3, x 4, respectively. Performing the necessary manipulations, Eqs. 23), 24), 44) and 45) are then rewritten as h 1 = 2 1 e 2) tan ϕ + 3e 2 sin ϕ cos ϕ 1 e 2 sin 2 ϕ 46a) 1 e 2 sin 2 ϕ ) sin ϕ cos ϕ h 2 = 1 e 2 46b) h 3 = 2 1 e 2) 1 e 2 sin 2 ϕ + 2e 2 sin 2 ϕ cos 2 ϕ 1 e2 sin 2 ϕ ) 2 cos2 ϕ 3e 2 1 e 2 sin 2 ϕ ) cos 2 ϕ sin 2 ϕ ) + 2e 2 sin 2 ϕ cos 2 ϕ 1 e2 sin 2 ϕ ) 2 h 4 = 46c) 1 e 2 sin 2 ϕ ) cos 2 ϕ sin 2 ϕ ) 2e 2 sin 2 ϕ cos 2 ϕ 1 e 2. 46d) Hence, the system of the four rst-order differential equations 43) can be solved on the interval [λ 0, λ 1 ] using a numerical integration method such as Runge-Kutta or a Taylor series see Butcher, 1987). The required initial conditions are described below. Equations 49) may give inaccurate results, since the evaluation of the arccosine suffers for arguments very close to 1 i.e., very small angles). We cannot overemphasize that Eqs. 48) and 49) serve only to provide an initial value for x 2, such that the accuracy of this value is of no real concern. This differentiates our approach from the efforts of Saito 1970), Vincenty 1975), and others in addressing the small angle problem. Furthermore, in our method one could instead use the very simple approximation formula for the initial value of x 2. x 2 = ϕ 0) 0 = ϕ 1 ϕ 0 50) λ 1 λ 0 Subsequently, in each iteration this value is corrected according to the method discussed in Section 3. Finally, the variables x 3 and x 4 always have initial values of 0 and 1, respectively. 5. Clairaut s constant The function f in Eq. 11) does not contain the independent variable λ explicitly. Therefore, along any geodesic ϕ = ϕλ) it holds that C is a constant van Brunt, 2004). f ϕ f ϕ = C 51) Equation 51), which is a special case of the Euler-Lagrange equation Eq. 12)), has been previously seen in the literature e.g.,
6 Journal of Geodetic Science 45 Struik, 1961; Karney and Deakin, 2010). In contrast, we make use of the complete expression, Eq. 18), and the constant C in Eq. 51) is computed immediately below. Substituting Eqs. 11) and 13) into Eq. 51) we obtain G = C. 52) Eϕ ) 2 + G However, for any curve on the ellipsoid, we can write ϕ = G E cot α. 53) Substituting Eq. 53) into Eq. 52) and using Eqs. 3) and 6) yields N cos ϕ sin α = C 54) which is the Clairaut equation. Hence, Eqs. 52) and 54) are equivalent and the Clairaut constant C can be computed by Eq. 52) at any value of the independent variable λ. Furthermore, since 0 < α < π, Eq. 54) implies that 0 < C < a. We also note that Eq. 52) involves the variables x 1 = ϕ and x 2 = ϕ, which are obtained by the numerical integration. In this way, one can check the accuracy of the numerical integration and subsequently compute the azimuths along the geodesic and the geodesic distance between two given points. 6. Azimuths and geodesic distance Traditionally, once the Clairaut constant C is known, the azimuth α at any point along the geodesic can be computed by solving Eq. 54) in α. This equation has two solutions ) C α = arcsin N cos ϕ ) C α = π arcsin N cos ϕ 55a) 55b) but this problem can be eliminated, allowing computation of the correct azimuth Sjöberg and Shirazian, 2012). Also, at the vertex point i.e., the point with maximum/minimum latitude and azimuth π/2), it holds that C = N cos ϕ. It is well-known that the derivative of arcsine is large for arguments near 1. Thus, small errors in C /N cos ϕ) due to factors such as computer rounding of oating point numbers, will produce large errors in the calculation of azimuth Thomas and Featherstone, 2005). Hence, using Eqs. 55) is likely to lead to inaccurate azimuths near the vertex. In order to avoid the previous problems, we compute the azimuth α at any point along the geodesic by solving Eq. 53), α = arccot ) E ϕ G 56) which gives π/2 α π/2. Since λ 0 λ 1, α 0. When the azimuth is negative, the correct azimuth is obtained as α = α + π. One should note that Eq. 56) involves the variables x 1 = ϕ and x 2 = ϕ, which are obtained by the numerical integration. Also, the forward azimuth α 01 at P 0 ϕ 0,λ 0 ) is α 01 = α 0 and the reverse azimuth at P 1 ϕ 1,λ 1 ) is α 10 = α 1 + π. Finally, the geodesic distance s between the two points P 0 ϕ 0,λ 0 ) and P 1 ϕ 1,λ 1 ), using Eqs. 10) and 11) is written as a de nite integral λ1 s = Eϕ ) 2 + Gdλ 57) λ 0 which can be computed by a numerical integration method such as the Newton Cotes formulas see, e.g., Hildebrand, 1974). 7. Numerical tests and comparisons In order to validate the algorithm which has been presented, the results were compared to those obtained using the Vincenty s inverse formulas. Because the problem is invariant under rotations around the z-axis, only starting points ϕ 0,λ 0 ) with λ 0 = 0 and ϕ 0 = 0, 30, 60 and 75 symmetry) were selected, as well as points ϕ 1,λ 1 ) with λ 1 = 5, 40, 80, 120, 160 and 170 and ϕ 1 = 75, 60, 30, 0, 30, 60, and 75. Note that, when ϕ 0 = 0 only the values ϕ 1 0 were used. Hence, in total 150 geodesics were tested. In our selection of the points and, especially, the maximum value of 170 for the longitude difference λ 1 λ 0 ), we followed the rationale of similar works, like Thomas and Featherstone 2005) and Sjöberg and Shirazian 2012), who limit their comparisons with Vincenty s method to geodesic distances up to km. All the numerical computations were carried out using the GRS80 ellipsoid Moritz, 1980), meaning that a = m and b = m. All algorithms were implemented in MATLAB. The system of four rst-order differential equations 43) was solved using the fourth-order Runge-Kutta numerical integration method. The number of steps n was selected as 8000 in order to cover all cases with sufficient accuracy. Thus, the maximum step size δλ corresponds to minutes of arc. The latitudes at λ 1 were required to converge with an accuracy ε = rad, which corresponds to approximately mm. As a result, in all cases, the iterative procedure reached convergence in three or four iterations. In any particular geodesic, the Clairaut constant C, which was computed using Eq. 52) at all values of the independent variable λ, had a maximum discrepancy 0.02 mm. Also, the geodesic distance s between the two points, given by Eq. 57), was computed by Simpson s rule i.e. the three point rule). For this study, Vincenty s algorithm Vincenty, 1975) was implemented with the requirement that the longitude differences were to converge with an accuracy rad mm. The results between our proposed method and Vincenty s inverse method show agreement to within seconds of arc for az-
7 46 Journal of Geodetic Science Table 1. Numerical tests and comparisons with Vincenty s inverse method. ϕ 0 ϕ 1 λ 1 C m) α 01 o ) α 10 o ) s m) α 01 ) α 10 ) s mm) 0 o 75 o 5 o o 75 o 170 o o 30 o 5 o o -30 o 170 o o -75 o 40 o o -30 o 120 o o 75 o 5 o o 30 o 170 o Table 2. Special cases. ϕ 0 ϕ 1 λ 1 C m) α 01 o ) α 10 o ) s m) α 01 ) α 10 ) s mm) 0 o 75 o 170 o o 30 o 170 o o -35 o 179 o o 0.5 o o imuth and 0.12 mm mean value of differences) in geodesic distanceup to 0.58 mm in the worst case). A small sample of results is presented in Table 1. The differences in the geodesic distances that appear in the last column of Table 1 are a consequence of our initial selection of the number of steps. We calculated the geodesics for a few special cases, using twice the usual number of steps instead of 8000). These results are presented in Table 2. For geodesics between near antipodal points, Saito 1970) pointed out that all methods based on conformal mapping are problematic and he proposed a method to face this antipodal problem. More recently, Karney 2013) has also also attempted to solve this problem. Since our method is not based on a conformal mapping, we tested it using two special cases and the results are also shown in Table 2. In these cases we used a longitude difference up to and our solution converges to the required accuracy after only 5 iterations. 8. Concluding remarks The method presented here describes the geodesic between two given points on the ellipsoid as a two-point boundary value problem. Using Taylor s theorem and a numerical integration method, this problem is replaced by an initial value problem. From its solution the coordinates and the azimuths at any point along the geodesic are determined. Also, the Clairaut constant is determined together with an accuracy check. The numerical tests show that the solutions practically agree with Vincenty s inverse solutions. Hence, this method can be used as an algorithm to plot the geodesic between two given points on the ellipsoid. As an independent method it can be used to validate Vincenty s inverse method and moreover is able to provide an accurate solution to the geodesic problem even in extreme cases, such as between points nearly antipodal to one another. The present method is universal in character and thus can be used to describe the geodesics on a triaxial ellipsoid. This presents an opportunity for use in many recent geodetic works dealing with the transformation problems on the triaxial ellipsoid, such as Feltens 2009) and Ligas 2012a, 2012b). Finally, by setting e = 0, the method is reduced to the sphere, i.e. showing that the geodesics on the sphere are obtained as a degenerate case. Acknowledgments The authors would like to thank Dr. G. Manoussakis, Department of Surveying Engineering, National Technical University of Athens, for his many suggestions and helpful comments on the various aspects of this research. We extend our appreciation to the reviewers of the original manuscript for their many constructive and useful remarks. References Butcher J. C., 1987, The numerical analysis of ordinary differential equations: Runge-Kutta and general linear methods, Wiley, New York. Deakin R. E. and Hunter M. N., 2008, Geometric Geodesy - Part A, Lecture Notes, School of Mathematical & Geospatial Sciences, RMIT University, Melbourne, Australia. Deakin R. E. and Hunter M. N., 2010, Geometric Geodesy - Part B, Lecture Notes, School of Mathematical & Geospatial Sciences, RMIT University, Melbourne, Australia. Feltens J., 2009, Vector method to compute the Cartesian X, Y, Z ) to geodetic ϕ, λ, h) transformation on a triaxial ellipsoid, J. Geod., 83, Fox L., 1990, The numerical solution of two-point bound-
8 Journal of Geodetic Science ary problems in ordinary differential equations, Dover, New York. Hildebrand F. B., 1974, Introduction to numerical analysis, 2nd ed., Dover, New York. Geod., 54, Rapp R. H., 1984, Geometric Geodesy - Part I, Department of Geodetic Science and Surveying, Ohio State University, Columbus, Ohio, USA. Jank W. and Kivioja L. A., 1980, Solution of the direct and inverse problems on reference ellipsoids by point-by-point integration using programmable pocket calculators, Surveying and Mapping, 15, Karney C. F. F. and Deakin R. E., 2010, F.W. Bessel 1825): The calculation of longitude and latitude from geodesic measurements, Astron. Nachr., 331, Karney C. F. F., 2013, Algorithms for geodesics, J. Geod., 87, Keller H. B., 1992, Numerical methods for two-point boundaryvalue problems, Dover, New York. Rapp R. H., 1993, Geometric Geodesy - Part II, Department of Geodetic Science and Surveying, Ohio State University, Columbus, Ohio, USA. Saito T., 1970, The computation of long geodesics on the ellipsoid by non-series expanding procedure, Bull. Geod., 98, Sjöberg L. E., 2007, Precise determination of the Clairaut constant in ellipsoidal geodesy, Surv. Rev., 39, Sjöberg L. E. and Shirazian M., 2012, Solving the direct and inverse geodetic problems on the ellipsoid by numerical integration, J. Surv. Eng., 138, Kivioja L. A., 1971, Computation of geodetic direct and indirect problems by computers accumulating increments from geodetic line elements, Bull. Geod., 99, Ligas M., 2012a, Cartesian to geodetic coordinates conversion on a triaxial ellipsoid, J. Geod., 86, Struik, D. J., 1961, Lectures on classical differential geometry, 2nd ed., Dover, New York. Thomas C. M. and Featherstone W. E., 2005, Validation of Vincenty s formulas for the geodesic using a new fourth-order extension of Kivioja s formula, J. Surv. Eng., 131, Ligas M., 2012b, Two modi ed algorithms to transform Cartesian to geodetic coordinates on a triaxial ellipsoid, Stud. Geoph. Geod., 56, Logan J. D., 2006, Applied mathematics, 3rd ed., Wiley- Interscience, New Jersey. van Brunt B., 2004, The calculus of variations, Springer- Verlag, New York. Vincenty T., 1975, Direct and inverse solutions of geodesics on the ellipsoid with application of nested equations, Surv. Rev., 23, Moritz H., 1980, Geodetic Reference System 1980, Bull.
Lecture # 7 Coordinate systems and georeferencing
Lecture # 7 Coordinate systems and georeferencing Coordinate Systems Coordinate reference on a plane Coordinate reference on a sphere Coordinate reference on a plane Coordinates are a convenient way of
More informationRECOMMENDATION ITU-R S.1257
Rec. ITU-R S.157 1 RECOMMENDATION ITU-R S.157 ANALYTICAL METHOD TO CALCULATE VISIBILITY STATISTICS FOR NON-GEOSTATIONARY SATELLITE ORBIT SATELLITES AS SEEN FROM A POINT ON THE EARTH S SURFACE (Questions
More informationOn Surfaces of Revolution whose Mean Curvature is Constant
On Surfaces of Revolution whose Mean Curvature is Constant Ch. Delaunay May 4, 2002 When one seeks a surface of given area enclosing a maximal volume, one finds that the equation this surface must satisfy
More informationMath 215 Project 1 (25 pts) : Using Linear Algebra to solve GPS problem
Due 11:55pm Fri. Sept. 28 NAME(S): Math 215 Project 1 (25 pts) : Using Linear Algebra to solve GPS problem 1 Introduction The age old question, Where in the world am I? can easily be solved nowadays by
More informationGauss and AGM. Burton Rosenberg. January 30, 2004
Gauss and AGM Burton Rosenberg January 3, 24 Introduction derivation of equation. what has it to do w/ the lemniscate agm properties of I elliptic integrals The Elliptic Integral of the First Kind Define
More informationAbout Doppler-Fizeau effect on radiated noise from a rotating source in cavitation tunnel
PROCEEDINGS of the 22 nd International Congress on Acoustics Signal Processing in Acoustics (others): Paper ICA2016-111 About Doppler-Fizeau effect on radiated noise from a rotating source in cavitation
More information10.1 Curves defined by parametric equations
Outline Section 1: Parametric Equations and Polar Coordinates 1.1 Curves defined by parametric equations 1.2 Calculus with Parametric Curves 1.3 Polar Coordinates 1.4 Areas and Lengths in Polar Coordinates
More informationOPTIMUM GEODETIC DATUM TRANSFORMATION TECHNIQUES FOR GPS SURVEYS IN EGYPT
Proceedings of Al-Azhar Engineering Sixth International Conference, Sept. 1-, 2000, Cairo, Egypt, Volume, pp. 09-1. OPTIMUM GEODETIC DATUM TRANSFORMATION TECHNIQUES FOR GPS SURVEYS IN EGYPT By Dr. Gomaa
More informationAn NGS Illustrated Guide to Geodesy for GIS Professionals
An NGS Illustrated Guide to Geodesy for GIS Professionals Michael Dennis, RLS, PE michael.dennis@noaa.gov Esri User Conference San Diego Convention Center July 14-18, 2014 San Diego, CA Why should we care
More informationGeodesy, Geographic Datums & Coordinate Systems
Geodesy, Geographic Datums & Coordinate Systems What is the shape of the earth? Why is it relevant for GIS? 1/23/2018 2-1 From Conceptual to Pragmatic Dividing a sphere into a stack of pancakes (latitude)
More informationLECTURE 7: POLYNOMIAL CONGRUENCES TO PRIME POWER MODULI
LECTURE 7: POLYNOMIAL CONGRUENCES TO PRIME POWER MODULI 1. Hensel Lemma for nonsingular solutions Although there is no analogue of Lagrange s Theorem for prime power moduli, there is an algorithm for determining
More informationExperiment O11e Optical Polarisation
Fakultät für Physik und Geowissenschaften Physikalisches Grundpraktikum Experiment O11e Optical Polarisation Tasks 0. During preparation for the laboratory experiment, familiarize yourself with the function
More informationLOCAL IONOSPHERIC MODELLING OF GPS CODE AND CARRIER PHASE OBSERVATIONS
Survey Review, 40, 309 pp.71-84 (July 008) LOCAL IONOSPHERIC MODELLING OF GPS CODE AND CARRIER PHASE OBSERVATIONS H. Nahavandchi and A. Soltanpour Norwegian University of Science and Technology, Division
More informationPost processing of multiple GPS receivers to enhance baseline accuracy
Michigan Technological University Digital Commons @ Michigan Tech Dissertations, Master's Theses and Master's Reports - Open Dissertations, Master's Theses and Master's Reports 2011 Post processing of
More informationExam 2 Summary. 1. The domain of a function is the set of all possible inputes of the function and the range is the set of all outputs.
Exam 2 Summary Disclaimer: The exam 2 covers lectures 9-15, inclusive. This is mostly about limits, continuity and differentiation of functions of 2 and 3 variables, and some applications. The complete
More informationReview guide for midterm 2 in Math 233 March 30, 2009
Review guide for midterm 2 in Math 2 March, 29 Midterm 2 covers material that begins approximately with the definition of partial derivatives in Chapter 4. and ends approximately with methods for calculating
More informationTHE CONDUCTANCE BANDWIDTH OF AN ELEC- TRICALLY SMALL ANTENNA IN ANTIRESONANT RANGES
Progress In Electromagnetics Research B, Vol. 24, 285 301, 2010 THE CONDUCTANCE BANDWIDTH OF AN ELEC- TRICALLY SMALL ANTENNA IN ANTIRESONANT RANGES O. B. Vorobyev Stavropol Institute of Radiocommunications
More informationRec. ITU-R P RECOMMENDATION ITU-R P PROPAGATION BY DIFFRACTION. (Question ITU-R 202/3)
Rec. ITU-R P.- 1 RECOMMENDATION ITU-R P.- PROPAGATION BY DIFFRACTION (Question ITU-R 0/) Rec. ITU-R P.- (1-1-1-1-1-1-1) The ITU Radiocommunication Assembly, considering a) that there is a need to provide
More informationCHAPTER 10 Conics, Parametric Equations, and Polar Coordinates
CHAPTER Conics, Parametric Equations, and Polar Coordinates Section. Conics and Calculus.................... Section. Plane Curves and Parametric Equations.......... Section. Parametric Equations and Calculus............
More informationApplications of Monte Carlo Methods in Charged Particles Optics
Sydney 13-17 February 2012 p. 1/3 Applications of Monte Carlo Methods in Charged Particles Optics Alla Shymanska alla.shymanska@aut.ac.nz School of Computing and Mathematical Sciences Auckland University
More informationLevel I Signal Modeling and Adaptive Spectral Analysis
Level I Signal Modeling and Adaptive Spectral Analysis 1 Learning Objectives Students will learn about autoregressive signal modeling as a means to represent a stochastic signal. This differs from using
More informationThe Compensating Polar Planimeter
The Compensating Polar Planimeter Description of a polar planimeter Standard operation The neutral circle How a compensating polar planimeter compensates Show and tell: actual planimeters References (Far
More informationWhat is an Intrinsically Straight Path?
What is an Intrinsically Straight Path? Dictionary.com: locally minimizes the distance. Equivalently, it is a path of minimal curvature. Yan-Bin Jia: The term geodesic comes from the science of geodesy,
More informationDyck paths, standard Young tableaux, and pattern avoiding permutations
PU. M. A. Vol. 21 (2010), No.2, pp. 265 284 Dyck paths, standard Young tableaux, and pattern avoiding permutations Hilmar Haukur Gudmundsson The Mathematics Institute Reykjavik University Iceland e-mail:
More informationMath 122: Final Exam Review Sheet
Exam Information Math 1: Final Exam Review Sheet The final exam will be given on Wednesday, December 1th from 8-1 am. The exam is cumulative and will cover sections 5., 5., 5.4, 5.5, 5., 5.9,.1,.,.4,.,
More informationMath 180 Chapter 6 Lecture Notes. Professor Miguel Ornelas
Math 180 Chapter 6 Lecture Notes Professor Miguel Ornelas 1 M. Ornelas Math 180 Lecture Notes Section 6.1 Section 6.1 Verifying Trigonometric Identities Verify the identity. a. sin x + cos x cot x = csc
More informationLimits on GNSS Performance at High Latitudes
Limits on GNSS Performance at High Latitudes Peter F. Swaszek, University of Rhode Island Richard J. Hartnett, U.S. Coast Guard Academy Kelly C. Seals, U.S. Coast Guard Academy Joseph D. Siciliano, U.S.
More informationTrigonometry. David R. Wilkins
Trigonometry David R. Wilkins 1. Trigonometry 1. Trigonometry 1.1. Trigonometric Functions There are six standard trigonometric functions. They are the sine function (sin), the cosine function (cos), the
More informationMAT187H1F Lec0101 Burbulla
Spring 17 What Is A Parametric Curve? y P(x, y) x 1. Let a point P on a curve have Cartesian coordinates (x, y). We can think of the curve as being traced out as the point P moves along it. 3. In this
More informationThe evaluation of the positioning accuracy of the EGNOS and DGPS systems based on the long-term measurements in the years
MARIUSZ SPECHT Department of Geodesy and Oceanography of Gdynia Maritime University mariuszspecht199@gmail.com Polish Cartographical Review Vol. 47, 015, no., pp. 99 108 DOI: 10.1515/pcr-015-0006 The evaluation
More informationRECOMMENDATION ITU-R S *
Rec. ITU-R S.1339-1 1 RECOMMENDATION ITU-R S.1339-1* Rec. ITU-R S.1339-1 SHARING BETWEEN SPACEBORNE PASSIVE SENSORS OF THE EARTH EXPLORATION-SATELLITE SERVICE AND INTER-SATELLITE LINKS OF GEOSTATIONARY-SATELLITE
More informationMAT01B1: Calculus with Polar coordinates
MAT01B1: Calculus with Polar coordinates Dr Craig 23 October 2018 My details: acraig@uj.ac.za Consulting hours: Monday 14h40 15h25 Thursday 11h30 12h55 Friday (this week) 11h20 12h25 Office C-Ring 508
More informationSOLUTIONS TO PROBLEM SET 5. Section 9.1
SOLUTIONS TO PROBLEM SET 5 Section 9.1 Exercise 2. Recall that for (a, m) = 1 we have ord m a divides φ(m). a) We have φ(11) = 10 thus ord 11 3 {1, 2, 5, 10}. We check 3 1 3 (mod 11), 3 2 9 (mod 11), 3
More informationCorresponding author: Rebecca Woodgate,
Correction of Teledyne Acoustic Doppler Current Profiler (ADCP) Bottom-Track Range Measurements for Instrument Pitch and Roll Rebecca A. Woodgate 1 and Alexander E. Holroyd 1 Applied Physics Laboratory,
More information2.1 Partial Derivatives
.1 Partial Derivatives.1.1 Functions of several variables Up until now, we have only met functions of single variables. From now on we will meet functions such as z = f(x, y) and w = f(x, y, z), which
More informationLecture 3: Geometrical Optics 1. Spherical Waves. From Waves to Rays. Lenses. Chromatic Aberrations. Mirrors. Outline
Lecture 3: Geometrical Optics 1 Outline 1 Spherical Waves 2 From Waves to Rays 3 Lenses 4 Chromatic Aberrations 5 Mirrors Christoph U. Keller, Leiden Observatory, keller@strw.leidenuniv.nl Lecture 3: Geometrical
More informationSolutions to the problems from Written assignment 2 Math 222 Winter 2015
Solutions to the problems from Written assignment 2 Math 222 Winter 2015 1. Determine if the following limits exist, and if a limit exists, find its value. x2 y (a) The limit of f(x, y) = x 4 as (x, y)
More informationArkansas Tech University MATH 2924: Calculus II Dr. Marcel B. Finan. Figure 50.1
50 Polar Coordinates Arkansas Tech University MATH 94: Calculus II Dr. Marcel B. Finan Up to this point we have dealt exclusively with the Cartesian coordinate system. However, as we will see, this is
More informationA Direct 2D Position Solution for an APNT-System
A Direct 2D Position Solution for an APNT-System E. Nossek, J. Dambeck and M. Meurer, German Aerospace Center (DLR), Institute of Communications and Navigation, Germany Technische Universität München (TUM),
More informationUniversity of California, Berkeley Department of Mathematics 5 th November, 2012, 12:10-12:55 pm MATH 53 - Test #2
University of California, Berkeley epartment of Mathematics 5 th November, 212, 12:1-12:55 pm MATH 53 - Test #2 Last Name: First Name: Student Number: iscussion Section: Name of GSI: Record your answers
More informationChannel Estimation for OFDM Systems in case of Insufficient Guard Interval Length
Channel Estimation for OFDM ystems in case of Insufficient Guard Interval Length Van Duc Nguyen, Michael Winkler, Christian Hansen, Hans-Peter Kuchenbecker University of Hannover, Institut für Allgemeine
More information10 Antenna gain, beam pattern, directivity
10 Antenna gain, beam pattern, directivity Adipoleantenna(oracloselyrelatedmonopoletobestudiedinLecture 18) is a near perfect radiator for purposes of broadcasting that is, sending waves of equal amplitudes
More informationCHAPTER 3. BASIC GEODESY
CHAPTER 3. BASIC GEODESY SECTION I. THE GEODETIC SYSTEM A geodetic system serves as a framework for determining coordinates on the Earth s surface with respect to a reference ellipsoid and the geoid. It
More informationCHAPTER 5 The Parallel Resonant Converter
CHAPTER 5 The Parallel Resonant Converter T he objective of this chapter is to describe the operation of the parallel resonant converter in detail. The concepts developed in chapter 3 are used to derive
More informationMATH Exam 2 Solutions November 16, 2015
MATH 1.54 Exam Solutions November 16, 15 1. Suppose f(x, y) is a differentiable function such that it and its derivatives take on the following values: (x, y) f(x, y) f x (x, y) f y (x, y) f xx (x, y)
More informationPredicting Low Earth Orbit Satellite Communications Quality and Visibility over Time
1 Predicting Low Earth Orbit Satellite Communications Quality and Visibility over Time John S. ilmore and Riaan Wolhuter Department of Electronic Engineering University of Stellenbosch 7600 Stellenbosch,
More informationTrigonometry. An Overview of Important Topics
Trigonometry An Overview of Important Topics 1 Contents Trigonometry An Overview of Important Topics... 4 UNDERSTAND HOW ANGLES ARE MEASURED... 6 Degrees... 7 Radians... 7 Unit Circle... 9 Practice Problems...
More informationMULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.
Trigonometry Final Exam Study Guide Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. The graph of a polar equation is given. Select the polar
More informationLecture 2: Geometrical Optics. Geometrical Approximation. Lenses. Mirrors. Optical Systems. Images and Pupils. Aberrations.
Lecture 2: Geometrical Optics Outline 1 Geometrical Approximation 2 Lenses 3 Mirrors 4 Optical Systems 5 Images and Pupils 6 Aberrations Christoph U. Keller, Leiden Observatory, keller@strw.leidenuniv.nl
More informationTHE SPECTRAL METHOD FOR PRECISION ESTIMATE OF THE CIRCLE ACCELERATOR ALIGNMENT
II/201 THE SPECTRAL METHOD FOR PRECISION ESTIMATE OF THE CIRCLE ACCELERATOR ALIGNMENT Jury Kirochkin Insitute for High Energy Physics, Protvino, Russia Inna Sedelnikova Moscow State Building University,
More informationTrigonometric identities
Trigonometric identities An identity is an equation that is satisfied by all the values of the variable(s) in the equation. For example, the equation (1 + x) = 1 + x + x is an identity. If you replace
More informationFUNCTIONS OF SEVERAL VARIABLES AND PARTIAL DIFFERENTIATION
FUNCTIONS OF SEVERAL VARIABLES AND PARTIAL DIFFERENTIATION 1. Functions of Several Variables A function of two variables is a rule that assigns a real number f(x, y) to each ordered pair of real numbers
More informationFiber Optic Communication Systems. Unit-04: Theory of Light. https://sites.google.com/a/faculty.muet.edu.pk/abdullatif
Unit-04: Theory of Light https://sites.google.com/a/faculty.muet.edu.pk/abdullatif Department of Telecommunication, MUET UET Jamshoro 1 Limitations of Ray theory Ray theory describes only the direction
More informationAmplitude, Reflection, and Period
SECTION 4.2 Amplitude, Reflection, and Period Copyright Cengage Learning. All rights reserved. Learning Objectives 1 2 3 4 Find the amplitude of a sine or cosine function. Find the period of a sine or
More informationReview Sheet for Math 230, Midterm exam 2. Fall 2006
Review Sheet for Math 230, Midterm exam 2. Fall 2006 October 31, 2006 The second midterm exam will take place: Monday, November 13, from 8:15 to 9:30 pm. It will cover chapter 15 and sections 16.1 16.4,
More informationCalculation of antenna radiation center using angular momentum
Calculation of antenna radiation center using angular momentum Fridén, Jonas; Kristensson, Gerhard Published in: 7th European Conference on Antennas and Propagation (EuCAP), 2013 2013 Link to publication
More informationNomograms for Synthesizing Crank Rocker Mechanism with a Desired Optimum Range of Transmission Angle
International Journal of Mining, Metallurgy & Mechanical Engineering (IJMMME Volume 3, Issue 3 (015 ISSN 30 4060 (Online Nomograms for Synthesizing Crank Rocker Mechanism with a Desired Optimum Range of
More informationTOPOLOGY, LIMITS OF COMPLEX NUMBERS. Contents 1. Topology and limits of complex numbers 1
TOPOLOGY, LIMITS OF COMPLEX NUMBERS Contents 1. Topology and limits of complex numbers 1 1. Topology and limits of complex numbers Since we will be doing calculus on complex numbers, not only do we need
More informationCivil Engineering Journal
Available online at www.civilejournal.org Civil Engineering Journal Vol. 2, No. 4, April, 2016 Marine Current Meter Calibration Using GNSS Receivers, a Comparison with Commercial Method Vahid Rezaali a
More informationii) When light falls on objects, it reflects the light and when the reflected light reaches our eyes then we see the objects.
Light i) Light is a form of energy which helps us to see objects. ii) When light falls on objects, it reflects the light and when the reflected light reaches our eyes then we see the objects. iii) Light
More informationUNIT Explain the radiation from two-wire. Ans: Radiation from Two wire
UNIT 1 1. Explain the radiation from two-wire. Radiation from Two wire Figure1.1.1 shows a voltage source connected two-wire transmission line which is further connected to an antenna. An electric field
More informationPREFACE. National Geographic Department would like to express our sincere thanks for your comments.
PREFACE According to the role of National Geographic Department on Prim Minister s Decree No 255 PM, dated August 16, 2005 regarding to Surveying, Aerial Photography and mapping activities in the territory
More informationDefinitions and claims functions of several variables
Definitions and claims functions of several variables In the Euclidian space I n of all real n-dimensional vectors x = (x 1, x,..., x n ) the following are defined: x + y = (x 1 + y 1, x + y,..., x n +
More informationP282 Two-point Paraxial Traveltime in Inhomogeneous Isotropic/Anisotropic Media - Tests of Accuracy
P8 Two-point Paraxial Traveltime in Inhomogeneous Isotropic/Anisotropic Media - Tests of Accuracy U. Waheed* (King Abdullah University of Science & Technology), T. Alkhalifah (King Abdullah University
More informationLower Bounds for the Number of Bends in Three-Dimensional Orthogonal Graph Drawings
ÂÓÙÖÒÐ Ó ÖÔ ÐÓÖØÑ Ò ÔÔÐØÓÒ ØØÔ»»ÛÛÛº ºÖÓÛÒºÙ»ÔÙÐØÓÒ»» vol.?, no.?, pp. 1 44 (????) Lower Bounds for the Number of Bends in Three-Dimensional Orthogonal Graph Drawings David R. Wood School of Computer Science
More informationLecture 2: Geometrical Optics. Geometrical Approximation. Lenses. Mirrors. Optical Systems. Images and Pupils. Aberrations.
Lecture 2: Geometrical Optics Outline 1 Geometrical Approximation 2 Lenses 3 Mirrors 4 Optical Systems 5 Images and Pupils 6 Aberrations Christoph U. Keller, Leiden Observatory, keller@strw.leidenuniv.nl
More informationTSUNAMI WAVE PROPAGATION ALONG WAVEGUIDES. Andrei G. Marchuk
TSUNAMI WAVE PROPAGATION ALONG WAVEGUIDES Andrei G. Marchuk Institute of Computational Mathematics and Mathematical Geophysics Siberian Division Russian Academy of Sciences, 639, Novosibirsk, RUSSIA mag@omzg.sscc.ru
More informationMath 148 Exam III Practice Problems
Math 48 Exam III Practice Problems This review should not be used as your sole source for preparation for the exam. You should also re-work all examples given in lecture, all homework problems, all lab
More information2.1 BASIC CONCEPTS Basic Operations on Signals Time Shifting. Figure 2.2 Time shifting of a signal. Time Reversal.
1 2.1 BASIC CONCEPTS 2.1.1 Basic Operations on Signals Time Shifting. Figure 2.2 Time shifting of a signal. Time Reversal. 2 Time Scaling. Figure 2.4 Time scaling of a signal. 2.1.2 Classification of Signals
More information10.3 Polar Coordinates
.3 Polar Coordinates Plot the points whose polar coordinates are given. Then find two other pairs of polar coordinates of this point, one with r > and one with r
More informationPermutation group and determinants. (Dated: September 19, 2018)
Permutation group and determinants (Dated: September 19, 2018) 1 I. SYMMETRIES OF MANY-PARTICLE FUNCTIONS Since electrons are fermions, the electronic wave functions have to be antisymmetric. This chapter
More information1111: Linear Algebra I
1111: Linear Algebra I Dr. Vladimir Dotsenko (Vlad) Lecture 7 Dr. Vladimir Dotsenko (Vlad) 1111: Linear Algebra I Lecture 7 1 / 8 Invertible matrices Theorem. 1. An elementary matrix is invertible. 2.
More informationMath 259 Winter Recitation Handout 9: Lagrange Multipliers
Math 259 Winter 2009 Recitation Handout 9: Lagrange Multipliers The method of Lagrange Multipliers is an excellent technique for finding the global maximum and global minimum values of a function f(x,
More information1. Measure angle in degrees and radians 2. Find coterminal angles 3. Determine the arc length of a circle
Pre- Calculus Mathematics 12 5.1 Trigonometric Functions Goal: 1. Measure angle in degrees and radians 2. Find coterminal angles 3. Determine the arc length of a circle Measuring Angles: Angles in Standard
More informationTHE TAYLOR EXPANSIONS OF tan x AND sec x
THE TAYLOR EXPANSIONS OF tan x AND sec x TAM PHAM AND RYAN CROMPTON Abstract. The report clarifies the relationships among the completely ordered leveled binary trees, the coefficients of the Taylor expansion
More informationThe Casey angle. A Different Angle on Perspective
A Different Angle on Perspective Marc Frantz Marc Frantz (mfrantz@indiana.edu) majored in painting at the Herron School of Art, where he received his.f.a. in 1975. After a thirteen-year career as a painter
More informationCut-off of Resonant Modes in Truncated Conical Cavities
Cut-off of Resonant Modes in Truncated Conical Cavities José J. Rodal, Ph.D. June 2015 Although the fact that a truncated conical cavity displays an absence of sharp cut-off frequencies has been remarked
More informationCHAPTER 2 GEODESY AND DATUMS IN NAVIGATION
CHAPTER 2 GEODESY AND DATUMS IN NAVIGATION GEODESY, THE BASIS OF CARTOGRAPHY 200. Definition Geodesy is the application of mathematics to model the size and shape of the physical earth, enabling us to
More informationPenetration of VLF Radio Waves through the Ionosphere
Penetration of VLF Radio Waves through the Ionosphere By Ken-ichi MAEDA and Hiroshi OYA Kyoto University, Kyoto, Japan (Read May 24; Received November 25, 1962) Abstract The rate of energy penetration
More informationAs the Planimeter s Wheel Turns
As the Planimeter s Wheel Turns December 30, 2004 A classic example of Green s Theorem in action is the planimeter, a device that measures the area enclosed by a curve. Most familiar may be the polar planimeter
More informationUnit 5 Investigating Trigonometry Graphs
Mathematics IV Frameworks Student Edition Unit 5 Investigating Trigonometry Graphs 1 st Edition Table of Contents INTRODUCTION:... 3 What s Your Temperature? Learning Task... Error! Bookmark not defined.
More informationA Planar Equiangular Spiral Antenna Array for the V-/W-Band
207 th European Conference on Antennas and Propagation (EUCAP) A Planar Equiangular Spiral Antenna Array for the V-/W-Band Paul Tcheg, Kolawole D. Bello, David Pouhè Reutlingen University of Applied Sciences,
More informationThe Mathematics of Geometrical and Physical Optics
Orestes N. Stavroudis The Mathematics of Geometrical and Physical Optics The fc-function and its Ramifications WILEY- VCH WILEY-VCH Verlag GmbH & Co. KGaA I Preliminaries 1 1 Fermat's Principle and the
More information2. Be able to evaluate a trig function at a particular degree measure. Example: cos. again, just use the unit circle!
Study Guide for PART II of the Fall 18 MAT187 Final Exam NO CALCULATORS are permitted on this part of the Final Exam. This part of the Final exam will consist of 5 multiple choice questions. You will be
More informationMinimal generating sets of Weierstrass semigroups of certain m-tuples on the norm-trace function field
Minimal generating sets of Weierstrass semigroups of certain m-tuples on the norm-trace function field Gretchen L. Matthews and Justin D. Peachey Abstract. The norm-trace function field is a generalization
More informationAppendix I EFFICIENT DERIVATION OF IMPEDANCE LOCI. 1.1 Adaptive Sampling Scheme
Appendix I EFFICIENT DERIVATION OF IMPEDANCE LOCI 1.1 Adaptive Sampling Scheme For a given system configuration and operating conditions, the impedance locus, described in Section 4.3.5, provides complete
More informationANNUAL OF NAVIGATION 16/2010
ANNUAL OF NAVIGATION 16/2010 STANISŁAW KONATOWSKI, MARCIN DĄBROWSKI, ANDRZEJ PIENIĘŻNY Military University of Technology VEHICLE POSITIONING SYSTEM BASED ON GPS AND AUTONOMIC SENSORS ABSTRACT In many real
More informationAnalysis on the Properties of a Permutation Group
International Journal of Theoretical and Applied Mathematics 2017; 3(1): 19-24 http://www.sciencepublishinggroup.com/j/ijtam doi: 10.11648/j.ijtam.20170301.13 Analysis on the Properties of a Permutation
More informationPATTERN DESIGN FOR NON-GEODESIC WINDING TOROIDAL PRESSURE VESSELS
16 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS PATTERN DESIGN FOR NON-GEODESIC WINDING TOROIDAL PRESSURE VESSELS Lei Zu*, Qin-Xiang He*, Qing-Qing Ni** * Dept. of Engineering Mechanics, Xi'an University
More informationFiltering and Data Cutoff in FSI Retrievals
Filtering and Data Cutoff in FSI Retrievals C. Marquardt, Y. Andres, L. Butenko, A. von Engeln, A. Foresi, E. Heredia, R. Notarpietro, Y. Yoon Outline RO basics FSI-type retrievals Spherical asymmetry,
More informationPHYSICS 140A : STATISTICAL PHYSICS HW ASSIGNMENT #1 SOLUTIONS
PHYSICS 40A : STATISTICAL PHYSICS HW ASSIGNMENT # SOLUTIONS () The information entropy of a distribution {p n } is defined as S n p n log 2 p n, where n ranges over all possible configurations of a given
More informationMAS336 Computational Problem Solving. Problem 3: Eight Queens
MAS336 Computational Problem Solving Problem 3: Eight Queens Introduction Francis J. Wright, 2007 Topics: arrays, recursion, plotting, symmetry The problem is to find all the distinct ways of choosing
More informationSeparation and Extraction of Short-Circuit Power Consumption in Digital CMOS VLSI Circuits
Separation and Extraction of Short-Circuit Power Consumption in Digital CMOS VLSI Circuits Atila Alvandpour, Per Larsson-Edefors, and Christer Svensson Div of Electronic Devices, Dept of Physics, Linköping
More informationPermutations with short monotone subsequences
Permutations with short monotone subsequences Dan Romik Abstract We consider permutations of 1, 2,..., n 2 whose longest monotone subsequence is of length n and are therefore extremal for the Erdős-Szekeres
More informationMATH 105: Midterm #1 Practice Problems
Name: MATH 105: Midterm #1 Practice Problems 1. TRUE or FALSE, plus explanation. Give a full-word answer TRUE or FALSE. If the statement is true, explain why, using concepts and results from class to justify
More informationA Novel Adaptive Algorithm for
A Novel Adaptive Algorithm for Sinusoidal Interference Cancellation H. C. So Department of Electronic Engineering, City University of Hong Kong Tat Chee Avenue, Kowloon, Hong Kong August 11, 2005 Indexing
More informationHigh Precision GNSS for Mapping & GIS Professionals
High Precision GNSS for Mapping & GIS Professionals Agenda Address your needs for GNSS knowledge. GNSS Basics Satellite Ranging Fundamentals (Code $ Carrier) Differential Corrections (Post Processed $
More informationSolutions 2: Probability and Counting
Massachusetts Institute of Technology MITES 18 Physics III Solutions : Probability and Counting Due Tuesday July 3 at 11:59PM under Fernando Rendon s door Preface: The basic methods of probability and
More informationChapter 3 Solution to Problems
Chapter 3 Solution to Problems 1. The telemetry system of a geostationary communications satellite samples 100 sensors on the spacecraft in sequence. Each sample is transmitted to earth as an eight-bit
More informationA Signal Space Theory of Interferences Cancellation Systems
A Signal Space Theory of Interferences Cancellation Systems Osamu Ichiyoshi Human Network for Better 21 Century E-mail: osamu-ichiyoshi@muf.biglobe.ne.jp Abstract Interferences among signals from different
More information