4.3 HW 3

  4.3.1 Problems listing
  4.3.2 Problem 1(a) (problem 5.2.3)
  4.3.3 Problem 1(b) (problem 5.2.4(iv))
  4.3.4 Problem 1(c) (problem 5.2.5)
  4.3.5 Problem 1(d) (problem 5.3.2)
  4.3.6 Problem 2(a) (problem 5.3.5)
  4.3.7 Problem 2(b) (problem 5.3.6)
  4.3.8 Problem 2(c)
  4.3.9 Problem 3
  4.3.10 Problem 4
  4.3.11 key solution for HW 3

4.3.1 Problems listing

PDF

PDF (letter size)
PDF (legal size)

4.3.2 Problem 1(a) (problem 5.2.3)

Solve for \(x\) and \(y\) given \[ \frac {2+3i}{6+7i}+\frac {2}{x+iy}=2+9i \] solution

Let \(z=x+iy\) be the complex number to solve for. The above becomes\begin {align*} \frac {2}{z} & =2+9i-\frac {2+3i}{6+7i}\\ & =2+9i-\frac {\left (2+3i\right ) \left (6-7i\right ) }{\left ( 6+7i\right ) \left (6-7i\right ) }\\ & =2+9i-\frac {12-14i+18i+21}{36+49}\\ & =2+9i-\frac {33+4i}{85}\\ & =\frac {85\left (2+9i\right ) -33-4i}{85}\\ & =\frac {170+765i-33-4i}{85}\\ & =\frac {137+761i}{85} \end {align*}

Therefore\begin {align*} \frac {2}{z} & =\frac {137+761i}{85}\\ z & =\frac {170}{137+761i}\\ & =\frac {170\left (137-761i\right ) }{\left (137+761i\right ) \left ( 137-761i\right ) }\\ & =\frac {23290-129370i}{597\,890}\\ & =\frac {23290}{597890}-\frac {129370}{597890}i\\ & =\frac {137}{3517}-\frac {761}{3517}i \end {align*}

But \(z=x+iy\). Hence\[ x+iy=\frac {137}{3517}-\frac {761}{3517}i \] Comparing real and imaginary parts shows that\begin {align*} x & =\frac {137}{3517}\\ y & =-\frac {761}{3517} \end {align*}

4.3.3 Problem 1(b) (problem 5.2.4(iv))

Find the real part, imaginary part, modulus, complex conjugate and inverse of the following (iv) \(\frac {1+\sqrt {2}i}{1-\sqrt {3}i}\)

solution

\begin {align*} \,z & =\frac {1+\sqrt {2}i}{1-\sqrt {3}i}\\ & =\frac {\left (1+\sqrt {2}i\right ) \left (1+\sqrt {3}i\right ) }{\left ( 1-\sqrt {3}i\right ) \left (1+\sqrt {3}i\right ) }\\ & =\frac {1+\sqrt {3}i+\sqrt {2}i-\sqrt {2}\sqrt {3}}{4}\\ & =\frac {1-\sqrt {6}}{4}+i\frac {\sqrt {3}+\sqrt {2}}{4} \end {align*}

Hence the real part is \(\frac {1-\sqrt {6}}{4}\) and the imaginary part is \(\frac {\sqrt {3}+\sqrt {2}}{4}\). Therefore we can now write\begin {align*} z & =x+iy\\ & =\left (\frac {1-\sqrt {6}}{4}\right ) +i\left (\frac {\sqrt {3}+\sqrt {2}}{4}\right ) \end {align*}

The modulus is \begin {align*} \left \vert z\right \vert & =\sqrt {x^{2}+y^{2}}\\ & =\sqrt {\left (\frac {1-\sqrt {6}}{4}\right ) ^{2}+\left (\frac {\sqrt {3}+\sqrt {2}}{4}\right ) ^{2}}\\ & =\sqrt {\frac {7}{16}-\frac {1}{8}\sqrt {6}+\frac {1}{8}\sqrt {6}+\frac {5}{16}}\\ & =\sqrt {\frac {3}{4}} \end {align*}

The complex conjugate of \(z\) is \(z^{\ast }\). Hence\begin {align*} z^{\ast } & =x-iy\\ & =\left (\frac {1-\sqrt {6}}{4}\right ) -i\left (\frac {\sqrt {3}+\sqrt {2}}{4}\right ) \end {align*}

The inverse is\begin {align*} \frac {1}{z} & =\frac {z^{\ast }}{zz^{\ast }}\\ & =\frac {z^{\ast }}{\left \vert z\right \vert ^{2}}\\ & =\frac {\left (\frac {1-\sqrt {6}}{4}\right ) -i\left (\frac {\sqrt {3}+\sqrt {2}}{4}\right ) }{\frac {3}{4}}\\ & =\frac {1-\sqrt {6}}{3}-i\frac {\sqrt {3}+\sqrt {2}}{3} \end {align*}

4.3.4 Problem 1(c) (problem 5.2.5)

Show that a polynomial with real coefficients has only real roots or complex roots that come in complex conjugate pairs.

solution

Let \[ p\relax (z) =a_{0}+a_{1}z+a_{2}z^{2}+\cdots +a_{n}z^{n}\] be polynomial in \(z\) where \(a_{i}\) are all real. We just need to show now that if \(\lambda \) is a root, then its complex conjugate \(\lambda ^{\ast }\) must also a root. If the root happened to be real, then its complex conjugate is itself. Hence nothing to do in this case. We only need to worry about the case when the root is complex and show that its complex conjugate must also be root.

Assuming \(\lambda \) is a root, then by definition of a root we have\begin {align} p\left (\lambda \right ) & =0\nonumber \\ & =a_{0}+a_{1}\lambda +a_{2}\lambda ^{2}+\cdots +a_{n}\lambda ^{n}\nonumber \\ & ={\displaystyle \sum \limits _{k=0}^{n}} a_{k}\lambda ^{k} \tag {1} \end {align}

Therefore, replacing \(\lambda \) by \(\lambda ^{\ast }\) on both sides of (1) gives\[ p\left (\lambda ^{\ast }\right ) ={\displaystyle \sum \limits _{k=0}^{n}} a_{k}\left (\lambda ^{\ast }\right ) ^{k}\] But \(\left (\lambda ^{\ast }\right ) ^{k}=\left (\lambda ^{k}\right ) ^{\ast }\) from complex numbers properties (equation 5.2.20 in book). The above becomes\[ p\left (\lambda ^{\ast }\right ) ={\displaystyle \sum \limits _{k=0}^{n}} a_{k}\left (\lambda ^{k}\right ) ^{\ast }\] Since \(a_{k}\) are real coefficients, then \(a_{k}^{\ast }=a_{k}\) and the above can be written as\[ p\left (\lambda ^{\ast }\right ) ={\displaystyle \sum \limits _{k=0}^{n}} \left (a_{k}\lambda ^{k}\right ) ^{\ast }\] Using property that \(A^{\ast }B^{\ast }=\left (AB\right ) ^{\ast }\) where \(A=a_{k},B=\lambda ^{k}\) in the above. Now we can move the complex conjugate outside the sum, using property that \(A^{\ast }+B^{\ast }=\left (A+B\right ) ^{\ast }\). Hence the above becomes\[ p\left (\lambda ^{\ast }\right ) =\left ( {\displaystyle \sum \limits _{k=0}^{n}} a_{k}\lambda ^{k}\right ) ^{\ast }\] But from (1), we know that \({\displaystyle \sum \limits _{k=0}^{n}} a_{k}\lambda ^{k}=0\), this is because \(\lambda \) is assumed to be a root. Therefore the above gives\begin {align*} p\left (\lambda ^{\ast }\right ) & =0^{\ast }\\ & =0 \end {align*}

The above shows that \(\lambda ^{\ast }\) is also a root if \(\lambda \) is a root. Therefore, the root can be either real, or complex. If the root is complex, its complex conjugate is also a root. A real root is just special case of complex root. QED.

