Taking the first derivatives of (9.12) w/r to the amplitude parameters and , we have
Taking the first derivative of (9.12) w/r to frequency , we get
We can also obtain the Fisher's Information matrix by taking the negative expected value of the second derivatives. From the above, we can obtain
All these quantities can be computed in the frequency domain with help of the FFT. Multiplication by matrix is the same as multiplying Fourier coefficient by . The ML estimation approach of Section 17.5 can be used. However, some parameters (the amplitudes ) can be determined in one step if the frequencies are assumed: