Comment by stagger87

If you change your Fourier transform approach into a time domain convolution, then you can generate it sample by sample. White noise -> FIR filter, this is a pretty simple and direct algorithm as far as DSP goes. The wiki on pink noise mentions this. I've done the FIR filter approach to simulate 1/f phase noise.