plotting phase transformations of a signal I would like to find the phase variations of the signal with respect to time axes:

S = sin(2*pi*100*t + pi/4) + cos(pi*500*t) + sin(2*pi*100*t + 5*pi/2);
My attempt to get the phase, see figure below: