Movatterモバイル変換


[0]ホーム

URL:


Jump to content
WikipediaThe Free Encyclopedia
Search

Spherical coordinate system

From Wikipedia, the free encyclopedia
(Redirected fromSpherical coordinates)

Coordinates comprising a distance and two angles
Thephysics convention. Spherical coordinates (r,θ,φ) as commonly used: (ISO80000-2:2019): radial distancer (slant distance to origin), polar angleθ (theta) (angle with respect to positive polar axis), and azimuthal angleφ (phi) (angle of rotation from the initial meridian plane).This is the convention followed in this article.

Inmathematics, aspherical coordinate system specifies a given point inthree-dimensional space by using a distance and two angles as its threecoordinates. These are

  • theradial distancer along the line connecting the point to a fixed point called theorigin;
  • thepolar angleθ between this radial line and a givenpolar axis;[a] and
  • theazimuthal angleφ, which is theangle of rotation of the radial line around the polar axis.[b]

(See graphic regarding the "physics convention".)

Once the radius is fixed, the three coordinates (r,θ,φ), known as a 3-tuple, provide a coordinate system on asphere, typically called thespherical polar coordinates.Theplane passing through the origin andperpendicular to the polar axis (where the polar angle is aright angle) is called thereference plane (sometimesfundamental plane).

Terminology

[edit]
The physics convention is followed in this article; (See both graphics re "physics convention" and re "mathematics convention".)

The radial distance from the fixed point of origin is also called theradius, orradial line, orradial coordinate. The polar angle may be calledinclination angle,zenith angle,normal angle, or thecolatitude. The user may choose to replace the inclination angle by itscomplement, theelevation angle (oraltitude angle), measured upward between the reference plane and the radial line—i.e., from the reference plane upward (towards to the positive z-axis) to the radial line. Thedepression angle is the negative of the elevation angle.(See graphic re the "physics convention"—not "mathematics convention".)

Both the use of symbols and the naming order of tuple coordinates differ among the several sources and disciplines. This article will use the ISO convention[1] frequently encountered inphysics, where the naming tuple gives the order as: radial distance, polar angle, azimuthal angle, or(r,θ,φ){\displaystyle (r,\theta ,\varphi )}. (See graphic re the "physics convention".) In contrast, the conventions in many mathematics books and texts give the naming order differently as: radial distance, "azimuthal angle", "polar angle", and(ρ,θ,φ){\displaystyle (\rho ,\theta ,\varphi )} or(r,θ,φ){\displaystyle (r,\theta ,\varphi )}—which switches the uses and meanings of symbols θ and φ. Other conventions may also be used, such asr for a radius from thez-axis that is not from the point of origin. Particular care must be taken to check the meaning of the symbols.

Themathematics convention. Spherical coordinates(r,θ,φ) as typically used: radial distancer, azimuthal angleθ, and polar angleφ. +The meanings ofθ andφ have been swapped—compared to thephysics convention. The 'south'-direction x-axis is depicted but the 'north'-direction x-axis is not. (As in physics,ρ (rho) is often used instead ofr to avoid confusion with the valuer in cylindrical and 2D polar coordinates.)

According to the conventions ofgeographical coordinate systems, positions are measured by latitude, longitude, and height (altitude). There are a number ofcelestial coordinate systems based on differentfundamental planes and with different terms for the various coordinates. The spherical coordinate systems used in mathematics normally useradians rather thandegrees; (note 90 degrees equalsπ2 radians). And these systems of themathematics convention may measure the azimuthal anglecounterclockwise (i.e., from the south directionx-axis, or 180°, towards the east directiony-axis, or +90°)—rather than measureclockwise (i.e., from the north direction x-axis, or 0°, towards the east direction y-axis, or +90°), as done in thehorizontal coordinate system.[2](See graphic re "mathematics convention".)

The spherical coordinate system of thephysics convention can be seen as a generalization of thepolar coordinate system inthree-dimensional space. It can be further extended to higher-dimensional spaces, and is then referred to as ahyperspherical coordinate system.

Definition

[edit]

To define a spherical coordinate system, one must designate anorigin point in space,O, and two orthogonal directions: thezenith reference direction and theazimuth reference direction. These choices determine a reference plane that is typically defined as containing the point of origin and thex– and y–axes, either of which may be designated as theazimuth reference direction. The reference plane is perpendicular (orthogonal) to the zenith direction, and typically is designated "horizontal" to the zenith direction's "vertical". The spherical coordinates of a pointP then are defined as follows:

  • Theradius orradial distance is theEuclidean distance from the originO toP.
  • Theinclination (orpolar angle) is the signed angle from the zenith reference direction to the line segmentOP. (Elevation may be used as the polar angle instead ofinclination; see below.)
  • Theazimuth (orazimuthal angle) is the signed angle measured from theazimuth reference direction to the orthogonal projection of the radial line segmentOP on the reference plane.

The sign of the azimuth is determined by designating the rotation that is thepositive sense of turning about the zenith. This choice is arbitrary, and is part of the coordinate system definition. (If the inclination is either zero or 180 degrees (=π radians), the azimuth is arbitrary. If the radius is zero, both azimuth and inclination are arbitrary.)

Theelevation is the signed angle from the x-y reference plane to the radial line segmentOP, where positive angles are designated as upward, towards the zenith reference.Elevation is 90 degrees (=π/2 radians)minus inclination. Thus, if the inclination is 60 degrees (=π/3 radians), then the elevation is 30 degrees (=π/6 radians).

