This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

The cocked hat

Imre Bárány, William Steiger, Sivan Toledo
Abstract.

We revisit the cocked hat – an old problem from navigation – and examine under what conditions its old solution is valid.

Key words and phrases:
Random rays, geometric probability, navigation
2000 Mathematics Subject Classification:
Primary 52A22, secondary 60D05

1. Introduction

Navigators used to plot on a map lines of position or lines of bearing, which are rays emanating from a landmark (e.g., a lighthouse or radio beacon) at a particular bearing (angle relative to north) that was estimated to be the direction from the landmark (which we also refer to as observation point) to the ship or plane. Two such rays usually intersect at a point, which the navigator would take as an estimate of the true position of the craft. Navigators were encouraged to plot three rays, to make position estimation more robust. The three rays normally created a triangle, called a cocked hat [6], as shown in Figure 1. The properties of the cocked hat were investigated thoroughly [1, 3, 4, 5, 8, 9, 10, 12], to help navigators interpret it and make good navigation decisions. The aim of this paper is to analyze the conditions under which an elegant property of the cocked hat holds. That property had been stated without a proof more than 80 years ago [8], proved informally (and essentially incorrectly) 70 years ago [9], and has been widely disseminated ever since [3, 5, 7, 11, 12], including in course material [11] and in a popular science book [7].

The property that we are interested in is the probability of the cocked hat containing the true position being 1/41/4. Under what conditions is this statement true?

Refer to caption
Figure 1. Three observation points P1,P2,P3P_{1},P_{2},P_{3}, the target FF, the three rays, and the cocked hat (shaded).

This claim first appeared in a 1938 navigation manual [8, page 166], without a proof and with only informal conditions on the error angles at the three landmarks, which we denote P1P_{1}, P2P_{2}, and P3P_{3} (see Figure 1). The error angles ϵ1\epsilon_{1}, ϵ2\epsilon_{2}, and ϵ3\epsilon_{3} are between the plotted rays, which we denote R1R_{1}, R2R_{2}, and R3R_{3} and the rays rir_{i} from PiP_{i} to the true position of the craft, which we denote by FF (see Figure 1). The informal conditions are that the errors are independent (the manual does not use this term, but this is what it means) and fairly small, around 1 degree. A 1947 article by Stansfield [9]111Stansfield developed the results published in the paper while serving in Operational Research Sections attached to the Royal Air Force Fighter Command and Coastal Command during World War II. cites the claim, gives more formal conditions for it, and sketches a proof. The conditions that Stansfield specified are remarkably weak: he claims that the result would hold if only two of the three errors have zero median. Stansfield writes that this assumption is equivalent to the following: “for two of the stations the observed bearings are equally likely to pass to the right or the left of the true position”. A 1951 article by Daniels [5] states Stansfield’s result in a more modern statistical language, saying that the cocked hat is a 25%25\% distribution-free confidence region; the term distribution free means that the result is not dependent on a particular error distribution, say Gaussian, but only on a parameter of the distribution, here the zero median222Daniels was a statistician and served as the president of the Royal Statistical Society from 1974 to 1975. His paper incorrectly states that the Admiralty Navigation Manual proves the 25%25\%-probability result; it does not; the first proof sketch appears in Stansfield’s paper. . Daniels then considers the case of nn landmarks and nn rays starting from there. The lines of these rays split the plane into finitely connected components, some of them bounded, some of them not. Daniels claims without proof a particular formula, 2n2n\frac{2n}{2^{n}}, for the probability that FF belongs to the union of the unbounded components. The 25%25\%-probability result was incorrectly extended again by Williams333Williams was a professional air navigator and served as president of the Royal Institute of Navigation from 1984-1987 [2]. in 1991. He claimed specific probabilities that the open regions around the cocked hat contain FF, again with only an informal specification of the assumptions and with only a sketch of the proof. Williams’s claims were shown to be false by Cook [3], using specific error distributions to which Williams answered with a witty (but scientifically wrong) rebuttal. Cook also repeated the claim that the probability of the cocked hat contains FF is 1/41/4.

Our aim in this paper is to show that the 25%25\%-probability result is valid only for error distributions that guarantee that the three rays intersect at three distinct points and form a triangle.

We note that the use of the cocked hat in navigation is today obsolete, having been replaced by estimation of confidence regions, usually circles or ellipses, by computer algorithms.

2. Generalizations to rays that do not intersect