4.3.5 Problem 1(d) (problem 5.3.2)

   4.3.5.1 First pair
   4.3.5.2 Second pair

For the following pairs of numbers, give their polar form, their complex conjugate, moduli, product, the quotient \(\frac {z_{1}}{z_{2}}\), and the complex conjugate of the quotient\begin {align*} z_{1} & =\frac {1+i}{\sqrt {2}}\qquad z_{2}=\sqrt {3}-i\\ z_{1} & =\frac {3+4i}{3-4i}\qquad z_{2}=\left (\frac {1+2i}{1-3i}\right ) ^{2} \end {align*}

solution

4.3.5.1 First pair

\[ z_{1}=\frac {1+i}{\sqrt {2}}\qquad z_{2}=\sqrt {3}-i \] The polar form of \(z\) is \(re^{i\theta }\) where \(r=\left \vert z\right \vert \) and \(\theta =\arctan \left (\frac {y}{x}\right ) \) when \(z=x+iy\). The first step is to write \(z=x+iy\)

For \(z_{1}\)

\begin {align*} z_{1} & =\frac {1+i}{\sqrt {2}}\\ & =\frac {1}{\sqrt {2}}+i\frac {1}{\sqrt {2}} \end {align*}

Hence \(x=\frac {1}{\sqrt {2}},y=\frac {1}{\sqrt {2}}\). Therefore \(\left \vert z_{1}\right \vert =\sqrt {x^{2}+y^{2}}=\sqrt {\frac {1}{2}+\frac {1}{2}}=1\). And \(\theta =\arctan \relax (1) =45^{0}\). Therefore in polar\begin {align*} z_{1} & =re^{i\theta }\\ & =e^{i\left (45^{0}\right ) }\\ & =e^{i\frac {\pi }{4}} \end {align*}

For \(z_{2}\) \[ z_{2}=\sqrt {3}-i \] Hence \(x=\sqrt {3},y=-1\). Therefore \(\left \vert z_{1}\right \vert =\sqrt {x^{2}+y^{2}}=\sqrt {3+1}=2\). And \(\theta =\arctan \left (\frac {-1}{\sqrt {3}}\right ) =-30^{0}\). Therefore in polar\begin {align*} z_{2} & =re^{i\theta }\\ & =2e^{i\left (-30^{0}\right ) }\\ & =2e^{-i\frac {\pi }{6}} \end {align*}

The complex conjugate is\begin {align*} z_{1}^{\ast } & =re^{-i\theta }\\ & =e^{-i\frac {\pi }{4}} \end {align*}

And\begin {align*} z_{2}^{\ast } & =re^{-i\theta }\\ & =2e^{i\frac {\pi }{6}} \end {align*}

And moduli is\begin {align*} \left \vert z_{1}\right \vert & =r\\ & =1 \end {align*}

And\begin {align*} \left \vert z_{2}\right \vert & =r\\ & =2 \end {align*}

And product\begin {align*} z_{1}z_{2} & =\left (r_{1}e^{i\theta _{1}}\right ) \left (r_{2}e^{i\theta _{2}}\right ) \\ & =r_{1}r_{2}e^{i\left (\theta _{1}+\theta _{2}\right ) } \end {align*}

But \(r_{1}=1,r_{2}=2,\theta _{1}=45^{0},\theta _{2}=-30^{0}\). The above becomes\begin {align*} z_{1}z_{2} & =2e^{i\left (45^{0}-30^{0}\right ) }\\ & =2e^{i\left (15^{0}\right ) }\\ & =2e^{i\frac {\pi }{12}} \end {align*}

And the quotient \(\frac {z_{1}}{z_{2}}\) is\begin {align*} \frac {z_{1}}{z_{2}} & =\frac {r_{1}e^{i\theta _{1}}}{r_{2}e^{i\theta _{2}}}\\ & =\frac {r_{1}}{r_{2}}e^{i\left (\theta _{1}-\theta _{2}\right ) } \end {align*}

But \(r_{1}=1,r_{2}=2,\theta _{1}=45^{0},\theta _{2}=-30^{0}\). The above becomes\begin {align*} \frac {z_{1}}{z_{2}} & =\frac {1}{2}e^{i\left (45^{0}+30^{0}\right ) }\\ & =\frac {1}{2}e^{i\left (75^{0}\right ) }\\ & =\frac {1}{2}e^{i\frac {5\pi }{12}} \end {align*}

And the complex conjugate of the quotient is\begin {align*} \left (\frac {z_{1}}{z_{2}}\right ) ^{\ast } & =\left (\frac {1}{2}e^{i\frac {5\pi }{12}}\right ) ^{\ast }\\ & =\frac {1}{2}e^{-i\frac {5\pi }{12}} \end {align*}

4.3.5.2 Second pair

\[ z_{1}=\frac {3+4i}{3-4i}\qquad z_{2}=\left (\frac {1+2i}{1-3i}\right ) ^{2}\] The polar form of \(z\) is \(re^{i\theta }\) where \(r=\left \vert z\right \vert \) and \(\theta =\arctan \left (\frac {y}{x}\right ) \) where \(z=x+iy\). Hence

For \(z_{1}\) \begin {align*} z_{1} & =\frac {3+4i}{3-4i}\\ & =\frac {\sqrt {3^{2}+4^{2}}e^{i\arctan \left (\frac {4}{3}\right ) }}{\sqrt {3^{2}+4^{2}}e^{i\arctan \left (-\frac {4}{3}\right ) }}\\ & =\frac {e^{i\arctan \left (\frac {4}{3}\right ) }}{e^{-i\arctan \left ( \frac {4}{3}\right ) }}\\ & =e^{i\arctan \left (\frac {4}{3}\right ) +\arctan \left (\frac {4}{3}\right ) }\\ & =e^{i\left (2\arctan \left (\frac {4}{3}\right ) \right ) }\\ & =e^{i\left (106.26^{0}\right ) } \end {align*}

For \(z_{2}\) \begin {align*} z_{2} & =\left (\frac {1+2i}{1-3i}\right ) ^{2}\\ & =\left (\frac {\left (1+2i\right ) \left (1+3i\right ) }{\left ( 1-3i\right ) \left (1+3i\right ) }\right ) ^{2}\\ & =\left (\frac {-5+5i}{10}\right ) ^{2}\\ & =\frac {25-25-50i}{100}\\ & =\frac {-1}{2}i \end {align*}

Hence \(x=0,y=-\frac {1}{2}\). Therefore \(\left \vert z_{1}\right \vert =\sqrt {x^{2}+y^{2}}=\sqrt {0+\frac {1}{4}}=\frac {1}{2}\). And \(\theta =\arctan \left (-\infty \right ) =-90^{0}\). Therefore in polar\begin {align*} z_{2} & =re^{i\theta }\\ & =\frac {1}{2}e^{i\left (-90^{0}\right ) }\\ & =\frac {1}{2}e^{-i\frac {\pi }{2}} \end {align*}

The complex conjugate is\begin {align*} z_{1}^{\ast } & =r_{1}e^{-i\theta _{1}}\\ & =e^{-i\left (2\arctan \left (\frac {4}{3}\right ) \right ) }\\ & =e^{i\left (-106.26^{0}\right ) } \end {align*}

And\begin {align*} z_{2}^{\ast } & =r_{2}e^{-i\theta _{2}}\\ & =\frac {1}{2}e^{i\frac {\pi }{2}} \end {align*}

And moduli is\begin {align*} \left \vert z_{1}\right \vert & =r_{1}\\ & =1 \end {align*}

And\begin {align*} \left \vert z_{2}\right \vert & =r_{2}\\ & =\frac {1}{2} \end {align*}

And product\begin {align*} z_{1}z_{2} & =\left (r_{1}e^{i\theta _{1}}\right ) \left (r_{2}e^{i\theta _{2}}\right ) \\ & =r_{1}r_{2}e^{i\left (\theta _{1}+\theta _{2}\right ) } \end {align*}