Inlinear algebra, thevector from the originO to the pointP is often called theposition vector ofP.

Conventions

[edit]

Several different conventions exist for representing spherical coordinates and prescribing the naming order of their symbols. The 3-tuple number set(r,θ,φ){\displaystyle (r,\theta ,\varphi )} denotes radial distance, the polar angle—"inclination", or as the alternative, "elevation"—and the azimuthal angle. It is the common practice within the physics convention, as specified byISO standard80000-2:2019, and earlier inISO 31-11 (1992).

As stated above, this article describes the ISO "physics convention"—unless otherwise noted.

However, some authors (including mathematicians) use the symbolρ (rho) for radius, or radial distance,φ for inclination (or elevation) andθ for azimuth—while others keep the use ofr for the radius; all which "provides a logical extension of the usual polar coordinates notation".[3] As to order, some authors list the azimuthbefore the inclination (or the elevation) angle. Some combinations of these choices result in aleft-handed coordinate system. The standard "physics convention" 3-tuple set(r,θ,φ){\displaystyle (r,\theta ,\varphi )} conflicts with the usual notation for two-dimensionalpolar coordinates and three-dimensionalcylindrical coordinates, whereθ is often used for the azimuth.[3]

Angles are typically measured indegrees (°) or inradians (rad), where 360° = 2π rad. The use of degrees is most common in geography, astronomy, and engineering, where radians are commonly used in mathematics andtheoretical physics. The unit for radial distance is usually determined by the context, as occurs in applications of the 'unit sphere', seeapplications.

When the system is used to designate physical three-space, it is customary to assign positive to azimuth angles measured in thecounterclockwise sense from the reference direction on the reference plane—as seen from the "zenith" side of the plane. This convention is used in particular for geographical coordinates, where the "zenith" direction isnorth and the positive azimuth (longitude) angles are measured eastwards from someprime meridian.

Major conventions
coordinates set ordercorresponding local geographical directions
(Z,X,Y)
right/left-handed
(r,θinc,φaz,right)(U,S,E)right
(r,φaz,right,θel)(U,E,N)right
(r,θel,φaz,right)(U,N,E)left

Note:Easting (E), Northing (N), Upwardness (U). In the case of(U,S,E) the localazimuth angle would be measuredcounterclockwise fromS toE.

Unique coordinates

[edit]

Any spherical coordinate triplet (or tuple)(r,θ,φ){\displaystyle (r,\theta ,\varphi )} specifies a single point of three-dimensional space. On the reverse view, any single point has infinitely many equivalent spherical coordinates. That is, the user can add or subtract any number of full turns to the angular measures without changing the angles themselves, and therefore without changing the point. It is convenient in many contexts to use negative radial distances, the convention being(r,θ,φ){\displaystyle (-r,\theta ,\varphi )}, which is equivalent to(r,θ+180,φ){\displaystyle (r,\theta {+}180^{\circ },\varphi )} or(r,90θ,φ+180){\displaystyle (r,90^{\circ }{-}\theta ,\varphi {+}180^{\circ })} for anyr,θ, andφ. Moreover,(r,θ,φ){\displaystyle (r,-\theta ,\varphi )} is equivalent to(r,θ,φ+180){\displaystyle (r,\theta ,\varphi {+}180^{\circ })}.

When necessary to define a unique set of spherical coordinates for each point, the user must restrict therange, aka interval, of each coordinate. A common choice is:

  • radial distance:r ≥ 0,
  • polar angle:0° ≤θ ≤ 180°, or0 rad ≤θπ rad,
  • azimuth :0° ≤φ < 360°, or0 rad ≤φ < 2π rad.

But instead of the interval[0°, 360°), the azimuthφ is typically restricted to thehalf-open interval(−180°, +180°], or(−π, +π ] radians, which is the standard convention for geographic longitude.

For the polar angleθ, the range (interval) for inclination is[0°, 180°], which is equivalent to elevation range (interval)[−90°, +90°]. In geography, the latitude is the elevation.

Even with these restrictions, if the polar angle (inclination) is 0° or 180°—elevation is −90° or +90°—then the azimuth angle is arbitrary; and ifr is zero, both azimuth and polar angles are arbitrary. To define the coordinates as unique, the user can assert the convention that (in these cases) the arbitrary coordinates are set to zero.

Plotting

[edit]

To plot any dot from its spherical coordinates(r,θ,φ), whereθ is inclination, the user would: mover units from the origin in the zenith reference direction (z-axis); then rotate by the amount of the azimuth angle (φ) about the originfrom the designatedazimuth reference direction, (i.e., either the x- or y-axis, seeDefinition, above); and then rotatefrom the z-axis by the amount of theθ angle.

Applications

[edit]
In themathematics convention: A globe showing aunit sphere, withtuple coordinates of pointP (red): its radial distancer (red, not labeled); its azimuthal angleθ (not labeled); and its polar angle ofinclinationφ (not labeled). The radial distance upward along thezenith–axis from the point of origin to the surface of the sphere is assigned the value unity, or 1. + In this image,r appears to equal 4/6, or .67, (of unity); i.e., four of the six 'nested shells' to the surface. The azimuth angleθ appears to equal positive 90°, as rotatedcounterclockwise from the azimuth-reference x–axis; and the inclinationφ appears to equal 30°, as rotated from the zenith–axis. (Note the 'full' rotation, or inclination, from the zenith–axis to the y–axis is 90°).