Two rays in the plane can intersect, but they can also fail to intersect. Lines of position plotted by navigators almost always intersected, because the error angles were small. Also, navigators were taught to choose landmarks so that no angle at the intersection is smaller than about 50 degrees – a small angle at the intersection implies ill conditioning (high sensitivity of the intersection point to bearing errors).

Stansfield’s formulation of the problem uses much more general assumptions on the errors, and no assumption about angles at the intersections. Stansfield, Daniels, and the authors that followed only require that the three errors ε1,ε2,ε3(π,π]\varepsilon_{1},\varepsilon_{2},\varepsilon_{3}\in(-\pi,\pi] are random, independent, and that the median of their distributions is zero. We replace the zero-median assumption by a consistent but slightly more general condition, namely

(2.1) Prob(ϵi<0)=Prob(ϵi>0)=12\mathrm{Prob}(\epsilon_{i}<0)=\mathrm{Prob}(\epsilon_{i}>0)=\frac{1}{2}

for every i[n]i\in[n], where [n][n] is a shorthand for the set {1,2,,n}\{1,2,\ldots,n\}. This means in particular that the target is never on RiR_{i} which is a necessity because if Prob(ϵi=0)>0\mathrm{Prob}(\epsilon_{i}=0)>0 were allowed, then Prob(FΔ)\mathrm{Prob}(F\in\Delta) could be close to one (e.g., if Prob(ϵi=0)\mathrm{Prob}(\epsilon_{i}=0) is close to one), implying that the 1/41/4 result does not hold in this case. We also consider the restriction of the errors to [π/2,π/2][-\pi/2,\pi/2].

Under these weak assumptions on the error distribution, the three rays might fail to form a triangle (the cocked hat). How can we formally express the 25%25\%-probability result when rays may fail to intersect? We propose four ways to express the result; the first three are fairly natural but are not sufficient for the 1/41/4 result, even under the restriction ϵi[π/2,π/2]\epsilon_{i}\in[-\pi/2,\pi/2]; the fourth is not particularly natural but is the only correct statement of the result.

Conjunction formulation. The probability that the three rays intersect at three points and that the triangle that they form contains FF is 1/41/4. In this formulation, we allow error distributions that could generate non-intersecting rays and we hope to prove that the probability that the rays intersect at fewer than three points or that the triangle does not contain FF is exactly 3/43/4. This is false.

Conditional probability formulation. The conditional probability that the triangle that the rays form contains FF, conditioned on the rays forming a triangle, is 1/41/4. In this formulation we again allow error distributions that generate non-intersecting rays, and we hope to prove that if the rays intersect at three points, then the probability that the triangle contains FF is 1/41/4. We do not care with what probability the rays fail to form a triangle. This again is false.

Lines formulation. We extend the rays rir_{i} to infinite lines i\ell_{i}, which always form a triangle, and we hope to show that the triangle that they form contains FF with probability 1/41/4. Here we must restrict ϵi[π/2,π/2]\epsilon_{i}\in[-\pi/2,\pi/2], otherwise the same line could appear both on the left and on the right of FF. Again, this claim is false.

Constrained distribution formulation. We assume that the distribution of errors is such that every pair of rays always intersects and we hope to show that the probability that the triangle contains FF is 1/41/4. We do not permit distributions under which two of the rays might fail to intersect. We show below that in this case Prob(FΔ)=14\mathrm{Prob}(F\in\Delta)=\frac{1}{4}.

We note that from the navigator’s perspective, the conditional probability is the most natural. You plot three rays. If they do not intersect at three points, you discard the measurements and try again, because you either picked bad observation points (e.g., two of them and your ship are almost collinear) or at least one of the bearings is way off. If they do intersect at three points, you want to know the (conditional) probability that the cocked hat contains FF. From the statistician’s perspective, any of the first three formulations makes sense. The fourth makes less statistical sense, because it is unusual to assume that independent error distributions satisfy some global structural constraint. In particular, it appears that Daniels may have believed that the lines formulation is correct, because he writes about geometrical lines in the plane, not about rays. He writes “a particular set of nn lines, no two of which are parallel, divides the plane in to 12(n2+n+2)\frac{1}{2}(n^{2}+n+2) polygons”.

3. Counterexamples

We now show that the Conjunction, Conditional probability, and Lines formulation are all false by giving counterexamples. Every example is a two-ray distribution that is concentrated on two rays Ri+R_{i}^{+} and RiR_{i}^{-}: Prob(Ri=Ri+)=Prob(Ri=Ri)=12\mathrm{Prob}(R_{i}=R_{i}^{+})=\mathrm{Prob}(R_{i}=R_{i}^{-})=\frac{1}{2}. This is no coincidence as we will see at the end of this section.

