Generating a phase modulation frequency spectrum with jitter

I would like to produce a frequency spectrum plot like this: