Is there a way to simulate time series data with a specific rolling mean and autocorrelation in R?

Question

I have an existing time series (1000 samples) and calculated the rolling mean using the filter() function in R, averaging across 30 samples each. The goal of th