Replace ordinate with its phase angle in radians.
g=phase(f)
This function replaces all ordinate values of f with the complex phase angle in radians. The phase angle will be in the range from -2*pi to 0.
imat_fn/abs, imat_fn/phased