If you want to perform a Phase Shift in time domain you need only to multiply your wavelet by exp(i*ang) where ang is the degree of the phase shift in radians.
Equation 2: (Ricker wavelet shifted by 45º) s.*exp(i*pi/4)
Noticed the usage of the dot product, because this is a pointwised multiplication.