p = quadtfdist(f, q);
For an input vector of length L, the kernel should be a L x L matrix. quadtfdist(f, q); computes a discrete quadratic time-frequency distribution.