Particles \(P\), of mass \(2\), and \(Q\), of mass \(1\), move along a line. Their distances from a fixed point are \(x_1\) and \(x_2\), respectively where \(x_2>x_1\,\). Each particle is subject to a repulsive force from the other of magnitude \(\displaystyle {2 \over z^3}\), where \(z = x_2-x_1 \,\). Initially, \(x_1=0\), \(x_2 = 1\), \(Q\) is at rest and \(P\) moves towards \(Q\) with speed 1. Show that \(z\) obeys the equation \(\displaystyle {\mathrm{d}^2 z \over \mathrm{d}t^2} = {3 \over z^3}\). By first writing \(\displaystyle {\mathrm{d}^2 z \over \mathrm{d}t^2} = v {\mathrm{d}v \over \mathrm{d}z} \,\), where \(\displaystyle v={\mathrm{d}z \over \mathrm{d}t}\,\), show that \(z=\sqrt{4t^2-2t+1}\,\). By considering the equation satisfied by \(2x_1+x_2\,\), find \(x_1\) and \(x_2\) in terms of \(t \,\).
Solution: \begin{align*} \text{N2}: && 2\ddot{x}_1 &= -\frac{2}{(x_2-x_1)^3}\\ \text{N2}: && \ddot{x}_2 &= \frac{2}{(x_2-x_1)^3}\\ \Rightarrow && \ddot{x}_2 - \ddot{x}_1 &= \frac{3}{(x_1-x_2)^3} \\ \Rightarrow && \frac{\d^2 z}{\d t^2} &= \frac{3}{z^3} \\ \Rightarrow && v \frac{\d v}{\d z} &= \frac{3}{z^3} \\ \Rightarrow && \int v \d v &= \int \frac{3}{z^3} \d z \\ \Rightarrow && \frac{v^2}{2} &= -\frac{3}{2}z^{-2} + C \\ \Rightarrow && v^2 &= -3 z^{-2} + C' \\ t=0,z=1,v=-1: && 1 &= -3+C \Rightarrow C = 4 \\ \Rightarrow && \frac{\d z}{\d t} &= -\sqrt{4-3z^{-2}} \\ \Rightarrow && \int \d t &= -\int \frac{1}{\sqrt{4-3z^{-2}}} \d z \\ \Rightarrow && t &= \int \frac{z}{\sqrt{4z^2-3}} \d z \\ \Rightarrow && t &= -\frac14\sqrt{4z^2-3} + C \\ t=0, z = 1: && 0 &= -\frac14+C \\ \Rightarrow && C &= \frac14\\ \Rightarrow && 4t -1 &= -\sqrt{4z^2-3} \\ \Rightarrow && 16t^2+1-8t &= 4z^2-3 \\ \Rightarrow && z &= \sqrt{4t^2-2t+1} \end{align*} \begin{align*} && 2\ddot{x}_1 + \ddot{x}_2 &= 0 \\ \Rightarrow && 2x_1+x_2 &= At + B \\ t = 0, v = -1: && 2x_1+x_2 &= -t+1 \\ \\ \Rightarrow && x_2-x_1 &= \sqrt{4t^2-2t+1}\\ && 2x_1+x_2 &= 1-t \\ \Rightarrow && x_1 &= \frac13 \left (1-t-\sqrt{4t^2-2t+1} \right) \\ && x_2 &= \frac13(1-t + \sqrt{4t^2-2t+1}) \end{align*} This method of considering the relative position and considering the motion of the centre of mass is extremely common for solving systems of particles problems.
A team of \(m\) players, numbered from \(1\) to \(m\), puts on a set of a \(m\) shirts, similarly numbered from \(1\) to \(m\). The players change in a hurry, so that the shirts are assigned to them randomly, one to each player. Let \(C_i\) be the random variable that takes the value \(1\) if player \(i\) is wearing shirt \(i\), and 0 otherwise. Show that \(\mathrm{E}\left(C_1\right)={1 \over m}\) and find \(\var \left(C_1\right)\) and \(\mathrm{Cov}\left(C_1 \, , \; C_2 \right) \,\). Let \(\, N = C_1 + C_2 + \cdots + C_m \,\) be the random variable whose value is the number of players who are wearing the correct shirt. Show that \(\mathrm{E}\left(N\right)= \var \left(N\right) = 1 \,\). Explain why a Normal approximation to \(N\) is not likely to be appropriate for any \(m\), but that a Poisson approximation might be reasonable. In the case \(m = 4\), find, by listing equally likely possibilities or otherwise, the probability that no player is wearing the correct shirt and verify that an appropriate Poisson approximation to \(N\) gives this probability with a relative error of about \(2\%\). [Use \(\e \approx 2\frac{72}{100} \,\).]
Solution: There are \(m!\) different ways of assigning the shirts, and in \((m-1)!\) of them player \(1\) gets their own shirt, ie \(\mathbb{E}(C_1) = \mathbb{P}(\text{player }1\text{ gets own shirt}) = \frac{(m-1)!}{m!} = \frac{1}{m}\). \(\var(C_1) = \mathbb{E}(C_1^2) - [\mathbb{E}(C_1)]^2 = \frac{1}{m} - \frac{1}{m^2} = \frac{m-1}{m^2}\). If we have two players, there are \((m-2)!\) ways they both get their own shirts, therefore \(\textrm{Cov}(C_1,C_2) = \mathbb{E}(C_1C_2) - \mathbb{E}(C_1)\mathbb{E}(C_2) = \frac{(m-2)!}{m!} - \frac{1}{m^2} = \frac{1}{m(m-1)} - \frac{1}{m^2} = \frac{m-m+1}{m^2(m-1)} = \frac{1}{m^2(m-1)}\). \begin{align*} \mathbb{E}(N) &= \mathbb{E}(C_1 + C_2 + \cdots + C_m) \\ &= \mathbb{E}(C_1) + \mathbb{E}(C_2) + \cdots + \mathbb{E}(C_m) \\ &= \frac{1}{m} + \frac{1}{m} +\cdots+ \frac1m \\ &= 1 \\ \\ \var(N) &= \sum_{r=1}^m \var(C_r) + 2\sum_{r=1}^{m-1} \sum_{s=2}^{m} \textrm{Cov}(C_r,C_s) \\ &= m \frac{m-1}{m^2} + 2 \frac{m(m-1)}{2}\frac{1}{m^2(m-1)} \\ &=\frac{m-1}{m} + \frac{1}{m} \\ &= 1 \end{align*} If we were to take a normal approximation, we would want to take \(N(1,1)\), but this would say things like \(-1\) is as likely as \(3\) shirts being correct, which is clearly a bad model. A Poisson is much more likely to be a sensible model as they have the same mean and variance as the parameter, and if \(m\) is large, the covariance between shirts is going to be very small, so it will appear similar to random events occurring. We can have \begin{align*} BADC \\ BCDA \\ BDAC \\ CADB \\ CDAB\\ CDBA \\ DABC\\ DCAB \\ DCBA \end{align*} Ie \(\frac{9}{24}\) ways to have no player wearing their own shirt with \(4\) players. \(Po(1)\) would say this probability is \(e^{-1}\), giving a relative error of: \begin{align*} \frac{e^{-1}-\frac{9}{24}}{\frac9{24}} &\approx \frac{\frac{100}{272} - \frac{9}{24}}{\frac9{24}} \\ &= -\frac{1}{51} \\ &\approx -2\% \end{align*}
A men's endurance competition has an unlimited number of rounds. In each round, a competitor has, independently, a probability \(p\) of making it through the round; otherwise, he fails the round. Once a competitor fails a round, he drops out of the competition; before he drops out, he takes part in every round. The grand prize is awarded to any competitor who makes it through a round which all the other remaining competitors fail; if all the remaining competitors fail at the same round the grand prize is not awarded. If the competition begins with three competitors, find the probability that:
Solution:
\textit{In this question, \(\Phi(z)\) is the cumulative distribution function of a standard normal random variable.} A random variable is known to have a Normal distribution with mean \(\mu\) and standard deviation either \(\sigma_0\) or \(\sigma_1\), where \(\sigma_0 < \sigma_1\,\). The mean, \(\overline{X}\), of a random sample of \(n\) values of \(X\) is to be used to test the hypothesis \(\mathrm{H}_0: \sigma = \sigma_0\) against the alternative \(\mathrm{H}_1: \sigma = \sigma_1\,\). Explain carefully why it is appropriate to use a two sided test of the form: accept \(\mathrm{H}_0\) if \phantom{} \(\mu - c < \overline{X} < \mu+c\,\), otherwise accept \(\mathrm{H}_1\). Given that the probability of accepting \(\mathrm{H}_1\) when \(\mathrm{H}_0\) is true is \(\alpha\), determine \(c\) in terms of \(n\), \(\sigma_0\) and \(z_{\alpha}\), where \(z_\alpha \) is defined by \(\ds\Phi(z_{\alpha}) = 1 - \tfrac{1}{2}\alpha\). The probability of accepting \(\mathrm{H}_0\) when \(\mathrm{H}_1\) is true is denoted by \(\beta\). Show that \(\beta\) is independent of \(n\). Given that \(\Phi(1.960)\approx 0.975\) and that \(\Phi(0.063) \approx 0.525\,\), determine, approximately, the minimum value of \(\ds \frac{\sigma_1}{\sigma_0}\) if \(\alpha\) and \(\beta\) are both to be less than \(0.05\,\).