Refer to caption
Figure 2. Two counterexamples.

In the first example FF is in the centroid of an equilateral triangle whose vertices are P1P_{1}, P2P_{2}, and P3P_{3}. Figure 2 (left) shows the two-ray error distributions. It is easy to see that Ri+R_{i}^{+} intersects neither Ri+1+R_{i+1}^{+} nor Ri+1R_{i+1}^{-} (subscripts are meant modulo 33). Therefore, if Ri+R_{i}^{+} is selected, then a cocked hat does not form. On the other hand, if R1R_{1}^{-}, R2R_{2}^{-}, and R3R_{3}^{-} are selected, then they form a cocked hat that contains FF. Therefore,

Prob(F| the rays form a cocked hat )\displaystyle\mathrm{Prob}\left(F\in\triangle|\text{ the rays form a cocked hat }\triangle\right) =\displaystyle= 1\displaystyle 1
Prob(the rays form a cocked hat  and F)\displaystyle\mathrm{Prob}\left(\text{the rays form a cocked hat }\triangle\text{ and }F\in\triangle\right) =\displaystyle= 18.\displaystyle\frac{1}{8}\;.

This shows that both the Conjunction formulation is false and that the Conditional probability formulation is false. Note that all the error angles have magnitude less than π/2\pi/2, so these formulations are false even with this restriction.

Figure 2 (right) shows another two-ray distribution. The error magnitudes are less than π/2\pi/2, actually as small as you wish. The true position FF lies outside all the triangles that the lines form, so the probability that the cocked hat (in the Lines formulation) contains FF is zero. We can move FF to the right by any amount and FΔF\notin\Delta will still hold. This example also shows that the conditional probability that a cocked hat formed by 33 rays contains FF can also be zero.

Refer to caption
Figure 3. The third counterexample.

The last example, given in Figure 3, shows that the probability that the triangle formed by the extension of the rays to lines contains FF can be 11. We again note that the error angles are bounded in magnitude by π/2\pi/2. In this example the three rays do not have three intersection points, so the cocked hat appears with probability zero. So this is another counterexample to the Conjunction formulation.

We close this section with a remark on two-ray distributions. The set of (Borel) probability distributions satisfying condition (2.1) is convex, and its extreme points are exactly the two-ray distributions, as one can easily check. Moreover Prob(FΔ)\mathrm{Prob}(F\in\Delta) is a linear function on the product of the distributions μ1,μ2,μ3\mu_{1},\mu_{2},\mu_{3} where μi\mu_{i} is the probability distribution of the ray RiR_{i}. Indeed, denoting by I(E)I(E) the indicator function of an event EE, we have

(3.1) Prob(FΔ)=I(FΔ)𝑑μ1𝑑μ2𝑑μ3,\mathrm{Prob}(F\in\Delta)=\int I(F\in\Delta)d\mu_{1}d\mu_{2}d\mu_{3},

a linear function of each μi\mu_{i}, so if it takes the value 14\frac{1}{4} on the two-ray distributions, then it takes the same value on all distributions satisfying (2.1). We will come back to such distributions in Section 5 again.

4. Intersecting rays

We now start the analysis when rays must intersect in pairs. We assume throughout that the n+1n+1 points P1,,Pn,FP_{1},\ldots,P_{n},F are in general position, so that no three are collinear and so that no other degeneracies arise.

We introduce some notation. We let XY\overrightarrow{XY} denote the ray emanating from XX in the direction of YY when X,YX,Y are distinct points in the plane; here we assume that XXYX\notin\overrightarrow{XY}. Thus ri=PiFr_{i}=\overrightarrow{P_{i}F} is the ray starting at PiP_{i} in the direction of the target FF, and i\ell_{i} is the line containing rir_{i}. From each PiP_{i} out goes a random ray RiR_{i} making a (signed) angle εi(π,π)\varepsilon_{i}\in(-\pi,\pi) with rir_{i}. Our basic assumption, besides (2.1), is that two random rays always intersect that is for distinct i,j[n]i,j\in[n]

(4.1) Prob(RiRj=)=0.\mathrm{Prob}(R_{i}\cap R_{j}=\emptyset)=0.

So ray RiR_{i} and RjR_{j} intersect almost surely but their intersection point is not PiP_{i} or PjP_{j} because of our convention that XXYX\notin\overrightarrow{XY} .