But \(r_{1}=1,r_{2}=\frac {1}{2},\theta _{1}=106.26^{0},\theta _{2}=-90^{0}\). The above becomes\begin {align*} z_{1}z_{2} & =\frac {1}{2}e^{i\left (106.26^{0}-90^{0}\right ) }\\ & =\frac {1}{2}e^{i\left (16.\,26^{0}\right ) } \end {align*}

And the quotient \(\frac {z_{1}}{z_{2}}\) is\begin {align*} \frac {z_{1}}{z_{2}} & =\frac {r_{1}e^{i\theta _{1}}}{r_{2}e^{i\theta _{2}}}\\ & =\frac {r_{1}}{r_{2}}e^{i\left (\theta _{1}-\theta _{2}\right ) } \end {align*}

But \(r_{1}=1,r_{2}=\frac {1}{2},\theta _{1}=106.26^{0},\theta _{2}=-90^{0}\). The above becomes\begin {align*} \frac {z_{1}}{z_{2}} & =2e^{i\left (106.26^{0}+90^{0}\right ) }\\ & =2e^{i\left (196.\,26^{0}\right ) } \end {align*}

And the complex conjugate of the quotient is\begin {align*} \left (\frac {z_{1}}{z_{2}}\right ) ^{\ast } & =\left (2e^{i\left ( 196.26^{0}\right ) }\right ) ^{\ast }\\ & =2e^{-i\left (196.26^{0}\right ) } \end {align*}

4.3.6 Problem 2(a) (problem 5.3.5)

Consider series \[ e^{i\theta }+e^{3i\theta }+\cdots +e^{\left (2n-1\right ) i\theta }\] Sum this geometric series, take the real and imaginary parts of both sides and show that \[ \cos \theta +\cos \left (3\theta \right ) +\cdots +\cos \left (\left ( 2n-1\right ) \theta \right ) =\frac {\sin \left (2n\theta \right ) }{2\sin \theta }\] And that a similar sum with sines adds up to \(\frac {\sin ^{2}\left ( n\theta \right ) }{\sin \theta }\)

solution

Let \begin {equation} S=e^{i\theta }+e^{3i\theta }+\cdots +e^{\left (2n-1\right ) i\theta }\tag {1} \end {equation} Then\begin {align} e^{2i\theta }S & =e^{2i\theta }\left (e^{i\theta }+e^{3i\theta }+\cdots +e^{\left (2n-1\right ) i\theta }\right ) \nonumber \\ & =e^{i3\theta }+e^{5i\theta }+\cdots +e^{\left (2n-1\right ) i\theta +2i\theta }\nonumber \\ & =e^{i3\theta }+e^{5i\theta }+\cdots +e^{\left (2n+1\right ) i\theta }\tag {2} \end {align}

Hence (2-1) gives\begin {align*} e^{2i\theta }S-S & =e^{\left (2n+1\right ) i\theta }-e^{i\theta }\\ S\left (e^{2i\theta }-1\right ) & =e^{\left (2n+1\right ) i\theta }-e^{i\theta }\\ S & =\frac {e^{\left (2n+1\right ) i\theta }-e^{i\theta }}{e^{2i\theta }-1} \end {align*}

Hence\begin {align*} S & =\frac {e^{i\theta }\left (e^{i2n\theta }-1\right ) }{e^{2i\theta }-1}\\ & =e^{i\theta }\frac {\left (e^{in\theta }\left (e^{in\theta }-e^{-in\theta }\right ) \right ) }{e^{i\theta }\left (e^{i\theta }-e^{-i\theta }\right ) }\\ & =\frac {e^{in\theta }\left (e^{in\theta }-e^{-in\theta }\right ) }{\left ( e^{i\theta }-e^{-i\theta }\right ) }\\ & =e^{in\theta }\frac {\left (e^{in\theta }-e^{-in\theta }\right ) }{\left ( e^{i\theta }-e^{-i\theta }\right ) }\\ & =e^{in\theta }\frac {\sin n\theta }{\sin \theta }\\ & =\cos \left (n\theta +i\sin n\theta \right ) \frac {\sin n\theta }{\sin \theta }\\ & =\frac {\cos \left (n\theta \right ) \sin \left (n\theta \right ) }{\sin \theta }+i\frac {\sin ^{2}\left (n\theta \right ) }{\sin \theta } \end {align*}

But \(\cos \left (n\theta \right ) \sin \left (n\theta \right ) =\frac {1}{2}\sin \left (2n\theta \right ) \,\). Therefore the above becomes\begin {equation} S=\frac {\sin \left (2n\theta \right ) }{2\sin \theta }+i\frac {\sin ^{2}\left ( n\theta \right ) }{\sin \theta }\nonumber \end {equation} Hence \begin {align} \operatorname {Re}\relax (S) & =\frac {\sin \left (2n\theta \right ) }{2\sin \theta }\tag {3}\\ \operatorname {Im}\relax (S) & =\frac {\sin ^{2}\left (n\theta \right ) }{\sin \theta }\tag {4} \end {align}

Now we look at the LHS. Since \(S=e^{i\theta }+e^{3i\theta }+\cdots +e^{\left ( 2n-1\right ) i\theta }\), then \begin {align} S & =\left (\cos \theta +i\sin \theta \right ) +\left (\cos 3\theta +i\sin 3\theta \right ) +\cdots +\left (\cos \left (2n-1\right ) \theta +i\sin \left ( 2n-1\right ) \theta \right ) \nonumber \\ & =\left (\cos \theta +\cos 3\theta +\cdots +\cos \left (2n-1\right ) \theta \right ) +i\left (\sin \theta +\sin 3\theta +\cdots +\sin \left ( 2n-1\right ) \theta \right ) \tag {5} \end {align}

Comparing (5) and (3,4) shows that\begin {align*} \cos \theta +\cos 3\theta +\cdots +\cos \left (2n-1\right ) \theta & =\operatorname {Re}\relax (S) \\ & =\frac {\sin \left (2n\theta \right ) }{2\sin \theta } \end {align*}

And\begin {align*} \sin \theta +\sin 3\theta +\cdots +\sin \left (2n-1\right ) \theta & =\operatorname {Im}\relax (S) \\ & =\frac {\sin ^{2}\left (n\theta \right ) }{\sin \theta } \end {align*}

Which is the result we are asked to show.

4.3.7 Problem 2(b) (problem 5.3.6)

   4.3.7.1 Part 1
   4.3.7.2 Part 2

(1) Consider De Moivre’s theorem, which states that \(\left (\cos \theta +i\sin \theta \right ) ^{n}=\cos n\theta +i\sin n\theta \). This follows from taking the \(nth\) power of both sides of Euler’s theorem. Find the formula for \(\cos 4\theta \) and \(\sin 4\theta \) in terms of \(\cos \theta \) and \(\sin \theta .\)

(2) Given \(e^{iA}e^{iB}=e^{i\left (A+B\right ) }\) deduce \(\cos \left ( A+B\right ) \) and \(\sin \left (A+B\right ) \)

solution

4.3.7.1 Part 1

