trigBasisExp function

Compute the design matrix for a trigonometric basis expansion