MATLAB: Generating a random bit stream in a frequency Band for PSK Modulation
You are attempting to simulate an RF signal directly at it’s carrier frequency which will require much more processing so is not recommended. You can achieve the same results by simulating everything at its equivalent complex baseband signal due to the linear properties of the underlying frequency translation. Equivalently everything in your simulation model (IF bands, channel models etc) can/should all be at complex baseband equivalent representation (the analytic signal).

That said, to do what you are trying to accomplish with either method, you would first create the complex modulation signal mapping your data bits to PSK symbols; for example for 8 PSK 3 bits map to every complex signal such as 000 could become $e^{j\pi/4}$, 111 could become $e^{j5\pi/4}$, etc- any mapping you choose but typically “grey coded” such that adjacent symbol locations only change by one bit.

This sequence which is then sampled only 1 sample per symbol then needs to be up sampled to support implementing pulse shaping (2 samples per symbol would be sufficient but in implementation more samples here can simplify later filtering). 

This now represents the complex baseband signal with pulse shaping and if 2 samples per symbol the represented spectrum from a PSD would extend from $-f_S/2$ to $+f_S/2$ where $f_S$ is the symbol rate. (Thus If oversampling more this frequency range can be extended further).

This shows the SIGNIFICANT advantage to doing all simulations at baseband - to now continue to to simulating at a carrier frequency at 900 MHz directly you would need to increase the sampling rate to be greater than the carrier frequency plus half the signal bandwidth plus some margin if completely a complex signal or twice that if simulating a real signal such that you want to model each carrier cycle. There is no reason to do this since everything that carrier is doing can be represented at complex baseband.

But if you still wanted to, against my advice, you would first interpolate the baseband waveform to this much higher rate and then multiply the signal by the carrier frequency ($e^{j\omega_c t}$ if complex or $\cos(\omega t)$ if modeling a real carrier).