Let \(n=4\), therefore, using De Moivre’s theorem gives\begin {equation} \left (\cos \theta +i\sin \theta \right ) ^{4}=\cos 4\theta +i\sin 4\theta \tag {1} \end {equation} We now expand the LHS of the above directly as follows\begin {equation} \left (\cos \theta +i\sin \theta \right ) ^{4}=\left (\cos \theta +i\sin \theta \right ) ^{2}\left (\cos \theta +i\sin \theta \right ) ^{2} \tag {2} \end {equation} But \[ \left (\cos \theta +i\sin \theta \right ) ^{2}=\cos ^{2}\theta -\sin ^{2}\theta +2i\cos \theta \sin \theta \] Substituting the above into (2) gives\begin {align*} \left (\cos \theta +i\sin \theta \right ) ^{4} & =\left (\cos ^{2}\theta -\sin ^{2}\theta +2i\cos \theta \sin \theta \right ) \left (\cos ^{2}\theta -\sin ^{2}\theta +2i\cos \theta \sin \theta \right ) \\ & =\cos ^{2}\theta \left (\cos ^{2}\theta -\sin ^{2}\theta +2i\cos \theta \sin \theta \right ) \\ & -\sin ^{2}\theta \left (\cos ^{2}\theta -\sin ^{2}\theta +2i\cos \theta \sin \theta \right ) \\ & +2i\cos \theta \sin \theta \left (\cos ^{2}\theta -\sin ^{2}\theta +2i\cos \theta \sin \theta \right ) \end {align*}

Expanding the RHS above more, then the above becomes\begin {align*} \left (\cos \theta +i\sin \theta \right ) ^{4} & =\left (\cos ^{4}\theta -\cos ^{2}\theta \sin ^{2}\theta +2i\cos ^{3}\theta \sin \theta \right ) \\ & -\left (\sin ^{2}\theta \cos ^{2}\theta -\sin ^{4}\theta +2i\cos \theta \sin ^{3}\theta \right ) \\ & +\left (2i\cos ^{3}\theta \sin \theta -2i\cos \theta \sin ^{3}\theta -4\cos ^{2}\theta \sin ^{2}\theta \right ) \end {align*}

Simplifying gives\begin {align} \left (\cos \theta +i\sin \theta \right ) ^{4} & =\cos ^{4}\theta -6\cos ^{2}\theta \sin ^{2}\theta +4i\cos ^{3}\theta \sin \theta +\sin ^{4}\theta -4i\cos \theta \sin ^{3}\theta \nonumber \\ & =\left (\cos ^{4}\theta +\sin ^{4}\theta -6\cos ^{2}\theta \sin ^{2}\theta \right ) +i\left (4\cos ^{3}\theta \sin \theta -4\cos \theta \sin ^{3}\theta \right ) \tag {3} \end {align}

Comparing the real and imaginary parts of (3) with the real and imaginary parts of (1)  shows that\begin {align*} \cos 4\theta & =\cos ^{4}\theta +\sin ^{4}\theta -6\cos ^{2}\theta \sin ^{2}\theta \\ \sin 4\theta & =4\cos ^{3}\theta \sin \theta -4\cos \theta \sin ^{3}\theta \end {align*}

4.3.7.2 Part 2

Given \[ e^{iA}e^{iB}=e^{i\left (A+B\right ) }\] Applying Euler’s formula \(e^{ix}=\cos x+i\sin x\), on both sides of the above results in\begin {align*} \left (\cos A+i\sin A\right ) \left (\cos B+i\sin B\right ) & =\cos \left ( A+B\right ) +i\sin \left (A+B\right ) \\ \cos A\cos B+i\cos A\sin B+i\sin A\cos B-\sin B\sin A & =\cos \left ( A+B\right ) +i\sin \left (A+B\right ) \\ \left (\cos A\cos B-\sin B\sin A\right ) +i\left (\cos A\sin B+\sin A\cos B\right ) & =\cos \left (A+B\right ) +i\sin \left (A+B\right ) \end {align*}

Comparing the real parts and the imaginary parts in the above shows that\[ \cos A\cos B-\sin B\sin A=\cos \left (A+B\right ) \] And\[ \cos A\sin B+\sin A\cos B=\sin \left (A+B\right ) \]

4.3.8 Problem 2(c)

Find \(\int _{0}^{\infty }xe^{-ax}\cos \left (kx\right ) dx\) using Euler’s formula.

solution

Let \[ I=\int _{0}^{\infty }xe^{-ax}\cos \left (kx\right ) dx \] Then, we replace \(\cos \left (kx\right ) \) by \(e^{ikx}\), evaluate the integral, and then take the real part of the result. Therefore\begin {align*} I & =\operatorname {Re}\left (\int _{0}^{\infty }xe^{-ax}e^{ikx}dx\right ) \\ & =\operatorname {Re}\left (\int _{0}^{\infty }xe^{x\left (-a+ik\right ) }dx\right ) \end {align*}

Integration by parts. Let \(u=x,du=dx\) and \(dv=e^{x\left (-a+ik\right ) },v=\frac {e^{x\left (-a+ik\right ) }}{-a+ik}\). The above now becomes\begin {align} I & =\operatorname {Re}\left (\left . uv\right \vert _{0}^{\infty }-\int _{0}^{\infty }vdu\right ) \nonumber \\ & =\operatorname {Re}\left (\frac {1}{-a+ik}\left . xe^{x\left (-a+ik\right ) }\right \vert _{0}^{\infty }-\int _{0}^{\infty }\frac {e^{x\left (-a+ik\right ) }}{-a+ik}dx\right ) \tag {1} \end {align}

But \[ \left . xe^{x\left (-a+ik\right ) }\right \vert _{0}^{\infty }=0 \] With the assumption that \(\operatorname {Re}\relax (a) >0\). To see this more clearly, let us write \(e^{x\left (-a+ik\right ) }=e^{-ax}e^{ikx}\). \(e^{ikx}\) is bounded since it is a complex exponential. So the contribution comes from \(e^{-ax}\). Hence when \(a>0\), and \(x\rightarrow \infty \) then the exponential will go to zero, and the whole term \(xe^{x\left (-a+ik\right ) }\rightarrow 0\), even though \(x\rightarrow \infty \), since exponential subdues any polynomial order. When \(x=0\), it is clear that \(xe^{x\left (-a+ik\right ) }=0\). Therefore (1) now simplifies to

\begin {align*} I & =\operatorname {Re}\left (-\int _{0}^{\infty }\frac {e^{x\left ( -a+ik\right ) }}{-a+ik}dx\right ) \\ & =\operatorname {Re}\left (-\frac {1}{-a+ik}\int _{0}^{\infty }e^{x\left ( -a+ik\right ) }dx\right ) \\ & =\operatorname {Re}\left (-\frac {1}{-a+ik}\left . \frac {e^{x\left ( -a+ik\right ) }}{-a+ik}\right \vert _{0}^{\infty }\right ) \\ & =\operatorname {Re}\left (-\frac {1}{\left (-a+ik\right ) ^{2}}\left . e^{x\left (-a+ik\right ) }\right \vert _{0}^{\infty }\right ) \end {align*}

But \(\left . e^{x\left (-a+ik\right ) }\right \vert _{0}^{\infty }=0-1=-1\). The above becomes\begin {align*} I & =\operatorname {Re}\left (\frac {1}{\left (-a+ik\right ) ^{2}}\right ) \\ & =\operatorname {Re}\left (\frac {1}{a^{2}-k^{2}-2aik}\right ) \\ & =\operatorname {Re}\left (\frac {\left (a^{2}-k^{2}+2aik\right ) }{\left ( a^{2}-k^{2}-2aik\right ) \left (a^{2}-k^{2}+2aik\right ) }\right ) \\ & =\operatorname {Re}\left (\frac {a^{2}-k^{2}+2aik}{\left (a^{2}-k^{2}\right ) ^{2}+4a^{2}k^{2}}\right ) \\ & =\operatorname {Re}\left (\frac {a^{2}-k^{2}}{\left (a^{2}-k^{2}\right ) ^{2}+4a^{2}k^{2}}+i\frac {2ak}{\left (a^{2}-k^{2}\right ) ^{2}+4a^{2}k^{2}}\right ) \\ & =\frac {a^{2}-k^{2}}{\left (a^{2}-k^{2}\right ) ^{2}+4a^{2}k^{2}} \end {align*}

