I am not aware of a specific phasor notation in MATLAB. Since they are a convenient way of representing complex numbers, leave them as complex variables for calculations and convert them to phasor notation when you output them:
z = 3 + j*5
phasorr = @(z) [abs(z) angle(z)]; % angle in radians
phasord = @(z) [abs(z) atan2d(imag(z),real(z))]; % angle in degrees
phasorr(z)
phasord(z)
I used atan2d in the degree version because in my experience, they do more accurate conversions than multiplying radian angles by (180/pi).