Further notations: hih_{i}^{-} resp. hi+h_{i}^{+} are the halfplanes bounded by i\ell_{i} with hih_{i}^{-} consisting of points XX such that the ray PiX\overrightarrow{P_{i}X} comes from a clockwise rotation from rir_{i} with angle less than π\pi, and hi+h_{i}^{+} is its complementary halfplane. When r,rr,r^{\prime} are two rays we denote by cone(Pi,r,r)\mathrm{cone}(P_{i},r,r^{\prime}) the cone whose apex is PiP_{i} and whose bounding rays are translated copies of rr and rr^{\prime}. Such a cone always has angle less than π\pi, because rr and rr^{\prime} will never have opposite directions.

Define Cij=cone(Pi,rj,PiPj)C_{ij}=\mathrm{cone}(P_{i},r_{j},\overrightarrow{P_{i}P_{j}}) for distinct i,j[n]i,j\in[n].

Lemma 4.1.

The cone CijC_{ij} contains rir_{i} and Prob(RiCij)=1\mathrm{Prob}(R_{i}\subset C_{ij})=1.

Proof. Assume first that PjhiP_{j}\in h_{i}^{-}. We define first the cones Cij=cone(Pi,ri,PiPj)C_{ij}^{-}=\mathrm{cone}(P_{i},r_{i},\overrightarrow{P_{i}P_{j}}) and Cij+=cone(Pi,ri,rj)C_{ij}^{+}=\mathrm{cone}(P_{i},r_{i},r_{j}), see Figure 4. Note that the angle of CijC_{ij}^{-} (resp. Cij+C_{ij}^{+}) is equal to the angle at PiP_{i} (and at FF) of the triangle with vertices Pi,Pj,FP_{i},P_{j},F. Then Cij=CijCij+C_{ij}=C_{ij}^{-}\cup C_{ij}^{+} because the angle of this cone is the sum of the angles of CijC_{ij}^{-} and Cij+C_{ij}^{+} so smaller than π\pi. Then riCijr_{i}\subset C_{ij} indeed as shown in Figure 4, left.

Refer to caption
Figure 4. Illustration for Lemma 4.1: the case PjhiP_{j}\in h_{i}^{-} on the left, the case Pjhi+P_{j}\in h_{i}^{+} on the right.

Suppose now that εi>0\varepsilon_{i}>0 which is the same as Rihi+R_{i}\subset h_{i}^{+}. If RiR_{i} does not lie in Cij+C_{ij}^{+}, then Rihi+Cij+R_{i}\subset h_{i}^{+}\setminus C_{ij}^{+}. The last set is a convex cone, disjoint from hjh_{j}^{-}, as they are separated by the line j\ell_{j}. So no RjR_{j} with εj<0\varepsilon_{j}<0 can intersect RiR_{i} contradicting (4.1). So RiCij+R_{i}\subset C_{ij}^{+}.

Let hh denote the halfplane containing FF and bounded by the line through PiP_{i} and PjP_{j}. Observe that by the previous argument RjhR_{j}\subset h because the complementary halfplane to hh is disjoint from Cij+C_{ij}^{+}, so RjR_{j} can intersect RiCij+R_{i}\subset C_{ij}^{+} only if it lies in hh.

Suppose next that εi<0\varepsilon_{i}<0. We show that RiCijR_{i}\subset C_{ij}^{-}. If not, then RihiCijR_{i}\subset h_{i}^{-}\setminus C_{ij}^{-}. The last set is a convex cone again, disjoint from hh, so RiRj=R_{i}\cap R_{j}=\emptyset for all RjR_{j} with εj<0\varepsilon_{j}<0 contradicting (4.1).

The argument for the case Pjhi+P_{j}\in h_{i}^{+} is symmetric (see Figure 4 right) but otherwise identical and is therefore omitted.∎

We remark here that Lemma 4.1 implies that the cone jiCij\bigcap_{j\neq i}C_{ij} is convex (that is, its angle is smaller than π\pi), it contains rir_{i}, and Prob(RijiCij)=1\mathrm{Prob}(R_{i}\subset\bigcap_{j\neq i}C_{ij})=1, of course only if n2n\geq 2. (For n=1n=1 condition (4.1) is void.) Define KiK_{i} as the smallest (with respect to inclusion) convex cone satisfying Prob(RiKi)=1\mathrm{Prob}(R_{i}\subset K_{i})=1. Note that KijiCijK_{i}\subset\bigcap_{j\neq i}C_{ij}. For later reference we state the following corollary.