Hence\begin {align*} \int _{0}^{\infty }xe^{-ax}\cos \left (kx\right ) dx & =\frac {a^{2}-k^{2}}{\left (a^{2}-k^{2}\right ) ^{2}+4a^{2}k^{2}}\\ & =\frac {a^{2}-k^{2}}{a^{4}+k^{4}-2a^{2}k^{2}+4a^{2}k^{2}}\\ & =\frac {a^{2}-k^{2}}{a^{4}+k^{4}+2a^{2}k^{2}}\\ & =\frac {a^{2}-k^{2}}{\left (a^{2}+k^{2}\right ) ^{2}}\qquad a>0 \end {align*}

4.3.9 Problem 3

   4.3.9.1 Appendix

Given the intensity pattern for the \(N\)-slit interference with separation \(d\) between adjacent slits, show that the pattern becomes that for the single-slit diffraction with slit width \(a\) when \(d\) goes to zero but with a fixed value of \(Nd=a\). (10 points)

Solution

Short version: In this version, The result for \(N\) slit \(\bar {I}_{N}\left ( \theta \right ) \) will be used as given in lecture notes without deriving it again, and will also use the single slit \(\,\bar {I}_{1}\left (\theta \right ) \) from the lecture notes, then show that \(\bar {I}_{N}\left (\theta \right ) \) becomes \(\,\bar {I}_{1}\left (\theta \right ) \) as \(d\rightarrow 0\) but with \(Nd=a\).

Here \(\bar {I}_{N}\left (\theta \right ) \) is the average intensity for \(N\) slits at location on the screen at angle \(\theta \) and similarly \(\bar {I}_{1}\left (\theta \right ) \) is the average intensity for one slit at same location on the screen at angle \(\theta .\) From lecture notes (lecture 3, pages 6,7) we have the expressions for \(\bar {I}_{N}\left (\theta \right ) ,\bar {I}_{1}\left (\theta \right ) \) given as\begin {align} \bar {I}_{N}\left (\theta \right ) & =\bar {I}\relax (0) \left ( \frac {\sin \left (\frac {N\pi d\sin \theta }{\lambda }\right ) }{N\sin \left ( \frac {\pi d\sin \theta }{\lambda }\right ) }\right ) ^{2}\tag {1}\\ \bar {I}_{1}\left (\theta \right ) & =\bar {I}\relax (0) \left ( \frac {\sin \left (\frac {\pi a\sin \theta }{\lambda }\right ) }{\frac {\pi a\sin \theta }{\lambda }}\right ) ^{2} \tag {2} \end {align}

Now we need to show that (1) gives same result as (2) when \(d\) goes to zero in the limit, but with a fixed value of \(Nd=a\). Replacing \(Nd=a\) in the numerator of (1) and taking the limit gives\begin {align} \lim _{d\rightarrow 0}\bar {I}_{N}\left (\theta \right ) & =\bar {I}\left ( 0\right ) \left (\lim _{d\rightarrow 0}\frac {\sin \left (\frac {\pi a\sin \left ( \theta \right ) }{\lambda }\right ) }{N\sin \left (\frac {\pi d\sin \left ( \theta \right ) }{\lambda }\right ) }\right ) ^{2}\nonumber \\ & =\bar {I}\relax (0) \left (\frac {\sin \left (\frac {\pi a\sin \left ( \theta \right ) }{\lambda }\right ) }{N\lim _{d\rightarrow 0}\sin \left (\frac {\pi d\sin \left (\theta \right ) }{\lambda }\right ) }\right ) ^{2} \tag {3} \end {align}

But \begin {equation} \lim _{d\rightarrow 0}\sin \left (\frac {\pi d\sin \left (\theta \right ) }{\lambda }\right ) \approx \frac {\pi d\sin \left (\theta \right ) }{\lambda }+\cdots \tag {4} \end {equation} In the above we used that \(\lim _{d\rightarrow 0}\sin \left (\frac {\pi d\sin \left (\theta \right ) }{\lambda }\right ) \approx \frac {\pi d\sin \left ( \theta \right ) }{\lambda }\). This comes from Taylor series expansion of \(\sin \) function, for small angle approximation by keeping only the linear term in the Taylor series expansion since \(\sin \relax (x) =x-\frac {x^{3}}{3!}+\frac {x^{5}}{5!}-\cdots \).  

Substituting (4) back into (3) gives\[ \lim _{d\rightarrow 0}\bar {I}_{N}\left (\theta \right ) =\bar {I}\left ( 0\right ) \left (\frac {\sin \left (\frac {\pi a\sin \left (\theta \right ) }{\lambda }\right ) }{N\frac {\pi d\sin \left (\theta \right ) }{\lambda }}\right ) ^{2}\] But \(Nd=a\). The above simplifies to\begin {equation} \lim _{d\rightarrow 0}\bar {I}_{N}\left (\theta \right ) =\bar {I}\left ( 0\right ) \left (\frac {\sin \left (\frac {\pi a\sin \left (\theta \right ) }{\lambda }\right ) }{\frac {\pi a\sin \left (\theta \right ) }{\lambda }}\right ) ^{2}\tag {5} \end {equation} Comparing (5) with (2) shows that are the same. Hence\[ \lim _{d\rightarrow 0}\bar {I}_{N}\left (\theta \right ) =\bar {I}_{1}\left ( \theta \right ) \] Which is what we are asked to show.

4.3.9.1 Appendix

Here, the derivation of\begin {equation} \bar {I}_{N}\left (\theta \right ) =\bar {I}\relax (0) \left (\frac {\sin \left (\frac {N\pi d\sin \theta }{\lambda }\right ) }{N\sin \left (\frac {\pi d\sin \theta }{\lambda }\right ) }\right ) ^{2}\tag {1} \end {equation} is given. First, let us consider a slit located at \(y_{n}\) relative to the origin as show in the diagram below

pict
Figure 4.7:Geometry for slit at location \(y_n\)

Therefore \begin {align*} r_{n} & =\sqrt {L^{2}+\left (y-y_{n}\right ) ^{2}}\\ & =\sqrt {L^{2}+\left (y^{2}+y_{n}^{2}-2yy_{n}\right ) }\\ & =\sqrt {L^{2}+y^{2}+y_{n}^{2}-2yy_{n}}\\ & =\sqrt {\left (L^{2}+y^{2}\right ) \left (1+\frac {y_{n}^{2}-2yy_{n}}{\left (L^{2}+y^{2}\right ) }\right ) }\\ & =\sqrt {\left (L^{2}+y^{2}\right ) }\sqrt {1+\frac {y_{n}^{2}-2yy_{n}}{\left ( L^{2}+y^{2}\right ) }}\\ & =\sqrt {\left (L^{2}+y^{2}\right ) }\sqrt {1-\frac {2yy_{n}}{\left ( L^{2}+y^{2}\right ) }+\frac {y_{n}^{2}}{\left (L^{2}+y^{2}\right ) }} \end {align*}

Since \(y_{n}\) is very small compared to \(\left (L^{2}+y^{2}\right ) \) and it is also of order 2, then we can ignore the term \(\frac {y_{n}^{2}}{\left ( L^{2}+y^{2}\right ) }\) above, giving\begin {align*} r_{n} & \approx \sqrt {\left (L^{2}+y^{2}\right ) }\sqrt {1-\frac {2yy_{n}}{\left (L^{2}+y^{2}\right ) }}\\ & =\sqrt {\left (L^{2}+y^{2}\right ) }\left (1-\frac {2yy_{n}}{\left ( L^{2}+y^{2}\right ) }\right ) ^{\frac {1}{2}} \end {align*}

But since \(y_{n}\) is very small compared to \(\left (L^{2}+y^{2}\right ) \), then the term \(\frac {2yy_{n}}{\left (L^{2}+y^{2}\right ) }\) is very small. So we can use \(\left (1+x\right ) ^{p}=1+px\) and ignore higher order terms. Hence the above becomes\begin {align*} r_{n} & \approx \sqrt {\left (L^{2}+y^{2}\right ) }\left (1+\frac {1}{2}\left (\frac {-2yy_{n}}{\left (L^{2}+y^{2}\right ) }\right ) \right ) \\ & =\sqrt {\left (L^{2}+y^{2}\right ) }-\frac {yy_{n}}{\sqrt {\left (L^{2}+y^{2}\right ) }}\\ & =r-y_{n}\frac {y}{r} \end {align*}

