I am trying to create a set of wavelets in frequency space--namely Cauchy
wavelets for an intensity analysis (von Tscharner, 2000). The wavelets are
defined by the following formula:
[(f/cf)^(cf*scale)]*[e^((-f/cf)+1)^(cf*scale)]
where *f *is frequency of length *n*, *cf* is center frequency (defined
below) and is an array of *j *columns and *n* row, and scale is a constant.
cf = (1/scale)*(j + q)^r
where *j *is the number of center frequencies (i.e., columns), and *q* and *
r* are constants.
I seem to be doing fine with most of the math, but the (*array*)^(cf*scale)
eludes me. It seems like it should be simple, but nothing is working. Any
suggestions??
[[alternative HTML version deleted]]