modeling stochastic sequence with determined mean, std and autocorrelation function

I have a problem to generate series of pseudo-stochastic numbers are representing some kind of distribution which are determined by their Mean, STD and autocorrelation function let say R=cos(at+b)*exp(-ct)