But \(\frac {y}{r}=\sin \theta \), therefore\begin {equation} r_{n}=r-y_{n}\sin \theta \tag {2} \end {equation} The electric field \(E_{n}\) measured at point \(\left (L,y\right ) \) due to slit at \(y_{n}\) is \[ E_{n}=E_{0}\sin \left (kr_{n}-\omega t\right ) \] Where \(k\) is the wave number \(k=\frac {2\pi }{\lambda }\).  Therefore for \(N\) slits, the total \(E\) is\begin {align} E & ={\displaystyle \sum \limits _{n=1}^{N}} E_{n}\nonumber \\ & ={\displaystyle \sum \limits _{n=1}^{N}} E_{0}\sin \left (kr_{n}-\omega t\right ) \nonumber \\ & =E_{0}\left (\operatorname {Im}{\displaystyle \sum \limits _{n=1}^{N}} e^{i\left (kr_{n}-\omega t\right ) }\right ) \nonumber \\ & =E_{0}\left (\operatorname {Im}{\displaystyle \sum \limits _{n=1}^{N}} e^{i\left (kr_{N}-\omega t\right ) }e^{ik\left (r_{n}-r_{N}\right ) }\right ) \nonumber \\ & =E_{0}\left (\operatorname {Im}\left (e^{i\left (kr_{N}-\omega t\right ) }{\displaystyle \sum \limits _{n=1}^{N}} e^{ik\left (r_{n}-r_{N}\right ) }\right ) \right ) \tag {3} \end {align}

But \begin {align} r_{n}-r_{N} & =\left (r-y_{n}\sin \theta \right ) -\left (r-y_{N}\sin \theta \right ) \nonumber \\ & =r-y_{n}\sin \theta -r+y_{N}\sin \theta \nonumber \\ & =\left (y_{N}-y_{n}\right ) \sin \theta \nonumber \\ & =\left (Nd-nd\right ) \sin \theta \nonumber \\ & =\left (N-n\right ) d\sin \theta \tag {4} \end {align}

Substituting (4) in (3) gives\[ E=E_{0}\left (\operatorname {Im}\left (e^{i\left (kr_{N}-\omega t\right ) }{\displaystyle \sum \limits _{n=1}^{N}} e^{ik\left (N-n\right ) d\sin \theta }\right ) \right ) \] Let \(m=N-n\). When \(n=1\) then \(m=N-1\). When \(n=N\) then \(m=0\). The above now becomes\begin {align*} E & =E_{0}\left (\operatorname {Im}\left (e^{i\left (kr_{N}-\omega t\right ) }{\displaystyle \sum \limits _{m=N-1}^{0}} e^{ikmd\sin \theta }\right ) \right ) \\ & =E_{0}\left (\operatorname {Im}\left (e^{i\left (kr_{N}-\omega t\right ) }{\displaystyle \sum \limits _{m=0}^{N-1}} e^{ikmd\sin \theta }\right ) \right ) \\ & =E_{0}\left (\operatorname {Im}\left (e^{i\left (kr_{N}-\omega t\right ) }{\displaystyle \sum \limits _{n=0}^{N-1}} e^{iknd\sin \theta }\right ) \right ) \end {align*}

Let \(\phi =kd\sin \theta \). The above becomes\begin {equation} E=E_{0}\left (\operatorname {Im}\left (e^{i\left (kr_{N}-\omega t\right ) }{\displaystyle \sum \limits _{n=0}^{N-1}} e^{in\phi }\right ) \right ) \tag {5} \end {equation} But \begin {align} {\displaystyle \sum \limits _{n=0}^{N-1}} e^{in\phi } & =\frac {1-e^{iN\phi }}{1-e^{i\phi }}\nonumber \\ & =\frac {e^{i\frac {N}{2}\phi }\left (e^{-i\frac {N}{2}\phi }-e^{i\frac {N}{2}\phi }\right ) }{e^{i\frac {\phi }{2}}\left (e^{-i\frac {\phi }{2}}-e^{i\frac {\phi }{2}}\right ) }\nonumber \\ & =\frac {e^{i\frac {N}{2}\phi }}{e^{i\frac {\phi }{2}}}\frac {\left (e^{-i\frac {N}{2}\phi }-e^{i\frac {N}{2}\phi }\right ) }{\left (e^{-i\frac {\phi }{2}}-e^{i\frac {\phi }{2}}\right ) }\nonumber \\ & =\frac {-e^{i\frac {N}{2}\phi }}{-e^{i\frac {\phi }{2}}}\frac {\left ( e^{i\frac {N}{2}\phi }-e^{-i\frac {N}{2}\phi }\right ) }{\left (e^{i\frac {\phi }{2}}-e^{-i\frac {\phi }{2}}\right ) }\nonumber \\ & =\frac {e^{i\frac {N}{2}\phi }}{e^{i\frac {\phi }{2}}}\frac {\sin \left (\frac {N}{2}\phi \right ) }{\sin \left (\frac {\phi }{2}\right ) }\nonumber \\ & =e^{i\frac {\left (N-1\right ) \phi }{2}}\frac {\sin \left (\frac {N}{2}\phi \right ) }{\sin \left (\frac {\phi }{2}\right ) }\tag {6} \end {align}

Substituting (6) in (5) gives\begin {align} E & =E_{0}\left (\operatorname {Im}\left (e^{i\left (kr_{N}-\omega t\right ) }e^{i\frac {\left (N-1\right ) \phi }{2}}\frac {\sin \left (\frac {N}{2}\phi \right ) }{\sin \left (\frac {\phi }{2}\right ) }\right ) \right ) \nonumber \\ & =E_{0}\left (\operatorname {Im}\left (e^{i\left (kr_{N}-\omega t+\frac {\left (N-1\right ) \phi }{2}\right ) }\frac {\sin \left (\frac {N}{2}\phi \right ) }{\sin \left (\frac {\phi }{2}\right ) }\right ) \right ) \tag {7} \end {align}

Let \[ t_{0}=kr_{N}+\frac {\left (N-1\right ) \phi }{2}\] Substituting this in (7) gives\begin {align} E & =E_{0}\left (\operatorname {Im}\left (e^{i\omega \left (t_{0}-t\right ) }\right ) \right ) \nonumber \\ & =E_{0}\frac {\sin \left (\frac {N}{2}\phi \right ) }{\sin \left (\frac {\phi }{2}\right ) }\left (\operatorname {Im}e^{i\omega \left (t_{0}-t\right ) }\right ) \nonumber \\ & =E_{0}\frac {\sin \left (\frac {N}{2}\phi \right ) }{\sin \left (\frac {\phi }{2}\right ) }\sin \left (\omega \left (t_{0}-t\right ) \right ) \tag {8} \end {align}

The electric field intensity is \begin {align*} I & =c\varepsilon _{0}E^{2}\\ & =c\varepsilon _{0}E_{0}^{2}\frac {\sin ^{2}\left (\frac {N}{2}\phi \right ) }{\sin ^{2}\left (\frac {\phi }{2}\right ) }\sin ^{2}\left (\omega \left ( t_{0}-t\right ) \right ) \end {align*}

The time (period) averaged intensity is therefore\begin {align*} I_{av} & =\frac {1}{T}\int _{0}^{T}Idt\\ & =\frac {1}{T}c\varepsilon _{0}E_{0}^{2}\frac {\sin ^{2}\left (\frac {N}{2}\phi \right ) }{\sin ^{2}\left (\frac {\phi }{2}\right ) }\int _{0}^{T}\sin ^{2}\left (\omega \left (t_{0}-t\right ) \right ) dt\\ & =\frac {1}{2}c\varepsilon _{0}E_{0}^{2}\frac {\sin ^{2}\left (\frac {N}{2}\phi \right ) }{\sin ^{2}\left (\frac {\phi }{2}\right ) } \end {align*}

