Polar coordinate system
In mathematics, the polar coordinate system specifies a given point in a plane by using a distance and an angle as its two coordinates. These are
- the point's distance from a reference point called the pole, and
- the point's direction from the pole relative to the direction of the polar axis, a ray drawn from the pole.
The distance from the pole is called the radial coordinate, radial distance or simply radius, and the angle is called the angular coordinate, polar angle, or azimuth.<ref name="brown">Template:Cite book</ref> The pole is analogous to the origin in a Cartesian coordinate system.
Polar coordinates are most appropriate in any context where the phenomenon being considered is inherently tied to direction and length from a center point in a plane, such as spirals. Planar physical systems with bodies moving around a central point, or phenomena originating from a central point, are often simpler and more intuitive to model using polar coordinates.
The polar coordinate system is extended to three dimensions in two ways: the cylindrical coordinate system adds a second distance coordinate, and the spherical coordinate system adds a second angular coordinate.
Grégoire de Saint-Vincent and Bonaventura Cavalieri independently introduced the system's concepts in the mid-17th century, though the actual term polar coordinates has been attributed to Gregorio Fontana in the 18th century. The initial motivation for introducing the polar system was the study of circular and orbital motion.
HistoryEdit
.jpg){kind=link}
The concepts of angle and radius were already used by ancient peoples of the first millennium BC. The Greek astronomer and astrologer Hipparchus (190–120 BC) created a table of chord functions giving the length of the chord for each angle, and there are references to his using polar coordinates in establishing stellar positions.<ref>{{#invoke:citation/CS1|citation |CitationClass=web }}</ref> In On Spirals, Archimedes describes the Archimedean spiral, a function whose radius depends on the angle. The Greek work, however, did not extend to a full coordinate system.
From the 8th century AD onward, astronomers developed methods for approximating and calculating the direction to Mecca (qibla)—and its distance—from any location on the Earth.<ref>Template:Cite book</ref> From the 9th century onward they were using spherical trigonometry and map projection methods to determine these quantities accurately. The calculation is essentially the conversion of the equatorial polar coordinates of Mecca (i.e. its longitude and latitude) to its polar coordinates (i.e. its qibla and distance) relative to a system whose reference meridian is the great circle through the given location and the Earth's poles and whose polar axis is the line through the location and its antipodal point.<ref>King (2005, p. 169). The calculations were as accurate as could be achieved under the limitations imposed by their assumption that the Earth was a perfect sphere.</ref>
There are various accounts of the introduction of polar coordinates as part of a formal coordinate system. The full history of the subject is described in Harvard professor Julian Lowell Coolidge's Origin of Polar Coordinates.<ref name="coolidge">Template:Cite journal</ref> Grégoire de Saint-Vincent and Bonaventura Cavalieri independently introduced the concepts in the mid-seventeenth century. Saint-Vincent wrote about them privately in 1625 and published his work in 1647, while Cavalieri published his in 1635 with a corrected version appearing in 1653. Cavalieri first used polar coordinates to solve a problem relating to the area within an Archimedean spiral. Blaise Pascal subsequently used polar coordinates to calculate the length of parabolic arcs.
In Method of Fluxions (written 1671, published 1736), Sir Isaac Newton examined the transformations between polar coordinates, which he referred to as the "Seventh Manner; For Spirals", and nine other coordinate systems.<ref>Template:Cite journal</ref> In the journal Acta Eruditorum (1691), Jacob Bernoulli used a system with a point on a line, called the pole and polar axis respectively. Coordinates were specified by the distance from the pole and the angle from the polar axis. Bernoulli's work extended to finding the radius of curvature of curves expressed in these coordinates.
The actual term polar coordinates has been attributed to Gregorio Fontana and was used by 18th-century Italian writers. The term appeared in English in George Peacock's 1816 translation of Lacroix's Differential and Integral Calculus.<ref>{{#invoke:citation/CS1|citation |CitationClass=web }}</ref><ref>Template:Cite book</ref> Alexis Clairaut was the first to think of polar coordinates in three dimensions, and Leonhard Euler was the first to actually develop them.<ref name="coolidge" />
ConventionsEdit
The radial coordinate is often denoted by r or ρ, and the angular coordinate by φ, θ, or t. The angular coordinate is specified as φ by ISO standard 31-11, now 80000-2:2019. However, in mathematical literature the angle is often denoted by θ instead.
Angles in polar notation are generally expressed in either degrees or radians (2[[pi|Template:Pi]] rad being equal to 360°). Degrees are traditionally used in navigation, surveying, and many applied disciplines, while radians are more common in mathematics and mathematical physics.<ref>Template:Cite book</ref>
The angle φ is defined to start at 0° from a reference direction, and to increase for rotations in either clockwise (cw) or counterclockwise (ccw) orientation. For example, in mathematics, the reference direction is usually drawn as a ray from the pole horizontally to the right, and the polar angle increases to positive angles for ccw rotations, whereas in navigation (bearing, heading) the 0°-heading is drawn vertically upwards and the angle increases for cw rotations. The polar angles decrease towards negative values for rotations in the respectively opposite orientations.
Uniqueness of polar coordinatesEdit
Adding any number of full turns (360°) to the angular coordinate does not change the corresponding direction. Similarly, any polar coordinate is identical to the coordinate with the negative radial component and the opposite direction (adding 180° to the polar angle). Therefore, the same point (r, φ) can be expressed with an infinite number of different polar coordinates Template:Nowrap and Template:Nowrap, where n is an arbitrary integer.<ref>{{#invoke:citation/CS1|citation |CitationClass=web }}</ref> Moreover, the pole itself can be expressed as (0, φ) for any angle φ.<ref>Template:Cite book</ref>
Where a unique representation is needed for any point besides the pole, it is usual to limit r to positive numbers (Template:Nowrap) and φ to either the interval Template:Closed-open or the interval Template:Open-closed, which in radians are Template:Closed-open or Template:Open-closed.<ref>Template:Cite book</ref> Another convention, in reference to the usual codomain of the arctan function, is to allow for arbitrary nonzero real values of the radial component and restrict the polar angle to Template:Open-closed. In all cases a unique azimuth for the pole (r = 0) must be chosen, e.g., φ = 0.
Converting between polar and Cartesian coordinatesEdit
The polar coordinates r and φ can be converted to the Cartesian coordinates x and y by using the trigonometric functions sine and cosine:
<math display=block>\begin{align}
x &= r \cos \varphi, \\ y &= r \sin \varphi.
\end{align}</math>
The Cartesian coordinates x and y can be converted to polar coordinates r and φ with r ≥ 0 and φ in the interval (−Template:Pi, Template:Pi] by:<ref>Template:Cite book</ref> <math display=block> \begin{align} r &= \sqrt{x^2 + y^2} = \operatorname{hypot}(x,y)\\ \varphi &= \operatorname{atan2}(y, x), \end{align}</math> where hypot is the Pythagorean sum and atan2 is a common variation on the arctangent function defined as <math display=block>\operatorname{atan2}(y, x) = \begin{cases}
\arctan\left(\frac{y}{x}\right) & \mbox{if } x > 0\\
\arctan\left(\frac{y}{x}\right) + \pi & \mbox{if } x < 0 \mbox{ and } y \ge 0\\
\arctan\left(\frac{y}{x}\right) - \pi & \mbox{if } x < 0 \mbox{ and } y < 0\\
\frac{\pi}{2} & \mbox{if } x = 0 \mbox{ and } y > 0\\
-\frac{\pi}{2} & \mbox{if } x = 0 \mbox{ and } y < 0\\
\text{undefined} & \mbox{if } x = 0 \mbox{ and } y = 0.
\end{cases}</math>
If r is calculated first as above, then this formula for φ may be stated more simply using the arccosine function: <math display=block>\varphi = \begin{cases}
\arccos\left(\frac{x}{r}\right) & \mbox{if } y \ge 0 \mbox{ and } r \neq 0 \\
-\arccos\left(\frac{x}{r}\right) & \mbox{if } y < 0 \\
\text{undefined} & \mbox{if } r = 0.
\end{cases}</math>
Complex numbersEdit
Every complex number can be represented as a point in the complex plane, and can therefore be expressed by specifying either the point's Cartesian coordinates (called rectangular or Cartesian form) or the point's polar coordinates (called polar form).
In polar form, the distance and angle coordinates are often referred to as the number's magnitude and argument respectively. Two complex numbers can be multiplied by adding their arguments and multiplying their magnitudes.
The complex number z can be represented in rectangular form as <math display="block">z = x + iy</math> where i is the imaginary unit, or can alternatively be written in polar form as <math display="block">z = r(\cos\varphi + i\sin\varphi)</math> and from there, by Euler's formula,<ref>Template:Cite book</ref> as <math display="block">z = re^{i\varphi} = r \exp i \varphi. </math> where e is Euler's number, and φ, expressed in radians, is the principal value of the complex number function arg applied to x + iy. To convert between the rectangular and polar forms of a complex number, the conversion formulae given above can be used. Equivalent are the Template:Math and angle notations: <math display="block"> z = r \operatorname\mathrm{cis} \varphi = r \angle \varphi .</math>
For the operations of multiplication, division, exponentiation, and root extraction of complex numbers, it is generally much simpler to work with complex numbers expressed in polar form rather than rectangular form. From the laws of exponentiation:
- Multiplication
- <math>r_0 e^{i\varphi_0}\, r_1 e^{i\varphi_1} = r_0 r_1 e^{i\left(\varphi_0 + \varphi_1\right)} </math>
- Division
- <math>\frac{r_0 e^{i\varphi_0}}{r_1 e^{i\varphi_1}} = \frac{r_0}{r_1}e^{i(\varphi_0 - \varphi_1)} </math>
- Exponentiation (De Moivre's formula)
- <math>\left(re^{i\varphi}\right)^n = r^n e^{in\varphi} </math>
- Root Extraction (Principal root)
- <math>\sqrt[n]{re^{i\varphi}} = \sqrt[n]{r} e^{i\varphi \over n} </math>
Polar equation of a curveEdit
The equation defining a plane curve expressed in polar coordinates is known as a polar equation. In many cases, such an equation can simply be specified by defining r as a function of φ. The resulting curve then consists of points of the form (r(φ), φ) and can be regarded as the graph of the polar function r. Note that, in contrast to Cartesian coordinates, the independent variable φ is the second entry in the ordered pair.
Different forms of symmetry can be deduced from the equation of a polar function r:
- If Template:Math the curve will be symmetrical about the horizontal (0°/180°) ray;
- If Template:Math it will be symmetric about the vertical (90°/270°) ray:
- If Template:Math it will be rotationally symmetric by α clockwise and counterclockwise about the pole.
Because of the circular nature of the polar coordinate system, many curves can be described by a rather simple polar equation, whereas their Cartesian form is much more intricate. Among the best known of these curves are the polar rose, Archimedean spiral, lemniscate, limaçon, and cardioid.
For the circle, line, and polar rose below, it is understood that there are no restrictions on the domain and range of the curve.
CircleEdit
The general equation for a circle with a center at <math>(r_0, \gamma)</math> and radius a is <math display="block">r^2 - 2 r r_0 \cos(\varphi - \gamma) + r_0^2 = a^2.</math>
This can be simplified in various ways, to conform to more specific cases, such as the equation <math display="block">r(\varphi)=a </math> for a circle with a center at the pole and radius a.<ref name="ping">{{#invoke:citation/CS1|citation |CitationClass=web }}</ref>
When Template:Math or the origin lies on the circle, the equation becomes <math display="block">r = 2 a\cos(\varphi - \gamma).</math>
In the general case, the equation can be solved for Template:Math, giving <math display="block">r = r_0 \cos(\varphi - \gamma) + \sqrt{a^2 - r_0^2 \sin^2(\varphi - \gamma)}</math> The solution with a minus sign in front of the square root gives the same curve.
LineEdit
Radial lines (those running through the pole) are represented by the equation <math display="block">\varphi = \gamma,</math> where <math>\gamma</math> is the angle of elevation of the line; that is, <math>\varphi = \arctan m</math>, where <math>m</math> is the slope of the line in the Cartesian coordinate system. The non-radial line that crosses the radial line <math>\varphi = \gamma</math> perpendicularly at the point <math>(r_0, \gamma)</math> has the equation <math display="block">r(\varphi) = r_0 \sec(\varphi - \gamma).</math>
Otherwise stated <math>(r_0, \gamma)</math> is the point in which the tangent intersects the imaginary circle of radius <math>r_0</math>
Polar roseEdit
.svg){kind=link}
A polar rose is a mathematical curve that looks like a petaled flower, and that can be expressed as a simple polar equation, <math display="block">r(\varphi) = a\cos\left(k\varphi + \gamma_0\right)</math>
for any constant γ0 (including 0). If k is an integer, these equations will produce a k-petaled rose if k is odd, or a 2k-petaled rose if k is even. If k is rational, but not an integer, a rose-like shape may form but with overlapping petals. Note that these equations never define a rose with 2, 6, 10, 14, etc. petals. The variable a directly represents the length or amplitude of the petals of the rose, while k relates to their spatial frequency. The constant γ0 can be regarded as a phase angle.
Archimedean spiralEdit
The Archimedean spiral is a spiral discovered by Archimedes which can also be expressed as a simple polar equation. It is represented by the equation <math display="block">r(\varphi) = a + b\varphi. </math> Changing the parameter a will turn the spiral, while b controls the distance between the arms, which for a given spiral is always constant. The Archimedean spiral has two arms, one for Template:Math and one for Template:Math. The two arms are smoothly connected at the pole. If Template:Math, taking the mirror image of one arm across the 90°/270° line will yield the other arm. This curve is notable as one of the first curves, after the conic sections, to be described in a mathematical treatise, and as a prime example of a curve best defined by a polar equation.
Conic sectionsEdit
A conic section with one focus on the pole and the other somewhere on the 0° ray (so that the conic's major axis lies along the polar axis) is given by: <math display="block">r = { \ell\over {1 - e \cos \varphi}}</math> where e is the eccentricity and <math>\ell</math> is the semi-latus rectum (the perpendicular distance at a focus from the major axis to the curve). If Template:Nowrap, this equation defines a hyperbola; if Template:Math, it defines a parabola; and if Template:Math, it defines an ellipse. The special case Template:Math of the latter results in a circle of the radius <math>\ell</math>.
QuadratrixEdit
{{#invoke:Labelled list hatnote|labelledList|Main article|Main articles|Main page|Main pages}}
A quadratrix in the first quadrant (x, y) is a curve with y = ρ sin θ equal to the fraction of the quarter circle with radius r determined by the radius through the curve point. Since this fraction is <math>\frac{2 r \theta}{\pi}</math>, the curve is given by <math>\rho (\theta) = \frac{2 r \theta}{\pi \sin \theta}</math>.<ref>N.H. Lucas, P.J. Bunt & J.D Bedient (1976) Historical Roots of Elementary Mathematics, page 113</ref>
Intersection of two polar curvesEdit
The graphs of two polar functions <math>r = f(\theta)</math> and <math>r = g(\theta)</math> have possible intersections of three types:
- In the origin, if the equations <math>f(\theta) = 0</math> and <math>g(\theta) = 0</math> have at least one solution each.
- All the points <math>[g(\theta_i),\theta_i]</math> where <math>\theta_i</math> are solutions to the equation <math>f(\theta+2k\pi)=g(\theta)</math> where <math>k</math> is an integer.
- All the points <math>[g(\theta_i),\theta_i]</math> where <math>\theta_i</math> are solutions to the equation <math>f(\theta+(2k+1)\pi)=-g(\theta)</math> where <math>k</math> is an integer.
CalculusEdit
Calculus can be applied to equations expressed in polar coordinates.<ref>{{#invoke:citation/CS1|citation |CitationClass=web }}</ref><ref>{{#invoke:citation/CS1|citation |CitationClass=web }}</ref>
The angular coordinate φ is expressed in radians throughout this section, which is the conventional choice when doing calculus.
Differential calculusEdit
Using Template:Math and Template:Math, one can derive a relationship between derivatives in Cartesian and polar coordinates. For a given function, u(x,y), it follows that (by computing its total derivatives) or <math display="block">\begin{align}
r \frac{du}{dr}
&= r \frac{\partial u}{\partial x} \cos\varphi +
r \frac{\partial u}{\partial y} \sin\varphi
= x \frac{\partial u}{\partial x} + y \frac{\partial u}{\partial y}, \\[2pt]
\frac{du}{d\varphi}
&= - \frac{\partial u}{\partial x} r \sin\varphi +
\frac{\partial u}{\partial y} r \cos\varphi
= -y \frac{\partial u}{\partial x} + x \frac{\partial u}{\partial y}.
\end{align}</math>
Hence, we have the following formula: <math display="block">\begin{align}
r \frac{d}{dr} &= x \frac{\partial}{\partial x} + y \frac{\partial}{\partial y} \\[2pt]
\frac{d}{d\varphi} &= -y \frac{\partial}{\partial x} + x \frac{\partial}{\partial y}.
\end{align}</math>
Using the inverse coordinates transformation, an analogous reciprocal relationship can be derived between the derivatives. Given a function u(r,φ), it follows that <math display="block">\begin{align}
\frac{du}{dx}
&= \frac{\partial u}{\partial r}\frac{\partial r}{\partial x} +
\frac{\partial u}{\partial \varphi}\frac{\partial \varphi}{\partial x}, \\[2pt]
\frac{du}{dy}
&= \frac{\partial u}{\partial r}\frac{\partial r}{\partial y} +
\frac{\partial u}{\partial \varphi}\frac{\partial \varphi}{\partial y},
\end{align}</math> or <math display="block">\begin{align}
\frac{du}{dx}
&= \frac{\partial u}{\partial r}\frac{x}{\sqrt{x^2+y^2}} -
\frac{\partial u}{\partial \varphi}\frac{y}{x^2+y^2} \\[2pt]
&= \cos \varphi \frac{\partial u}{\partial r} -
\frac{1}{r} \sin\varphi \frac{\partial u}{\partial \varphi}, \\[2pt]
\frac{du}{dy}
&= \frac{\partial u}{\partial r}\frac{y}{\sqrt{x^2+y^2}} +
\frac{\partial u}{\partial \varphi}\frac{x}{x^2+y^2} \\[2pt]
&= \sin\varphi \frac{\partial u}{\partial r} +
\frac{1}{r} \cos\varphi \frac{\partial u}{\partial \varphi}.
\end{align}</math>
Hence, we have the following formulae: <math display="block">\begin{align}
\frac{d}{dx}
&= \cos \varphi \frac{\partial}{\partial r} -
\frac{1}{r} \sin\varphi \frac{\partial}{\partial \varphi} \\[2pt]
\frac{d}{dy}
&= \sin \varphi \frac{\partial}{\partial r} +
\frac{1}{r} \cos\varphi \frac{\partial}{\partial \varphi}.
\end{align}</math>
To find the Cartesian slope of the tangent line to a polar curve r(φ) at any given point, the curve is first expressed as a system of parametric equations. <math display="block">\begin{align}
x &= r(\varphi)\cos\varphi \\ y &= r(\varphi)\sin\varphi
\end{align}</math>
Differentiating both equations with respect to φ yields <math display="block">\begin{align}
\frac{dx}{d\varphi} &= r'(\varphi)\cos\varphi - r(\varphi)\sin\varphi \\[2pt]
\frac{dy}{d\varphi} &= r'(\varphi)\sin\varphi + r(\varphi)\cos\varphi.
\end{align}</math>
Dividing the second equation by the first yields the Cartesian slope of the tangent line to the curve at the point Template:Nowrap: <math display="block">\frac{dy}{dx} = \frac{r'(\varphi)\sin\varphi + r(\varphi)\cos\varphi}{r'(\varphi)\cos\varphi-r(\varphi)\sin\varphi}.</math>
For other useful formulas including divergence, gradient, and Laplacian in polar coordinates, see curvilinear coordinates.
Integral calculus (arc length)Edit
The arc length (length of a line segment) defined by a polar function is found by the integration over the curve r(φ). Let L denote this length along the curve starting from points A through to point B, where these points correspond to φ = a and φ = b such that Template:Math. The length of L is given by the following integral <math display="block">L = \int_a^b \sqrt{ \left[r(\varphi)\right]^2 + \left[ {\tfrac{dr(\varphi) }{ d\varphi }} \right] ^2 } d\varphi</math>
Integral calculus (area)Edit
Let R denote the region enclosed by a curve r(φ) and the rays φ = a and φ = b, where Template:Nowrap. Then, the area of R is <math display="block">\frac12\int_a^b \left[r(\varphi)\right]^2\, d\varphi.</math>
This result can be found as follows. First, the interval Template:Nowrap is divided into n subintervals, where n is some positive integer. Thus Δφ, the angle measure of each subinterval, is equal to Template:Math (the total angle measure of the interval), divided by n, the number of subintervals. For each subinterval i = 1, 2, ..., n, let φi be the midpoint of the subinterval, and construct a sector with the center at the pole, radius r(φi), central angle Δφ and arc length r(φi)Δφ. The area of each constructed sector is therefore equal to <math display="block">\left[r(\varphi_i)\right]^2 \pi \cdot \frac{\Delta \varphi}{2\pi} = \frac{1}{2}\left[r(\varphi_i)\right]^2 \Delta \varphi.</math> Hence, the total area of all of the sectors is <math display="block">\sum_{i=1}^n \tfrac12r(\varphi_i)^2\,\Delta\varphi.</math>
As the number of subintervals n is increased, the approximation of the area improves. Taking Template:Nowrap, the sum becomes the Riemann sum for the above integral.
A mechanical device that computes area integrals is the planimeter, which measures the area of plane figures by tracing them out: this replicates integration in polar coordinates by adding a joint so that the 2-element linkage effects Green's theorem, converting the quadratic polar integral to a linear integral.
GeneralizationEdit
Using Cartesian coordinates, an infinitesimal area element can be calculated as dA = dx dy. The substitution rule for multiple integrals states that, when using other coordinates, the Jacobian determinant of the coordinate conversion formula has to be considered: <math display="block">J = \det \frac{\partial(x, y)}{\partial(r, \varphi)}
= \begin{vmatrix}
\frac{\partial x}{\partial r} & \frac{\partial x}{\partial \varphi} \\[2pt]
\frac{\partial y}{\partial r} & \frac{\partial y}{\partial \varphi}
\end{vmatrix}
= \begin{vmatrix}
\cos\varphi & -r\sin\varphi \\
\sin\varphi & r\cos\varphi
\end{vmatrix}
= r\cos^2\varphi + r\sin^2\varphi = r.
</math>
Hence, an area element in polar coordinates can be written as <math display="block">dA = dx\,dy\ = J\,dr\,d\varphi = r\,dr\,d\varphi.</math>
Now, a function, that is given in polar coordinates, can be integrated as follows: <math display="block">\iint_R f(x, y)\, dA = \int_a^b \int_0^{r(\varphi)} f(r, \varphi)\,r\,dr\,d\varphi.</math>
Here, R is the same region as above, namely, the region enclosed by a curve r(φ) and the rays φ = a and φ = b. The formula for the area of R is retrieved by taking f identically equal to 1.
.svg){kind=link}
A more surprising application of this result yields the Gaussian integral: <math display="block">\int_{-\infty}^\infty e^{-x^2} \, dx = \sqrt\pi.</math>
Vector calculusEdit
Vector calculus can also be applied to polar coordinates. For a planar motion, let <math>\mathbf{r}</math> be the position vector Template:Math, with r and Template:Math depending on time t.
Template:AnchorWe define an orthonormal basis with three unit vectors: radial, transverse, and normal directions. The radial direction is defined by normalizing <math>\mathbf{r}</math>: <math display="block">\hat{\mathbf{r}} = (\cos(\varphi), \sin(\varphi))</math> Radial and velocity directions span the plane of the motion, whose normal direction is denoted <math>\hat{\mathbf{k}}</math>: <math display="block">\hat{\mathbf{k}} = \hat{\mathbf{v}} \times \hat{\mathbf{r}}.</math> The transverse direction is perpendicular to both radial and normal directions: <math display="block">\hat \boldsymbol \varphi = (-\sin(\varphi), \cos(\varphi)) = \hat{\mathbf{k}} \times \hat{\mathbf{r}} \ ,</math>
Then <math display="block">\begin{align}
\mathbf{r} &= (x,\ y) = r(\cos\varphi,\ \sin\varphi) = r \hat{\mathbf{r}}\ , \\[1.5ex]
\dot{\mathbf{r}} &= \left(\dot{x},\ \dot{y}\right) = \dot{r}(\cos\varphi,\ \sin\varphi) + r\dot{\varphi}(-\sin\varphi,\ \cos\varphi) = \dot{r}\hat{\mathbf{r}} + r\dot{\varphi}\hat{\boldsymbol{\varphi}}\ ,\\[1.5ex]
\ddot{\mathbf{r}} &= \left(\ddot{x},\ \ddot{y}\right) \\[1ex]
&= \ddot{r}(\cos\varphi,\ \sin\varphi) + 2\dot{r}\dot{\varphi}(-\sin\varphi,\ \cos\varphi) + r\ddot{\varphi}(-\sin\varphi,\ \cos\varphi) - r\dot{\varphi}^2(\cos\varphi,\ \sin\varphi) \\[1ex]
&= \left(\ddot{r} - r\dot{\varphi}^2\right) \hat{\mathbf{r}} + \left(r\ddot{\varphi} + 2\dot{r}\dot{\varphi}\right) \hat{\boldsymbol{\varphi}} \\[1ex]
&= \left(\ddot{r} - r\dot{\varphi}^2\right) \hat{\mathbf{r}} + \frac{1}{r}\; \frac{d}{dt} \left(r^2\dot{\varphi}\right) \hat{\boldsymbol{\varphi}}.
\end{align}</math>
This equation can be obtain by taking derivative of the function and derivatives of the unit basis vectors.
For a curve in 2D where the parameter is <math>\theta</math> the previous equations simplify to: <math display="block">\begin{aligned} \mathbf{r} &= r(\theta) \hat \mathbf{e}_r \\[1ex] \frac {d\mathbf{r}}{d\theta} &= \frac {dr} {d\theta} \hat \mathbf{e}_r + r \hat \mathbf{e}_\theta\\[1ex] \frac {d^2\mathbf{r}}{d\theta^2} &= \left(\frac {d^2 r} {d\theta^2}-r\right) \hat \mathbf{e}_r + \frac {dr} {d\theta} \hat \mathbf{e}_\theta \end{aligned}</math>
Centrifugal and Coriolis termsEdit
Template:See also Template:Multiple image
The term <math>r\dot\varphi^2</math> is sometimes referred to as the centripetal acceleration, and the term <math>2\dot r \dot\varphi</math> as the Coriolis acceleration. For example, see Shankar.<ref name="Shankar">Template:Cite book</ref>
Note: these terms, that appear when acceleration is expressed in polar coordinates, are a mathematical consequence of differentiation; they appear whenever polar coordinates are used. In planar particle dynamics these accelerations appear when setting up Newton's second law of motion in a rotating frame of reference. Here these extra terms are often called fictitious forces; fictitious because they are simply a result of a change in coordinate frame. That does not mean they do not exist, rather they exist only in the rotating frame.
Co-rotating frameEdit
For a particle in planar motion, one approach to attaching physical significance to these terms is based on the concept of an instantaneous co-rotating frame of reference.<ref name="Taylor">For the following discussion, see Template:Cite book</ref> To define a co-rotating frame, first an origin is selected from which the distance r(t) to the particle is defined. An axis of rotation is set up that is perpendicular to the plane of motion of the particle, and passing through this origin. Then, at the selected moment t, the rate of rotation of the co-rotating frame Ω is made to match the rate of rotation of the particle about this axis, dφ/dt. Next, the terms in the acceleration in the inertial frame are related to those in the co-rotating frame. Let the location of the particle in the inertial frame be (r(t), φ(t)), and in the co-rotating frame be (r′(t), φ′(t)). Because the co-rotating frame rotates at the same rate as the particle, dφ′/dt = 0. The fictitious centrifugal force in the co-rotating frame is mrΩ2, radially outward. The velocity of the particle in the co-rotating frame also is radially outward, because dφ′/dt = 0. The fictitious Coriolis force therefore has a value −2m(dr/dt)Ω, pointed in the direction of increasing φ only. Thus, using these forces in Newton's second law we find: <math display="block">\mathbf{F} + \mathbf{F}_\text{cf} + \mathbf{F}_\text{Cor} = m\ddot{\mathbf{r}} \, , </math> where over dots represent derivatives with respect to time, and F is the net real force (as opposed to the fictitious forces). In terms of components, this vector equation becomes: <math display="block">\begin{align}
F_r + mr\Omega^2 &= m\ddot{r} \\
F_\varphi - 2m\dot{r}\Omega &= mr\ddot{\varphi} \ ,
\end{align}</math> which can be compared to the equations for the inertial frame: <math display="block">\begin{align}
F_r &= m\ddot{r} - mr\dot{\varphi}^2 \\
F_\varphi &= mr\ddot{\varphi} + 2m\dot{r}\dot{\varphi} \ .
\end{align}</math>
This comparison, plus the recognition that by the definition of the co-rotating frame at time t it has a rate of rotation Ω = dφ/dt, shows that we can interpret the terms in the acceleration (multiplied by the mass of the particle) as found in the inertial frame as the negative of the centrifugal and Coriolis forces that would be seen in the instantaneous, non-inertial co-rotating frame.
For general motion of a particle (as opposed to simple circular motion), the centrifugal and Coriolis forces in a particle's frame of reference commonly are referred to the instantaneous osculating circle of its motion, not to a fixed center of polar coordinates. For more detail, see centripetal force.
Differential geometryEdit
In the modern terminology of differential geometry, polar coordinates provide coordinate charts for the differentiable manifold Template:Math, the plane minus the origin. In these coordinates, the Euclidean metric tensor is given by<math display="block">ds^2 = dr^2 + r^2 d\theta^2.</math>This can be seen via the change of variables formula for the metric tensor, or by computing the differential forms dx, dy via the exterior derivative of the 0-forms Template:Math, Template:Math and substituting them in the Euclidean metric tensor Template:Math.
Template:Collapse top Let <math>p_1=(x_1,y_1)=(r_1,\theta_1)</math>, and <math>p_2=(x_2,y_2)=(r_2,\theta_2)</math> be two points in the plane given by their cartesian and polar coordinates. Then
- <math>ds^2=dx^2+dy^2=(x_2-x_1)^2+(y_2-y_1)^2.</math>
Since <math>dx^2=(r_2\cos\theta_2-r_1\cos\theta_1)^2</math>, and <math>dy^2=(r_2\sin\theta_2-r_1\sin\theta_1)^2</math>, we get that
- <math>ds^2=r_2^2\cos^2\theta_2-2r_1r_2\cos\theta_1\cos\theta_2+r_1^2\cos^2\theta_1+r_2^2\sin^2\theta_2-2r_1r_2\sin\theta_1\sin\theta_2+r_1^2\sin^2\theta_1=</math>
- <math>r_2^2(\cos^2\theta_2+\sin^2\theta_2)+r_1^2(\cos^2\theta_1+\sin^2\theta_1)-2r_1r_2(\cos\theta_1\cos\theta_2+\sin\theta_1\sin\theta_2)=</math>
- <math>r_1^2+r_2^2-2r_1r_2(1-1+\cos\theta_1\cos\theta_2+\sin\theta_1\sin\theta_2)=</math>
- <math>(r_2-r_1)^2+2r_1r_2(1-\cos\theta_1\cos\theta_2-\sin\theta_1\sin\theta_2).</math>
Now we use the trigonometric identity <math>\cos(\theta_2-\theta_1)=\cos\theta_1\cos\theta_2+\sin\theta_1\sin\theta_2</math> to proceed:
- <math>ds^2=dr^2+2r_1r_2(1-\cos d\theta).</math>
If the radial and angular quantities are near to each other, and therefore near to a common quantity <math>r</math> and <math>\theta</math>, we have that <math>r_1r_2\approx r^2</math>. Moreover, the cosine of <math>d\theta</math> can be approximated with the Taylor series of the cosine up to linear terms:
- <math>\cos d\theta\approx1-\frac{d\theta^2}{2},</math>
so that <math>1-\cos d\theta\approx\frac{d\theta^2}{2}</math>, and <math>2r_1r_2(1-\cos d\theta)\approx2r^2\frac{d\theta^2}{2}=r^2d\theta^2</math>. Therefore, around an infinitesimally small domain of any point,
- <math>ds^2=dr^2+r^2d\theta^2,</math>
as stated.
An orthonormal frame with respect to this metric is given by<math display="block">e_r = \frac{\partial}{\partial r}, \quad e_\theta = \frac1r \frac{\partial}{\partial \theta},</math>with dual coframe<math display="block">e^r = dr, \quad e^\theta = r d\theta.</math>The connection form relative to this frame and the Levi-Civita connection is given by the skew-symmetric matrix of 1-forms<math display="block">{\omega^i}_j = \begin{pmatrix} 0 & -d\theta \\ d\theta & 0\end{pmatrix}</math>and hence the curvature form Template:Math vanishes. Therefore, as expected, the punctured plane is a flat manifold.
Extensions in 3DEdit
The polar coordinate system is extended into three dimensions with two different coordinate systems, the cylindrical and spherical coordinate system.
ApplicationsEdit
Polar coordinates are two-dimensional and thus they can be used only where point positions lie on a single two-dimensional plane. They are most appropriate in any context where the phenomenon being considered is inherently tied to direction and length from a center point. For instance, the examples above show how elementary polar equations suffice to define curves—such as the Archimedean spiral—whose equation in the Cartesian coordinate system would be much more intricate. Moreover, many physical systems—such as those concerned with bodies moving around a central point or with phenomena originating from a central point—are simpler and more intuitive to model using polar coordinates. The initial motivation for the introduction of the polar system was the study of circular and orbital motion.
Edit
Polar coordinates are used often in navigation as the destination or direction of travel can be given as an angle and distance from the object being considered. For instance, aircraft use a slightly modified version of the polar coordinates for navigation. In this system, the one generally used for any sort of navigation, the 0° ray is generally called heading 360, and the angles continue in a clockwise direction, rather than counterclockwise, as in the mathematical system. Heading 360 corresponds to magnetic north, while headings 90, 180, and 270 correspond to magnetic east, south, and west, respectively.<ref>{{#invoke:citation/CS1|citation |CitationClass=web }}</ref> Thus, an aircraft traveling 5 nautical miles due east will be traveling 5 units at heading 90 (read zero-niner-zero by air traffic control).<ref>{{#invoke:citation/CS1|citation |CitationClass=web }}</ref>
ModelingEdit
Systems displaying radial symmetry provide natural settings for the polar coordinate system, with the central point acting as the pole. A prime example of this usage is the groundwater flow equation when applied to radially symmetric wells. Systems with a radial force are also good candidates for the use of the polar coordinate system. These systems include gravitational fields, which obey the inverse-square law, as well as systems with point sources, such as radio antennas.
Radially asymmetric systems may also be modeled with polar coordinates. For example, a microphone's pickup pattern illustrates its proportional response to an incoming sound from a given direction, and these patterns can be represented as polar curves. The curve for a standard cardioid microphone, the most common unidirectional microphone, can be represented as Template:Nowrap at its target design frequency.<ref>Template:Cite book</ref> The pattern shifts toward omnidirectionality at lower frequencies.
See alsoEdit
- Curvilinear coordinates
- List of common coordinate transformations
- Log-polar coordinates
- Polar decomposition
- Unit circle
ReferencesEdit
General referencesEdit
External linksEdit
- Template:Springer
- Coordinate Converter — converts between polar, Cartesian and spherical coordinates
- Polar Coordinate System Dynamic Demo
Template:Orthogonal coordinate systems Template:Authority control