How to plot the 2d complex function with amplitude and phase

I have a 2d complex function as follows

psi=sqrt(x^2+y^2)*exp(-(x^2+y^2)/16)*exp(i*atan2(y,x));

How to realize the following figure by using Matlab?

In this image, the color on the ring presents the changes of phase of the function from small to big.

Answer 1

[Cris LaPierre]

Perhaps you would be interested in this File Exchange submission:

https://www.mathworks.com/matlabcentral/fileexchange/66850-complex-function-visualization-amplitude-phase

[Qiu Xu]

Thanks for your helps