But \(\phi =kd\sin \theta \) and \(k=\frac {2\pi }{\lambda }\), then the above becomes\begin {align*} I\left (\theta \right ) _{av} & =\frac {1}{2}c\varepsilon _{0}E_{0}^{2}\frac {\sin ^{2}\left (\frac {N\pi d\sin \theta }{\lambda }\right ) }{\sin ^{2}\left (\frac {\pi d\sin \theta }{\lambda }\right ) }\\ & =\frac {1}{2}c\varepsilon _{0}E_{0}^{2}\left (\frac {\sin \left (\frac {N\pi d\sin \theta }{\lambda }\right ) }{\sin \left (\frac {\pi d\sin \theta }{\lambda }\right ) }\right ) ^{2} \end {align*}

At \(\theta =0\), we have\begin {align*} I\relax (0) _{av} & =\lim _{\theta \rightarrow 0}\frac {1}{2}c\varepsilon _{0}E_{0}^{2}\frac {\sin ^{2}\left (\frac {N\pi d\sin \theta }{\lambda }\right ) }{\sin ^{2}\left (\frac {\pi d\sin \theta }{\lambda }\right ) }\\ & =N^{2}\frac {1}{2}c\varepsilon _{0}E_{0}^{2} \end {align*}

Hence \begin {align*} \frac {I\left (\theta \right ) _{av}}{I\relax (0) _{av}} & =\frac {\frac {1}{2}c\varepsilon _{0}E_{0}^{2}\left (\frac {\sin \left (\frac {N\pi d\sin \theta }{\lambda }\right ) }{\sin \left (\frac {\pi d\sin \theta }{\lambda }\right ) }\right ) ^{2}}{N^{2}\frac {1}{2}c\varepsilon _{0}E_{0}^{2}}\\ & =\frac {1}{N^{2}}\left (\frac {\sin \left (\frac {N\pi d\sin \theta }{\lambda }\right ) }{\sin \left (\frac {\pi d\sin \theta }{\lambda }\right ) }\right ) ^{2}\\ & =\left (\frac {\sin \left (\frac {N\pi d\sin \theta }{\lambda }\right ) }{N\sin \left (\frac {\pi d\sin \theta }{\lambda }\right ) }\right ) ^{2} \end {align*}

Therefore\[ I\left (\theta \right ) _{av}=I\relax (0) _{av}\left (\frac {\sin \left (\frac {N\pi d\sin \theta }{\lambda }\right ) }{N\sin \left (\frac {\pi d\sin \theta }{\lambda }\right ) }\right ) ^{2}\] Which is the formula used in the earlier derivation.

4.3.10 Problem 4

   4.3.10.1 Part 1
   4.3.10.2 Part 2
   4.3.10.3 Part 3

(1) Find the roots \(z_{n}\left (n=1,2,\cdots ,N\right ) \) of the complex equation \(z^{N}=1\). (2) Find \(S_{N}=\sum _{n=1}^{N}z_{n}\) and give a geometric interpretation of the result. (3) Note that \(1-z^{N}=\left ( 1-z\right ) \left (1+z+z^{2}+\cdots +z^{N-1}\right ) \). Relate this result and the roots \(z_{n}\) to the conditions for destructive interference among \(N\) slits.

solution

4.3.10.1 Part 1

\begin {align*} Z^{N} & =1\\ Z & =1^{\frac {1}{N}} \end {align*}

But \(1=e^{i\left (2\pi \right ) }\) and the above becomes\begin {align*} Z & =\left (e^{i\left (2\pi \right ) }\right ) ^{\frac {1}{N}}\\ Z_{n} & =\left (\cos \left (2\pi +\left (2\pi \right ) n\right ) +i\sin \left (2\pi +\left (2\pi \right ) n\right ) \right ) ^{\frac {1}{N}}\qquad n=0,1,2,\cdots ,N-1 \end {align*}

Since \(\cos \) and \(\sin \) are periodic with period \(2\pi \). Using De Moivre’s theorem the above becomes\begin {align*} Z_{n} & =\left (\cos \left (\frac {2\pi }{N}+\frac {n}{N}\left (2\pi \right ) \right ) +i\sin \left (\frac {2\pi }{N}+\frac {n}{N}\left (2\pi \right ) \right ) \right ) \\ & =e^{i\left (\frac {2\pi }{N}+\frac {n}{N}\left (2\pi \right ) \right ) }\\ & =e^{i\left (\frac {2\pi \left (n+1\right ) }{N}\right ) }\qquad n=0,1,2,\cdots ,N-1 \end {align*}

Which is the same as\[ Z_{n}=e^{i\left (\frac {2\pi n}{N}\right ) }\qquad n=1,2,\cdots ,N \] For an example, let \(N=3\). Therefore we have 3 roots, given by \(n=1,2,3\). They are\begin {align*} Z_{1} & =e^{i\left (\frac {2\pi }{3}\right ) }=e^{i(120^{0})}\\ Z_{2} & =e^{i\left (\frac {2\pi \relax (2) }{3}\right ) }=e^{i\left ( \frac {4\pi }{3}\right ) }=e^{i\left (240^{0}\right ) }\\ Z_{3} & =e^{i\left (\frac {2\pi \relax (3) }{3}\right ) }=e^{i\left ( 2\pi \right ) }=e^{i360^{0}} \end {align*}

The roots are \(120^{0}\) degrees apart on the unit circle. First root has phase \(0^{0}\) (or \(360^{0}\)), second at \(120^{0}\) and the third at \(240^{0}\). There are only 3 unique roots, since after that, they repeat. Here is a diagram showing the roots for \(N=3\) for illustration. The root with phase \(0^{0}\) is the real root \(1\) since \(e^{i0^{0}}=1\), the other two roots are complex valued, and complex conjugate of each others.

pict
Figure 4.8:Roots of \(z^N\) for case of \(N=3\)

There are only 3 unique roots, since after moving around the unit circle once, the roots repeat.

4.3.10.2 Part 2

\begin {equation} S_{N}=\sum _{n=1}^{N}z_{n} \tag {1} \end {equation} It is assumed that \(z_{n}\) above are all the roots of \(Z^{N}\) from part(a), even though the problem did not say that. Hence all roots have same modulus. But differ by the phase as found in part 1.

Let \(z=x+iy=e^{i\theta }\) where \(r=\sqrt {x^{2}+y^{2}}\) and \(\theta =\arctan \left (\frac {y}{x}\right ) \). The above becomes\begin {align*} S_{N} & =z_{1}+z_{2}+\cdots +z_{N}\\ & =re^{i\theta _{1}}+re^{i\theta _{2}}+\cdots +re^{i\theta _{N}}\\ & =r\left (e^{i\theta _{1}}+e^{i\theta _{2}}+\cdots +e^{i\theta _{N}}\right ) \end {align*}

But \(r=1\), hence\[ S_{N}=e^{i\theta _{1}}+e^{i\theta _{2}}+\cdots +e^{i\theta _{N}}\] From part 1, we found that\begin {equation} \theta _{n}=\frac {2\pi n}{N}\qquad n=1,2,3,\cdots ,N \tag {2} \end {equation} Using (2) in (1), now the sum can be written as\begin {equation} S_{N}=\sum _{n=1}^{N}e^{i\frac {2\pi n}{N}} \tag {3} \end {equation} If \(N=1\), then the sum is just \(e^{i2\pi }=1\). But if \(N>1\) then to find the partial sum, let \begin {align} S_{N} & =e^{i\frac {2\pi }{N}}+e^{i\left (\frac {4\pi }{N}\right ) }+e^{i\frac {6\pi }{N}}+\cdots +e^{i\frac {N\left (2\pi \right ) }{N}}\tag {4}\\ e^{i\frac {2\pi }{N}}S_{N} & =e^{i\frac {4\pi }{N}}+e^{i\frac {6\pi }{N}}+\cdots +e^{i\frac {N\left (2\pi \right ) +2\pi }{N}} \tag {5} \end {align}

