r = drihaczekdist(f);
drihaczekdist(f) computes a discrete Rihaczek distribution of vector f. The discrete Rihaczek distribution is computed by
where \(k, l=0,\ldots,L-1\) and \(c\) is the Fourier transform of \(f\).
WARNING: The quadratic time-frequency distributions are highly redundant. For an input vector of length L, the quadratic time-frequency distribution will be a \(L \times L\) matrix. If f is multichannel (\(L\times W\) matrix), the resulting distributions are stacked along the third dimension such that the result is \(L\times L \times W\) cube.