Corollary 4.1.

Under conditions (2.1) and (4.1) KiK_{i} is a convex cone, riKir_{i}\subset K_{i} and Prob(RiKi)=1\mathrm{Prob}(R_{i}\subset K_{i})=1 for every i[n]i\in[n].

Theorem 4.1.

Under conditions (2.1) and (4.1)

Prob(FΔ)=14.\mathrm{Prob}(F\in\Delta)=\frac{1}{4}.

Proof. Set T=conv{P1,P2,P3,F}T=\mathrm{conv}\{P_{1},P_{2},P_{3},F\}, the convex hull of P1,P2,P3P_{1},P_{2},P_{3}, and FF. We will have to consider three cases separately: when TT is a triangle with FF inside TT (Case 1), when TT is a triangle with FF a vertex of TT (Case 2), and when TT is a quadrilateral (Case 3).

Case 1. Define Ci=cone(Pi,PiPi1,PiPi+1)C_{i}=\mathrm{cone}(P_{i},\overrightarrow{P_{i}P_{i-1}},\overrightarrow{P_{i}P_{i+1}}) for i=1,2,3i=1,2,3 where the subscripts are taken mod 33, see Figure 5 left.

We claim that RiCiR_{i}\subset C_{i} for all ii. By symmetry it suffices to show this for i=2i=2. By Lemma 4.1 R2C21C23R_{2}\subset C_{21}\cap C_{23}. So it is enough to check that C2=C21C23C_{2}=C_{21}\cap C_{23}, and this is evident: the rays bounding C2C_{2} are P2P1\overrightarrow{P_{2}P_{1}} (which bounds C21C_{21}) and P2P3\overrightarrow{P_{2}P_{3}} (which bounds C23C_{23}).

Refer to caption
Figure 5. The cone C2C_{2} in Case 1 (left) and 2 (right).
Refer to caption
Figure 6. Illustration for the proof of Theorem 4.1.

We can now finish the proof of the theorem in Case 1. There are 88 sub-cases with equal probabilities that correspond to the signs of ϵ1\epsilon_{1}, ϵ2\epsilon_{2}, and ϵ3\epsilon_{3}, as shown in Figure 6. Only in two of them, namely when all ϵi\epsilon_{i}s have the same sign, we have FF\in\triangle, so the probability of this event is 1/41/4.

Case 2. We assume (by symmetry) that P2P_{2} is inside the triangle TT. We define the cones C1=cone(P1,r2,P1P2)C_{1}=\mathrm{cone}(P_{1},r_{2},\overrightarrow{P_{1}P_{2}}), C2=cone(P2,r1,r3)C_{2}=\mathrm{cone}(P_{2},r_{1},r_{3}), and C3=cone(P3,r2,P3P2)C_{3}=\mathrm{cone}(P_{3},r_{2},\overrightarrow{P_{3}P_{2}}) and we claim that RiCiR_{i}\subset C_{i} for all ii. From Lemma 4.1 we have that R2C21C23R_{2}\subset C_{21}\cap C_{23}. The bounding rays of C2C_{2} are a translate of r2r_{2} (bounding C21C_{21}) and a translate of r3r_{3} (bounding C23C_{23}), so C2=C21C23C_{2}=C_{21}\cap C_{23} (see Figure 5 right).

The cases i=1i=1 and 33 are symmetric and very simple. We only consider i=1i=1. Again, by Lemma 4.1 R1C12R_{1}\subset C_{12} and then C1=C12C_{1}=C_{12} implying R1C1R_{1}\subset C_{1}.

Again there are 8 subcases, corresponding to the 8 possible sign patterns of ε1,ε2,ε3\varepsilon_{1},\varepsilon_{2},\varepsilon_{3}. It is easy to see that FΔF\in\Delta in exactly two of them.

Case 3. We assume again by symmetry that the segment P2FP_{2}F is a diagonal of the quadrilateral TT. Define cones C1=cone(P1,r2,P1P3)C_{1}=\mathrm{cone}(P_{1},r_{2},\overrightarrow{P_{1}P_{3}}), C2=cone(P2,r1,r3)C_{2}=\mathrm{cone}(P_{2},r_{1},r_{3}), and C3=cone(P3,r2,P3P2)C_{3}=\mathrm{cone}(P_{3},r_{2},\overrightarrow{P_{3}P_{2}}). We claim again that RiCiR_{i}\subset C_{i} for all ii. The proof is similar to the previous ones using Lemma 4.1 and is omitted here. Again, FΔF\in\Delta in exactly two out of the 8 cases.∎