Just as the two-dimensionalCartesian coordinate system is useful—has a wide set of applications—on a planar surface, a two-dimensional spherical coordinate system is useful on the surface of a sphere. For example, one sphere that is described inCartesian coordinates with the equationx2 +y2 +z2 =c2 can be described inspherical coordinates by the simple equationr =c. (In this system—shown here in the mathematics convention—the sphere is adapted as aunit sphere, where the radius is set to unity and then can generally be ignored, see graphic.)

This (unit sphere) simplification is also useful when dealing with objects such asrotational matrices. Spherical coordinates are also useful in analyzing systems that have some degree of symmetry about a point, including:volume integrals inside a sphere; the potential energy field surrounding a concentrated mass or charge; or global weather simulation in a planet's atmosphere.

The output pattern of the industrialloudspeaker shown here uses spherical polar plots taken at six frequencies

Three dimensional modeling ofloudspeaker output patterns can be used to predict their performance. A number of polar plots are required, taken at a wide selection of frequencies, as the pattern changes greatly with frequency. Polar plots help to show that many loudspeakers tend toward omnidirectionality at lower frequencies.

An important application of spherical coordinates provides for theseparation of variables in twopartial differential equations—theLaplace and theHelmholtz equations—that arise in many physical problems. The angular portions of the solutions to such equations take the form ofspherical harmonics. Another application isergonomic design, wherer is the arm length of a stationary person and the angles describe the direction of the arm as it reaches out. The spherical coordinate system is also commonly used in 3Dgame development to rotate the camera around the player's position[4]

In geography

[edit]
Main article:Geographic coordinate system
See also:ECEF

Instead of inclination, thegeographic coordinate system uses elevation angle (orlatitude), in the range (akadomain)−90° ≤φ ≤ 90° and rotated north from theequator plane. Latitude (i.e.,the angle of latitude) may be eithergeocentric latitude, measured (rotated) from the Earth's center—and designated variously byψ,q,φ′,φc,φg—orgeodetic latitude, measured (rotated) from the observer'slocal vertical, and typically designatedφ.The polar angle (inclination), which is 90° minus the latitude and ranges from 0 to 180°, is calledcolatitude in geography.

The azimuth angle (orlongitude) of a given position on Earth, commonly denoted byλ, is measured in degrees east or west from some conventional referencemeridian (most commonly theIERS Reference Meridian); thus its domain (or range) is−180° ≤λ ≤ 180° and a given reading is typically designated "East" or "West". For positions on theEarth or other solidcelestial body, the reference plane is usually taken to be the plane perpendicular to theaxis of rotation.

Instead of the radial distancer geographers commonly usealtitude above or below some local reference surface (vertical datum), which, for example, may be themean sea level. When needed, the radial distance can be computed from the altitude by adding theradius of Earth, which is approximately 6,360 ± 11 km (3,952 ± 7 miles).

However, modern geographical coordinate systems are quite complex, and the positions implied by these simple formulae may be inaccurate by several kilometers. The precise standard meanings oflatitude, longitude andaltitude are currently defined by theWorld Geodetic System (WGS), and take into account the flattening of the Earth at the poles (about 21 km or 13 miles) and many other details.

Planetary coordinate systems use formulations analogous to the geographic coordinate system.

In astronomy

[edit]

