Taking the first derivatives of (8.12) w/r to the amplitude parameters and , we have

Taking the first derivative of (8.12) w/r to frequency , we get

where is the matrix

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 16.5 can be used. However, some parameters (the amplitudes ) can be determined in one step if the frequencies are assumed:

So, it is only necessary to iterate on the frequencies.

Baggenstoss 2017-05-19