5. Daniels’ statement

We assume now that there are n3n\geq 3 observation points P1,,PnP_{1},\ldots,P_{n} plus the target point FF and that these n+1n+1 points are in general position. A random ray RiR_{i} starts at each PiP_{i} satisfying conditions (2.1) and (4.1). The lines of the rays RiR_{i} split the plane into connected components, let UU denote the union of the 2n2n unbounded components. Here comes Daniels’ statement.

Theorem 5.1.

Under conditions (2.1) and (4.1)

Prob(FU)=2n2n.\mathrm{Prob}(F\in U)=\frac{2n}{2^{n}}.

The case n=2n=2 is trivial and not interesting. The case n=3n=3 is just Theorem 4.1. We note that condition (4.1) is a necessity, even for n=3n=3 as the counterexamples in Section 3 show.

We are going to prove this theorem under the assumption that each RiR_{i} is a two-ray distribution, that is, Prob(Ri=Ri+)=Prob(Ri=Ri)=12\mathrm{Prob}(R_{i}=R_{i}^{+})=\mathrm{Prob}(R_{i}=R_{i}^{-})=\frac{1}{2} and explain, after the proof, how this special case implies the theorem. We also assume that the 2n2n rays Ri+,RiR_{i}^{+},R_{i}^{-}, together with the points P1,,Pn,FP_{1},\ldots,P_{n},F are in general position. This is not a serious restriction because the general case of two-ray distributions follows from this by a routine limiting argument.

Proof. To simplify the writing, we set Di=cone(Pi,Ri+,Ri)D_{i}=\mathrm{cone}(P_{i},R_{i}^{+},R_{i}^{-}), which is equivalent to Di=conv(Ri+Ri)D_{i}=\mathrm{conv}(R_{i}^{+}\cup R_{i}^{-}). Lemma 4.1 implies that riDir_{i}\subset D_{i} for every i[n]i\in[n]. Let SS be a circle centered at FF such that SDiS\subset D_{i} for every i[n]i\in[n]. Observe that for distinct i,j[n]i,j\in[n], the intersection DiDjD_{i}\cap D_{j} is a convex quadrilateral containing SS and of course FF, see Figure 7 left. This follows from condition (4.1): both Ri+R_{i}^{+} and RiR_{i}^{-} intersect both Rj+R_{j}^{+} and RjR_{j}^{-} and the four intersection points are the vertices of DiDjD_{i}\cap D_{j} which is then a convex quadrilateral.

Refer to caption
Figure 7. The intersection DiDjD_{i}\cap D_{j} and the translated cone DiD_{i}^{*}.

Let Li+L_{i}^{+} (resp. LiL_{i}^{-}) denote the line of the ray Ri+R_{i}^{+} (and RiR_{i}^{-}). For a selection δ1,,δn{1,1}\delta_{1},\ldots,\delta_{n}\in\{1,-1\} of signs the lines L1δ1,,LnδnL_{1}^{\delta_{1}},\ldots,L_{n}^{\delta_{n}} split the plane into finitely many connected components. We are going to show that out of the 2n2^{n} possible selections there are exactly 2n2n for which FF lies in an unbounded component.

We reduce this statement to another one about arcs on the unit circle. First comes a simpler reduction. Translate each cone DiD_{i} into a new (and actually unique) position DiD_{i}^{*} so that its rays touch the circle SS (see Figure 7 right). Let Qi+,Mi+Q_{i}^{+},M_{i}^{+} (resp. Qi,MiQ_{i}^{-},M_{i}^{-}) be the translated copies of Ri+,Li+R_{i}^{+},L_{i}^{+} (and Ri,LiR_{i}^{-},L_{i}^{-}). Note that DiDjD_{i}^{*}\cap D_{j}^{*} is again a convex quadrilateral.

We claim next that for a fixed selection δ1,,δn\delta_{1},\ldots,\delta_{n} of signs, FF lies in an unbounded component for the lines L1δ1,,LnδnL_{1}^{\delta_{1}},\ldots,L_{n}^{\delta_{n}} if and only if it lies in the corresponding unbounded component for the lines M1δ1,,MnδnM_{1}^{\delta_{1}},\ldots,M_{n}^{\delta_{n}}. This is simple. The point FF lies in an unbounded component for the lines LiδiL_{i}^{\delta_{i}} if and only if there is a halfline RR starting at FF and disjoint from each LiδiL_{i}^{\delta_{i}} which happens if and only if RR is disjoint from the lines MiδiM_{i}^{\delta_{i}} as well.

Assume now that SS is the unit circle. Let ai+a_{i}^{+} (resp. aia_{i}^{-}) be the points where Mi+M_{i}^{+} (and MiM_{i}^{-}) touch SS, and let AiA_{i} be the shorter arc on SS between ai+a_{i}^{+} and aia_{i}^{-}, see Figure  7 right. It is clear that ai+a_{i}^{+} and aia_{i}^{-} are not opposite points on SS so AiA_{i} is welldefined. These arcs completely determine DiD_{i}^{*}. They satisfy the conditions

  • (i)

    each AiA_{i} is shorter than π\pi, and

  • (ii)

    AiAjA_{i}\cup A_{j} is an arc in SS longer than π\pi for all i,j[n],iji,j\in[n],i\neq j.

The latter condition follows from the fact that DiDjD_{i}^{*}\cap D_{j}^{*} is a convex quadrilateral.

We call a selection δ1,,δn\delta_{1},\ldots,\delta_{n} special if it gives an unbounded component containing FF. We claim that a selection is special if and only if the points a1δ1,,anδna_{1}^{\delta_{1}},\ldots,a_{n}^{\delta_{n}} lie on an arc of SS shorter than π\pi. This is also simple. If there is such an arc, call it II and let QQ be the centre point of the complementary arc SIS\setminus I. The ray FQ\overrightarrow{FQ} avoids every line MiδiM_{i}^{\delta_{i}}. If there is no such arc, then the connected component containing FF (and SS) is bounded as one can check easily. Therefore it suffices to prove the following lemma.

Lemma 5.1.

Under the above conditions there are exactly 2n2n special selections.

Proof. For a special selection δ=(δ1,,δn)\delta=(\delta_{1},\ldots,\delta_{n}) let I(δ)I(\delta) denote the shortest arc on SS containing every aiδia_{i}^{\delta_{i}}, i[n]i\in[n]. Thus I(δ)I(\delta) is the shorter arc between points aiδia_{i}^{\delta_{i}} and ajδja_{j}^{\delta_{j}} for some distinct i,j[n]i,j\in[n], and they are the endpoints of I(δ)I(\delta).

Claim 5.1.

Each ai+a_{i}^{+} (and aia_{i}^{-}) is the endpoint of I(δ)I(\delta) for exactly two special selections δ\delta.

It suffices to prove this claim since it implies Lemma 5.1 and then Theorem 5.1 as well.

Proof  of the claim. It is enough to work with a1+a_{1}^{+}. Using the notation on Figure 8 we assume that a1a_{1}^{-} is from the halfcircle S+S^{+} so A1S+A_{1}\subset S^{+}.

Refer to caption
Figure 8. The definition of S+S^{+} and SS^{-}.

Define X={a1+,a1,,an+,an}X=\{a_{1}^{+},a_{1}^{-},\ldots,a_{n}^{+},a_{n}^{-}\} and Y=X{a1+,a1}Y=X\setminus\{a_{1}^{+},a_{1}^{-}\}. Observe first that S+S^{+} can’t contain any AiA_{i}, i>1i>1 as otherwise A1,AiS+A_{1},A_{i}\subset S^{+} contradicting (ii). Moreover, SS^{-} can’t contain two arcs Ai,AjA_{i},A_{j} with distinct i,j>1i,j>1 because of (ii) again. It follows then that |S+Y|=n1|S^{+}\cap Y|=n-1 or n2n-2.

Case 1 when |S+Y|=n1|S^{+}\cap Y|=n-1. Then |SY|=n1|S^{-}\cap Y|=n-1 as well and S+S^{+} contains exactly one element from each pair {ai+,ai}\{a_{i}^{+},a_{i}^{-}\}, i>1i>1, and then so does SS^{-}. This gives exactly two special selection δ\delta and ε\varepsilon with I(δ)S+I(\delta)\subset S^{+} and I(ε)SI(\varepsilon)\subset S^{-}, with a1+a_{1}^{+} an endpoint of both.

