Generate the coefficients in the equations that determine the phase constraints.

The result will be a structure containing:

A

A matrix [Nf x Nb], where Nf is the number of nonnegative frequencies in the Fourier transform and Nb is the number of basis functions (i.e., EOFs). A[k,j] = (a_k^i a_k^j), where i is the index of the basis function used for the reference (conventionally 2, which is the first spatial basis function. Don't use the time-like basis function, i=1, since it will generally have small coefficients).

i

The index of the reference basis function.

phase_eqn_coef(Fx, i = 2)

Arguments

Fx

The Fourier transform of the EOF projection coefficients.
i

The index of the basis function to use for the reference. Don't use EOF-0 (i=1).

Details

These equations are given as equation XX in the paper.