(4-5) gives\begin {align*} S_{N}-e^{i\frac {2\pi }{N}}S_{N} & =e^{i\frac {2\pi }{N}}-e^{i\frac {N\left ( 2\pi \right ) +2\pi }{N}}\\ S_{N}\left (1-e^{i\frac {2\pi }{N}}\right ) & =e^{i\frac {2\pi }{N}}-e^{i\frac {N\left (2\pi \right ) +2\pi }{N}}\\ S_{N} & =\frac {e^{i\frac {2\pi }{N}}-e^{i\frac {N\left (2\pi \right ) +2\pi }{N}}}{1-e^{i\frac {2\pi }{N}}} \end {align*}

But \(e^{i\frac {\left (N+1\right ) (2\pi )}{N}}=e^{i\frac {N(2\pi )+2\pi }{N}}=e^{i2\pi }e^{i\frac {2\pi }{N}}=e^{i\frac {2\pi }{N}}\). The above becomes\begin {align*} S_{N} & =\frac {e^{i\frac {2\pi }{N}}-e^{i\frac {2\pi }{N}}}{1-e^{i\frac {2\pi }{N}}}\\ & =0 \end {align*}

Therefore, the final result is \[ S_{N}=\left \{ \begin {array} [c]{ccc}1 & & N=1\\ 0 & & N>1 \end {array} \right . \] For geometric interpretation. Each root \(z_{n}\) is a unit vector, where the angle between each root is the same. it is \(\frac {2\pi }{N}\). Looking at each root as a vector in the complex plane, these vectors originate from the origin and end up at the unit circle, each with phase which is \(\frac {2\pi }{N}\) more than the vector just to the right of it as we go anticlockwise around the circle. The first vector starts with phase 0.

The sum \(\sum _{n=1}^{N}z_{n}\) is therefore the a vector sum of these \(N\) root. The easiest way to see that this sum is zero geometrically, is to add these vectors, by putting each vector tail, at the tip of the previous vector.  To illustrate this, we will look at the case of \(N=3\) where the angle between each vector is \(120^{0}\). This is because \(\frac {2\pi }{3}=120^{0}\). Using this method to add the roots gives this

pict
Figure 4.9:Geometric interpretation of adding the roots. Example for \(N=3\)

The above generalizes for any \(N\). If the vector sum using the tail to tip method gives a closed shape which in this case ends up back at the origin, then the vector sum is zero.

4.3.10.3 Part 3

Looking at the Electric field \(E\) at an observation point at angle \(\theta \) we obtain the following diagram

pict
Figure 4.10:Contribution of \(E\) from each slit

In the above, the \(E\) contribution from slit \(1\) was normalized to be \(E_{1}=1\). Therefore, the contribution of \(E_{2}\) from the second slit will have a phase shift relative to the first slit. This is given by \(d\sin \theta \) as seen in the diagram. For each addition slit, the phase will increase by \(d\sin \theta \). Hence the \(E_{3}\) will have phase of \(2d\sin \theta \) and so on until the last slit \(N\) which will have phase shift of \(\left (N-1\right ) d\sin \theta \).

Therefore we see that electric field at the observation point is the sum of all \(E_{n}\) from each slit, and given by\begin {align} E & =E_{1}+E_{2}+\cdots +E_{N}\nonumber \\ & =1+e^{ikd\sin \theta }+e^{ik(2d\sin \theta )}+e^{ik(3d\sin \theta )}+\cdots +e^{ik((N-1)d\sin \theta )} \tag {1} \end {align}

Now, from lecture notes, we are given the conditions for minima (i.e. destructive interference) as \begin {equation} d\sin \theta =\frac {k}{N}\lambda \qquad k=\pm 1,\pm 2,\cdots \tag {2} \end {equation} Substituting (2) into (1) gives\[ E=1+e^{ik\left (\frac {k}{N}\lambda \right ) }+e^{ik(2\left (\frac {k}{N}\lambda \right ) )}+e^{ik(3\left (\frac {k}{N}\lambda \right ) )}+\cdots +e^{ik((N-1)\left (\frac {k}{N}\lambda \right ) )}\] Replacing the first \(k\) in each term by \(\frac {2\pi }{\lambda }\) since \(k\) is wave number, then the above becomes\begin {align*} E & =1+e^{i\left (\frac {2\pi }{\lambda }\right ) \left (\frac {k}{N}\lambda \right ) }+e^{i\left (\frac {2\pi }{\lambda }\right ) (2\left (\frac {k}{N}\lambda \right ) )}+e^{i\left (\frac {2\pi }{\lambda }\right ) (3\left ( \frac {k}{N}\lambda \right ) )}+\cdots +e^{i\left (\frac {2\pi }{\lambda }\right ) ((N-1)\left (\frac {k}{N}\lambda \right ) )}\\ & =1+e^{i\frac {2\pi k}{N}}+e^{i2\left (\frac {2\pi k}{N}\right ) }+e^{i3\left (\frac {2\pi k}{N}\right ) }+\cdots +e^{i\left (N-1\right ) \left (\frac {2\pi k}{N}\right ) } \end {align*}

Let \(\phi =\frac {2\pi }{N}\,\). The above becomes\begin {equation} E=1+e^{ik\phi }+e^{i2k\phi }+e^{i3k\phi }+\cdots +e^{ik\left (N-1\right ) \phi } \tag {3} \end {equation} Comparing the above to the result obtain in part 1 we found that the sum of the roots for \(Z^{n}=1\) to be\begin {align} S_{N} & =z_{0}+z_{1}+\cdots +z_{N-1}\nonumber \\ & =e^{i\theta _{0}}+e^{i\theta _{1}}+\cdots +e^{i\theta _{N-1}} \tag {4} \end {align}

Where \begin {equation} \theta _{n}=\frac {2\pi n}{N}\qquad n=0,1,2,\cdots ,N-1 \tag {5} \end {equation} Hence (4) becomes\begin {equation} S_{N}=1+e^{i\frac {2\pi }{N}}+e^{i2\frac {2\pi }{N}}+\cdots +e^{i\left ( N-1\right ) \frac {2\pi }{N}} \tag {6} \end {equation} Therefore, for each different \(k\) Eq(3) is the same as (6). So (3) can be written as\begin {equation} E=1+z+z^{2}+\cdots +z^{N-1} \tag {7} \end {equation} Where now \(z=e^{ik\phi }\) with \(\phi =\frac {2\pi }{N}\). But we know that \begin {equation} \left (1-z^{N}\right ) =\left (1-z\right ) \left (1+z+z^{2}+\cdots +z^{N-1}\right ) \tag {8} \end {equation} But \(\left (1-z^{N}\right ) =0\) since \(1\) is root of \(z^{N}\). Hence the above becomes\[ 0=\left (1-z\right ) \left (1+z+z^{2}+\cdots +z^{N-1}\right ) \] Since \(z\neq 1\) (unless \(\frac {2\pi }{N}k\) happened to be exact multiple of \(2\pi \)), then we conclude that \(1+z+z^{2}+\cdots +z^{N-1}\) must be zero. This implies that\begin {align*} E & =1+e^{ik\phi }+e^{i2k\phi }+e^{i3k\phi }+\cdots +e^{ik\left (N-1\right ) \phi }\\ & =0 \end {align*}

Under the condition of destructive interference. This says the total Electric field from the \(N\) slits will vanish at the observation point when destructive interference condition is applied. Which is what we are asked to show.

4.3.11 key solution for HW 3

PDF