Case 2 when |S+Y|=n2|S^{+}\cap Y|=n-2. Then S+S^{+} contains no I(δ)I(\delta) with δ\delta special, |SY|=n|S^{-}\cap Y|=n and so AiSA_{i}\subset S^{-} for a unique i>1i>1. This gives again two special selections δ\delta and ε\varepsilon where a1+a_{1}^{+} is the endpoint of I(δ)I(\delta) and I(ε)I(\varepsilon). In fact δ\delta and ε\varepsilon coincide except at position ii: δj=εj\delta_{j}=\varepsilon_{j} for all j[n]j\in[n] but j=ij=i and δ1=ε1=1\delta_{1}=\varepsilon_{1}=1.∎

We explain now how the case of two-ray distributions implies Theorem 5.1, or rather give a sketch of this and leave the technical details to the interested reader. Assume each ray RiR_{i} follows a generic distribution μi\mu_{i} for all i[n]i\in[n] still satisfying conditions (2.1) and (4.1). Note that by Corollary 4.1, Prob(RiKi)=1\mathrm{Prob}(R_{i}\subset K_{i})=1. Using this one can check that every μi\mu_{i} can be approximated with high precision by a convex combination of two-ray distributions, each having Ri+,RiKiR_{i}^{+},R_{i}^{-}\subset K_{i}. One has to show as well that this approximation can be chosen so that RiδiRjδjR_{i}^{\delta_{i}}\cap R_{j}^{\delta_{j}}\neq\emptyset for distinct i,j[n]i,j\in[n] and for every choice of signs δi,δj\delta_{i},\delta_{j}. As in (3.1), Prob(FU)\mathrm{Prob}(F\in U) is a linear function of the underlying distributions μi\mu_{i}, and this linear function equals 2n/2n2n/2^{n} on the product of two-ray distributions. Therefore this linear function equals 2n/2n2n/2^{n} on any convex combination of products of two-ray distributions and consequently Prob(FU)\mathrm{Prob}(F\in U) must be equal to 2n/2n2n/2^{n} on the product of the μi\mu_{i}s.


Acknowledgements. The first author was partially supported by Hungarian National Research grants 131529, 131696, and KKP-133819. The last author was partially supported by grants 965/15, 863/15, and 1919/19 from the Israel Science Foundation and by a grant from the Israeli Ministry of Science and Technology.


References

  • [1] E. W. Anderson. The treatment of navigational errors. Journal of Navigation, 5(2):103–124, 1952.
  • [2] J. Charnley. J. E. D. Williams. Journal of Navigation, 46(2):290–292, 1993. An obituary.
  • [3] I. Cook. Random cocked hats. Journal of Navigation, 46(1):132–135, 1993. A response by J. E. D. Williams appears on pages 135–137.
  • [4] C. H. Cotter. The cocked hat. Journal of Navigation, 14(2):223–230, 1961.
  • [5] H. E. Daniels. The theory of position finding. Journal of the Royal Statistical Society, Series B (Methodological), 13(2):186–207, 1951.
  • [6] I.C.B. Dear and P. Kemp, editors. The Oxford Companion to Ships and the Sea. Oxford University Press, 2nd edition, 2006.
  • [7] M. Denny. The Science of Navigation: From Dead Reckoning to GPS. The Johns Hopkins University Press, 2012.
  • [8] His Majesty’s Navigation School. Admiralty Navigation Manual, volume 3. His Magesty’s Stationary Office, 1938.
  • [9] R. G. Stansfield. Statistical theory of D.F. fixing. Journal of the Institution of Electrical Engineers Part IIIA: Radiocommunication, 94(15):762–770, 1947.
  • [10] R. G. Stuart. Probabilities in a Gaussian cocked hat. Journal of Navigation, 72(6):1496–1512, 2019.
  • [11] The Open University. M245 Probability and statistic: Item 15, Findings One’s Bearings. UK, 1984. A video that is part of course materials, available at https:// www.open.ac.uk/library/digital-archive/video:FOUM190W.
  • [12] J. E. D. Williams. The cocked hat. Journal of Navigation, 44(2):269–271, 1991.

Imre Bárány
Rényi Institute of Mathematics,
13-15 Reáltanoda Street, Budapest, 1053 Hungary
barany.imre@renyi.hu
and
Department of Mathematics
University College London
Gower Street, London, WC1E 6BT, UK

William Steiger
Department of Computer Science, Rutgers University
110 Frelinghuysen Road,Piscataway, NJ 08854-8019, USA
steiger@cs.rutgers.edu

Sivan Toledo
Blavatnik School of Computer Science, Tel Aviv University,
Tel Aviv 6997801, Israel
stoledo@mail.tau.ac.il