A series ofastronomical coordinate systems are used to measure the elevation angle from severalfundamental planes. These reference planes include: the observer'shorizon, thegalactic equator (defined by the rotation of theMilky Way), thecelestial equator (defined byEarth's rotation), the plane of theecliptic (defined by Earth's orbit around theSun), and the plane of the earthterminator (normal to the instantaneous direction to theSun).

Coordinate system conversions

[edit]
See also:List of common coordinate transformations § To spherical coordinates

As the spherical coordinate system is only one of many three-dimensional coordinate systems, there exist equations for converting coordinates between the spherical coordinate system and others.

Cartesian coordinates

[edit]

The spherical coordinates of a point in the ISO convention (i.e. for physics: radiusr, inclinationθ, azimuthφ) can be obtained from itsCartesian coordinates(x,y,z) by the formulae

r=x2+y2+z2θ=arccoszx2+y2+z2=arccoszr={arctanx2+y2zif z>0π+arctanx2+y2zif z<0+π2if z=0 and x2+y20undefinedif x=y=z=0φ=sgn(y)arccosxx2+y2={arctan(yx)if x>0,arctan(yx)+πif x<0 and y0,arctan(yx)πif x<0 and y<0,+π2if x=0 and y>0,π2if x=0 and y<0,undefinedif x=0 and y=0.{\displaystyle {\begin{aligned}r&={\sqrt {x^{2}+y^{2}+z^{2}}}\\\theta &=\arccos {\frac {z}{\sqrt {x^{2}+y^{2}+z^{2}}}}=\arccos {\frac {z}{r}}={\begin{cases}\arctan {\frac {\sqrt {x^{2}+y^{2}}}{z}}&{\text{if }}z>0\\\pi +\arctan {\frac {\sqrt {x^{2}+y^{2}}}{z}}&{\text{if }}z<0\\+{\frac {\pi }{2}}&{\text{if }}z=0{\text{ and }}{\sqrt {x^{2}+y^{2}}}\neq 0\\{\text{undefined}}&{\text{if }}x=y=z=0\\\end{cases}}\\\varphi &=\operatorname {sgn}(y)\arccos {\frac {x}{\sqrt {x^{2}+y^{2}}}}={\begin{cases}\arctan({\frac {y}{x}})&{\text{if }}x>0,\\\arctan({\frac {y}{x}})+\pi &{\text{if }}x<0{\text{ and }}y\geq 0,\\\arctan({\frac {y}{x}})-\pi &{\text{if }}x<0{\text{ and }}y<0,\\+{\frac {\pi }{2}}&{\text{if }}x=0{\text{ and }}y>0,\\-{\frac {\pi }{2}}&{\text{if }}x=0{\text{ and }}y<0,\\{\text{undefined}}&{\text{if }}x=0{\text{ and }}y=0.\end{cases}}\end{aligned}}}

Theinverse tangent denoted inφ = arctany/x must be suitably defined, taking into account the correct quadrant of(x,y), as done in the equations above. See the article onatan2.

Alternatively, the conversion can be considered as two sequentialrectangular to polar conversions: the first in the Cartesianxy plane from(x,y) to(R,φ), whereR is the projection ofr onto thexy-plane, and the second in the CartesianzR-plane from(z,R) to(r,θ). The correct quadrants forφ andθ are implied by the correctness of the planar rectangular to polar conversions.

These formulae assume that the two systems have the same origin, that the spherical reference plane is the Cartesianxy plane, thatθ is inclination from thez direction, and that the azimuth angles are measured from the Cartesianx axis (so that they axis hasφ = +90°). Ifθ measures elevation from the reference plane instead of inclination from the zenith the arccos above becomes an arcsin, and thecosθ andsinθ below become switched.

Conversely, the Cartesian coordinates may be retrieved from the spherical coordinates (radiusr,inclinationθ,azimuthφ), wherer[0, ∞),θ[0,π],φ[0, 2π), byx=rsinθcosφ,y=rsinθsinφ,z=rcosθ.{\displaystyle {\begin{aligned}x&=r\sin \theta \,\cos \varphi ,\\y&=r\sin \theta \,\sin \varphi ,\\z&=r\cos \theta .\end{aligned}}}

Cylindrical coordinates

[edit]
Main article:Cylindrical coordinate system

Cylindrical coordinates (axialradiusρ,azimuthφ,elevationz) may be converted into spherical coordinates (central radiusr,inclinationθ,azimuthφ), by the formulas

r=ρ2+z2,θ=arctanρz=arccoszρ2+z2,φ=φ.{\displaystyle {\begin{aligned}r&={\sqrt {\rho ^{2}+z^{2}}},\\\theta &=\arctan {\frac {\rho }{z}}=\arccos {\frac {z}{\sqrt {\rho ^{2}+z^{2}}}},\\\varphi &=\varphi .\end{aligned}}}

Conversely, the spherical coordinates may be converted into cylindrical coordinates by the formulae

ρ=rsinθ,φ=φ,z=rcosθ.{\displaystyle {\begin{aligned}\rho &=r\sin \theta ,\\\varphi &=\varphi ,\\z&=r\cos \theta .\end{aligned}}}

These formulae assume that the two systems have the same origin and same reference plane, measure the azimuth angleφ in the same senses from the same axis, and that the spherical angleθ is inclination from the cylindricalz axis.

Ellipsoidal coordinates

[edit]
See also:Ellipsoidal coordinates

It is also possible to deal with ellipsoids in Cartesian coordinates by using a modified version of the spherical coordinates.

Let P be an ellipsoid specified by the level set

ax2+by2+cz2=d.{\displaystyle ax^{2}+by^{2}+cz^{2}=d.}

The modified spherical coordinates of a point in P in the ISO convention (i.e. for physics:radiusr,inclinationθ,azimuthφ) can be obtained from itsCartesian coordinates(x,y,z) by the formulae

x=1arsinθcosφ,y=1brsinθsinφ,z=1crcosθ,r2=ax2+by2+cz2.{\displaystyle {\begin{aligned}x&={\frac {1}{\sqrt {a}}}r\sin \theta \,\cos \varphi ,\\y&={\frac {1}{\sqrt {b}}}r\sin \theta \,\sin \varphi ,\\z&={\frac {1}{\sqrt {c}}}r\cos \theta ,\\r^{2}&=ax^{2}+by^{2}+cz^{2}.\end{aligned}}}

An infinitesimal volume element is given by

dV=|(x,y,z)(r,θ,φ)|drdθdφ=1abcr2sinθdrdθdφ=1abcr2drdΩ.{\displaystyle \mathrm {d} V=\left|{\frac {\partial (x,y,z)}{\partial (r,\theta ,\varphi )}}\right|\,dr\,d\theta \,d\varphi ={\frac {1}{\sqrt {abc}}}r^{2}\sin \theta \,\mathrm {d} r\,\mathrm {d} \theta \,\mathrm {d} \varphi ={\frac {1}{\sqrt {abc}}}r^{2}\,\mathrm {d} r\,\mathrm {d} \Omega .}

The square-root factor comes from the property of thedeterminant that allows a constant to be pulled out from a column:

|kabckdefkghi|=k|abcdefghi|.{\displaystyle {\begin{vmatrix}ka&b&c\\kd&e&f\\kg&h&i\end{vmatrix}}=k{\begin{vmatrix}a&b&c\\d&e&f\\g&h&i\end{vmatrix}}.}

Integration and differentiation in spherical coordinates

[edit]
Unit vectors in spherical coordinates

The following equations (Iyanaga 1977) assume that the colatitudeθ is the inclination from the positivez axis, as in thephysics convention discussed.

Theline element for an infinitesimal displacement from(r,θ,φ) to(r + dr,θ + dθ,φ + dφ) isdr=drr^+rdθθ^+rsinθdφφ^,{\displaystyle \mathrm {d} \mathbf {r} =\mathrm {d} r\,{\hat {\mathbf {r} }}+r\,\mathrm {d} \theta \,{\hat {\boldsymbol {\theta }}}+r\sin {\theta }\,\mathrm {d} \varphi \,\mathbf {\hat {\boldsymbol {\varphi }}} ,}wherer^=sinθcosφx^+sinθsinφy^+cosθz^,θ^=cosθcosφx^+cosθsinφy^sinθz^,φ^=sinφx^+cosφy^{\displaystyle {\begin{aligned}{\hat {\mathbf {r} }}&=\sin \theta \cos \varphi \,{\hat {\mathbf {x} }}+\sin \theta \sin \varphi \,{\hat {\mathbf {y} }}+\cos \theta \,{\hat {\mathbf {z} }},\\{\hat {\boldsymbol {\theta }}}&=\cos \theta \cos \varphi \,{\hat {\mathbf {x} }}+\cos \theta \sin \varphi \,{\hat {\mathbf {y} }}-\sin \theta \,{\hat {\mathbf {z} }},\\{\hat {\boldsymbol {\varphi }}}&=-\sin \varphi \,{\hat {\mathbf {x} }}+\cos \varphi \,{\hat {\mathbf {y} }}\end{aligned}}}are the local orthogonalunit vectors in the directions of increasingr,θ, andφ, respectively,and,ŷ, and are the unit vectors in Cartesian coordinates. The linear transformation to this right-handed coordinate triplet is arotation matrix,R=(sinθcosφsinθsinφcosθcosθcosφcosθsinφsinθsinφcosφ0).{\displaystyle R={\begin{pmatrix}\sin \theta \cos \varphi &\sin \theta \sin \varphi &{\hphantom {-}}\cos \theta \\\cos \theta \cos \varphi &\cos \theta \sin \varphi &-\sin \theta \\-\sin \varphi &\cos \varphi &{\hphantom {-}}0\end{pmatrix}}.}

This gives the transformation from the Cartesian to the spherical, the other way around is given by its inverse.Note: the matrix is anorthogonal matrix, that is, its inverse is simply itstranspose.

The Cartesian unit vectors are thus related to the spherical unit vectors by:[x^y^z^]=[sinθcosφcosθcosφsinφsinθsinφcosθsinφcosφcosθsinθ0][r^θ^φ^]{\displaystyle {\begin{bmatrix}\mathbf {\hat {x}} \\\mathbf {\hat {y}} \\\mathbf {\hat {z}} \end{bmatrix}}={\begin{bmatrix}\sin \theta \cos \varphi &\cos \theta \cos \varphi &-\sin \varphi \\\sin \theta \sin \varphi &\cos \theta \sin \varphi &{\hphantom {-}}\cos \varphi \\\cos \theta &-\sin \theta &{\hphantom {-}}0\end{bmatrix}}{\begin{bmatrix}{\boldsymbol {\hat {r}}}\\{\boldsymbol {\hat {\theta }}}\\{\boldsymbol {\hat {\varphi }}}\end{bmatrix}}}

The general form of the formula to prove the differential line element, is[5]dr=irxidxi=i|rxi|rxi|rxi|dxi=i|rxi|dxix^i,{\displaystyle \mathrm {d} \mathbf {r} =\sum _{i}{\frac {\partial \mathbf {r} }{\partial x_{i}}}\,\mathrm {d} x_{i}=\sum _{i}\left|{\frac {\partial \mathbf {r} }{\partial x_{i}}}\right|{\frac {\frac {\partial \mathbf {r} }{\partial x_{i}}}{\left|{\frac {\partial \mathbf {r} }{\partial x_{i}}}\right|}}\,\mathrm {d} x_{i}=\sum _{i}\left|{\frac {\partial \mathbf {r} }{\partial x_{i}}}\right|\,\mathrm {d} x_{i}\,{\hat {\boldsymbol {x}}}_{i},}that is, the change inr{\displaystyle \mathbf {r} } is decomposed into individual changes corresponding to changes in the individual coordinates.

To apply this to the present case, one needs to calculate howr{\displaystyle \mathbf {r} } changes with each of the coordinates. In the conventions used,r=[rsinθcosφrsinθsinφrcosθ],x1=r,x2=θ,x3=φ.{\displaystyle \mathbf {r} ={\begin{bmatrix}r\sin \theta \,\cos \varphi \\r\sin \theta \,\sin \varphi \\r\cos \theta \end{bmatrix}},x_{1}=r,x_{2}=\theta ,x_{3}=\varphi .}

Thus,rr=[sinθcosφsinθsinφcosθ]=r^,rθ=[rcosθcosφrcosθsinφrsinθ]=rθ^,rφ=[rsinθsinφrsinθcosφ0]=rsinθφ^.{\displaystyle {\frac {\partial \mathbf {r} }{\partial r}}={\begin{bmatrix}\sin \theta \,\cos \varphi \\\sin \theta \,\sin \varphi \\\cos \theta \end{bmatrix}}=\mathbf {\hat {r}} ,\quad {\frac {\partial \mathbf {r} }{\partial \theta }}={\begin{bmatrix}r\cos \theta \,\cos \varphi \\r\cos \theta \,\sin \varphi \\-r\sin \theta \end{bmatrix}}=r\,{\hat {\boldsymbol {\theta }}},\quad {\frac {\partial \mathbf {r} }{\partial \varphi }}={\begin{bmatrix}-r\sin \theta \,\sin \varphi \\{\hphantom {-}}r\sin \theta \,\cos \varphi \\0\end{bmatrix}}=r\sin \theta \,\mathbf {\hat {\boldsymbol {\varphi }}} .}

The desired coefficients are the magnitudes of these vectors:[5]|rr|=1,|rθ|=r,|rφ|=rsinθ.{\displaystyle \left|{\frac {\partial \mathbf {r} }{\partial r}}\right|=1,\quad \left|{\frac {\partial \mathbf {r} }{\partial \theta }}\right|=r,\quad \left|{\frac {\partial \mathbf {r} }{\partial \varphi }}\right|=r\sin \theta .}

Thesurface element spanning fromθ toθ + dθ andφ toφ + dφ on a spherical surface at (constant) radiusr is thendSr=rθ×rφdθdφ=|rθ^×rsinθφ^|dθdφ=r2sinθdθdφ .{\displaystyle \mathrm {d} S_{r}=\left\|{\frac {\partial {\mathbf {r} }}{\partial \theta }}\times {\frac {\partial {\mathbf {r} }}{\partial \varphi }}\right\|\mathrm {d} \theta \,\mathrm {d} \varphi =\left|r{\hat {\boldsymbol {\theta }}}\times r\sin \theta {\boldsymbol {\hat {\varphi }}}\right|\mathrm {d} \theta \,\mathrm {d} \varphi =r^{2}\sin \theta \,\mathrm {d} \theta \,\mathrm {d} \varphi ~.}

Thus the differentialsolid angle isdΩ=dSrr2=sinθdθdφ.{\displaystyle \mathrm {d} \Omega ={\frac {\mathrm {d} S_{r}}{r^{2}}}=\sin \theta \,\mathrm {d} \theta \,\mathrm {d} \varphi .}

The surface element in a surface of polar angleθ constant (a cone with vertex at the origin) isdSθ=rsinθdφdr.{\displaystyle \mathrm {d} S_{\theta }=r\sin \theta \,\mathrm {d} \varphi \,\mathrm {d} r.}

The surface element in a surface of azimuthφ constant (a vertical half-plane) isdSφ=rdrdθ.{\displaystyle \mathrm {d} S_{\varphi }=r\,\mathrm {d} r\,\mathrm {d} \theta .}

Thevolume element spanning fromr tor + dr,θ toθ + dθ, andφ toφ + dφ is specified by thedeterminant of theJacobian matrix ofpartial derivatives,J=(x,y,z)(r,θ,φ)=(sinθcosφrcosθcosφrsinθsinφsinθsinφrcosθsinφrsinθcosφcosθrsinθ0),{\displaystyle J={\frac {\partial (x,y,z)}{\partial (r,\theta ,\varphi )}}={\begin{pmatrix}\sin \theta \cos \varphi &r\cos \theta \cos \varphi &-r\sin \theta \sin \varphi \\\sin \theta \sin \varphi &r\cos \theta \sin \varphi &{\hphantom {-}}r\sin \theta \cos \varphi \\\cos \theta &-r\sin \theta &{\hphantom {-}}0\end{pmatrix}},}namelydV=|(x,y,z)(r,θ,φ)|drdθdφ=r2sinθdrdθdφ=r2drdΩ .{\displaystyle \mathrm {d} V=\left|{\frac {\partial (x,y,z)}{\partial (r,\theta ,\varphi )}}\right|\,\mathrm {d} r\,\mathrm {d} \theta \,\mathrm {d} \varphi =r^{2}\sin \theta \,\mathrm {d} r\,\mathrm {d} \theta \,\mathrm {d} \varphi =r^{2}\,\mathrm {d} r\,\mathrm {d} \Omega ~.}

Thus, for example, a functionf(r,θ,φ) can be integrated over every point inR3 by thetriple integral02π0π0f(r,θ,φ)r2sinθdrdθdφ .{\displaystyle \int \limits _{0}^{2\pi }\int \limits _{0}^{\pi }\int \limits _{0}^{\infty }f(r,\theta ,\varphi )r^{2}\sin \theta \,\mathrm {d} r\,\mathrm {d} \theta \,\mathrm {d} \varphi ~.}

Thedel operator in this system leads to the following expressions for thegradient andLaplacian for scalar fields,f=frr^+1rfθθ^+1rsinθfφφ^,2f=1r2r(r2fr)+1r2sinθθ(sinθfθ)+1r2sin2θ2fφ2=(2r2+2rr)f+1r2sinθθ(sinθθ)f+1r2sin2θ2φ2f ,{\displaystyle {\begin{aligned}\nabla f&={\partial f \over \partial r}{\hat {\mathbf {r} }}+{1 \over r}{\partial f \over \partial \theta }{\hat {\boldsymbol {\theta }}}+{1 \over r\sin \theta }{\partial f \over \partial \varphi }{\hat {\boldsymbol {\varphi }}},\\[8pt]\nabla ^{2}f&={1 \over r^{2}}{\partial \over \partial r}\left(r^{2}{\partial f \over \partial r}\right)+{1 \over r^{2}\sin \theta }{\partial \over \partial \theta }\left(\sin \theta {\partial f \over \partial \theta }\right)+{1 \over r^{2}\sin ^{2}\theta }{\partial ^{2}f \over \partial \varphi ^{2}}\\[8pt]&=\left({\frac {\partial ^{2}}{\partial r^{2}}}+{\frac {2}{r}}{\frac {\partial }{\partial r}}\right)f+{1 \over r^{2}\sin \theta }{\partial \over \partial \theta }\left(\sin \theta {\frac {\partial }{\partial \theta }}\right)f+{\frac {1}{r^{2}\sin ^{2}\theta }}{\frac {\partial ^{2}}{\partial \varphi ^{2}}}f~,\\[8pt]\end{aligned}}}And it leads to the following expressions for thedivergence andcurl ofvector fields,

A=1r2r(r2Ar)+1rsinθθ(sinθAθ)+1rsinθAφφ,{\displaystyle \nabla \cdot \mathbf {A} ={\frac {1}{r^{2}}}{\partial \over \partial r}\left(r^{2}A_{r}\right)+{\frac {1}{r\sin \theta }}{\partial \over \partial \theta }\left(\sin \theta A_{\theta }\right)+{\frac {1}{r\sin \theta }}{\partial A_{\varphi } \over \partial \varphi },}×A=1rsinθ[θ(Aφsinθ)Aθφ]r^+1r[1sinθArφr(rAφ)]θ^+1r[r(rAθ)Arθ]φ^,{\displaystyle {\begin{aligned}\nabla \times \mathbf {A} ={}&{\frac {1}{r\sin \theta }}\left[{\partial \over \partial \theta }\left(A_{\varphi }\sin \theta \right)-{\partial A_{\theta } \over \partial \varphi }\right]{\hat {\mathbf {r} }}\\[4pt]&{}+{\frac {1}{r}}\left[{1 \over \sin \theta }{\partial A_{r} \over \partial \varphi }-{\partial \over \partial r}\left(rA_{\varphi }\right)\right]{\hat {\boldsymbol {\theta }}}\\[4pt]&{}+{\frac {1}{r}}\left[{\partial \over \partial r}\left(rA_{\theta }\right)-{\partial A_{r} \over \partial \theta }\right]{\hat {\boldsymbol {\varphi }}},\end{aligned}}}

Further, the inverse Jacobian in Cartesian coordinates isJ1=(xryrzrxzr2x2+y2yzr2x2+y2(x2+y2)r2x2+y2yx2+y2xx2+y20).{\displaystyle J^{-1}={\begin{pmatrix}{\dfrac {x}{r}}&{\dfrac {y}{r}}&{\dfrac {z}{r}}\\\\{\dfrac {xz}{r^{2}{\sqrt {x^{2}+y^{2}}}}}&{\dfrac {yz}{r^{2}{\sqrt {x^{2}+y^{2}}}}}&{\dfrac {-\left(x^{2}+y^{2}\right)}{r^{2}{\sqrt {x^{2}+y^{2}}}}}\\\\{\dfrac {-y}{x^{2}+y^{2}}}&{\dfrac {x}{x^{2}+y^{2}}}&0\end{pmatrix}}.}Themetric tensor in the spherical coordinate system isg=JTJ{\displaystyle g=J^{T}J}.

Distance in spherical coordinates

[edit]

In spherical coordinates, given two points withφ being the azimuthal coordinater=(r,θ,φ),r=(r,θ,φ){\displaystyle {\begin{aligned}{\mathbf {r} }&=(r,\theta ,\varphi ),\\{\mathbf {r} '}&=(r',\theta ',\varphi ')\end{aligned}}}The distance between the two points can be expressed as[6]D=r2+r22rr(sinθsinθcos(φφ)+cosθcosθ){\displaystyle {\begin{aligned}{\mathbf {D} }&={\sqrt {r^{2}+r'^{2}-2rr'(\sin {\theta }\sin {\theta '}\cos {(\varphi -\varphi ')}+\cos {\theta }\cos {\theta '})}}\end{aligned}}}

Kinematics

[edit]

In spherical coordinates, the position of a point or particle (although better written as atriple(r,θ,φ){\displaystyle (r,\theta ,\varphi )}) can be written as[7]r=rr^.{\displaystyle \mathbf {r} =r\mathbf {\hat {r}} .}Its velocity is then[7]v=drdt=r˙r^+rθ˙θ^+rφ˙sinθφ^{\displaystyle \mathbf {v} ={\frac {\mathrm {d} \mathbf {r} }{\mathrm {d} t}}={\dot {r}}\mathbf {\hat {r}} +r\,{\dot {\theta }}\,{\hat {\boldsymbol {\theta }}}+r\,{\dot {\varphi }}\sin \theta \,\mathbf {\hat {\boldsymbol {\varphi }}} }and its acceleration is[7]a=dvdt=+(r¨rθ˙2rφ˙2sin2θ)r^+(rθ¨+2r˙θ˙rφ˙2sinθcosθ)θ^+(rφ¨sinθ+2r˙φ˙sinθ+2rθ˙φ˙cosθ)φ^{\displaystyle {\begin{aligned}\mathbf {a} ={}&{\frac {\mathrm {d} \mathbf {v} }{\mathrm {d} t}}\\[1ex]={}&{\hphantom {+}}\;\left({\ddot {r}}-r\,{\dot {\theta }}^{2}-r\,{\dot {\varphi }}^{2}\sin ^{2}\theta \right)\mathbf {\hat {r}} \\&{}+\left(r\,{\ddot {\theta }}+2{\dot {r}}\,{\dot {\theta }}-r\,{\dot {\varphi }}^{2}\sin \theta \cos \theta \right){\hat {\boldsymbol {\theta }}}\\&{}+\left(r{\ddot {\varphi }}\,\sin \theta +2{\dot {r}}\,{\dot {\varphi }}\,\sin \theta +2r\,{\dot {\theta }}\,{\dot {\varphi }}\,\cos \theta \right){\hat {\boldsymbol {\varphi }}}\end{aligned}}}

The angular momentum isL=r×p=r×mv=mr2(φ˙sinθθ^+θ˙φ^){\displaystyle \mathbf {L} =\mathbf {r} \times \mathbf {p} =\mathbf {r} \times m\mathbf {v} =mr^{2}\left(-{\dot {\varphi }}\sin \theta \,\mathbf {\hat {\boldsymbol {\theta }}} +{\dot {\theta }}\,{\hat {\boldsymbol {\varphi }}}\right)}Wherem{\displaystyle m} is mass. In the case of a constantφ or elseθ =π/2, this reduces tovector calculus in polar coordinates.

The corresponding angular momentum operator then follows from the phase-space reformulation of the above,L=i r×=i(θ^sin(θ)ϕϕ^θ).{\displaystyle \mathbf {L} =-i\hbar ~\mathbf {r} \times \nabla =i\hbar \left({\frac {\hat {\boldsymbol {\theta }}}{\sin(\theta )}}{\frac {\partial }{\partial \phi }}-{\hat {\boldsymbol {\phi }}}{\frac {\partial }{\partial \theta }}\right).}

The torque is given as[7]τ=dLdt=r×F=m(2rr˙φ˙sinθ+r2φ¨sinθ+2r2θ˙φ˙cosθ)θ^+m(r2θ¨+2rr˙θ˙r2φ˙2sinθcosθ)φ^{\displaystyle \mathbf {\tau } ={\frac {\mathrm {d} \mathbf {L} }{\mathrm {d} t}}=\mathbf {r} \times \mathbf {F} =-m\left(2r{\dot {r}}{\dot {\varphi }}\sin \theta +r^{2}{\ddot {\varphi }}\sin {\theta }+2r^{2}{\dot {\theta }}{\dot {\varphi }}\cos {\theta }\right){\hat {\boldsymbol {\theta }}}+m\left(r^{2}{\ddot {\theta }}+2r{\dot {r}}{\dot {\theta }}-r^{2}{\dot {\varphi }}^{2}\sin \theta \cos \theta \right){\hat {\boldsymbol {\varphi }}}}

Thekinetic energy is given as[7]Ek=12m[(r˙)2+(rθ˙)2+(rφ˙sinθ)2]{\displaystyle E_{k}={\frac {1}{2}}m\left[\left({\dot {r}}\right)^{2}+\left(r{\dot {\theta }}\right)^{2}+\left(r{\dot {\varphi }}\sin \theta \right)^{2}\right]}

See also

[edit]

Notes

[edit]
  1. ^Anoriented line, so the polar angle is anoriented angle reckoned from the polar axis maindirection, not its opposite direction.
  2. ^If the polar axis is made to coincide with positivez-axis, the azimuthal angleφ may be calculated as the angle between either of thex-axis ory-axis and theorthogonal projection of the radial line onto the referencex-y-plane — which isorthogonal to thez-axis and passes through the fixed point of origin, completing a three-dimensionalCartesian coordinate system.

References

[edit]
  1. ^"ISO 80000-2:2019 Quantities and units – Part 2: Mathematics".ISO. 19 May 2020. pp. 20–21. Item no. 2-17.3. Retrieved2020-08-12.
  2. ^Duffett-Smith, P and Zwart, J, p. 34.
  3. ^abEric W. Weisstein (2005-10-26)."Spherical Coordinates".MathWorld. Retrieved2010-01-15.
  4. ^"Video Game Math: Polar and Spherical Notation".Academy of Interactive Entertainment (AIE). Retrieved2022-02-16.
  5. ^ab"Line element (dl) in spherical coordinates derivation/diagram".Stack Exchange. October 21, 2011.
  6. ^"Distance between two points in spherical coordinates".
  7. ^abcdeReed, Bruce Cameron (2019).Keplerian ellipses : the physics of the gravitational two-body problem. Morgan & Claypool Publishers, Institute of Physics. San Rafael [California] (40 Oak Drive, San Rafael, CA, 94903, US).ISBN 978-1-64327-470-6.OCLC 1104053368.{{cite book}}: CS1 maint: location (link) CS1 maint: location missing publisher (link)

Bibliography

[edit]

External links

[edit]
Two dimensional
Three dimensional
Retrieved from "https://en.wikipedia.org/w/index.php?title=Spherical_coordinate_system&oldid=1303772995"
Categories:
Hidden categories:

[8]ページ先頭

©2009-2025 Movatter.jp